TIFF文件裁剪

python2.7环境
话不多说,直接上代码

import arcpy
import os

arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput = 1
arcpy.env.overwriteOutput=True
arcpy.env.workspace = "G:\\RF\\01_result\\tif" #shp文件所在地址
shps = arcpy.ListFiles("*.shp")

arcpy.env.workspace = "G:\\RF\\01_result\\tif" #tif文件所在地址
rasters = arcpy.ListRasters("*", "tiff")

for shp in shps:

	mask = shp

	print(mask)

	catchment_name = 'china'
	#

	for raster in rasters:
		print (raster)

		if os.path.exists(r"G:\\RF\\01_result\\01_hu_result_tif"):#输出文件所在地址
			out = (r"G:\\RF\\01_result\\01_hu_result_tif\\" + os.path.splitext(raster)[0] + ".tif")
		else:
			os.makedirs(r"G:\\RF\\01_result\\01_hu_result_tif")#输出文件所在地址
			out = (r"G:\\RF\\01_result\\01_hu_result_tif\\" + os.path.splitext(raster)[0] + ".tif")

		arcpy.gp.ExtractByMask_sa(raster, mask, out)

		print("clip_"+raster+"  has done")
	print(catchment_name+"_Q_has done")

print("All done")


`
## tif转csv,CSV转tif,csv文件操作,tif文件操作,nc文件操作,机器学习XGBOOST,RF等看作者名字


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值