[MapProxy] **SPAM** WMS EPSG:102067 reprojection problem

[Martin Kokeš]

Hello list,

I would like to setup mapproxy as temporary tile proxy for a WMS server. This WMS service (http://geoportal.cuzk.cz/WMS_ORTOFOTO_PUB/service.svc/get?service=WMS&request=getCapabilities) have a source data in EPSG:102067. Built-in WMS reprojection to EPSG:3857 / EPSG:900913 however fails with shifted coordinates by approx. 120 m NE.

So I want to get WMS tiles with EPSG:102067, reproject and cache it on my server in EPSG:900913 and then serve as tiles. However can't get correctly reprojected files.

Instead a request like this:
http://geoportal.cuzk.cz/WMS_ORTOFOTO_PUB/service.svc/get?SERVICE=WMS&REQUEST=GetMap&VERSION=1.3.0&LAYERS=GR_ORTFOTORGB&FORMAT=image%2Fjpeg&WIDTH=1218&HEIGHT=707&BBOX=-647115.8771133041,-1067104.2684050582,-643930.8657550423,-1065255.497444228&CRS=EPSG%3A102067&STYLES=

I get in source-requests.log request like this:
http://geoportal.cuzk.cz/WMS_ORTOFOTO_PUB/service.svc/get?styles=&format=image%2Fpng&height=1184
&bbox=-22072429.9911,-40399557.4804,-22058115.041,-40385241.7224&layers=GR_ORTFOTORGB&crs=EPSG%3A102067&service=WMS&request=Ge
tMap&width=1184&version=1.3.0

As you see, bbox is somehow wrong.

I have /usr/share/proj/epsg with correct EPSG:102067 definition:
# S-JTSK Krovak East North
<102067> +proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=30.28813972222222 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +pm=greenwich +units=m +no_defs <>

One reference point transformation is correct:
echo "-868208.53 -1095793.57 512.30" | cs2cs +init=esri:102067 \
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:4326

12d48'25.16"E	49d27'8.146"N 559.261

My mapproxy.yaml is here: http://pastebin.com/zBn2U3iS

Anyone can look at this and give me some advice, what I should check?

Thanks

Martin