<div dir="ltr"><br><div class="gmail_extra">Guillaume, </div><div class="gmail_extra"><br></div><div class="gmail_extra">Yes, I believe this is exactly the issue. Good luck!</div><div class="gmail_extra"><br></div><div class="gmail_extra">
Best regards,</div><div class="gmail_extra">Frank<br><br><div class="gmail_quote">On Thu, May 16, 2013 at 5:38 AM, Guillaume Sueur <span dir="ltr"><<a href="mailto:no-reply@neogeo-online.net" target="_blank">no-reply@neogeo-online.net</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I think I undestood. mapgdal.c retrieves the mapserver projection definition and translates it to proj4 string, and then ask GDAL for a WKS. I think it should be using ESPG code when available to allow GDAL to retrieve more information from its CRS database than the proj4 string allows it to.<br>
<br>
Gonna try to hack it a bit...<br>
<br>
<br>
Le 16/05/13 13:29, Thomas Gratier a écrit :<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Hello,<br>
<br>
Just a guess: the Mapserver projection file comes maybe from GEOS and<br>
GDAL retrieve the wrong proj (deduce from the mapfile) after (although<br>
with your analysis it doesn't seem so) >><br>
<a href="http://mapserver.org/faq.html#where-do-i-find-my-epsg-code" target="_blank">http://mapserver.org/faq.html#<u></u>where-do-i-find-my-epsg-code</a><br>
<br>
Regards<br>
<br>
Thomas Gratier<br>
<br>
<br>
2013/5/16 Guillaume Sueur <<a href="mailto:no-reply@neogeo-online.net" target="_blank">no-reply@neogeo-online.net</a><br></div>
<mailto:<a href="mailto:no-reply@neogeo-online.net" target="_blank">no-reply@neogeo-<u></u>online.net</a>>><div class="im"><br>
<br>
Hi list(s),<br>
<br>
Sorry for cross-posting gdal and mapserver dev lists, but I don't<br>
know exactly where my problem comes from.<br>
<br>
When using a SHAPEZIP OUTPUTFORMAT with MapServer, I retrieve a prj<br>
file which is not exactly the correct projection definition (for<br>
EPSG:3946) :<br>
<br></div>
PROJCS["Lambert_Conformal___<u></u>Conic",GEOGCS["GCS_GRS 1980(IUGG,<br>
1980)",DATUM["D_unknown",__<u></u>SPHEROID["GRS80",6378137,298._<u></u>_257222101]],PRIMEM["<u></u>Greenwich"__,0],UNIT["Degree",<u></u>0.__017453292519943295]],__<u></u>PROJECTION["Lambert_Conformal_<u></u>__Conic"],PARAMETER["standard_<u></u>__parallel_1",45.25],<u></u>PARAMETER["__standard_<u></u>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],<u></u>UNIT[__"Meter",1]]<div class="im">
<br>
<br>
I took a look at the code in mapogroutput.c and noticed it was all<br>
driven by OGR.<br>
<br>
So I did the same extract to a EPSG:3946 shapefile with ogr2ogr. It<br>
gives the correct result :<br>
<br></div>
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_<u></u>of___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]]<div class="im">
<br>
<br>
which makes the shapefile's projection understandable by QGIS.<br>
<br>
So, the question is : where does the prj definition generated by<br>
mapserver come from ? Why is it wrong and how to fix this ?<br>
<br>
Best regards<br>
<br>
Guillaume<br></div>
______________________________<u></u>___________________<br>
mapserver-dev mailing list<br>
<a href="mailto:mapserver-dev@lists.osgeo.org" target="_blank">mapserver-dev@lists.osgeo.org</a> <mailto:<a href="mailto:mapserver-dev@lists.osgeo.org" target="_blank">mapserver-dev@lists.<u></u>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>
<<a href="http://lists.osgeo.org/mailman/listinfo/mapserver-dev" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/mapserver-dev</a><u></u>><br>
<br>
<br>
</blockquote><div class="HOEnZb"><div class="h5">
______________________________<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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>---------------------------------------+--------------------------------------<br>I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>and watch the world go round - Rush | Geospatial Software Developer<br>
</div></div>