#include <rtkFourDROOSTERConeBeamReconstructionFilter.h>

Inheritance diagram for rtk::FourDROOSTERConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >:

Collaboration diagram for rtk::FourDROOSTERConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >:

Public Types
using	AddFilterType = itk::AddImageFilter< VolumeSeriesType, VolumeSeriesType >

using	AverageOutOfROIFilterType = std::conditional_t< std::is_same_v< VolumeSeriesType, CPUVolumeSeriesType >, AverageOutOfROIImageFilter< VolumeSeriesType, VolumeType >, CudaAverageOutOfROIImageFilter >

using	BackProjectionType = typename Superclass::BackProjectionType

using	CovariantVectorForSpatialGradient = itk::CovariantVector< typename VolumeSeriesType::ValueType, VolumeSeriesType::ImageDimension - 1 >

using	CovariantVectorForTemporalGradient = itk::CovariantVector< typename VolumeSeriesType::ValueType, 1 >

using	CPUVolumeSeriesType = typename itk::Image< typename VolumeSeriesType::PixelType, VolumeSeriesType::ImageDimension >

using	DVFImageType = typename VolumeSeriesType::template RebindImageType< DVFVectorType, VolumeSeriesType::ImageDimension - 1 >

using	DVFSequenceImageType = typename VolumeSeriesType::template RebindImageType< DVFVectorType, VolumeSeriesType::ImageDimension >

using	DVFVectorType = CovariantVectorForSpatialGradient

using	ForwardProjectionType = typename Superclass::ForwardProjectionType

using	FourDCGFilterType = rtk::FourDConjugateGradientConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >

using	Pointer = itk::SmartPointer< Self >

using	ResampleFilterType = itk::ResampleImageFilter< VolumeType, VolumeType >

using	Self = FourDROOSTERConeBeamReconstructionFilter

using	SpatialGradientImageType = typename VolumeSeriesType::template RebindImageType< CovariantVectorForSpatialGradient, VolumeSeriesType::ImageDimension >

using	SpatialTVDenoisingFilterType = rtk::TotalVariationDenoiseSequenceImageFilter< VolumeSeriesType >

using	SpatialWaveletsDenoisingFilterType = rtk::DaubechiesWaveletsDenoiseSequenceImageFilter< VolumeSeriesType >

using	SubtractFilterType = itk::SubtractImageFilter< VolumeSeriesType, VolumeSeriesType >

using	Superclass = rtk::IterativeConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >

using	TemporalGradientImageType = typename VolumeSeriesType::template RebindImageType< CovariantVectorForTemporalGradient, VolumeSeriesType::ImageDimension >

using	TemporalL0DenoisingFilterType = rtk::LastDimensionL0GradientDenoisingImageFilter< VolumeSeriesType >

using	TemporalTVDenoisingFilterType = std::conditional_t< std::is_same_v< VolumeSeriesType, CPUVolumeSeriesType >, TotalVariationDenoisingBPDQImageFilter< VolumeSeriesType, TemporalGradientImageType >, CudaLastDimensionTVDenoisingImageFilter >

using	ThresholdFilterType = itk::ThresholdImageFilter< VolumeSeriesType >

using	TNVDenoisingFilterType = rtk::TotalNuclearVariationDenoisingBPDQImageFilter< VolumeSeriesType, SpatialGradientImageType >

using	UnwarpSequenceFilterType = rtk::UnwarpSequenceImageFilter< VolumeSeriesType, DVFSequenceImageType, VolumeType, DVFImageType >

using	VolumeType = ProjectionStackType

using	WarpSequenceFilterType = rtk::WarpSequenceImageFilter< VolumeSeriesType, DVFSequenceImageType, VolumeType, DVFImageType >

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

using	BackProjectionPointerType = typename BackProjectionFilterType::Pointer

enum	BackProjectionType

using	ConstPointer = itk::SmartPointer< const Self >

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

using	ForwardProjectionPointerType = typename ForwardProjectionFilterType::Pointer

enum	ForwardProjectionType

using	Pointer = itk::SmartPointer< Self >

using	Self = IterativeConeBeamReconstructionFilter

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

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

using	VolumeType = ProjectionStackType

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

virtual int	GetCG_iterations ()

virtual bool	GetComputeInverseWarpingByConjugateGradient ()

virtual bool	GetCudaConjugateGradient ()

virtual float	GetGammaTNV ()

