14# pragma warning (disable: 5055)
21 template<
class GeodType>
22 int PolygonAreaT<GeodType>::transit(
real lon1,
real lon2) {
34 lon12 > 0 && ((lon1 < 0 && lon2 >= 0) ||
36 (lon1 > 0 && lon2 == 0)) ? 1 :
38 (lon12 < 0 && lon1 >= 0 && lon2 < 0 ? -1 : 0);
49 template<
class GeodType>
50 int PolygonAreaT<GeodType>::transitdirect(
real lon1,
real lon2) {
64 return ( (lon2 >= 0 && lon2 <
Math::td ? 0 : 1) -
65 (lon1 >= 0 && lon1 <
Math::td ? 0 : 1) );
68 template<
class GeodType>
75 _earth.GenInverse(_lat1, _lon1, lat, lon, _mask,
76 s12, t, t, t, t, t, S12);
80 _crossings += transit(_lon1, lon);
82 _lat1 = lat; _lon1 = lon;
87 template<
class GeodType>
90 real lat, lon, S12, t;
91 _earth.GenDirect(_lat1, _lon1, azi,
false, s, _mask,
92 lat, lon, t, t, t, t, t, S12);
96 _crossings += transitdirect(_lon1, lon);
98 _lat1 = lat; _lon1 = lon;
103 template<
class GeodType>
105 real& perimeter, real& area)
const
115 perimeter = _perimetersum();
118 _earth.GenInverse(_lat1, _lon1, _lat0, _lon0, _mask,
119 s12, t, t, t, t, t, S12);
120 perimeter = _perimetersum(s12);
123 int crossings = _crossings + transit(_lon1, _lon0);
124 AreaReduce(tempsum, crossings, reverse, sign);
125 area =
real(0) + tempsum();
129 template<
class GeodType>
131 bool reverse,
bool sign,
132 real& perimeter, real& area)
const
140 perimeter = _perimetersum();
141 real tempsum = _polyline ? 0 : _areasum();
142 int crossings = _crossings;
143 unsigned num = _num + 1;
144 for (
int i = 0; i < (_polyline ? 1 : 2); ++i) {
146 _earth.GenInverse(i == 0 ? _lat1 : lat, i == 0 ? _lon1 : lon,
147 i != 0 ? _lat0 : lat, i != 0 ? _lon0 : lon,
148 _mask, s12, t, t, t, t, t, S12);
152 crossings += transit(i == 0 ? _lon1 : lon,
153 i != 0 ? _lon0 : lon);
160 AreaReduce(tempsum, crossings, reverse, sign);
161 area =
real(0) + tempsum;
165 template<
class GeodType>
167 bool reverse,
bool sign,
168 real& perimeter, real& area)
const
176 unsigned num = _num + 1;
177 perimeter = _perimetersum() + s;
181 real tempsum = _areasum();
182 int crossings = _crossings;
184 real lat, lon, s12, S12, t;
185 _earth.GenDirect(_lat1, _lon1, azi,
false, s, _mask,
186 lat, lon, t, t, t, t, t, S12);
188 crossings += transitdirect(_lon1, lon);
189 _earth.GenInverse(lat, lon, _lat0, _lon0, _mask,
190 s12, t, t, t, t, t, S12);
193 crossings += transit(lon, _lon0);
196 AreaReduce(tempsum, crossings, reverse, sign);
197 area =
real(0) + tempsum;
201 template<
class GeodType>
204 bool reverse,
bool sign)
const {
206 if (crossings & 1) area += (area < 0 ? 1 : -1) * _area0/2;
209 if (!reverse) area *= -1;
214 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.