0

パッケージ 'maptools' のデータによって提供される世界地図の国境に従ってトリミングしたいラスター グリッドがあります。crop()いくつかの再検索を行うことで、関数を使用してから関数を使用する必要があることがわかりましたmask()が、エラーメッセージが表示されます。

これが私のコードです:

# load the worldmap SpatialPolygonDataFrame
library(maptools)
data(wrld_simpl)
ll=CRS("+init=epsg:4326")
world<-spTransform(wrld_simpl, ll)
ext<-extent(-10.417,31.917,34.083,71.083)
# get only region covering europe 
world@bbox<-as.matrix(ext)


# create the origin in WGS84 CRS
ll = CRS("+init=epsg:4326") 
origin = SpatialPoints(cbind(7,40),ll) 
# create the raster grid 
grid = raster() 
# define extent and resolution
e = extent(5,18,40,50)
extent(grid)<-e
res(grid)=c(0.08,0.08) 
# fill the raster values
grid[] = runif(1:20250)

# crop the grid with worldmap to erase sea areas
test<-crop(grid,world)
rasterize(world,test)
Found 246 region(s) and 3768 polygon(s)
Error in if (x1a > rxmx) { : argument is of length zero
4

1 に答える 1

0

あなたが思うように、トリミングは国のポリゴンにトリミングされません。実際、crop(grid, world)ワールドの範囲はグリッドの範囲よりも大きいため、実際には何もトリミングされません。プロットすることで確認できますtest。と同じgridです。

代わりに、次のgrid[] = runif(1:ncell(grid))ことができた後：

# we're going to make 3 plots
par(mfrow=c(1,3))

# crop the Worldmap to grid 
worldcrop<-crop(world, grid)
plot(worldcrop)

# rasterize output, give cells value of NAME(seas are NA)
worldcropr = rasterize(worldcrop,grid, field='NAME', fun='first')
plot(worldcropr)

# mask random grid by worldcropr
gridmask = mask(x=grid, mask=worldcropr)
plot(gridmask)

# number of land pixels:
sum(!is.na(gridmask@data@values))

# number of sea pixels:
sum(is.na(gridmask@data@values))

私はそれがあなたが求めているものだと信じています。 トリミングされたグリッド

于 2014-05-23T13:00:56.677 に答える