ArcGis Python脚本——根据接图表批量裁切分幅影像 使用ArcMap做一个1:5000标准分幅图并编号

年前写了一个用渔网工具制作图幅接图表的文章,链接在这里:

使用ArcMap做一个1:5000标准分幅图并编号

本文提供一个使用ArcMap利用接图表图斑裁切一幅影像为多幅的方法。

第一步,将接图表拆分,使之每一个图斑为一个要素类(Shapefile文件)

拆分方法很多,提供两个参考:

1、可以使用Python脚本

2、可以使用ArcToolbox——分析工具——split

拆分的结果是每个影像框图斑为一个Shapefile文件,所有Shapefile文件生成在一个文件夹下。

第二步,使用Python脚本遍历每个Shapefile文件,以其范围extent裁切影像

 1 import arcpy
 2 #设置工作空间,即裁切参照shp文件所在的目录
 3 arcpy.env.workspace=r'A:/1_5000/00/shp'
 4 #获取工作空间中的Shapefile文件列表
 5 featureClasses=arcpy.ListFeatureClasses()
 6 #待裁切影像文件全路径,这里是.img格式
 7 inRasterPath=r'A:/1_5000/GF2-2199394-20170221/GF2219939420170221F.img'
 8 #裁切出的影像文件输出文件夹路径
 9 outFolderPath=r'A:/1_5000/00/tif/'
10 for fc in featureClasses:
11     arcpy.Clip_management(inRasterPath,'#',outFolderPath+fc[:-4]+'.tif',fc,0,"ClippingGeometry")
12     #arcpy.Clip_management函数的用法可以参考官方帮助
13     #outFolderPath+fc[:-4]+'.tif'为输出文件的全路径,tif格式,用切割参照shp的名字命名,fc[:-4]即截断名字最后4个字符,即.shp

arcpy.Clip_management函数帮助页 http://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/data-management-toolbox/clip.htm

猜你喜欢

转载自www.cnblogs.com/yzhyingcool/p/10486352.html