/// Set the geodetic elevation from the argument given in feet
void setElevationFt(double elevation);
- // Compare two geodetic positions for equality
+ /// Compare two geodetic positions for equality
bool operator == ( const SGGeod & other ) const;
+ /// check the Geod contains sane values (finite, inside appropriate
+ /// ranges for lat/lon)
+ bool isValid() const;
private:
/// This one is private since construction is not unique if you do
/// not know the units of the arguments. Use the factory methods for
_elevation == other._elevation;
}
+inline
+bool
+SGGeod::isValid() const
+{
+ if (SGMiscd::isNaN(_lon))
+ return false;
+ if (SGMiscd::isNaN(_lat))
+ return false;
+#ifdef SG_GEOD_NATIVE_DEGREE
+ return (_lon >= -180.0) && (_lon <= 180.0) &&
+ (_lat >= -90.0) && (_lat <= 90.0);
+#else
+ return (_lon >= -SGD_PI) && (_lon <= SGD_PI) &&
+ (_lat >= -SGD_PI_2) && (_lat <= SGD_PI_2);
+#endif
+}
+
/// Output to an ostream
template<typename char_type, typename traits_type>
inline