Python多种方式实现批量裁剪

最新推荐文章于 2024-08-15 11:41:50 发布

xinxgcq

最新推荐文章于 2024-08-15 11:41:50 发布

阅读量1.4k

点赞数

分类专栏： Python 文章标签： Python 批量裁剪 Mask Clip

本文链接：https://blog.csdn.net/qq_38281942/article/details/107030606

版权

Python 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

Python多种方式实现批量裁剪

栅格裁剪（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()

栅格裁剪（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()

矢量裁剪（Clip_analysis）

#矢量裁剪（Clip_analysis）
#Author: Changqing Guo
# coding: utf-8
import arcpy
from  arcpy import env

from arcpy.sa import *

arcpy.env.workspace = "C:/Downloads"
arcpy.env.overwriterOutput = True

outfiles = "E:/01data/ruoergai"
infiles = "C:/Downloads/ld_dian/ld_dian.shp"
clipfiles = arcpy.ListFiles("*.shp")

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

    outputfiles = outfiles + "/" + filename

    arcpy.Clip_analysis(infiles,inclipfiles,outputfiles + ".shp","")

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