flightsimulatorplatfrom/pilotclient
1
0
Fork 0
mirror of https://github.com/swift-project/pilotclient.git synced 2026-03-22 14:55:36 +08:00
Code Issues Packages Projects Releases Wiki Activity

add greatCircleDistance

Browse Source 
refs #81
This commit is contained in:
Klaus Basan
2013-12-10 19:32:02 +00:00
committed by Mathew Sutcliffe
parent ded6fc012e
commit 0518e680c3
2 changed files with 35 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
src/blackmisc/coordinategeodetic.cpp
View File

@@ -68,5 +68,31 @@ namespace BlackMisc
}
/*
* Great circle distance
*/
PhysicalQuantities::CLength greatCircleDistance(const ICoordinateGeodetic &coordinate1, const ICoordinateGeodetic &coordinate2)
{
// same coordinate results in 0 distance
if (coordinate1.latitude() == coordinate2.latitude() && coordinate1.longitude() == coordinate2.longitude())
{
return CLength(0, CLengthUnit::NM());
}
// first, prelimary distance calculation
// http://www.geodatasource.com/developers/c
double dist;
double lon1rad = coordinate1.longitude().value(CAngleUnit::rad());
double lon2rad = coordinate2.longitude().value(CAngleUnit::rad());
double lat1rad = coordinate1.latitude().value(CAngleUnit::rad());
double lat2rad = coordinate2.latitude().value(CAngleUnit::rad());
double theta = lon1rad - lon2rad;
dist = qSin(lat1rad) * qSin(lat2rad) + qCos(lat1rad) * qCos(lat2rad) * cos(theta);
dist = qAcos(dist);
dist = CMath::rad2deg(dist);
dist = dist * 60; // dist in NM
return CLength(qAbs(dist), CLengthUnit::NM());
}
} // namespace
} // namespace

10
src/blackmisc/coordinategeodetic.h
View File

@@ -55,10 +55,18 @@ namespace BlackMisc
};
/*!
* \brief Great circle distance between points
* \param coordinate1
* \param coordinate2
* \return
*/
BlackMisc::PhysicalQuantities::CLength greatCircleDistance(const ICoordinateGeodetic &coordinate1, const ICoordinateGeodetic &coordinate2);
/*!
* \brief Geodetic coordinate
*/
class CCoordinateGeodetic : public CValueObject
class CCoordinateGeodetic : public CValueObject, public ICoordinateGeodetic
{
private:
BlackMisc::Geo::CLatitude m_latitude; //!< Latitude