Python ArcPy批量掩膜/重采样大量遥感影像

科技资讯 投稿 11800 0 评论

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_manag

编程笔记 » Python ArcPy批量掩膜/重采样大量遥感影像

赞同 (54) or 分享 (0)
游客 发表我的评论   换个身份
取消评论

表情
(0)个小伙伴在吐槽