Want Python code for what Calculate Geonetry tool does for xy to meters

AOrlich

Hello

I have points with XY coordinates in GCS84 and need to convert them to meters before extracting raster cell values from the points. I know I can do this with the Calculate Geometry window from the drop-down menu off the field name in the attribute table, but I'd like to add this step into a longer python code after adding the new field.

Should be easy, right? The Calculate Geometry window asks for the following input (I've entered mine):
Property: Y Coordinate of Point
Coordinate System: Use coordinate system of the data frame (PCS:Polar Stereographic)
Units: Meters

Thanks in advance for any ideas or help!
~Alice

JakeSkinner

Hi Alice,

You can do this using the UpdateCursor.  You can specify a coordinate system that's in meters.  Ex:

fc = "Points"

sr = arcpy.SpatialReference(3032)

with arcpy.da.UpdateCursor(fc, ["SHAPE@XY", "XCOORD", "YCOORD"], "", sr) as cursor:
    for row in cursor:
        row[1] = row[0][0]
        row[2] = row[0][1]
        cursor.updateRow(row)

del row, cursor

You can find the code of your coordinate system within the data frame properties > Coordinate System tab:

[ATTACH=CONFIG]30312[/ATTACH]

screen.png