CoordsAlgorithms.h

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
/***********************************************************************************/
/*  Copyright 2009 WSL Institute for Snow and Avalanche Research    SLF-DAVOS      */
/***********************************************************************************/
/* This file is part of MeteoIO.
    MeteoIO is free software: you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    MeteoIO is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU Lesser General Public License for more details.
 
    You should have received a copy of the GNU Lesser General Public License
    along with MeteoIO.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef COORDSALGORITHMS_H
#define COORDSALGORITHMS_H
 
#include <string>
 
namespace mio {
class CoordsAlgorithms {
public:
        enum ELLIPSOIDS_NAMES {
                E_WGS84, 
                E_WGS72, 
                E_GRS80, 
                E_AIRY, 
                E_INTL1924, 
                E_CLARKE1880, 
                E_CLARKE1866, 
                E_GRS67, 
                E_EVEREST1830, 
                E_IERS2003, 
                E_KRASSOVSKY, 
                NONE 
        };
 
        //Handling of lat/lon
        static double dms_to_decimal(const std::string& dms);
        static std::string decimal_to_dms(const double& decimal);
        static void parseLatLon(const std::string& coordinates, double& lat, double& lon);
        static std::string printLatLon(const double& latitude, const double& longitude);
        static double lat_degree_lenght(const double& latitude);
        static double lon_degree_lenght(const double& latitude);
        
        //EPSG helper methods
        static bool isUTM(const int& epsg);
        static short int str_to_EPSG(const std::string& coordsystem, const std::string& coordparam);
        static void EPSG_to_str(const int& epsg, std::string& coordsystem, std::string& coordparam);
 
        //handling of rotated coordinates and datum conversions
        static void rotatedToTrueLatLon(const double& lat_N, const double& lon_N, const double& lat_rot, const double& lon_rot, double &lat_true, double &lon_true);
        static void trueLatLonToRotated(const double& lat_N, const double& lon_N, const double& lat_true, const double& lon_true, double &lat_rot, double &lon_rot);
        static void Molodensky(const double& lat_in, const double& lon_in, const double& alt_in, const ELLIPSOIDS_NAMES& ellipsoid_in,
                                               double &lat_out, double &lon_out, double &alt_out, const ELLIPSOIDS_NAMES& ellipsoid_out,
                                               const double& delta_x=0., const double& delta_y=0., const double& delta_z=0.);
 
        //handling of distances on a sphere
        static double cosineDistance(const double& lat1, const double& lon1, const double& lat2, const double& lon2, double& alpha);
        static void cosineInverse(const double& lat_ref, const double& lon_ref, const double& distance, const double& bearing, double& lat, double& lon);
        static double VincentyDistance(const double& lat1, const double& lon1, const double& lat2, const double& lon2, double& alpha);
        static void VincentyInverse(const double& lat_ref, const double& lon_ref, const double& distance, const double& bearing, double& lat, double& lon);
 
        //Coordinates conversions
        static void WGS84_to_CH1903(const double& lat_in, const double& long_in, double& east_out, double& north_out);
        static void CH1903_to_WGS84(const double& east_in, const double& north_in, double& lat_out, double& long_out);
        static void WGS84_to_UTM(const double& lat_in, double long_in, const std::string& coordparam, double& east_out, double& north_out);
        static void UTM_to_WGS84(double east_in, double north_in, const std::string& coordparam, double& lat_out, double& long_out);
        static void WGS84_to_UPS(const double& lat_in, const double& long_in, const std::string& coordparam, double& east_out, double& north_out);
        static void UPS_to_WGS84(const double& east_in, const double& north_in, const std::string& coordparam, double& lat_out, double& long_out);
        static void WGS84_to_PROJ(const double& lat_in, const double& long_in, const std::string& coordparam, double& east_out, double& north_out);
        static void PROJ_to_WGS84(const double& east_in, const double& north_in, const std::string& coordparam, double& lat_out, double& long_out);
 
        static int getUTMZone(const double& latitude, const double& longitude, std::string& zone_out);
        static void parseUTMZone(const std::string& zone_info, char& zoneLetter, short int& zoneNumber);
 
 private:
        struct ELLIPSOID {
                double a;
                double b;
        };
        static const struct ELLIPSOID ellipsoids[12];
};
} //end namespace
 
#endif