博文
Python+ArcPy: Extract by Mask (Spatial Analyst)
|||
原帖ENVI+IDL：Clip a Raster by a shapefile出现一个Bug，当shapefile文件覆盖范围大于Raster数据时程序出现错误，裁剪不成功，替代方案是基于ArcGIS的ArcPy再编辑一个Python版的裁剪程序。
代码参考：Extract by Mask (Spatial Analyst)，运行流畅，成功输出预期结果。
两点需要特别注意：
[1] 代码似乎不能定义裁剪背景的数值，默认是-32768，在ArcGIS中显示为NoData。
[2] 注意原始栅格文件的背景可能会被包含在裁剪结果之中，图 1上方黑色的区域即为原始栅格文件的无效区域，后期的操作应对这种区域加以处理，避免异常数值对分析结果的干扰。
图 1
附上测试数据及代码（
pyclip.rar）。
https://wap.sciencenet.cn/blog-1148346-877265.html
上一篇：Maths: The Negative Values of R2
下一篇：Matlab: Mosaicking Images
扫一扫，分享此博文