使用Python脚本批量裁切栅格

对栅格的裁切，我们通常使用裁切（数据管理-栅格-栅格处理）或按掩膜提取（空间分析-提取分析）来裁切，裁切的矢量要素通常是一个要素图层或Shape文件。如果要进行批量处理，可以使用ToolBox中的批量处理工具。但是，有时我们选择的裁切矢量要素不是一个要素图层，而是要素图层中的一个要素，如按标准图幅裁切栅格影像，以前的做法是将要素一个一个的导出来，再进行批量处理。快速导出方法：ArcGIS按字段属性分割文件。

现在，可以使用Python脚本来批量裁切，在网上找了位大神的原始代码学习一下：

#-*- encoding:UTF-8 -*-# Tool Name: 遍历要素裁剪栅格
# Source Name: extractByMask.py
# Version: ArcGIS 10.1
# Author: GIScloud
# Created: 2012/11/1import 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()

将脚本导入ToolBox中（导入方法：ArcGIS使用Python脚本工具），执行报错：UnicodeEncodeError: 'ascii' codec can't encode characters in position 32-34: ordinal not in range(128)问题，使用的是ulipad编译器。

看样子是编码的问题，查了一下解决方案：在开头加上代码设置默认编辑

import sys
reload(sys)
sys.setdefaultencoding( "utf-8" )

能够实现部分裁切了，但最后一个出了问题，看消息输出的文件名中怎么多了一个句号？

查看属性表，没有问题。将输出字段"北川羌族自治县"改为"北川羌族自治"能够正确输出，如果改成"北川羌族自治北川羌族自治"，又报相同的错误。说明应该就是输出名称长度的问题。

还有一个问题，输出的不同栅格边界接连有点小问题，如果避免这种情况，只能更改系统变量，输出更小的单元格大小。

参考：http://bbs.esrichina-bj.cn/ESRI/viewthread.php?tid=121642

http://blog.sina.com.cn/s/blog_64a3795a01018vyp.html

转载于:https://www.cnblogs.com/liweis/p/4252097.html

使用Python脚本批量裁切栅格相关推荐

python脚本批量生成数据
在平时的工作中,经常会遇到造数据,特别是性能测试的时候更是需要大量的数据.如果一条条的插入数据库或者一条条的创建数据,效率未免有点低.如何快速的造大量的测试数据呢?在不熟悉存储过程的情况下,今天给大家 ...
Python脚本批量重命名图片
Python脚本批量重命名图片 # -*- coding: utf-8 -*- """ Created on Tue Mar 23 13:34:59 2021@autho ...
使用python脚本批量修改vc工程文件
作者:朱金灿来源:clever101的专栏为什么大多数人学不会人工智能编程?>>> 问题描述 vc工程文件中的很多属性设置,比如输出目录.中间目录和输出文件名等等.使用vs ...
arcgis使用python脚本批量裁剪影像
功能描述: 对于含有多个要素的矢量文件shp.栅格影像raster,按照shp中的各要素范围,逐个对raster进行提取,并保存到文件夹中效果如图所示: 主要思路: 1)获取矢量图层.栅格图层 2) ...
python脚本批量生成雷达回波数据集
批量生成雷达回波数据集目前,没有找到探地雷达回波图的开源数据集,鉴于实验成本限制,采用gprMax v3.1.5仿真生成. 这里参照两个博主: 链接1:https://blog.csdn.net/g ...
在blender中使用python脚本批量复制平移生成模型
本案例需求: 从基本的建筑单元按照字形平面布局生成综合建筑体. 先在blender中用手工制作好一个建筑单元(名称定为"cube"),然后在blender中打开一个"Te ...
python编写arcgis脚本_ArcGis Python脚本——批量添加字段
插个广告,制作ArcGIS的Tool工具学习下面的教程就对了: 以下是正文: 先看如何增加一个字段函数:arcpy.AddField_management 语法:AddFields_manageme ...
python脚本批量登录crt_Python实现批量新建SecureCRT Session
最近因为工作需要,我需要在ssh的时候保存几千台网关的session,工作量相当大(也就是ssh的时候需要记住用户名和密码,然后还要再session选项中录入enable密码,相当繁琐),而且设备的用 ...
python脚本--批量下载GNSS数据
文章目录 1. 前言 2. 思路 3. Talk is cheap 4. 总结 1. 前言应某一瓜皮要求,为其写一个用于批量下载GNSS数据的脚本.恰好前几天研究过如何在教务系统中自动签到,便有了这 ...

使用Python脚本批量裁切栅格

使用Python脚本批量裁切栅格相关推荐

最新文章

热门文章