LatLong-UTMconversion.h

Go to the documentation of this file.

/*
 * $Id: LatLong-UTMconversion.h 17606 2014-09-26 01:09:51Z predoehl $
 */

#ifndef LATLONGCONV_H_INCLUDED_UOFARIZONA_VISION
#define LATLONGCONV_H_INCLUDED_UOFARIZONA_VISION

#include <topo_cpp/layer.h>
#include <vector>

namespace kjb
{
namespace TopoFusion
{

void LLtoUTM(
    int ReferenceEllipsoid,
    const double Lat,
    const double Long,
    pt& utm
);

#if 0
void UTMtoLL(
    int ReferenceEllipsoid,
    const pt& utm,
    double& Lat,
    double& Long
);
#endif

void UTMtoGPXLL(
    int ReferenceEllipsoid,
    const pt& utm,
    double& Lat,
    double& Long
);


void utm_to_lat_long(
    int ReferenceEllipsoid,
    const pt& utm,
    double& Lat,
    double& Long
);


void utm_to_lat_long(
    int,
    const std::vector<pt>&,
    std::vector<double>*,
    std::vector<double>* 
);


double getNewEasting(
    const pt& utm,
    int new_zone
);


char zone_of(int ReferenceEllipspoid, const pt&);


#if 0
struct Debug_spy_utm_ll_conversion {
    double eccPrimeSquared, e1, mu, phi1Rad, N1, T1, C1, R1, D, Lat, Long;
};

extern Debug_spy_utm_ll_conversion spy[2];
#endif

struct Ellipsoid
{
    int id;                     
    const char* ellipsoidName;  
    double EquatorialRadius;    
    double eccentricitySquared; 
    double oneOverFlattening;   

    Ellipsoid(int Id, const char* name, double radius, double ecc, double flat)
    :   id(Id),
        ellipsoidName(name),
        EquatorialRadius(radius),
        eccentricitySquared(ecc),
        oneOverFlattening(flat)
    {}
};

enum ELLIPSOID_ID {
    ELLIPSOID_Airy = 1,
    ELLIPSOID_Australian_National,
    ELLIPSOID_Bessel_1841,
    ELLIPSOID_Bessel_1841_Nambia,
    ELLIPSOID_Clarke_1866,              
    ELLIPSOID_Clarke_1880,
    ELLIPSOID_Everest,
    ELLIPSOID_Fischer_1960_Mercury,
    ELLIPSOID_Fischer_1968,
    ELLIPSOID_GRS_1967,
    ELLIPSOID_GRS_1980,                 
    ELLIPSOID_Helmert_1906,
    ELLIPSOID_Hough,
    ELLIPSOID_International,
    ELLIPSOID_Krassovsky,
    ELLIPSOID_Modified_Airy,
    ELLIPSOID_Modified_Everest,
    ELLIPSOID_Modified_Fischer_1960,
    ELLIPSOID_South_American_1969,
    ELLIPSOID_WGS_60,
    ELLIPSOID_WGS_66,
    ELLIPSOID_WGS_72,
    ELLIPSOID_WGS_84,                   
    ELLIPSOID_END_OF_LIST,

    ELLIPSOID_START_OF_LIST = ELLIPSOID_Airy
};


int validate_ellipsoid_table();


#if 0
struct datum
{
    int id;                     
    const char* dName;          
    int ellipsoid;              
    int dx, dy, dz;             // offsets relative to WGS-84

    datum( int Id, const char* name, int el, int Dx, int Dy, int Dz )
    :   id( Id ),
        dName( name ),
        ellipsoid( el ),
        dx( Dx ),
        dy( Dy ),
        dz( Dz )
    {}
};
#endif


} // end namespace TopoFusion
} // end namespace kjb


#endif