#include <rtkADMMTotalVariationConeBeamReconstructionFilter.h>

Inheritance diagram for rtk::ADMMTotalVariationConeBeamReconstructionFilter< TOutputImage >:

Collaboration diagram for rtk::ADMMTotalVariationConeBeamReconstructionFilter< TOutputImage >:

Public Types
using	AddGradientsFilterType = itk::AddImageFilter< GradientImageType >

using	BackProjectionFilterPointer = typename BackProjectionFilterType::Pointer

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

using	BackProjectionType = typename Superclass::BackProjectionType

using	CGOperatorFilterType = rtk::ADMMTotalVariationConjugateGradientOperator< TOutputImage >

using	ConjugateGradientFilterType = rtk::ConjugateGradientImageFilter< TOutputImage >

using	CPUImageType = itk::Image< typename TOutputImage::PixelType, TOutputImage::ImageDimension >

using	DisplacedDetectorFilterType = rtk::DisplacedDetectorImageFilter< TOutputImage >

using	ForwardProjectionFilterPointer = typename ForwardProjectionFilterType::Pointer

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

using	ForwardProjectionType = typename Superclass::ForwardProjectionType

using	GatingWeightsFilterType = rtk::MultiplyByVectorImageFilter< TOutputImage >

typedef std::conditional< std::is_same< TOutputImage, CPUImageType >::value, itk::Image< VectorPixelType, TOutputImage::ImageDimension >, itk::CudaImage< VectorPixelType, TOutputImage::ImageDimension > >::type	GradientImageType

using	ImageDivergenceFilterType = rtk::BackwardDifferenceDivergenceImageFilter< GradientImageType, TOutputImage >

using	ImageGradientFilterType = ForwardDifferenceGradientImageFilter< TOutputImage, typename TOutputImage::ValueType, typename TOutputImage::ValueType, GradientImageType >

using	MultiplyGradientFilterType = itk::MultiplyImageFilter< GradientImageType >

using	MultiplyVolumeFilterType = itk::MultiplyImageFilter< TOutputImage >

using	Pointer = itk::SmartPointer< Self >

using	Self = ADMMTotalVariationConeBeamReconstructionFilter

typedef rtk::SoftThresholdTVImageFilter< GradientImageType >	SoftThresholdTVFilterType

using	SubtractGradientsFilterType = itk::SubtractImageFilter< GradientImageType >

using	SubtractVolumeFilterType = itk::SubtractImageFilter< TOutputImage >

using	Superclass = IterativeConeBeamReconstructionFilter< TOutputImage, TOutputImage >

using	VectorPixelType = itk::CovariantVector< typename TOutputImage::ValueType, TOutputImage::ImageDimension >

Public Types inherited from rtk::IterativeConeBeamReconstructionFilter< TOutputImage, TOutputImage >
using	BackProjectionFilterType = rtk::BackProjectionImageFilter< TOutputImage, VolumeType >

using	BackProjectionPointerType = typename BackProjectionFilterType::Pointer

enum	BackProjectionType

using	ConstPointer = itk::SmartPointer< const Self >

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

using	ForwardProjectionPointerType = typename ForwardProjectionFilterType::Pointer

enum	ForwardProjectionType

using	Pointer = itk::SmartPointer< Self >

using	Self = IterativeConeBeamReconstructionFilter

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

using	TClipImageType = itk::Image< double, VolumeType::ImageDimension >

using	VolumeType = TOutputImage

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

virtual float	GetAL_iterations ()

virtual float	GetAlpha ()

virtual float	GetBeta ()

virtual float	GetCG_iterations ()

const char *	GetNameOfClass () const override

virtual void	SetAL_iterations (float _arg)

virtual void	SetAlpha (float _arg)

virtual void	SetBeta (float _arg)

void	SetBetaForCurrentIteration (int iter)

virtual void	SetCG_iterations (float _arg)

void	SetGatingWeights (std::vector< float > weights)

virtual void	SetGeometry (ThreeDCircularProjectionGeometry *_arg)


virtual void	SetDisableDisplacedDetectorFilter (bool _arg)

virtual bool	GetDisableDisplacedDetectorFilter ()

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

const char *	GetNameOfClass () const override

virtual void	SetForwardProjectionFilter (ForwardProjectionType fwtype)

