利用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 bandcountarcpy.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")

博客来源引用：地理遥感生态网www.gisrs.cn

同时，地理遥感生态网www.gisrs.cn上分享了很多地理遥感领域的科学数据（土地利用数据、npp净初级生产力数据数据、NDVI数据、气象数据（降雨量、气温、蒸散量、辐射、湿度、日照时数、风速、水汽压数据）、径流量数据、夜间灯光数据、统计年鉴、道路网、POI兴趣点数据、GDP分布、人口密度分布、三级流域矢量边界、地质灾害分布数据、土壤类型、土壤质地、土壤有机质、土壤PH值、土壤质地、土壤侵蚀、植被类型、自然保护区分布、建筑轮廓分布等等地理数据，以及关于gis、遥感从方面的操作教程）。

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

利用ArcGIS Python批量拼接裁剪遥感影像(arcpy batch processing)
本篇文章将说明如何利用ArcGIS 10.1自带的Python IDLE进行遥感影像的批量拼接与裁剪. 1.运行环境:ArcGIS10.1 (安装传送门).Python IDLE 2.数据来源:地理空 ...
使用PhotoShop拼接遥感影像
如何使用PhotoShop拼接遥感影像,并进行调色处理只是记录一下流程,不会贴详细的步骤首先要知道的重要的步骤,感觉还是很繁琐只是记录一下流程,不会贴详细的步骤首先要知道的因为大众版的ps ...
python/gdal处理遥感影像（读取、投影转换、裁剪、建立图像金字塔等）
python/gdal处理遥感影像(读取.投影转换.裁剪.建立图像金字塔等) gdal库简单介绍 python使用gdal 一.安装python环境二.安装gdal库三.使用gdal处理遥感影像 ...
基于Python(Arcpy)批量对遥感影像MK与SEN检验，结果为两张TIF图（MK与SEN）
def mkSen_cal(source_folder,outpt_folder):arcpy.CheckOutExtension("spatial")arcpy.env.work ...
Python批量实现栅格影像可视化
任务描述: 对栅格影像进行可视化并保存图片数据准备: 文件夹(工作空间)内含有需要批量可视化的tif影像工具准备: Python环境操作: 替换代码中的路径.根据需要调整可视化图片格式 # -* ...
Arcgis中怎么进行遥感影像的标准分幅
第一步:91地图上下载图幅框第二步:KML转图层--转换工具. 第三步:选择自己研究区域的图框--按位置选择. 第四步:分割图框--分割工具. 注:id为各小图框的标准图幅编号.如下所示: 第五步: ...
使用Python批量拼接图片
前言当需要将多张图像拼接成一张更大的图像时,通常会用到图片拼接技术.这种技术在许多领域中都有广泛的应用,例如计算机视觉.图像处理.卫星图像.地理信息系统等等.在实际应用中,拼接图像可以用于创建全景图 ...
python可视化DEM遥感影像（tif格式）||xarray使用
1.利用xarray导入tif格式的DEM影像,并让其可视化. 参考博客:主要文章和博客1及博客2. 代码如下: 首先导入相关的包: import numpy as np import matplot ...
python批量拼接两个文件夹相同名字的图片
结果如上图所示: 为了直观的观察图片效果:可以将图片拼接在一起:代码如下 import cv2 import numpy as np import pandas as pd import os, ha ...
用Python批量下载DACC的MODIS数据
本人初次尝试用Python批量下载DACC的MODIS数据,记下步骤,提醒自己,数据还在下载,成功是否未知,等待结果中...... 若有大佬发现步骤有不对之处,望指出,不胜感激. 1.下载Python ...

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

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

最新文章

热门文章