1 // runways.hxx -- a simple class to manage airport runway info
3 // Written by Curtis Olson, started August 2000.
5 // Copyright (C) 2000 Curtis L. Olson - curt@flightgear.org
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.
28 #include <math.h> // fabs()
29 #include <stdio.h> // sprintf()
31 #include <simgear/compiler.h>
33 #include <simgear/debug/logstream.hxx>
34 #include <simgear/misc/sgstream.hxx>
37 #include STL_FUNCTIONAL
38 #include STL_ALGORITHM
40 #include "runways.hxx"
42 SG_USING_NAMESPACE(std);
45 FGRunway::FGRunway() {
49 FGRunway::~FGRunway() {
53 FGRunways::FGRunways( const string& file ) {
54 // open the specified database readonly
55 storage = new c4_Storage( file.c_str(), false );
57 if ( !storage->Strategy().IsValid() ) {
58 SG_LOG( SG_GENERAL, SG_ALERT, "Cannot open file: " << file );
62 vRunway = new c4_View;
64 storage->GetAs("runway[ID:S,Rwy:S,Longitude:F,Latitude:F,Heading:F,Length:F,Width:F,SurfaceFlags:S,End1Flags:S,End2Flags:S]");
70 // search for the specified apt id
71 bool FGRunways::search( const string& aptid, FGRunway* r ) {
72 c4_StringProp pID ("ID");
73 c4_StringProp pRwy ("Rwy");
74 c4_FloatProp pLon ("Longitude");
75 c4_FloatProp pLat ("Latitude");
76 c4_FloatProp pHdg ("Heading");
77 c4_FloatProp pLen ("Length");
78 c4_FloatProp pWid ("Width");
79 c4_StringProp pSurf ("SurfaceFlags");
80 c4_StringProp pEnd1 ("End1Flags");
81 c4_StringProp pEnd2 ("End2Flags");
83 int index = vRunway->Find(pID[aptid.c_str()]);
84 // cout << "index = " << index << endl;
90 next_index = index + 1;
92 c4_RowRef row = vRunway->GetAt(index);
94 r->id = (const char *) pID(row);
95 r->rwy_no = (const char *) pRwy(row);
96 r->lon = (double) pLon(row);
97 r->lat = (double) pLat(row);
98 r->heading = (double) pHdg(row);
99 r->length = (double) pLen(row);
100 r->width = (double) pWid(row);
101 r->surface_flags = (const char *) pSurf(row);
102 r->end1_flags = (const char *) pEnd1(row);
103 r->end2_flags = (const char *) pEnd2(row);
109 // search for the specified apt id and runway no
110 bool FGRunways::search( const string& aptid, const string& rwyno, FGRunway* r )
112 c4_StringProp pID ("ID");
113 c4_StringProp pRwy ("Rwy");
114 c4_FloatProp pLon ("Longitude");
115 c4_FloatProp pLat ("Latitude");
116 c4_FloatProp pHdg ("Heading");
117 c4_FloatProp pLen ("Length");
118 c4_FloatProp pWid ("Width");
119 c4_StringProp pSurf ("SurfaceFlags");
120 c4_StringProp pEnd1 ("End1Flags");
121 c4_StringProp pEnd2 ("End2Flags");
123 int index = vRunway->Find(pID[aptid.c_str()]);
124 // cout << "index = " << index << endl;
130 c4_RowRef row = vRunway->GetAt(index);
131 string rowid = (const char *) pID(row);
132 string rowrwyno = (const char *) pRwy(row);
133 while ( rowid == aptid ) {
134 next_index = index + 1;
136 if ( rowrwyno == rwyno ) {
137 r->id = (const char *) pID(row);
138 r->rwy_no = (const char *) pRwy(row);
139 r->lon = (double) pLon(row);
140 r->lat = (double) pLat(row);
141 r->heading = (double) pHdg(row);
142 r->length = (double) pLen(row);
143 r->width = (double) pWid(row);
144 r->surface_flags = (const char *) pSurf(row);
145 r->end1_flags = (const char *) pEnd1(row);
146 r->end2_flags = (const char *) pEnd2(row);
152 c4_RowRef row = vRunway->GetAt(index);
153 string rowid = (const char *) pID(row);
154 string rowrwyno = (const char *) pRwy(row);
161 FGRunway FGRunways::search( const string& aptid ) {
168 // search for the specified id
169 bool FGRunways::next( FGRunway* r ) {
170 c4_StringProp pID ("ID");
171 c4_StringProp pRwy ("Rwy");
172 c4_FloatProp pLon ("Longitude");
173 c4_FloatProp pLat ("Latitude");
174 c4_FloatProp pHdg ("Heading");
175 c4_FloatProp pLen ("Length");
176 c4_FloatProp pWid ("Width");
177 c4_StringProp pSurf ("SurfaceFlags");
178 c4_StringProp pEnd1 ("End1Flags");
179 c4_StringProp pEnd2 ("End2Flags");
181 int size = vRunway->GetSize();
182 // cout << "total records = " << size << endl;
184 int index = next_index;
185 // cout << "index = " << index << endl;
187 if ( index == -1 || index >= size ) {
191 next_index = index + 1;
193 c4_RowRef row = vRunway->GetAt(index);
194 r->id = (const char *) pID(row);
195 r->rwy_no = (const char *) pRwy(row);
196 r->lon = (double) pLon(row);
197 r->lat = (double) pLat(row);
198 r->heading = (double) pHdg(row);
199 r->length = (double) pLen(row);
200 r->width = (double) pWid(row);
201 r->surface_flags = (const char *) pSurf(row);
202 r->end1_flags = (const char *) pEnd1(row);
203 r->end2_flags = (const char *) pEnd2(row);
209 // Return the runway closest to a given heading
210 bool FGRunways::search( const string& aptid, const int tgt_hdg,
215 double found_dir = 0.0;
217 if ( !search( aptid, &r ) ) {
218 SG_LOG( SG_GENERAL, SG_ALERT,
219 "Failed to find " << aptid << " in database." );
224 double min_diff = 360.0;
226 while ( r.id == aptid ) {
228 diff = tgt_hdg - r.heading;
229 while ( diff < -180.0 ) { diff += 360.0; }
230 while ( diff > 180.0 ) { diff -= 360.0; }
232 // SG_LOG( SG_GENERAL, SG_INFO,
233 // "Runway " << r.rwy_no << " heading = " << r.heading <<
234 // " diff = " << diff );
235 if ( diff < min_diff ) {
243 diff = tgt_hdg - r.heading - 180.0;
244 while ( diff < -180.0 ) { diff += 360.0; }
245 while ( diff > 180.0 ) { diff -= 360.0; }
247 // SG_LOG( SG_GENERAL, SG_INFO,
248 // "Runway -" << r.rwy_no << " heading = " <<
249 // r.heading + 180.0 <<
250 // " diff = " << diff );
251 if ( diff < min_diff ) {
260 // SG_LOG( SG_GENERAL, SG_INFO, "closest runway = " << runway->rwy_no
261 // << " + " << found_dir );
264 double heading = runway->heading + found_dir;
265 while ( heading >= 360.0 ) { heading -= 360.0; }
266 runway->heading = heading;
267 //cout << runway->heading << " " << runway->rwy_no
268 // << " " << runway->id << endl;
273 // Return the runway number of the runway closest to a given heading
274 string FGRunways::search( const string& aptid, const int tgt_hdg ) {
277 double found_dir = 0.0;
279 if ( !search( aptid, &r ) ) {
280 SG_LOG( SG_GENERAL, SG_ALERT,
281 "Failed to find " << aptid << " in database." );
286 double min_diff = 360.0;
288 while ( r.id == aptid ) {
290 diff = tgt_hdg - r.heading;
291 while ( diff < -180.0 ) { diff += 360.0; }
292 while ( diff > 180.0 ) { diff -= 360.0; }
294 // SG_LOG( SG_GENERAL, SG_INFO,
295 // "Runway " << r.rwy_no << " heading = " << r.heading <<
296 // " diff = " << diff );
297 if ( diff < min_diff ) {
304 diff = tgt_hdg - r.heading - 180.0;
305 while ( diff < -180.0 ) { diff += 360.0; }
306 while ( diff > 180.0 ) { diff -= 360.0; }
308 // SG_LOG( SG_GENERAL, SG_INFO,
309 // "Runway -" << r.rwy_no << " heading = " <<
310 // r.heading + 180.0 <<
311 // " diff = " << diff );
312 if ( diff < min_diff ) {
321 // SG_LOG( SG_GENERAL, SG_INFO, "closest runway = " << r.rwy_no
322 // << " + " << found_dir );
324 // cout << "In search, rn = " << rn << endl;
325 if ( found_dir == 180 ) {
326 int irn = atoi(rn.c_str());
331 char buf[4]; // 2 chars + string terminator + 1 for safety
332 sprintf(buf, "%i", irn);
341 FGRunways::~FGRunways( void ) {
347 FGRunwaysUtil::FGRunwaysUtil() {
352 int FGRunwaysUtil::load( const string& file ) {
356 runways.erase( runways.begin(), runways.end() );
358 sg_gzifstream in( file );
359 if ( !in.is_open() ) {
360 SG_LOG( SG_GENERAL, SG_ALERT, "Cannot open file: " << file );
364 // skip first line of file
366 in.getline( tmp, 2048 );
368 // read in each line of the file
374 while ( in.get(c) && c != '\0' ) {
378 } else if ( c == 'R' ) {
381 runways.push_back(r);
391 while ( ! in.eof() ) {
397 } else if ( c == 'R' ) {
400 // cout << apt_id << " " << r.rwy_no << endl;
401 runways.push_back(r);
414 // save the data in gdbm format
415 bool FGRunwaysUtil::dump_mk4( const string& file ) {
416 // open database for writing
417 c4_Storage storage( file.c_str(), true );
419 // need to do something about error handling here!
421 // define the properties
422 c4_StringProp pID ("ID");
423 c4_StringProp pRwy ("Rwy");
424 c4_FloatProp pLon ("Longitude");
425 c4_FloatProp pLat ("Latitude");
426 c4_FloatProp pHdg ("Heading");
427 c4_FloatProp pLen ("Length");
428 c4_FloatProp pWid ("Width");
429 c4_StringProp pSurf ("SurfaceFlags");
430 c4_StringProp pEnd1 ("End1Flags");
431 c4_StringProp pEnd2 ("End2Flags");
433 // Start with an empty view of the proper structure.
435 storage.GetAs("runway[ID:S,Rwy:S,Longitude:F,Latitude:F,Heading:F,Length:F,Width:F,SurfaceFlags:S,End1Flags:S,End2Flags:S]");
439 iterator current = runways.begin();
440 const_iterator end = runways.end();
441 while ( current != end ) {
442 // add each runway record
443 pID (row) = current->id.c_str();
444 pRwy (row) = current->rwy_no.c_str();
445 pLon (row) = current->lon;
446 pLat (row) = current->lat;
447 pHdg (row) = current->heading;
448 pLen (row) = current->length;
449 pWid (row) = current->width;
450 pSurf (row) = current->surface_flags.c_str();
451 pEnd1 (row) = current->end1_flags.c_str();
452 pEnd2 (row) = current->end2_flags.c_str();
458 // commit our changes
466 // search for the specified id
468 FGRunwaysUtil::search( const string& id, FGRunway* a ) const
470 const_iterator it = runways.find( FGRunway(id) );
471 if ( it != runways.end() )
484 FGRunwaysUtil::search( const string& id ) const
487 this->search( id, &a );
493 FGRunwaysUtil::~FGRunwaysUtil( void ) {