landsat パッケージのslopeasp()関数を使用して勾配の値を取得するために、SpatialPointsDataframeをSpatialGridDataframeに変換しようとしています。SpatialPointsDataframe をラスターに変換して変換しようとしましたが、これも機能せず、エラーが発生しました。
as(pts, "spatialGridDataFrame") のエラー: "SpatialPointsDataFrame" を "spatialGridDataFrame" に強制するためのメソッドまたはデフォルトはありません</p>
以下は私のコードです:
#create fake data
x <- rep(rep(seq(12,36,0.5),41))
y <-rep(seq(32,52,0.5), each=49)
z <- rnorm(2009, 26.5, 44.0)
pts <- as.data.frame(matrix(c(x,y,z), ncol=3, byrow=FALSE))
colnames(pts)=c("x", "y", "z")
# create a SpatialPointsDataFrame
coordinates(pts) = ~x+y
# create an empty raster object to the extent of the points
rast <- raster(ext=extent(pts), resolution=250)
# rasterize irregular points
rasOut<-rasterize(pts, rast, pts$z, fun = mean)
#attempt from spdf to sgdf
sgdf <- as(pts, 'spatialGridDataFrame')
#atempt from raster to sgdf
sgdf <- as(rasOut, 'spatialGridDataFrame')
SpatialGridDataFrameを取得する方法について何か提案はありますか? ありがとう!