source: ntrip/trunk/BNC/src/rinex/dopplot.cpp@ 4767

// Part of BNC, a utility for retrieving decoding and
// converting GNSS data streams from NTRIP broadcasters.
//
// Copyright (C) 2007
// German Federal Agency for Cartography and Geodesy (BKG)
// http://www.bkg.bund.de
// Czech Technical University Prague, Department of Geodesy
// http://www.fsv.cvut.cz
//
// Email: euref-ip@bkg.bund.de
//
// 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, version 2.
//
// 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., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.

/* -------------------------------------------------------------------------
 * BKG NTRIP Client
 * -------------------------------------------------------------------------
 *
 * Class:      t_dopPlot
 *
 * Purpose:    Plot PDOP, GDOP, and number of satellites
 *
 * Author:     L. Mervart
 *
 * Created:    09-Sep-2012
 *
 * Changes:    
 *
 * -----------------------------------------------------------------------*/

#include <qwt_scale_draw.h>
#include <qwt_text.h>
#include <qwt_legend.h>
#include <qwt_plot_marker.h>

#include "dopplot.h"
#include "reqcanalyze.h"

//
//////////////////////////////////////////////////////////////////////////////
class t_scaleDrawTime : public QwtScaleDraw {
 public:
  t_scaleDrawTime() {}
  virtual QwtText label(double mjdX24) const {
    bncTime epoTime; epoTime.setmjd(mjdX24/24.0);
    return QwtText(epoTime.timestr(0,':').c_str());
  }
};

// Constructor
//////////////////////////////////////////////////////////////////////////////
t_dopPlot::t_dopPlot(QWidget* parent, t_obsStat* obsStat) 
: QwtPlot(parent) {

  setCanvasBackground(QColor(Qt::white));

  // Axes
  // ----
  setAxisScaleDraw(QwtPlot::xBottom, new t_scaleDrawTime());
  setAxisLabelRotation(QwtPlot::xBottom, -10.0);
  setAxisLabelAlignment(QwtPlot::xBottom, Qt::AlignLeft | Qt::AlignBottom);

  enableAxis(QwtPlot::yRight);

  QwtText textPDOP("PDOP");
  QFont   fontPDOP = textPDOP.font();
  fontPDOP.setPointSize(fontPDOP.pointSize()*0.8);
  textPDOP.setFont(fontPDOP);
  textPDOP.setColor(Qt::red);
  setAxisTitle(QwtPlot::yRight, textPDOP);
  setAxisScale(QwtPlot::yRight, 0,   6);

  QwtText textNumSat("# Sat");
  QFont   fontNumSat = textNumSat.font();
  fontNumSat.setPointSize(fontNumSat.pointSize()*0.8);
  textNumSat.setFont(fontNumSat);
  textNumSat.setColor(Qt::blue);
  setAxisTitle(QwtPlot::yLeft,  textNumSat);
  double maxNumSat = 20.0;
  if (obsStat) {
    for (int ii = 0; ii < obsStat->_numSat.size(); ii++) {
      if (maxNumSat < obsStat->_numSat[ii]) {
        maxNumSat = obsStat->_numSat[ii] + 5;
      } 
    }
  }
  setAxisScale(QwtPlot::yLeft,  0,  maxNumSat);

  // Legend
  // ------
  QwtLegend* legend = new QwtLegend;
  insertLegend(legend, QwtPlot::RightLegend);

  // Curves
  // ------
  if (obsStat) {

    QwtPlotCurve* curveNumSat = new QwtPlotCurve(textNumSat);
    curveNumSat->setXAxis(QwtPlot::xBottom);
    curveNumSat->setYAxis(QwtPlot::yLeft);
    curveNumSat->setSamples(obsStat->_mjdX24, obsStat->_numSat);
    curveNumSat->setPen(QPen(textNumSat.color()));
    curveNumSat->attach(this);

    QwtPlotCurve* curvePDOP = new QwtPlotCurve(textPDOP);
    curvePDOP->setXAxis(QwtPlot::xBottom);
    curvePDOP->setYAxis(QwtPlot::yRight);
    curvePDOP->setSamples(obsStat->_mjdX24, obsStat->_PDOP);
    curvePDOP->setPen(QPen(textPDOP.color()));
    curvePDOP->attach(this);
  }
  
  // Important !!!
  // -------------
  replot();
}