1 // apt_loader.cxx -- a front end loader of the apt.dat file. This loader
2 // populates the runway and basic classes.
4 // Written by Curtis Olson, started August 2000.
6 // Copyright (C) 2000 Curtis L. Olson - http://www.flightgear.org/~curt
8 // This program is free software; you can redistribute it and/or
9 // modify it under the terms of the GNU General Public License as
10 // published by the Free Software Foundation; either version 2 of the
11 // License, or (at your option) any later version.
13 // This program is distributed in the hope that it will be useful, but
14 // WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // General Public License for more details.
18 // You should have received a copy of the GNU General Public License
19 // along with this program; if not, write to the Free Software
20 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
29 #include "apt_loader.hxx"
31 #include <simgear/compiler.h>
33 #include <stdlib.h> // atof(), atoi()
34 #include <string.h> // memchr()
35 #include <ctype.h> // isspace()
37 #include <simgear/constants.h>
38 #include <simgear/debug/logstream.hxx>
39 #include <simgear/misc/sgstream.hxx>
40 #include <simgear/misc/strutils.hxx>
41 #include <simgear/structure/exception.hxx>
42 #include <simgear/bucket/newbucket.hxx>
47 #include "runways.hxx"
48 #include "pavement.hxx"
50 # include <ATCDCL/commlist.hxx>
52 #include <ATC/atcutils.hxx>
59 static FGPositioned::Type fptypeFromRobinType(int aType)
62 case 1: return FGPositioned::AIRPORT;
63 case 16: return FGPositioned::SEAPORT;
64 case 17: return FGPositioned::HELIPORT;
66 SG_LOG(SG_GENERAL, SG_ALERT, "unsupported type:" << aType);
67 throw sg_range_exception("Unsupported airport type", "fptypeFromRobinType");
85 void parseAPT(const string &aptdb_file, FGCommList *comm_list)
87 sg_gzifstream in( aptdb_file );
89 if ( !in.is_open() ) {
90 SG_LOG( SG_GENERAL, SG_ALERT, "Cannot open file: " << aptdb_file );
98 unsigned int line_id = 0;
99 unsigned int line_num = 0;
101 while ( ! in.eof() ) {
102 in.getline(tmp, 2048);
103 line = tmp; // string copy, ack
106 if ( !line.size() || isspace(tmp[0])) {
110 if (line.size() >= 3) {
111 char *p = (char *)memchr(tmp, ' ', 3);
117 if ( tmp[0] == 'I' ) {
118 // First line, indicates IBM (i.e. DOS line endings I
121 // move past this line and read and discard the next line
122 // which is the version and copyright information
123 in.getline(tmp, 2048);
124 // vector<string> vers_token = simgear::strutils::split( tmp );
125 if ( strlen(tmp) > 4 ) {
126 char *p = (char *)memchr(tmp, ' ', 4);
130 SG_LOG( SG_GENERAL, SG_INFO, "Data file version = " << tmp );
131 } else if ( line_id == 1 /* Airport */ ||
132 line_id == 16 /* Seaplane base */ ||
133 line_id == 17 /* Heliport */ ) {
134 parseAirportLine(simgear::strutils::split(line));
135 } else if ( line_id == 10 ) { // Runway v810
136 parseRunwayLine810(simgear::strutils::split(line));
137 } else if ( line_id == 100 ) { // Runway v850
138 parseRunwayLine850(simgear::strutils::split(line));
139 } else if ( line_id == 101 ) { // Water Runway v850
140 parseWaterRunwayLine850(simgear::strutils::split(line));
141 } else if ( line_id == 102 ) { // Helipad v850
142 parseHelipadLine850(simgear::strutils::split(line));
143 } else if ( line_id == 18 ) {
144 // beacon entry (ignore)
145 } else if ( line_id == 14 ) {
146 // control tower entry
147 vector<string> token(simgear::strutils::split(line));
149 double lat = atof( token[1].c_str() );
150 double lon = atof( token[2].c_str() );
151 double elev = atof( token[3].c_str() );
152 tower = SGGeod::fromDegFt(lon, lat, elev + last_apt_elev);
154 } else if ( line_id == 19 ) {
155 // windsock entry (ignore)
156 } else if ( line_id == 20 ) {
157 // Taxiway sign (ignore)
158 } else if ( line_id == 21 ) {
159 // lighting objects (ignore)
160 } else if ( line_id == 15 ) {
161 // custom startup locations (ignore)
162 } else if ( line_id == 0 ) {
164 } else if ( line_id == 50 ) {
166 parseATISLine(comm_list, simgear::strutils::split(line));
168 } else if ( line_id >= 51 && line_id <= 56 ) {
169 // other frequency entries (ignore)
170 } else if ( line_id == 110 ) {
172 parsePavementLine850(simgear::strutils::split(line, 0, 4));
173 } else if ( line_id >= 111 && line_id <= 114 ) {
175 parsePavementNodeLine850(line_id, simgear::strutils::split(line));
176 } else if ( line_id >= 115 && line_id <= 116 ) {
177 // other pavement nodes (ignore)
178 } else if ( line_id == 120 ) {
180 } else if ( line_id == 130 ) {
182 } else if ( line_id == 99 ) {
183 SG_LOG( SG_GENERAL, SG_DEBUG, "End of file reached" );
185 SG_LOG( SG_GENERAL, SG_ALERT,
186 "Unknown line(#" << line_num << ") in file: " << line );
195 vector<string> token;
196 double rwy_lat_accum;
197 double rwy_lon_accum;
198 double last_rwy_heading;
202 string last_apt_name;
203 double last_apt_elev;
205 string last_apt_info;
206 string last_apt_type;
207 string pavement_ident;
210 vector<FGRunwayPtr> runways;
211 vector<FGTaxiwayPtr> taxiways;
212 vector<FGPavementPtr> pavements;
216 if (last_apt_id.empty()) {
221 SG_LOG(SG_GENERAL, SG_ALERT, "ERROR: No runways for " << last_apt_id
226 double lat = rwy_lat_accum / (double)rwy_count;
227 double lon = rwy_lon_accum / (double)rwy_count;
230 // tower height hard coded for now...
231 const float tower_height = 50.0f;
232 // make a little off the heading for 1 runway airports...
233 float fudge_lon = fabs(sin(last_rwy_heading * SGD_DEGREES_TO_RADIANS)) * .003f;
234 float fudge_lat = .003f - fudge_lon;
235 tower = SGGeod::fromDegFt(lon + fudge_lon, lat + fudge_lat, last_apt_elev + tower_height);
238 SGGeod pos(SGGeod::fromDegFt(lon, lat, last_apt_elev));
239 FGAirport* apt = new FGAirport(last_apt_id, pos, tower, last_apt_name, false,
240 fptypeFromRobinType(atoi(last_apt_type.c_str())));
241 apt->setRunwaysAndTaxiways(runways, taxiways, pavements);
244 void parseAirportLine(const vector<string>& token)
246 const string& id(token[4]);
247 double elev = atof( token[1].c_str() );
252 last_apt_elev = elev;
257 for ( unsigned int i = 5; i < token.size() - 1; ++i ) {
258 last_apt_name += token[i];
259 last_apt_name += " ";
261 last_apt_name += token[token.size() - 1];
262 last_apt_type = token[0];
264 // clear runway list for start of next airport
270 void parseRunwayLine810(const vector<string>& token)
272 double lat = atof( token[1].c_str() );
273 double lon = atof( token[2].c_str() );
274 rwy_lat_accum += lat;
275 rwy_lon_accum += lon;
278 const string& rwy_no(token[3]);
280 double heading = atof( token[4].c_str() );
281 double length = atoi( token[5].c_str() );
282 double width = atoi( token[8].c_str() );
284 last_rwy_heading = heading;
286 int surface_code = atoi( token[10].c_str() );
287 SGGeod pos(SGGeod::fromDegFt(lon, lat, last_apt_elev));
289 if (rwy_no[0] == 'x') {
291 FGTaxiway* t = new FGTaxiway(rwy_no, pos, heading, length, width, surface_code);
292 taxiways.push_back(t);
295 string rwy_displ_threshold = token[6];
297 = simgear::strutils::split( rwy_displ_threshold, "." );
298 double displ_thresh1 = atof( displ[0].c_str() );
299 double displ_thresh2 = atof( displ[1].c_str() );
301 string rwy_stopway = token[7];
303 = simgear::strutils::split( rwy_stopway, "." );
304 double stopway1 = atof( stop[0].c_str() );
305 double stopway2 = atof( stop[1].c_str() );
307 FGRunway* rwy = new FGRunway(NULL, rwy_no, pos, heading, length,
308 width, displ_thresh1, stopway1, surface_code, false);
309 runways.push_back(rwy);
311 FGRunway* reciprocal = new FGRunway(NULL, FGRunway::reverseIdent(rwy_no),
312 pos, heading + 180.0, length, width,
313 displ_thresh2, stopway2, surface_code, true);
315 runways.push_back(reciprocal);
317 rwy->setReciprocalRunway(reciprocal);
318 reciprocal->setReciprocalRunway(rwy);
322 void parseRunwayLine850(const vector<string>& token)
324 double width = atof( token[1].c_str() ) * SG_METER_TO_FEET;
325 int surface_code = atoi( token[2].c_str() );
327 double lat_1 = atof( token[9].c_str() );
328 double lon_1 = atof( token[10].c_str() );
329 SGGeod pos_1(SGGeod::fromDegFt(lon_1, lat_1, 0.0));
330 rwy_lat_accum += lat_1;
331 rwy_lon_accum += lon_1;
334 double lat_2 = atof( token[18].c_str() );
335 double lon_2 = atof( token[19].c_str() );
336 SGGeod pos_2(SGGeod::fromDegFt(lon_2, lat_2, 0.0));
337 rwy_lat_accum += lat_2;
338 rwy_lon_accum += lon_2;
341 double length, heading_1, heading_2, dummy;
342 SGGeodesy::inverse( pos_1, pos_2, heading_1, heading_2, length );
344 SGGeodesy::direct( pos_1, heading_1, length / 2.0, pos, dummy );
345 length *= SG_METER_TO_FEET;
347 last_rwy_heading = heading_1;
349 const string& rwy_no_1(token[8]);
350 const string& rwy_no_2(token[17]);
351 if ( rwy_no_1.size() == 0 || rwy_no_2.size() == 0 )
354 double displ_thresh1 = atof( token[11].c_str() );
355 double displ_thresh2 = atof( token[20].c_str() );
357 double stopway1 = atof( token[12].c_str() );
358 double stopway2 = atof( token[21].c_str() );
360 FGRunway* rwy = new FGRunway(NULL, rwy_no_1, pos, heading_1, length,
361 width, displ_thresh1, stopway1, surface_code, false);
362 runways.push_back(rwy);
364 FGRunway* reciprocal = new FGRunway(NULL, rwy_no_2,
365 pos, heading_2, length, width,
366 displ_thresh2, stopway2, surface_code, true);
367 runways.push_back(reciprocal);
369 rwy->setReciprocalRunway(reciprocal);
370 reciprocal->setReciprocalRunway(rwy);
373 void parseWaterRunwayLine850(const vector<string>& token)
375 double width = atof( token[1].c_str() ) * SG_METER_TO_FEET;
377 double lat_1 = atof( token[4].c_str() );
378 double lon_1 = atof( token[5].c_str() );
379 SGGeod pos_1(SGGeod::fromDegFt(lon_1, lat_1, 0.0));
380 rwy_lat_accum += lat_1;
381 rwy_lon_accum += lon_1;
384 double lat_2 = atof( token[7].c_str() );
385 double lon_2 = atof( token[8].c_str() );
386 SGGeod pos_2(SGGeod::fromDegFt(lon_2, lat_2, 0.0));
387 rwy_lat_accum += lat_2;
388 rwy_lon_accum += lon_2;
391 double length, heading_1, heading_2, dummy;
392 SGGeodesy::inverse( pos_1, pos_2, heading_1, heading_2, length );
394 SGGeodesy::direct( pos_1, heading_1, length / 2.0, pos, dummy );
396 last_rwy_heading = heading_1;
398 const string& rwy_no_1(token[3]);
399 const string& rwy_no_2(token[6]);
401 FGRunway* rwy = new FGRunway(NULL, rwy_no_1, pos, heading_1, length,
402 width, 0.0, 0.0, 13, false);
403 runways.push_back(rwy);
405 FGRunway* reciprocal = new FGRunway(NULL, rwy_no_2,
406 pos, heading_2, length, width,
408 runways.push_back(reciprocal);
410 rwy->setReciprocalRunway(reciprocal);
411 reciprocal->setReciprocalRunway(rwy);
414 void parseHelipadLine850(const vector<string>& token)
416 double length = atof( token[5].c_str() ) * SG_METER_TO_FEET;
417 double width = atof( token[6].c_str() ) * SG_METER_TO_FEET;
419 double lat = atof( token[2].c_str() );
420 double lon = atof( token[3].c_str() );
421 SGGeod pos(SGGeod::fromDegFt(lon, lat, 0.0));
422 rwy_lat_accum += lat;
423 rwy_lon_accum += lon;
426 double heading = atof( token[4].c_str() );
428 last_rwy_heading = heading;
430 const string& rwy_no(token[1]);
431 int surface_code = atoi( token[7].c_str() );
433 FGRunway* rwy = new FGRunway(NULL, rwy_no, pos, heading, length,
434 width, 0.0, 0.0, surface_code, false);
435 runways.push_back(rwy);
438 void parsePavementLine850(const vector<string>& token)
440 if ( token.size() >= 5 ) {
441 pavement_ident = token[4];
442 if ( !pavement_ident.empty() && pavement_ident[pavement_ident.size()-1] == '\r' )
443 pavement_ident.erase( pavement_ident.size()-1 );
445 pavement_ident = "xx";
449 void parsePavementNodeLine850(int num, const vector<string>& token)
451 double lat = atof( token[1].c_str() );
452 double lon = atof( token[2].c_str() );
453 SGGeod pos(SGGeod::fromDegFt(lon, lat, 0.0));
456 if ( !pavement_ident.empty() ) {
457 pvt = new FGPavement( pavement_ident, pos );
458 pavements.push_back( pvt );
461 pvt = pavements.back();
463 if ( num == 112 || num == 114 ) {
464 double lat_b = atof( token[3].c_str() );
465 double lon_b = atof( token[4].c_str() );
466 SGGeod pos_b(SGGeod::fromDegFt(lon_b, lat_b, 0.0));
467 pvt->addBezierNode(pos, pos_b, num == 114);
469 pvt->addNode(pos, num == 113);
473 void parseATISLine(FGCommList *comm_list, const vector<string>& token)
475 if ( rwy_count <= 0 ) {
476 SG_LOG( SG_GENERAL, SG_ALERT,
477 "No runways; skipping AWOS for " + last_apt_id);
480 // This assumes/requires that any code-50 line (ATIS or AWOS)
481 // applies to the preceding code-1 line (airport ID and name)
482 // and that a full set of code-10 lines (runway descriptors)
483 // has come between the code-1 and code-50 lines.
484 // typical code-50 lines:
487 // This code parallels code found in "operator>>" in ATC.hxx;
488 // FIXME: unify the code.
491 a.geod = SGGeod::fromDegFt(rwy_lon_accum / (double)rwy_count,
492 rwy_lat_accum / (double)rwy_count, last_apt_elev);
493 a.range = 50; // give all ATISs small range
494 a.ident = last_apt_id;
495 a.name = last_apt_name;
496 // short int representing tens of kHz:
497 a.freq = atoi(token[1].c_str());
498 if (token[2] == "ATIS") a.type = ATIS;
499 else a.type = AWOS; // ASOS same as AWOS
501 // generate cartesian coordinates
502 a.cart = SGVec3d::fromGeod(a.geod);
503 comm_list->commlist_freq[a.freq].push_back(a);
505 SGBucket bucket(a.geod);
506 int bucknum = bucket.gen_index();
507 comm_list->commlist_bck[bucknum].push_back(a);
511 SG_LOG( SG_GENERAL, SG_ALERT,
512 "Loaded ATIS/AWOS for airport: " << a.ident
513 << " lat: " << a.geod.getLatitudeDeg()
514 << " lon: " << a.geod.getLongitudeDeg()
515 << " freq: " << a.freq
516 << " type: " << a.type );
523 // Load the airport data base from the specified aptdb file. The
524 // metar file is used to mark the airports as having metar available
526 bool fgAirportDBLoad( const string &aptdb_file,
527 FGCommList *comm_list, const std::string &metar_file )
531 ld.parseAPT(aptdb_file, comm_list);
533 // Load the metar.dat file and update apt db with stations that
537 sg_gzifstream metar_in( metar_file );
538 if ( !metar_in.is_open() ) {
539 SG_LOG( SG_GENERAL, SG_ALERT, "Cannot open file: " << metar_file );
545 if ( ident == "#" || ident == "//" ) {
548 FGAirport* apt = FGAirport::findByIdent(ident);
555 SG_LOG(SG_GENERAL, SG_INFO, "[FINISHED LOADING]");