ArcPy实现栅格图像的创建与遥感影像的批量拼接
本文介绍基于 Python 语言 arcpy 模块,实现栅格影像图层建立与多幅遥感影像数据批量拼接(Mosaic)的操作。
首先,相关操作所需具体代码如下: import os import arcpy file_path="G:/Postgraduate/LAI_Glass_RTlab/A2018161_Dif/DRT/" out_file_path="G:/Postgraduate/LAI_Glass_RTlab/A2018161_Dif/DRT/" out_file_name="Global.tif" file_name_list=os.listdir(file_path) tif_file_path=file_path+file_name_list[0] cell_size_x=arcpy.GetRasterProperties_management(tif_file_path,"CELLSIZEX") cell_size=cell_size_x.getOutput(0) value_type=arcpy.GetRasterProperties_management(tif_file_path,"VALUETYPE") describe=arcpy.Describe(tif_file_path) spatial_reference=describe.spatialReference arcpy.CreateRasterDataset_management(out_file_path,out_file_name,cell_size,"16_BIT_SIGNED", spatial_reference,"1") out_file=out_file_path+out_file_name for file in file_name_list: file_path_name=file_path+file print(file_path_name) arcpy.Mosaic_management([file_path_name],out_file)
其中, file_path 为存放有多景初始遥感影像的路径格式为.tif 栅格文件(如果不是.tif 格式,例如是.hdf 等文件,需首先进行文件格式的转换);out_file_path 为拼接后所得结果栅格图层的存放路径;out_file_name 为拼接后所得结果栅格图层的文件名称,其可选格式有很多,如下图所示。
在这里,我们默认所得拼接结果图层为一个(也就是file_path 文件夹中全部的待处理遥感影像最终全拼接在一起);如果大家需要使得拼接结果图层是多幅(也就是file_path 文件夹中待处理遥感影像依据区域、时间等分为很多不同的部分,每一部分拼接在一起),可以参考Python中GDAL读取处理栅格影像并结合质量评估QA波段筛选掩膜图像,利用其中的循环方式实现需求。
随后,通过os.listdir() 函数获取file_path 路径下的栅格文件,并存储于file_name_list 列表中。
接下来需要创建一个新的栅格图层。之所以要进行这一步骤,是因为本文后期选择用arcpy.Mosaic_management() 函数进行栅格的批量拼接,因此需要首先创建一个新的、空的栅格图层作为拼接的基准。如果大家的需求不是批量拼接栅格数据,而是单纯想利用arcpy 进行新栅格的创建,那就只看这一部分的代码即可。
在这里,我们选择用file_path 路径下的第一个栅格数据(下称"第一栅格")作为新栅格图层中各项属性(例如像素边长、像素数据格式等)的依据。首先,arcpy.GetRasterProperties_management() 函数获取第一栅格的像素x 边边长;因为一般栅格数据中像素都是正方形,因此我们就通过cell_size=cell_size_x.getOutput(0) 将第一栅格的像素x 边边长作为新栅格图层像素x 边与y 边二者的边长。再利用arcpy.GetRasterProperties_management() 函数获取第一栅格的数据格式;最后利用中间变量describe 获取第一栅格的空间参考信息。
完成以上步骤后,将已获取的第一栅格的各类信息通过函数arcpy.CreateRasterDataset_management() 带入新栅格中。在这里需要注意:尽可能在将要拼接时选择新栅格为"16_BIT_SIGNED" 及以下的数据格式(具体数据格式类别如下图),且将file_path 路径下待拼接的栅格数据的数据格式也全部修改为这一格式;否则可能会由于数据量大而导致拼接过程极慢。我之前就是由于选用了32 bit float 格式的栅格数据进行拼接,导致全球范围的MODIS一个植被产品数据拼接花了将近一天的时间。如果大家的栅格像素数据包含小数,可以通过乘上一个缩放系数的方式进行数据整数化。
代码最后的一个for 循环,就是遍历file_name_list 中的各个栅格数据,并通过arcpy.Mosaic_management() 函数加以拼接即可。
以上,便完成了本次批量拼接的操作。这里还有一点需要注意:由于arcpy 模块的限制,如果大家的Python版本是3.0 及以上,往往不能直接运行上述代码,最好是在ArcMap的Python运行框或其对应IDLE(如下图所示)中运行。
七夕了,怎么能把女友拍成女神今天终于是七夕了这么一个情侣狂欢的节日跟对象出去玩,拍照是必不可少的了之前也有很多朋友来问每次出去拍照,女友对照片都不满意那怎样才能把女友最美的一面展现出来呢首先,想好拍出美丽的照
不走出这三大误区,你永远是摄影新手!现在摄影已经成为了大众化的娱乐项目,门槛越来越低,加上手机拍照系统越来越好,很多人以为买了相机以后就会拍的比手机好,结果发现买了以后拍的也不是那么好,就放到一边吃灰去了。。按道理相
G29到底值不值得买?看看大家怎么说近年来,电竞行业通过互联网的传播,拥有了更多的熟知者和爱好者,各种类型的游戏层出不穷,手游,端游和模拟游戏的爱好者基数都在增长。其中有许多朋友就是赛车游戏的爱好者,无论是出于对儿时
G29神力科莎游戏体验上次和大家分享了G29在尘埃4里的使用体验,整体来说使用感是不错的,不知道在神力科莎里G29的表现会如何呢,我们一起来看看吧!神力科莎是一款模拟赛车游戏,较强的拟真性优秀的物理引擎
鼠标怎么选?适合自己的最重要二十世纪以来,无论是办公还是娱乐,人们都越来越离不开电脑了,自然而然的,也越来越离不开鼠标。生活中,许多人对外设产品不了解,所以总是买不到适合自己的产品,使用之后不是握着不舒服就是
纠结学习python还是c?看看这篇漫画吧而我们仅与Python打交道。简单来说Python可以自降身份,给其他语言开发的系统打个下手,提高他们的开发效率。我们把租来的车和司机看做一个整体,他们指代的是CC开发的系统,而我
计算机计数只要一推出那就已将落后了只要一推出那就已将落后了计算机计数,只要一推出那就已将落后了。学不会新技术或者效率不高,其实是很危险的。之前对于学习没有什么概念,就像是随意翻本书,看看书名,然后就随意的放在书架上
笑出腹肌!程序员从不撒谎,但注释却会大牛总说代码就是最好的注释。所以,世上本没有注释,代码写多了,也便有了注释。1这是一个被代码耽误的诗人2来一份1987年的代码看看3产品经理要对此负责4不敢看,也不敢问5Nike官
2021靳埭强设计奖报名开始!截稿至11月7日2021靳埭强设计奖如期而至!二十二年来,大赛一直秉承创立时的初衷,立足中国传统文化,探讨现代设计,希望借此推动中国艺术设计教育的发展,亦希望引导并鼓励参赛者不断进取,创新设计省思
中秋晚会丨中秋圆中秋缘团圆就现在中秋缘我们总说下次再约过几天再约中秋圆中秋缘团圆就!现!在!中秋晚会人间天上歌起舞乐飞旋为避免人员聚集奇酷教育以班级为单位举办中秋晚会歌王舞王齐聚她们的婀娜他们的动感天涯海角共赏表
剪不断理还乱英特尔十一代处理器广告惊现MacBook众所周知,苹果公司去年宣布了基于AppleSiliconM1处理器的Macbook和iMac电脑,并且确认会在两年内完成旗下所有电脑产品切换到AppleSilicon处理器,这意味