#include <itkLandweberDeconvolutionImageFilter.h>

Inheritance diagram for itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >:

Collaboration diagram for itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >:

Detailed Description

template<class TInputImage, class TKernelImage = TInputImage, class TOutputImage = TInputImage, class TInternalPrecision = double>
class itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >

Deconvolve an image using the Landweber deconvolution algorithm.

This filter implements the Landweber deconvolution algorthm as defined in Bertero M and Boccacci P, "Introduction to Inverse Problems in Imaging", 1998. The algorithm assumes that the input image has been formed by a linear shift-invariant system with a known kernel.

The Landweber algorithm converges to a solution that minimizes the sum of squared errors $||f \otimes h - g||$ where $f$ is the estimate of the unblurred image, $\otimes$ is the convolution operator, $h$ is the blurring kernel, and $g$ is the blurred input image. As such, it is best suited for images that have zero-mean Gaussian white noise.

This is the base implementation of the Landweber algorithm. It may produce results with negative values. For a version of this algorithm that enforces a positivity constraint on each intermediate solution, see ProjectedLandweberDeconvolutionImageFilter.

This code was adapted from the Insight Journal contribution:

"Deconvolution: infrastructure and reference algorithms" by Gaetan Lehmann http://hdl.handle.net/10380/3207

Author: Gaetan Lehmann, Biologie du Developpement et de la Reproduction, INRA de Jouy-en-Josas, France; Cory Quammen, The University of North Carolina at Chapel Hill

See Also: IterativeDeconvolutionImageFilter; RichardsonLucyDeconvolutionImageFilter; ProjectedLandweberDeconvolutionImageFilter

Definition at line 101 of file itkLandweberDeconvolutionImageFilter.h.

Public Types
typedef SmartPointer< const Self >	ConstPointer

typedef TInputImage	InputImageType

typedef Superclass::InternalComplexImagePointerType	InternalComplexImagePointerType

typedef Superclass::InternalComplexImageType	InternalComplexImageType

typedef Superclass::InternalComplexType	InternalComplexType

typedef Superclass::InternalImagePointerType	InternalImagePointerType

typedef Superclass::InternalImageType	InternalImageType

typedef TKernelImage	KernelImageType

typedef TOutputImage	OutputImageType

typedef SmartPointer< Self >	Pointer

typedef LandweberDeconvolutionImageFilter	Self

typedef IterativeDeconvolutionImageFilter < TInputImage, TKernelImage, TOutputImage, TInternalPrecision >	Superclass

Public Types inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
typedef SmartPointer< const Self >	ConstPointer

typedef TInputImage	InputImageType

typedef Superclass::InternalComplexImagePointerType	InternalComplexImagePointerType

typedef Superclass::InternalComplexImageType	InternalComplexImageType

typedef Superclass::InternalComplexType	InternalComplexType

typedef Superclass::InternalImagePointerType	InternalImagePointerType

typedef Superclass::InternalImageType	InternalImageType

typedef TKernelImage	KernelImageType

typedef TOutputImage	OutputImageType

typedef SmartPointer< Self >	Pointer

typedef IterativeDeconvolutionImageFilter	Self

typedef FFTConvolutionImageFilter < TInputImage, TKernelImage, TOutputImage, TInternalPrecision >	Superclass

Public Types inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
typedef Superclass::BoundaryConditionPointerType	BoundaryConditionPointerType

typedef Superclass::BoundaryConditionType	BoundaryConditionType

typedef SmartPointer< const Self >	ConstPointer

typedef TInputImage	InputImageType

typedef InputImageType::IndexType	InputIndexType

typedef InputImageType::PixelType	InputPixelType

typedef InputImageType::RegionType	InputRegionType

typedef InputImageType::SizeType	InputSizeType

typedef InternalComplexImageType::Pointer	InternalComplexImagePointerType

