Я пытаюсь использовать функцию кригинга для создания вертикальных карт химических параметров на разрезе океана, и мне трудно начать.
Мои данные выглядят так:
horiz=rep(1:5, 5)
depth=runif(25)
value = horiz+runif(25)/5
df <- data.frame(horiz, depth, value)
Функция autoKrige
в пакете automap
выглядит так, как будто она должна выполнять эту работу за меня, но она принимает объект класса SpatialPointsDataFrame
. Насколько я могу судить, функция spTransform
в пакете rgdal
создает SpatialPointsDataFrame
объектов, но есть две проблемы:
Двоичные файлы OSX этого недоступны из CRAN, и моя копия RStudio, работающая на OXS 10.7, похоже, не может ее установить, и
Эта функция, по-видимому, работает с данными широты/долготы и корректирует значения расстояния для кривизны Земли. Поскольку я имею дело с вертикальной плоскостью (и короткими расстояниями, масштабом в сотни метров), я не хочу исправлять свои расстояния.
В R здесь есть отличное обсуждение кригинга, но из-за проблем, перечисленных выше, я не совсем понимаю, как применить это к моей конкретной проблеме.
Мне нужна матрица или кадр данных, описывающий сетку точек с интерполированными значениями для моих химических параметров, которые я затем могу построить (в идеале с помощью ggplot2). Я подозреваю, что решение моей проблемы значительно проще, чем я это представляю.