Tile nummers WMTS rest service Arcgis

Datasets
1

Ik probeer via de volgende dienst van ArcGis afbeeldingen op te vragen:
https://tiles.arcgis.com/tiles/tz8zJgIoWvr9xGSh/arcgis/rest/services/03_16k_atl20_2/MapServer?f=html&cacheKey=95768509ab583bf2
Het lukt me echter niet de juiste tile nummers te genereren. Onderstaande link volgens mij Deventer moeten opleveren:
https://tiles.arcgis.com/tiles/tz8zJgIoWvr9xGSh/arcgis/rest/services/03_16k_atl20_2/MapServer/tile/4/8/11
Bij een zoom level van 4 zou nederland in 16x16 vakken zijn opgedeeld. De “extent” die wordt aangegeven is:

>   "xmin": -5249.903000000864,
>   "ymin": 306846.1979999989,
>   "xmax": 278026.08999999985,
>   "ymax": 640598.2080000006,
>   "spatialReference": 
>     "wkid": 28992,
>     "latestWkid": 28992

Als ik dit vlak opdeel in 16 bij 16 vlakken dan zou Deventer in vak Y=8 en X=11 moeten liggen. Zoek ik het handmatig uit dan blijkt het echter dat Deventer in vak Y=8 en X=9 ligt.

Ik heb alles gecontroleerd en veel verklaringen proberen te bedenken, maar helaas tot nu toe zonder resultaat.

2

Ik heb niet echt de tijd om in je vraag te duiken, maar misschien ligt de oplossing in de gratis Qgis-plugin tiles_xyz die op basis van een WMS een vectorlaag aanmaakt met de tegels:

https://plugins.qgis.org/plugins/tiles_xyz/

1 like
4

Uiteindelijk de oplossingen gevonden:
Origin:
“x”: -285401.92
“y”: 903401.92
For this service (origin at top-left, tile size 512 px):

tileX = floor( (X - origin_x) / (tile_size * resolution[z]) )
tileY = floor( (origin_y - Y) / (tile_size * resolution[z]) )