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 49

#define GEODESIC_VERSION_PATCH 0

#define GEODESIC_VERSION_NUM(a,b,c) ((((a) * 10000 + (b)) * 100) + (c))

#define GEODESIC_VERSION \
 GEODESIC_VERSION_NUM(GEODESIC_VERSION_MAJOR, \
                      GEODESIC_VERSION_MINOR, \
                      GEODESIC_VERSION_PATCH)

#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 salp1;               
    double calp1;               
    double a13;                 
    double s13;                 
    double b, c2, f1, salp0, calp0, k2,
      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_direct(const struct geod_geodesic* g,
                   double lat1, double lon1, double azi1, double s12,
                   double* plat2, double* plon2, double* pazi2);

  double geod_gendirect(const struct geod_geodesic* g,
                        double lat1, double lon1, double azi1,
                        unsigned flags, double s12_a12,
                        double* plat2, double* plon2, double* pazi2,
                        double* ps12, double* pm12, double* pM12, double* pM21,
                        double* pS12);

  void geod_inverse(const struct geod_geodesic* g,
                    double lat1, double lon1, double lat2, double lon2,
                    double* ps12, double* pazi1, double* pazi2);

  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);

  void geod_lineinit(struct geod_geodesicline* l,
                     const struct geod_geodesic* g,
                     double lat1, double lon1, double azi1, unsigned caps);

  void geod_directline(struct geod_geodesicline* l,
                       const struct geod_geodesic* g,
                       double lat1, double lon1, double azi1, double s12,
                       unsigned caps);

  void geod_gendirectline(struct geod_geodesicline* l,
                          const struct geod_geodesic* g,
                          double lat1, double lon1, double azi1,
                          unsigned flags, double s12_a12,
                          unsigned caps);

  void geod_inverseline(struct geod_geodesicline* l,
                       const struct geod_geodesic* g,
                       double lat1, double lon1, double lat2, double lon2,
                       unsigned caps);

  void geod_position(const struct geod_geodesicline* l, double s12,
                     double* plat2, double* plon2, double* pazi2);

  double geod_genposition(const struct geod_geodesicline* l,
                          unsigned flags, double s12_a12,
                          double* plat2, double* plon2, double* pazi2,
                          double* ps12, double* pm12,
                          double* pM12, double* pM21,
                          double* pS12);

  void geod_setdistance(struct geod_geodesicline* l, double s13);

  void geod_gensetdistance(struct geod_geodesicline* l,
                           unsigned flags, double s13_a13);

  void geod_polygon_init(struct geod_polygon* p, int polylinep);

  void geod_polygon_clear(struct geod_polygon* p);

  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         
  };

  enum geod_flags {
    GEOD_NOFLAGS      = 0U,     
    GEOD_ARCMODE      = 1U<<0,  
    GEOD_LONG_UNROLL  = 1U<<15  
  };

#if defined(__cplusplus)
}
#endif

#endif