0

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を取得する方法について何か提案はありますか? ありがとう!

4

1 に答える 1