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
本文出自 “记笔记” 博客，请务必保留此出处http://dressame.blog.51cto.com/12551541/1918863
原文：http://dressame.blog.51cto.com/12551541/1918863