animaMultiT2EPGRelaxometryCostFunction.h

Go to the documentation of this file.

#pragma once

#include <itkSingleValuedCostFunction.h>
#include "AnimaRelaxometryExport.h"

namespace anima
{
    
class ANIMARELAXOMETRY_EXPORT MultiT2EPGRelaxometryCostFunction :
public itk::SingleValuedCostFunction
{
public:
    typedef MultiT2EPGRelaxometryCostFunction Self;
    typedef itk::SingleValuedCostFunction   Superclass;
    typedef itk::SmartPointer<Self>         Pointer;
    typedef itk::SmartPointer<const Self>   ConstPointer;

    itkNewMacro(Self);

    itkTypeMacro(MultiT2EPGRelaxometryCostFunction, Superclass);

    typedef Superclass::MeasureType MeasureType;
    typedef Superclass::DerivativeType DerivativeType;
    typedef Superclass::ParametersType ParametersType;
    typedef std::vector < std::complex <double> > ComplexVectorType;
    typedef std::vector <ComplexVectorType> MatrixType;

    virtual MeasureType GetValue(const ParametersType & parameters) const ITK_OVERRIDE;
    virtual void GetDerivative(const ParametersType & parameters, DerivativeType & derivative) const ITK_OVERRIDE {} // No derivative

    itkSetMacro(EchoSpacing, double)
    itkSetMacro(ExcitationFlipAngle, double)

    void SetT2RelaxometrySignals(ParametersType &relaxoSignals) {m_T2RelaxometrySignals = relaxoSignals;}
    void SetT2FlipAngles(std::vector <double> & flipAngles) {m_T2FlipAngles = flipAngles;}

    itkSetMacro(T1Value, double)
    itkSetMacro(M0Value, double)

    void SetT2Values(std::vector <double> &values) {m_T2Values = values;}
    void SetT2Weights(ParametersType &weights) {m_T2Weights = weights;}

    unsigned int GetNumberOfParameters() const ITK_OVERRIDE
    {
        return 1;
    }

protected:
    MultiT2EPGRelaxometryCostFunction()
    {
        m_T1Value = 1;
        m_M0Value = 1;

        m_EchoSpacing = 1;
    }

    virtual ~MultiT2EPGRelaxometryCostFunction() {}

private:
    MultiT2EPGRelaxometryCostFunction(const Self&); //purposely not implemented
    void operator=(const Self&); //purposely not implemented

    double m_EchoSpacing;
    ParametersType m_T2RelaxometrySignals;

    double m_ExcitationFlipAngle;
    std::vector <double> m_T2FlipAngles;

    std::vector <double> m_T2Values;
    ParametersType m_T2Weights;

    double m_T1Value, m_M0Value;
};
    
} // end namespace anima