↖ | ↑ | ↗ | ||
← | N 79 |
← 56.41 m → | N 79 |
→ |
↑ 56.38 m ↓ |
↑ 56.38 m ↓ |
|||
N 79 |
← 56.41 m → 3 180 m² |
N 79 |
||
↙ | ↓ | ↘ |
Die Berechnung
-
Aus der URL werden die Kachel-Parameter übernommen:
Name Parameter Wert Wertebereich Erklärung zum Wertebereich Vergrößerungsstufe tz
17 0… Keine obere Grenze. Es sollte aber Kacheln für die gewählte Stufe geben. Kachel-X tx
16570 0…2zoom-1 Kachel-Y ty
16020 0…2zoom-1 -
Aus der Kartenposition x=0.126422882080078 y=0.122226715087891 und der
Vergrößerungsstufe zoom=17 berechnen wir die Kachelnummer:
Name Formel Berechnung Ergebnis Kachel-X (tx) floor (x × 2zoom) floor (0.126422882080078 × 217)
floor (0.126422882080078 × 131072)
floor (16570.5)tx = 16570 Kachel-Y (ty) floor (y × 2zoom) floor (0.122226715087891 × 217)
floor (0.122226715087891 × 131072)
floor (16020.5)ty = 16020 Kachel-Pfad (ti) "zoom/tx/tz" 17 / 16570 / 16020 ti = "17/16570/16020" -
Anzeige der Kachel
und Kaffeepause.https://a.tile.openstreetmap.org/17/16570/16020.png
-
Wir berechnen aus der Kachelnummer die Koordinaten der linken oberen Kachelecke.
Name Formel Berechnung Ergebnis X-Position (x) tx ÷ 2tz 16570 ÷ 217
16570 ÷ 131072x = 0.126419067382812 Y-Position (y) ty ÷ 2tz 16020 ÷ 217
16020 ÷ 131072y = 0.122222900390625 Länge (Λ)
(Merkator)+(x × 2 - 1) × π +(0.126419067382812 × 2 - 1) × π
-0.747161865234375 × 3.1415926535Λ = -2.34727823 Breite (Φ)
(Merkator)-(y × 2 - 1) × π -(0.122222900390625 × 2 - 1) × π
0.75555419921875 × 3.1415926535Φ = 2.3736435215867 Länge (λ) Λ (unverändert) -2.34727823} λ = -2.34727823} Breite (φ) 2×atan(exp(Φ))-π/2 2×atan(exp(2.3736435215867))-π/2
2×atan(10.7364395554107)-π/2
2×1.47792352467942-π/2
2.95584704935883-1.57079632675φ = 1.38505072 Länge in Grad λ ÷ π × 180° -2.34727823} ÷ 3.1415926535 × 180° lon = -134.489136° Breite in Grad φ ÷ π × 180° 1.38505072 ÷ 3.1415926535 × 180° lat = 79.357561° -
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 16570 KachelY 16020 -2.34727823 1.38505072 -134.489136 79.357561 Oben rechts KachelX + 1 16571 KachelY 16020 -2.34723029 1.38505072 -134.486389 79.357561 Unten links KachelX 16570 KachelY + 1 16021 -2.34727823 1.38504187 -134.489136 79.357054 Unten rechts KachelX + 1 16571 KachelY + 1 16021 -2.34723029 1.38504187 -134.486389 79.357054 -
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(1.38505072-1.38504187) × R
8.84999999994918e-06 × 6371000dl = 56.3833499996762m Rechte Seite abs(φOR-φUR) × R abs(1.38505072-1.38504187) × R
8.84999999994918e-06 × 6371000dr = 56.3833499996762m Obere Seite abs(λOL-λOR) × cos(φOL) × R abs(-2.34727823--2.34723029) × cos(1.38505072) × R
4.79399999999686e-05 × 0.184679366276907 × 6371000do = 56.4058321078183m Untere Seite abs(λUL-λUR) × cos(φUL) × R abs(-2.34727823--2.34723029) × cos(1.38504187) × R
4.79399999999686e-05 × 0.184688064039281 × 6371000du = 56.4084886283279m -
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(λ1-λ2)× abs(sinφ1-sinφ2)× R² abs(λ1-λ2)× abs(sin(1.38505072)-sin(1.38504187))× R²
abs(λ1-λ2)×abs(0.184679366276907-0.184688064039281)× R²
abs(-2.34723029--2.34727823)×8.69776237444331e-06× R²
4.79399999999686e-05×8.69776237444331e-06× 6371000²
4.79399999999686e-05×8.69776237444331e-06× 40589641000000 ar = 3180.42466541328m²