SpecFitSettingsWidgetRadio.qo.h

Go to the documentation of this file.

//# Copyright (C) 2005
//# Associated Universities, Inc. Washington DC, USA.
//#
//# This library is free software; you can redistribute it and/or modify it
//# under the terms of the GNU Library General Public License as published by
//# the Free Software Foundation; either version 2 of the License, or (at your
//# option) any later version.
//#
//# This library 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 Library General Public
//# License for more details.
//#
//# You should have received a copy of the GNU Library General Public License
//# along with this library; if not, write to the Free Software Foundation,
//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
//#
//# Correspondence concerning AIPS++ should be addressed as follows:
//#        Internet email: aips2-request@nrao.edu.
//#        Postal address: AIPS++ Project Office
//#                        National Radio Astronomy Observatory
//#                        520 Edgemont Road
//#                        Charlottesville, VA 22903-2475 USA
//#
#ifndef SPECFITSETTINGSWIDGETRADIO_QO_H
#define SPECFITSETTINGSWIDGETRADIO_QO_H

#include <QWidget>
#include <QProgressDialog>
#include <display/QtPlotter/SpecFitSettingsWidgetRadio.ui.h>
#include <display/QtPlotter/ProfileTaskFacilitator.h>
#include <display/QtPlotter/GaussianEstimateDialog.qo.h>
#include <casa/Containers/Record.h>
#include <display/Display/DisplayCoordinateSystem.h>

namespace casa {

        class ImageProfileFitter;
        class SpectralList;
        class SpectralElement;
        class SpecFit;
        class SpecFitThread;

        class SpecFitSettingsWidgetRadio : public QWidget, public ProfileTaskFacilitator {
                Q_OBJECT

        public:
                SpecFitSettingsWidgetRadio(QWidget *parent = 0);
                ~SpecFitSettingsWidgetRadio();
                void setUnits( QString units );
                void setDisplayYUnits( const QString& units );
                void setImageYUnits( const QString& units );
                void setCurveName( const QString& curveName );
                void addCurveName( const QString& curveName );
                void setRange(double start, double end );
                void reset();
                void pixelsChanged( int pixX, int pixY );

        signals:
                void gaussEstimateCountChanged( int count );

        private slots:
                void polyFitChanged( int state );
                void gaussCountChanged( int count );
                void fitRatioChanged( int count );
                void clean();
                void specLineFit();
                void setOutputLogFile();
                void viewOutputLogFile();
                void saveOutputChanged( int state );
                void specFitEstimateSpecified(double xValue,double yValue, bool centerPeak);
                void fitDone( bool newData = true );
                void cancelFit();
                void specifyGaussianEstimates();
                void gaussianEstimatesChanged();

        private:

                bool _constructFitter( std::shared_ptr<const casacore::ImageInterface<float> >& image,
                                        const casacore::String& region, const casacore::Record* const &regionPtr, const casacore::String& box,
                                        const casacore::String& chans, const casacore::String& stokes, const casacore::String& mask, const int axis,
                                        const casacore::uInt ngauss, const SpectralList& spectralList );
                void setEstimateValue( int row, int col, double val );
                bool isValidEstimate( QString& peakStr, QString& centerStr,
                                      QString& fwhmStr, QString& fixedStr, int rowIndex );
                SpectralList buildSpectralList( int nGauss, bool& validList );
                bool isValidFitSpecification( int gaussCount, bool polyFit );
                void setCanvas( QtCanvas* canvas );
                void doFit( float startVal, float endVal, uint gaussCount, bool fitPoly, int polyN );
                casacore::String getChannels( float startVal, float endVal, const casacore::Vector<float>& specValues ) const;
                void getFitBounds( float& startVal, float& endVal ) const;
                int getFitCount(int& startChannelIndex, int& endChannelIndex );
                void setFitEstimate( int row, double xValue, double yValue, bool centerPeak );
                void clearEstimates();
                void clear();
                void emptyCurveList();
                void resolveOutputLogFile( );
                void getConversion( const casacore::String& unitStr, bool& velocity, bool& wavelength ) const ;
                bool processFitResults(casacore::Vector<float>& xValues, casacore::Vector<float>& xValuesPix, const casacore::String& yUnit);
                void getEstimateStrings( int index, QString& peakStr, QString& centerStr, QString& fwhmStr ) const;
                bool processFitResultGaussian( const SpectralElement* solution,
                                               int index, QList<SpecFit*>& curves, const casacore::String& yUnit);
                bool processFitResultPolynomial( const SpectralElement* solution,
                                                 QList<SpecFit*>& curves);
                bool isInRange( double xValue ) const;
                QString settingsToString() const;
                double toPixels( double val, casacore::String units, casacore::SpectralCoordinate& coord ) const;
                void drawCurves( int pixelX, int pixelY );

                enum TableHeaders {PEAK,CENTER,FWHM,FIXED,END_COLUMN};
                Ui::SpecFitSettingsWidgetRadio ui;
                ImageProfileFitter* fitter;
                SpecFitThread* specFitThread;
                QProgressDialog progressDialog;
                GaussianEstimateDialog gaussEstimateDialog;
                casacore::Record regionRecord;


                QString outputLogPath;
                QString displayYUnits;
                QString imageYUnits;
                QList<QList<SpecFit*> > curveList;
                int SUM_FIT_INDEX;
                bool viewerSettingsWritten;
                volatile bool fitCancelled;
        };
}
#endif // SPECFITSETTINGSWIDGETRADIO_H