Python?ArcPy批量掩膜、重采樣大量遙感影像的操作

更新時間：2023年03月23日 08:28:44 作者：瘋狂學習GIS

這篇文章主要介紹了Python?ArcPy批量掩膜、重采樣大量遙感影像,本文介紹基于Python中ArcPy模塊，對大量柵格遙感影像文件進行批量掩膜與批量重采樣的操作，需要的朋友可以參考下

本文介紹基于Python中ArcPy模塊，對大量柵格遙感影像文件進行批量掩膜與批量重采樣的操作。

首先，我們來明確一下本文的具體需求?，F(xiàn)有一個存儲有大量.tif格式遙感影像的文件夾；且其中除了.tif格式的遙感影像文件外，還具有其它格式的文件。

我們希望，依據(jù)一個已知的面要素矢量圖層文件，對上述文件夾中的全部.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()函數(shù)，獲取路徑下原有的全部.tif格式的圖像文件，并存放于tif_file_name中；隨后，遍歷tif_file_path路徑下全部.tif格式圖像文件（即遍歷tif_file_name），并利用ExtractByMask()函數(shù)進行掩膜操作；其次，對于掩膜好的圖層，在其原有文件名后添加"_Mask.tif"后綴，作為新文件的文件名。

對全部圖像文件完成掩膜操作后，我們繼續(xù)進行重采樣操作。和前述代碼思路類似，我們依然還是先遍歷文件，并在其原有文件名后添加"_Re.tif"后綴，作為新文件的文件名；隨后，利用Resample_management()函數(shù)進行重采樣。其中，1000表示重采樣的空間分辨率，在這里單位為米；"BILINEAR"表示用雙線性插值的方法完成重采樣。

以上便是本次操作的全部代碼；我們這里選擇在 IDLE (Python GUI) 中運行代碼。運行完畢，得到的一個結(jié)果文件如下圖；可以看到，遙感影像已經(jīng)完成了掩膜，且空間分辨率已經(jīng)為1000 m。