virtual float	GetGammaTVSpace ()

virtual float	GetGammaTVTime ()

virtual const ThreeDCircularProjectionGeometry *	GetGeometry () const

virtual int	GetL0_iterations ()

virtual float	GetLambdaL0Time ()

virtual int	GetMainLoop_iterations ()

const char *	GetNameOfClass () const override

virtual bool	GetPerformL0TemporalDenoising ()

virtual bool	GetPerformMotionMask ()

virtual bool	GetPerformPositivity ()

virtual bool	GetPerformTNVDenoising ()

virtual bool	GetPerformTVSpatialDenoising ()

virtual bool	GetPerformTVTemporalDenoising ()

virtual bool	GetPerformWarping ()

virtual bool	GetPerformWaveletsSpatialDenoising ()

virtual float	GetPhaseShift ()

virtual float	GetSoftThresholdWavelets ()

virtual int	GetTV_iterations ()

virtual bool	GetUseNearestNeighborInterpolationInWarping ()

virtual void	SetCG_iterations (int _arg)

virtual void	SetComputeInverseWarpingByConjugateGradient (bool _arg)

virtual void	SetCudaConjugateGradient (bool _arg)

virtual void	SetGammaTNV (float _arg)

virtual void	SetGammaTVSpace (float _arg)

virtual void	SetGammaTVTime (float _arg)

virtual void	SetGeometry (const ThreeDCircularProjectionGeometry *_arg)

virtual void	SetL0_iterations (int _arg)

virtual void	SetLambdaL0Time (float _arg)

virtual void	SetMainLoop_iterations (int _arg)

virtual void	SetPerformL0TemporalDenoising (bool _arg)

virtual void	SetPerformMotionMask (bool _arg)

virtual void	SetPerformPositivity (bool _arg)

virtual void	SetPerformTNVDenoising (bool _arg)

virtual void	SetPerformTVSpatialDenoising (bool _arg)

virtual void	SetPerformTVTemporalDenoising (bool _arg)

virtual void	SetPerformWarping (bool _arg)

virtual void	SetPerformWaveletsSpatialDenoising (bool _arg)

virtual void	SetPhaseShift (float _arg)

virtual void	SetSignal (const std::vector< double > signal)

virtual void	SetSoftThresholdWavelets (float _arg)

virtual void	SetTV_iterations (int _arg)

virtual void	SetUseNearestNeighborInterpolationInWarping (bool _arg)

virtual void	SetWeights (const itk::Array2D< float > _arg)


void	SetInputVolumeSeries (const VolumeSeriesType *VolumeSeries)

VolumeSeriesType::ConstPointer	GetInputVolumeSeries ()


void	SetInputProjectionStack (const ProjectionStackType *Projection)

ProjectionStackType::Pointer	GetInputProjectionStack ()


void	SetMotionMask (const VolumeType *mask)

VolumeType::Pointer	GetMotionMask ()


void	SetDisplacementField (const DVFSequenceImageType *DVFs)

void	SetInverseDisplacementField (const DVFSequenceImageType *DVFs)

DVFSequenceImageType::Pointer	GetDisplacementField ()

DVFSequenceImageType::Pointer	GetInverseDisplacementField ()


virtual void	SetDisableDisplacedDetectorFilter (bool _arg)

virtual bool	GetDisableDisplacedDetectorFilter ()


virtual void	SetUseCudaCyclicDeformation (bool _arg)

virtual bool	GetUseCudaCyclicDeformation ()


virtual unsigned int	GetNumberOfLevels ()

virtual void	SetNumberOfLevels (unsigned int _arg)


virtual unsigned int	GetOrder ()

virtual void	SetOrder (unsigned int _arg)

Public Member Functions inherited from rtk::IterativeConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >
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< VolumeSeriesType, ProjectionStackType >
static Pointer	New ()

Protected Member Functions
	FourDROOSTERConeBeamReconstructionFilter ()

void	GenerateData () override

void	GenerateInputRequestedRegion () override

void	GenerateOutputInformation () override

void	VerifyInputInformation () const override

void	VerifyPreconditions () const override

	~FourDROOSTERConeBeamReconstructionFilter () override=default

Protected Member Functions inherited from rtk::IterativeConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >
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
AddFilterType::Pointer	m_AddFilter

AverageOutOfROIFilterType::Pointer	m_AverageOutOfROIFilter

int	m_CG_iterations

bool	m_ComputeInverseWarpingByConjugateGradient

