ArcGis Python腳本——根據接圖表批量裁切分幅影像

插個廣告，制作ArcGIS的Tool工具學習下面的教程就對了：

零基礎學習Python制作ArcGIS自定義工具觀看鏈接

《零基礎學習Python制作ArcGIS自定義工具》課程簡介

---

 

年前寫了一個用漁網工具制作圖幅接圖表的文章，鏈接在這里：

使用ArcMap做一個1:5000標准分幅圖並編號

本文提供一個使用ArcMap利用接圖表圖斑裁切一幅影像為多幅的方法。

第一步，將接圖表拆分，使之每一個圖斑為一個要素類（Shapefile文件）

拆分方法很多，提供兩個參考：

1、可以使用Python腳本

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

拆分的結果是每個影像框圖斑為一個Shapefile文件，所有Shapefile文件生成在一個文件夾下。

第二步，使用Python腳本遍歷每個Shapefile文件，以其幾何或外接矩形裁切影像

import arcpy
#設置工作空間，即裁切參照shp文件所在的目錄
arcpy.env.workspace=r'A:/1_5000/00/shp'
#獲取工作空間中的Shapefile文件列表
featureClasses=arcpy.ListFeatureClasses()
#待裁切影像文件全路徑，這里是.img格式
inRasterPath=r'A:/1_5000/GF2-2199394-20170221/GF2219939420170221F.img'
#裁切出的影像文件輸出文件夾路徑
outFolderPath=r'A:/1_5000/00/tif/'
for fc in featureClasses:
    arcpy.Clip_management(inRasterPath,'#',outFolderPath+fc[:-4]+'.tif',fc,0,"ClippingGeometry")
    #arcpy.Clip_management函數的用法可以參考官方幫助
    #outFolderPath+fc[:-4]+'.tif'為輸出文件的全路徑，tif格式，用切割參照shp的名字命名，fc[:-4]即截斷名字最后4個字符，即.shp

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

arcpy.Clip_management函數幫助頁 http://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/data-management-toolbox/clip.htm