geodesic.h

Go to the documentation of this file.

 
#if !defined(GEODESIC_H)
#define GEODESIC_H 1
 
#define GEODESIC_VERSION_MAJOR 1
#define GEODESIC_VERSION_MINOR 32
#define GEODESIC_VERSION_PATCH 0
 
#if defined(__cplusplus)
extern "C" {
#endif
 
struct geod_geodesic {
    double a; 
    double f; 
    double f1, e2, ep2, n, b, c2, etol2;
    double A3x[6], C3x[15], C4x[21];
};
 
struct geod_geodesicline {
    double lat1; 
    double lon1; 
    double azi1; 
    double a;    
    double f;    
    double b, c2, f1, salp0, calp0, k2, salp1, calp1, ssig1, csig1, dn1, stau1, ctau1, somg1, comg1, A1m1, A2m1, A3c,
        B11, B21, B31, A4, B41;
    double C1a[6 + 1], C1pa[6 + 1], C2a[6 + 1], C3a[6], C4a[6];
    unsigned caps; 
};
 
struct geod_polygon {
    double lat; 
    double lon; 
    double lat0;
    double lon0;
    double A[2];
    double P[2];
    int polyline;
    int crossings;
    unsigned num; 
};
 
void geod_init(struct geod_geodesic* g, double a, double f);
 
void geod_lineinit(struct geod_geodesicline* l, const struct geod_geodesic* g, double lat1, double lon1, double azi1,
                   unsigned caps);
 
void geod_direct(const struct geod_geodesic* g, double lat1, double lon1, double azi1, double s12, double* plat2,
                 double* plon2, double* pazi2);
 
void geod_inverse(const struct geod_geodesic* g, double lat1, double lon1, double lat2, double lon2, double* ps12,
                  double* pazi1, double* pazi2);
 
void geod_position(const struct geod_geodesicline* l, double s12, double* plat2, double* plon2, double* pazi2);
 
double geod_gendirect(const struct geod_geodesic* g, double lat1, double lon1, double azi1, int arcmode, double s12_a12,
                      double* plat2, double* plon2, double* pazi2, double* ps12, double* pm12, double* pM12,
                      double* pM21, double* pS12);
 
double geod_geninverse(const struct geod_geodesic* g, double lat1, double lon1, double lat2, double lon2, double* ps12,
                       double* pazi1, double* pazi2, double* pm12, double* pM12, double* pM21, double* pS12);
 
double geod_genposition(const struct geod_geodesicline* l, int arcmode, double s12_a12, double* plat2, double* plon2,
                        double* pazi2, double* ps12, double* pm12, double* pM12, double* pM21, double* pS12);
 
void geod_polygon_init(struct geod_polygon* p, int polylinep);
 
void geod_polygon_addpoint(const struct geod_geodesic* g, struct geod_polygon* p, double lat, double lon);
 
void geod_polygon_addedge(const struct geod_geodesic* g, struct geod_polygon* p, double azi, double s);
 
unsigned geod_polygon_compute(const struct geod_geodesic* g, const struct geod_polygon* p, int reverse, int sign,
                              double* pA, double* pP);
 
unsigned geod_polygon_testpoint(const struct geod_geodesic* g, const struct geod_polygon* p, double lat, double lon,
                                int reverse, int sign, double* pA, double* pP);
 
unsigned geod_polygon_testedge(const struct geod_geodesic* g, const struct geod_polygon* p, double azi, double s,
                               int reverse, int sign, double* pA, double* pP);
 
void geod_polygonarea(const struct geod_geodesic* g, double lats[], double lons[], int n, double* pA, double* pP);
 
enum geod_mask {
    GEOD_NONE = 0U,                                    
    GEOD_LATITUDE = 1U << 7 | 0U,                      
    GEOD_LONGITUDE = 1U << 8 | 1U << 3,                
    GEOD_AZIMUTH = 1U << 9 | 0U,                       
    GEOD_DISTANCE = 1U << 10 | 1U << 0,                
    GEOD_DISTANCE_IN = 1U << 11 | 1U << 0 | 1U << 1,   
    GEOD_REDUCEDLENGTH = 1U << 12 | 1U << 0 | 1U << 2, 
    GEOD_GEODESICSCALE = 1U << 13 | 1U << 0 | 1U << 2, 
    GEOD_AREA = 1U << 14 | 1U << 4,                    
    GEOD_ALL = 0x7F80U | 0x1FU                         
};
 
#if defined(__cplusplus)
}
#endif
 
#endif