通过掩码提取的R保留所有原始数量的单元格

时间:2016-04-27 22:24:44

标签: r raster

我试图在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)

enter image description here

#extract raster based on the shape file
dat<-mask(rast$hd85tn501, shape)
plot(dat)

enter image description here

正如您所看到的,尽管代码有效,但我的最终栅格太小了。如果我输入dat,它仍会显示单元格总数为1944000.我想要的是保留并仅显示我想要提取的单元格而不是整个1944000单元格。

编辑:我也试过这个,但它给了我一些我不需要的额外网格(mask能够做到的)。

raster::crop: plot(crop(rast$hd85tn501, shape))

enter image description here

0 个答案:

没有答案