449 lines
10 KiB
Python
449 lines
10 KiB
Python
from __future__ import absolute_import
|
|
from __future__ import division
|
|
from __future__ import print_function
|
|
from __future__ import unicode_literals
|
|
|
|
from __future__ import absolute_import
|
|
from rowers.utils import lbstoN
|
|
import numpy as np
|
|
|
|
import pandas as pd
|
|
from scipy import optimize
|
|
from django.utils import timezone
|
|
|
|
from math import log10
|
|
|
|
nometrics = [
|
|
'originalvelo',
|
|
# 'cumdist',
|
|
'strokes_slsh_min',
|
|
' WorkPerStroke (joules)',
|
|
' activityIdx',
|
|
' lapIdx',
|
|
# ' pointIdx',
|
|
' WorkoutType',
|
|
' IntervalType',
|
|
' WorkoutState',
|
|
' RowingState',
|
|
' WorkoutDurationType',
|
|
' WorkoutIntervalCount',
|
|
'ergpace',
|
|
'ref',
|
|
'id',
|
|
'deltat',
|
|
'workoutid',
|
|
'totalangle',
|
|
'hr_bottom',
|
|
'Position',
|
|
'Extensions',
|
|
'GPS Speed',
|
|
'Split (IMP)',
|
|
'Speed (IMP)',
|
|
'driveenergy',
|
|
'Distance/Stroke (IMP)',
|
|
'Distance/Stroke (GPS)',
|
|
'Distance (IMP)',
|
|
'equivergpower',
|
|
'cum_dist.1'
|
|
]
|
|
|
|
rowingmetrics = (
|
|
('time',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Time',
|
|
'ax_min': 0,
|
|
'ax_max': 1e5,
|
|
'mode':'both',
|
|
'type': 'basic',
|
|
'group':'basic',
|
|
'sigfigs': -1}),
|
|
|
|
('hr',{
|
|
'numtype':'integer',
|
|
'null':True,
|
|
'verbose_name': 'Heart Rate (bpm)',
|
|
'ax_min': 100,
|
|
'ax_max': 200,
|
|
'mode':'both',
|
|
'type': 'basic',
|
|
'group':'athlete',
|
|
'sigfigs':0,}),
|
|
|
|
('pace',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Pace (/500m)',
|
|
'ax_min': 1.0e3*210,
|
|
'ax_max': 1.0e3*75,
|
|
'mode':'both',
|
|
'type': 'basic',
|
|
'sigfigs': 1,
|
|
'group': 'basic'}),
|
|
|
|
('velo',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Boat Speed (m/s)',
|
|
'ax_min': 0,
|
|
'ax_max': 8,
|
|
'default': 0,
|
|
'mode':'both',
|
|
'sigfigs': 1,
|
|
'type':'pro',
|
|
'group': 'basic'}),
|
|
|
|
# ('powerhr',{
|
|
# 'numtype':'float',
|
|
# 'null':True,
|
|
# 'verbose_name': 'Power Heart Rate Efficiency (J/beat)',
|
|
# 'ax_min':0,
|
|
# 'ax_max':300,
|
|
# 'mode':'both',
|
|
# 'type':'pro',
|
|
# 'group':'athlete'}),
|
|
|
|
('spm',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Stroke Rate (spm)',
|
|
'ax_min': 15,
|
|
'ax_max': 45,
|
|
'mode':'both',
|
|
'sigfigs': 1,
|
|
'type': 'basic',
|
|
'group':'basic'}),
|
|
|
|
('driveenergy',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Work Per Stroke (J)',
|
|
'ax_min': 0,
|
|
'ax_max': 1000,
|
|
'mode':'both',
|
|
'sigfigs': 0,
|
|
'type': 'pro',
|
|
'group':'forcepower'}),
|
|
|
|
('power',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Power (W)',
|
|
'ax_min': 0,
|
|
'ax_max': 600,
|
|
'mode':'both',
|
|
'sigfigs': 0,
|
|
'type': 'basic',
|
|
'group':'forcepower'}),
|
|
|
|
('averageforce',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Average Drive Force (N)',
|
|
'ax_min': 0,
|
|
'ax_max': 1200,
|
|
'mode':'both',
|
|
'sigfigs': 0,
|
|
'type': 'pro',
|
|
'group':'forcepower'}),
|
|
|
|
('peakforce',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Peak Drive Force (N)',
|
|
'ax_min': 0,
|
|
'ax_max': 1500,
|
|
'sigfigs': 0,
|
|
'mode':'both',
|
|
'type': 'pro',
|
|
'group':'forcepower'}),
|
|
|
|
('drivelength',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Drive Length (m)',
|
|
'ax_min': 0,
|
|
'ax_max': 2.5,
|
|
'mode':'rower',
|
|
'sigfigs': 2,
|
|
'type': 'pro',
|
|
'group':'stroke'}),
|
|
|
|
('forceratio',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Average/Peak Force Ratio',
|
|
'ax_min': 0,
|
|
'ax_max': 1,
|
|
'sigfigs': 2,
|
|
'mode':'both',
|
|
'type': 'pro',
|
|
'group': 'forcepower'}),
|
|
|
|
('distance',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Interval Distance (m)',
|
|
'ax_min': 0,
|
|
'ax_max': 1e5,
|
|
'sigfigs': 0,
|
|
'mode':'both',
|
|
'type': 'basic',
|
|
'group':'basic'}),
|
|
|
|
('cumdist',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Cumulative Distance (m)',
|
|
'ax_min': 0,
|
|
'ax_max': 1e5,
|
|
'mode':'both',
|
|
'sigfigs': 0,
|
|
'type': 'basic',
|
|
'group':'basic'}),
|
|
|
|
('drivespeed',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Drive Speed (m/s)',
|
|
'ax_min': 0,
|
|
'ax_max': 4,
|
|
'default': 0,
|
|
'sigfigs': 2,
|
|
'mode':'both',
|
|
'type': 'pro',
|
|
'group': 'stroke'}),
|
|
|
|
|
|
('catch',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Catch Angle (degrees)',
|
|
'ax_min': -40,
|
|
'ax_max': -75,
|
|
'default': 0,
|
|
'sigfigs': 0,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group': 'stroke'}),
|
|
|
|
('slip',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Slip (degrees)',
|
|
'ax_min': 0,
|
|
'ax_max': 20,
|
|
'default': 0,
|
|
'sigfigs': 1,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group': 'stroke'}),
|
|
|
|
('finish',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Finish Angle (degrees)',
|
|
'ax_min': 20,
|
|
'ax_max': 55,
|
|
'default': 0,
|
|
'sigfigs': 0,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group': 'stroke'}),
|
|
|
|
('wash',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Wash (degrees)',
|
|
'ax_min': 0,
|
|
'ax_max': 30,
|
|
'default': 0,
|
|
'sigfigs': 1,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group': 'stroke'}),
|
|
|
|
('peakforceangle',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Peak Force Angle',
|
|
'ax_min': -50,
|
|
'ax_max': 50,
|
|
'default': 0,
|
|
'sigfigs': 0,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group':'stroke'}),
|
|
|
|
|
|
('totalangle',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Drive Length (deg)',
|
|
'ax_min': 40,
|
|
'ax_max': 140,
|
|
'default': 0,
|
|
'sigfigs': 0,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group':'stroke'}),
|
|
|
|
|
|
('effectiveangle',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Effective Drive Length (deg)',
|
|
'ax_min': 40,
|
|
'ax_max': 140,
|
|
'default': 0,
|
|
'sigfigs': 0,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group':'stroke'}),
|
|
|
|
('rhythm',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Stroke Rhythm',
|
|
'ax_min': 20,
|
|
'ax_max': 55,
|
|
'default': 1.0,
|
|
'sigfigs': 0,
|
|
'mode':'erg',
|
|
'type': 'pro',
|
|
'group':'stroke'}),
|
|
|
|
('efficiency',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'OTW Efficiency (%)',
|
|
'ax_min': 0,
|
|
'ax_max': 110,
|
|
'default': 0,
|
|
'sigfigs': 0,
|
|
'mode':'water',
|
|
'type': 'pro',
|
|
'group':'forcepower'}),
|
|
|
|
('distanceperstroke',{
|
|
'numtype':'float',
|
|
'null':True,
|
|
'verbose_name': 'Distance per Stroke (m)',
|
|
'ax_min': 0,
|
|
'ax_max': 15,
|
|
'default': 0,
|
|
'sigfigs': 1,
|
|
'mode':'both',
|
|
'type': 'basic',
|
|
'group':'basic'}),
|
|
|
|
)
|
|
|
|
metricsgroups = list(set([d['group'] for n,d in rowingmetrics]))
|
|
|
|
dtypes = {}
|
|
|
|
for name,d in rowingmetrics:
|
|
if d['numtype'] == 'float':
|
|
dtypes[name] = float
|
|
elif d['numtype'] == 'int':
|
|
dtypes[name] = int
|
|
|
|
axesnew = [
|
|
(name,d['verbose_name'],d['ax_min'],d['ax_max'],d['type']) for name,d in rowingmetrics
|
|
]
|
|
|
|
axes = tuple(axesnew+[('None','<None>',0,1,'basic')])
|
|
|
|
axlabels = {ax[0]:ax[1] for ax in axes}
|
|
|
|
yaxminima = {ax[0]:ax[2] for ax in axes}
|
|
|
|
yaxmaxima = {ax[0]:ax[3] for ax in axes}
|
|
|
|
defaultfavoritecharts = (
|
|
{
|
|
'yparam1':'pace',
|
|
'yparam2':'spm',
|
|
'xparam':'time',
|
|
'plottype':'line',
|
|
'workouttype':'both',
|
|
'reststrokes':True,
|
|
'notes':"""This chart shows your pace and stroke rate versus time. """,
|
|
},
|
|
{
|
|
'yparam1':'pace',
|
|
'yparam2':'hr',
|
|
'xparam':'time',
|
|
'plottype':'line',
|
|
'workouttype':'both',
|
|
'reststrokes':True,
|
|
'notes':"""This chart shows your pace and heart rate versus time.
|
|
Heart rate values will be shown only when it is in your data, i.e.
|
|
in case you recorded your heart rate during your workout""",
|
|
},
|
|
{
|
|
'yparam1':'distanceperstroke',
|
|
'yparam2':'hr',
|
|
'xparam':'time',
|
|
'plottype':'line',
|
|
'workouttype':'otw',
|
|
'reststrokes':True,
|
|
'notes':"""This chart shows the Distance covered per stroke, and your
|
|
heart rate versus time. """,
|
|
},
|
|
{
|
|
'yparam1':'driveenergy',
|
|
'yparam2':'hr',
|
|
'xparam':'time',
|
|
'plottype':'line',
|
|
'workouttype':'ote',
|
|
'reststrokes':True,
|
|
'notes':"""This chart shows the Work per Stroke and your heart rate
|
|
plotted versus time. """,
|
|
},
|
|
{
|
|
'yparam1':'distanceperstroke',
|
|
'yparam2':'None',
|
|
'xparam':'spm',
|
|
'plottype':'scatter',
|
|
'workouttype':'otw',
|
|
'reststrokes':True,
|
|
'notes':"""This chart shows the Distance per Stroke versus stroke
|
|
stroke rate. You should see a steady decline of the Distance per Stroke
|
|
as you increase stroke rate. Typical values are > 10m for steady state
|
|
dropping to 8m for race pace in the single.""",
|
|
},
|
|
{
|
|
'yparam1':'driveenergy',
|
|
'yparam2':'None',
|
|
'xparam':'spm',
|
|
'plottype':'scatter',
|
|
'workouttype':'ote',
|
|
'reststrokes':True,
|
|
'notes':"""This chart shows the Work per Stroke versus Stroke Rate.
|
|
This value should be fairly constant across all stroke rates.""",
|
|
},
|
|
)
|
|
|
|
|
|
|
|
def calc_trimp(df,sex,hrmax,hrmin,hrftp):
|
|
if sex == 'male':
|
|
f = 1.92
|
|
else:
|
|
f = 1.67
|
|
|
|
dt = df['time'].diff()/6.e4
|
|
|
|
hrr = (df['hr']-hrmin)/(hrmax-hrmin)
|
|
hrrftp = (hrftp-hrmin)/float(hrmax-hrmin)
|
|
trimp1hr = 60*hrrftp*0.64*np.exp(f*hrrftp)
|
|
|
|
trimpdata = dt*hrr*0.64*np.exp(f*hrr)
|
|
trimp = trimpdata.sum()
|
|
|
|
hrtss = 100*trimp/trimp1hr
|
|
|
|
return trimp,hrtss
|