1
我有autoKrig功能的問題,我試圖在這裏做重複的例子:與autoKrige功能投影錯誤
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
我得到同樣的錯誤與我的原始數據。我感謝任何幫助。
謝謝,但由於錯誤顯示我的輸入和輸出都有相同的投影。 –
@ N16您的數據在'longlat'投影中,幫助文件詳細描述爲「非投影系統(即緯度 - 經度)」。 – ialm
所以在這種情況下,你推薦什麼?我的數據集是全球性的。 –