Die Mathematik von Kartenkacheln interaktiv erklärt

Tele ← 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 → Makro
13 / **4252** / **2736**
↖	↑			↗
←	`N 51.179343° E 6.855469°`	← 3 063.27 m →	`N 51.179343° E 6.899414°`	→
	↑ 3 064.13 m ↓		↑ 3 064.13 m ↓
	`N 51.151786° E 6.855469°`	← 3 065.10 m → 9 389 064 m²	`N 51.151786° E 6.899414°`
↙	↓			↘

Die Berechnung

Aus der URL werden die geographischen Koordinaten übernommen:

Name	Parameter	Wert	Wertebereich	Erklärung zum Wertebereich
Länge in Grad	`lon`	6.8655°	-180°…+180°	Die Merkatortransformation läßt die Länge unverändert, vergrößert aber die Breite. Der Wertebereich der Breite ist so gewählt, dass die Transformation den erlaubten Bereich von Länge und Breite auf eine quadratische Fläche abbildet. Wenn jemand Sie fragen sollte, wie das Limit für die Breite berechnet wird: der ist, Taschenrechner flugs auf Degrees eingestellt, E hoch PI, davon der Arcus Tangens mal zwei, das Ganze minus 90°. Mit PHP: php -r 'echo 2atan(exp(pi())) 180/pi() - 90;', mit JavaScript 2Math.atan(Math.exp(Math.PI)) 180/Math.PI - 90.
Breite in Grad	`lat`	51.1578°	-85°…+85°
Vergrößerungsstufe	`zoom`	13	0…	Keine obere Grenze. Es sollte aber Kacheln für die gewählte Stufe geben.

Länge und Breite werden in das Bogenmaß umgerechnet (Multiplikation mit π/180°):

Name	Formel	Berechnung	Ergebnis
Länge (λ)	Länge in Grad × π ÷ 180°	6.8655 × 3.1415926535 ÷ 180	λ = 0.119825579792246
Breite (φ)	Breite in Grad × π ÷ 180°	51.1578 × 3.1415926535 ÷ 180	φ = 0.892872048051235

Länge und Breite werden durch die Mercator-Projektion transformiert:

Name	Formel	Berechnung	Ergebnis
Länge (Λ)	λ (unverändert)	0.119825579792246}	Λ = 0.119825579792246}
Breite (Φ)	log (tan φ + sec φ) = log (tan φ + 1/cos φ)	log (tan 0.892872048051235 + 1/cos 0.892872048051235) = log 1.24187500346969 + 1/0.627177646342629) = log 1.24187500346969 + 1.59444458174088) = log 2.83631958521056	Φ = 1.04250729104121

Länge und Breite werden in das Karten-Koordinatensystem transformiert.
Dieses beginnt links oben mit [0,0] und endet rechts unten mit [1,1]:

Name	Formel	Berechnung	Ergebnis
X-Position (x)	(1 + Λ/π) ÷ 2	(1 + 0.119825579792246/3.1415926535) ÷ 2	x = 0.519070833333333
Y-Position (y)	(1 - Φ/π) ÷ 2	(1 - 1.04250729104121/3.1415926535) ÷ 2	y = 0.334079811416708

Aus der Kartenposition x=0.519070833333333 y=0.334079811416708 und der Vergrößerungsstufe zoom=13 berechnen wir die Kachelnummer:

Name	Formel	Berechnung	Ergebnis
Kachel-X (tx)	floor (x × 2^zoom)	floor (0.519070833333333 × 2¹³) floor (0.519070833333333 × 8192) floor (4252.22826666667)	tx = 4252
Kachel-Y (ty)	floor (y × 2^zoom)	floor (0.334079811416708 × 2¹³) floor (0.334079811416708 × 8192) floor (2736.78181512567)	ty = 2736
Kachel-Pfad (ti)	"zoom/tx/tz"	13 / 4252 / 2736	ti = "13/4252/2736"

Anzeige der Kachel https://a.tile.openstreetmap.org/13/4252/2736.png und Kaffeepause.

Aus den Nachkommastellen 0.228266666666968 und 0.7818151256688 sowie der Kachelgröße in Pixeln tilesize=256 berechnen wir die Position des Markers:

