查看原文
其他

利用Python在ArcGIS中批量裁剪影像(带重叠度)「附源代码下载」

王双炎 GIS前沿 2022-03-16


之前在在CAD平台上做过将分幅图框外扩一个插件「C# | 批量将CAD图幅网格外扩生成新图框(附源代码下载)」,将50*50网格进行外扩一定距离,使用外扩后的新图框进行裁剪图形或影像,这样就会有重叠部分、方便工作过程中使用。

下面我们外扩后的新图框结合Arcpy对影像进行裁剪。

工作流程

1. CAD网格和文字注记转成shp文件;
2. 对转换完成后的网格面由文字注记点进行赋值,将图幅号赋值到面中,用于命名影像; 
3. 使用ArcPy对影像进行裁剪;
4. 叠加裁剪后的影像和矢量数据检查。

如下图CAD图形所示,白色图形为标准50*50分幅图框,红色图形是白色图形外扩10米后的新图框。
 

CAD分幅图框

将CAD图形标准分幅图框、外扩后的图框和文字注记分别转成shp文件,并把注记赋值到面属性中。CAD注记转shp后图幅号在【HTEXT】中。
 

文字转shp

外扩后的图框转shp,并新建一个图幅号字段。
 

网格转shp

点向面赋值,使用[ArcToolbox]—>[分析工具] —>[叠加分析] —>[空间连接]工具。
 

赋值前
 

空间连接

空间连接后结果,使用字段计算器将【HTBXT】字段值赋值到【TFH】上。
 

空间连接后结果

创建ArcPy脚本
 

ArcPy裁剪

脚本执行过程
 

裁剪后结果

使用ArcMap加载裁剪后的影像,可以看出在不同的图幅之间、裁剪后的影像保留了一部分重叠的部分。
 

裁剪后影像

进入本公众号后台

发送数字

5211

获取下载链接


 - END -


Think LiDAR in Python
C# | CAD批量导出多段线节点坐标(附源代码下载)
Python学的好,数据处理快又好!ArcPy入门到精通(手册可下载)
基于Python的卫星影像数据处理与分析
Python在GIS中的应用实例


戳原文,更有料!

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存