Re: [postgis-users] Importing a shape file with NAD83 / Conus Albers
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 > 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 >>> 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
Re: [postgis-users] Importing a shape file with NAD83 / Conus Albers
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 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 >> 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
Re: [postgis-users] Importing a shape file with NAD83 / Conus Albers
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 > 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
Re: [postgis-users] Importing a shape file with NAD83 / Conus Albers
Try running MorphFromESRI on the srs before generating the proj4? P On Mon, Aug 15, 2011 at 2:25 PM, Javier de la Torre 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] Importing a shape file with NAD83 / Conus Albers
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
