Линейная интерполяция с функцией Криге в R - PullRequest
Новая
       12

Линейная интерполяция с функцией Криге в R

0 голосов
/ 15 апреля 2019

Я делаю регрессионный кригинг в R и уже сделал обычный кригинг с остатками, который работал абсолютно нормально.Теперь я хочу выполнить линейную регрессию с тремя переменными прогнозирования, используя krige. 

Slope <- readGDAL("Slope.tiff")
proj4string(Slope) <- CRS(paste("+init=epsg:",32632,sep=""))
SRTM <- readGDAL("SRTM_NiSa.tiff")
proj4string(SRTM) <- CRS(paste("+init=epsg:",32632,sep=""))
Wetness <- readGDAL("Wetness.tiff")
proj4string(Wetness) <- CRS(paste("+init=epsg:",32632,sep=""))

...

VarioTon <- variogram(Lucas@data$TonKorrigiertFinal~1, LucasTransformiert, width = 5000)
plot(VarioTon, type = "b", main = "Experimentelles Variogramm (Ton)")
vmTon <- vgm(100, "Exp", 12000, 10)
vmfTon <- fit.variogram(VarioTon, vmTon)

step_ton <- step(model_ton)
geo_data@data$step_tonfit <- fitted(step_ton)
geo_data@data$step_tonres <- residuals(step_ton)

Residuals_Ton <- krige(geo_data@data$step_tonres~1, 
                       locations = geo_data, newdata = SRTM, model= vmfTon)

krige(geo_data@data$step_tonfit ~ Slope+SRTM+Wetness, locations = geo_data, newdata = Slope)

Это приводит к сообщению об ошибке «Недопустимый тип (S4) для переменной« SRTM »».SRTM - это формальный класс SpatialGridDataframe, остальные - Large SpatialGridDataFrames, хотя все они имеют одинаковый экстент и были открыты с помощью readGDAL.Выход из SRTM приводит к появлению сообщения об ошибке «объект не является матрицей».Только использование наклона приводит к тому же сообщению об ошибке.Дополнительные предупреждения об ошибках касаются разного количества строк.

...