|
- AreaSum0=sum([row[0] for row in arcpy.da.UpdateCursor('DLTB0',u'椭球面积')]) #计算老DLTB总面椭球面积
- AreaSum1=sum([row[0] for row in arcpy.da.UpdateCursor('DLTB1',u'椭球面积')]) #计算分割后新DLTB总面椭球面积
- AreaAD=AreaSum1-AreaSum0 #计算要平差的面积,可正可负
- arcpy.AddField_management("DLTB1",'Scale',"DOUBLE") #添加比例scale字段
- arcpy.AddField_management("DLTB1",u'椭球面积_平差',"DOUBLE") #添加用于存放平差后面积字段
- arcpy.CalculateField_management("DLTB1","Scale","!Shape_Area!/AreaSum1","PYTHON_9.3") #计算scale字段,字段值是dltb1每个要素与总面积的比例
- arcpy.CalculateField_management("DLTB1",u"椭球面积_平差","!Scale!*AreaAD+!椭球面积!","PYTHON_9.3") #将AreaAD按照比例分别平差到各个图斑
复制代码
你好,大神。我操作了一下,怎么不行啊,说我SQL语句错误是啥意思,这段代码是VB脚本 还是 pyhone 请教了。这段代码如何使用啊,请教请教。
|
|