#include <itkLandweberDeconvolutionImageFilter.h>

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

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

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 >	Superclass
Public Member Functions
virtual ::itk::LightObject::Pointer	CreateAnother (void) const
virtual const char *	GetNameOfClass () const

virtual void	SetAlpha (double _arg)
virtual double	GetAlpha ()
Static Public Member Functions
static Pointer	New ()
Protected Types
typedef Superclass::FFTFilterType	FFTFilterType
typedef Superclass::IFFTFilterType	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 ()
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

Detailed Description

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

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:: Gaeten 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.