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., 675 Mass Ave, Cambridge, MA 02139, USA.
24 #include <simgear/math/polar3d.hxx>
25 #include <simgear/math/sg_geodesy.hxx>
27 #include "waypoint.hxx"
31 SGWayPoint::SGWayPoint( const double lon, const double lat, const double alt,
32 const modetype m, const string s ) {
42 SGWayPoint::~SGWayPoint() {
46 // Calculate course and distances. For WGS84 and SPHERICAL
47 // coordinates lat, lon, and course are in degrees, alt and distance
48 // are in meters. For CARTESIAN coordinates x = lon, y = lat. Course
49 // is in degrees and distance is in what ever units x and y are in.
50 void SGWayPoint::CourseAndDistance( const double cur_lon,
53 double *course, double *distance ) const {
54 if ( mode == WGS84 ) {
56 geo_inverse_wgs_84( cur_alt, cur_lat, cur_lon, target_lat, target_lon,
57 course, &reverse, distance );
58 } else if ( mode == SPHERICAL ) {
59 Point3D current( cur_lon * SGD_DEGREES_TO_RADIANS, cur_lat * SGD_DEGREES_TO_RADIANS, 0.0 );
60 Point3D target( target_lon * SGD_DEGREES_TO_RADIANS, target_lat * SGD_DEGREES_TO_RADIANS, 0.0 );
61 calc_gc_course_dist( current, target, course, distance );
62 *course = 360.0 - *course * SGD_RADIANS_TO_DEGREES;
63 } else if ( mode == CARTESIAN ) {
64 double dx = target_lon - cur_lon;
65 double dy = target_lat - cur_lat;
66 *course = -atan2( dy, dx ) * SGD_RADIANS_TO_DEGREES - 90;
67 while ( *course < 0 ) {
70 while ( *course > 360.0 ) {
73 *distance = sqrt( dx * dx + dy * dy );
77 // Calculate course and distances between two waypoints
78 void SGWayPoint::CourseAndDistance( const SGWayPoint &wp,
79 double *course, double *distance ) const {
80 CourseAndDistance( wp.get_target_lon(),