ArcGis Python脚本

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

使用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

# code source： https://www.cnblogs.com/yzhyingcool/
# QQ：975601416

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