Sports & Health Lab · Automat Labs

Performance,
quantified.

Four signals of the body — motion, sleep, temperature and biometrics — measured with sensors, pulled together by automated pipelines, and analysed with the same metrics used to score surgical skill. The goal: know what actually makes you faster and healthier, and how much of it is the athlete, the equipment or the conditions.

Motion Sleep Temperature Biometrics

01Motion analysis

Motion, scored like surgical skill.

In a study published with KU Leuven in Scientific Reports (Nature Portfolio), we tracked a fetal surgeon’s instruments in 6D and computed 39 motion metrics to test whether 3D vision improves skill. The same mathematics reads a badminton stroke, a running stride or a pedal stroke: how far, how fast, how smooth, how consistent.

Origin · fetal surgery

Method published in Scientific Reports (Nature Portfolio) with KU Leuven, 2023 — 3D vs. 2D simulated fetoscopy for spina bifida repair: a quantitative motion analysis. Ahmad, Weiler, Joyeux et al.

Read on nature.com → · doi:10.1038/s41598-023-47531-9

1.1

Three sports, three kinds of load.

Where the sensors go, what counts as one movement, and which metrics matter most.

01 · Segment · stroke

Badminton

Explosive, intermittent, and decided by technique. Measure the racket and the feet, not just the heart.

IMU in the racket handle and on the wrist. The non-racket arm plays the part of the study’s assistant hand: it stabilises, and its motion counts too.

Measure

  • Racket IMU — swing speed, stroke type & count
  • Insole pressure — footwork, split-step timing
  • Heart rate — rally load & recovery

Key metrics

  • Smax at contact · stroke time T
  • LDLJ & SPARC per stroke type
  • NSP — hesitation in the swing
  • SPARC on ω(t) — pronation
02 · Segment · stride

Running — trail & marathon

Long efforts where pacing, economy and heat decide the result. Trail adds climbing and technical ground.

IMU on the shoe or ankle and at the sacrum. One stride from foot strike to foot strike; a marathon is thousands of repetitions of the same segment.

Measure

  • HR / HRV · pace · cadence · ground contact
  • Running power · VAM on climbs (trail)
  • Pacing drift & fuelling (marathon)

Key metrics

  • σS stride to stride
  • LDLJ drift with fatigue
  • D along the vertical — oscillation
  • TOR on technical trail descents
03 · Segment · revolution

Cycling — gravel, MTB & road

Power is the headline number, but the terrain decides how much of it reaches the road.

IMU on the crank or pedal and at the knee. One crank revolution per segment, for gravel, MTB and road alike.

Measure

  • Power · cadence · HR · speed
  • Vibration & suspension travel (gravel, MTB)
  • Aero position & pacing (road)

Key metrics

  • σS of pedal speed — uniformity
  • NAP — dead spots
  • κ & TOR of the knee path
  • Bike sway out of the saddle
1.2

A movement is a curve in space.

Every sample is a point p(t). Its path, bounding volume and speed along the way are the raw material of every metric below.

Δx 1552 mmΔy 250 mmΔz 828 mmp(t₀) · startSmax 11.6 m/sRACKET · HEAD IMUsimulated · 200 Hz · T = 0.45 s011.6 m/sspeed v(t)bounding box · P = 2.45 m

Badminton · overhead clear

IMU at the racket head. One stroke from backswing to follow-through: the arc sets the path length, the speed peaks at contact.

colour
speed v(t), blue → ochre
box
extent Δx · Δy · Δz
shadow
floor projection
data
simulated, seeded, 200 Hz
Simulated signal — metrics are computed live
1.3

From raw sensor to clean segments.

The preprocessing pipeline of the study, unchanged. Only the reference frames and the segment boundaries change with the sport.

a

Format

Raw logs — EM-tracker files in the study, FIT/IMU logs here — into one table per sensor on one clock.

b

Transform

All sensors into one reference frame: the operating-table corner there; court, pelvis or bike frame here.

c

Prune

Gaps > 1 s removed, shorter gaps bridged with a quadratic spline, then report how much is real data:

