Transform XY coordinate to WGS84 lat, long in a spatially enabled dataframe

Question

Hi,

As part of my python script, I have a spatially enabled data frame with ~1000 points in OSGB (EPSG:27700) X,Y which i would like to convert to WGS84 lat, long.  Is there a way, using the ArcGIS API for Python, to transform the OSGB XY coordinates to Lat/Long?  Or to transform the SHAPE field of the SDF from OSGB to WGS84?

Thanks in advance!

markjones6 · Answer

Hi, sorry for taking a while to reply, but thank you for your advice, it led me to a solution that's worked for my purposes.

What I did was to create a dictionary from my x,y columns and use the following to transform from OSGB to WGS84:

gis = GIS()
wgs84_coords_dict = project(geometries = osgb_coords_dict, in_sr=27700, out_sr=4326, transformation=1314)

This creates a new coordinate dictionary - wgs84_coords_dict, which can then by merged with the remainder of the dataframe.

MehdiPira1 · Answer

@markjones6 ,

Then the code snippet below is the way to follow:

from arcgis import geometry

# access the SHAPE field in your SDF and create a list
shape_list=[row.rings[0] for row in rows.SHAPE]
# do the projection on the list
for shape in shape_list:
    result = geometry.project(geometries=shape, in_sr=3857, out_sr=4326)
    print(result)

Bear in mind that this is based on the format below:

if you couldn't get it to work, I need to see a sample of your SDF to check the format.

markjones6 · Answer

This is really interesting @MehdiPira1 , thank you!

Is it possible to do something similar where the SDF doesn't exist as a feature layer? Basically I would like to work with my data only in my notebook as i have a number of data analysis to do before I publish out a "final" dataset to a feature layer.

I was building on the work in this post - ArcGIS API for Python - geometry.project as below, but it would be great rather than transform the coordinates in this way, we could transform the SHAPE to lat/long without the need for an "existing" dataset or feature layer.

input_geoms = [{"x": -17568824.55, "y": 2428377.35}, {"x": -17568456.88, "y": 2428431.352}]
result = project(geometries=input_geoms, in_sr=3857, out_sr=4326)
result