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

[Frank Warmerdam]

Guillaume,

Yes, I believe this is exactly the issue.  Good luck!

Best regards,
Frank

On Thu, May 16, 2013 at 5:38 AM, Guillaume Sueur <no-reply at neogeo-online.net
> wrote:

> 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.
>
> Gonna try to hack it a bit...
>
>
> Le 16/05/13 13:29, Thomas Gratier a écrit :
>
>> 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<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
>> <mailto:no-reply at neogeo-**online.net <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 <mailto:mapserver-dev at lists.**osgeo.org<mapserver-dev at lists.osgeo.org>
>> >
>>     http://lists.osgeo.org/__**mailman/listinfo/mapserver-dev<http://lists.osgeo.org/__mailman/listinfo/mapserver-dev>
>>     <http://lists.osgeo.org/**mailman/listinfo/mapserver-dev<http://lists.osgeo.org/mailman/listinfo/mapserver-dev>
>> **>
>>
>>
>>  ______________________________**_________________
> 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>
>



-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam,
warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Software Developer
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-dev/attachments/20130516/a67c4466/attachment.html>