时间管理系统 网站开发,百度关键词优化推广,网站集群建设合同,Wordpress请求接口数据1 ARCGIS #xff08;推荐#xff01;#xff01;#xff01;-速度很快#xff09;
利用Polygon to Raster 注意#xff1a;由于我们想要的mask有效值是1#xff0c;在进行转换的时候#xff0c;注意设置转换字段【Value field】 【Value field】通过编辑shp文件属性表…1 ARCGIS 推荐-速度很快
利用Polygon to Raster 注意由于我们想要的mask有效值是1在进行转换的时候注意设置转换字段【Value field】 【Value field】通过编辑shp文件属性表进行设置【设置为1】 可通过Environments设置经纬度范围 得到栅格数据后导出栅格 通过设置为NODATA0将目标矢量区域外的格网值设置为0 可以自定义栅格分辨率或者栅格大小 这里我生成的是tif数据得到最终的mask栅格数据方便后期在Python或者MATLAB中使用
2 MATLAB 方法速度太慢
利用函数shaperead和inpolygon
clear;clc;
# 1 读取shp文件
China1shaperead(省.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(mask0) nan;学习链接
【matlab利用shp文件制作mask白化文件】