Re: Mapserver bumped from 7.0.7 to 7.4.1 leads to wrong projection of features

Jaap Donker <[email protected]> Tue, 19 Nov 2019 14:02:45 +0100
Newsgroups gmane.comp.gis.mapserver.user
Message-ID <CAC=rA=OjMFrDEXuyaswYR=mhTV6X29Y8MiWEsnZff_OgNWWA8w@mail.gmail.com>
--===============0554351809228778901==
Content-Type: multipart/alternative; boundary="000000000000cb14ad0597b2ae04"

--000000000000cb14ad0597b2ae04
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

Thanks for your response.

In our frontend we use the following projection string for EPSG:28992:
+proj=3Dsterea +lat_0=3D52.15616055555555 +lon_0=3D5.38763888888889
+k=3D0.9999079 +x_0=3D155000 +y_0=3D463000 +ellps=3Dbessel
+towgs84=3D565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725
+units=3Dm +no_defs

I hope these parts of the map file are informative enough:

MAP
  PROJECTION
    'init=3Depsg:28992'
  END

  OUTPUTFORMAT
    NAME "GTiff"
    DRIVER GDAL/GTiff
    MIMETYPE "image/tiff"
    IMAGEMODE INT16
    EXTENSION "tif"
  END

  OUTPUTFORMAT
    NAME "GTiff-FLOAT"
    DRIVER GDAL/GTiff
    MIMETYPE "image/tiff"
    IMAGEMODE FLOAT32
    EXTENSION "tif"
  END

  OUTPUTFORMAT
    NAME "geojson"
    DRIVER "OGR/GEOJSON"
    MIMETYPE "application/json; subtype=3Dgeojson; charset=3Dutf-8"
    FORMATOPTION "STORAGE=3Dmemory"
    FORMATOPTION "FORM=3DSIMPLE"
  END
...
  LAYER
    METADATA
      "wfs_srs"            "EPSG:4326"
      "wfs_enable_request" "*"
      "wfs_getfeature_formatlist" "geojson"
      "gml_include_items"  "all"
      "wfs_maxfeatures"    "100"  # Max features to load
    END
    TYPE POLYGON
    STATUS ON
    PROJECTION
      "init=3Depsg:4326"
    END
  END
END

Kind regards,

Jaap


On Thu, 14 Nov 2019 at 18:53, Rahkonen Jukka (MML) <
[email protected]> wrote:

> Hi,
>
>
>
> Please add relevant parts of your mapfile. Obviously all that relates to
> projections and an example of the client side requests as well if you are
> changing the coordinate system on the fly.
>
>
>
> -Jukka Rahkonen-
>
>
>
> *L=C3=A4hett=C3=A4j=C3=A4:* mapserver-users <mapserver-users-bounces@list=
s.osgeo.org> *Puolesta
> *Jaap Donker
> *L=C3=A4hetetty:* torstai 14. marraskuuta 2019 17.51
> *Vastaanottaja:* [email protected]
> *Aihe:* [mapserver-users] Mapserver bumped from 7.0.7 to 7.4.1 leads to
> wrong projection of features
>
>
>
> Hi,
>
>
>
> Upon updating mapserver, projections for the features we use in our webap=
p
> were off.
>
> We updated from version 7.0.7 to version 7.4.1. I'd like to stress that
> nothing else changed in setup during testing.
>
> In the images below you can see the situation before (left) and after
> (right). The blue lines are polygons representing seamlines for the
> orthomosaic of aerials we are using.
>
> This change in projection is troubling.
>
>
>
> Original situation:
>
> https://pasteboard.co/IGFtgdY.png
>
>
>
> New situation:
>
> https://pasteboard.co/IGFtxfs.png
>
>
>
> Going through the backlog I see that there's some changes in PROJ4, but
> I'm a bit at a loss as to how best to investigate this issue.
>
>
>
> Any help or guidance would be very much appreciated.
>
>
>
> Regards,
>
>
>
> Jaap
>

--000000000000cb14ad0597b2ae04
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr"><div dir=3D"ltr"><pre>Thanks for your response.=20

In our frontend we use the following projection string for EPSG:28992:
+proj=3Dsterea +lat_0=3D52.15616055555555 +lon_0=3D5.38763888888889 +k=3D0.=
9999079 +x_0=3D155000 +y_0=3D463000 +ellps=3Dbessel +towgs84=3D565.417,50.3=
319,465.552,-0.398957,0.343988,-1.8774,4.0725 +units=3Dm +no_defs

I hope these parts of the map file are informative enough:

MAP
  PROJECTION
    &#39;init=3Depsg:28992&#39;
  END

  OUTPUTFORMAT
    NAME &quot;GTiff&quot;
    DRIVER GDAL/GTiff
    MIMETYPE &quot;image/tiff&quot;
    IMAGEMODE INT16
    EXTENSION &quot;tif&quot;
  END

  OUTPUTFORMAT
    NAME &quot;GTiff-FLOAT&quot;
    DRIVER GDAL/GTiff
    MIMETYPE &quot;image/tiff&quot;
    IMAGEMODE FLOAT32
    EXTENSION &quot;tif&quot;
  END

  OUTPUTFORMAT
    NAME &quot;geojson&quot;
    DRIVER &quot;OGR/GEOJSON&quot;
    MIMETYPE &quot;application/json; subtype=3Dgeojson; charset=3Dutf-8&quo=
