1 // waypoint.cxx -- Class to hold data and return info relating to a waypoint
3 // Written by Curtis Olson, started September 2000.
5 // Copyright (C) 2000 Curtis L. Olson - curt@hfrl.umn.edu
7 // This program is free software; you can redistribute it and/or
8 // modify it under the terms of the GNU General Public License as
9 // published by the Free Software Foundation; either version 2 of the
10 // License, or (at your option) any later version.
12 // This program is distributed in the hope that it will be useful, but
13 // WITHOUT ANY WARRANTY; without even the implied warranty of
14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 // General Public License for more details.
17 // You should have received a copy of the GNU General Public License
18 // along with this program; if not, write to the Free Software
19 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
24 # include <simgear_config.h>
27 #include <simgear/math/polar3d.hxx>
28 #include <simgear/math/sg_geodesy.hxx>
30 #include "waypoint.hxx"
34 SGWayPoint::SGWayPoint( const double lon, const double lat, const double alt,
35 const modetype m, const string s, const string n ) {
46 SGWayPoint::~SGWayPoint() {
50 // Calculate course and distances. For WGS84 and SPHERICAL
51 // coordinates lat, lon, and course are in degrees, alt and distance
52 // are in meters. For CARTESIAN coordinates x = lon, y = lat. Course
53 // is in degrees and distance is in what ever units x and y are in.
54 void SGWayPoint::CourseAndDistance( const double cur_lon,
57 double *course, double *dist ) const {
58 if ( mode == WGS84 ) {
60 geo_inverse_wgs_84( cur_alt, cur_lat, cur_lon, target_lat, target_lon,
61 course, &reverse, dist );
62 } else if ( mode == SPHERICAL ) {
63 Point3D current( cur_lon * SGD_DEGREES_TO_RADIANS, cur_lat * SGD_DEGREES_TO_RADIANS, 0.0 );
64 Point3D target( target_lon * SGD_DEGREES_TO_RADIANS, target_lat * SGD_DEGREES_TO_RADIANS, 0.0 );
65 calc_gc_course_dist( current, target, course, dist );
66 *course = 360.0 - *course * SGD_RADIANS_TO_DEGREES;
67 } else if ( mode == CARTESIAN ) {
68 double dx = target_lon - cur_lon;
69 double dy = target_lat - cur_lat;
70 *course = -atan2( dy, dx ) * SGD_RADIANS_TO_DEGREES - 90;
71 while ( *course < 0 ) {
74 while ( *course > 360.0 ) {
77 *dist = sqrt( dx * dx + dy * dy );
81 // Calculate course and distances between two waypoints
82 void SGWayPoint::CourseAndDistance( const SGWayPoint &wp,
83 double *course, double *dist ) const {
84 CourseAndDistance( wp.get_target_lon(),