-
Notifications
You must be signed in to change notification settings - Fork 19
Description
I'm converting an Albers projection from WKT to Proj4. My WKT is:
PROJCS["NAD_1983_Albers",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.2572221010042,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",29.5],PARAMETER["standard_parallel_2",45.5],PARAMETER["latitude_of_center",23],PARAMETER["longitude_of_center",-96],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]
Then this line:
pycrs.parse.from_ogc_wkt(wkt).to_proj4()
Gives:
'+proj=aea +datum=NAD83 +ellps=GRS80 +a=6378137.0 +rf=298.2572221010042 +pm=0 +lat_1=29.5 +lat_2=45.5 +lonc=-96.0 +x_0=0.0 +y_0=0.0 +units=m +axis=enu +no_defs'
As far as I am aware, +lonc should be +lon_0