vi≜Tmotion,iTi
d

Smooth

Butterworth low-pass, fc = 10 Hz, n = 5, then an exponential filter:

yk=αxk+(1−α)yk−1,α=0.05
e

Annotate & debug

Segments marked from video: surgical steps there; strokes, strides, revolutions here. Interpolated stretches beyond physical limits (6 m/s² at a surgeon’s hand) are rejected; the limit is set per sensor.

Source: Scientific Reports (Nature Portfolio) 13:20951 (2023) — Methods · Motion data processing.

1.4

Position, then three derivatives.

Speed, acceleration and jerk are the first, second and third derivative of position. Smoothness metrics live in the third.

position · x y z 𝒑(t)=(xyz)⊤,p=x2+y2+z2
speed 𝒗(t)=d𝒑dt,v=x˙2+y˙2+z˙2
acceleration 𝒂(t)=d2𝒑dt2,a=‖𝒂(t)‖
jerk 𝒋(t)=d3𝒑dt3,j=‖𝒋(t)‖
1.5

Metric bench: smooth vs intermittent.

Slide from one fluent movement to one built from hesitant sub-movements with tremor. The speed peaks, the spectrum and the metrics respond as they did between 2D and 3D surgery.

k = 0.40 Simulated
αS̄ · NSP(α/100) Smax · MAPR13.09.76.53.20.00.0 s0.1 s0.2 s0.3 s0.5 sv(t) · m/s
V̄ = 0.05ωc = 6.5 Hz0246810 Hz1.000.750.500.250.00V̂(ω) / V̂(0)SPARC = −1.53
P2.45mvs smooth ±0 %
Smax11.57m/svs smooth −5 %
S̄5.39m/svs smooth +1 %
σS3.73m/svs smooth −15 %
NSPα=12vs smooth +100 %
MAPRα=1077%vs smooth +8 %
LDLJ7.68vs smooth +3 %
SPARC−1.53vs smooth −9 %
TOR1.000vs smooth ±0 %

NSP counts speed maxima above α·S̄ (α = 1). MAPR: share of time above 10 % of Smax. LDLJ as in the paper, ln|DLJ| — higher means jerkier. SPARC closer to 0 means smoother. Change shown vs the smooth (k = 0) movement.

1.6

The metric atlas.

Every definition used in the study (Scientific Reports (Nature Portfolio) 13:20951, Table 1), plus curvature, turning angle and tortuosity. Index i is a movement segment — a surgical step there, a stroke, stride or crank revolution here.

T · P · V · D · TOR

Time & space

How long, how far, how much room. → stroke arc · stride loop · pedal circle, and how much the envelope wanders.

Total time
Ti≜∫ti,startti,enddt

Shorter in experts; the primary outcome of the study.

Path length
Pi≜∫0Tivi(t)dt

Distance travelled. Significant 3D vs 2D.

Volume
Vi=convhull(𝒑i(t))

Convex hull of all positions: the space a movement claims.

Depth perception
Di≜∫0Ti(dz′dt)2dt

Path along the tool axis z′ — here the racket shaft or the tibia.

Tortuosity
TORi=1+a,a=−dlnL(s)dlns

Divider method: path length L vs divider step s. 1 = straight, 2 = maximally convoluted.

S · NSP · SPR · MAPR

Speed

How fast, how even, how often it surges. → racket speed at contact · ankle swing speed · crank uniformity (σ → 0 = even pedalling).

Max · mean · consistency
Si,max≜max(vi),S¯i≜mean(vi),σSi≜std(vi)

Peak, average and spread of the speed profile.

Number of speed peaks
NSPi,α≜|{vi(t)∣dvidt=0,d2vidt2<0,vi>αS¯i}|

Each peak is a sub-movement or a correction. Significant 3D vs 2D.

Speed peak rate
SPRα,i=NSPi,αTi

Peaks per second, less dependent on duration.

Movement arrest period ratio
MAPRi,α≜Ti,vi>(α/100)Si,maxTi

Share of time actually moving vs pausing.

