[MapProxy] Placer dans la boîte de réception Plus 1 sur 552 Tout imprimer Dans une nouvelle fenêtre Reprojecting to EPSG:4326 cuts data at prime meridian

[Marc Monnerat]

Hello,

We are using MapProxy 1.8.0, the problem is seen both running under Apache
2.2 and `mapproxy-util serve-develop`.

The source data are WMTS tiles in EPSG:21781 projection. When reprojecting
to WGS1984, either as WMS or WMTS, some data are missing. This pattern is
not shown with other projection, like EPSG:3857 (webmercator).

The "minimal" config reproducing the error is in attachement (test.yaml).  On
this screeshot, we see that the reprojected data (in green) are not
displayed west of the prime origin projection prime meridian (600'000m, at
Bern).
The light grey data are show for orientation only (source is another WMS
service).



[image: Images intégrées 1]

I did the following requests against mapproxy-util serve-develop --debug -b
127.0.0.1:9001 mapproxy/test.yaml:


*WMS*

*Both the EPSG:4326 and EPSG:4258 have truncated images west of Bern*


http://localhost:9001/service?LAYERS=ch.swisstopo.pixelkarte-farbe_20151231&FORMAT=image%2Fjpeg&SRS=EPSG%3A4258&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&BBOX=3.2209693748826,44.022010038859,12.747055982558,48.753510009559&WIDTH=1208&HEIGHT=600


http://localhost:9001/service?LAYERS=ch.swisstopo.pixelkarte-farbe_20151231&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&BBOX=3.796635204651,44.219155870071,13.322721812327,48.950655840771&WIDTH=1208&HEIGHT=600



*In the opposite, the image in EPSG:3857 is fine:*


http://localhost:9001/service?LAYERS=ch.swisstopo.pixelkarte-farbe_20151231&FORMAT=image%2Fjpeg&SRS=EPSG%3A3857&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&BBOX=368677.88243586,5636562.1333018,1499245.0927522,6198102.1384258&WIDTH=1208&HEIGHT=600




*WMTS*


*As observed with WMS, data west of Bern are missing, as in these two tiles
(EPSG:4326 and EPSG:4258)*


http://localhost:9001/wmts/1.0.0/ch.swisstopo.pixelkarte-farbe/default/20151231/epsg_4326/8/266/61.jpeg


http://localhost:9001/wmts/1.0.0/ch.swisstopo.pixelkarte-farbe/default/20151231/epsg_4258/8/266/61.jpeg

*In the opposite, the image in EPSG:3857 is fine:*


http://localhost:9001/wmts/1.0.0/ch.swisstopo.pixelkarte-farbe/default/20151231/epsg_3857/8/133/90.jpeg


So if I shift the meridian to the west, form the official definition:

<21781> +proj=somerc +lat_0=46.95240555555556 *+lon_0=7.439583333333333*
 +k_0=1 *+x_0=600000* +y_0=200000 +ellps=bessel
+towgs84=674.4,15.1,405.3,0,0,0,0 +units=m +no_defs  <>

into (something very unprecise)

<21781> +proj=somerc +lat_0=46.95240555555556 *+lon_0=6.0* +k_0=1
*+x_0=450000 *+y_0=200000 +ellps=bessel +towgs84=674.4,15.1,405.3,0,0,0,0
+units=m +no_defs  <>

everything is fine, but that not the goal.

I hope, it is a bit clearer


Thank you for any hints

Marc
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapproxy/attachments/20150904/979fd12c/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: missing_west.jpg
Type: image/jpeg
Size: 45499 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/mapproxy/attachments/20150904/979fd12c/attachment-0001.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.yaml
Type: application/octet-stream
Size: 4502 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/mapproxy/attachments/20150904/979fd12c/attachment-0001.obj>