geodetic.hpp

Go to the documentation of this file.

#pragma once
 
#include "../centers.hpp"
#include "../ffi_core.hpp"
#include "../frames.hpp"
 
#include <qtty/qtty.hpp>
 
#include <ostream>
 
namespace siderust {
namespace cartesian {
template <typename C, typename F, typename U> struct Position;
}
 
struct Geodetic {
  qtty::Degree lon;   
  qtty::Degree lat;   
  qtty::Meter height; 
 
  Geodetic() : lon(qtty::Degree(0)), lat(qtty::Degree(0)), height(qtty::Meter(0)) {}
 
  Geodetic(qtty::Degree lon_, qtty::Degree lat_, qtty::Meter h = qtty::Meter(0))
      : lon(lon_), lat(lat_), height(h) {}
 
  Geodetic(double lon_deg, double lat_deg, double height_m = 0.0)
      : lon(qtty::Degree(lon_deg)), lat(qtty::Degree(lat_deg)), height(qtty::Meter(height_m)) {}
 
  siderust_geodetic_t to_c() const { return {lon.value(), lat.value(), height.value()}; }
 
  static Geodetic from_c(const siderust_geodetic_t &c) {
    return Geodetic(c.lon_deg, c.lat_deg, c.height_m);
  }
 
  template <typename U = qtty::Meter>
  cartesian::Position<centers::Geocentric, frames::ECEF, U> to_cartesian() const;
};
 
// ============================================================================
// Stream operators
// ============================================================================
 
inline std::ostream &operator<<(std::ostream &os, const Geodetic &geo) {
  return os << "Geodetic (lon=" << geo.lon << ", lat=" << geo.lat << ", h=" << geo.height << ')';
}
 
} // namespace siderust