t;
    FORMATOPTION &quot;STORAGE=3Dmemory&quot;
    FORMATOPTION &quot;FORM=3DSIMPLE&quot;
  END
...
  LAYER
    METADATA
      &quot;wfs_srs&quot;            &quot;EPSG:4326&quot;
      &quot;wfs_enable_request&quot; &quot;*&quot;
      &quot;wfs_getfeature_formatlist&quot; &quot;geojson&quot;
      &quot;gml_include_items&quot;  &quot;all&quot;
      &quot;wfs_maxfeatures&quot;    &quot;100&quot;  # Max features to loa=
d
    END
    TYPE POLYGON
    STATUS ON
    PROJECTION
      &quot;init=3Depsg:4326&quot;
    END
  END
END

Kind regards,

Jaap</pre></div><br><div class=3D"gmail_quote"><div dir=3D"ltr" class=3D"gm=
ail_attr">On Thu, 14 Nov 2019 at 18:53, Rahkonen Jukka (MML) &lt;<a href=3D=
"mailto:[email protected]" target=3D"_blank">jukka.rahkon=
[email protected]</a>&gt; wrote:<br></div><blockquote class=3D"gmail_=
quote" style=3D"margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,=
204);padding-left:1ex">





<div lang=3D"FI">
<div>
<p class=3D"MsoNormal"><span>Hi,<u></u><u></u></span></p>
<p class=3D"MsoNormal"><span><u></u>=C2=A0<u></u></span></p>
<p class=3D"MsoNormal"><span lang=3D"EN-US">Please add relevant parts of yo=
ur mapfile. Obviously all that relates to projections and an example of the=
 client side requests as well if you are changing the coordinate system on =
the
 fly.<u></u><u></u></span></p>
<p class=3D"MsoNormal"><span lang=3D"EN-US"><u></u>=C2=A0<u></u></span></p>
<p class=3D"MsoNormal"><span lang=3D"EN-US">-Jukka Rahkonen-<u></u><u></u><=
/span></p>
<p class=3D"MsoNormal"><span lang=3D"EN-US"><u></u>=C2=A0<u></u></span></p>
<p class=3D"MsoNormal"><b>L=C3=A4hett=C3=A4j=C3=A4:</b> mapserver-users &lt=
;<a href=3D"mailto:[email protected]" target=3D"_blan=
k">[email protected]</a>&gt;
<b>Puolesta </b>Jaap Donker<br>
<b>L=C3=A4hetetty:</b> torstai 14. marraskuuta 2019 17.51<br>
<b>Vastaanottaja:</b> <a href=3D"mailto:[email protected]" ta=
rget=3D"_blank">[email protected]</a><br>
<b>Aihe:</b> [mapserver-users] Mapserver bumped from 7.0.7 to 7.4.1 leads t=
o wrong projection of features<u></u><u></u></p>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
<div>
<div>
<p class=3D"MsoNormal">Hi,<u></u><u></u></p>
</div>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">Upon updating mapserver, projections for the feature=
s we use in our webapp were off.<u></u><u></u></p>
</div>
<p class=3D"MsoNormal">We updated from version 7.0.7 to version 7.4.1. I&#3=
9;d like to stress that nothing else changed in setup during testing.<br>
<br>
In the images below you can see the situation before (left) and after (righ=
t). The blue lines are polygons representing seamlines for the orthomosaic =
of aerials we are using.<u></u><u></u></p>
<div>
<p class=3D"MsoNormal">This change in projection is troubling.<u></u><u></u=
></p>
</div>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">Original situation:<u></u><u></u></p>
</div>
<p class=3D"MsoNormal"><a href=3D"https://pasteboard.co/IGFtgdY.png" target=
=3D"_blank">https://pasteboard.co/IGFtgdY.png</a>=C2=A0=C2=A0
<u></u><u></u></p>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">New situation:<u></u><u></u></p>
</div>
<div>
<p class=3D"MsoNormal"><a href=3D"https://pasteboard.co/IGFtxfs.png" target=
=3D"_blank">https://pasteboard.co/IGFtxfs.png</a><u></u><u></u></p>
</div>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">Going through the backlog I see that there&#39;s som=
e changes in PROJ4, but I&#39;m a bit at a loss as to how best to investiga=
te this issue.<u></u><u></u></p>
</div>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">Any help or guidance would be very much appreciated.=
<u></u><u></u></p>
</div>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">Regards,<u></u><u></u></p>
</div>
<div>
<p class=3D"MsoNormal"><u></u>=C2=A0<u></u></p>
</div>
<div>
<p class=3D"MsoNormal">Jaap <u></u><u></u></p>
</div>
</div>
</div>
</div>

</blockquote></div></div>

--000000000000cb14ad0597b2ae04--

--===============0554351809228778901==
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline

X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX18KbWFwc2VydmVy
LXVzZXJzIG1haWxpbmcgbGlzdAptYXBzZXJ2ZXItdXNlcnNAbGlzdHMub3NnZW8ub3JnCmh0dHBz
Oi8vbGlzdHMub3NnZW8ub3JnL21haWxtYW4vbGlzdGluZm8vbWFwc2VydmVyLXVzZXJz

--===============0554351809228778901==--