ITK
5.0.0
Insight Segmentation and Registration Toolkit
|
#include <itkLandweberDeconvolutionImageFilter.h>
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 where is the estimate of the unblurred image, is the convolution operator, is the blurring kernel, and 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 https://hdl.handle.net/10380/3207
Definition at line 101 of file itkLandweberDeconvolutionImageFilter.h.
Public Types | |
using | ConstPointer = SmartPointer< const Self > |
using | InputImageType = TInputImage |
using | InternalComplexImagePointerType = typename Superclass::InternalComplexImagePointerType |
using | InternalComplexImageType = typename Superclass::InternalComplexImageType |
using | InternalComplexType = typename Superclass::InternalComplexType |
using | InternalImagePointerType = typename Superclass::InternalImagePointerType |
using | InternalImageType = typename Superclass::InternalImageType |
using | KernelImageType = TKernelImage |
using | OutputImageType = TOutputImage |
using | Pointer = SmartPointer< Self > |
using | Self = LandweberDeconvolutionImageFilter |
using | Superclass = IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > |
Public Types inherited from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > | |
using | ConstPointer = SmartPointer< const Self > |
using | InputImageType = TInputImage |
using | InternalComplexImagePointerType = typename Superclass::InternalComplexImagePointerType |
using | InternalComplexImageType = typename Superclass::InternalComplexImageType |
using | InternalComplexType = typename Superclass::InternalComplexType |
using | InternalImagePointerType = typename Superclass::InternalImagePointerType |
using | InternalImageType = typename Superclass::InternalImageType |
using | KernelImageType = TKernelImage |
using | OutputImageType = TOutputImage |
using | Pointer = SmartPointer< Self > |
using | Self = IterativeDeconvolutionImageFilter |
using | Superclass = FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > |
Public Types inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > | |
using | BoundaryConditionPointerType = typename Superclass::BoundaryConditionPointerType |
using | BoundaryConditionType = typename Superclass::BoundaryConditionType |
using | ConstPointer = SmartPointer< const Self > |
using | InputImageType = TInputImage |
using | InputIndexType = typename InputImageType::IndexType |
using | InputPixelType = typename InputImageType::PixelType |
using | InputRegionType = typename InputImageType::RegionType |
using | InputSizeType = typename InputImageType::SizeType |
using | InternalComplexImagePointerType = typename InternalComplexImageType::Pointer |
using | InternalComplexImageType = Image< InternalComplexType, TInputImage::ImageDimension > |
using | InternalComplexType = std::complex< TInternalPrecision > |
using | InternalImagePointerType = typename InternalImageType::Pointer |
using | InternalImageType = Image< TInternalPrecision, TInputImage::ImageDimension > |
using | KernelImageType = TKernelImage |
using | KernelIndexType = typename KernelImageType::IndexType |
using | KernelPixelType = typename KernelImageType::PixelType |
using | KernelRegionType = typename KernelImageType::RegionType |
using | KernelSizeType = typename KernelImageType::SizeType |
using | OutputImageType = TOutputImage |
using | OutputIndexType = typename OutputImageType::IndexType |
using | OutputPixelType = typename OutputImageType::PixelType |
using | OutputRegionType = typename OutputImageType::RegionType |
using | OutputSizeType = typename OutputImageType::SizeType |
using | Pointer = SmartPointer< Self > |
using | Self = FFTConvolutionImageFilter |
using | SizeValueType = typename InputSizeType::SizeValueType |
using | Superclass = ConvolutionImageFilterBase< TInputImage, TKernelImage, TOutputImage > |
Public Types inherited from itk::ConvolutionImageFilterBase< TInputImage, TKernelImage, TOutputImage > | |
using | BoundaryConditionPointerType = BoundaryConditionType * |
using | BoundaryConditionType = ImageBoundaryCondition< TInputImage > |
using | ConstPointer = SmartPointer< const Self > |
using | DefaultBoundaryConditionType = ZeroFluxNeumannBoundaryCondition< TInputImage > |
using | InputImageType = TInputImage |
using | InputIndexType = typename InputImageType::IndexType |
using | InputPixelType = typename InputImageType::PixelType |
using | InputRegionType = typename InputImageType::RegionType |
using | InputSizeType = typename InputImageType::SizeType |
using | KernelImageType = TKernelImage |
using | KernelIndexType = typename KernelImageType::IndexType |
using | KernelPixelType = typename KernelImageType::PixelType |
using | KernelRegionType = typename KernelImageType::RegionType |
using | KernelSizeType = typename KernelImageType::SizeType |
using | OutputImageType = TOutputImage |
using | OutputIndexType = typename OutputImageType::IndexType |
using | OutputPixelType = typename OutputImageType::PixelType |
enum | OutputRegionModeType { SAME = 0, VALID } |
using | OutputRegionType = typename OutputImageType::RegionType |
using | OutputSizeType = typename OutputImageType::SizeType |
using | Pointer = SmartPointer< Self > |
using | Self = ConvolutionImageFilterBase |
using | SizeValueType = typename InputSizeType::SizeValueType |
using | Superclass = ImageToImageFilter< TInputImage, TOutputImage > |
Public Types inherited from itk::ImageToImageFilter< TInputImage, TOutputImage > | |
using | ConstPointer = SmartPointer< const Self > |
using | InputImageConstPointer = typename InputImageType::ConstPointer |
using | InputImagePixelType = typename InputImageType::PixelType |
using | InputImagePointer = typename InputImageType::Pointer |
using | InputImageRegionType = typename InputImageType::RegionType |
using | InputImageType = TInputImage |
using | OutputImagePixelType = typename Superclass::OutputImagePixelType |
using | OutputImageRegionType = typename Superclass::OutputImageRegionType |
using | Pointer = SmartPointer< Self > |
using | Self = ImageToImageFilter |
using | Superclass = ImageSource< TOutputImage > |
Public Types inherited from itk::ImageSource< TOutputImage > | |
using | ConstPointer = SmartPointer< const Self > |
using | DataObjectIdentifierType = Superclass::DataObjectIdentifierType |
using | DataObjectPointer = DataObject::Pointer |
using | DataObjectPointerArraySizeType = Superclass::DataObjectPointerArraySizeType |
using | OutputImagePixelType = typename OutputImageType::PixelType |
using | OutputImagePointer = typename OutputImageType::Pointer |
using | OutputImageRegionType = typename OutputImageType::RegionType |
using | OutputImageType = TOutputImage |
using | Pointer = SmartPointer< Self > |
using | Self = ImageSource |
using | Superclass = ProcessObject |
Public Types inherited from itk::ProcessObject | |
using | ConstPointer = SmartPointer< const Self > |
using | DataObjectIdentifierType = DataObject::DataObjectIdentifierType |
using | DataObjectPointer = DataObject::Pointer |
using | DataObjectPointerArray = std::vector< DataObjectPointer > |
using | DataObjectPointerArraySizeType = DataObjectPointerArray::size_type |
using | MultiThreaderType = MultiThreaderBase |
using | NameArray = std::vector< DataObjectIdentifierType > |
using | Pointer = SmartPointer< Self > |
using | Self = ProcessObject |
using | Superclass = Object |
Public Types inherited from itk::Object | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = Object |
using | Superclass = LightObject |
Public Types inherited from itk::LightObject | |
using | ConstPointer = SmartPointer< const Self > |
using | Pointer = SmartPointer< Self > |
using | Self = LightObject |
Static Public Member Functions | |
static Pointer | New () |
Static Public Member Functions inherited from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > | |
static Pointer | New () |
Static Public Member Functions inherited from itk::Object | |
static bool | GetGlobalWarningDisplay () |
static void | GlobalWarningDisplayOff () |
static void | GlobalWarningDisplayOn () |
static Pointer | New () |
static void | SetGlobalWarningDisplay (bool flag) |
Static Public Member Functions inherited from itk::LightObject | |
static void | BreakOnError () |
static Pointer | New () |
Private Attributes | |
double | m_Alpha |
IFFTFilterType::Pointer | m_IFFTFilter |
LandweberFilterType::Pointer | m_LandweberFilter |
InternalComplexImagePointerType | m_TransformedInput |
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::ConstPointer = SmartPointer< const Self > |
Definition at line 114 of file itkLandweberDeconvolutionImageFilter.h.
|
protected |
Definition at line 153 of file itkLandweberDeconvolutionImageFilter.h.
|
protected |
Definition at line 154 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::InputImageType = TInputImage |
Other useful type alias.
Definition at line 117 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::InternalComplexImagePointerType = typename Superclass::InternalComplexImagePointerType |
Definition at line 126 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::InternalComplexImageType = typename Superclass::InternalComplexImageType |
Definition at line 125 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::InternalComplexType = typename Superclass::InternalComplexType |
Definition at line 124 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::InternalImagePointerType = typename Superclass::InternalImagePointerType |
Definition at line 123 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::InternalImageType = typename Superclass::InternalImageType |
Internal types used by the FFT filters.
Definition at line 122 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::KernelImageType = TKernelImage |
Definition at line 118 of file itkLandweberDeconvolutionImageFilter.h.
|
private |
Definition at line 171 of file itkLandweberDeconvolutionImageFilter.h.
|
private |
Definition at line 166 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::OutputImageType = TOutputImage |
Definition at line 119 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::Pointer = SmartPointer< Self > |
Definition at line 113 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::Self = LandweberDeconvolutionImageFilter |
Standard type alias.
Definition at line 108 of file itkLandweberDeconvolutionImageFilter.h.
using itk::LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >::Superclass = IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > |
Definition at line 112 of file itkLandweberDeconvolutionImageFilter.h.
|
protected |
|
overrideprotected |
|
virtual |
Create an object from an instance, potentially deferring to a factory. This method allows you to create an instance of an object that is exactly the same type as the referring object. This is useful in cases where an object has been cast back to a base class.
Reimplemented from itk::FFTConvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >.
Reimplemented in itk::ProjectedLandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >, and itk::ProjectedIterativeDeconvolutionImageFilter< LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > >.
|
overrideprotectedvirtual |
Runs after iterating.
Reimplemented from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >.
|
virtual |
Set/get relaxation factor.
|
virtual |
Runtime information support.
Reimplemented from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >.
Reimplemented in itk::ProjectedLandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >, and itk::ProjectedIterativeDeconvolutionImageFilter< LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > >.
|
overrideprotectedvirtual |
Runs before iterating .
Reimplemented from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >.
Reimplemented in itk::ProjectedIterativeDeconvolutionImageFilter< LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > >.
|
overrideprotectedvirtual |
Performs an iteration.
Implements itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >.
Reimplemented in itk::ProjectedIterativeDeconvolutionImageFilter< LandweberDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision > >.
|
static |
Method for creation through the object factory.
|
overrideprotectedvirtual |
Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes.
Reimplemented from itk::IterativeDeconvolutionImageFilter< TInputImage, TKernelImage, TOutputImage, TInternalPrecision >.
|
virtual |
Set/get relaxation factor.
|
private |
Definition at line 159 of file itkLandweberDeconvolutionImageFilter.h.
|
private |
Definition at line 174 of file itkLandweberDeconvolutionImageFilter.h.
|
private |
Definition at line 173 of file itkLandweberDeconvolutionImageFilter.h.
|
private |
Definition at line 161 of file itkLandweberDeconvolutionImageFilter.h.