typedef Image < InternalComplexType, TInputImage::ImageDimension >	InternalComplexImageType

typedef std::complex < TInternalPrecision >	InternalComplexType

typedef InternalImageType::Pointer	InternalImagePointerType

typedef Image < TInternalPrecision, TInputImage::ImageDimension >	InternalImageType

typedef TKernelImage	KernelImageType

typedef KernelImageType::IndexType	KernelIndexType

typedef KernelImageType::PixelType	KernelPixelType

typedef KernelImageType::RegionType	KernelRegionType

typedef KernelImageType::SizeType	KernelSizeType

typedef TOutputImage	OutputImageType

typedef OutputImageType::IndexType	OutputIndexType

typedef OutputImageType::PixelType	OutputPixelType

typedef OutputImageType::RegionType	OutputRegionType

typedef OutputImageType::SizeType	OutputSizeType

typedef SmartPointer< Self >	Pointer

typedef FFTConvolutionImageFilter	Self

typedef InputSizeType::SizeValueType	SizeValueType

typedef ConvolutionImageFilterBase < TInputImage, TKernelImage, TOutputImage >	Superclass

Public Types inherited from itk::ConvolutionImageFilterBase< TInputImage, TKernelImage, TOutputImage >
typedef BoundaryConditionType *	BoundaryConditionPointerType

typedef ImageBoundaryCondition < TInputImage >	BoundaryConditionType

typedef SmartPointer< const Self >	ConstPointer

typedef ZeroFluxNeumannBoundaryCondition < TInputImage >	DefaultBoundaryConditionType

typedef TInputImage	InputImageType

typedef InputImageType::IndexType	InputIndexType

typedef InputImageType::PixelType	InputPixelType

typedef InputImageType::RegionType	InputRegionType

typedef InputImageType::SizeType	InputSizeType

typedef TKernelImage	KernelImageType

typedef KernelImageType::IndexType	KernelIndexType

typedef KernelImageType::PixelType	KernelPixelType

typedef KernelImageType::RegionType	KernelRegionType

typedef KernelImageType::SizeType	KernelSizeType

typedef TOutputImage	OutputImageType

typedef OutputImageType::IndexType	OutputIndexType

typedef OutputImageType::PixelType	OutputPixelType

enum	OutputRegionModeType { SAME = 0, VALID }

typedef OutputImageType::RegionType	OutputRegionType

typedef OutputImageType::SizeType	OutputSizeType

typedef SmartPointer< Self >	Pointer

typedef ConvolutionImageFilterBase	Self

typedef InputSizeType::SizeValueType	SizeValueType

typedef ImageToImageFilter < TInputImage, TOutputImage >	Superclass

Public Types inherited from itk::ImageToImageFilter< TInputImage, TOutputImage >
typedef SmartPointer< const Self >	ConstPointer

typedef InputImageType::ConstPointer	InputImageConstPointer

typedef InputImageType::PixelType	InputImagePixelType

typedef InputImageType::Pointer	InputImagePointer

typedef InputImageType::RegionType	InputImageRegionType

typedef TInputImage	InputImageType

typedef Superclass::OutputImagePixelType	OutputImagePixelType

typedef Superclass::OutputImageRegionType	OutputImageRegionType

typedef SmartPointer< Self >	Pointer

typedef ImageToImageFilter	Self

typedef ImageSource< TOutputImage >	Superclass

Public Types inherited from itk::ImageSource< TOutputImage >
typedef SmartPointer< const Self >	ConstPointer

typedef Superclass::DataObjectIdentifierType	DataObjectIdentifierType

typedef DataObject::Pointer	DataObjectPointer

typedef Superclass::DataObjectPointerArraySizeType	DataObjectPointerArraySizeType

typedef OutputImageType::PixelType	OutputImagePixelType

typedef OutputImageType::Pointer	OutputImagePointer

typedef OutputImageType::RegionType	OutputImageRegionType

