EllipsoidTangentPlane.h

#pragma once
 
#include <CesiumGeometry/Plane.h>
#include <CesiumGeospatial/Ellipsoid.h>
#include <CesiumGeospatial/Library.h>
 
#include <glm/fwd.hpp>
#include <glm/vec2.hpp>
#include <glm/vec3.hpp>
 
namespace CesiumGeospatial {

class CESIUMGEOSPATIAL_API EllipsoidTangentPlane final {
public:
  EllipsoidTangentPlane(
      const glm::dvec3& origin,
      const Ellipsoid& ellipsoid CESIUM_DEFAULT_ELLIPSOID);

  EllipsoidTangentPlane(
      const glm::dmat4& eastNorthUpToFixedFrame,
      const Ellipsoid& ellipsoid CESIUM_DEFAULT_ELLIPSOID);

  const Ellipsoid& getEllipsoid() const noexcept { return this->_ellipsoid; }

  const glm::dvec3& getOrigin() const noexcept { return this->_origin; }

  const glm::dvec3& getXAxis() const noexcept { return this->_xAxis; }

  const glm::dvec3& getYAxis() const noexcept { return this->_yAxis; }

  const glm::dvec3& getZAxis() const noexcept {
    return this->_plane.getNormal();
  }

  const CesiumGeometry::Plane& getPlane() const noexcept {
    return this->_plane;
  }

  glm::dvec2
  projectPointToNearestOnPlane(const glm::dvec3& cartesian) const noexcept;
 
private:
  static glm::dmat4 computeEastNorthUpToFixedFrame(
      const glm::dvec3& origin,
      const Ellipsoid& ellipsoid CESIUM_DEFAULT_ELLIPSOID);
 
  Ellipsoid _ellipsoid;
  glm::dvec3 _origin;
  glm::dvec3 _xAxis;
  glm::dvec3 _yAxis;
  CesiumGeometry::Plane _plane;
};
 
} // namespace CesiumGeospatial