Я использую пакет geoR
для пространственной интерполяции осадков. Должен сказать, что я новичок в геостатистике. Благодаря некоторым видеоурокам на YouTube, я понял (ну, я так думаю) теорию вариограммы. Насколько я понимаю, количество пар должно уменьшаться с увеличением расстояния задержки. Например, если мы рассматриваем участок длиной 100 м (скажем, поперечное сечение русла реки длиной 100 м), количество пар для 5-метрового лага равно 20, а количество пар для 10-метрового лага равно 10 и так далее. Но меня немного смущает вывод функции variog
в geoR
package. Пример приведен ниже
mydata
X Y a
[1,] 415720 432795 2.551415
[2,] 415513 432834 2.553177
[3,] 415325 432740 2.824652
[4,] 415356 432847 2.751844
[5,] 415374 432858 2.194091
[6,] 415426 432774 2.598897
[7,] 415395 432811 2.699066
[8,] 415626 432762 2.916368
это мой набор данных, где a
- моя переменная (интенсивность дождя), а x, y
- координаты точек. Расчет вариграммы показан ниже.
geodata=as.geodata(data,header=TRUE)
variogram=variog(geodata,coords=geodata$coords,data=geodata$data)
variogram[1:3]
$u
[1] 46.01662 107.37212 138.04987 199.40537 291.43861 352.79411
$v
[1] 0.044636453 0.025991469 0.109742986 0.029081575 0.006289056 0.041963076
$n
[1] 3 8 3 3 3 2
куда
u: вектор с расстояниями.
v: вектор с оценочными значениями вариограммы на расстояниях, указанных в u.
n: количество пар в каждой ячейке
В соответствии с этим количество пар (n) имеет случайный образец, тогда как соответствующее расстояние запаздывания (u) увеличивается. Мне трудно это понять. Кто-нибудь может объяснить, что происходит? Также мы высоко ценим любые предложения / советы по улучшению расчета вариограммы для этого приложения (пространственная интерполяция интенсивности дождя), поскольку я новичок в геостатистике. Заранее спасибо.