fly/flightgear
1
0
Fork 0
Code Activity

Initial revision.

Browse source
This commit is contained in:
curt 1998-09-01 19:34:32 +00:00
parent c9d460db06
commit 96a0774d3f
4 changed files with 592 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

61
GenAirports/Makefile.am Normal file
View file

@ -0,0 +1,61 @@
#---------------------------------------------------------------------------
# Makefile
#
# Written by Curtis Olson, started January 1998.
#
# Copyright (C) 1998 Curtis L. Olson - curt@me.umn.edu
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
#
# $Id$
# (Log is kept at end of this file)
#---------------------------------------------------------------------------
bin_PROGRAMS = genapts
genapts_SOURCES = area.cxx area.hxx main.cxx
genapts_LDADD = \
$(top_builddir)/Lib/Bucket/libBucket.a \
$(top_builddir)/Lib/zlib/libz.a
INCLUDES += -I$(top_builddir) -I$(top_builddir)/Lib
#---------------------------------------------------------------------------
# $Log$
# Revision 1.1 1998/09/01 19:34:32 curt
# Initial revision.
#
# Revision 1.2 1998/07/30 23:49:18 curt
# Removed libtool support.
#
# Revision 1.1 1998/07/20 12:54:53 curt
# Whoops, need to commit Makefile.am, not Makefile.
#
# Revision 1.2 1998/04/14 02:25:59 curt
# Code reorganizations. Added a Lib/ directory for more general libraries.
#
# Revision 1.1 1998/04/08 22:54:57 curt
# Adopted Gnu automake/autoconf system.
#
# Revision 1.2 1998/01/21 02:55:46 curt
# Incorporated new make system from Bob Kuehne <rpk@sgi.com>.
#
# Revision 1.1 1998/01/15 02:45:25 curt
# Initial revision.
#

232
GenAirports/area.cxx Normal file
View file

