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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
|
/**
* \file NETGeographicLib/Rhumb.cpp
* \brief Implementation for NETGeographicLib::Rhumb and NETGeographicLib::RhumbLine class
*
* NETGeographicLib is copyright (c) Scott Heiman (2013)
* GeographicLib is Copyright (c) Charles Karney (2010-2012)
* <charles@karney.com> and licensed under the MIT/X11 License.
* For more information, see
* http://geographiclib.sourceforge.net/
**********************************************************************/
#include "stdafx.h"
#include "GeographicLib/Rhumb.hpp"
#include "Rhumb.h"
#include "NETGeographicLib.h"
using namespace NETGeographicLib;
//*****************************************************************************
Rhumb::!Rhumb(void)
{
if ( m_pRhumb != NULL )
{
delete m_pRhumb;
m_pRhumb = NULL;
}
}
//*****************************************************************************
Rhumb::Rhumb(double a, double f, bool exact)
{
try
{
m_pRhumb = new GeographicLib::Rhumb( a, f, exact );
}
catch ( GeographicLib::GeographicErr& err )
{
throw gcnew GeographicErr( err.what() );
}
catch ( std::bad_alloc )
{
throw gcnew System::Exception("Failed to allocate memory for a Rhumb.");
}
}
//*****************************************************************************
void Rhumb::Direct(double lat1, double lon1, double azi12, double s12,
[System::Runtime::InteropServices::Out] double% lat2,
[System::Runtime::InteropServices::Out] double% lon2)
{
double ilat2, ilon2;
m_pRhumb->Direct( lat1, lon1, azi12, s12, ilat2, ilon2 );
lat2 = ilat2;
lon2 = ilon2;
}
//*****************************************************************************
void Rhumb::Inverse(double lat1, double lon1, double lat2, double lon2,
[System::Runtime::InteropServices::Out] double% s12,
[System::Runtime::InteropServices::Out] double% azi12)
{
double is12, iazi12;
m_pRhumb->Inverse( lat1, lon1, lat2, lon2, is12, iazi12 );
s12 = is12;
azi12 = iazi12;
}
//*****************************************************************************
RhumbLine^ Rhumb::Line(double lat1, double lon1, double azi12)
{
return gcnew RhumbLine( new GeographicLib::RhumbLine(m_pRhumb->Line( lat1, lon1, azi12 )) );
}
//*****************************************************************************
double Rhumb::MajorRadius::get() { return m_pRhumb->MajorRadius(); }
//*****************************************************************************
double Rhumb::Flattening::get() { return m_pRhumb->Flattening(); }
//*****************************************************************************
Rhumb^ Rhumb::WGS84()
{
return gcnew Rhumb( GeographicLib::Constants::WGS84_a(),
GeographicLib::Constants::WGS84_f(), false );
}
//*****************************************************************************
// RhumbLine functions
//*****************************************************************************
RhumbLine::!RhumbLine(void)
{
if ( m_pRhumbLine != NULL )
{
delete m_pRhumbLine;
m_pRhumbLine = NULL;
}
}
//*****************************************************************************
RhumbLine::RhumbLine( GeographicLib::RhumbLine* pRhumbLine )
{
if ( pRhumbLine == NULL )
throw gcnew System::Exception("Invalid pointer in RhumbLine constructor.");
m_pRhumbLine = pRhumbLine;
}
//*****************************************************************************
void RhumbLine::Position(double s12,
[System::Runtime::InteropServices::Out] double% lat2,
[System::Runtime::InteropServices::Out] double% lon2)
{
double ilat2, ilon2;
m_pRhumbLine->Position( s12, ilat2, ilon2 );
lat2 = ilat2;
lon2 = ilon2;
}
//*****************************************************************************
double RhumbLine::Latitude::get()
{
return m_pRhumbLine->Latitude();
}
//*****************************************************************************
double RhumbLine::Longitude::get()
{
return m_pRhumbLine->Longitude();
}
//*****************************************************************************
double RhumbLine::Azimuth::get()
{
return m_pRhumbLine->Azimuth();
}
//*****************************************************************************
double RhumbLine::MajorRadius::get()
{
return m_pRhumbLine->MajorRadius();
}
//*****************************************************************************
double RhumbLine::Flattening::get()
{
return m_pRhumbLine->Flattening();
}
|
