Naar aanleiding van een vraag over de formule voor de andere kant op, constateerde ik dat er een fout in bovenstaande formule voor de y-coördinaat zit. De juiste formules beide kanten op staan hieronder.
Van ZXY van een tile naar xy-coördinaten in RD:
t = (903401.92-22598.08)*0.5^Z 'tegelgrootte in meters'
x = X*t-285401.92 'x-coördinaat in meters van de linkeronderhoek van de tile'
y = -(Y+1)*t+903401.92 'y-coördinaat in meters van de linkeronderhoek van de tile'
De andere kant op, van xy-coördinaten in RD naar XY-index van een tile bij zoomniveau Z:
t = (903401.92-22598.08)*0.5^Z 'tegelgrootte in meters'
X = floor((285401.92+x)/t) 'rij-index van de tile'
Y = floor((903401.92-y)/t) 'kolomindex van de tile'
Waar bij de getallen in bovenstaande formules uit het document over de richtlijn komen:
903401.92 'y (links)boven bij zoomniveau 0'
22598.08 'y (rechts)onder bij zoomniveau 0'
-285401.92 'x links(boven) bij zoomniveau 0'