[gdal-dev] [mapserver-dev] Wrong .prj file with SHAPEZIP Outputformat

Thomas Gratier osgeo.mailinglist at gmail.com
Thu May 16 04:29:10 PDT 2013


Hello,

Just a guess: the Mapserver projection file comes maybe from GEOS and GDAL
retrieve the wrong proj (deduce from the mapfile) after (although with your
analysis it doesn't seem so) >>
http://mapserver.org/faq.html#where-do-i-find-my-epsg-code

Regards

Thomas Gratier


2013/5/16 Guillaume Sueur <no-reply at neogeo-online.net>

> Hi list(s),
>
> Sorry for cross-posting gdal and mapserver dev lists, but I don't know
> exactly where my problem comes from.
>
> When using a SHAPEZIP OUTPUTFORMAT with MapServer, I retrieve a prj file
> which is not exactly the correct projection definition (for EPSG:3946) :
>
> PROJCS["Lambert_Conformal_**Conic",GEOGCS["GCS_GRS 1980(IUGG,
> 1980)",DATUM["D_unknown",**SPHEROID["GRS80",6378137,298.**
> 257222101]],PRIMEM["Greenwich"**,0],UNIT["Degree",0.**
> 017453292519943295]],**PROJECTION["Lambert_Conformal_**
> Conic"],PARAMETER["standard_**parallel_1",45.25],PARAMETER["**
> standard_parallel_2",46.75],**PARAMETER["latitude_of_origin"**
> ,46],PARAMETER["central_**meridian",3],PARAMETER["false_**
> easting",1700000],PARAMETER["**false_northing",5200000],UNIT[**"Meter",1]]
>
> I took a look at the code in mapogroutput.c and noticed it was all driven
> by OGR.
>
> So I did the same extract to a EPSG:3946 shapefile with ogr2ogr. It gives
> the correct result :
>
> PROJCS["RGF93_Lambert_93",**GEOGCS["GCS_RGF93",DATUM["D_**
> RGF_1993",SPHEROID["GRS_1980",**6378137,298.257222101]],**
> PRIMEM["Greenwich",0],UNIT["**Degree",0.017453292519943295]]**
> ,PROJECTION["Lambert_**Conformal_Conic"],PARAMETER["**
> standard_parallel_1",49],**PARAMETER["standard_parallel_**
> 2",44],PARAMETER["latitude_of_**origin",46.5],PARAMETER["**
> central_meridian",3],**PARAMETER["false_easting",**
> 700000],PARAMETER["false_**northing",6600000],UNIT["**Meter",1]]
>
> which makes the shapefile's projection understandable by QGIS.
>
> So, the question is : where does the prj definition generated by mapserver
> come from ? Why is it wrong and how to fix this ?
>
> Best regards
>
> Guillaume
> ______________________________**_________________
> mapserver-dev mailing list
> mapserver-dev at lists.osgeo.org
> http://lists.osgeo.org/**mailman/listinfo/mapserver-dev<http://lists.osgeo.org/mailman/listinfo/mapserver-dev>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20130516/55e1bad8/attachment.html>


More information about the gdal-dev mailing list