2015年10月9日写的,从QQ空间转过来。
裁切感兴趣区域外部影像
如图所示,1.shp为矢量图,1.tif为栅格图,我们尝试裁切矢量图的外面栅格区域。

步骤如下:
1.提取多边形所有节点的坐标。
用要素节点转点工具进行转换,提取出所有节点;
2.计算所有点的横纵坐标,用添加XY坐标工具;
3.导出坐标进excel中,用CONCATENATE命令按照代码格式编写文本;
arcpy.Point(40427531.4426999,4106332.48579999),
arcpy.Point(40427960.0685999,4106210.7773),
arcpy.Point(40427970.652, 4106067.902),
arcpy.Point(40427965.3602999,4105776.85969999),
arcpy.Point(40427912.4434999,4105618.1094),
arcpy.Point(40427579.0678, 4105766.27639999),
arcpy.Point(40427351.5257,4105988.52679999),
arcpy.Point(40427531.4426999,4106332.48579999),
4.在python中写入代码
>>> import arcpy
from arcpy import env
from arcpy.sa import *
polyPoints = [arcpy.Point(40427531.4426999,4106332.48579999),
arcpy.Point(40427960.0685999, 4106210.7773),
arcpy.Point(40427970.652, 4106067.902),
arcpy.Point(40427965.3602999, 4105776.85969999),
arcpy.Point(40427912.4434999, 4105618.1094),
arcpy.Point(40427579.0678, 4105766.27639999),
arcpy.Point(40427351.5257, 4105988.52679999),
arcpy.Point(40427531.4426999, 4106332.48579999)]
env.workspace = "D:\tttt" //工作目录
extPolygonOut = ExtractByPolygon("1.tif",polyPoints, "OUTSIDE") //裁切的影像
extPolygonOut.save("D:\tttt\out")
5.运行后,如下