typedef TOutputImage	OutputImageType

typedef SmartPointer< Self >	Pointer

typedef ImageSource	Self

typedef ProcessObject	Superclass

Public Types inherited from itk::ProcessObject
typedef SmartPointer< const Self >	ConstPointer

typedef DataObject::DataObjectIdentifierType	DataObjectIdentifierType

typedef DataObject::Pointer	DataObjectPointer

typedef std::vector < DataObjectPointer >	DataObjectPointerArray

typedef DataObjectPointerArray::size_type	DataObjectPointerArraySizeType

typedef std::vector < DataObjectIdentifierType >	NameArray

typedef SmartPointer< Self >	Pointer

typedef ProcessObject	Self

typedef Object	Superclass

Public Types inherited from itk::Object
typedef SmartPointer< const Self >	ConstPointer

typedef SmartPointer< Self >	Pointer

typedef Object	Self

typedef LightObject	Superclass

Public Types inherited from itk::LightObject
typedef SmartPointer< const Self >	ConstPointer

typedef SmartPointer< Self >	Pointer

typedef LightObject	Self

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

virtual const char *	GetNameOfClass () const


virtual void	SetAlpha (double _arg)

virtual double	GetAlpha ()

Public Member Functions inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
virtual const InternalImageType *	GetCurrentEstimate () const

virtual unsigned int	GetIteration () const

virtual InternalImageType *	GetModifiableCurrentEstimate ()

virtual void	SetNumberOfIterations (unsigned int _arg)

virtual unsigned int	GetNumberOfIterations ()

void	SetStopIteration (bool stop)

virtual bool	GetStopIteration () const

Static Public Member Functions
static Pointer	New ()

Protected Types
typedef Superclass::FFTFilterType	FFTFilterType

typedef Superclass::IFFTFilterType	IFFTFilterType

Protected Types inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
typedef Superclass::FFTFilterType	FFTFilterType

typedef Superclass::IFFTFilterType	IFFTFilterType

Protected Types inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
typedef RealToHalfHermitianForwardFFTImageFilter < InternalImageType, InternalComplexImageType >	FFTFilterType

typedef HalfHermitianToRealInverseFFTImageFilter < InternalComplexImageType, InternalImageType >	IFFTFilterType

Protected Member Functions
virtual void	Finish (ProgressAccumulator *progress, float progressWeight)

virtual void	Initialize (ProgressAccumulator *progress, float progressWeight, float iterationProgressWeight)

virtual void	Iteration (ProgressAccumulator *progress, float iterationProgressWeight)

	LandweberDeconvolutionImageFilter ()

virtual void	PrintSelf (std::ostream &os, Indent indent) const

virtual	~LandweberDeconvolutionImageFilter ()

Protected Member Functions inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
virtual void	GenerateData ()

virtual void	GenerateInputRequestedRegion ()

	IterativeDeconvolutionImageFilter ()

virtual	~IterativeDeconvolutionImageFilter ()

Protected Member Functions inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
void	CropOutput (InternalImageType paddedOutput, ProgressAccumulator progress, float progressWeight)

	FFTConvolutionImageFilter ()

InputSizeType	GetPadLowerBound () const

InputSizeType	GetPadSize () const

bool	GetXDimensionIsOdd () const

void	PadInput (const InputImageType input, InternalImagePointerType &paddedInput, ProgressAccumulator progress, float progressWeight)

void	PrepareInput (const InputImageType input, InternalComplexImagePointerType &preparedInput, ProgressAccumulator progress, float progressWeight)

void	PrepareInputs (const InputImageType input, const KernelImageType kernel, InternalComplexImagePointerType &preparedInput, InternalComplexImagePointerType &preparedKernel, ProgressAccumulator *progress, float progressWeight)

void	PrepareKernel (const KernelImageType kernel, InternalComplexImagePointerType &preparedKernel, ProgressAccumulator progress, float progressWeight)

