16 template <
class GeodType>
17 int PolygonAreaT<GeodType>::transit(
real lon1,
real lon2) {
29 lon12 > 0 && ((lon1 < 0 && lon2 >= 0) ||
31 (lon1 > 0 && lon2 == 0)) ? 1 :
33 (lon12 < 0 && lon1 >= 0 && lon2 < 0 ? -1 : 0);
44 template <
class GeodType>
45 int PolygonAreaT<GeodType>::transitdirect(
real lon1,
real lon2) {
60 return ( (lon2 >= 0 && lon2 <
Math::td ? 0 : 1) -
61 (lon1 >= 0 && lon1 <
Math::td ? 0 : 1) );
64 template <
class GeodType>
71 _earth.GenInverse(_lat1, _lon1, lat, lon, _mask,
72 s12, t, t, t, t, t, S12);
76 _crossings += transit(_lon1, lon);
78 _lat1 = lat; _lon1 = lon;
83 template <
class GeodType>
86 real lat, lon, S12, t;
87 _earth.GenDirect(_lat1, _lon1, azi,
false, s, _mask,
88 lat, lon, t, t, t, t, t, S12);
92 _crossings += transitdirect(_lon1, lon);
94 _lat1 = lat; _lon1 = lon;
99 template <
class GeodType>
101 real& perimeter, real& area)
const
111 perimeter = _perimetersum();
114 _earth.GenInverse(_lat1, _lon1, _lat0, _lon0, _mask,
115 s12, t, t, t, t, t, S12);
116 perimeter = _perimetersum(s12);
119 int crossings = _crossings + transit(_lon1, _lon0);
120 AreaReduce(tempsum, crossings, reverse, sign);
121 area =
real(0) + tempsum();
125 template <
class GeodType>
127 bool reverse,
bool sign,
128 real& perimeter, real& area)
const
136 perimeter = _perimetersum();
137 real tempsum = _polyline ? 0 : _areasum();
138 int crossings = _crossings;
139 unsigned num = _num + 1;
140 for (
int i = 0; i < (_polyline ? 1 : 2); ++i) {
142 _earth.GenInverse(i == 0 ? _lat1 : lat, i == 0 ? _lon1 : lon,
143 i != 0 ? _lat0 : lat, i != 0 ? _lon0 : lon,
144 _mask, s12, t, t, t, t, t, S12);
148 crossings += transit(i == 0 ? _lon1 : lon,
149 i != 0 ? _lon0 : lon);
156 AreaReduce(tempsum, crossings, reverse, sign);
157 area =
real(0) + tempsum;
161 template <
class GeodType>
163 bool reverse,
bool sign,
164 real& perimeter, real& area)
const
172 unsigned num = _num + 1;
173 perimeter = _perimetersum() + s;
177 real tempsum = _areasum();
178 int crossings = _crossings;
180 real lat, lon, s12, S12, t;
181 _earth.GenDirect(_lat1, _lon1, azi,
false, s, _mask,
182 lat, lon, t, t, t, t, t, S12);
184 crossings += transitdirect(_lon1, lon);
185 _earth.GenInverse(lat, lon, _lat0, _lon0, _mask,
186 s12, t, t, t, t, t, S12);
189 crossings += transit(lon, _lon0);
192 AreaReduce(tempsum, crossings, reverse, sign);
193 area =
real(0) + tempsum;
197 template <
class GeodType>
198 template <
typename T>
200 bool reverse,
bool sign)
const {
202 if (crossings & 1) area += (area < 0 ? 1 : -1) * _area0/2;
205 if (!reverse) area *= -1;
210 else if (area <= -_area0/2)
#define GEOGRAPHICLIB_EXPORT
GeographicLib::Math::real real
Header for GeographicLib::PolygonAreaT class.
static T AngNormalize(T x)
static T AngDiff(T x, T y, T &e)
void AddPoint(real lat, real lon)
unsigned Compute(bool reverse, bool sign, real &perimeter, real &area) const
unsigned TestPoint(real lat, real lon, bool reverse, bool sign, real &perimeter, real &area) const
void AddEdge(real azi, real s)
unsigned TestEdge(real azi, real s, bool reverse, bool sign, real &perimeter, real &area) const
Namespace for GeographicLib.