// FGGround - a class to provide ground control at larger airports. // // Written by David Luff, started March 2002. // // Copyright (C) 2002 David C. Luff - david.luff@nottingham.ac.uk // // 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. #include #include #include #include #include #include
#include #include STL_FSTREAM #include "ground.hxx" #include "ATCmgr.hxx" #include "ATCutils.hxx" #include "ATCdisplay.hxx" #include "AILocalTraffic.hxx" SG_USING_STD(ifstream); SG_USING_STD(cout); node::node() { } node::~node() { for(unsigned int i=0; i < arcs.size(); ++i) { delete arcs[i]; } } // Make sure that a_path.cost += distance is safe from the moment it's created. a_path::a_path() { cost = 0; } FGGround::FGGround() { ATCmgr = globals->get_ATC_mgr(); display = false; networkLoadOK = false; ground_traffic.erase(ground_traffic.begin(), ground_traffic.end()); ground_traffic_itr = ground_traffic.begin(); // Init the property nodes - TODO - need to make sure we're getting surface winds. wind_from_hdg = fgGetNode("/environment/wind-from-heading-deg", true); wind_speed_knots = fgGetNode("/environment/wind-speed-kts", true); } FGGround::FGGround(string id) { ATCmgr = globals->get_ATC_mgr(); display = false; networkLoadOK = false; ground_traffic.erase(ground_traffic.begin(), ground_traffic.end()); ground_traffic_itr = ground_traffic.begin(); ident = id; // Init the property nodes - TODO - need to make sure we're getting surface winds. wind_from_hdg = fgGetNode("/environment/wind-from-heading-deg", true); wind_speed_knots = fgGetNode("/environment/wind-speed-kts", true); } FGGround::~FGGround() { } void FGGround::ParseRwyExits(node* np, char* es) { char* token; char estr[20]; strcpy(estr, es); const char delimiters[] = "-"; token = strtok(estr, delimiters); while(token != NULL) { int i = atoi(token); //cout << "token = " << token << endl; //cout << "rwy number = " << i << endl; //runways[(atoi(token))].exits.push_back(np); runways[i].exits.push_back(np); //cout << "token = " << token << '\n'; token = strtok(NULL, delimiters); } } // Load the ground logical network of the current instances airport // Return true if successfull. // TODO - currently the file is assumed to reside in the base/ATC directory. // This might change to something more thought out in the future. // NOTE - currently it is assumed that all nodes are loaded before any arcs. // It won't work ATM if this doesn't hold true. bool FGGround::LoadNetwork() { node* np; arc* ap; Gate* gp; int gateCount = 0; // This is used to allocate gateID's from zero upwards // This may well change in the future - probably to reading in the real-world // gate numbers from file. ifstream fin; SGPath path = globals->get_fg_root(); //string taxiPath = "ATC/" + ident + ".taxi"; string taxiPath = "ATC/KEMT.taxi"; // FIXME - HARDWIRED FOR TESTING path.append(taxiPath); SG_LOG(SG_ATC, SG_INFO, "Trying to read taxiway data for " << ident << "..."); //cout << "Trying to read taxiway data for " << ident << "..." << endl; fin.open(path.c_str(), ios::in); if(!fin) { SG_LOG(SG_ATC, SG_ALERT, "Unable to open taxiway data input file " << path.c_str()); //cout << "Unable to open taxiway data input file " << path.c_str() << endl; return(false); } char ch; char buf[30]; while(!fin.eof()) { fin >> buf; // Node, arc, or [End]? //cout << "Read in ground network element type = " << buf << endl; if(!strcmp(buf, "[End]")) { // TODO - maybe make this more robust to spelling errors by just looking for '[' SG_LOG(SG_ATC, SG_INFO, "Done reading " << path.c_str() << endl); break; } else if(!strcmp(buf, "N")) { // Node np = new node; np->struct_type = NODE; fin >> buf; np->nodeID = atoi(buf); fin >> buf; np->pos.setlon(atof(buf)); fin >> buf; np->pos.setlat(atof(buf)); fin >> buf; np->pos.setelev(atof(buf)); fin >> buf; // node type if(!strcmp(buf, "J")) { np->type = JUNCTION; } else if(!strcmp(buf, "T")) { np->type = TJUNCTION; } else if(!strcmp(buf, "H")) { np->type = HOLD; } else { SG_LOG(SG_ATC, SG_ALERT, "**** ERROR ***** Unknown node type in taxi network...\n"); delete np; return(false); } fin >> buf; // rwy exit information - gets parsed later - FRAGILE - will break if buf is reused. // Now the name fin >> ch; // strip the leading " off np->name = ""; while(1) { fin.unsetf(ios::skipws); fin >> ch; if((ch == '"') || (ch == 0x0A)) { break; } // we shouldn't need the 0x0A but it makes a nice safely in case someone leaves off the " np->name += ch; } fin.setf(ios::skipws); network.push_back(np); // FIXME - fragile - replies on buf not getting modified from exits read to here // see if we also need to push it onto the runway exit list //cout << "strlen(buf) = " << strlen(buf) << endl; if(strlen(buf) > 2) { //cout << "Calling ParseRwyExits for " << buf << endl; ParseRwyExits(np, buf); } } else if(!strcmp(buf, "A")) { ap = new arc; ap->struct_type = ARC; fin >> buf; ap->n1 = atoi(buf); fin >> buf; ap->n2 = atoi(buf); fin >> buf; if(!strcmp(buf, "R")) { ap->type = RUNWAY; } else if(!strcmp(buf, "T")) { ap->type = TAXIWAY; } else { SG_LOG(SG_ATC, SG_ALERT, "**** ERROR ***** Unknown arc type in taxi network...\n"); delete ap; return(false); } // directed? fin >> buf; if(!strcmp(buf, "Y")) { ap->directed = true; } else if(!strcmp(buf, "N")) { ap->directed = false; } else { SG_LOG(SG_ATC, SG_ALERT, "**** ERROR ***** Unknown arc directed value in taxi network - should be Y/N !!!\n"); delete ap; return(false); } // Now the name ap->name = ""; while(1) { fin.unsetf(ios::skipws); fin >> ch; ap->name += ch; if((ch == '"') || (ch == 0x0A)) { break; } // we shouldn't need the 0x0A but it makes a nice safely in case someone leaves off the " } fin.setf(ios::skipws); ap->distance = (int)dclGetHorizontalSeparation(network[ap->n1]->pos, network[ap->n2]->pos); //cout << "Distance = " << ap->distance << '\n'; network[ap->n1]->arcs.push_back(ap); network[ap->n2]->arcs.push_back(ap); } else if(!strcmp(buf, "G")) { gp = new Gate; gp->struct_type = NODE; gp->type = GATE; fin >> buf; gp->nodeID = atoi(buf); fin >> buf; gp->pos.setlon(atof(buf)); fin >> buf; gp->pos.setlat(atof(buf)); fin >> buf; gp->pos.setelev(atof(buf)); fin >> buf; // gate type - ignore this for now fin >> buf; // gate heading gp->heading = atoi(buf); // Now the name gp->name = ""; while(1) { fin.unsetf(ios::skipws); fin >> ch; gp->name += ch; if((ch == '"') || (ch == 0x0A)) { break; } // we shouldn't need the 0x0A but it makes a nice safely in case someone leaves off the " } fin.setf(ios::skipws); gp->id = gateCount; // Warning - this will likely change in the future. gp->used = false; network.push_back(gp); gates[gateCount] = gp; gateCount++; } else { // Something has gone seriously pear-shaped SG_LOG(SG_ATC, SG_ALERT, "********* ERROR - unknown ground network element type... aborting read of " << path.c_str() << '\n'); return(false); } fin >> skipeol; } return(true); } void FGGround::Init() { display = false; untowered = false; // Figure out which is the active runway - TODO - it would be better to have ground call tower // for runway operation details, but at the moment we can't guarantee that tower control at a given airport // will be initialised before ground so we can't do that. DoRwyDetails(); //cout << "In FGGround::Init, active rwy is " << activeRwy << '\n'; ortho.Init(rwy.threshold_pos, rwy.hdg); networkLoadOK = LoadNetwork(); } void FGGround::Update(double dt) { // Each time step, what do we need to do? // We need to go through the list of outstanding requests and acknowedgements // and process at least one of them. // We need to go through the list of planes under our control and check if // any need to be addressed. // We need to check for planes not under our control coming within our // control area and address if necessary. // Lets take the example of a plane which has just contacted ground // following landing - presumably requesting where to go? // First we need to establish the position of the plane within the logical network. // Next we need to decide where its going. if(ground_traffic.size()) { if(ground_traffic_itr == ground_traffic.end()) { ground_traffic_itr = ground_traffic.begin(); } //Process(*ground_traffic_itr); GroundRec* g = *ground_traffic_itr; if(g->taxiRequestOutstanding) { double responseTime = 10.0; // seconds - this should get more sophisticated at some point if(g->clearanceCounter > responseTime) { // DO CLEARANCE // TODO - move the mechanics of making up the transmission out of the main Update(...) routine. string trns = ""; trns += g->plane.callsign; trns += " taxi holding point runway "; // TODO - add the holding point name // eg " taxi holding point G2 runway " trns += ConvertRwyNumToSpokenString(activeRwy); if(display) { globals->get_ATC_display()->RegisterSingleMessage(trns, 0); } g->planePtr->RegisterTransmission(1); // cleared to taxi g->clearanceCounter = 0.0; g->taxiRequestOutstanding = false; } else { g->clearanceCounter += (dt * ground_traffic.size()); } } else if(((FGAILocalTraffic*)(g->planePtr))->AtHoldShort()) { // That's a hack - eventually we should monitor actual position // HACK ALERT - the automatic cast to AILocalTraffic has to go once we have other sorts working!!!!! FIXME TODO // NOTE - we don't need to do the contact tower bit unless we have separate tower and ground string trns = g->plane.callsign; trns += " contact Tower "; double f = globals->get_ATC_mgr()->GetFrequency(ident, TOWER); char buf[10]; sprintf(buf, "%f", f); trns += buf; if(display) { globals->get_ATC_display()->RegisterSingleMessage(trns, 0); } g->planePtr->RegisterTransmission(2); // contact tower delete *ground_traffic_itr; ground_traffic.erase(ground_traffic_itr); ground_traffic_itr = ground_traffic.begin(); } ++ground_traffic_itr; } } // Figure out which runways are active. // For now we'll just be simple and do one active runway - eventually this will get much more complex // Copied from FGTower - TODO - it would be better to implement this just once, and have ground call tower // for runway operation details, but at the moment we can't guarantee that tower control at a given airport // will be initialised before ground so we can't do that. void FGGround::DoRwyDetails() { //cout << "GetRwyDetails called" << endl; // Based on the airport-id and wind get the active runway SGPath path( globals->get_fg_root() ); path.append( "Airports" ); path.append( "runways.mk4" ); FGRunways r_ways( path.c_str() ); //wind double hdg = wind_from_hdg->getDoubleValue(); double speed = wind_speed_knots->getDoubleValue(); hdg = (speed == 0.0 ? 270.0 : hdg); //cout << "Heading = " << hdg << '\n'; FGRunway runway; bool rwyGood = r_ways.search(ident, int(hdg), &runway); if(rwyGood) { activeRwy = runway.rwy_no; rwy.rwyID = runway.rwy_no; SG_LOG(SG_ATC, SG_INFO, "In FGGround, active runway for airport " << ident << " is " << activeRwy); // Get the threshold position double other_way = runway.heading - 180.0; while(other_way <= 0.0) { other_way += 360.0; } // move to the +l end/center of the runway //cout << "Runway center is at " << runway.lon << ", " << runway.lat << '\n'; Point3D origin = Point3D(runway.lon, runway.lat, aptElev); Point3D ref = origin; double tshlon, tshlat, tshr; double tolon, tolat, tor; rwy.length = runway.length * SG_FEET_TO_METER; geo_direct_wgs_84 ( aptElev, ref.lat(), ref.lon(), other_way, rwy.length / 2.0 - 25.0, &tshlat, &tshlon, &tshr ); geo_direct_wgs_84 ( aptElev, ref.lat(), ref.lon(), runway.heading, rwy.length / 2.0 - 25.0, &tolat, &tolon, &tor ); // Note - 25 meters in from the runway end is a bit of a hack to put the plane ahead of the user. // now copy what we need out of runway into rwy rwy.threshold_pos = Point3D(tshlon, tshlat, aptElev); Point3D takeoff_end = Point3D(tolon, tolat, aptElev); //cout << "Threshold position = " << tshlon << ", " << tshlat << ", " << aptElev << '\n'; //cout << "Takeoff position = " << tolon << ", " << tolat << ", " << aptElev << '\n'; rwy.hdg = runway.heading; // Set the projection for the local area based on this active runway ortho.Init(rwy.threshold_pos, rwy.hdg); rwy.end1ortho = ortho.ConvertToLocal(rwy.threshold_pos); // should come out as zero rwy.end2ortho = ortho.ConvertToLocal(takeoff_end); } else { SG_LOG(SG_ATC, SG_ALERT, "Help - can't get good runway in FGTower!!"); activeRwy = "NN"; } } // Return a random gate ID of an unused gate. // Two error values may be returned and must be checked for by the calling function: // -2 signifies that no gates exist at this airport. // -1 signifies that all gates are currently full. int FGGround::GetRandomGateID() { // Check that this airport actually has some gates!! if(!gates.size()) { return(-2); } gate_vec_type gateVec; int num = 0; int thenum; int ID; gatesItr = gates.begin(); while(gatesItr != gates.end()) { if((gatesItr->second)->used == false) { gateVec.push_back(gatesItr->second); num++; } ++gatesItr; } // Check that there are some unused gates! if(!gateVec.size()) { return(-1); } // Randomly select one from the list sg_srandom_time(); thenum = (int)(sg_random() * gateVec.size()); ID = gateVec[thenum]->id; return(ID); } // Return a pointer to an unused gate node Gate* FGGround::GetGateNode() { int id = GetRandomGateID(); if(id < 0) { return(NULL); } else { return(gates[id]); } } // WARNING - This is hardwired to my prototype logical network format // and will almost certainly change when Bernie's stuff comes on-line. node* FGGround::GetThresholdNode(string rwyID) { // For now go through all the nodes and parse their names // Maybe in the future we'll map threshold nodes by ID //cout << "Size of network is " << network.size() << '\n'; for(unsigned int i=0; ipath.push_back(A); pathPtr->cost = 0; pathMap[A->nodeID] = pathPtr; bool solution_found = false; // Flag to indicate that at least one candidate path has been found int solution_cost = -1; // Cost of current best cost solution. -1 indicates no solution found yet. a_path solution_path; node* nPtr; // nPtr is used to point to the node we are currently working with while(nodesLeft.size()) { //cout << "\n*****nodesLeft*****\n"; //for(unsigned int i=0; inodeID << '\n'; //} //cout << "*******************\n\n"; nPtr = *nodesLeft.begin(); // Thought - definate optimization possibilities here in the choice of which nodes we process first. nodesLeft.erase(nodesLeft.begin()); //cout << "Walking out from node " << nPtr->nodeID << '\n'; for(unsigned int i=0; iarcs.size(); ++i) { //cout << "ARC TO " << ((nPtr->arcs[i]->n1 == nPtr->nodeID) ? nPtr->arcs[i]->n2 : nPtr->arcs[i]->n1) << '\n'; } if((solution_found) && (solution_cost <= pathMap[nPtr->nodeID]->cost)) { // Do nothing - we've already found a solution and this partial path is already more expensive } else { // This path could still be better than the current solution - check it out for(unsigned int i=0; i<(nPtr->arcs.size()); i++) { // Map the new path against the end node, ie. *not* the one we just started with. unsigned int end_nodeID = ((nPtr->arcs[i]->n1 == nPtr->nodeID) ? nPtr->arcs[i]->n2 : nPtr->arcs[i]->n1); //cout << "end_nodeID = " << end_nodeID << '\n'; //cout << "pathMap size is " << pathMap.size() << '\n'; if(end_nodeID == nPtr->nodeID) { //cout << "Circular arc!\n"; // Then its a circular arc - don't bother!! //nPtr->arcs.erase(nPtr->arcs.begin() + i); } else { // see if the end node is already in the map or not if(pathMap.find(end_nodeID) == pathMap.end()) { //cout << "Not in the map" << endl;; // Not in the map - easy! pathPtr = new a_path; pathsCreated++; *pathPtr = *pathMap[nPtr->nodeID]; // *copy* the path pathPtr->path.push_back(nPtr->arcs[i]); pathPtr->path.push_back(network[end_nodeID]); pathPtr->cost += nPtr->arcs[i]->distance; pathMap[end_nodeID] = pathPtr; nodesLeft.push_back(network[end_nodeID]); // By definition this can't be in the list already, or // it would also have been in the map and hence OR'd with this one. if(end_nodeID == B->nodeID) { //cout << "Solution found!!!" << endl; // Since this node wasn't in the map this is by definition the first solution solution_cost = pathPtr->cost; solution_path = *pathPtr; solution_found = true; } } else { //cout << "Already in the map" << endl; // In the map - not so easy - need to get rid of an arc from the higher cost one. //cout << "Current cost of node " << end_nodeID << " is " << pathMap[end_nodeID]->cost << endl; int newCost = pathMap[nPtr->nodeID]->cost + nPtr->arcs[i]->distance; //cout << "New cost is of node " << nPtr->nodeID << " is " << newCost << endl; if(newCost >= pathMap[end_nodeID]->cost) { // No need to do anything. //cout << "Not doing anything!" << endl; } else { delete pathMap[end_nodeID]; pathsCreated--; pathPtr = new a_path; pathsCreated++; *pathPtr = *pathMap[nPtr->nodeID]; // *copy* the path pathPtr->path.push_back(nPtr->arcs[i]); pathPtr->path.push_back(network[end_nodeID]); pathPtr->cost += nPtr->arcs[i]->distance; pathMap[end_nodeID] = pathPtr; // We need to add this node to the list-to-do again to force a recalculation // onwards from this node with the new lower cost to node cost. nodesLeft.push_back(network[end_nodeID]); if(end_nodeID == B->nodeID) { //cout << "Solution found!!!" << endl; // Need to check if there is a previous better solution if((solution_cost < 0) || (pathPtr->cost < solution_cost)) { solution_cost = pathPtr->cost; solution_path = *pathPtr; solution_found = true; } } } } } } } } // delete all the paths before returning shortest_path_map_iterator spItr = pathMap.begin(); while(spItr != pathMap.end()) { if(spItr->second != NULL) { delete spItr->second; --pathsCreated; } ++spItr; } //cout << "pathsCreated = " << pathsCreated << '\n'; if(pathsCreated > 0) { SG_LOG(SG_ATC, SG_ALERT, "WARNING - Possible memory leak in FGGround::GetShortestPath\n\ Please report to flightgear-devel@flightgear.org\n"); } //cout << (solution_found ? "Result: solution found\n" : "Result: no solution found\n"); return(solution_path.path); // TODO - we really ought to have a fallback position incase a solution isn't found. } // Randomly or otherwise populate some of the gates with parked planes // (might eventually be done by the AIMgr if and when lots of AI traffic is generated) // Return a list of exits from a given runway // It is up to the calling function to check for non-zero size of returned array before use node_array_type FGGround::GetExits(string rwyID) { // FIXME - get a 07L or similar in here and we're stuffed!!! return(runways[atoi(rwyID.c_str())].exits); } void FGGround::RequestDeparture(PlaneRec plane, FGAIEntity* requestee) { // For now we'll just automatically clear all planes to the runway hold. // This communication needs to be delayed 20 sec or so from receiving the request. // Even if display=false we still need to start the timer in case display=true when communication starts. // We also need to bear in mind we also might have other outstanding communications, although for now we'll punt that issue! // FIXME - sort the above! // HACK - assume that anything requesting departure is new for now - FIXME LATER GroundRec* g = new GroundRec; g->plane = plane; g->planePtr = requestee; g->taxiRequestOutstanding = true; g->clearanceCounter = 0; g->cleared = false; g->incoming = false; // TODO - need to handle the next 3 as well //Point3D current_pos; //node* destination; //node* last_clearance; ground_traffic.push_back(g); } #if 0 void FGGround::NewArrival(plane_rec plane) { // What are we going to do here? // We need to start a new ground_rec and add the plane_rec to it // We need to decide what gate we are going to clear it to. // Then we need to add clearing it to that gate to the pending transmissions queue? - or simply transmit? // Probably simply transmit for now and think about a transmission queue later if we need one. // We might need one though in order to add a little delay for response time. ground_rec* g = new ground_rec; g->plane_rec = plane; g->current_pos = ConvertWGS84ToXY(plane.pos); g->node = GetNode(g->current_pos); // TODO - might need to sort out node/arc here AssignGate(g); g->cleared = false; ground_traffic.push_back(g); NextClearance(g); } void FGGround::NewContact(plane_rec plane) { // This is a bit of a convienience function at the moment and is likely to change. if(at a gate or apron) NewDeparture(plane); else NewArrival(plane); } void FGGround::NextClearance(ground_rec &g) { // Need to work out where we can clear g to. // Assume the pilot doesn't need progressive instructions // We *should* already have a gate or holding point assigned by the time we get here // but it wouldn't do any harm to check. // For now though we will hardwire it to clear to the final destination. } void FGGround::AssignGate(ground_rec &g) { // We'll cheat for now - since we only have the user's aircraft and a couple of airports implemented // we'll hardwire the gate! // In the long run the logic of which gate or area to send the plane to could be somewhat non-trivial. } #endif //0