<div dir="ltr"><div><div><div>Hello,<br><br></div>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) >> <a href="http://mapserver.org/faq.html#where-do-i-find-my-epsg-code">http://mapserver.org/faq.html#where-do-i-find-my-epsg-code</a><br>

<br></div>Regards<br><br></div>Thomas Gratier<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/5/16 Guillaume Sueur <span dir="ltr"><<a href="mailto:no-reply@neogeo-online.net" target="_blank">no-reply@neogeo-online.net</a>></span><br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi list(s),<br>
<br>
Sorry for cross-posting gdal and mapserver dev lists, but I don't know exactly where my problem comes from.<br>
<br>
When using a SHAPEZIP OUTPUTFORMAT with MapServer, I retrieve a prj file which is not exactly the correct projection definition (for EPSG:3946) :<br>
<br>
PROJCS["Lambert_Conformal_<u></u>Conic",GEOGCS["GCS_GRS 1980(IUGG, 1980)",DATUM["D_unknown",<u></u>SPHEROID["GRS80",6378137,298.<u></u>257222101]],PRIMEM["Greenwich"<u></u>,0],UNIT["Degree",0.<u></u>017453292519943295]],<u></u>PROJECTION["Lambert_Conformal_<u></u>Conic"],PARAMETER["standard_<u></u>parallel_1",45.25],PARAMETER["<u></u>standard_parallel_2",46.75],<u></u>PARAMETER["latitude_of_origin"<u></u>,46],PARAMETER["central_<u></u>meridian",3],PARAMETER["false_<u></u>easting",1700000],PARAMETER["<u></u>false_northing",5200000],UNIT[<u></u>"Meter",1]]<br>


<br>
I took a look at the code in mapogroutput.c and noticed it was all driven by OGR.<br>
<br>
So I did the same extract to a EPSG:3946 shapefile with ogr2ogr. It gives the correct result :<br>
<br>
PROJCS["RGF93_Lambert_93",<u></u>GEOGCS["GCS_RGF93",DATUM["D_<u></u>RGF_1993",SPHEROID["GRS_1980",<u></u>6378137,298.257222101]],<u></u>PRIMEM["Greenwich",0],UNIT["<u></u>Degree",0.017453292519943295]]<u></u>,PROJECTION["Lambert_<u></u>Conformal_Conic"],PARAMETER["<u></u>standard_parallel_1",49],<u></u>PARAMETER["standard_parallel_<u></u>2",44],PARAMETER["latitude_of_<u></u>origin",46.5],PARAMETER["<u></u>central_meridian",3],<u></u>PARAMETER["false_easting",<u></u>700000],PARAMETER["false_<u></u>northing",6600000],UNIT["<u></u>Meter",1]]<br>


<br>
which makes the shapefile's projection understandable by QGIS.<br>
<br>
So, the question is : where does the prj definition generated by mapserver come from ? Why is it wrong and how to fix this ?<br>
<br>
Best regards<br>
<br>
Guillaume<br>
______________________________<u></u>_________________<br>
mapserver-dev mailing list<br>
<a href="mailto:mapserver-dev@lists.osgeo.org" target="_blank">mapserver-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/mapserver-dev" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/mapserver-dev</a><br>
</blockquote></div><br></div>