ProjectionAEA.h

Go to the documentation of this file.

// This file is part of the ExaHyPE2 project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
#include <cstddef>
 
// TODO: Using std=c++23 constexpr would be possible
 
namespace projection::AEA {
  namespace ellipsoid {
    /*
     * @brief
     * Projects geographic coordinates (longitude, latitude) onto a planar surface
     * using the Albers Equal-Area Conic projection for the WGS_1984 ellipsoid.
     *
     * This implementation is based on the equations provided in USGS Professional Paper 1395.
     * It supports standard parallels and is suitable for regional-scale mapping where
     * area preservation is important.
     *
     * @see https://doi.org/10.3133/pp1395
     */
    template <auto EllipsoidA, auto EllipsoidE, auto StdParallel1Deg, auto StdParallel2Deg, auto Lon0Deg, auto Lat0Deg>
    static inline auto forward(const auto lonDeg, const auto latDeg);
 
    template <
      auto EllipsoidA,
      auto EllipsoidE,
      auto StdParallel1Deg,
      auto StdParallel2Deg,
      auto Lon0Deg,
      auto Lat0Deg,
      auto MaxIter   = 10,
      auto Tolerance = 1e-12>
    static inline auto inverse(const auto x, const auto y);
  } // namespace ellipsoid
 
  namespace sphere {
    template <auto Lon0Deg, auto Lat0Deg, auto StdParallel1Deg, auto StdParallel2Deg, auto EarthRadius>
    static inline auto forward(const auto lon, const auto lat);
 
    template <auto Lon0Deg, auto Lat0Deg, auto StdParallel1Deg, auto StdParallel2Deg, auto EarthRadius>
    static inline auto inverse(const auto x, const auto y);
  } // namespace sphere
} // namespace projection::AEA
 
#include "ProjectionAEA.cpph"