void	ProduceOutput (InternalComplexImageType paddedOutput, ProgressAccumulator progress, float progressWeight)

void	TransformPaddedInput (const InternalImageType paddedInput, InternalComplexImagePointerType &transformedInput, ProgressAccumulator progress, float progressWeight)

	~FFTConvolutionImageFilter ()

Protected Member Functions inherited from itk::ConvolutionImageFilterBase< TInputImage, TKernelImage, TOutputImage >
	ConvolutionImageFilterBase ()

void	GenerateOutputInformation ()

OutputRegionType	GetValidRegion () const

virtual void	VerifyInputInformation ()

	~ConvolutionImageFilterBase ()

Protected Member Functions inherited from itk::ImageToImageFilter< TInputImage, TOutputImage >
virtual void	CallCopyInputRegionToOutputRegion (OutputImageRegionType &destRegion, const InputImageRegionType &srcRegion)

virtual void	CallCopyOutputRegionToInputRegion (InputImageRegionType &destRegion, const OutputImageRegionType &srcRegion)

	ImageToImageFilter ()

	~ImageToImageFilter ()

void	PushBackInput (const DataObject *input)

void	PushFrontInput (const DataObject *input)

Protected Member Functions inherited from itk::ImageSource< TOutputImage >
virtual void	AfterThreadedGenerateData ()

virtual void	AllocateOutputs ()

virtual void	BeforeThreadedGenerateData ()

virtual const ImageRegionSplitterBase *	GetImageRegionSplitter (void) const

	ImageSource ()

virtual unsigned int	SplitRequestedRegion (unsigned int i, unsigned int num, OutputImageRegionType &splitRegion)

virtual void	ThreadedGenerateData (const OutputImageRegionType &outputRegionForThread, ThreadIdType threadId)

virtual	~ImageSource ()

Protected Member Functions inherited from itk::ProcessObject
virtual void	AddInput (DataObject *input)

virtual void	AddOutput (DataObject *output)

virtual void	CacheInputReleaseDataFlags ()

virtual void	GenerateOutputRequestedRegion (DataObject *output)

virtual const DataObjectPointerArraySizeType &	GetNumberOfRequiredOutputs ()

bool	IsIndexedInputName (const DataObjectIdentifierType &) const

bool	IsIndexedOutputName (const DataObjectIdentifierType &) const

	itkLegacyMacro (virtual void RemoveOutput(DataObject *output))

	itkLegacyMacro (void SetNumberOfOutputs(DataObjectPointerArraySizeType num))

	itkLegacyMacro (virtual void RemoveInput(DataObject *input))

	itkLegacyMacro (void SetNumberOfInputs(DataObjectPointerArraySizeType num))

DataObjectPointerArraySizeType	MakeIndexFromInputName (const DataObjectIdentifierType &name) const

DataObjectPointerArraySizeType	MakeIndexFromOutputName (const DataObjectIdentifierType &name) const

DataObjectIdentifierType	MakeNameFromInputIndex (DataObjectPointerArraySizeType idx) const

DataObjectIdentifierType	MakeNameFromOutputIndex (DataObjectPointerArraySizeType idx) const

	ProcessObject ()

virtual void	PropagateResetPipeline ()

virtual void	ReleaseInputs ()

virtual void	RemoveInput (const DataObjectIdentifierType &key)

virtual void	RemoveInput (DataObjectPointerArraySizeType)

virtual void	RemoveOutput (const DataObjectIdentifierType &key)

virtual void	RemoveOutput (DataObjectPointerArraySizeType idx)

virtual void	RestoreInputReleaseDataFlags ()

virtual void	SetInput (const DataObjectIdentifierType &key, DataObject *input)

virtual void	SetNthInput (DataObjectPointerArraySizeType num, DataObject *input)

virtual void	SetNthOutput (DataObjectPointerArraySizeType num, DataObject *output)

