Arcpy批量点统计

最新推荐文章于 2024-07-12 16:16:27 发布

蒋户川柯基

最新推荐文章于 2024-07-12 16:16:27 发布

阅读量140

点赞数 1

文章标签： python Powered by 金山文档

本文链接：https://blog.csdn.net/m0_46663205/article/details/129747570

版权

选择模板栅格捕捉导出结果

# -*- coding: utf-8 -*-
#批量点统计

import arcpy
import os
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "D:/Desktop/diantongjishppro"
env.snapRaster = "D:/Desktop/soc060.tif"
outspace = "D:\\Desktop\\diantongjitiff\\"


list_shapefiles = arcpy.ListFiles("*shp")

for shp_name in list_shapefiles:

    print shp_name
    # Set local variables
    inPointFeatures = shp_name
    field = "count"
    cellSize = 1000
    neighborhood = NbrCircle(5000, "MAP")

    arcpy.CheckOutExtension("Spatial")

    outPointStatistics = PointStatistics(inPointFeatures, field, cellSize,neighborhood, "SUM")

    outPointStatistics.save(outspace+shp_name.split('.')[0]+'.tif')
    print ("         ---- OK！ ----         ")

print("   ---   All project is OK!   ---   ")