首页 > 编程语言 > 详细

ArcGis Python脚本——根据接图表批量裁切分幅影像

时间:2019-03-06 22:35:01      阅读:519      评论:0      收藏:0      [点我收藏+]

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

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

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

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

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

1、可以使用Python脚本

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

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

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

 1 import arcpy
 2 #设置工作空间,即裁切参照shp文件所在的目录
 3 arcpy.env.workspace=rA:/1_5000/00/shp
 4 #获取工作空间中的Shapefile文件列表
 5 featureClasses=arcpy.ListFeatureClasses()
 6 #待裁切影像文件全路径,这里是.img格式
 7 inRasterPath=rA:/1_5000/GF2-2199394-20170221/GF2219939420170221F.img
 8 #裁切出的影像文件输出文件夹路径
 9 outFolderPath=rA:/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

 

ArcGis Python脚本——根据接图表批量裁切分幅影像

原文:https://www.cnblogs.com/yzhyingcool/p/10486352.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!