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();
