您的位置:首页 > 其它

提取矢量图层外部的栅格图

2017-10-09 18:56 393 查看
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.运行后,如下



内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: