def _readPointsFromActivity(activity, workdir, kwargs): try: spatialFilterDistance = kwargs["spatialFilterDistance"] except: spatialFilterDistance = None # Parse the XML document with open(workdir + os.sep + activity["GPX File"], "r") as f: soup = BeautifulSoup(f.read()) # Read all GPS trackpoints trkpts = soup.find_all("trkpt") # Convert to GPSPoint structures points = [GPSPoint.trackpointToGPSPoint(trkpt, activity) for trkpt in trkpts] if spatialFilterDistance is not None: points = GPSPoint.removeWithinRange(points, spatialFilterDistance) return points
def processPoints(workdir, points): # Get unique activity types activities = rkIO.parseIndexFile(workdir) types = set([a['Type'] for a in activities]) types = [i for i in types] # Filter points into activity types filteredpoints = [] for t in types: filteredpoints.append([point for i,point in enumerate(points) if activities[point.activity]['Type'] == t]) points = [] for i in range(len(types)): # Remove points within a specified distance print("Filtering type '" + types[i] + "'") n0 = len(filteredpoints[i]) filteredpoints[i] = GPSPoint.removeWithinRange(filteredpoints[i], 20, verbose=True) n1 = len(filteredpoints[i]) print('Pass 1 kept ' + str(n1) + '/' + str(n0)) points.extend(filteredpoints[i]) return points
#r"C:\Users\wkulp\Documents\GitHub\RunkeeperMapper\data" print('Working dir: ' + workdir) rkIO.loadArchive(workdir, 'data' + os.sep + 'data-2013-10-09-1728.zip') # Read the master activity list activities = rkIO.parseIndexFile(workdir) # Parse each GPX file into points print('Processing GPS track files') nprocs = cpu_count() print('Using ' + str(nprocs) + " CPU cores...") #points = rkIO.sp_readAllPoints(activities, workdir, spatialFilterDistance = 20) points = rkIO.mp_readAllPoints(nprocs, activities, workdir, spatialFilterDistance = 20) GPSPoint.savePoints(workdir, points, 'rawpoints.csv') points = processPoints(workdir, points) GPSPoint.savePoints(workdir, points, 'processedpoints.csv') points = GPSPoint.loadPoints(workdir, 'processedpoints.csv') writeOutput(workdir, points) print('Done!') finally: #shutil.rmtree(workdir) pass import pdb;pdb.set_trace();