快捷導(dǎo)航

Python?ArcPy實(shí)現(xiàn)批量拼接長(zhǎng)時(shí)間序列柵格圖像

更新時(shí)間：2023年03月16日 08:23:35 作者：瘋狂學(xué)習(xí)GIS

這篇文章主要介紹了如何基于Python中ArcPy模塊，對(duì)大量不同時(shí)相的柵格遙感影像按照其成像時(shí)間依次執(zhí)行批量拼接的方法，感興趣的可以了解一下

本文介紹基于Python中ArcPy模塊，對(duì)大量不同時(shí)相的柵格遙感影像按照其成像時(shí)間依次執(zhí)行批量拼接的方法。

在前期的文章Python arcpy創(chuàng)建柵格、批量拼接?xùn)鸥?/a>中，我們介紹了利用Python實(shí)現(xiàn)柵格遙感影像批量拼接的方法；但這篇文章實(shí)現(xiàn)的操作是將某個(gè)保存路徑下全部的柵格圖像文件加以拼接，換句話說，是對(duì)不同空間位置的同一時(shí)相的若干圖像加以拼接，拼接結(jié)果就只有一景大的圖像。而在實(shí)踐中，我們經(jīng)常還會(huì)需要對(duì)不同空間位置的不同時(shí)相的圖像分別加以拼接，拼接結(jié)果是很多景不同時(shí)相的大的圖像。那么，這種需求該怎么實(shí)現(xiàn)呢？

首先，我們來(lái)明確一下本文的具體需求?，F(xiàn)有一個(gè)存儲(chǔ)有大量.tif格式遙感影像的文件夾，其中每一個(gè)遙感影像的文件名中都包含有該圖像的成像時(shí)間，如下圖所示。

我們希望，對(duì)于同一天成像的遙感影像進(jìn)行拼接——例如，上圖中具有2001年第185天成像的遙感影像10幅，每一幅都是這一天在不同空間位置的成像；同時(shí)有2001年第193天成像的遙感影像10幅。我們希望首先將第185天成像的10幅遙感影像加以拼接，隨后再對(duì)第193天成像的10幅遙感影像加以拼接，以此類推。在遙感影像整體數(shù)量較少時(shí)，我們或許還可以逐一手動(dòng)拼接；而當(dāng)圖像數(shù)量很多時(shí)，就需要借助代碼來(lái)實(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格式的圖像文件，并截取第一個(gè)文件的部分文件名，從而獲取其成像時(shí)間；接下來(lái)，做好創(chuàng)建一個(gè)新的柵格文件的準(zhǔn)備，這一部分代碼的含義在本文開頭提及的那一篇文章 Python arcpy創(chuàng)建柵格、批量拼接?xùn)鸥?/a>中已有提及，這里就不再贅述。

接下來(lái)，遍歷tif_file_path路徑下全部.tif格式圖像文件。其中，我們通過一個(gè)簡(jiǎn)單的判斷語(yǔ)句，來(lái)確定某一成像時(shí)間的遙感影像是否已經(jīng)讀取完畢——如果已經(jīng)讀取完畢，例如假如第185天成像的10幅遙感影像都已經(jīng)遍歷過了，那么就對(duì)這十景遙感影像加以拼接；如果還沒有讀取完畢，例如假如第185天成像的10幅遙感影像目前僅遍歷到了第8幅，那么就不拼接，繼續(xù)往下遍歷。

這里相信大家也看到了為什么我們要在前期先將文件夾中的文件按照“名稱”排序——是為了保證同一成像時(shí)間的所有遙感影像都排列在一起，遍歷時(shí)只要遇到一個(gè)新的成像時(shí)間，程序就知道上一個(gè)成像時(shí)間的所有圖像都已經(jīng)遍歷完畢了，就可以將上一個(gè)成像時(shí)間的所有柵格圖像加以拼接。

最后，通過tif_file==tif_file_name[len(tif_file_name)-1]這個(gè)判斷，來(lái)確認(rèn)是否目前已經(jīng)遍歷到文件夾中的最后一個(gè)圖像文件。如果是的話，就需要將當(dāng)前成像時(shí)間的所有圖像進(jìn)行拼接，并完成代碼的運(yùn)行。

在 IDLE (Python GUI) 中運(yùn)行代碼。代碼運(yùn)行完畢后，我們開看一下結(jié)果文件夾?？梢钥吹?，其中的圖像已經(jīng)是按照成像時(shí)間，分別完成拼接后的結(jié)果了。

至此，大功告成。

到此這篇關(guān)于Python ArcPy實(shí)現(xiàn)批量拼接長(zhǎng)時(shí)間序列柵格圖像的文章就介紹到這了,更多相關(guān)Python ArcPy拼接?xùn)鸥駡D像內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！