Changeset 4563 in ntrip


Ignore:
Timestamp:
Aug 28, 2012, 1:44:43 PM (12 years ago)
Author:
mervart
Message:
 
Location:
trunk/BNC/src/rinex
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/rinex/reqcanalyze.cpp

    r4562 r4563  
    244244  // Analyze the Multipath
    245245  // ---------------------
    246   QVector<t_polarPoint*>* dataMP1 = new QVector<t_polarPoint*>;
    247   QVector<t_polarPoint*>* dataMP2 = new QVector<t_polarPoint*>;
     246  QVector<t_polarPoint*>* dataMP1  = new QVector<t_polarPoint*>;
     247  QVector<t_polarPoint*>* dataMP2  = new QVector<t_polarPoint*>;
     248  QVector<t_polarPoint*>* dataSNR1 = new QVector<t_polarPoint*>;
     249  QVector<t_polarPoint*>* dataSNR2 = new QVector<t_polarPoint*>;
    248250
    249251  QMapIterator<QString, t_satStat> it(_satStat);
     
    252254    QString          prn     = it.key();
    253255    const t_satStat& satStat = it.value();
    254     analyzeMultipath(prn, satStat, xyz, obsFile->interval(), dataMP1, dataMP2);
     256    analyzeMultipathAndSNR(prn, satStat, xyz, obsFile->interval(),
     257                           dataMP1, dataMP2, dataSNR1, dataSNR2);
    255258  }
    256259
    257260  emit displayGraph(obsFile->fileName(), "MP1", dataMP1, "MP2", dataMP2, 2.0);
     261  emit displayGraph(obsFile->fileName(), "SNR11", dataSNR1, "SNR2", dataSNR2, 9.0);
    258262
    259263  if (_log) {
     
    341345// 
    342346////////////////////////////////////////////////////////////////////////////
    343 void t_reqcAnalyze::analyzeMultipath(const QString& prn,
    344                                      const t_satStat& satStat,
    345                                      const ColumnVector& xyz,
    346                                      double obsInterval,
    347                                      QVector<t_polarPoint*>* dataMP1,
    348                                      QVector<t_polarPoint*>* dataMP2) {
     347void t_reqcAnalyze::analyzeMultipathAndSNR(const QString& prn,
     348                                           const t_satStat& satStat,
     349                                           const ColumnVector& xyz,
     350                                           double obsInterval,
     351                                           QVector<t_polarPoint*>* dataMP1,
     352                                           QVector<t_polarPoint*>* dataMP2,
     353                                           QVector<t_polarPoint*>* dataSNR1,
     354                                           QVector<t_polarPoint*>* dataSNR2) {
    349355
    350356  const int chunkStep = int( 30.0 / obsInterval); // chunk step (30 sec) 
     
    359365    double mean1    = 0.0;
    360366    double mean2    = 0.0;
     367    double SNR1     = 0.0;
     368    double SNR2     = 0.0;
     369
    361370
    362371    for (int ii = 0; ii < numEpo; ii++) {
     
    365374      mean1 += anaObs->_MP1;
    366375      mean2 += anaObs->_MP2;
     376
     377      SNR1 = anaObs->_SNR1;
     378      SNR2 = anaObs->_SNR2;
    367379 
    368380      // Check Slip
     
    430442    // Add new Point
    431443    // -------------
    432     (*dataMP1) << (new t_polarPoint(az, zen, MP1));
    433     (*dataMP2) << (new t_polarPoint(az, zen, MP2));
     444    (*dataMP1)  << (new t_polarPoint(az, zen, MP1));
     445    (*dataMP2)  << (new t_polarPoint(az, zen, MP2));
     446    (*dataSNR1) << (new t_polarPoint(az, zen, SNR1));
     447    (*dataSNR2) << (new t_polarPoint(az, zen, SNR2));
    434448
    435449    if (_log) {
  • trunk/BNC/src/rinex/reqcanalyze.h

    r4558 r4563  
    9595
    9696  void analyzeFile(t_rnxObsFile* obsFile);
    97   void analyzeMultipath(const QString& prn,
    98                         const t_satStat& satStat,
    99                         const ColumnVector& xyz,
    100                         double obsInterval,
    101                         QVector<t_polarPoint*>* dataMP1,
    102                         QVector<t_polarPoint*>* dataMP2);
     97  void analyzeMultipathAndSNR(const QString& prn,
     98                              const t_satStat& satStat,
     99                              const ColumnVector& xyz,
     100                              double obsInterval,
     101                              QVector<t_polarPoint*>* dataMP1,
     102                              QVector<t_polarPoint*>* dataMP2,
     103                              QVector<t_polarPoint*>* dataSNR1,
     104                              QVector<t_polarPoint*>* dataSNR2);
    103105
    104106  QString                  _logFileName;
Note: See TracChangeset for help on using the changeset viewer.