RdTransformAds_bLogFile.cc

Go to the documentation of this file.

#include "RdTransformAds_bLogFile.h"
#include "AirplaneUtil.h"
#include "Ads_bDataBase.h"
#include "Airplane.h"
#include "Flight.h"

#include <limits> // setprecision
#include <iomanip>  // setprecision

using namespace evt;
using namespace revt;
using namespace std;
using namespace utl;
using namespace fwk;
using namespace RdAirplane;
using namespace boost;
namespace fs = filesystem;

RdTransformAds_bLogFile::RdTransformAds_bLogFile() :
  _outputFile_(),
  _ads_bDataBase_(NULL) {}

RdTransformAds_bLogFile::~RdTransformAds_bLogFile() {}

VModule::ResultFlag RdTransformAds_bLogFile::Init() {
  
  _globalCoordinateSystem_ = det::Detector::GetInstance().GetReferenceCoordinateSystem();
  string outputDirectory;
  
  CentralConfig* centralConfig = CentralConfig::GetInstance();
  Branch timingTopBranch =  centralConfig->GetTopBranch("RdTransformAds_bLogFile");
  timingTopBranch.GetChild("OutputDirectoy").GetData(outputDirectory);

  //reading the ads-b data that is  grouped by airplane and by flight over the AERA field
  _ads_bDataBase_ = Ads_bDataBase::getSharedDataBase();
  INFO("Assuming a time shift of the capture time of the ADS-B signal and the real position measurement on the airplane of " + lexical_cast<string>(_ads_bDataBase_->getTimeShift().GetInterval()/ms) + "ms.");
  
  _outputFile_ = fs::path(outputDirectory);
  if(!fs::exists(_outputFile_)) {
    ERROR("The output directory does not exist: " + _outputFile_.string());
    return eFailure;
  }
  _outputFile_+="/adsbData.dat";

  if(fs::exists(_outputFile_)) {
    ERROR("The output file already exists: " + _outputFile_.string());
    return eFailure;
  }
  INFO("Data conversion started.");

  ofstream outputFileHandle;
  try {
    outputFileHandle.open(_outputFile_.c_str());
    if(!outputFileHandle.is_open()) {
      ERROR("Could not write data.");
      return eFailure;
    }
    outputFileHandle
      << "Flight Id| "
      << "Time (Unix-Timestamp)| "
      << "Latitude (UTM) (degree)| "
      << "Longitude (UTM) (degree)| "
      << "Pressure Altitude (UTM) (m)| "
      << "Corrected Pressure Altitude (m)| "
      << "GPS Altitude (m)| "
      << endl
      << setprecision(numeric_limits<long double>::digits10 + 1);
    const Ads_bDataBase::AirplaneMap& airplanes = _ads_bDataBase_->getAirplanes();
    //iterate over all airplanes
    for(Ads_bDataBase::AirplaneMap::const_iterator airplaneIter = airplanes.begin();
        airplaneIter != airplanes.end();
        ++airplaneIter) {
      const Airplane* airplane = airplaneIter->second;
      const Airplane::FlightList& flights = airplane->getFlights();
      for(Airplane::FlightList::const_iterator flightIter = flights.begin();
          flightIter != flights.end();
          ++flightIter) {
        const Flight* flight = *flightIter;
        const Flight::Ads_bEventSet& ads_bEvents = flight->getAdsbEvents();
        // iterate over all fly-overs
        for(Flight::Ads_bEventSet::iterator ads_bEventIter = ads_bEvents.begin();
            ads_bEventIter != ads_bEvents.end();
            ++ads_bEventIter) {
          const Ads_bEvent* ads_bEvent = *ads_bEventIter;
          Ads_bEvent::AltitudeType altitudeType = ads_bEvent->getAltitudeType();
          double pressureAltitude = altitudeType == Ads_bEvent::eGPSAltitude ? -1 : ads_bEvent->getAltitude(Ads_bEvent::ePressureAltitude);
          double correctedPressureAltitude = altitudeType == Ads_bEvent::eGDASCorrectedPressureAltitude ? ads_bEvent->getAltitude(Ads_bEvent::eGDASCorrectedPressureAltitude) : -1;
          double gpsAltitude = altitudeType == Ads_bEvent::eGPSAltitude ? ads_bEvent->getAltitude(Ads_bEvent::eGPSAltitude) : -1;
          const UTMPoint& utmCoordinates = ads_bEvent->getCoordinates();
          tuple<double, double, double> geodeticCoordinates = utmCoordinates.GetGeodeticCoordinates();
          double latitude = geodeticCoordinates.get<0>();
          double longitude = geodeticCoordinates.get<1>();

          outputFileHandle
            << flight->getId() << " "
            << UTCDateTime(ads_bEvent->getCaptureTime()).GetUnixSecond() << " "
            << latitude/degree << " "
            << longitude/degree << " "
            << pressureAltitude/m << " "
            << correctedPressureAltitude/m << " "
            << gpsAltitude/m
            << endl;

        }
      }
    }
  } catch(...) {
    outputFileHandle.close();
    throw;
  }
  outputFileHandle.close();
  INFO("Data conversion finished. (Aaaaand Offline might have killed itself. If so hit Ctr+C)");
  return eSuccess;
}

VModule::ResultFlag RdTransformAds_bLogFile::Run(Event& event) {
 event.HasREvent(); // suppress unused variable warning
 return eSuccess;
}


VModule::ResultFlag RdTransformAds_bLogFile::Finish() {
  return eSuccess;
}