[postgis-users] Importing a shape file with NAD83 / Conus Albers

[Javier de la Torre]

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
> <jatorre at 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 at postgis.refractions.net
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> 
>> 
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users