利用ArcGIS Python批量拼接遥感影像(arcpy batch processing)

 本篇文章将说明如何利用ArcGIS 10.1自带的Python IDLE进行遥感影像的批量拼接与裁剪。

1.运行环境：ArcGIS10.1 (安装传送门)、Python IDLE

2.数据来源：地理空间数据云 GDEMV2 30M分辨率数字高程数据

3.解决问题：制作山西省的DEM影像

如下图所示，以30M分辨率数字高程数据为例，影像皆是固定范围的经纬度保存在其服务器上，外在表现以小幅正方形影像。如果手动进行拼接，工作量会非常大且容易出错。
 

1.拼接影像

开始菜单-ArcGIS-IDLE（Python GUI） 打开IDLE。File-New Window, 文件以 .py结尾，如MosaicToNewRasters.py

import arcpy
import os
 
#指定工作目录，即存放影像的目录
arcpy.env.workspace = r"E:\arcpyData\shanxi\N34_N35"
 
#指定该工作空间下的一副影像为基础影像，为后面的参数提取做准备
base = "ASTGTM2_N34E109_dem.tif"
 
#以下一段代码是为执行拼接做参数准备
out_coor_system = arcpy.Describe(base).spatialReference #获取坐标系统
dataType = arcpy.Describe(base).DataType 
piexl_type = arcpy.Describe(base).pixelType 
cellwidth = arcpy.Describe(base).meanCellWidth #获取栅格单元的的宽度
bandcount = arcpy.Describe(base).bandCount #获取bandCount
 
#打印一些信息
print out_coor_system.name
print dataType
print piexl_type
print cellwidth
print bandcount
 
arcpy.CheckOutExtension("Spatial")
 
#提取待拼接影像的文件名，且中间以;隔开，例如：a.tif;b.tif;c.tif
rasters = []
for ras in arcpy.ListRasters("*dem.tif"):    #for循环，对wrokspace下的所有以dem.tif结尾的影像进行过滤
    rasters.append(ras)
 
ras_list = ";".join(rasters)                 #字符串拼接
 
#打印出来，看看什么结果吧
print ras_list
 
 
#指定输出文件夹
outFolder = r"E:\arcpyData\shanxi\N34_N35\n34_n35_out"
 
#执行拼接操作
arcpy.MosaicToNewRaster_management(ras_list, outFolder, "shanxi_n34_n35_dem_data.tif", out_coor_system, "16_BIT_SIGNED", cellwidth, bandcount, "LAST", "FIRST")
 

《遥感影像数据》共享方法如下：

（1）人员，限定为关注小编的用户。

（2）各类项目（包括各类科研项目）申请本数据扔享受免费政策，但需向本号捐赠一定数量的硬盘才能获取。

（3）捐赠硬盘可免留言获取数据。