Hi all,
I've been trying to set up the GDAL WMS driver to work with an ArcGIS tile
service, as described at
https://gdal.org/drivers/raster/wms.html#arcgis-rest-api
Everything works fine for Web Mercator services, but I'm not able to get this
to work with a different projection.
I've used the following to create the XML service configuration file:
gdal_translate
https://utility.arcgis.com/usrsvcs/servers/6a6319c6b6bc4876b742eb907642da2e/rest/services/MapGeniePremiumITM/MapServer?f=json
wms.xml -of WMS
This outputs:
<GDAL_WMS>
<Service name="TMS">
<ServerUrl>https://utility.arcgis.com/usrsvcs/servers/6a6319c6b6bc4876b742eb907642da2e/rest/services/MapGeniePremiumITM/MapServer/tile/${z}/${y}/${x}</ServerUrl>
</Service>
<DataWindow>
<UpperLeftX>-5022200.00000000</UpperLeftX>
<UpperLeftY>4821100.00000000</UpperLeftY>
<LowerRightX>-4751266.12479892</LowerRightX>
<LowerRightY>4550166.12479892</LowerRightY>
<TileLevel>12</TileLevel>
<TileCountX>1</TileCountX>
<YOrigin>top</YOrigin>
</DataWindow>
<Projection>EPSG:2157</Projection>
<BlockSizeX>256</BlockSizeX>
<BlockSizeY>256</BlockSizeY>
<Cache/>
</GDAL_WMS>
The following should then produce valid output (a tiny test area), and does
when pointing at an XML config file for the Web Mercator services:
SET CPL_DEBUG=ON
gdal_translate -of PNG -outsize 512 512 -projwin 611540 640119 611640
640019 -projwin_srs "EPSG:2157" wms.xml output.png
For the above XML file the output image is empty. I've played around with all
options including setting the DataWindow to the full extent referenced in the
service JSON file with no luck:
<UpperLeftX>399999.908983152</UpperLeftX>
<LowerRightY>499999.966133265</LowerRightY>
<LowerRightX>810000.120650242</LowerRightX>
<UpperLeftY>999999.966133265</UpperLeftY>
Changing to the world bounds at least makes requests, but the X and Y values
are way off:
<UpperLeftX>-20037508.34</UpperLeftX>
<UpperLeftY>20037508.34</UpperLeftY>
<LowerRightX>20037508.34</LowerRightX>
<LowerRightY>-20037508.34</LowerRightY>
Requesting... /tile/12/2948/2113
Actual tiles are more in the range /tile/12/61727/83181
Setting TileX and TileY can increase these values, but I'm not sure how (or if)
I should be calculating them.
Some questions:
- Are custom projections and DataWindows even supported?
- Are the values output to the XML config correct, or are they simply a
template from which to start with?
- Are all the values for the XML in the service JSON file, or are there
additional calculations to make?
Any advice greatly appreciated,
Seth
--
web:https://geographika.net & https://mapserverstudio.net
twitter: @geographika
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
