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

R 中国陆地蒙板