1

この例は再現可能で、最初の部分 (大きなポリゴンを含む) は機能しますが、2 番目の部分は失敗します。

library(raster)

USA.altitude <- getData('alt', country='USA')
lower48 <- USA.alt[[1]]

Srs1 = Polygons(list(Polygon(cbind(c(-96,-95,-95,-96),c(40,40,41,40)))), "s1")
Srs2 = Polygons(list(Polygon(cbind(c(-97,-96,-96,-97),c(40,40,41,40)))), "s2")

spdf.large <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs1,Srs2)), 
                              data.frame( z=1:2, row.names=c("s1","s2") ) )
plot(spdf.large)
usa.average.elev<-extract(lower48, spdf.large, fun=mean)

わーい!はい、大丈夫ですよね?小さなポリゴンで:

Srs3 = Polygons(list(Polygon(cbind(c(-95.01,-94.99,-94.99,-95.01),c(40,40,40.01,40)))), "s3")
Srs4 = Polygons(list(Polygon(cbind(c(-94.01,-93.99,-93.99,-94.01),c(40,40,40.01,40)))), "s4")


spdf.small <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs3,Srs4)), 
                              data.frame( z=1:2, row.names=c("s3","s4") ) )
plot(spdf.small)
usa.average.elev<-extract(lower48, spdf.small, fun=mean)

エラーで失敗します:

pmax(ceiling(10/dim(r)[1:2]), c(1, 1)) のエラー: オブジェクト 'r' が見つかりません

4

0 に答える 0