A · NAP · APR · IAV

Acceleration

Force in disguise: starts, stops, impacts. → lunges and split-steps · foot strike · dead spots at 12 and 6 o’clock.

Max · mean · consistency
Ai,max≜max(ai),A¯i≜mean(ai),σAi≜std(ai)

Peak, average and spread of acceleration.

Number of acceleration peaks
NAPi,α≜|{ai(t)∣daidt=0,d2aidt2<0,ai>αA¯i}|

Bursts of force. Significant 3D vs 2D.

Acceleration peak rate
APRi,α=NAPi,αTi

Bursts per second.

Integral of acceleration
IAVi≜∫0Tiai(t)dt

Total “effort” of changing velocity.

Sm · SAL · SPARC · LDLJ

Smoothness

Where skill shows: the third derivative and the spectrum. → the fluency of a stroke · form breaking down in the last 10 km · a round pedal stroke.

Smoothness
Smi≜1Ti12∫0Ti[(d3xdt3)2+(d3ydt3)2+(d3zdt3)2]dt

Root jerk; not dimensionless, so it scales with amplitude and duration.

Spectral arc length
SALi≜−∫0ωc(1ωc)2+(dV^i(ω)dω)2dω,V^i(ω)=Vi(ω)Vi(0)

Arc length of the normalised speed spectrum: a smooth movement has a short, simple spectrum.

SPARC · adaptive cut-off
ωc≜min{ωcmax,min{ω∣V^i(r)<V¯,∀r>ω}}

SAL with a cut-off chosen per movement; independent of temporal scaling. Closer to 0 = smoother.

Dimensionless · log jerk
DLJi≜(Ti)3vi,peak2∫0Tiji(t)2dt LDLJi≜ln|DLJi|

Jerk made dimensionless. Significant 3D vs 2D.

κ · TA

Shape

How straight, how sharply it turns. → the arc of a smash vs a drive · knee tracking in the frontal plane · cornering lines.

Curvature
κi(t)≜‖𝒗i(t)×𝒂i(t)‖vi3(t)

≈ 0 for straight movements; summarised by mean, max or spread.

Turning angle
TAi(t)≜arccos(𝒖i,t−1·𝒖i,t+1‖𝒖i,t−1‖‖𝒖i,t+1‖)

With ut−1 = p(t) − p(t−1): angle between consecutive steps.

q · ω

Orientation

The study’s open lead: the tracker also logs rotation, and almost no study uses it. → forearm pronation and racket-face angle · foot-strike angle · ankling through the pedal stroke.

Quaternion
𝒒(t)≜(q4q1q2q3)⊤,‖𝒒‖=1

Orientation without gimbal lock, straight from the sensor.

Angular speed
ω(t)=2‖𝒒˙(t)⊗𝒒*(t)‖

Feed ω(t) instead of v(t) into NSP, SPARC or LDLJ: rotation gets its own smoothness score.

Scroll or drag the strip · hover a formula for its name and meaning.

02Sleep

Recovery is where training pays off.

Sleep stages from a wrist or ring sensor, nocturnal heart-rate variability, and movement during the night — scored with the same peak detection as the motion toolbox.

Sleep efficiency
SE=TSTTIB×100%

Total sleep time over time in bed.

Wake after sleep onset
WASO=∑e>SOLΔte[se=wake]

Minutes awake after sleep-onset latency SOL; 30-second epochs e.

Night movement
NAPα(awrist),MAPRα(vwrist)

Restlessness, counted with the motion metrics above.

WakeREMN1N2N3wrist movement23:001:003:005:007:00simulated · SE = 94 %
Hypnogram · 23:00–07:00 · 30-s epochsSimulated

03Body temperature

Heat is a hidden load.

Core and skin temperature drive fatigue. How do layers, fabrics and fit change skin and core temperature, sweat rate and heart-rate drift — in heat and in cold?

Mean skin temperature · Ramanathan
T¯sk=0.3Tchest+0.3Tarm+0.2Tthigh+0.2Tleg

Four skin sensors, one weighted number.

