Как мне оценить вариограмму для данных о глобусе?

Question

У меня есть данные с соответствующими долготами и широтами. Как получить вариограмму для этих данных на основе расстояний между точками большого круга?

Этот простой пример содержит все данные на экваторе:

require(geoR)

long <- seq(-179, 180)
x <- sin(pi * long / 180) + rnorm(length(long))
V <- variog(data=x, coords=cbind(long, 0))
# variog: computing omnidirectional variogram
plot(V)

Первая и последняя точки на самом деле находятся на расстоянии всего лишь 1 градуса, но моя наивная попытка приводит к variog мнению, что они разделены на 359 градусов.

Answer 1

Вы должны использовать полувариограмму от nmle . Позволяет указать матрицу расстояний, которую вы можете самостоятельно разработать для себя.

Answer 2

Из поста на R-sig-geo (список рассылки, посвященный пространственным данным в R) я, кажется, помню, что в R нет готовых к работе функций, которые поддерживают расстояния большого круга:

http://r -sig-geo.2731867.n2.nabble.com / Great-Circle-distance-in-Automap-Gstat-td6863940.html

Мое предложение будет проецировать вашданные, а затем выполнить интерполяцию на проецируемых данных.