Arcpy批量裁剪影像

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

贾村拓哉

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

阅读量152

点赞数

文章标签： arcgis

本文链接：https://blog.csdn.net/Jschao/article/details/129258015

版权

Arcpy批量裁剪影像，可以使用shp或者栅格影像裁剪

# -*- coding: utf-8 -*-
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "D:/LST_8day_eartdata/Night"  # 待裁剪tif影像存放路径
outputpath = "D:/LST_8day_eartdata/tif"  # 输出结果影像存放路径

# Set local variables
rasters = arcpy.ListRasters("*", "tif")
inMaskData = "D:/LST_8day_eartdata/mask.tif"  # shp或者栅格裁剪文件存放路径及名称

for raster in rasters:
    # 去掉后缀和文件名中的'.'
    i = 0
    file_name = raster.split('.')[i]
    i = i + 1
    print(file_name)
    # Check out the ArcGIS Spatial Analyst extension license
    arcpy.CheckOutExtension("Spatial")
    # Execute ExtractByMask
    outExtractByMask = ExtractByMask(raster, inMaskData)
    # Save the output
    outExtractByMask.save(outputpath + '/' + file_name + '_C.tif')

print('finished!')