I am trying to convert some lat and lon to x and y with projection Lambert conformal conic. I am trying the Python package
pyproj for that but getting the wrong values. The issue is that the coordinate reference system is not available in EPSG codes, I built a custom CRS based on the details here. I think that's where I am doing something wrong.
Below is the code I tried-
import pyproj proj_daymet = "+proj=lcc +lat_0=42.5 +lon_0=-100 +lat_1=25 +lat_2=60 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" #my custom CRS geog_crs = pyproj.crs.GeographicCRS() lat = 31.177685 lon = -87.103554 # the transform command below gives incorrect values pyproj.transform(geog_crs,proj_daymet , lat, lon) # it outputs: 87465941.89322367, 6696423.32671964, the correct answer is 1192729, 1114943