GeodeticUtils.h

#pragma once
 
#include "../FoundationBase.h"
 
#include <glm/vec2.hpp>
#include <glm/vec3.hpp>
#include <glm/mat3x3.hpp>
 
COGSFOUNDATION_API void LLtoUTMZone(int&         zoneNumber,
                                    bool&        northernHemisphere,
                                    const double latitude,
                                    const double longitude,
                                    const int    referenceEllipsoid = 23);
 
COGSFOUNDATION_API void LLtoUTM(double&      easting,
                                double&      northing,
                                const double latitude,
                                const double longitude,
                                int          zoneNumber,
                                bool         northernHemisphere,
                                const int    referenceEllipsoid = 23);
 
COGSFOUNDATION_API glm::dvec3 LLtoECEF(double latitude, double longitude, double altitude = 0.0);
 
COGSFOUNDATION_API glm::dvec2 calcUTMVectorBetweenCoordinates(double lat1, double long1, double lat2, double long2, int zone1 = -1, int zone2 = -1);
 
COGSFOUNDATION_API float calcBearingBetweenCoordinates(double lat1, double long1, double lat2, double long2, int zone = -1);
 
namespace Cogs {
 
  class COGSFOUNDATION_API LLToENUConverter {
  public:
                          LLToENUConverter(double latitude, double longitude, double altitude = 0.0);
                          LLToENUConverter(const LLToENUConverter&) = default;
                          LLToENUConverter(LLToENUConverter&&) = default;
 
            glm::dvec3    convert(double latitude, double longitude, double altitude = 0.0) const;
            glm::dvec2    calcVectorBetweenCoordinates(double lat1, double long1, double lat2, double long2) const;
 
  private:
            glm::dvec3    ecefRef;
            glm::dmat3x3  ecefToENUMatrix;
  };
}