ForwardProjectionType	GetForwardProjectionFilter ()

virtual void	SetBackProjectionFilter (BackProjectionType bptype)

BackProjectionType	GetBackProjectionFilter ()

void	SetAttenuationMap (const VolumeType *attenuationMap)

VolumeType::ConstPointer	GetAttenuationMap ()

void	SetInferiorClipImage (const TClipImageType *inferiorClipImage)

TClipImageType::ConstPointer	GetInferiorClipImage ()

void	SetSuperiorClipImage (const TClipImageType *superiorClipImage)

TClipImageType::ConstPointer	GetSuperiorClipImage ()

virtual double	GetSigmaZero ()

virtual void	SetSigmaZero (double _arg)

virtual double	GetAlphaPSF ()

virtual void	SetAlphaPSF (double _arg)

virtual double	GetStepSize () const

virtual void	SetStepSize (double _arg)

Static Public Member Functions
static Pointer	New ()

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

Protected Member Functions
	ADMMTotalVariationConeBeamReconstructionFilter ()

void	GenerateData () override

void	VerifyInputInformation () const override

void	VerifyPreconditions () const override

	~ADMMTotalVariationConeBeamReconstructionFilter () override=default


void	GenerateInputRequestedRegion () override

void	GenerateOutputInformation () override

Protected Member Functions inherited from rtk::IterativeConeBeamReconstructionFilter< TOutputImage, TOutputImage >
virtual BackProjectionPointerType	InstantiateBackProjectionFilter (int bptype)

BackProjectionPointerType	InstantiateCudaBackProjection ()

BackProjectionPointerType	InstantiateCudaBackProjection ()

ForwardProjectionPointerType	InstantiateCudaForwardProjection ()

ForwardProjectionPointerType	InstantiateCudaForwardProjection ()

BackProjectionPointerType	InstantiateCudaRayCastBackProjection ()

BackProjectionPointerType	InstantiateCudaRayCastBackProjection ()

virtual ForwardProjectionPointerType	InstantiateForwardProjectionFilter (int fwtype)

BackProjectionPointerType	InstantiateJosephBackAttenuatedProjection ()

BackProjectionPointerType	InstantiateJosephBackAttenuatedProjection ()

ForwardProjectionPointerType	InstantiateJosephForwardAttenuatedProjection ()

ForwardProjectionPointerType	InstantiateJosephForwardAttenuatedProjection ()

BackProjectionPointerType	InstantiateZengBackProjection ()

BackProjectionPointerType	InstantiateZengBackProjection ()

ForwardProjectionPointerType	InstantiateZengForwardProjection ()

ForwardProjectionPointerType	InstantiateZengForwardProjection ()

	IterativeConeBeamReconstructionFilter ()

	~IterativeConeBeamReconstructionFilter () override=default

Protected Attributes
AddGradientsFilterType::Pointer	m_AddGradientsFilter

BackProjectionImageFilter< TOutputImage, TOutputImage >::Pointer	m_BackProjectionFilter

BackProjectionImageFilter< TOutputImage, TOutputImage >::Pointer	m_BackProjectionFilterForConjugateGradient

CGOperatorFilterType::Pointer	m_CGOperator

ConjugateGradientFilterType::Pointer	m_ConjugateGradientFilter

bool	m_DisableDisplacedDetectorFilter

DisplacedDetectorFilterType::Pointer	m_DisplacedDetectorFilter

ImageDivergenceFilterType::Pointer	m_DivergenceFilter

ForwardProjectionImageFilter< TOutputImage, TOutputImage >::Pointer	m_ForwardProjectionFilter

std::vector< float >	m_GatingWeights

GatingWeightsFilterType::Pointer	m_GatingWeightsFilter

ImageGradientFilterType::Pointer	m_GradientFilter1

ImageGradientFilterType::Pointer	m_GradientFilter2

bool	m_IsGated

MultiplyVolumeFilterType::Pointer	m_MultiplyFilter

SoftThresholdTVFilterType::Pointer	m_SoftThresholdFilter

SubtractGradientsFilterType::Pointer	m_SubtractFilter1

SubtractGradientsFilterType::Pointer	m_SubtractFilter2

SubtractVolumeFilterType::Pointer	m_SubtractVolumeFilter

