Python?ArcPy實(shí)現(xiàn)批量拼接長時間序列柵格圖像
本文介紹基于Python中ArcPy模塊,對大量不同時相的柵格遙感影像按照其成像時間依次執(zhí)行批量拼接的方法。
首先,我們來明確一下本文的具體需求?,F(xiàn)有一個存儲有大量.tif
格式遙感影像的文件夾,其中每一個遙感影像的文件名中都包含有該圖像的成像時間,如下圖所示。
我們希望,對于同一天成像的遙感影像進(jìn)行拼接——例如,上圖中具有2001
年第185
天成像的遙感影像10
幅,每一幅都是這一天在不同空間位置的成像;同時有2001
年第193
天成像的遙感影像10
幅。我們希望首先將第185
天成像的10
幅遙感影像加以拼接,隨后再對第193
天成像的10
幅遙感影像加以拼接,以此類推。在遙感影像整體數(shù)量較少時,我們或許還可以逐一手動拼接;而當(dāng)圖像數(shù)量很多時,就需要借助代碼來實(shí)現(xiàn)了。
明確了需求后,我們就可以開始具體的操作。首先,本文所需用到的代碼如下。
# -*- coding: utf-8 -*- """ Created on Fri Apr 15 13:21:55 2022 @author: fkxxgis """ import os import arcpy tif_file_path="E:/LST/Data/NDVI/02_TIFF/" out_file_path="E:/LST/Data/NDVI/03_Mosaic/" arcpy.env.workspace=tif_file_path tif_file_name=arcpy.ListRasters("*","tif") tif_file_date=tif_file_name[0][1:8] one_day_tif_list=[] tif_file_example_path=tif_file_path+tif_file_name[0] cell_size_x=arcpy.GetRasterProperties_management(tif_file_example_path,"CELLSIZEX") cell_size=cell_size_x.getOutput(0) value_type=arcpy.GetRasterProperties_management(tif_file_example_path,"VALUETYPE") describe=arcpy.Describe(tif_file_example_path) spatial_reference=describe.spatialReference for tif_file in tif_file_name: if tif_file[1:8]==tif_file_date: one_day_tif_list.append(tif_file) tif_file_temp=tif_file if tif_file==tif_file_name[len(tif_file_name)-1]: out_file_name=tif_file[1:8]+".tif" 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 tif_file_new in one_day_tif_list: arcpy.Mosaic_management([tif_file_path+tif_file_new],out_file) else: out_file_name=tif_file_temp[1:8]+".tif" 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 tif_file_new in one_day_tif_list: arcpy.Mosaic_management([tif_file_path+tif_file_new],out_file) one_day_tif_list=[] one_day_tif_list.append(tif_file) tif_file_date=tif_file[1:8]
其中,tif_file_path
是原有拼接前遙感圖像的保存路徑,out_file_path
是我們新生成的拼接后遙感影像的保存路徑。
在這里,我們需要首先在資源管理器中,將tif_file_path
路徑下的各文件以“名稱”排序的方式進(jìn)行排序;隨后,利用arcpy.ListRasters()
函數(shù),獲取路徑下原有的全部.tif
格式的圖像文件,并截取第一個文件的部分文件名,從而獲取其成像時間;接下來,做好創(chuàng)建一個新的柵格文件的準(zhǔn)備,這一部分代碼的含義在本文開頭提及的那一篇文章Python arcpy創(chuàng)建柵格、批量拼接?xùn)鸥?/a>中已有提及,這里就不再贅述。
接下來,遍歷tif_file_path
路徑下全部.tif
格式圖像文件。其中,我們通過一個簡單的判斷語句,來確定某一成像時間的遙感影像是否已經(jīng)讀取完畢——如果已經(jīng)讀取完畢,例如假如第185
天成像的10幅遙感影像都已經(jīng)遍歷過了,那么就對這十景遙感影像加以拼接;如果還沒有讀取完畢,例如假如第185
天成像的10幅遙感影像目前僅遍歷到了第8幅,那么就不拼接,繼續(xù)往下遍歷。
這里相信大家也看到了為什么我們要在前期先將文件夾中的文件按照“名稱”排序——是為了保證同一成像時間的所有遙感影像都排列在一起,遍歷時只要遇到一個新的成像時間,程序就知道上一個成像時間的所有圖像都已經(jīng)遍歷完畢了,就可以將上一個成像時間的所有柵格圖像加以拼接。
最后,通過tif_file==tif_file_name[len(tif_file_name)-1]
這個判斷,來確認(rèn)是否目前已經(jīng)遍歷到文件夾中的最后一個圖像文件。如果是的話,就需要將當(dāng)前成像時間的所有圖像進(jìn)行拼接,并完成代碼的運(yùn)行。
在 IDLE (Python GUI) 中運(yùn)行代碼。代碼運(yùn)行完畢后,我們開看一下結(jié)果文件夾??梢钥吹剑渲械膱D像已經(jīng)是按照成像時間,分別完成拼接后的結(jié)果了。
至此,大功告成。
到此這篇關(guān)于Python ArcPy實(shí)現(xiàn)批量拼接長時間序列柵格圖像的文章就介紹到這了,更多相關(guān)Python ArcPy拼接?xùn)鸥駡D像內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
使用python BeautifulSoup庫抓取58手機(jī)維修信息
這篇文章主要介紹了一個使用python抓取58手機(jī)的精準(zhǔn)商家信息,使用BeautifulSoup API的方法2013-11-11python學(xué)習(xí)之SpaCy庫的高級特性詳解
在之前的文章中,我們介紹了SpaCy庫的一些基本概念和功能,在這篇文章中,我們將深入學(xué)習(xí)一些更高級的特性,包括詞向量、依賴性解析、和自定義組件2023-07-07

python 讀入多行數(shù)據(jù)的實(shí)例

python爬蟲入門教程之點(diǎn)點(diǎn)美女圖片爬蟲代碼分享

python矩陣轉(zhuǎn)換為一維數(shù)組的實(shí)例

Python創(chuàng)建對稱矩陣的方法示例【基于numpy模塊】

python基于socket進(jìn)行端口轉(zhuǎn)發(fā)實(shí)現(xiàn)后門隱藏的示例