计算当前地图文档的每个数据框每个图层的面积,基于数据源的坐标系
import arcpy
import os
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
#设置输出文件
outputfile="D:\\aaa\\output.txt"
#获取当前地图文档下的数据框
mxd = arcpy.mapping.MapDocument('current')
dfs = arcpy.mapping.ListDataFrames(mxd)
#打开输出文件
f = open(outputfile,'w')
#遍历数据框
for df in dfs:
#f.write("%s\n"%df.name)
#遍历当前数据框的每个图层,获取其数据源并计算面积
lyrs= arcpy.mapping.ListLayers(mxd,"",df)
for lyr in lyrs:
#获取图层的数据源,并计算每个数据的面积
fc = lyr.dataSource
fields=["SHAPE@AREA"]
#计算面积
areas=[]
with arcpy.da.SearchCursor(fc, fields) as cursor:
for row in cursor:
areas.append(row[0])
#面积写入到文件,每行的格式为,制表符分割:
#数据框名称-图层名称-面积
f.write("%s\t%s\t%s\n"%(df.name,lyr.name,sum(areas)))
#关闭结果文件
f.close()