nmelone/GoldenCheetah
1
0
Fork 0
mirror of https://github.com/GoldenCheetah/GoldenCheetah.git synced 2026-02-15 17:09:56 +00:00
Code Issues Packages Projects Releases Wiki Activity
Files
a38696e25d3b0756aceab41b406d6b3718bcdb76
GoldenCheetah/src/FileIO/FixDerivePower.cpp
Mark Liversedge a38696e25d Data Processor 'On Save' and Snippets 
.. a data processor can now be configured to run on save

.. in addition, the data processor is now passed the operation
   that is triggering it: "Manual", "ADD", "UPDATE", "DELETE",
   "IMPORT".

.. a new data processor has been added to write a JSON file to
   the snippet directory (added to the athlete structure).

.. this is so we can output a snippet every time an activity is
   added, updated or deleted (i.e. runs "on save").

.. these JSON snippets can be used to update external datastores
   where coaches use other analytic software but don't want to
   manually sync GC changes with there external stores.
2016-11-12 22:36:44 +00:00

294 lines
12 KiB
C++
Raw Blame History

/*
* Copyright (c) 2010 Mark Liversedge (liversedge@gmail.com)
*
* This program is free software; you can redistribute it and/or modify it
* under the terms of the GNU General Public License as published by the Free
* Software Foundation; either version 2 of the License, or (at your option)
* any later version.
*
* This program is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for
* more details.
*
* You should have received a copy of the GNU General Public License along
* with this program; if not, write to the Free Software Foundation, Inc., 51
* Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include "DataProcessor.h"
#include "LTMOutliers.h"
#include "Settings.h"
#include "Units.h"
#include "HelpWhatsThis.h"
#include <algorithm>
#include <QVector>
#ifndef MATHCONST_PI
#define MATHCONST_PI 3.141592653589793238462643383279502884L /* pi */
#endif
// Config widget used by the Preferences/Options config panes
class FixDerivePower;
class FixDerivePowerConfig : public DataProcessorConfig
{
Q_DECLARE_TR_FUNCTIONS(FixDerivePowerConfig)
friend class ::FixDerivePower;
protected:
QVBoxLayout *layoutV;
QHBoxLayout *layout;
QHBoxLayout *windLayout;
QLabel *bikeWeightLabel;
QDoubleSpinBox *bikeWeight;
QLabel *crrLabel;
QDoubleSpinBox *crr;
QLabel *windSpeedLabel;
QDoubleSpinBox *windSpeed;
QLabel *windHeadingLabel;
QDoubleSpinBox *windHeading;
public:
FixDerivePowerConfig(QWidget *parent) : DataProcessorConfig(parent) {
HelpWhatsThis *help = new HelpWhatsThis(parent);
parent->setWhatsThis(help->getWhatsThisText(HelpWhatsThis::MenuBar_Edit_EstimatePowerValues));
layout = new QHBoxLayout(this);
layout->setContentsMargins(0,0,0,0);
setContentsMargins(0,0,0,0);
bikeWeightLabel = new QLabel(tr("Bike Weight (kg)"));
crrLabel = new QLabel(tr("Crr"));
windSpeedLabel = new QLabel(tr("Wind (kph)"));
windHeadingLabel = new QLabel(tr(", heading"));
bikeWeight = new QDoubleSpinBox();
bikeWeight->setMaximum(99.9);
bikeWeight->setMinimum(0);
bikeWeight->setSingleStep(0.1);
bikeWeight->setDecimals(1);
crr = new QDoubleSpinBox();
crr->setMaximum(0.0999);
crr->setMinimum(0);
crr->setSingleStep(0.0001);
crr->setDecimals(4);
windSpeed = new QDoubleSpinBox();
windSpeed->setMaximum(99.9);
windSpeed->setMinimum(0);
windSpeed->setSingleStep(0.1);
windSpeed->setDecimals(1);
windHeading = new QDoubleSpinBox();
windHeading->setMaximum(180.0);
windHeading->setMinimum(-179.0);
windHeading->setSingleStep(1);
windHeading->setDecimals(0);
layout->addWidget(bikeWeightLabel);
layout->addWidget(bikeWeight);
layout->addWidget(crrLabel);
layout->addWidget(crr);
layout->addWidget(windSpeedLabel);
layout->addWidget(windSpeed);
layout->addWidget(windHeadingLabel);
layout->addWidget(windHeading);
layout->addStretch();
}
//~FixDerivePowerConfig() {} // deliberately not declared since Qt will delete
// the widget and its children when the config pane is deleted
QString explain() {
return(QString(tr("Derive estimated power data based on "
"speed/elevation/weight etc\n\n"
"Bike Weight parameter is added to athlete's "
"weight to compound total mass, it should "
"include apparel, shoes, etc\n\n"
"CRR parameter is the coefficient of rolling "
"resistance, it depends on tires and surface\n"
"wind speed shall be indicated in kph\n"
"wind heading (origin) unit is degrees "
"from -179 to +180 (-90=W, 0=N, 90=E, 180=S)\n"
"Note: if the ride file already contain wind data\n"
" it will be overridden if wind is entered manually")));
}
void readConfig() {
double MBik = appsettings->value(NULL, GC_DPDP_BIKEWEIGHT, "9.5").toDouble();
bikeWeight->setValue(MBik);
double Crr = appsettings->value(NULL, GC_DPDP_CRR, "0.0031").toDouble();
crr->setValue(Crr);
windSpeed->setValue(0.0);
windHeading->setValue(0.0);
}
void saveConfig() {
appsettings->setValue(GC_DPDP_BIKEWEIGHT, bikeWeight->value());
appsettings->setValue(GC_DPDP_CRR, crr->value());
}
};
// RideFile Dataprocessor -- Derive estimated power data based on
// speed/elevation/weight etc
//
class FixDerivePower : public DataProcessor {
Q_DECLARE_TR_FUNCTIONS(FixDerivePower)
public:
FixDerivePower() {}
~FixDerivePower() {}
// the processor
bool postProcess(RideFile *, DataProcessorConfig* config, QString op);
// the config widget
DataProcessorConfig* processorConfig(QWidget *parent) {
return new FixDerivePowerConfig(parent);
}
// Localized Name
QString name() {
return (tr("Estimate Power Values"));
}
};
static bool FixDerivePowerAdded = DataProcessorFactory::instance().registerProcessor(QString("Estimate Power Values"), new FixDerivePower());
bool
FixDerivePower::postProcess(RideFile *ride, DataProcessorConfig *config=0, QString op="")
{
// get settings
double MBik; // Bike weight kg
double CrV; // Coefficient of Rolling Resistance
double windSpeed; // wind speed
double windHeading; //wind direction
if (config == NULL) { // being called automatically
MBik = appsettings->value(NULL, GC_DPDP_BIKEWEIGHT, "9.5").toDouble();
CrV = appsettings->value(NULL, GC_DPDP_CRR, "0.0031").toDouble();
windSpeed = 0.0;
windHeading = 0.0;
} else { // being called manually
MBik = ((FixDerivePowerConfig*)(config))->bikeWeight->value();
CrV = ((FixDerivePowerConfig*)(config))->crr->value();
windSpeed = ((FixDerivePowerConfig*)(config))->windSpeed->value(); // kph
windHeading = ((FixDerivePowerConfig*)(config))->windHeading->value() / 180 * MATHCONST_PI; // rad
}
// Do nothing for swims and runs
if (ride->isSwim() || ride->isRun()) return false;
// if its already there do nothing !
if (ride->areDataPresent()->watts) return false;
// no dice if we don't have alt and speed
if (!ride->areDataPresent()->alt || !ride->areDataPresent()->kph) return false;
// Power Estimation Constants (mostly constants...)
double hRider = ride->getHeight(); //Height in m
double M = ride->getWeight(); //Weight kg
double T = 15; //Temp degC in not in ride data
double W = 0; // headwind (from records or based on wind parameters entered manually)
double bearing = 0.0; //cyclist direction used to compute headwind
double cCad=.002;
double afCd = 0.62;
double afSin = 0.89;
double afCm = 1.025;
double afCdBike = 1.2;
double afCATireV = 1.1;
double afCATireH = 0.9;
double afAFrame = 0.048;
double ATire = 0.031;
double CrEff = CrV;
double adipos = sqrt(M/(hRider*750));
double CwaBike = afCdBike * (afCATireV * ATire + afCATireH * ATire + afAFrame);
//qDebug()<<"CwaBike="<<CwaBike<<", afCdBike="<<afCdBike<<", afCATireV="<<afCATireV<<", ATire="<<ATire<<", afCATireH="<<afCATireH<<", afAFrame="<<afAFrame;
// apply the change
ride->command->startLUW("Estimate Power");
if (ride->areDataPresent()->slope && ride->areDataPresent()->alt
&& ride->areDataPresent()->km) {
for (int i=0; i<ride->dataPoints().count(); i++) {
RideFilePoint *p = ride->dataPoints()[i];
// compute bearing in order to calculate headwind
if (i>=1)
{
RideFilePoint *prevPoint = ride->dataPoints()[i-1];
// ensure a movement occurred and valid lat/lon in order to compute cyclist direction
if ( (prevPoint->lat != p->lat || prevPoint->lon != p->lon )
&& (prevPoint->lat != 0 || prevPoint->lon != 0 )
&& (p->lat != 0 || p->lon != 0 ) )
bearing = atan2(cos(p->lat)*sin(p->lon - prevPoint->lon),
cos(prevPoint->lat)*sin(p->lat)-sin(prevPoint->lat)*cos(p->lat)*cos(p->lon - prevPoint->lon));
}
// else keep previous bearing (or 0 at the beginning)
// wind parameters to be considered
if (windSpeed || windHeading) // if wind parameters were entered manually then use it and override rideFile headwind
W = cos(bearing - windHeading) * windSpeed * 0.27777777777778; // Absolute wind speed relative to cyclist orientation (m/s)
else if (ride->areDataPresent()->headwind) // otherwise use headwind from rideFile records (typ. weather forecast included in FIT files)
W = (p->headwind - p->kph) * 0.27777777777778; // Compute absolute wind speed relative to cyclist orientation (m/s)
else
W = 0.0; //otherwise assume no wind
// Estimate Power if not in data
double cad = ride->areDataPresent()->cad ? p->cad : 85.00;
if (cad > 0) {
if (ride->areDataPresent()->temp) T = p->temp;
double Slope = atan(p->slope * .01);
double V = p->kph * 0.27777777777778; // Cyclist speed m/s
double CrDyn = 0.1 * cos(Slope);
double CwaRider, Ka;
double Frg = 9.81 * (MBik + M) * (CrEff * cos(Slope) + sin(Slope));
double vw=V+W; // Wind speed against cyclist = cyclist speed + wind speed
CwaRider = (1 + cad * cCad) * afCd * adipos * (((hRider - adipos) * afSin) + adipos);
Ka = 176.5 * exp(-p->alt * .0001253) * (CwaRider + CwaBike) / (273 + T);
//qDebug()<<"acc="<<p->kphd<<" , V="<<V<<" , m="<<M<<" , Pa="<<(p->kphd > 1 ? 1 : p->kphd*V*M);
double watts = (afCm * V * (Ka * (vw * vw) + Frg + V * CrDyn))+(p->kphd > 1 ? 1 : p->kphd*V*M);
ride->command->setPointValue(i, RideFile::watts, watts > 0 ? (watts > 1000 ? 1000 : watts) : 0);
// qDebug() << "watts = "<<p->watts;
// qDebug() << " " << afCm * V * Ka * (vw * vw) << " = afCm(=" << afCm << ") * V(=" << V << ") * Ka(="<<Ka<<") * (vw^2(=" << V+W << "^2))";
// qDebug() << " " << afCm * V * Frg << " = afCm * V * Frg(=" << Frg << ")";
// qDebug() << " " << afCm * V * V * CrDyn << " = afCm * V^2 * CrDyn(=" << CrDyn << ")";
// qDebug() << " " << p->kphd*V*M << " = kphd(=" << p->kphd << ") * V * M(=" << M << ")";
// qDebug() << " Ka="<<Ka<<", CwaRi="<<CwaRider<<", slope="<<p->slope<<", v="<<p->kph<<" Cwa="<<(CwaRider + CwaBike);
} else {
ride->command->setPointValue(i, RideFile::watts, 0);
}
}
int smoothPoints = 3;
// initialise rolling average
double rtot = 0;
for (int i=smoothPoints; i>0 && ride->dataPoints().count()-i >=0; i--) {
rtot += ride->dataPoints()[ride->dataPoints().count()-i]->watts;
}
// now run backwards setting the rolling average
for (int i=ride->dataPoints().count()-1; i>=smoothPoints; i--) {
double here = ride->dataPoints()[i]->watts;
ride->dataPoints()[i]->watts = rtot / smoothPoints;
if (ride->dataPoints()[i]->watts<0) ride->dataPoints()[i]->watts = 0;
rtot -= here;
rtot += ride->dataPoints()[i-smoothPoints]->watts;
}
ride->setDataPresent(ride->watts, true);
}
ride->command->endLUW();
return true;
}
Reference in New Issue View Git Blame Copy Permalink