# Projecting single point in arcpy

577
2
01-27-2022 01:18 PM
Occasional Contributor III

I have a script that takes in lot, long values (in decimal degrees) and determines whether the points are contained within a polygon in another feature class.  I assume I have to create a point object and get it into the same coordinate system as the polygon, but the projection comes back empty (no points).  I don't fully understand all the ins and outs of projections so hopefully I'm doing something obviously incorrect here.

``````x = -88.155675
y = 43.223869
f_cs = 'WGS 1984 Web Mercator (auxiliary sphere)'
t_cs = 'NAD 1983 UTM Zone 16N'
pt = arcpy.Point(x,y)
pg = arcpy.PointGeometry(pt, arcpy.SpatialReference(f_cs))
pgp = pg.projectAs(t_cs, xform)
print (pgp.firstPoint)``````

Tags (1)
1 Solution

Accepted Solutions
MVP Esteemed Contributor

you x, y aren't in web mercator,

``````x = -88.155675
y = 43.223869
sr = arcpy.SpatialReference(4326)
t_cs = 'NAD 1983 UTM Zone 16N'
pt = arcpy.Point(x,y)
pg = arcpy.PointGeometry(pt, sr)
srto = arcpy.SpatialReference('NAD 1983 UTM Zone 16N')
pgp = pg.projectAs(srto, xform)

pgp.firstPoint
<Point (406145.6122366393, 4786322.940660958, #, #)>

pt
<Point (-88.155675, 43.223869, #, #)>``````

seems to work

... sort of retired...
2 Replies
MVP Esteemed Contributor

you x, y aren't in web mercator,

``````x = -88.155675
y = 43.223869
sr = arcpy.SpatialReference(4326)
t_cs = 'NAD 1983 UTM Zone 16N'
pt = arcpy.Point(x,y)
pg = arcpy.PointGeometry(pt, sr)
srto = arcpy.SpatialReference('NAD 1983 UTM Zone 16N')
pgp = pg.projectAs(srto, xform)

pgp.firstPoint
<Point (406145.6122366393, 4786322.940660958, #, #)>

pt
<Point (-88.155675, 43.223869, #, #)>``````

seems to work

... sort of retired...
Occasional Contributor III

That's it! I must have tried a hundred combinations and at some point I was using 4326 but must have had something else off at that point.  Thanks Dan!