GeographicProjection.h

#pragma once
 
#include "Ellipsoid.h"
#include "GlobeRectangle.h"
#include "Library.h"
 
#include <CesiumUtility/Math.h>
 
#include <glm/vec2.hpp>
#include <glm/vec3.hpp>
 
namespace CesiumGeospatial {
 
class Cartographic;
 
class CESIUMGEOSPATIAL_API GeographicProjection final {
public:
  static constexpr GlobeRectangle MAXIMUM_GLOBE_RECTANGLE = GlobeRectangle(
      -CesiumUtility::Math::OnePi,
      -CesiumUtility::Math::PiOverTwo,
      CesiumUtility::Math::OnePi,
      CesiumUtility::Math::PiOverTwo);
 
  static constexpr CesiumGeometry::Rectangle computeMaximumProjectedRectangle(
      const Ellipsoid& ellipsoid CESIUM_DEFAULT_ELLIPSOID) noexcept {
    const double longitudeValue =
        ellipsoid.getMaximumRadius() * CesiumUtility::Math::OnePi;
    const double latitudeValue =
        ellipsoid.getMaximumRadius() * CesiumUtility::Math::PiOverTwo;
    return CesiumGeometry::Rectangle(
        -longitudeValue,
        -latitudeValue,
        longitudeValue,
        latitudeValue);
  }
 
  GeographicProjection(
      const Ellipsoid& ellipsoid CESIUM_DEFAULT_ELLIPSOID) noexcept;
 
  const Ellipsoid& getEllipsoid() const noexcept { return this->_ellipsoid; }
 
  glm::dvec3 project(const Cartographic& cartographic) const noexcept;
 
  CesiumGeometry::Rectangle
  project(const CesiumGeospatial::GlobeRectangle& rectangle) const noexcept;
 
  Cartographic unproject(const glm::dvec2& projectedCoordinates) const noexcept;
 
  Cartographic unproject(const glm::dvec3& projectedCoordinates) const noexcept;
 
  CesiumGeospatial::GlobeRectangle
  unproject(const CesiumGeometry::Rectangle& rectangle) const noexcept;
 
  bool operator==(const GeographicProjection& rhs) const noexcept {
    return this->_ellipsoid == rhs._ellipsoid;
  };
 
  bool operator!=(const GeographicProjection& rhs) const noexcept {
    return !(*this == rhs);
  };
 
private:
  Ellipsoid _ellipsoid;
  double _semimajorAxis;
  double _oneOverSemimajorAxis;
};
 
} // namespace CesiumGeospatial