#include <rtkADMMWaveletsConjugateGradientOperator.h>

Inheritance diagram for rtk::ADMMWaveletsConjugateGradientOperator< TOutputImage >:

Collaboration diagram for rtk::ADMMWaveletsConjugateGradientOperator< TOutputImage >:

Public Types
using	AddFilterType = itk::AddImageFilter< TOutputImage >

using	BackProjectionFilterPointer = typename BackProjectionFilterType::Pointer

using	BackProjectionFilterType = rtk::BackProjectionImageFilter< TOutputImage, TOutputImage >

using	DisplacedDetectorFilterType = rtk::DisplacedDetectorImageFilter< TOutputImage >

using	ForwardProjectionFilterPointer = typename ForwardProjectionFilterType::Pointer

using	ForwardProjectionFilterType = rtk::ForwardProjectionImageFilter< TOutputImage, TOutputImage >

using	MultiplyFilterType = itk::MultiplyImageFilter< TOutputImage >

using	Pointer = itk::SmartPointer< Self >

using	Self = ADMMWaveletsConjugateGradientOperator

using	Superclass = ConjugateGradientOperator< TOutputImage >

Public Types inherited from rtk::ConjugateGradientOperator< TOutputImage >
using	Pointer = itk::SmartPointer< Self >

using	Self = ConjugateGradientOperator

using	Superclass = itk::ImageToImageFilter< TOutputImage, TOutputImage >

Public Member Functions
virtual ::itk::LightObject::Pointer	CreateAnother () const

void	SetBackProjectionFilter (BackProjectionFilterType *_arg)

virtual void	SetBeta (float _arg)

void	SetForwardProjectionFilter (ForwardProjectionFilterType *_arg)

void	SetGeometry (ThreeDCircularProjectionGeometry *_arg)


virtual const char *	GetNameOfClass () const


virtual void	SetDisableDisplacedDetectorFilter (bool _arg)

virtual bool	GetDisableDisplacedDetectorFilter ()

Public Member Functions inherited from rtk::ConjugateGradientOperator< TOutputImage >
virtual ::itk::LightObject::Pointer	CreateAnother () const

virtual void	SetX (const TOutputImage *OutputImage)

Static Public Member Functions
static Pointer	New ()

Static Public Member Functions inherited from rtk::ConjugateGradientOperator< TOutputImage >
static Pointer	New ()

Protected Member Functions
	ADMMWaveletsConjugateGradientOperator ()

void	GenerateData () override

void	VerifyInputInformation () const override

	~ADMMWaveletsConjugateGradientOperator () override=default


void	GenerateInputRequestedRegion () override

void	GenerateOutputInformation () override

Protected Member Functions inherited from rtk::ConjugateGradientOperator< TOutputImage >
	ConjugateGradientOperator ()

	~ConjugateGradientOperator () override=default

Protected Attributes
AddFilterType::Pointer	m_AddFilter

BackProjectionFilterPointer	m_BackProjectionFilter

float	m_Beta

bool	m_DisableDisplacedDetectorFilter

DisplacedDetectorFilterType::Pointer	m_DisplacedDetectorFilter

ForwardProjectionFilterPointer	m_ForwardProjectionFilter

MultiplyFilterType::Pointer	m_MultiplyFilter

MultiplyFilterType::Pointer	m_ZeroMultiplyProjectionFilter

MultiplyFilterType::Pointer	m_ZeroMultiplyVolumeFilter

Detailed Description

template<typename TOutputImage>
class rtk::ADMMWaveletsConjugateGradientOperator< TOutputImage >

Implements the operator A used in the conjugate gradient step of ADMM reconstruction with wavelets regularization.

This filter implements the operator A used in the conjugate gradient step of a reconstruction method based on compressed sensing. The method attempts to find the f that minimizes || Rf -p ||_2^2 + alpha * || W(f) ||_1, with R the forward projection operator, p the measured projections, and W the Daubechies wavelets transform. Details on the method and the calculations can be found in

Mory, C., B. Zhang, V. Auvray, M. Grass, D. Schafer, F. Peyrin, S. Rit, P. Douek, and L. Boussel. "ECG-Gated C-Arm Computed Tomography Using L1 Regularization." In Proceedings of the 20th European Signal Processing Conference (EUSIPCO), 2728-32, 2012.

This filter takes in input f and outputs R_t R f + beta * W_t W f. The Daubechies wavelets being orthogonal, W_t happens to be the inverse of W, and therefore the filter outputs R_t R f + beta * f

Test:: rtkadmmWaveletstest.cxx

Author: Cyril Mory

Definition at line 98 of file rtkADMMWaveletsConjugateGradientOperator.h.