Python?ArcPy批量掩膜、重采樣大量遙感影像的操作方法是什么

今天小編給大家分享一下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。

以上就是“Python ArcPy批量掩膜、重采樣大量遙感影像的操作方法是什么”這篇文章的所有內容，感謝各位的閱讀！相信大家閱讀完這篇文章都有很大的收獲，小編每天都會為大家更新不同的知識，如果還想學習更多的知識，請關注億速云行業資訊頻道。