geodetictools.h

Go to the documentation of this file.

/* ============================================================
 *
 * This file is a part of digiKam project
 * https://www.digikam.org
 *
 * Date        : 2008-05-05
 * Description : Geodetic tools
 *
 * Copyright (C) 2008-2011 by Marcel Wiesweg <marcel dot wiesweg at gmx dot de>
 * Copyright (C) 2015-2022 by Gilles Caulier <caulier dot gilles at gmail dot 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, 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.
 *
 * ============================================================ */
 
#ifndef DIGIKAM_GEODETIC_TOOLS_H
#define DIGIKAM_GEODETIC_TOOLS_H
 
// C++ includes
 
#include <cmath>
 
// Qt includes
 
#include <QString>
#include <QPointF>
 
// Local includes
 
#include "digikam_export.h"
 
namespace Digikam
{
 
namespace Coordinates
{
 
 
inline double toRadians(double deg)
{
    return (deg * M_PI / 180.0);
}
 
inline double toRadiansFactor()
{
    return (M_PI / 180.0);
}
 
inline double toDegrees(double rad)
{
    return (rad * 180.0 / M_PI);
}
 
inline double toDegreesFactor()
{
    return (180.0 / M_PI);
}
 
} // namespace Coordinates
 
// ------------------------------------------------------------------------------------------------
 
class DIGIKAM_EXPORT Ellipsoid
{
 
public:
 
    static Ellipsoid WGS84();
 
    static Ellipsoid GRS80();
 
    static Ellipsoid INTERNATIONAL_1924();
 
    static Ellipsoid CLARKE_1866();
 
    static Ellipsoid SPHERE();
 
    static Ellipsoid createEllipsoid(const QString& name,
                                     double semiMajorAxis,
                                     double semiMinorAxis);
 
    static Ellipsoid createFlattenedSphere(const QString& name,
                                           double semiMajorAxis,
                                           double inverseFlattening);
 
    double semiMajorAxis() const;
 
    double semiMinorAxis() const;
 
    double eccentricity() const;
 
    double inverseFlattening() const;
 
    bool isIvfDefinitive() const;
 
    bool isSphere() const;
 
    double orthodromicDistance(double x1, double y1, double x2, double y2);
 
    double radiusOfCurvature(double latitude);
 
protected:
 
    Ellipsoid(const QString& name,
              double semiMajorAxis,
              double semiMinorAxis,
              double inverseFlattening,
              bool ivfDefinitive);
    Ellipsoid(const QString& name,
              double radius,
              bool ivfDefinitive);
 
protected:
 
    QString name;
 
    double m_semiMajorAxis;
 
    double m_semiMinorAxis;
 
    double m_inverseFlattening;
 
    bool   m_ivfDefinitive;
 
    bool   m_isSphere;
};
 
// ------------------------------------------------------------------------------------------------
 
class DIGIKAM_EXPORT GeodeticCalculator
{
public:
 
    explicit GeodeticCalculator(const Ellipsoid& e = Ellipsoid::WGS84());
 
    Ellipsoid ellipsoid() const;
 
    void setStartingGeographicPoint(double longitude, double latitude);
 
    void setDestinationGeographicPoint(double longitude, double latitude);
 
    bool destinationGeographicPoint(double* longitude, double* latitude);
    QPointF destinationGeographicPoint();
 
    void setDirection(double azimuth, double distance);
 
    double azimuth();
 
    double orthodromicDistance();
 
    bool checkOrthodromicDistance();
 
    bool computeDestinationPoint();
 
    double meridianArcLength(double latitude1, double latitude2);
 
    double meridianArcLengthRadians(double P1, double P2);
 
    bool computeDirection();
 
protected:
 
    double castToAngleRange(const double alpha);
 
    bool checkLatitude(double* latitude);
 
    bool checkLongitude(double* longitude);
 
    bool checkAzimuth(double* azimuth);
 
    bool checkOrthodromicDistance(const double distance);
 
    double TOLERANCE_0, TOLERANCE_1, TOLERANCE_2, TOLERANCE_3;
 
    double TOLERANCE_CHECK;
 
    Ellipsoid m_ellipsoid;
 
    double m_semiMajorAxis;
 
    double m_semiMinorAxis;
 
    double m_eccentricitySquared;
 
    double m_maxOrthodromicDistance;
 
    double A, B, C, D, E, F;
 
    double fo, f, f2, f3, f4;
 
    double T1, T2, T4, T6;
 
    double a01, a02, a03, a21, a22, a23, a42, a43, a63;
 
    double m_lat1, m_long1;
 
    double m_lat2, m_long2;
 
    double m_distance, m_azimuth;
 
    bool m_destinationValid;
 
    bool m_directionValid;
};
 
} // namespace Digikam
 
#endif // DIGIKAM_GEODETIC_TOOLS_H