src/Navaids/fixlist.cxx

   1 // fixlist.cxx -- fix list management class
   2 //
   3 // Written by Curtis Olson, started April 2000.
   4 //
   5 // Copyright (C) 2000  Curtis L. Olson - curt@flightgear.org
   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., 675 Mass Ave, Cambridge, MA 02139, USA.
  20 //
  21 // $Id$
  22
  23
  24 #ifdef HAVE_CONFIG_H
  25 #  include <config.h>
  26 #endif
  27
  28 #include <simgear/debug/logstream.hxx>
  29 #include <simgear/misc/sgstream.hxx>
  30 #include <simgear/math/sg_geodesy.hxx>
  31
  32 #include "fixlist.hxx"
  33
  34
  35 FGFixList *current_fixlist;
  36
  37
  38 // Constructor
  39 FGFixList::FGFixList( void ) {
  40 }
  41
  42
  43 // Destructor
  44 FGFixList::~FGFixList( void ) {
  45 }
  46
  47
  48 // load the navaids and build the map
  49 bool FGFixList::init( SGPath path ) {
  50
  51     fixlist.erase( fixlist.begin(), fixlist.end() );
  52
  53     sg_gzifstream in( path.str() );
  54     if ( !in.is_open() ) {
  55         SG_LOG( SG_GENERAL, SG_ALERT, "Cannot open file: " << path.str() );
  56         exit(-1);
  57     }
  58
  59     // toss the first two lines of the file
  60     in >> skipeol;
  61     in >> skipeol;
  62
  63     // read in each remaining line of the file
  64
  65 #ifdef __MWERKS__
  66     char c = 0;
  67     while ( in.get(c) && c != '\0' ) {
  68         in.putback(c);
  69 #else
  70     while ( ! in.eof() ) {
  71 #endif
  72
  73         FGFix fix;
  74         in >> fix;
  75         if ( fix.get_lat() > 95.0 ) {
  76             break;
  77         }
  78
  79         /* cout << "ident=" << fix.get_ident()
  80              << ", lat=" << fix.get_lat()
  81              << ", lon=" << fix.get_lon() << endl; */
  82
  83         fixlist[fix.get_ident()] = fix;
  84         in >> skipcomment;
  85     }
  86
  87     return true;
  88 }
  89
  90
  91 // query the database for the specified fix, lon and lat are in
  92 // degrees, elev is in meters
  93 bool FGFixList::query( const string& ident, FGFix *fix ) {
  94     *fix = fixlist[ident];
  95     if ( ! fix->get_ident().empty() ) {
  96         return true;
  97     } else {
  98         return false;
  99     }
 100 }
 101
 102
 103 // query the database for the specified fix, lon and lat are in
 104 // degrees, elev is in meters
 105 bool FGFixList::query_and_offset( const string& ident, double lon, double lat,
 106                                   double elev, FGFix *fix, double *heading,
 107                                   double *dist )
 108 {
 109     *fix = fixlist[ident];
 110     if ( fix->get_ident().empty() ) {
 111         return false;
 112     }
 113
 114     double az1, az2, s;
 115     geo_inverse_wgs_84( elev, lat, lon,
 116                         fix->get_lat(), fix->get_lon(),
 117                         &az1, &az2, &s );
 118     // cout << "  dist = " << s << endl;
 119     *heading = az2;
 120     *dist = s;
 121     return true;
 122 }