ArcGIS中的工具均提供Python函数,方便支持批量处理。
以掩膜裁剪(提取分析/Extract by Mask)为例,
![50d4e2ce2dc619e997959bece3f4d46c.png](https://i-blog.csdnimg.cn/blog_migrate/862ec7791c660cef396101bdc1f871d3.jpeg)
Python资源共享群:484031800
可以查找到帮助文档“按掩膜提取”(源1)(源2),路径“工具参考/工具/Spatial Analyst 工具箱/按掩膜提取”函数为:Extract by Mask
Extract by Mask的基本语法是:
ExtractByMask (in_raster, in_mask_data)
in_raster:提取像元的输入栅格。
in_mask_data:用于定义提取区域的输入掩膜数据。它可以是栅格,也可以是要素数据集。
代码示例
以下以批量读取裁剪某文件夹下的tif数据为例(特别注意Python2脚本循环语句格式的空格不能用Tab,否则会报错Indent;
另外在使用不同的矢量裁剪同一影像时,为保证裁剪后范围一致需设置范围为输入栅格范围。例如下面的:
rd = arcpy.sa.Raster(inRaster)
# Set the extent environment as the raster, very important for clip with different vector