Python：arcpy批处理

最新推荐文章于 2023-06-22 21:52:14 发布

Pria.GIS

最新推荐文章于 2023-06-22 21:52:14 发布

阅读量288

点赞数

分类专栏：自用文章标签： python

本文链接：https://blog.csdn.net/qq_39733892/article/details/124943841

版权

自用专栏收录该内容

2 篇文章 1 订阅

订阅专栏

这篇博客介绍了如何使用Python的arcpy库进行批量GIS数据处理，包括按单一矢量文件批量裁剪多个栅格数据、按多个矢量文件裁剪一个栅格数据以及将多个栅格图像拼接到一起。这些操作对于地理信息系统分析和地图制作非常实用。

摘要由CSDN通过智能技术生成

1. 批量裁剪（one shp）

按照一个矢量文件批量裁剪多个raster数据：

import os
import arcpy
import glob

filepath=r'F:\RasterData'
Clppath=r'F:\Clp_RasterData'
mask=r'F:\SHP\rec.shp'

#筛选tif文件
files=glob.glob(os.path.join(filepath,'*.tif'))
for f in files:
	outname=os.path.join(Clppath, os.path.basename(f).split(".")[0] + ".tif")
	out_extract=arcpy.sa.ExtractByMask(f, mask)
	out_extract.save(outname)
print('OVER')

2. 批量裁剪（multi-shp）

按照多个矢量文件裁剪一个raster数据：
主要区别在于文件的搜索

import os
import arcpy
import glob

filepath=r'F:\ratserdata\DEM.tif'
shppath=r'F:\SHP'
respath=r'F:\SplitDEM'

shps=glob.glob(os.path.join(shppath, "*.shp"))
for s in shps:
	outname = os.path.join(respath, os.path.basename(s).split(".")[0] + ".tif")
	out_extract = arcpy.sa.ExtractByMask(filepath, s)
	out_extract.save(outname)
print('OVER')

3.批量拼接

将所有图像都拼接到第一个图像中去：

import os
import arcpy
import glob

filepath=r'F:\ratserdata'
files=glob.glob(os.path.join(filepath, "*.tif"))
r0=files[0]

for f in files[1:]:
	# f向r0拼接，留在r0上
	#拼接方法还有：LAST，FIRST，MAXIMUM，MINIMUM
	arcpy.Mosaic_management(f,r0,'BLEND')
print('OVER')