Я пытаюсь вычислить площадь многоугольника в моем скрипте Python. Я создаю новый многоугольник из слияния двух вместе, и я хотел бы добавить область полученного многоугольника к полю в выходном файле. Многоугольник хранится в обычном шейп-файле и проецируется. Площадь желательно в единицах карты.
Я бы подумал, что это довольно распространенная и простая задача, но, несмотря на большое количество Google, я до сих пор не смог найти работающие решения.
Я планировал использовать arcpy.updateCursor
для вставки значения после его вычисления (на этом этапе в FC есть только одна особенность), поэтому проще всего будет вернуть его в виде переменной. Любое альтернативное решение, которое выполняет ту же задачу (получение значения площади в правильном поле), также будет работать.
Я также попробовал Калькулятор поля от Python. Модифицированный со страниц справки я думал, что следующее будет работать, но пока не повезло.
arcpy.AddField_management(tempPgs, "Shape_area", 'DOUBLE')
exp = "float(!SHAPE.AREA!.split())"
arcpy.CalculateField_management(tempPgs, "Shape_area", exp)
Запуск ArcGIS Basic 10.1 SP1 с Python 2.7 в Windows 7.
Соответствующие части моего текущего кода выглядят так:
#/.../
arcpy.Copy_management(inpgs, outpgs)
arcpy.AddField_management(outpgs, 'Shape_area', 'LONG')
fields = AM.FieldLst(outpgs)
#/.../
# Identify and search for shapes smaller than minimum area
where1 = '"' + 'Shape_Area' + '" < ' + str(msz)
polyrows = arcpy.SearchCursor(inpgs, where1)
for prow in polyrows:
grd1 = prow.GridID # GridID on the current polygon
grd2 = nDD.get(grd1) # GridID on the polygon downstream
# Update features
if grd2
geometry1 = prow.Shape
geometry2 = geometryDictionary[grd2]
# Update temporary features
arcpy.Merge_management([geometry1, geometry2], tempMerged)
arcpy.Dissolve_management(tempMerged, tempPgs)
fds = AM.FieldLst(tempPgs)
for field in fields[2:]:
arcpy.AddField_management(tempPgs, field, 'DOUBLE')
for fd in fds[2:]:
arcpy.DeleteField_management(tempPgs, fd)
exp = "float(!SHAPE.AREA!.split())"
arcpy.CalculateField_management(tempPgs, "Shape_area", exp)
# Append them to output FC
try:
arcpy.Append_management(tempPgs, outpgs, "TEST")
except arcgisscripting.ExecuteError:
arcpy.Append_management(tempPgs, outpgs, "NO_TEST")
elif ...
else ...
SHAPE@AREA
как часть вашего курсора, чтобы прочитать область; но структура кода зависит от того, находится ли ваша область в тех единицах, которые вы хотите выписать.