新民网站建设价格咨询巨量算数

1 ARCGIS （推荐！！！-速度很快）

1. 利用Polygon to Raster
   注意：由于我们想要的mask有效值是1，在进行转换的时候，注意设置转换字段【Value field】
   【Value field】通过编辑shp文件属性表进行设置【设置为1】
   
   可通过Environments设置经纬度范围
   
2. 得到栅格数据后，导出栅格
   
   通过设置为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白化文件】