R_номер пар для каждого лага в вариограмме
Я использую geoR
пакет для пространственной интерполяции осадков. Должен сказать, что я новичок в геостатистике. Благодаря некоторым видеоурокам на YouTube я понял (ну, я так думаю) теорию вариограммы. Насколько я понимаю, количество пар должно уменьшаться с увеличением лагов. Например, если мы рассмотрим отрезок длиной 100 м (скажем, поперечное сечение русла длиной 100 м), то количество пар для отставания в 5 м составляет 20, а количество пар для отставания в 10 м - 10 и так далее. Но я немного запутался с выходом из variog
функция в geoR
пакет. Пример приведен ниже
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
где
и: вектор с расстояниями.
v: вектор с оцененными значениями вариограммы на расстояниях, указанных в и.
n: количество пар в каждой корзине
В соответствии с этим количество пар (n) имеет случайный характер, тогда как соответствующее расстояние запаздывания (u) увеличивается. Мне трудно это понять. Кто-нибудь может объяснить, что происходит? Кроме того, любые предложения / советы по улучшению расчета вариограммы для этого приложения (пространственная интерполяция интенсивности осадков) высоко ценится, поскольку я новичок в геостатистике. Заранее спасибо.
1 ответ
На линейном разрезе 100 м с регулярным интервалом 5 м между наблюдениями, если у вас будет 20 пар с отставанием в 5 м, у вас будет 19 пар с отставанием в 10 м. Эта идея не распространяется на ваши данные, потому что они распределены нерегулярно и распределены по двум измерениям. Для нерегулярно распределенных данных у вас очень мало пар точек для очень коротких расстояний. Совет для получения более качественной вариограммы заключается в работе с большим набором данных: геостатистика начинает становиться интересной с 30 наблюдениями, а с более чем 100 наблюдениями.