在Zonal Statistics工具中使用重叠的区

 Howto: 在Zonal Statistics工具中使用重叠的区

文章编号： 39516

软件：  ArcGIS - ArcEditor 9.3, 9.3.1, 10 ArcGIS - ArcInfo 9.3, 9.3.1, 10 ArcGIS - ArcView 9.3, 9.3.1, 10

操作系统： N/A

摘要：

提供的脚本用于对相互重叠的区进行Zonal Statistics计算。

当使用多边形要素作为Zonal Statistics工具的区时，该工具会将多边形要素栅格化。如果多边形要素存在的重叠，栅格化会导致数据丢失。

Zonal Statistics中的这一缺陷，可以通过使用python逐一处理单个要素来解决。

内容：

# 导入系统模块

import arcgisscripting, os

# 创建Geoprocessor对象

gp = arcgisscripting.create(9.3)

gp.OverWriteOutput = 1

# 检查许可

gp.CheckOutExtension("Spatial")

#设置默认工作空间

gp.workspace = "C:/Data"

# 设置变量

raster = "measurements.tif"

inFC = "zone_features.shp"

outTable = "zonal_stats.dbf"

DBF_dir = "dbf" #folder for temp dbf outputs

if not os.path.exists(DBF_dir):

os.mkdir(DBF_dir)

#循环处理要素

inRows = gp.searchcursor(inputFC)

inRows.reset()

inRow = inRows.next()

while inRow:

# Zone_ID is the assumed name of the zone field

print "{0} is being processed".format(inRow.Zone_ID)

lyr = "Zone {0}".format(inRow.Zone_ID)

gp.MakeFeatureLayer_management(inFC, lyr, '"FID" = {0}'.format(inRow.FID))

tempTable = "{0}/DBF_{1}.dbf".format(DBF_dir, inRow.Test)

gp.ZonalStatisticsAsTable_sa(lyr, "Zone_ID", raster, tempTable, "DATA")

inRow = inRows.next()

#设置临时DBF文件的存储路径

gp.workspace = DBF_dir

TableList = gp.ListTables("*")

#合并临时DBF文件

gp.Merge(TableList, outTable)