Method cartesianToCartographic

CesiumGeospatial/src/Ellipsoid.cpp:47–63

}

std::optional<Cartographic>
Ellipsoid::cartesianToCartographic(const glm::dvec3& cartesian) const noexcept {
  std::optional<glm::dvec3> p = this->scaleToGeodeticSurface(cartesian);

  if (!p) {
    return std::optional<Cartographic>();
  }

  const glm::dvec3 n = this->geodeticSurfaceNormal(p.value());
  const glm::dvec3 h = cartesian - p.value();

  const double longitude = glm::atan(n.y, n.x);
  const double latitude = glm::asin(n.z);
  const double height = Math::sign(glm::dot(h, cartesian)) * glm::length(h);

  return Cartographic(longitude, latitude, height);
}

std::optional<glm::dvec3>
Ellipsoid::scaleToGeodeticSurface(const glm::dvec3& cartesian) const noexcept {