[Talk-de] Distanz-Berechnungen aus Koordinaten

[Florian Heer]

Hi!

Sascha Silbe schrieb:
> On Fri, Sep 19, 2008 at 03:42:11AM +0200, Florian Heer wrote:
>
>> http://www.openstreetmap.org/user/florianheer/traces/160678
> Habs gerade mal kurz von Python durchrechnen lassen (mittels 
> Vincenty), komme auf 8427.5424101584249m (lat/lon korrekt) bzw. 
> 11603.755080105428m (lat/lon vertauscht). Keine Filterung (967 Punkte, 
> 966 Segmente), keine Glättung, eingebauter Fließkommatyp (IEEE 754 
> double precision). Maximale Einzeldistanz ist ca. 37.88m, minimale >0 
> ist ca. 18.4cm.

Mit kurz-/lang-Auswertung:
Exp-080803-8,53km-.gpx is 10.421032628337267 km long (simple)
shortest: 0.14556823651036202m longest: 43.1536557160521m
Exp-080803-8,53km-.gpx is 10.432387536300029 km long (Vincenty)
shortest: 0.1456457768180979m longest: 43.155324462589306m


Lustig, hab gerade mal bei mir Lat/Lon getauscht:
Exp-080803-8,53km-.gpx is 14.463829250546146 km long (simple)
shortest: 0.18494348788838302m longest: 66.67561313341307m
Exp-080803-8,53km-.gpx is 14.39166351008208 km long (Vincenty)
shortest: 0.18420196663375732m longest: 66.26012533075907m

Da tauchen wenigstens die ca 18.4cm auf.

Nochmal zur Klärung, was ich da aufsummiere, bzw. was ich da mache, 
vielleicht ist da ja schon ein Fehler:
Parsen des Tracks auf _einfachste_ Art und Weise: alle trkpt werden mit 
Lat (double) / Lon (double) / Time (DateTime) in eine MySQL-DB 
geschrieben. Select auf alle Punkte, sortiert nach Zeit, aufsteigend.
Dann Distanz zwischen Punkt und Vorgänger berechnen und aufsummieren. 
Und nein: ich habe nicht vergessen, dass ich beim ersten Punkt noch 
keinen Vorgänger habe :-)

>> Sofern OSM da keine Glättungen macht, aber ich kanns mir nicht 
>> wirklich vorstellen...
> Da fällt mir noch was ein: Welche Datenformate hast Du benutzt und 
> wieviele Nachkommastellen wurden da jeweils gespeichert? Das GPX sieht 
> aus, als wären beim Upload nur 5 Stellen hinter dem Komma angegeben 
> worden, das ist eine Auflösung von 1.1m (Latitude) bzw. 0.7m 
> (Longitude). 94 der von mir berechneten Werte sind kleiner als 2*1.1m, 
> mit einer Summe von 91.8m. Kein großer Fehler (dürfte im Bereich um 1% 
> sein), aber kleinere Unterschiede könnte es erklären.

Das "Original" des Tracks ist NoniGPSPlot, von dort nach GPX exportiert. 
Dieses GPX wurde ungeändert nach OSM hochgeladen und bildet gleichzeitig 
die Basis für meinen lokalen Import. In dem Track sind etwas mehr 
Nachkommastellen (die auch korrekt in meine DB geparst wurden, ich habs 
jetzt extra nochmal geprüft, so langsam bin ich verwundert...), Auszug:

      <trkpt lat="50.7982019999" lon="7.3804293333">
        <ele>215.2</ele>
        <time>2008-08-03T18:44:24Z</time>
        <course>0.0</course>
        <speed>0.0</speed>
      </trkpt>
      <trkpt lat="50.7982020512" lon="7.3804462821">
        <ele>215.1</ele>
        <time>2008-08-03T18:44:24Z</time>
        <course>0.0</course>
        <speed>0.0</speed>
      </trkpt>
      <trkpt lat="50.7982020512" lon="7.3804462821">
        <ele>215.1</ele>
        <time>2008-08-03T18:44:31Z</time>
        <course>0.0</course>
        <speed>0.0</speed>
      </trkpt>


> So, jetzt haben wir die Daten mit 3 Implementierungen (Noni-GPS, 
> Deine, meine) durchrechnen lassen und kommen auf 3 unterschiedliche 
> Werte. Prima. ;)

Hehe, ich guck mir morgen mal genau an, was bei raus kommt, wenn ich 
Distanzen berechne, die ich mir bei Google Earth mit dem Lineal vorgebe.

> Bei Deinen 10.42km würde ich ja fast auf lat/lon vertauscht tippen. 
> Das GPX deckt sich recht gut mit der Karte, dürfte also korrekt sein.

Das GPX ist an sich recht okay, aber lat/lon Tausch... Selbst wenn ich 
den Test nicht jetzt gerade gemacht hätte würde es mich wundern, da, wie 
schon erwähnt, bei den meisten der Tracks die Distanzangabe von Noni mit 
meinen Berechnungen harmoniert. Aber wo Du mit Deinen Werten zumindest 
zwischen die 8.53 (Noni) und 8.3 (Kompf-upload) gekommen bist, muss wohl 
doch etwas bei mir schief liegen.

> In welcher Form liegen denn die Rohdaten vor?
GPX bzw. Nonis eigenes .ngt-Format.


Viele Grüße, Florian.