#coding:utf-8

裁剪太行山区域的降水量和气象数据

import sys
import os, string
import arcpy
import arcpy.sa
import pdb
from arcpy import env
from arcpy.sa import *

env.overwriteOutput=True #可以覆盖
env.workspace=r’D:/cp45_pt/’# 读取的栅格数据集所在文件夹
rasters=arcpy.ListRasters()
dsc=arcpy.Describe(r’D:/CANKAO/dem1.tif’) #找一个坐标系
dsc.spatialReference
print rasters
mask_fea=‘D:/cankao/taihangshan.shp’
aa=arcpy.CheckOutExtension(“Spatial”)
print 1.0
for ra in rasters:
outRas=“D:/OutPut/pt/”+ra+“clip.tif”
arcpy.DefineProjection_management(ra,dsc.spatialReference)
# 给定义一个坐标系，使用常用的WGS_1984坐标系
ores=arcpy.sa.ExtractByMask(ra,mask_fea)
ores.save(outRas)
print ra
print ‘work have done’

转载自：https://blog.csdn.net/weixin_43187378/article/details/88687816