本文介紹基於Python中ArcPy模塊，對大量柵格遙感影像文件進行批量掩膜與批量重採樣的操作。 首先，我們來明確一下本文的具體需求。現有一個存儲有大量.tif格式遙感影像的文件夾；且其中除了.tif格式的遙感影像文件外，還具有其它格式的文件。 我們希望，依據一個已知的面要素矢量圖層文件，對上述文件 ...
本文介紹基於Python中ArcPy模塊，對大量柵格遙感影像文件進行批量掩膜與批量重採樣的操作。
首先，我們來明確一下本文的具體需求。現有一個存儲有大量.tif格式遙感影像的文件夾；且其中除了.tif格式的遙感影像文件外，還具有其它格式的文件。
我們希望，依據一個已知的面要素矢量圖層文件，對上述文件夾中的全部.tif格式遙感影像進行掩膜，並對掩膜後的遙感影像文件再分別加以批量重採樣，使得其空間解析度為1000 m。
明確了需求後，我們就可以開始具體的操作。首先，本文所需用到的代碼如下。
# -*- coding: utf-8 -*-
"""
Created on Fri Apr 15 16:44:26 2022
@author: fkxxgis
"""
import arcpy
from arcpy.sa import *
tif_file_path="E:/LST/Data/NDVI/03_Mosaic/"
shp_file="E:/LST/Data/Region/YellowRiver_nineprovince.shp"
out_file_path="E:/LST/Data/NDVI/04_Mask/"
resample_file_path="E:/LST/Data/NDVI/05_Resample/"
arcpy.env.workspace=tif_file_path
arcpy.env.extent=shp_file
tif_file_name=arcpy.ListRasters("*","tif")
for tif_file in tif_file_name:
mask_result=ExtractByMask(tif_file,shp_file)
mask_result_path=out_file_path+"/"+tif_file.strip(".tif")+"_Mask.tif"
mask_result.save(mask_result_path)
arcpy.env.workspace=out_file_path
tif_file_name=arcpy.ListRasters("*","tif")
for tif_file in tif_file_name:
resample_file_name=tif_file.strip(".tif")+"_Re.tif"
arcpy.Resample_management(tif_file,resample_file_path+resample_file_name,
1000,"BILINEAR")
其中，tif_file_path是原有掩膜前遙感圖像的保存路徑，shp_file是已知面要素矢量圖層文件的保存路徑，out_file_path是我們新生成的掩膜後遙感影像的保存路徑，resample_file_path則是最終重採樣後遙感影像的保存路徑。
在這裡，我們首先利用arcpy.ListRasters()函數，獲取路徑下原有的全部.tif格式的圖像文件，並存放於tif_file_name中；隨後，遍歷tif_file_path路徑下全部.tif格式圖像文件（即遍歷tif_file_name），並利用ExtractByMask()函數進行掩膜操作；其次，對於掩膜好的圖層，在其原有文件名後添加"_Mask.tif"尾碼，作為新文件的文件名。
對全部圖像文件完成掩膜操作後，我們繼續進行重採樣操作。和前述代碼思路類似，我們依然還是先遍歷文件，併在其原有文件名後添加"_Re.tif"尾碼，作為新文件的文件名；隨後，利用Resample_management()函數進行重採樣。其中，1000表示重採樣的空間解析度，在這裡單位為米；"BILINEAR"表示用雙線性插值的方法完成重採樣。
以上便是本次操作的全部代碼；我們這裡選擇在 IDLE (Python GUI) 中運行代碼。運行完畢，得到的一個結果文件如下圖；可以看到，遙感影像已經完成了掩膜，且空間解析度已經為1000 m。
至此，大功告成。
