diff options
Diffstat (limited to 'articles/2011-06-14_darstellen_von_gps_daten_mit_gnuplot.md')
| -rw-r--r-- | articles/2011-06-14_darstellen_von_gps_daten_mit_gnuplot.md | 32 |
1 files changed, 13 insertions, 19 deletions
diff --git a/articles/2011-06-14_darstellen_von_gps_daten_mit_gnuplot.md b/articles/2011-06-14_darstellen_von_gps_daten_mit_gnuplot.md index cc97add..28c95de 100644 --- a/articles/2011-06-14_darstellen_von_gps_daten_mit_gnuplot.md +++ b/articles/2011-06-14_darstellen_von_gps_daten_mit_gnuplot.md @@ -4,19 +4,18 @@ Bei meiner letzten Wanderung in den Schweizer Alpen habe ich spaßeshalber das N Die Daten der Messpunkte sind im XML als `trkpt`-Tags gespeichert. Enthalten sind jeweils der Längen- und Breitengrad, die Uhrzeit, der Modus (3d / 2d), die Höhe über Null und die Anzahl der zur Positionsbestimmung genutzten Satelliten. Aussehen tut das ganze dann z.B. so: -~~~ +```xml <trkpt lat="47.320591" lon="9.329439"> <time>2011-06-12T07:57:39Z</time> <fix>3d</fix> <ele>870</ele> <sat>6</sat> </trkpt> -~~~ -{: .language-xml} +``` Diese Daten lassen sich nun sehr einfach Verarbeiten – ich habe das Python `xml.dom.minidom` Modul verwendet. Um die Positionen einfacher verwenden zu können, werden sie mit dieser Funktion in Listenform gebracht: -~~~ +```python def getPositions(xml): doc = minidom.parse(xml) node = doc.documentElement @@ -29,12 +28,11 @@ def getPositions(xml): pos["ele"] = int(TrkPt.getElementsByTagName("ele")[0].childNodes[0].nodeValue) positions.append(pos) return positions -~~~ -{: .language-python} +``` Aus dieser Liste kann ich jetzt schon einige Kennzahlen ziehen: -~~~ +```python def printStats(gpxPositions): highEle = gpxPositions[0]["ele"] lowEle = gpxPositions[0]["ele"] @@ -48,8 +46,7 @@ def printStats(gpxPositions): print "Lowest elevation: " + str(lowEle) print "Highest elevation: " + str(highEle) print "Height difference: " + str(eleDiv) -~~~ -{: .language-python} +``` Die Kennzahlen für meine Testdaten wären: @@ -61,21 +58,20 @@ Die Kennzahlen für meine Testdaten wären: Da die Daten ja, wie schon im Titel angekündigt, mit gnuplot dargestellt werden sollen werden sie mit dieser Funktion in für gnuplot lesbares CSV gebracht: -~~~ +```python def printCsv(gpxPositions): separator = ';' for pos in gpxPositions: print pos["lat"] + separator + pos["lon"] + separator + str(pos["ele"]) -~~~ -{: .language-python} +``` ## Plotten mit gnuplot -{: .full .clear} +{.full .clear} Eine solche, dreidimensionale Ausgabe der GPS Daten zu erzeugen ist mit der `splot`-Funktion sehr einfach. -~~~ +```sh #!/usr/bin/gnuplot set terminal png size 1280,1024 set output "output.png" @@ -90,16 +86,14 @@ set datafile separator ';' splot "/home/adrian/projects/gpxplot/wanderung_120611.csv" with impulses lt 3 lw 1 splot "/home/adrian/projects/gpxplot/wanderung_120611.csv" with lines lw 2 unset multiplot -~~~ -{: .language-sh} +``` Mit `set terminal png size 1280,1024` und `set output "output.png"` werden zuerst das Ausgabemedium, die Größe und der Dateiname der Ausgabe definiert. Dannach aktiviert `set multiplot` den gnuplot-Modus, bei dem mehrere Plots in einer Ausgabe angezeigt werden können. Dieses Verhalten brauchen wir hier, um sowohl die Strecke selbst als rote Line, als auch die zur Verdeutlichung verwendeten blauen Linien gleichzeitig anzuzeigen. Mit `set [y,x,z]range` werden die Außengrenzen des zu plottenden Bereichs gesetzt. Dies ließe sich natürlich auch über ein Script automatisch erledigen. Als Nächstes wird mit `set view 28,272,1,1` die Blickrichtung und Skalierung definiert. `set ticslevel 0` sorgt dafür, dass die Z-Achse direkt auf der Grundebene beginnt. Um ein Gitter auf der Grundfläche anzuzeigen, gibt es `set grid`. Als letztes werden jetzt die zwei Plots mit `splot` gezeichnet. Die Angaben hinter `with` steuern hierbei das Aussehen der Linien. -Falls jemand den Artikel mit meinen Daten nachvollziehen möchte - das GPX-File kann hier heruntergeladen werden: - [2011-06-12.gpx](https://static.kummerlaender.eu/media/2011-06-12.gpx) +Falls jemand den Artikel mit meinen Daten nachvollziehen möchte - das GPX-File kann hier heruntergeladen werden: [2011-06-12.gpx](https://static.kummerlaender.eu/media/2011-06-12.gpx) Zum Schluss hier noch ein Blick vom Weg auf den Kronberg Richtung Jakobsbad im Appenzell: -{: .full} +{.full} |