Datum.h

Go to the documentation of this file.

#ifndef DATUM_H_
#define DATUM_H_

#include <math.h>
#include <cassert>

#include "AuvMath.h"

#include "../../Lib/geotrans/dt_cc/utm/utm.h"
#include "../../Lib/geotrans/dt_cc/tranmerc/tranmerc.h"

static const double UTM_ZONE_WIDTH = M_PI / 30.;

class Wgs84
{
public:

    static long UtmToLatLon( const double northings, const double eastings,
                             const unsigned int zone, const bool northernHemi,
                             double &lat, double &lon )
    {
        double myLat = 0.0, myLon = 0.0;
        long myZone = zone;
        long error = Convert_UTM_To_Geodetic( myZone, northernHemi ? 'N' : 'S',
                                              eastings, northings, &myLat, &myLon );

        lat = myLat;
        lon = myLon;
        return error;
    }

    static long LatLonToUtm( const double lat, const double lon,
                             double &northings, double &eastings,
                             unsigned int &zone, bool &northernHemi )
    {
        long myZone = 0;
        double myNorthings = 0.0, myEastings = 0.0;
        char myHemi = ' ';
        long error = Convert_Geodetic_To_UTM( lat, lon,
                                              &myZone, &myHemi,
                                              &myEastings, &myNorthings );
        zone = myZone;
        eastings = myEastings;
        northings = myNorthings;
        northernHemi = ( myHemi == 'N' );
        return error;
    }

    static const char *UtmErrorToString( long error )
    {
        switch( error )
        {
        case UTM_NO_ERROR:
            return "UTM_NO_ERROR";
        case UTM_LAT_ERROR:
            return "UTM_LAT_ERROR";
        case UTM_LON_ERROR:
            return "UTM_LON_ERROR";
        case UTM_EASTING_ERROR:
            return "UTM_EASTING_ERROR";
        case UTM_NORTHING_ERROR:
            return "UTM_NORTHING_ERROR";
        case UTM_ZONE_ERROR:
            return "UTM_ZONE_ERROR";
        case UTM_HEMISPHERE_ERROR:
            return "UTM_HEMISPHERE_ERROR";
        case UTM_ZONE_OVERRIDE_ERROR:
            return "UTM_ZONE_OVERRIDE_ERROR";
        case UTM_A_ERROR:
            return "UTM_A_ERROR";
        case UTM_INV_F_ERROR:
            return "UTM_INV_F_ERROR";
        default:
            return "(unknown)";
        }
    }

private:

    Wgs84()
    {}
    ;
};

#endif /*DATUM_H_*/