planimeter.c

Go to the documentation of this file.

#include <stdio.h>
#include "geodesic.h"

#if defined(_MSC_VER)
/* Squelch warnings about scanf */
#  pragma warning (disable: 4996)
#endif

int main() {
  double a = 6378137, f = 1/298.257223563; /* WGS84 */
  double lat, lon, A, P;
  int n;
  struct geod_geodesic g;
  struct geod_polygon p;
  geod_init(&g, a, f);
  geod_polygon_init(&p, 0);

  while (scanf("%lf %lf", &lat, &lon) == 2)
    geod_polygon_addpoint(&g, &p, lat, lon);
  n = geod_polygon_compute(&g, &p, 0, 1, &A, &P);
  printf("%d %.8f %.3f\n", n, P, A);
  return 0;
}