[postgis-users] Projection questions

Paul Ramsey pramsey at refractions.net
Wed Jun 22 15:59:14 PDT 2005


Hi Debbie,
You would first INSERT a new row into the SPATIAL_REF_SYS table, with 
the WKT, the PROJ4 and other info, along with whatever SRID you like.
Then run shp2pgsql using the -s flag and your new made-up SRID. Having 
done that, you will find that the transform() function works as advertised.
If you are never planning on using transform(), you can just ignore your 
SRID, but that's bad data management practice (what if someone else 
comes along and wants to use your data).
Yours,
Paul

Pagurek,Debbie [NCR] wrote:

> I have just received a shapefile (polygons) that is in an Albers
> projection and from what I can see it does not have an associated EPSG
> code.
> My question is this:
> Do I have to reproject the shapefile to meet an EPSG code (e.g.42304)
> before loading the data into postgis? Or can I load it using shp2pgsql
> (what would I set the -s to?) and then somehow set up a new SRID that
> defines the projection below? How would I add a new SRID to the
> SPATIAL_REF_SYS Table?
> In the end, I would like to have this in postgis and serve it as WMS/WFS
> using Mapserver. 
> 
> Just trying to find the easiest way to deal with this non-EPSG
> projection.
> 
> Here is the projection in WKT SRS representation:
> Layer name: pda_al83_polygon
> Geometry: Polygon
> Feature Count: 5013
> Extent: (3668545.746945, 620284.375175) - (9021698.683472,
> 5271375.948747)
> Layer SRS WKT:
> PROJCS["NAD_1983_Albers",
>     GEOGCS["GCS_North_American_1983",
>         DATUM["North_American_Datum_1983",
>             SPHEROID["GRS_1980",6378137.0,298.257222101]],
>         PRIMEM["Greenwich",0.0],
>         UNIT["Degree",0.0174532925199433]],
>     PROJECTION["Albers_Conic_Equal_Area"],
>     PARAMETER["False_Easting",6200000.0],
>     PARAMETER["False_Northing",3000000.0],
>     PARAMETER["longitude_of_center",-91.86666666666666],
>     PARAMETER["Standard_Parallel_1",49.0],
>     PARAMETER["Standard_Parallel_2",77.0],
>     PARAMETER["latitude_of_center",63.390675],
>     UNIT["Meter",1.0]]
> 
> Thanks in advance for advice that anyone can offer.
> 
> D. Pagurek
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users




More information about the postgis-users mailing list