WebMercatorProjection.h

#pragma once
 
#include "Ellipsoid.h"
#include "GlobeRectangle.h"
#include "Library.h"
 
#include <glm/vec2.hpp>
#include <glm/vec3.hpp>
 
namespace CesiumGeospatial {
 
class Cartographic;
 
class CESIUMGEOSPATIAL_API WebMercatorProjection final {
public:
  static const double MAXIMUM_LATITUDE;
 
  static const GlobeRectangle MAXIMUM_GLOBE_RECTANGLE;
 
  static constexpr CesiumGeometry::Rectangle computeMaximumProjectedRectangle(
      const Ellipsoid& ellipsoid CESIUM_DEFAULT_ELLIPSOID) noexcept {
    const double value =
        ellipsoid.getMaximumRadius() * CesiumUtility::Math::OnePi;
    return CesiumGeometry::Rectangle(-value, -value, value, value);
  }
 
  WebMercatorProjection(
      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;
 
  static double mercatorAngleToGeodeticLatitude(double mercatorAngle) noexcept;
 
  static double geodeticLatitudeToMercatorAngle(double latitude) noexcept;
 
  bool operator==(const WebMercatorProjection& rhs) const noexcept {
    return this->_ellipsoid == rhs._ellipsoid;
  };
 
  bool operator!=(const WebMercatorProjection& rhs) const noexcept {
    return !(*this == rhs);
  };
 
private:
  Ellipsoid _ellipsoid;
  double _semimajorAxis;
  double _oneOverSemimajorAxis;
};
 
} // namespace CesiumGeospatial