If one uses arcpy to calculate areas for a field in metric units on a polygon shapefile in a geographic coordinate system it will work, and give you plausible results. According to the documentation it uses 'a geodesic algorithm' but I don't understand what this means. How is it producing these results in square kilometers in the following code?
result_file = settings.get_boundary_file_results_location(boundary_short_name) area_field_name = settings.NEW_FIELD_NAMES['Area square kilometers'] arcpy.AddField_management(result_file, area_field_name, "FLOAT") expression = "{0}".format("!SHAPE.area@SQUAREKILOMETERS!") arcpy.CalculateField_management(result_file, area_field_name, expression, "PYTHON_9.3") And the properties of the shapefile:
أكثر...
result_file = settings.get_boundary_file_results_location(boundary_short_name) area_field_name = settings.NEW_FIELD_NAMES['Area square kilometers'] arcpy.AddField_management(result_file, area_field_name, "FLOAT") expression = "{0}".format("!SHAPE.area@SQUAREKILOMETERS!") arcpy.CalculateField_management(result_file, area_field_name, expression, "PYTHON_9.3") And the properties of the shapefile:
أكثر...