Changeset 6653 in ntrip


Ignore:
Timestamp:
Mar 12, 2015, 1:26:34 PM (7 years ago)
Author:
stuerze
Message:

sinex tro file support added,
troposphere results in overall ppp logfile added

Location:
trunk/BNC/src
Files:
2 added
8 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/PPP/pppClient.cpp

    r6463 r6653  
    336336    copy(&_filter->Q().data()[0], &_filter->Q().data()[6], _output->_covMatrix);
    337337
     338    _output->_trp0     = t_tropo::delay_saast(_staRover->xyzApr(), M_PI/2.0);
     339    _output->_trp      = _filter->trp();
     340    _output->_trpStdev = _filter->trpStdev();
     341
    338342    _output->_numSat     = _filter->numSat();
    339343    _output->_pDop       = _filter->PDOP();
  • trunk/BNC/src/PPP/pppFilter.h

    r6011 r6653  
    2828  double PDOP() const {return _dop.P;}
    2929  double GDOP() const {return _dop.G;}
     30  double trp() const {
     31    const std::vector<t_pppParam*>& par = _parlist->params();
     32    for (unsigned ii = 0; ii < par.size(); ++ii) {
     33      if (par[ii]->type() == t_pppParam::trp) {
     34        return x()[ii];
     35      }
     36    }
     37    return 0.0;
     38  };
     39  double trpStdev() const {
     40    const std::vector<t_pppParam*>& par = _parlist->params();
     41    for (unsigned ii = 0; ii < par.size(); ++ii) {
     42      if (par[ii]->type() == t_pppParam::trp) {
     43        return sqrt(Q()[ii][ii]);
     44      }
     45    }
     46    return 0.0;
     47  };
    3048
    3149 private:
     
    5068    double G;
    5169  };
     70
    5271  t_irc processSystem(const std::vector<t_lc::type>& LCs,
    5372                      const std::vector<t_pppSatObs*>& obsVector);
     
    6079
    6180  void cmpDOP(const std::vector<t_pppSatObs*>& obsVector);
     81
     82  void getTrp(void);
    6283
    6384  void predictCovCrdPart(const SymmetricMatrix& QFltOld);
  • trunk/BNC/src/PPP_free/pppClient.cpp

    r6634 r6653  
    145145    output->_numSat      = _filter->numSat();
    146146    output->_pDop        = _filter->PDOP();
     147    output->_trp0        = _filter->delay_saast(M_PI/2.0);
     148    output->_trp         = _filter->trp();
    147149  }
    148150  else {
  • trunk/BNC/src/PPP_free/pppFilter.h

    r6634 r6653  
    169169    return 0.0;
    170170  }
     171  double trpStdev() const {
     172    for (unsigned ii = 0; ii < _params.size(); ++ii) {
     173      t_pppParam* pp = _params[ii];
     174      if (pp->type() == t_pppParam::TROPO) {
     175        return sqrt(Q()[ii][ii]);
     176      }
     177    }
     178    return 0.0;
     179  }
    171180  double Glonass_offset() const {
    172181    for (int ii = 0; ii < _params.size(); ++ii) {
  • trunk/BNC/src/pppInclude.h

    r6465 r6653  
    3131  double       _covMatrix[6];
    3232  double       _neu[3]; 
     33  double       _trp0;
     34  double       _trp;
     35  double       _trpStdev;
    3336  int          _numSat;       
    3437  double       _pDop;         
  • trunk/BNC/src/pppRun.cpp

    r6608 r6653  
    5353#include "bncsettings.h"
    5454#include "bncoutf.h"
     55#include "bncsinextro.h"
    5556#include "rinex/rnxobsfile.h"
    5657#include "rinex/rnxnavfile.h"
     
    159160    }
    160161    else {
    161       snxtroFileSkl.replace("${STATION}", roverName);
    162     }
    163     _snxtroFile = new bncoutf(snxtroFileSkl, "1 day", 0);
     162      snxtroFileSkl.replace("${STATION}", roverName.left(4));
     163    }
     164    int samplSnxTro = settings.value("PPP/snxtroSampl").toInt();
     165    _snxtroFile = new bncSinexTro(_opt, snxtroFileSkl, "1 day", samplSnxTro);
    164166  }
    165167}
     
    272274          << " Y = "  << setprecision(4) << output._xyzRover[1]
    273275          << " Z = "  << setprecision(4) << output._xyzRover[2]
    274           << " NEU: " << setprecision(4) << output._neu[0]     
    275           << " " << setprecision(4) << output._neu[1]     
    276           << " " << setprecision(4) << output._neu[2];
     276          << " NEU: " << showpos << setw(8) << setprecision(4) << output._neu[0]
     277          << " "      << showpos << setw(8) << setprecision(4) << output._neu[1]
     278          << " "      << showpos << setw(8) << setprecision(4) << output._neu[2]
     279          << " TRP: " << showpos << setw(8) << setprecision(4) << output._trp0
     280          << " "      << showpos << setw(8) << setprecision(4) << output._trp;
    277281    }
    278282
     
    294298
    295299    if (_snxtroFile && output._epoTime.valid()) {
    296      // _snxtroFile->write(output._epoTime.gpsw(), output._epoTime.gpssec(),QString("TEST"));
     300      _snxtroFile->write(staID, int(output._epoTime.gpsw()), output._epoTime.gpssec(),
     301                  output._trp0 + output._trp, output._trpStdev);
    297302    }
    298303
  • trunk/BNC/src/pppRun.h

    r6608 r6653  
    1414class t_corrFile;
    1515class bncoutf;
     16class bncSinexTro;
    1617
    1718namespace BNC_PPP {
     
    7273  bncoutf*               _logFile;
    7374  bncoutf*               _nmeaFile;
    74   bncoutf*               _snxtroFile;
     75  bncSinexTro*           _snxtroFile;
    7576};
    7677
  • trunk/BNC/src/pppWidgets.cpp

    r6646 r6653  
    213213  _minEle->setValue(settings.value(_minEle->objectName()).toInt());
    214214  _corrWaitTime->setValue(settings.value(_corrWaitTime->objectName()).toInt());
     215  _snxtroSampl->setValue(settings.value(_snxtroSampl->objectName()).toInt());
    215216
    216217  // Table with stations
     
    264265  settings.setValue(_nmeaPort    ->objectName(), _nmeaPort    ->text());
    265266  settings.setValue(_snxtroFile  ->objectName(), _snxtroFile  ->text());
     267  settings.setValue(_snxtroSampl ->objectName(), _snxtroSampl ->value());
    266268  settings.setValue(_lcGPS       ->objectName(), _lcGPS       ->currentText());
    267269  settings.setValue(_lcGLONASS   ->objectName(), _lcGLONASS   ->currentText());
Note: See TracChangeset for help on using the changeset viewer.