examples/example-PolarStereographic.cpp
Go to the documentation of this file.
1 // Example of using the GeographicLib::PolarStereographic class
2 
3 #include <iostream>
4 #include <iomanip>
5 #include <exception>
7 
8 using namespace std;
9 using namespace GeographicLib;
10 
11 int main() {
12  try {
13  PolarStereographic proj(Constants::WGS84_a(), Constants::WGS84_f(),
14  Constants::UPS_k0());
15  // Alternatively:
16  // const PolarStereographic& proj = PolarStereographic::UPS();
17  bool northp = true;
18  {
19  // Sample forward calculation
20  double lat = 61.2, lon = -149.9; // Anchorage
21  double x, y;
22  proj.Forward(northp, lat, lon, x, y);
23  cout << x << " " << y << "\n";
24  }
25  {
26  // Sample reverse calculation
27  double x = -1637e3, y = 2824e3;
28  double lat, lon;
29  proj.Reverse(northp, x, y, lat, lon);
30  cout << lat << " " << lon << "\n";
31  }
32  }
33  catch (const exception& e) {
34  cerr << "Caught exception: " << e.what() << "\n";
35  return 1;
36  }
37 }
Scalar * y
static const double lat
Definition: Half.h:150
void Forward(bool northp, real lat, real lon, real &x, real &y, real &gamma, real &k) const
void Reverse(bool northp, real x, real y, real &lat, real &lon, real &gamma, real &k) const
Namespace for GeographicLib.
Array< double, 1, 3 > e(1./3., 0.5, 2.)
Polar stereographic projection.
static const double lon
Header for GeographicLib::PolarStereographic class.
set noclip points set clip one set noclip two set bar set border lt lw set xdata set ydata set zdata set x2data set y2data set boxwidth set dummy x


gtsam
Author(s):
autogenerated on Sat May 8 2021 02:42:02