Это продолжение этого вопроса . Теперь мне удалось построить вариограммы данных, используя ответы данный вопрос Код, который я использую (где im
- это переменная, содержащая данные dput
, указанные в https://gist.github.com/2780792is
library(gstat)
point_data <- as(im, 'SpatialPointsDataFrame')
v <- variogram(band1 ~ 1, data = point_data)
plot(v)
Это дает мне следующий график:
Как видите, числа на оси расстояний варьируются от 0 до примерно 150. Однако единицами данных являются градусы:
> head(coordinates(point_data))
x y
[1,] -1.849353 52.06165
[2,] -1.759728 52.06165
[3,] -1.401227 52.06165
[4,] -1.311602 52.06165
[5,] -1.221977 52.06165
[6,] -1.132352 52.06165
Учитывая это, каковы единицы измерения расстояния по оси X? Я ожидаю, что они будут в тех же единицах, что и координаты. Для них не имеет смысла быть метрами, я думаю, они могут быть километрами.
Я не могу найти ничего об этом в руководстве по gstat, поэтому я предположил, что они использовали базовые единицы измерения данных, но, похоже, это не так.
Есть идеи?