MultiplyGradientFilterType::Pointer	m_ZeroMultiplyGradientFilter

MultiplyVolumeFilterType::Pointer	m_ZeroMultiplyVolumeFilter

Protected Attributes inherited from rtk::IterativeConeBeamReconstructionFilter< TOutputImage, TOutputImage >
double	m_AlphaPSF

BackProjectionType	m_CurrentBackProjectionConfiguration

ForwardProjectionType	m_CurrentForwardProjectionConfiguration

std::default_random_engine	m_DefaultRandomEngine

double	m_SigmaZero

double	m_StepSize

Private Attributes
unsigned int	m_AL_iterations

float	m_Alpha

float	m_Beta

unsigned int	m_CG_iterations

ThreeDCircularProjectionGeometry::Pointer	m_Geometry

Additional Inherited Members
Protected Types inherited from rtk::IterativeConeBeamReconstructionFilter< TOutputImage, TOutputImage >
using	CPUImageType = typename itk::Image< typename TOutputImage ::PixelType, TOutputImage ::ImageDimension >

using	EnableCudaScalarAndVectorType = typename std::enable_if< !std::is_same_v< CPUImageType, ImageType > &&std::is_same_v< typename itk::PixelTraits< typename ImageType::PixelType >::ValueType, float > &&(itk::PixelTraits< typename ImageType::PixelType >::Dimension==1\|\|itk::PixelTraits< typename ImageType::PixelType >::Dimension==2\|\|itk::PixelTraits< typename ImageType::PixelType >::Dimension==3)>::type

using	DisableCudaScalarAndVectorType = typename std::enable_if< std::is_same_v< CPUImageType, ImageType >\|\|!std::is_same_v< typename itk::PixelTraits< typename ImageType::PixelType >::ValueType, float >\|\|(itk::PixelTraits< typename ImageType::PixelType >::Dimension !=1 &&itk::PixelTraits< typename ImageType::PixelType >::Dimension !=2 &&itk::PixelTraits< typename ImageType::PixelType >::Dimension !=3)>::type

using	EnableCudaScalarType = typename std::enable_if< !std::is_same_v< CPUImageType, ImageType > &&std::is_same_v< typename itk::PixelTraits< typename ImageType::PixelType >::ValueType, float > &&itk::PixelTraits< typename ImageType::PixelType >::Dimension==1 >::type

using	DisableCudaScalarType = typename std::enable_if< std::is_same_v< CPUImageType, ImageType >\|\|!std::is_same_v< typename itk::PixelTraits< typename ImageType::PixelType >::ValueType, float >\|\|itk::PixelTraits< typename ImageType::PixelType >::Dimension !=1 >::type

using	EnableVectorType = typename std::enable_if< itk::PixelTraits< typename ImageType::PixelType >::Dimension !=1 >::type

using	DisableVectorType = typename std::enable_if< itk::PixelTraits< typename ImageType::PixelType >::Dimension==1 >::type

Detailed Description

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

Implements the ADMM reconstruction with total variation regularization.

This filter implements a reconstruction method based on compressed sensing. The method attempts to find the f that minimizes || sqrt(D) (Rf -p) ||_2^2 + alpha * TV(f), with R the forward projection operator, p the measured projections, D the displaced detector weighting matrix, and TV the total variation. Details on the method and the calculations can be found on page 48 of

Mory, C. "Cardiac C-Arm Computed Tomography", PhD thesis, 2014. https://hal.inria.fr/tel-00985728/document

\( f_{k+1} \) is obtained by linear conjugate solving the following:

\[ ( R^T R + \beta \nabla^T \nabla ) f = R^T p + \beta \nabla^T ( g_k + d_k ) \]

\( g_{k+1} \) is obtained by soft thresholding:

\[ g_{k+1} = ST_{ \frac{\alpha}{2 \beta} } ( \nabla f_{k+1} - d_k ) \]

\( d_{k+1} \) is a simple subtraction:

\[ d_{k+1} = g_{k+1} - ( \nabla f_{k+1} - d_k) \]

Test:: rtkadmmtotalvariationtest.cxx

Author: Cyril Mory

Definition at line 135 of file rtkADMMTotalVariationConeBeamReconstructionFilter.h.