↖ | ↑ | ↗ | ||
← | N 51 |
← 191.55 m → | N 51 |
→ |
↑ 191.51 m ↓ |
↑ 191.51 m ↓ |
|||
N 51 |
← 191.56 m → 36 685 m² |
N 51 |
||
↙ | ↓ | ↘ |
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 2*atan(exp(pi())) * 180/pi() - 90;'
, mit JavaScript2*Math.atan(Math.exp(Math.PI)) * 180/Math.PI - 90
.Breite in Grad lat
51.1578° -85°…+85° Vergrößerungsstufe zoom
17 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=17 berechnen wir die Kachelnummer:
Name Formel Berechnung Ergebnis Kachel-X (tx) floor (x × 2zoom) floor (0.519070833333333 × 217)
floor (0.519070833333333 × 131072)
floor (68035.6522666667)tx = 68035 Kachel-Y (ty) floor (y × 2zoom) floor (0.334079811416708 × 217)
floor (0.334079811416708 × 131072)
floor (43788.5090420107)ty = 43788 Kachel-Pfad (ti) "zoom/tx/tz" 17 / 68035 / 43788 ti = "17/68035/43788" -
Anzeige der Kachel
und Kaffeepause.https://a.tile.openstreetmap.org/17/68035/43788.png
-
Aus den Nachkommastellen 0.652266666671494 und 0.509042010700796
sowie der Kachelgröße in Pixeln tilesize=256
berechnen wir die Position des Markers:
Name Formel Berechnung Ergebnis Marker-X (px) floor (fract (x × 2zoom) * tilesize) floor (fract (0.519070833333333 × 217) × 256)
floor (fract (0.519070833333333 × 131072) × 256)
floor (fract (68035.6522666667) × 256)
floor (0.652266666671494 × 256)
floor (166.980266667902)px = 166 Marker-Y (py) floor (fract (y × 2zoom) * tilesize) floor (fract (0.334079811416708 × 217) × 256)
floor (fract (0.334079811416708 × 131072) × 256)
floor (fract (43788.5090420107) × 256)
floor (0.509042010700796 × 256)
floor (130.314754739404)py = 130 -
Wir berechnen aus der Kachelnummer die Koordinaten der linken oberen Kachelecke.
Name Formel Berechnung Ergebnis X-Position (x) tx ÷ 2tz 68035 ÷ 217
68035 ÷ 131072x = 0.519065856933594 Y-Position (y) ty ÷ 2tz 43788 ÷ 217
43788 ÷ 131072y = 0.334075927734375 Länge (Λ)
(Merkator)+(x × 2 - 1) × π +(0.519065856933594 × 2 - 1) × π
0.0381317138671875 × 3.1415926535Λ = 0.11979431 Breite (Φ)
(Merkator)-(y × 2 - 1) × π -(0.334075927734375 × 2 - 1) × π
0.33184814453125 × 3.1415926535Φ = 1.04253169293698 Länge (λ) Λ (unverändert) 0.11979431} λ = 0.11979431} Breite (φ) 2×atan(exp(Φ))-π/2 2×atan(exp(1.04253169293698))-π/2
2×atan(2.83638879762991)-π/2
2×1.23184183951212-π/2
2.46368367902425-1.57079632675φ = 0.89288735 Länge in Grad λ ÷ π × 180° 0.11979431} ÷ 3.1415926535 × 180° lon = 6.863708° Breite in Grad φ ÷ π × 180° 0.89288735 ÷ 3.1415926535 × 180° lat = 51.158677° -
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 68035 KachelY 43788 0.11979431 0.89288735 6.863708 51.158677 Oben rechts KachelX + 1 68036 KachelY 43788 0.11984225 0.89288735 6.866455 51.158677 Unten links KachelX 68035 KachelY + 1 43789 0.11979431 0.89285729 6.863708 51.156954 Unten rechts KachelX + 1 68036 KachelY + 1 43789 0.11984225 0.89285729 6.866455 51.156954 -
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.89288735-0.89285729) × R
3.00599999999429e-05 × 6371000dl = 191.512259999636m Rechte Seite abs(φOR-φUR) × R abs(0.89288735-0.89285729) × R
3.00599999999429e-05 × 6371000dr = 191.512259999636m Obere Seite abs(λOL-λOR) × cos(φOL) × R abs(0.11979431-0.11984225) × cos(0.89288735) × R
4.79399999999963e-05 × 0.627165727944857 × 6371000do = 191.552556560182m Untere Seite abs(λUL-λUR) × cos(φUL) × R abs(0.11979431-0.11984225) × cos(0.89285729) × R
4.79399999999963e-05 × 0.62718914096981 × 6371000du = 191.559707500654m -
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(0.89288735)-sin(0.89285729))× R²
abs(λ1-λ2)×abs(0.627165727944857-0.62718914096981)× R²
abs(0.11984225-0.11979431)×2.34130249524567e-05× R²
4.79399999999963e-05×2.34130249524567e-05× 6371000²
4.79399999999963e-05×2.34130249524567e-05× 40589641000000 ar = 36685.3477647586m²