Physiological strain index · Moran
PSI=5Tc,t−Tc,039.5−Tc,0+5HRt−HR0180−HR0

Heat strain on a 0–10 scale from core temperature and heart rate.

Cardiac drift
ΔHRdrift=HR¯2nd half−HR¯1st halfHR¯1st half|P¯=const

Heart-rate rise at constant power or pace: the price of heat.

39.538.537.536.535.534.533.532.519016514011590Tc coreT̄sk mean skinHR · heart rate°Cbpmsimulated · 60 min · PSI = 7.10102030405060 min
60-min run in the heat · core, skin, HRSimulated

04Biometrics

Heart, load and energy.

The classic numbers — heart rate, heart-rate variability, training load and energy expenditure — computed the same way every time, so sessions compare across weeks, sports and equipment.

Heart-rate reserve · Karvonen
HRtarget=HRrest+x·(HRmax−HRrest)

Zones as a fraction x of the reserve.

HRV · RMSSD
RMSSD=1N−1∑i=1N−1(RRi+1−RRi)2

Beat-to-beat variability: a recovery signal.

Training impulse · Banister
TRIMP=∑tΔt·ht·0.64e1.92ht,ht=HRt−HRrestHRmax−HRrest

Internal load, weighted towards high intensity (coefficients 0.86 and 1.67 for women).

Energy from power
Emet=1η∫0TP(t)dt,η≈0.24⇒kcal≈kJmech

Gross efficiency turns power-meter work into energy burned.

811 msR1888 msR2840 msR3788 msR4774 msR5839 msR6ECG · RRᵢsimulated · RMSSD = 56 ms
ECG · R-R intervals · RMSSDSimulated

05Equipment & environment

How much is the athlete, and how much is the kit?

The research focus of the lab: isolating the effect of equipment, materials and conditions — with the same controlled 2D-vs-3D mindset, now applied to a racket, a shoe or a tyre.

M

Materials

Rackets, strings, shoes and shuttles. What do stiffness, weight balance and foam actually change in the data — and in the smoothness of the movement?

  • Racket frame layup & stiffness
  • String tension
  • Shoe foam & stack · shuttle type
K

Mechanics & terrain

The bike as a machine: geometry, tyres and drivetrain on asphalt, gravel and trail.

  • Tyre width & pressure vs rolling resistance
  • Compliance & vibration per terrain
  • Geometry · drivetrain efficiency

06Metrics

What gets measured.

MetricPillarSensorSportsWhat it quantifies
Path, speed & peaksMotionIMUAllP, Smax, NSP per stroke, stride or revolution
SmoothnessMotionIMUAllLDLJ and SPARC — fluency, fatigue, technique
OrientationMotionIMU (quaternion)Badminton, cyclingPronation, racket face, ankling
Foot pressureMotionPressure insolesBadminton, runningFootwork, landing pattern, left/right balance
Sleep stages & efficiencySleepWrist / ringAllTST, SE, WASO, night movement
Skin & core temperatureTemperatureTemperature sensorsAllHeat strain (PSI), clothing and material effects
Heart rate & HRVBiometricsChest strapAllInternal load, TRIMP, recovery, cardiac drift
Power & energyBiometricsPower meterCycling, runningExternal work, kcal — the reference for every comparison
Pace, speed & elevationBiometricsGPS · barometerRunning, cyclingPacing, climbing rate (VAM), terrain sections
Vibration & tyre pressureEquipmentFrame accelerometer · pressureGravel, MTBTerrain roughness, rolling resistance vs grip

07Project log

Work in the open.

Planned
Motion-metrics toolbox: MATLAB → Python
motionPython
Tyre pressure vs terrain on gravel
cyclingmechanics
Clothing layers vs skin temperature
runningtemperature
Marathon pacing drift analysis
runningHR
In progress
Racket IMU — stroke detection
badmintonsensors
Automated FIT/GPX ingestion pipeline
Pythondata
Published

Train, test, compare.

Club, coach, equipment brand or curious athlete — questions, ideas or a sensor to test? Get in touch.