BatchClipRasterByElement批量裁切影像

BatchClipRasterByElement批量裁切影像

上一篇 / 下一篇  2013-04-22 21:22:56 / 个人分类：Python

查看( 24 ) / 评论( 0 ) / 评分( 0 / 0 )

# 中文编码问题
# -*- coding: gbk -*-
#Author : RYQ
#CreateTime: 2013年4月18日
#工具说明：通过面状图层中得每个面状要素将栅格影像进行分割，并保存到GDB中。

print "正在导入所需模块，请稍侯........"
import arcpy
import os
from arcpy import env

#需要修改以下三个参数
#1.输入工作空间
env.workspace = r"G:\21Python\4.18\ENVI.gdb"
#2.输入需要裁剪得栅格影像名称
inRaster="浉河流域Raster"
#3.输入裁剪面状要素名称
inPolygonFeatures = "浉河流域1"

lyrFile = "Tmplyr"
arcpy.MakeFeatureLayer_management(inPolygonFeatures, lyrFile)
rows = arcpy.SearchCursor(inPolygonFeatures)
for row in rows:
    intValue = str(row.getValue("OBJECTID"))
    whereClause = "OBJECTID = "+ intValue
    arcpy.SelectLayerByAttribute_management(lyrFile,"NEW_SELECTION",whereClause)
    utRaster=inRaster + "r_" + intValue
    arcpy.Clip_management(inRaster,"#",outRaster,lyrFile,"#","ClippingGeometry")