python批量裁剪矢量数据_ARCGIS\python批量裁剪栅格数据

#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’

转载自：/weixin_43187378/article/details/88687816