convert.h

Go to the documentation of this file.

/*
  Copyright (c) Mircea Neacsu (2014-2025) Licensed under MIT License.
  This file is part of MLIB project. See LICENSE file for full license terms.
*/
 
#pragma once
 
#include <math.h>
#include <numbers>
#include <string>
#include <cmath>
 
#ifndef M_PI
#define M_PI std::numbers::pi
#endif
 
// Useful constants

constexpr double D2R = M_PI / 180.;

constexpr double A_WGS84 = 6378137.000;

constexpr double F_WGS84 = 0.003352810664747;

constexpr double F1_WGS84 = 298.257223563;
 
#define NM2M     1852.                     
#define MPS2KNOT (3600. / NM2M)            
#define USFOOT2M (1200. / 3937.)           

constexpr double MAS = M_PI / (180 * 3600000.);
 

constexpr double operator"" _ftUS (long double ftus)
{
  return ftus * USFOOT2M;
}
 
constexpr double operator"" _ftUS (unsigned long long ftus)
{
  return ftus * USFOOT2M;
}


constexpr double operator"" _nmi (long double nmi)
{
  return nmi * NM2M;
}
 
constexpr double operator"" _nmi (unsigned long long nmi)
{
  return nmi * NM2M;
}

 
// ================== Degrees conversion functions ==========================
 
// Decimal degrees---------------------------------------------------------

constexpr double DEG (double dd)
{
  return dd * D2R;
}

constexpr double D2rad (double dd)
{
  return dd * D2R;
}

constexpr double rad2D (double r)
{
  return r / D2R;
}
 
// Degrees, minutes --------------------------------------------------------

constexpr double DM (double dd, double mm)
{
  return (dd + mm / 60.) * D2R;
}

constexpr double DM2deg (double ddmm)
{
  int sign = (ddmm >= 0) ? 1 : -1;
  if (ddmm < 0)
    ddmm = -ddmm;
  int deg = (int)(ddmm / 100.);
  ddmm -= (double)deg * 100.;
  return sign * (deg + ddmm / 60.);
}

constexpr double deg2DM (double dd)
{
  int deg = (int)dd;
  return (dd - deg) * 60. + deg * 100.;
}

constexpr double rad2DM (double rad)
{
  return deg2DM (rad / D2R);
}

constexpr double DM2rad (double val)
{
  return DM2deg (val) * D2R;
}
 
// Degrees, minutes, seconds ------------------------------------------------

constexpr double DMS (double dd, double mm, double ss)
{
  return (dd + mm / 60. + ss / 3600.) * D2R;
}

constexpr double DMS2deg (double dms)
{
  int sign = (dms >= 0) ? 1 : -1;
  if (dms < 0)
    dms = -dms;
  int deg = (int)(dms / 10000.);
  dms -= (double)deg * 10000;
  int min = (int)(dms / 100.);
  dms -= (double)min * 100.;
  return sign * (deg + min / 60. + dms / 3600.);
}

constexpr double DMS2rad (double dms)
{
  return DMS2deg (dms) * D2R;
}
 
// =================== Angle units literal operators =========================

constexpr double operator"" _deg (long double deg)
{
  return deg * D2R;
}
 
constexpr double operator"" _deg (unsigned long long deg)
{
  return deg * D2R;
}


constexpr double operator"" _arcmin (long double min)
{
  return min/60. * D2R;
}
 
constexpr double operator"" _arcmin (unsigned long long min)
{
  return min/60. * D2R;
}


constexpr double operator ""_dm (long double val)
{
  return DM2rad (val);
}
 
constexpr double operator ""_dm (unsigned long long val)
{
  return DM2rad ((double)val);
}


constexpr double operator ""_dms (long double val)
{
  return DMS2rad (val);
}
 
constexpr double operator ""_dms (unsigned long long val)
{
  return DMS2rad ((double)val);
}


constexpr double operator"" _arcsec (long double sec)
{
  return sec/3600. * D2R;
}
 
constexpr double operator"" _arcsec (unsigned long long sec)
{
  return sec/3600. * D2R;
}


constexpr double operator"" _mas (long double mas)
{
  return mas * MAS;
}
 
constexpr double operator"" _mas (unsigned long long mas)
{
  return mas * MAS;
}

 
 
namespace mlib {
double deg_reduce (double value);

enum deg_fmt
{
  degrees = 0, 
  minutes = 1, 
  seconds = 2 
};

std::string degtoa (double degrees, deg_fmt format, bool latitude, int precision);

double atodeg (const std::string& str);

template <typename T>
std::pair<T, T> sincos (T val)
{
  return {std::sin (val), std::cos (val)};
}
} // namespace mlib