<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi all,<div><br></div><div>I got a shapefile with a .prj file on it. I pasted the contents to <a href="http://prj2epsg.org/">http://prj2epsg.org</a> and got that this <a href="http://prj2epsg.org/epsg/5070">http://prj2epsg.org/epsg/5070</a></div><div><br></div><div>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 <a href="http://spatialreference.org/ref/sr-org/7178/">http://spatialreference.org/ref/sr-org/7178/</a></div><div><br></div><div>But it does not seem to include the proj4text data needed for postgis to do transformations.</div><div><br></div><div>I then tried loading it in OGR with Python doing this:</div><div><br></div><div><div>>>> srs = osr.SpatialReference()</div><div>>>> 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"]]'</div><div>>>> srs.ImportFromWkt([wkt])</div></div><div>>>>srs.ExportToProj4()</div><div><br></div><div>But I get "ERROR 6: No translation for Albers Equal Area to PROJ.4 format is known."</div><div><br></div><div>Anybody knows how can I get the PROJ4TEXT from this?</div><div><br></div><div>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.</div><div><br></div><div><br></div><div>Thanks for any advice,</div><div><br></div><div>Javier.</div><div><br></div></body></html>