再議使用Python批量裁切柵格

標籤：

曾經寫過《使用Python指令碼批量裁切柵格》，但今天又遇到這個情況則發現了問題。我們遇到的實際問題往往是有一個需要裁剪的影像（大塊的），另外有一個向量面，現在需要按向量面每一個要素進行裁剪，無奈arcgis裡的工具無法方便地做到。只能自己寫工具，這次使用了clip而不是ExtractByMask，因為ExtractByMask有很多限制！

下面是工具的操作樣本：按每一個要素進行裁剪柵格，輸出柵格以選擇的欄位命名，前提是欄位的每個值是唯一的。

下面是訊息輸入和裁剪向量表的屬性工作表：

下面是Python原始碼

# ---------------------------------------------------------------------------# Purpose : ClipRasterByFeature# Author ：gisweis# Date ：2015.7.21# Version : ArcGIS 10.1# Email :[email protected]# Notes :# --------------------------------------------------------------------------- import sysreload(sys)sys.setdefaultencoding( "utf-8" ) import arcpyimport string try:raster = arcpy.GetParameterAsText(0) #clip rasterclip_feat = arcpy.GetParameterAsText(1) #clip featureclassfield = arcpy.GetParameterAsText(2) #name fieldoutworkspace = arcpy.GetParameterAsText(3) #output wsouttype = arcpy.GetParameterAsText(4) #output ws total = int(arcpy.GetCount_management(clip_feat).getOutput(0))count= 1for row in arcpy.SearchCursor(clip_feat):mask=row.getValue("Shape")extent=str(mask.extent.XMin)+" " +str(mask.extent.YMin)+" " +str(mask.extent.XMax)+" " +str(mask.extent.YMax)outPath=outworkspace+"\\"+str(row.getValue(field)+outtype)arcpy.AddMessage("chipping: " + str(row.getValue(field)) + "...count:"+str(total)+"\\"+str(count))arcpy.Clip_management(raster,extent,outPath,mask,"0","ClippingGeometry")count=count+1except arcpy.ExecuteError:    print arcpy.GetMessages()

　　

 

再議使用Python批量裁切柵格