@ -0,0 +1,232 @@
// area.c -- routines to assist with inserting "areas" into FG terrain
//
// Written by Curtis Olson, started March 1998.
//
// Copyright (C) 1998 Curtis L. Olson - curt@me.umn.edu
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//
// $Id$
// (Log is kept at end of this file)
//
#include <math.h>
#include <stdio.h>
#include <Include/fg_constants.h>
#include "area.hxx"
// calc new x, y for a rotation
double rot_x(double x, double y, double theta) {
return ( x * cos(theta) + y * sin(theta) );
}
// calc new x, y for a rotation
double rot_y(double x, double y, double theta) {
return ( -x * sin(theta) + y * cos(theta) );
}
// calc new lon/lat given starting lon/lat, and offset radial, and
// distance. NOTE: distance is specified in meters (and converted
// internally to radians)
point2d calc_lon_lat( point2d orig, point2d offset ) {
point2d result;
// printf("calc_lon_lat() offset.theta = %.2f offset.dist = %.2f\n",
// offset.theta, offset.dist);
offset.dist *= METER_TO_NM * NM_TO_RAD;
result.lat = asin( sin(orig.lat) * cos(offset.dist) +
cos(orig.lat) * sin(offset.dist) * cos(offset.theta) );
if ( cos(result.lat) < FG_EPSILON ) {
result.lon = orig.lon; // endpoint a pole
} else {
result.lon =
fmod(orig.lon - asin( sin(offset.theta) * sin(offset.dist) /
cos(result.lat) ) + FG_PI, FG_2PI) - FG_PI;
}
return(result);
}
point2d cart_to_polar_2d(point2d in) {
point2d result;
result.dist = sqrt( in.x * in.x + in.y * in.y );
result.theta = atan2(in.y, in.x);
return(result);
}
list < point2d >
batch_cart_to_polar_2d( list < point2d > in_list)
{
list < point2d > out_list;
list < point2d > :: iterator current;
list < point2d > :: iterator last;
point2d p;
current = in_list.begin();
last = in_list.end();
while ( current != last ) {
p = cart_to_polar_2d( *current );
out_list.push_back(p);
current++;
}
return out_list;
}
// given a set of 2d coordinates relative to a center point, and the
// lon, lat of that center point (specified in degrees), as well as a
// potential orientation angle, generate the corresponding lon and lat
// of the original 2d verticies.
list < point2d >
gen_area(point2d origin, double angle, list < point2d > cart_list)
{
list < point2d > rad_list;
list < point2d > result_list;
list < point2d > :: iterator current;
list < point2d > :: iterator last;
point2d origin_rad, p;
origin_rad.lon = origin.lon * DEG_TO_RAD;
origin_rad.lat = origin.lat * DEG_TO_RAD;
// convert to polar coordinates
rad_list = batch_cart_to_polar_2d(cart_list);
/*
// display points
printf("converted to polar\n");
current = rad_list.begin();
last = rad_list.end();
while ( current != last ) {
printf("(%.2f, %.2f)\n", current->theta, current->dist);
current++;
}
printf("\n");
*/
// rotate by specified angle
// printf("Rotating points by %.2f\n", angle);
current = rad_list.begin();
last = rad_list.end();
while ( current != last ) {
current->theta -= angle;
while ( current->theta > FG_2PI ) {
current->theta -= FG_2PI;
}
// printf("(%.2f, %.2f)\n", current->theta, current->dist);
current++;
}
// printf("\n");
// find actual lon,lat of coordinates
// printf("convert to lon, lat relative to %.2f %.2f\n",
// origin.lon, origin.lat);
current = rad_list.begin();
last = rad_list.end();
while ( current != last ) {
p = calc_lon_lat(origin_rad, *current);
// printf("(%.8f, %.8f)\n", p.lon, p.lat);
result_list.push_back(p);
current++;
}
// printf("\n");
return result_list;
}
// generate an area for a runway
list < point2d >
gen_runway_area( double lon, double lat, double heading,
double length, double width)
{
list < point2d > result_list;
list < point2d > tmp_list;
list < point2d > :: iterator current;
list < point2d > :: iterator last;
point2d p;
point2d origin;
double l, w;
int i;
/*
printf("runway: lon = %.2f lat = %.2f hdg = %.2f len = %.2f width = %.2f\n",
lon, lat, heading, length, width);
*/
origin.lon = lon;
origin.lat = lat;
l = length / 2.0;
w = width / 2.0;
// generate untransformed runway area vertices
p.x = l; p.y = w; tmp_list.push_back(p);
p.x = l; p.y = -w; tmp_list.push_back(p);
p.x = -l; p.y = -w; tmp_list.push_back(p);
p.x = -l; p.y = w; tmp_list.push_back(p);
/*
// display points
printf("Untransformed, unrotated runway\n");
current = tmp_list.begin();
last = tmp_list.end();
while ( current != last ) {
printf("(%.2f, %.2f)\n", current->x, current->y);
current++;
}
printf("\n");
*/
// rotate, transform, and convert points to lon, lat
result_list = gen_area(origin, heading, tmp_list);
/*
// display points
printf("Results in radians.\n");
current = result_list.begin();
last = result_list.end();
while ( current != last ) {
printf("(%.8f, %.8f)\n", current->lon, current->lat);
current++;
}
printf("\n");
*/
return result_list;
}
// $Log$
// Revision 1.1 1998/09/01 19:34:33 curt
// Initial revision.
//
// Revision 1.1 1998/07/20 12:54:05 curt
// Initial revision.
//
//

63
GenAirports/area.hxx Normal file
View file

@ -0,0 +1,63 @@
// area.h -- routines to assist with inserting "areas" into FG terrain
//
// Written by Curtis Olson, started February 1998.
//
// Copyright (C) 1998 Curtis L. Olson - curt@me.umn.edu
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//
// $Id$
// (Log is kept at end of this file)
//
#ifndef _AREA_H
#define _AREA_H
#include <list>
typedef struct {
union {
double x;
double dist;
double lon;
};
union {
double y;
double theta;
double lat;
};
} point2d;
// generate an area for a runway
list < point2d >
gen_runway_area( double lon, double lat, double heading,
double length, double width);
#endif // _AREA_H
// $Log$
// Revision 1.1 1998/09/01 19:34:33 curt
// Initial revision.
//
// Revision 1.1 1998/07/20 12:54:05 curt
// Initial revision.
//
//

236
GenAirports/main.cxx Normal file
View file

