栅格矢量数据的裁剪

kelly1250230225

已于 2022-06-23 20:35:10 修改

阅读量979

点赞数

文章标签：前端大数据

于 2022-06-23 20:34:11 首次发布

本文链接：https://blog.csdn.net/kelly1250230225/article/details/125434456

版权

栅格裁剪（ExtractByMask）：同一栅格文件，不同矢量文件

#栅格裁剪（ExtractByMask）：同一栅格文件，不同矢量文件
#Author: Changqing Guo
# coding: utf-8
import arcpy
from  arcpy import env

from arcpy.sa import *

arcpy.CheckOutExtension("Spatial")
arcpy.env.workspace = "D:/Test/shp"
arcpy.env.overwriterOutput = True

outfiles = "D:/Test/out"
infiles = "D:/Test/Tif/00_dem.tif"
clipfiles = arcpy.ListFiles("*.shp")

for filename in clipfiles:
    print("Processing:" + filename)
    clipputfiles = arcpy.env.workspace + "/" + filename

    outputfiles = outfiles + "/" + filename[:-4]

    outExtractByMask = ExtractByMask(infiles, clipputfiles)
    outExtractByMask.save(outputfiles + ".tif")


print "***OVER***"
print arcpy.GetMessages()

栅格裁剪（ExtractByMask）：同一矢量文件，不同栅格文件

#栅格裁剪（ExtractByMask）：同一矢量文件，不同栅格文件
#Author: Changqing Guo
# coding: utf-8
import arcpy
from  arcpy import env

from arcpy.sa import *

arcpy.CheckOutExtension("Spatial")
arcpy.env.workspace = "Z:/0/01/01"
arcpy.env.overwriterOutput = True

outfiles = "Z:/0/01/05caijian/02et"
clipfiles = "Z://New_Shapefile.shp"
infiles = arcpy.ListFiles("*.tif")

for filename in infiles:
    print("Processing:" + filename)
    inputfiles = arcpy.env.workspace + "/" + filename

    outputfiles = outfiles + "/" + filename

    outExtractByMask = ExtractByMask(inputfiles, clipfiles)
    outExtractByMask.save(outputfiles + ".tif")

print "***OVER***"
print arcpy.GetMessages()