How can I convert a projection that's not part of spatial_ref_sys?

Posted by Summer on Stack Overflow See other posts from Stack Overflow or by Summer
Published on 2010-04-22T22:25:44Z Indexed on 2010/04/22 22:33 UTC
Read the original article Hit count: 496

Filed under:
|
|
|

Hi,

I'm importing shapefiles into a Postgres+PostGIS database.

Here's my usual procedure:
* Find an srid in the spatial_ref_sys table where srtext appears to match the shapefile's .prj file
* Upload the data into a new table using the shp2pgsql utility, specifying the srid using the -s flag
* Add the new table to my main geometry table, and on the way convert to an srid of 4269 (the Census standard projection) using ST_Transform

Unfortunately, the spatial_ref_sys table doesn't include Mississippi state's standard projection. The contents of their .prj file is as follows, where I've bolded the parts I usually try to match:

PROJCS["mstm",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",1300000.0],PARAMETER["Central_Meridian",-89.75],PARAMETER["Scale_Factor",0.9998335],PARAMETER["Latitude_Of_Origin",32.5],UNIT["Meter",1.0]]

I eventually found the ogr2ogr utility, and especially with the "peace and joy" promises, I decided to give it a try. I tried this command:

ogr2ogr -update -f "PostgreSQL" PG:"Connection details" 
    "File name.shp" -t_srs EPSG:4269 -nln Table_Name

I am now getting the error "Terminating translation prematurely after failed translation of layer" -- which seems to indicate that ogr2ogr is not going to be the savior I imagined in getting arbitrary .prj files neatly into the 4269 projection.

Any ideas about what to do?

© Stack Overflow or respective owner

Related posts about gis

Related posts about postgresql