1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
// Boost.Geometry
// Copyright (c) 2018 Adeel Ahmad, Islamabad, Pakistan.
// Contributed and/or modified by Adeel Ahmad, as part of Google Summer of Code 2018 program.
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
// Formula example - Show how to use Karney's direct method.
#include <boost/geometry.hpp>
#include <boost/geometry/formulas/karney_direct.hpp>
using namespace boost::geometry;
int main()
{
double lon1_deg = 0.;
double lat1_deg = 73.114273316483;
double distance_m = 19992866.6147806;
double azi12_deg = 78.154765899661;
// Create an alias of the formula.
typedef formula::karney_direct<double, true, true, true, true, 8> karney_direct;
// Structure to hold the resulting values.
formula::result_direct<double> result;
// WGS-84 spheroid.
srs::spheroid<double> spheroid(6378137.0, 6356752.3142451793);
result = karney_direct::apply(lon1_deg, lat1_deg, distance_m, azi12_deg, spheroid);
return 0;
}
|