利用shp文件构建mask【MATLAB和ARCGIS】两种方法
1 ARCGIS （推荐！！！-速度很快）
- 利用Polygon to Raster
注意：由于我们想要的mask有效值是1，在进行转换的时候，注意设置转换字段【Value field】
【Value field】通过编辑shp文件属性表进行设置【设置为1】
可通过Environments设置经纬度范围
- 得到栅格数据后，导出栅格
通过设置为NODATA=0将目标矢量区域外的格网值设置为0
可以自定义栅格分辨率或者栅格大小
这里，我生成的是tif数据，得到最终的mask栅格数据，方便后期在Python或者MATLAB中使用
2 MATLAB 方法（速度太慢）
利用函数shaperead和inpolygon
clear;clc;
# 1 读取shp文件
China1=shaperead('省.shp');
# 2 设置整体区域【一般是一个大格网】的经纬度
mask = ones(128,128);
lat = 30.8:0.05:37.15;
lon = 110.35:0.05:116.7;
# 注意：mask的行列数量需要和lat[lon]保持一致
[latt,lonn] = meshgrid(lat,lon); # 生成经纬度格网点
# 3 mask圈定
# idx1为正整数
inshp = inpolygon(lonn,latt,China1(idx1).X,China1(idx1).Y);
mask = mask.*inshp; # 点乘
mask(mask==0) = nan;
学习链接
- 【matlab利用shp文件制作mask白化文件】