@ -0,0 +1,236 @@
// main.cxx -- main loop
//
// Written by Curtis Olson, started March 1998.
//
// Copyright (C) 1998 Curtis L. Olson - curt@me.umn.edu
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//
// $Id$
// (Log is kept at end of this file)
//
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#ifdef HAVE_STDLIB_H
#include <stdlib.h>
#endif
#include <list>
#include <stdio.h>
#include <string.h>
#include <string>
#include <Bucket/bucketutils.h>
#include <Include/fg_constants.h>
#include <Include/fg_zlib.h>
#include "area.hxx"
// process and airport + runway list
void process_airport( string last_airport, list < string > & runway_list ) {
list < point2d > rwy_list, apt_list;
list < point2d > :: iterator current;
list < point2d > :: iterator last;
string line_str;
double lon, lat;
int len, width, hdg, label_hdg, elev;
char codes[10];
char side;
printf( "(apt) %s", last_airport.c_str() );
while ( runway_list.size() ) {
line_str = runway_list.front();
runway_list.pop_front();
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,
(double)len * FEET_TO_METER,
(double)width * FEET_TO_METER );
// add rwy_list to apt_list
current = rwy_list.begin();
last = rwy_list.end();
while ( current != last ) {
apt_list.push_back(*current);
current++;
}
}
printf("Final results in degrees\n");
current = apt_list.begin();
last = apt_list.end();
while ( current != last ) {
// printf( "(%.4f, %.4f)\n",
printf( "%.5f %.5f\n",
current->lon * RAD_TO_DEG,
current->lat * RAD_TO_DEG );
current++;
}
printf("\n");
}
// reads the apt_full file and extracts and processes the individual
// 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;
*/
if ( argc != 3 ) {
printf("Usage %s <apt_file> <work dir>\n", argv[0]);
exit(0);
}
apt_path = argv[1];
gz_apt_path = apt_path + ".gz";
// 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());
}
}
while ( fggets(f, line, 250) != NULL ) {
// printf("%s", line);
if ( strlen(line) == 0 ) {
// empty, skip
} else if ( line[0] == '#' ) {
// comment, skip
} else if ( line[0] == '\t' ) {
// runway entry
runway_list.push_back(line);
} else {
// start of airport record
airport = line;
if ( last_airport.length() ) {
// process previous record
process_airport(last_airport, runway_list);
}
last_airport = airport;
}
}
if ( last_airport.length() ) {
// process previous record
process_airport(last_airport, runway_list);
}
fgclose(f);
return(1);
}
#if 0
// P13 (Globe, AZ)
// lon = -110.6642442;
// lat = 33.3528903;
// heading = 102.0 * DEG_TO_RAD;
// length = 1769;
// width = 23;
// KANE
lon = -93.2113889;
lat = 45.145;
elevation = 912 * FEET_TO_METER;
heading = 270.0 * DEG_TO_RAD;
length = 1220;
width = 23;
gen_runway_area( lon * DEG_TO_RAD, lat * DEG_TO_RAD,
heading, length, width, nodes, &count );
fgBucketFind(lon, 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);
sprintf(path, "%s/Scenery/%s", argv[1], base);
sprintf(command, "mkdir -p %s\n", path);
system(command);
sprintf(exfile, "%s/%ld.node.ex", path, index);
sprintf(file, "%s/%ld.poly", path, index);
printf( "extra node file = %s\n", exfile);
printf( "poly file = %s\n", file);
// output extra nodes
if ( (fd = fopen(exfile, "w")) == NULL ) {
printf("Cannot open file: %s\n", exfile);
exit(-1);
}
fprintf(fd, "%d 2 0 0\n", count);
for ( i = 0; i < count; i++ ) {
fprintf( fd, "%d %.2f %.2f %.2f\n", i + 1,
nodes[i].lon * RAD_TO_ARCSEC, nodes[i].lat * RAD_TO_ARCSEC,
elevation);
}
fclose(fd);
// output poly
if ( (fd = fopen(file, "w")) == NULL ) {
printf("Cannot open file: %s\n", file);
exit(-1);
}
// output empty node list
fprintf(fd, "0 2 0 0\n");
// output segments
fprintf(fd, "%d 0\n", count);
for ( i = 0; i < count - 1; i++ ) {
fprintf( fd, "%d %d %d\n", i + 1, i + 1, i + 2 );
}
fprintf( fd, "%d %d %d\n", count, count, 1 );
// output hole center
fprintf( fd, "1\n");
fprintf( fd, "1 %.2f %.2f\n", lon * 3600.0, lat * 3600);
fclose(fd);
#endif
// $Log: main.c,v
//