與autoKrige功能投影錯誤

Question

我有autoKrig功能的問題，我試圖在這裏做重複的例子：

library(automap) 
library(raster) 
library(dismo) 

bio <- getData("worldclim", var="bio", res=10) 

bio1 <- raster(bio, layer=1) 
bio12 <- raster(bio, layer=12) 
predictors <- stack(bio1, bio12) 

bg <- randomPoints(bio1, 50) 

data <- extract(predictors, bg) 
data <- cbind(bg,data) 
data <- data.frame(data) 

coordinates(data)=~x+y 
proj4string(data) = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0") 
kg <- autoKrige(bio1~bio12, data, new_data=predictors) 

這將導致：

Error in autoKrige(bio1 ~ bio12, data, new_data = predictors) : 
    Either input_data or new_data is in LongLat, please reproject. 
    input_data: +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
    new_data: +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 

我得到同樣的錯誤與我的原始數據。我感謝任何幫助。

Answer 1

如果您閱讀幫助文件，它會告訴您爲什麼它會拋出該錯誤。

autoKrige對input_data 和new_data的座標系執行一些檢查。如果其中一個或兩個都是NA，則分配其他投影 。如果他們有不同的預測，則會產生錯誤。 如果其中一個或兩個都有非投影系統（即緯度 - 經度），則會引發錯誤。由於'gstat在地理座標中使用 球面距離，但是 通常的變異函數模型通常不是非負確定球體，並且沒有合適的模型可用'（Edzer Pebesma on r -sig-GEO）。

看起來您需要在調用autoKrige之前投影您的數據。