simgear/route/waypoint.cxx

   1 // waypoint.cxx -- Class to hold data and return info relating to a waypoint
   2 //
   3 // Written by Curtis Olson, started September 2000.
   4 //
   5 // Copyright (C) 2000  Curtis L. Olson  - curt@hfrl.umn.edu
   6 //
   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.
  11 //
  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.
  16 //
  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.
  20 //
  21 // $Id$
  22
  23 #ifdef HAVE_CONFIG_H
  24 #  include <simgear_config.h>
  25 #endif
  26
  27 #include <simgear/math/polar3d.hxx>
  28 #include <simgear/math/sg_geodesy.hxx>
  29
  30 #include "waypoint.hxx"
  31
  32
  33 // Constructor
  34 SGWayPoint::SGWayPoint( const double lon, const double lat, const double alt,
  35                         const modetype m, const string& s, const string& n ) {
  36     target_lon = lon;
  37     target_lat = lat;
  38     target_alt = alt;
  39     mode = m;
  40     id = s;
  41     name = n;
  42 }
  43
  44
  45 // Destructor
  46 SGWayPoint::~SGWayPoint() {
  47 }
  48
  49
  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,
  55                                     const double cur_lat,
  56                                     const double cur_alt,
  57                                     double *course, double *dist ) const {
  58     if ( mode == WGS84 ) {
  59         double reverse;
  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 ) {
  72             *course += 360.0;
  73         }
  74         while ( *course > 360.0 ) {
  75             *course -= 360.0;
  76         }
  77         *dist = sqrt( dx * dx + dy * dy );
  78     }
  79 }
  80
  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(),
  85                        wp.get_target_lat(),
  86                        wp.get_target_alt(),
  87                        course, dist );
  88 }