#include <config.h>
#endif
+#include <Include/compiler.h>
+
#ifdef HAVE_STDLIB_H
#include <stdlib.h>
#endif
#include <list>
#include <stdio.h>
#include <string.h>
-#include <string>
+#include STL_STRING
-#include <Bucket/bucketutils.h>
+#include <Bucket/newbucket.hxx>
+#include <Debug/logstream.hxx>
#include <Include/fg_constants.h>
-#include <Include/fg_zlib.h>
+#include <Misc/fgstream.hxx>
+#include <Polygon/index.hxx>
#include "area.hxx"
#include "convex_hull.hxx"
+// write out airport data
+void write_airport( int p_index, list_container hull_list, FGBucket b,
+ const string& root, const bool cut_and_keep ) {
+ char tile_name[256], poly_index[256];
+
+ long int b_index = b.gen_index();
+ string base = b.gen_base_path();
+ string path = root + "/Scenery/" + base;
+ string command = "mkdir -p " + path;
+ system( command.c_str() );
+
+ sprintf(tile_name, "%ld", b_index);
+ string aptfile = path + "/" + tile_name;
+
+ sprintf( poly_index, "%d", p_index );
+ aptfile += ".";
+ aptfile += poly_index;
+ cout << "apt file = " << aptfile << endl;
+
+ FILE *fd;
+ if ( (fd = fopen(aptfile.c_str(), "a")) == NULL ) {
+ cout << "Cannot open file: " << aptfile << endl;
+ exit(-1);
+ }
+
+ // polygon type
+ if ( cut_and_keep ) {
+ fprintf( fd, "AirportKeep\n" );
+ } else {
+ fprintf( fd, "AirportIgnore\n" );
+ }
+
+ // number of contours
+ fprintf( fd, "1\n" );
+
+ // size of first contour
+ fprintf( fd, "%d\n", hull_list.size() );
+
+ // write contour (polygon) points
+ list_iterator current = hull_list.begin();
+ list_iterator last = hull_list.end();
+ for ( ; current != last ; ++current ) {
+ fprintf( fd, "%.7f %.7f\n", (*current).lon, (*current).lat );
+ }
+
+ fclose(fd);
+}
+
+
// process and airport + runway list
-void process_airport( string last_airport, list < string > & runway_list,
+void process_airport( string airport, list < string > & runway_list,
const string& root ) {
list_container rwy_list, apt_list, hull_list;
list_iterator current, last;
- string line_str;
- double lon, lat;
- int len, width, hdg, label_hdg, elev;
- char codes[10];
- char side;
- point2d average;
- double sum_x, sum_y;
+ // parse main airport information
+ int elev;
+
+ cout << airport << endl;
+ string apt_type = airport.substr(0, 1);
+ string apt_code = airport.substr(2, 4);
+ string apt_lat = airport.substr(7, 10);
+ string apt_lon = airport.substr(18, 11);
+ string apt_elev = airport.substr(30, 5);
+ sscanf( apt_elev.c_str(), "%d", &elev );
+ string apt_use = airport.substr(36, 1);
+ string apt_twr = airport.substr(37, 1);
+ string apt_bldg = airport.substr(38, 1);
+ string apt_name = airport.substr(40);
- FILE *fd;
- fgBUCKET b;
- long int index;
- char base[256], tmp[256];
- string path, command, exfile, file, aptfile;
- int i, count;
+ /*
+ cout << " type = " << apt_type << endl;
+ cout << " code = " << apt_code << endl;
+ cout << " lat = " << apt_lat << endl;
+ cout << " lon = " << apt_lon << endl;
+ cout << " elev = " << apt_elev << " " << elev << endl;
+ cout << " use = " << apt_use << endl;
+ cout << " twr = " << apt_twr << endl;
+ cout << " bldg = " << apt_bldg << endl;
+ cout << " name = " << apt_name << endl;
+ */
- printf( "(apt) %s", last_airport.c_str() );
+ // parse runways and generate the vertex list
+ string rwy_str;
+ double lon, lat, hdg;
+ int len, width;
list < string >::iterator last_runway = runway_list.end();
for ( list < string >::iterator current_runway = runway_list.begin();
current_runway != last_runway ; ++current_runway ) {
- line_str = (*current_runway);
- printf( "%s", line_str.c_str() );
-
- sscanf( line_str.c_str(), "%lf %lf %d %d %d %s %d %c %d\n",
- &lon, &lat, &len, &width, &hdg, codes, &label_hdg,
- &side, &elev );
-
- rwy_list = gen_runway_area( lon, lat, (double)hdg * DEG_TO_RAD,
+ rwy_str = (*current_runway);
+
+ cout << rwy_str << endl;
+ string rwy_no = rwy_str.substr(2, 4);
+ string rwy_lat = rwy_str.substr(6, 10);
+ sscanf( rwy_lat.c_str(), "%lf", &lat);
+ string rwy_lon = rwy_str.substr(17, 11);
+ sscanf( rwy_lon.c_str(), "%lf", &lon);
+ string rwy_hdg = rwy_str.substr(29, 7);
+ sscanf( rwy_hdg.c_str(), "%lf", &hdg);
+ string rwy_len = rwy_str.substr(36, 7);
+ sscanf( rwy_len.c_str(), "%d", &len);
+ string rwy_width = rwy_str.substr(43, 4);
+ sscanf( rwy_width.c_str(), "%d", &width);
+ string rwy_sfc = rwy_str.substr(47, 4);
+ string rwy_end1 = rwy_str.substr(52, 6);
+ string rwy_end2 = rwy_str.substr(59, 6);
+
+ /*
+ cout << " no = " << rwy_no << endl;
+ cout << " lat = " << rwy_lat << " " << lat << endl;
+ cout << " lon = " << rwy_lon << " " << lon << endl;
+ cout << " hdg = " << rwy_hdg << " " << hdg << endl;
+ cout << " len = " << rwy_len << " " << len << endl;
+ cout << " width = " << rwy_width << " " << width << endl;
+ cout << " sfc = " << rwy_sfc << endl;
+ cout << " end1 = " << rwy_end1 << endl;
+ cout << " end2 = " << rwy_end2 << endl;
+ */
+
+ rwy_list = gen_runway_area( lon, lat, hdg * DEG_TO_RAD,
(double)len * FEET_TO_METER,
(double)width * FEET_TO_METER );
}
}
- printf("Runway points in degrees\n");
- current = apt_list.begin();
- last = apt_list.end();
- for ( ; current != last; ++current ) {
- // printf( "(%.4f, %.4f)\n",
- printf( "%.5f %.5f\n", current->lon, current->lat );
+ if ( apt_list.size() == 0 ) {
+ cout << "no runway points generated" << endl;
+ return;
}
- printf("\n");
+
+ // printf("Runway points in degrees\n");
+ // current = apt_list.begin();
+ // last = apt_list.end();
+ // for ( ; current != last; ++current ) {
+ // printf( "%.5f %.5f\n", current->lon, current->lat );
+ // }
+ // printf("\n");
// generate convex hull
hull_list = convex_hull(apt_list);
- // find average center point of convex hull
- count = hull_list.size();
+ // get next polygon index
+ int index = poly_index_next();
+ // find average center, min, and max point of convex hull
+ point2d average, min, max;
+ double sum_x, sum_y;
+ int count = hull_list.size();
current = hull_list.begin();
last = hull_list.end();
sum_x = sum_y = 0.0;
+ min.x = min.y = 200.0;
+ max.x = max.y = -200.0;
for ( ; current != last; ++current ) {
// printf("return = %.6f %.6f\n", (*current).x, (*current).y);
sum_x += (*current).x;
sum_y += (*current).y;
- }
+ if ( (*current).x < min.x ) { min.x = (*current).x; }
+ if ( (*current).y < min.y ) { min.y = (*current).y; }
+ if ( (*current).x > max.x ) { max.x = (*current).x; }
+ if ( (*current).y > max.y ) { max.y = (*current).y; }
+ }
average.x = sum_x / count;
average.y = sum_y / count;
- // find bucket based on average center point of hull list.
- // Eventually we'll need to handle cases where the area crosses
- // bucket boundaries.
- fgBucketFind( average.lon, average.lat, &b);
- printf( "Bucket = lon,lat = %d,%d x,y index = %d,%d\n",
- b.lon, b.lat, b.x, b.y);
-
- index = fgBucketGenIndex(&b);
- fgBucketGenBasePath(&b, base);
- path = root + "/Scenery/" + base;
- command = "mkdir -p " + path;
- system( command.c_str() );
-
- sprintf(tmp, "%ld", index);
- exfile = path + "/" + tmp + ".node.ex";
- file = path + "/" + tmp + ".poly";
- aptfile = path + "/" + tmp + ".apt";
- printf( "extra node file = %s\n", exfile.c_str() );
- printf( "poly file = %s\n", file.c_str() );
- printf( "apt file = %s\n", aptfile.c_str() );
-
- // output exclude nodes
- printf("Output exclude nodes\n");
- if ( (fd = fopen(exfile.c_str(), "w")) == NULL ) {
- printf("Cannot open file: %s\n", exfile.c_str());
- exit(-1);
- }
-
- fprintf( fd, "%d 2 0 0\n", count );
-
- current = hull_list.begin();
- last = hull_list.end();
- i = 1;
- for ( ; current != last ; ++current ) {
- // printf( "(%.4f, %.4f)\n",
- fprintf( fd, "%d %.2f %.2f %.2f\n", i,
- (*current).lon * 3600.0, (*current).lat * 3600.0,
- (double)elev * FEET_TO_METER );
- ++i;
- }
- fclose(fd);
-
- // output poly
- printf("Output poly\n");
- if ( (fd = fopen(file.c_str(), "w")) == NULL ) {
- printf("Cannot open file: %s\n", file.c_str());
- exit(-1);
- }
-
- // output empty node list
- fprintf(fd, "0 2 0 0\n");
-
- // output segments
- fprintf( fd, "%d 0\n", count );
- for ( i = 1; i < count; i++ ) {
- fprintf( fd, "%d %d %d\n", i, i, i + 1 );
- }
- fprintf( fd, "%d %d %d\n", count, count, 1 );
-
- // output hole center
- fprintf( fd, "1\n");
- fprintf( fd, "1 %.2f %.2f\n", average.x * 3600.0, average.y * 3600);
-
- fclose(fd);
-
- // output "apt" file
- printf("Output airport\n");
- if ( (fd = fopen(aptfile.c_str(), "w")) == NULL ) {
- printf("Cannot open file: %s\n", aptfile.c_str());
- exit(-1);
- }
-
- // write main airport identifier
- fprintf(fd, "a %s", last_airport.c_str() );
-
- // write perimeter polygon
- current = hull_list.begin();
- last = hull_list.end();
- for ( ; current != last ; ++current ) {
- fprintf( fd, "p %.7f %.7f %.2f\n", (*current).lon, (*current).lat,
- (double)elev * FEET_TO_METER );
- }
+ // find buckets for center, min, and max points of convex hull.
+ // note to self: self, you should think about checking for runways
+ // that span the data line
+ FGBucket b(average.lon, average.lat);
+ FGBucket b_min(min.x, min.y);
+ FGBucket b_max(max.x, max.y);
+ cout << "Bucket center = " << b << endl;
+ cout << "Bucket min = " << b_min << endl;
+ cout << "Bucket max = " << b_max << endl;
+
+ if ( b_min == b_max ) {
+ write_airport( index, hull_list, b, root, true );
+ } else {
+ FGBucket b_cur;
+ int dx, dy, i, j;
+
+ fgBucketDiff(b_min, b_max, &dx, &dy);
+ cout << "airport spans tile boundaries" << endl;
+ cout << " dx = " << dx << " dy = " << dy << endl;
+
+ if ( (dx > 1) || (dy > 1) ) {
+ cout << "somethings really wrong!!!!" << endl;
+ exit(-1);
+ }
- // write runway info
- for ( list < string >::iterator current_runway = runway_list.begin();
- current_runway != last_runway ; ++current_runway ) {
- line_str = (*current_runway);
- line_str = line_str.substr(1, line_str.size());
- fprintf(fd, "r %s", line_str.c_str() );
+ for ( j = 0; j <= dy; j++ ) {
+ for ( i = 0; i <= dx; i++ ) {
+ b_cur = fgBucketOffset(min.x, min.y, i, j);
+ if ( b_cur == b ) {
+ write_airport( index, hull_list, b_cur, root, true );
+ } else {
+ write_airport( index, hull_list, b_cur, root, false );
+ }
+ }
+ }
+ // string answer; cin >> answer;
}
-
- fclose(fd);
}
// airport records
int main( int argc, char **argv ) {
list < string > runway_list;
- string apt_path, gz_apt_path;
string airport, last_airport;
- fgFile f;
- char line[256];
- /*
- fgBUCKET b;
- point2d nodes[4];
- char base[256], path[256], command[256], file[256], exfile[256];
- double lon, lat, elevation, heading;
- double length, width;
- long int index;
- */
+ string line;
+ char tmp[256];
+
+ fglog().setLogLevels( FG_ALL, FG_DEBUG );
if ( argc != 3 ) {
- printf("Usage %s <apt_file> <work dir>\n", argv[0]);
- exit(0);
+ FG_LOG( FG_GENERAL, FG_ALERT,
+ "Usage " << argv[0] << " <apt_file> <work_dir>" );
+ exit(-1);
}
- apt_path = argv[1];
- gz_apt_path = apt_path + ".gz";
+ // make work directory
+ string work_dir = argv[2];
+ string command = "mkdir -p " + work_dir;
+ system( command.c_str() );
- // first try "path.gz"
- if ( (f = fgopen(gz_apt_path.c_str(), "rb")) == NULL ) {
- // next try "path"
- if ( (f = fgopen(apt_path.c_str(), "rb")) == NULL ) {
- printf( "Cannot open file: %s\n", apt_path.c_str());
- }
+ // initialize persistant polygon counter
+ string counter_file = work_dir + "/polygon.counter";
+ poly_index_init( counter_file );
+
+ fg_gzifstream in( argv[1] );
+ if ( !in ) {
+ FG_LOG( FG_GENERAL, FG_ALERT, "Cannot open file: " << argv[1] );
+ exit(-1);
}
- while ( fggets(f, line, 250) != NULL ) {
- // printf("%s", line);
- if ( strlen(line) == 0 ) {
+ // throw away the first 3 lines
+ in.getline(tmp, 256);
+ in.getline(tmp, 256);
+ in.getline(tmp, 256);
+
+ last_airport = "";
+
+ while ( ! in.eof() ) {
+ in.getline(tmp, 256);
+ line = tmp;
+ // cout << line << endl;
+
+ if ( line.length() == 0 ) {
// empty, skip
} else if ( line[0] == '#' ) {
// comment, skip
- } else if ( line[0] == '\t' ) {
- // runway entry
- runway_list.push_back(line);
- } else {
+ } else if ( (line[0] == 'A') || (line[0] == 'S') ) {
// start of airport record
airport = line;
runway_list.erase(runway_list.begin(), runway_list.end());
last_airport = airport;
+ } else if ( line[0] == 'R' ) {
+ // runway entry
+ runway_list.push_back(line);
+ } else if ( line == "99" ) {
+ // end of file
+ break;
+ } else {
+ FG_LOG( FG_GENERAL, FG_ALERT,
+ "Unknown line in file" << endl << line );
+ exit(-1);
}
}
process_airport(last_airport, runway_list, argv[2]);
}
- fgclose(f);
-
- return(1);
+ return 0;
}
// $Log$
+// Revision 1.8 1999/03/01 15:35:26 curt
+// Fixed bug in output format generated.
+//
+// Revision 1.7 1999/02/25 21:32:49 curt
+// Modified to adhere to new polygon naming convention, and also to read the
+// new Robin Peel aiport format.
+//
+// Revision 1.6 1999/02/11 01:10:51 curt
+// Start of scenery revamp project.
+//
+// Revision 1.5 1998/09/17 18:40:43 curt
+// Debug message tweaks.
+//
// Revision 1.4 1998/09/09 20:59:56 curt
// Loop construct tweaks for STL usage.
// Output airport file to be used to generate airport scenery on the fly
projects / flightgear.git / blobdiff