Arcpy专题：批量裁剪影像

最新推荐文章于 2024-06-18 16:30:18 发布

敏敏先生

最新推荐文章于 2024-06-18 16:30:18 发布

阅读量584

点赞数

分类专栏： ArcGIS arcpy 遥感文章标签：深度学习人工智能

本文链接：https://blog.youkuaiyun.com/weixin_45929124/article/details/124783601

版权

ArcGIS 同时被 3 个专栏收录

4 篇文章

订阅专栏

遥感

3 篇文章

订阅专栏

arcpy

1 篇文章

订阅专栏

这是一个使用Python和arcpy库批量裁剪遥感影像的脚本。脚本根据指定的shapefile文件，对指定文件夹下的影像数据进行裁剪，并将结果保存到输出文件夹中。主要涉及的空间分析操作是Extract By Mask。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

# -*- coding: utf-8 -*-
# @Time : 2020/12/10 19:16
# @Author : ZhangMin
import arcpy
import glob
import os

'''
批量裁剪：
    对输入文件夹下的数据，根据shp文件，进行批量裁剪，结果存放于输出文件夹
    
    需要修改：
       

'''
arcpy.CheckOutExtension('Spatial')

# 输出文件夹
outfile = r"H:\大论文\福建省\2008-2012Landsat5&7\30m分辨率修正版、\按潮区裁剪的影像".decode('utf-8')
# shp文件夹
shpmask = r"H:\大论文\福建省\1998-2002Landsat4&5\30m修正版\新潮区划分".decode('utf-8')
# 影像路径
rasters = r'H:\大论文\福建省\2008-2012Landsat5&7\30m分辨率修正版、\潮滩频率\20082012潮滩范围频率.tif'.decode('utf-8')

files = []

for shp in os.listdir(shpmask):
    if os.path.splitext(shp)[1] == '.shp':
        sourcefile = os.path.join(shpmask, shp)
        files.append(sourcefile)
        # print(sourcefile)

        # 截断文件名
        (filepath, tempfilename) = os.path.split(sourcefile)
        (filename, extension) = os.path.splitext(tempfilename)
        print(filename)
        # 输出路径+文件名
        outname = os.path.join(outfile, filename + ".tif")
        print(outname)

        # 裁剪
        out_extract = arcpy.sa.ExtractByMask(rasters, sourcefile)
        out_extract.save(outname)