mirror of
https://github.com/swift-project/pilotclient.git
synced 2026-03-22 23:05:36 +08:00
85 lines
3.6 KiB
C++
85 lines
3.6 KiB
C++
/* Copyright (C) 2015
|
||
* swift project Community / Contributors
|
||
*
|
||
* This file is part of swift project. It is subject to the license terms in the LICENSE file found in the top-level
|
||
* directory of this distribution. No part of swift project, including this file, may be copied, modified, propagated,
|
||
* or distributed except according to the terms contained in the LICENSE file.
|
||
*/
|
||
|
||
//! \file
|
||
//! \ingroup sampleblackmiscquantities
|
||
|
||
#include "samplesgeo.h"
|
||
#include "blackmisc/aviation/altitude.h"
|
||
#include "blackmisc/geo/coordinategeodetic.h"
|
||
#include "blackmisc/geo/earthangle.h"
|
||
#include "blackmisc/geo/latitude.h"
|
||
#include "blackmisc/geo/longitude.h"
|
||
#include "blackmisc/pq/physicalquantity.h"
|
||
#include "blackmisc/pq/units.h"
|
||
#include "blackmisc/stringutils.h"
|
||
|
||
#include <QString>
|
||
#include <QTextStream>
|
||
#include <QVector3D>
|
||
|
||
using namespace BlackMisc;
|
||
using namespace BlackMisc::Aviation;
|
||
using namespace BlackMisc::PhysicalQuantities;
|
||
using namespace BlackMisc::Geo;
|
||
|
||
namespace BlackSample
|
||
{
|
||
int CSamplesGeo::samples(QTextStream &out)
|
||
{
|
||
// comparions https://www.fcc.gov/general/degrees-minutes-seconds-tofrom-decimal-degrees
|
||
const int digits = 12;
|
||
const QString latStr("N 48° 7′ 6.3588"); // 48.118433
|
||
const QString lngStr("E 16° 33′ 39.924");
|
||
out.setRealNumberPrecision(digits);
|
||
CCoordinateGeodetic geo = CCoordinateGeodetic::fromWgs84(latStr, lngStr, CAltitude(100, CAltitude::MeanSeaLevel, CLengthUnit::m()));
|
||
CLatitude lat = CLatitude::fromWgs84(latStr);
|
||
CLongitude lng = CLongitude::fromWgs84(lngStr);
|
||
|
||
// check if conversions to xyz have messed something up
|
||
CLatitude deltaLat = geo.latitude() - lat;
|
||
CLongitude deltaLng = geo.longitude() - lng;
|
||
|
||
out << latStr << " " << lngStr << endl;
|
||
out <<
|
||
lat.value(CAngleUnit::deg()) << " " << lat.value(CAngleUnit::sexagesimalDeg()) << " " <<
|
||
lng.value(CAngleUnit::deg()) << " " << lng.value(CAngleUnit::sexagesimalDeg()) << endl;
|
||
out <<
|
||
geo.latitude().value(CAngleUnit::deg()) << " " << geo.latitude().value(CAngleUnit::sexagesimalDeg()) << " " <<
|
||
geo.longitude().value(CAngleUnit::deg()) << " " << geo.longitude().value(CAngleUnit::sexagesimalDeg()) << endl;
|
||
|
||
out << deltaLat.valueRoundedWithUnit(digits) << " " << deltaLng.valueRoundedWithUnit(digits) << endl;
|
||
|
||
// equal test
|
||
out << "Equal? " <<
|
||
BlackMisc::boolToYesNo(lat == geo.latitude()) << " " <<
|
||
BlackMisc::boolToYesNo(lng == geo.longitude()) << endl;
|
||
|
||
// check if conversions to xyz have messed something up
|
||
QVector3D geoVector = geo.normalVector();
|
||
CCoordinateGeodetic geo2(geoVector);
|
||
deltaLat = geo2.latitude() - lat;
|
||
deltaLng = geo2.longitude() - lng;
|
||
out << deltaLat.valueRoundedWithUnit(digits) << " " << deltaLng.valueRoundedWithUnit(digits) << endl;
|
||
|
||
// Heading/bearing of same values
|
||
CAngle bearing = geo.calculateBearing(geo);
|
||
CLength distance = geo.calculateGreatCircleDistance(geo);
|
||
out << bearing.valueRoundedWithUnit(CAngleUnit::deg(), 2) << " " << distance.valueRoundedWithUnit(CLengthUnit::m(), 2) << endl;
|
||
|
||
const CCoordinateGeodetic nullCoordinate;
|
||
bearing = geo.calculateBearing(nullCoordinate);
|
||
distance = geo.calculateGreatCircleDistance(nullCoordinate);
|
||
out << bearing.valueRoundedWithUnit(CAngleUnit::deg(), 2) << " " << distance.valueRoundedWithUnit(CLengthUnit::m(), 2) << endl;
|
||
|
||
// bye
|
||
out << "-----------------------------------------------" << endl;
|
||
return 0;
|
||
}
|
||
} // namespace
|