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

[Javier de la Torre]

Thanks Paul,

Now, it is a pity there is no way to automate this to a decent level.

Is there nobody working on a generic shapefile importer that can accept any sort of shapefile and try to load it using the .prj file?

I will keep trying :)

Javier.

On Aug 15, 2011, at 5:55 PM, Paul Ramsey wrote:

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