void	SetNumberOfIndexedInputs (DataObjectPointerArraySizeType num)

void	SetNumberOfIndexedOutputs (DataObjectPointerArraySizeType num)

virtual void	SetNumberOfRequiredOutputs (DataObjectPointerArraySizeType _arg)

virtual void	SetOutput (const DataObjectIdentifierType &key, DataObject *output)

virtual void	SetPrimaryInput (DataObject *input)

virtual void	SetPrimaryOutput (DataObject *output)

virtual void	VerifyPreconditions ()

	~ProcessObject ()

DataObject *	GetInput (const DataObjectIdentifierType &key)

const DataObject *	GetInput (const DataObjectIdentifierType &key) const

DataObject *	GetInput (DataObjectPointerArraySizeType idx)

const DataObject *	GetInput (DataObjectPointerArraySizeType idx) const

DataObject *	GetPrimaryInput ()

const DataObject *	GetPrimaryInput () const

virtual void	SetPrimaryInputName (const DataObjectIdentifierType &key)

virtual const char *	GetPrimaryInputName (void) const

virtual void	SetNumberOfRequiredInputs (DataObjectPointerArraySizeType)

virtual const DataObjectPointerArraySizeType &	GetNumberOfRequiredInputs ()

bool	AddRequiredInputName (const DataObjectIdentifierType &)

bool	AddRequiredInputName (const DataObjectIdentifierType &, DataObjectPointerArraySizeType idx)

bool	RemoveRequiredInputName (const DataObjectIdentifierType &)

bool	IsRequiredInputName (const DataObjectIdentifierType &) const

void	SetRequiredInputNames (const NameArray &)

DataObject *	GetOutput (const DataObjectIdentifierType &key)

const DataObject *	GetOutput (const DataObjectIdentifierType &key) const

virtual void	SetPrimaryOutputName (const DataObjectIdentifierType &key)

virtual const char *	GetPrimaryOutputName (void) const

DataObject *	GetOutput (DataObjectPointerArraySizeType idx)

const DataObject *	GetOutput (DataObjectPointerArraySizeType idx) const

DataObject *	GetPrimaryOutput ()

const DataObject *	GetPrimaryOutput () const

Protected Member Functions inherited from itk::Object
	Object ()

bool	PrintObservers (std::ostream &os, Indent indent) const

virtual void	SetTimeStamp (const TimeStamp &time)

virtual	~Object ()

Protected Member Functions inherited from itk::LightObject
virtual LightObject::Pointer	InternalClone () const

	LightObject ()

virtual void	PrintHeader (std::ostream &os, Indent indent) const

virtual void	PrintTrailer (std::ostream &os, Indent indent) const

virtual	~LightObject ()

Private Types
typedef TernaryFunctorImageFilter < InternalComplexImageType, InternalComplexImageType, InternalComplexImageType, InternalComplexImageType, LandweberFunctor >	LandweberFilterType

typedef Functor::LandweberMethod < InternalComplexType, InternalComplexType, InternalComplexType, InternalComplexType >	LandweberFunctor

Private Member Functions
	LandweberDeconvolutionImageFilter (const Self &)

void	operator= (const Self &)

Private Attributes
double	m_Alpha

IFFTFilterType::Pointer	m_IFFTFilter

LandweberFilterType::Pointer	m_LandweberFilter

InternalComplexImagePointerType	m_TransformedInput

Additional Inherited Members
Static Public Attributes inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
static const unsigned int	ImageDimension = TInputImage::ImageDimension

Static Protected Member Functions inherited from itk::ImageSource< TOutputImage >
static const ImageRegionSplitterBase *	GetGlobalDefaultSplitter ()

static ITK_THREAD_RETURN_TYPE	ThreaderCallback (void *arg)

Protected Attributes inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
InternalImagePointerType	m_CurrentEstimate

InternalComplexImagePointerType	m_TransferFunction