bool	m_CudaConjugateGradient

bool	m_DimensionsProcessedForTVSpace [VolumeSeriesType::ImageDimension]

bool	m_DimensionsProcessedForTVTime [VolumeSeriesType::ImageDimension]

bool	m_DisableDisplacedDetectorFilter

itk::ImageToImageFilter< VolumeSeriesType, VolumeSeriesType >::Pointer	m_DownstreamFilter

FourDCGFilterType::Pointer	m_FourDCGFilter

float	m_GammaTNV

float	m_GammaTVSpace

float	m_GammaTVTime

rtk::ThreeDCircularProjectionGeometry::ConstPointer	m_Geometry

WarpSequenceFilterType::Pointer	m_InverseWarp

int	m_L0_iterations

TemporalL0DenoisingFilterType::Pointer	m_L0DenoisingTime

float	m_LambdaL0Time

int	m_MainLoop_iterations

unsigned int	m_NumberOfLevels

unsigned int	m_Order

bool	m_PerformL0TemporalDenoising

bool	m_PerformMotionMask

bool	m_PerformPositivity

bool	m_PerformTNVDenoising

bool	m_PerformTVSpatialDenoising

bool	m_PerformTVTemporalDenoising

bool	m_PerformWarping

bool	m_PerformWaveletsSpatialDenoising

float	m_PhaseShift

ThresholdFilterType::Pointer	m_PositivityFilter

ResampleFilterType::Pointer	m_ResampleFilter

std::vector< double >	m_Signal

float	m_SoftThresholdWavelets

SubtractFilterType::Pointer	m_SubtractFilter

TNVDenoisingFilterType::Pointer	m_TNVDenoising

int	m_TV_iterations

SpatialTVDenoisingFilterType::Pointer	m_TVDenoisingSpace

TemporalTVDenoisingFilterType::Pointer	m_TVDenoisingTime

UnwarpSequenceFilterType::Pointer	m_Unwarp

bool	m_UseCudaCyclicDeformation

bool	m_UseNearestNeighborInterpolationInWarping

WarpSequenceFilterType::Pointer	m_Warp

SpatialWaveletsDenoisingFilterType::Pointer	m_WaveletsDenoisingSpace

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

BackProjectionType	m_CurrentBackProjectionConfiguration

ForwardProjectionType	m_CurrentForwardProjectionConfiguration

std::default_random_engine	m_DefaultRandomEngine

double	m_SigmaZero

double	m_StepSize

Additional Inherited Members
Protected Types inherited from rtk::IterativeConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >
using	CPUImageType = typename itk::Image< typename ProjectionStackType::PixelType, ProjectionStackType::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 VolumeSeriesType, typename ProjectionStackType>
class rtk::FourDROOSTERConeBeamReconstructionFilter< VolumeSeriesType, ProjectionStackType >

Implements 4D RecOnstructiOn using Spatial and TEmporal Regularization (short 4D ROOSTER)

See the reference paper: "Cardiac C-arm computed tomography using a 3D + time ROI reconstruction method with spatial and temporal regularization" by Mory et al.

4D ROOSTER reconstruction consists in performing 4D Conjugate Gradient reconstruction, then applying several regularization steps :

Replacing all negative values by zero
Averaging along time where no movement is expected
Applying total variation denoising in space
Applying wavelets denoising in space
Applying total variation denoising in time
Applying gradient's L0 norm denoising in time
Applying total nuclear variation denoising and starting over as many times as the number of main loop iterations desired.

If both the displacement vector fields to a reference phase and from a reference phase are provided, 4D ROOSTER performs the denoising in time the following way:

each 3D volume of the sequence is warped to the reference phase using the first DVF
denoising in time is applied on the warped sequence
the difference sequence between the warped-then-denoised sequence and the warped sequence is computed
that difference sequence is warped from the reference phase using the second DVF, and added to the output of spatial denoising

If only the displacement vector field to a reference phase is provided, 4D ROOSTER performs total variation denoising in time the following way:

each 3D volume of the sequence is warped to the reference phase using the first DVF
denoising in time is applied on the warped sequence
the warped-then-denoised sequence is warped from the reference phase by an iterative procedure based on conjugate gradient. This significantly increases computation time.

Test:: rtkfourdroostertest.cxx

Author: Cyril Mory

Definition at line 207 of file rtkFourDROOSTERConeBeamReconstructionFilter.h.