??本文介紹基于Python中whitebox模塊，對大量長時間序列柵格遙感影像的每一個像元進行忽略NoData值的多時序平均值求取，
??在文章Python ArcPy批量計算多時相遙感影像的各像元平均值中，我們介紹了基于Python中Arcpy模塊實作多時相遙感影像資料的平均值求取方法，但是這一方法具有一個問題，即對于任意一個像元，只要該像元在任意一個時相的影像中是無效值（即為NoData），那么該像元在最終求出的平均值結果圖中像素值也將會是無效值NoData，這就導致在我們最終計算得到的平均值結果圖層中，具有很多空白區域（像素值為NoData的區域），
??為了解決這一問題，這里我們再介紹一種基于Python中另一個地理空間資料分析庫——whitebox，實作多時像遙感影像資料逐像元平均值的求取方法，
??首先，需要下載并安裝whitebox這一模塊，如果大家電腦中已經有了Anaconda環境，就可以直接按照Python地理分析庫whitebox在Anaconda中的配置這篇文章中介紹的方法下載、安裝whitebox，
??本文要實作的需求和文章Python ArcPy批量計算多時相遙感影像的各像元平均值中的一致，這里就不再贅述，本文所需用到的代碼如下，
# -*- coding: utf-8 -*-
"""
Created on Sun Apr 17 15:04:29 2022
@author: fkxxgis
"""
import glob
from whitebox import WhiteboxTools
tif_file_path="E:/LST/Data/MODIS/test/"
average_file_path="E:/LST/Data/MODIS/06_Average/"
wbt=WhiteboxTools()
wbt.work_dir=tif_file_path
tif_file_name=glob.glob(tif_file_path+"*.tif")
tif_file_year=tif_file_name[0][-18:-14]
one_year_tif_list=[]
for tif_file in tif_file_name:
if tif_file[-18:-14]==tif_file_year:
one_year_tif_list.append(tif_file)
tif_file_temp=tif_file
if tif_file==tif_file_name[len(tif_file_name)-1]:
wbt.average_overlay(inputs=';'.join(one_year_tif_list),
output=average_file_path+tif_file_year+"_Ave.tif")
else:
wbt.average_overlay(inputs=';'.join(one_year_tif_list),
output=average_file_path+tif_file_year+"_Ave.tif")
one_year_tif_list=[]
one_year_tif_list.append(tif_file)
tif_file_year=tif_file[-18:-14]
??其中，tif_file_path是原有計算平均值前遙感影像的保存路徑，average_file_path是我們新生成的求取平均值后遙感影像的保存路徑，也就是結果保存路徑，
??上述代碼的整體思路其實和文章Python ArcPy批量計算多時相遙感影像的各像元平均值這篇文章是非常類似的，首先，同樣需要在資源管理器中，將tif_file_path路徑下的各檔案以“名稱”排序的方式進行排序；隨后，利用arcpy.ListRasters()函式，獲取路徑下原有的全部.tif格式的影像檔案，并截取第一個檔案的部分檔案名，從而獲取其成像時間的具體年份，
??接下來，遍歷tif_file_path路徑下全部.tif格式影像檔案，其中，我們通過一個簡單的判斷陳述句if tif_file[0:4]==tif_file_year:，來確定某一年的遙感影像是否已經讀取完畢——如果已經讀取完畢，例如假如2001年成像的8幅遙感影像都已經遍歷過了，那么就對這8景遙感影像加以逐像元的平均值求取，并開始對下一個年份（即2005年）成像的遙感影像繼續加以計算；如果還沒有讀取完畢，例如假如2001年成像的8幅遙感影像目前僅遍歷到了第5幅，那么就不求平均值，繼續往下遍歷，直到遍歷完2001年成像的8幅遙感影像，
??這里相信大家也看到了為什么我們要在前期先將檔案夾中的檔案按照“名稱”排序——是為了保證同一年成像的所有遙感影像都排列在一起，遍歷時只要遇到一個新的年份，程式就知道上一個年份的所有影像都已經遍歷完畢了，就可以將上一個年份的所有柵格影像加以平均值求取，
??本文代碼和前期博客中代碼不一樣的部分就在于，這里是用到whitebox模塊而非arcpy模塊來實作同一年份遙感影像的逐像元平均值求取，在這里，wbt.average_overlay()函式就是我們實作這一步驟的關鍵，其中inputs引數表示需要進行平均值計算的同一年份的所有遙感影像，output表示求取平均值后得到的結果影像，
??最后，通過if tif_file==tif_file_name[len(tif_file_name)-1]:這個判斷，來確認是否目前已經遍歷到檔案夾中的最后一個影像檔案，如果是的話，就需要將當前成像年份的所有影像進行平均值的求取，并宣告代碼完成運行，
??這里需要注意，由于我們在此沒有用到arcpy模塊，因此代碼也就不一定非要在 IDLE (Python GUI) 中運行了，常見的編譯器都可以運行，在代碼運行程序中，還可以看到具體運行情況與進度，
??代碼運行完畢后，即可得到求解平均值后的結果圖層，
??最后還有一個問題——在我用這一代碼進行實踐后發現，如果計算平均值前的圖層具有兩個或兩個以上的波段，那么得到的結果圖層整體看還好，如下圖所示，
??但放大后，會發現得到的結果呈現出如下所示的條帶狀，
??而如果計算平均值前的圖層僅具有一個波段的話，就不會出現這種問題；如下圖所示，
??因此，大家在使用本文的代碼對大量長時間序列柵格遙感影像的每一個像元進行忽略Nodata值的多時序平均值求取時，一定注意輸入圖層要僅含有一個波段；否則結果就會出現條帶狀的錯誤，
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/552481.html
標籤：其他
上一篇：序列檢測器（兩種設計方法和四種檢測模式|verilog代碼|Testbench|仿真結果）
下一篇：返回列表