#include <itkLandweberDeconvolutionImageFilter.h>

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

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

List of all members.

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 >
Public Types inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
typedef Superclass::BoundaryConditionPointerType	BoundaryConditionPointerType
typedef Superclass::BoundaryConditionType	BoundaryConditionType
typedef InputImageType::IndexType	InputIndexType
typedef InputImageType::PixelType	InputPixelType
typedef InputImageType::RegionType	InputRegionType
typedef InputImageType::SizeType	InputSizeType
typedef KernelImageType::IndexType	KernelIndexType
typedef KernelImageType::PixelType	KernelPixelType
typedef KernelImageType::RegionType	KernelRegionType
typedef KernelImageType::SizeType	KernelSizeType
typedef OutputImageType::IndexType	OutputIndexType
typedef OutputImageType::PixelType	OutputPixelType
typedef OutputImageType::RegionType	OutputRegionType
typedef OutputImageType::SizeType	OutputSizeType
typedef InputSizeType::SizeValueType	SizeValueType
Public Types inherited from itk::ConvolutionImageFilterBase< TInputImage, TKernelImage, TOutputImage >
typedef ZeroFluxNeumannBoundaryCondition < TInputImage >	DefaultBoundaryConditionType
enum	OutputRegionModeType { SAME = 0, VALID }
Public Types inherited from itk::ImageToImageFilter< TInputImage, TOutputImage >
typedef InputImageType::ConstPointer	InputImageConstPointer
typedef InputImageType::PixelType	InputImagePixelType
typedef InputImageType::Pointer	InputImagePointer
typedef InputImageType::RegionType	InputImageRegionType
typedef Superclass::OutputImagePixelType	OutputImagePixelType
typedef Superclass::OutputImageRegionType	OutputImageRegionType
Public Types inherited from itk::ImageSource< TOutputImage >
typedef Superclass::DataObjectIdentifierType	DataObjectIdentifierType
typedef DataObject::Pointer	DataObjectPointer
typedef Superclass::DataObjectPointerArraySizeType	DataObjectPointerArraySizeType
typedef OutputImageType::Pointer	OutputImagePointer
Public Types inherited from itk::ProcessObject
typedef std::vector < DataObjectPointer >	DataObjectPointerArray
typedef std::vector < DataObjectIdentifierType >	NameArray
Public Types inherited from itk::Object
Public Types inherited from itk::LightObject

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 ()
virtual unsigned int	GetIteration () const
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 >
Protected Types inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >

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 ()
	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)
bool	AddRequiredInputName (const DataObjectIdentifierType &)
virtual void	CacheInputReleaseDataFlags ()
virtual void	GenerateOutputRequestedRegion (DataObject *output)
virtual const DataObjectPointerArraySizeType &	GetNumberOfRequiredInputs ()
virtual const DataObjectPointerArraySizeType &	GetNumberOfRequiredOutputs ()
bool	IsIndexedName (const DataObjectIdentifierType &) const
bool	IsRequiredInputName (const DataObjectIdentifierType &) const
	itkLegacyMacro (virtual void RemoveInput(DataObject *input))
	itkLegacyMacro (virtual void RemoveOutput(DataObject *output))
	itkLegacyMacro (void SetNumberOfInputs(DataObjectPointerArraySizeType num))
	itkLegacyMacro (void SetNumberOfOutputs(DataObjectPointerArraySizeType num))
DataObjectPointerArraySizeType	MakeIndexFromName (const DataObjectIdentifierType &) const
DataObjectIdentifierType	MakeNameFromIndex (DataObjectPointerArraySizeType) 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)
bool	RemoveRequiredInputName (const DataObjectIdentifierType &)
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	SetNumberOfRequiredInputs (DataObjectPointerArraySizeType)
virtual void	SetNumberOfRequiredOutputs (DataObjectPointerArraySizeType _arg)
virtual void	SetOutput (const DataObjectIdentifierType &key, DataObject *output)
virtual void	SetPrimaryInput (DataObject *input)
virtual void	SetPrimaryOutput (DataObject *output)
void	SetRequiredInputNames (const NameArray &)
virtual void	VerifyPreconditions ()
	~ProcessObject ()
DataObject *	GetInput (const DataObjectIdentifierType &key)
const DataObject *	GetInput (const DataObjectIdentifierType &key) const
DataObject *	GetInput (DataObjectPointerArraySizeType)
const DataObject *	GetInput (DataObjectPointerArraySizeType idx) const
DataObject *	GetPrimaryInput ()
const DataObject *	GetPrimaryInput () const
DataObject *	GetOutput (const DataObjectIdentifierType &key)
const DataObject *	GetOutput (const DataObjectIdentifierType &key) 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
Protected Attributes inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >
InternalImagePointerType	m_CurrentEstimate
InternalComplexImagePointerType	m_TransferFunction

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.