1
0
mirror of https://github.com/tumic0/GPXSee.git synced 2024-10-06 23:03:22 +02:00
GPXSee/src/map/angularunits.cpp
2018-10-11 18:19:35 +02:00

94 lines
1.8 KiB
C++

#include "common/coordinates.h"
#include "common/util.h"
#include "angularunits.h"
static double sDMS2deg(double val)
{
double angle;
const char *decimal;
QString qstr(QString::number(qAbs(val), 'f', 7));
QByteArray ba = qstr.toLatin1();
const char *str = ba.constData();
decimal = strrchr(str, '.');
int deg = str2int(str, decimal - str);
int min = str2int(decimal + 1, 2);
int sec = str2int(decimal + 3, 2);
int f = str2int(decimal + 5, 3);
angle = deg + min/60.0 + sec/3600.0 + (f/1000.0)/3600.0;
return (val < 0) ? -angle : angle;
}
static double deg2sDMS(double val)
{
double aval = qAbs(val);
int deg = aval;
double r1 = aval - deg;
int min = r1 * 60.0;
double r2 = r1 - (min / 60.0);
int sec = r2 * 3600.0;
double r3 = r2 - (sec / 3600.0);
int f = (int)(r3 * 3600.0 * 1000.0);
QString str(QString("%1.%2%3%4").arg(deg).arg(min, 2, 10, QChar('0'))
.arg(sec, 2, 10, QChar('0')).arg(f, 3, 10, QChar('0')));
return (val < 0) ? -str.toDouble() : str.toDouble();
}
AngularUnits::AngularUnits(int code) : _code(code)
{
switch (code) {
case 9101:
_f = 180.0 / M_PI;
break;
case 9102:
case 9107:
case 9108:
case 9110:
case 9122:
_f = 1.0;
break;
case 9103:
_f = 1 / 60.0;
break;
case 9104:
_f = 1 / 3600.0;
break;
case 9105:
_f = 180.0 / 200.0;
break;
case 9106:
_f = 180.0 / 200.0;
break;
case 9109:
_f = 180.0 / (M_PI * 1000000.0);
break;
default:
_f = NAN;
}
}
double AngularUnits::toDegrees(double val) const
{
return (_code == 9110) ? sDMS2deg(val) : val * _f;
}
double AngularUnits::fromDegrees(double val) const
{
return (_code == 9110) ? deg2sDMS(val) : val / _f;
}
#ifndef QT_NO_DEBUG
QDebug operator<<(QDebug dbg, const AngularUnits &au)
{
dbg.nospace() << "AngularUnits(" << deg2rad(au._f) << ")";
return dbg.space();
}
#endif // QT_NO_DEBUG