#include <qxrdpolarintensityfitter.h>

Public Member Functions
	QxrdPolarIntensityFitter (QcepDoubleImageDataPtr data, int colNum)

	~QxrdPolarIntensityFitter ()

int	fit ()

double	fittedAverage ()

double	fittedAmplitude ()

Private Attributes
QcepDoubleImageDataPtr	m_Data

int	m_ColNum

double	m_FittedAverage

double	m_FittedAmplitude

Detailed Description

Definition at line 7 of file qxrdpolarintensityfitter.h.

Constructor & Destructor Documentation

QxrdPolarIntensityFitter::QxrdPolarIntensityFitter	(	QcepDoubleImageDataPtr	data,
		int	colNum
	)

Definition at line 6 of file qxrdpolarintensityfitter.cpp.

                                                                                           :
   m_Data(data),
   m_ColNum(colNum),
   m_FittedAverage(0.0),
   m_FittedAmplitude(0.0)
 {
 }

QxrdPolarIntensityFitter::~QxrdPolarIntensityFitter ( )

Definition at line 14 of file qxrdpolarintensityfitter.cpp.

15 {

16 }

Member Function Documentation

int QxrdPolarIntensityFitter::fit ( )

Definition at line 28 of file qxrdpolarintensityfitter.cpp.

References m_ColNum, m_Data, m_FittedAmplitude, and m_FittedAverage.

Referenced by QxrdPolarNormalization::executeCol().

 {
   if (m_Data) {
     int nRows = m_Data->get_Height();
 
     double vStart = m_Data->get_VStart();
     double vStep  = m_Data->get_VStep();
 
     double sumn = 0, sumaa = 0, sumab = 0, sumbb = 0, sumva = 0, sumvb = 0;
 
     for (int i=0; i<nRows; i++) {
       double x = vStart + i*vStep;
       double a = 1;
       double b = cos(2.0*x*M_PI/180.0);
       double v = m_Data->value(m_ColNum, i);
 
       if (v == v) {
         sumn  += 1;
         sumaa += a*a;
         sumab += a*b;
         sumbb += b*b;
 
         sumva += v*a;
         sumvb += v*b;
       }
     }
 
     if (sumn > 5) {
       QMatrix4x4 m; // Initialized as unit matrix
 
       m(0,0) = sumaa;
       m(1,0) = sumab;
       m(0,1) = sumab;
       m(1,1) = sumbb;
 
       bool invertible;
 
       QMatrix4x4 inv = m.inverted(&invertible);
 
       if (invertible) {
         m_FittedAverage   = inv(0,0)*sumva + inv(0,1)*sumvb;
         m_FittedAmplitude = inv(1,0)*sumva + inv(1,1)*sumvb;
       } else {
         m_FittedAverage   = sumva/sumaa;
         m_FittedAmplitude = 0;
       }
     }
   }
 
   return true;
 }