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

[Paul Ramsey]

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