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
