arcpy计算当前地图文档的每个数据框每个图层的面积

计算当前地图文档的每个数据框每个图层的面积，基于数据源的坐标系

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()