Very odd. Well, the proj4text you want is: +proj=aea +lon_0=-96 +lat_0=23 +lat_1=29.5 +lat_2=45.5 +datum=NAD83
But, why OGR cannot produce that? Don't know. P On Mon, Aug 15, 2011 at 2:38 PM, Javier de la Torre <jato...@vizzuality.com> wrote: > 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 > _______________________________________________ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users
