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