import gpx without gpsbabel about pygohome HOT 4 CLOSED

Thank you, I have actually originally started with gpxpy and utm (I'm already using utm to convert back to latlon to plot in ipyleaflet), but it was comparatively much slower (speaking about 60k+ points, which represents over 900km of tracks) when loading all raw GPX files each time:

def read_segments():
    for path in sorted(Path(".").glob("gpx/*.gpx")):
        gpx = gpxpy.parse(path.read_text())
        for track in gpx.tracks:
            for segment in track.segments:
                if len(segment.points) < 5:
                    continue
                yield [
                    (point.time_difference(segment.points[0]),
                     utm.from_latlon(point.latitude, point.longitude, 32, "U")[:2])
                    for point in segment.points
                ]
segments = list(read_segments())

UTM system is perfect for local geometry and distances, but maybe with GeoPandas the expensive conversion could be skipped.

The performance problem is also the conversion of ALL points, which could be solved by some simplifying methods of gpxpy.

Or some sort of database/cache for all converted points. And we're back at CSVs, which offer the user a stable and reviewable step in data processing.

from pygohome.

Another option is to run gpsbabel in the Jupyter notebook with Unix commands (!)

https://intro.syzygy.ca/unix-tricks/#unix-and-magic-in-a-notebook

(Would have to run with ‘wget’ or ‘git clone’ to get local ‘gpsbabel’ binary.)

from pygohome.

utm can do NumPy arrays, it is just too shy about it. That's why I added it to its docs: Turbo87/utm#50
Converting 1M points in a loop takes 90 seconds, a 1M numpy array needs 0.26s. Now it's an option again.

from pygohome.

gpxpy and utm used to load GPX files in 1a30a69
Thanks @philshem

from pygohome.