1、從http://www.gadm.org/country 下載China和Taiwan的R格式的level 0地圖
CHN_adm0.rds, TWN_adm0.rds
2、蒙板
r <- raster(nrow=1201,ncol=1401,ext=extent(70,140,0,60))
chn_mask <- rasterize(readRDS("e:/CHN_adm0.rds"),r)
twn_mask <- rasterize(readRDS("e:/TWN_adm0.rds"),r)
#或操作,合併中國大陸和臺灣
all_mask <- twn_mask | chn_mask
3、業務相關
產品柵格先行後列,從左向右(70到140),從下向上排列(0到60)
#反轉行序
flip_all_mask <- flip(all_mask,'y')
#轉爲矩陣
mask_matrix <- as.matrix(flip_all_mask)
4、蒙板使用
#柵格產品加載爲矩陣m
m[is.na(mask_matrix)] <- NA