作者:不言不语都可以温柔_631 | 来源:互联网 | 2023-09-08 09:06
ArcGIS自带的ZonalStatisticsAsTable在分区统计栅格面积时很好用,但是这个工具只能针对栅格,不能针对矢量。最近有这个需求,所以就做了一个,基本都是ArcPy
ArcGIS 自带的 Zonal Statistics As Table 在分区统计栅格面积时很好用,但是这个工具只能针对栅格,不能针对矢量。
最近有这个需求,所以就做了一个,基本都是 ArcPy 的自带函数,输入的参数和 Zonal Statistics As Table 高度类似,在此不作详细解释,请君查看代码。
下附:
import arcpy,os,pandas,shutil
def write2DListToCsv(listItem,fileName):
fileObject = open(fileName, 'w')
for i in range(len(listItem)):
for j in range(len(listItem[i])):
fileObject.write(str(listItem[i][j]))
if j fileObject.write(",")
if i fileObject.write("\n")
fileObject.close()
def main():
maskShp = r"E:\xxx\xxx.shp"
valueField = "xxxID"
polygon = r"E:\xxx\xxx.shp"
dstFile = r"E:\xxx\xxx.shp"
tempDir = "E:\\xxx\\xxx\\temp\\"
masksDir = tempDir+"masks\\"
shutil.rmtree(tempDir)
os.mkdir(tempDir)
os.mkdir(masksDir)
table = [[valueField,"AREA"]]
arcpy.Split_analysis(maskShp,maskShp,valueField,masksDir)
for mask in os.listdir(masksDir):
if not os.path.splitext(mask)[1] == ".shp":
continue
print(mask)
value = mask[0:-4]
tempFile = value+".shp"
arcpy.Clip_analysis(polygon,masksDir+mask,tempDir+tempFile)
featureCount = arcpy.GetCount_management(tempDir+tempFile)[0]
if int(featureCount)>0:
arcpy.CalculateAreas_stats(tempDir+tempFile, tempDir+"AREA_"+tempFile)
arcpy.TableToTable_conversion(tempDir+"AREA_"+value+".dbf",tempDir,value+".csv")
areaSum = pandas.read_csv(tempDir+value+".csv")["F_AREA"].sum()
print(areaSum)
table.append([value,areaSum])
print(table)
write2DListToCsv(table,dstFile)
if __name__ == "__main__":
main()