我试图在R中使用掩码函数提取。这就是我所做的:
library(raster)
# download my shape file
shape <- getData('GADM', country='IND', level=0)
#download my raster file
rast<-getData('CMIP5', var='tmin', res=10, rcp=85, model='HD', year=50)
plot(rast)
#extract raster based on the shape file
dat<-mask(rast$hd85tn501, shape)
plot(dat)
正如您所看到的,尽管代码有效,但我的最终栅格太小了。如果我输入dat
,它仍会显示单元格总数为1944000.我想要的是保留并仅显示我想要提取的单元格而不是整个1944000单元格。
编辑:我也试过这个,但它给了我一些我不需要的额外网格(mask
能够做到的)。
raster::crop: plot(crop(rast$hd85tn501, shape))