Name	Formel	Berechnung	Ergebnis
Marker-X (px)	floor (fract (x × 2^zoom) * tilesize)	floor (fract (0.519070833333333 × 2¹³) × 256) floor (fract (0.519070833333333 × 8192) × 256) floor (fract (4252.22826666667) × 256) floor (0.228266666666968 × 256) floor (58.4362666667439)	px = 58
Marker-Y (py)	floor (fract (y × 2^zoom) * tilesize)	floor (fract (0.334079811416708 × 2¹³) × 256) floor (fract (0.334079811416708 × 8192) × 256) floor (fract (2736.78181512567) × 256) floor (0.7818151256688 × 256) floor (200.144672171213)	py = 200

Wir berechnen aus der Kachelnummer die Koordinaten der linken oberen Kachelecke.

Name	Formel	Berechnung	Ergebnis
X-Position (x)	tx ÷ 2^tz	4252 ÷ 2¹³ 4252 ÷ 8192	x = 0.51904296875
Y-Position (y)	ty ÷ 2^tz	2736 ÷ 2¹³ 2736 ÷ 8192	y = 0.333984375
Länge (Λ) (Merkator)	+(x × 2 - 1) × π	+(0.51904296875 × 2 - 1) × π 0.0380859375 × 3.1415926535	Λ = 0.11965050
Breite (Φ) (Merkator)	-(y × 2 - 1) × π	-(0.333984375 × 2 - 1) × π 0.33203125 × 3.1415926535	Φ = 1.04310693573242
Länge (λ)	Λ (unverändert)	0.11965050}	λ = 0.11965050}
Breite (φ)	2×atan(exp(Φ))-π/2	2×atan(exp(1.04310693573242))-π/2 2×atan(2.8380208792274)-π/2 2×1.23202218538612-π/2 2.46404437077225-1.57079632675	φ = 0.89324804
Länge in Grad	λ ÷ π × 180°	0.11965050} ÷ 3.1415926535 × 180°	lon = 6.855469°
Breite in Grad	φ ÷ π × 180°	0.89324804 ÷ 3.1415926535 × 180°	lat = 51.179343°

Die Koordinaten der anderen Ecken werden aus den Kachelnummern der Nachbarkacheln berechnet:

Ecke	tx		ty		λ	φ	Länge in Grad	Breite in Grad
Oben links	KachelX	4252	KachelY	2736	0.11965050	0.89324804	6.855469	51.179343
Oben rechts	KachelX + 1	4253	KachelY	2736	0.12041749	0.89324804	6.899414	51.179343
Unten links	KachelX	4252	KachelY + 1	2737	0.11965050	0.89276709	6.855469	51.151786
Unten rechts	KachelX + 1	4253	KachelY + 1	2737	0.12041749	0.89276709	6.899414	51.151786

Die Kantenlängen werden aus den Koordinaten der Eckpunkte sowie dem Erdradius R berechnet:

Name	Formel	Berechnung	Ergebnis
Linke Seite	abs(φ_OL-φ_UL) × R	abs(0.89324804-0.89276709) × R 0.000480950000000036 × 6371000	dl = 3064.13245000023m
Rechte Seite	abs(φ_OR-φ_UR) × R	abs(0.89324804-0.89276709) × R 0.000480950000000036 × 6371000	dr = 3064.13245000023m
Obere Seite	abs(λ_OL-λ_OR) × cos(φ_OL) × R	abs(0.11965050-0.12041749) × cos(0.89324804) × R 0.000766989999999995 × 0.626884750821521 × 6371000	do = 3063.26812849266m
Untere Seite	abs(λ_UL-λ_UR) × cos(φ_UL) × R	abs(0.11965050-0.12041749) × cos(0.89276709) × R 0.000766989999999995 × 0.627259392220409 × 6371000	du = 3065.09881117449m

Die Fläche wird aus den Koordinaten von linker oberer und rechter unterer Ecke sowie dem Erdradius R berechnet:

Name	Formel	Berechnung	Ergebnis
Fläche	abs(λ₁-λ₂)×abs(sinφ₁-sinφ₂)×R²	abs(λ₁-λ₂)×abs(sin(0.89324804)-sin(0.89276709))×R² abs(λ₁-λ₂)×abs(0.626884750821521-0.627259392220409)×R² abs(0.12041749-0.11965050)×0.00037464139888832×R² 0.000766989999999995×0.00037464139888832×6371000² 0.000766989999999995×0.00037464139888832×40589641000000	ar = 9389064.18365333m²