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,
32 const modetype m, const string s ) {
40 SGWayPoint::SGWayPoint() {
41 SGWayPoint( 0.0, 0.0, WGS84, "" );
46 SGWayPoint::~SGWayPoint() {
50 // Calculate course and distances. For WGS84 and SPHERICAL
51 // coordinates lat, lon, and course are in degrees and distance in
52 // meters. For CARTESIAN coordinates x = lon, y = lat. Course is in
53 // degrees and distance is in what ever units x and y are in.
54 void SGWayPoint::CourseAndDistance( const double cur_lon, const double cur_lat,
55 double *course, double *distance ) {
56 if ( mode == WGS84 ) {
58 geo_inverse_wgs_84( 0.0, cur_lat, cur_lon, target_lat, target_lon,
59 course, &reverse, distance );
60 } else if ( mode == SPHERICAL ) {
61 Point3D current( cur_lon * DEG_TO_RAD, cur_lat * DEG_TO_RAD, 0.0 );
62 Point3D target( target_lon * DEG_TO_RAD, target_lat * DEG_TO_RAD, 0.0 );
63 calc_gc_course_dist( current, target, course, distance );
64 *course = 360.0 - *course * RAD_TO_DEG;
65 } else if ( mode == CARTESIAN ) {
66 double dx = target_lon - cur_lon;
67 double dy = target_lat - cur_lat;
68 *course = -atan2( dy, dx ) * RAD_TO_DEG - 90;
69 while ( *course < 0 ) {
72 while ( *course > 360.0 ) {
75 *distance = sqrt( dx * dx + dy * dy );