I am looking to transform the CRS of the geometry from my initial shapefile data to standard lat-long coordinates - espg 4236.
Here's the geopandas dataframe:
I wrote some script below that transforms the geometry to the espg CRS.
geodata.crs
{'proj': 'lcc',
'lat_1': 37.06666666666667,
'lat_2': 38.43333333333333,
'lat_0': 36.5,
'lon_0': -120.5,
'x_0': 2000000,
'y_0': 500000.0000000002,
'ellps': 'GRS80',
'units': 'us-ft',
'no_defs': True}
pointsEB is point datatype:
POINT (6091769.359598021 2090180.343520682)
# Function to transform CRS projection
inProj = Proj(geodata.crs)
outProj = Proj(init='epsg:4326')
pointsEBt = []
for pt in pointsEB:
x,y = transform(inProj, outProj, pt.x, pt.y)
pointsEBt.append(Point(x,y))
Transformed Point:
pointsEBt[1].x, pointsEBt[1].y
(-68.42841348959915, 40.16978821064519) --> Located in Antartica.
My question: There seems to be something going on with the conversion to lat-long CRS. I used this as reference and "+init=epsg:4326" seems to be the format. However, when I plot a transformed point on the map the location is not what I expect it to be, far from it.