Python 批量 Arcgis 通过矢量数据裁剪栅格数据

    网上已经有批量裁剪的源代码,但是运行问题有很多,经过改动,成功运行,希望对他人有所启发。如果有问题,欢迎留言,希望能够回复帮忙。以下是从网上找到的源码

#-*- encoding:UTF-8 -*-

  

# Tool Name: 遍历要素裁剪栅格

# Source Name: extractByMask.py

# Version: ArcGIS 10.1

# Author: GIScloud

# Created: 2012/11/1

  

import arcpy

import string

from arcpy.sa import *

  

try:

    raster = arcpy.GetParameterAsText(0) #要裁剪的栅格

    clip_feat = arcpy.GetParameterAsText(1) #裁剪要素类

    field = arcpy.GetParameterAsText(2) #命名字段

    outworkspace = arcpy.GetParameterAsText(3) #命名字段裁剪后输出目录

  

    for row in arcpy.SearchCursor(clip_feat):

        mask=row.getValue("Shape")

        outPath=outworkspace+"\\"+str(row.getValue(field))

        outExtractByMask = ExtractByMask(raster,mask)

        outExtractByMask.save(outPath)

  

except arcpy.ExecuteError:

    print arcpy.GetMessages()

此为源码,资源来自互联网,链接为:https://www.cnblogs.com/liweis/p/4252097

运行过程中发现有很多问题,总结如下:

1)报错0100240,不能将运行好的文件保存到某路径。解决办法,更改environment中的workspace。此外,看看你的属性字段是不是有空格,字符是不是超过了13个,是不是有特殊符号,例如’之类,都会引起报错,如果这些都检查没有问题,再看看导入到arcgis工具箱中的脚本是不是有缩进的问题,因为python是由缩进控制作用域,很可能引起代码不循环,不判断等等。更改后代码如下:

#-*- encoding:UTF-8 -*-
  
# Tool Name: maskproperty
# Source Name: extractByMask.py
# Version: ArcGIS 10.1
# Author: GIScloud
# Created: 2012/11/1
  
import sys
reload(sys)
sys.setdefaultencoding( "utf-8" )
import arcpy
from arcpy import env
import string
from arcpy.sa import *

arcpy.CheckOutExtension("spatial")
env.workspace = r"G:\image_bach\origal_data"#enveriment setting
  
try:
   raster = arcpy.GetParameterAsText(0)                                       #input roster
   clip_feat = arcpy.GetParameterAsText(1)                                  #input vector
   field = arcpy.GetParameterAsText(2)                                         #file 字段
   outworkspace = arcpy.GetParameterAsText(3)                         #outputdectory
  
   for row in arcpy.SearchCursor(clip_feat):
      mask=row.getValue("Shape")
      outpart=str(row.getValue(field)) 
      n=len(outpart)
      if n <= 13:
        outPath=outworkspace+"\\"+''.join(outpart.split())       #由于python的文件命名不能超过13个字符,判断是否超了
      else:
        noblankstr=''.join(outpart.split())                                  #由于命名字段可能存在空格,去掉空格
        outPath=outworkspace+"\\"+noblankstr[0:13]             #如果字段的字符超过13个,截取前13个
      outExtractByMask = ExtractByMask(raster,mask)
      outExtractByMask.save(outPath)
except arcpy.ExecuteError:
   print arcpy.GetMessages()

    运行时候是将此工具配置到arcgis中自己新建的script中运行的。

运行结果:

 此为运行结果,文件名是属性字段命名的,如果有空格,将空格去掉作为文件名。

    生成文件夹,以属性字段命名,生成的格式是arcgis的栅格文件。

  此为原创,转载请注明出处。 

转载自:https://blog.csdn.net/qq_42237611/article/details/81204832

You may also like...