You mean: srs.MorphFromESRI() srs.ExportToProj4()
That still produce srs.ExportToProj4() ERROR 6: No translation for Albers Equal Area to PROJ.4 format is known. On Aug 15, 2011, at 5:28 PM, Paul Ramsey wrote: > Try running MorphFromESRI on the srs before generating the proj4? > P > > On Mon, Aug 15, 2011 at 2:25 PM, Javier de la Torre > <jato...@vizzuality.com> wrote: >> Hi all, >> I got a shapefile with a .prj file on it. I pasted the contents >> to http://prj2epsg.org and got that this http://prj2epsg.org/epsg/5070 >> This srs does not seem to be in PostGIS so I am trying to add it. Now, the >> first thing I tried was sending it to spatialreference and I got >> this http://spatialreference.org/ref/sr-org/7178/ >> But it does not seem to include the proj4text data needed for postgis to do >> transformations. >> I then tried loading it in OGR with Python doing this: >>>>> srs = osr.SpatialReference() >>>>> wkt = 'PROJCS["NAD83 / Conus Albers", GEOGCS["NAD83", DATUM["North >>>>> American Datum 1983", SPHEROID["GRS 1980", 6378137.0, 298.257222101, >>>>> AUTHORITY["EPSG","7019"]], TOWGS84[1.0, 1.0, -1.0, 0.0, 0.0, 0.0, >>>>> 0.0], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich", 0.0, >>>>> AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], >>>>> AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], >>>>> AUTHORITY["EPSG","4269"]], PROJECTION["Albers Equal Area", >>>>> AUTHORITY["EPSG","9822"]], PARAMETER["central_meridian", -96.0], >>>>> PARAMETER["latitude_of_origin", 23.0], PARAMETER["standard_parallel_1", >>>>> 29.5], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", >>>>> 0.0], PARAMETER["standard_parallel_2", 45.5], UNIT["m", 1.0], >>>>> AXIS["Easting", EAST], AXIS["Northing", NORTH], >>>>> AUTHORITY["EPSG","5070"]]' >>>>> srs.ImportFromWkt([wkt]) >>>>> srs.ExportToProj4() >> But I get "ERROR 6: No translation for Albers Equal Area to PROJ.4 format is >> known." >> Anybody knows how can I get the PROJ4TEXT from this? >> Ideally what I would like is to inspect automatially a shapefile coming with >> a .prj file, try to find if it is in PostGIS and if not register the proj >> and import the data using shp2psql. >> >> Thanks for any advice, >> Javier. >> >> _______________________________________________ >> postgis-users mailing list >> postgis-users@postgis.refractions.net >> http://postgis.refractions.net/mailman/listinfo/postgis-users >> >> > _______________________________________________ > postgis-users mailing list > postgis-users@postgis.refractions.net > http://postgis.refractions.net/mailman/listinfo/postgis-users _______________________________________________ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users
