ITKShotNoiseImage

=================

Group (Subgroup)

ITKImageProcessing (ITKImageProcessing)

Description

Alter an image with shot noise.

The shot noise follows a Poisson distribution:

\par \f$ I = N(I_0) \f$

\par where \f$ N(I_0) \f$ is a Poisson-distributed random variable of mean \f$ I_0 \f$ . The noise is thus dependent on the pixel intensities in the image.

The intensities in the image can be scaled by a user provided value to map pixel values to the actual number of particles. The scaling can be seen as the inverse of the gain used during the acquisition. The noisy signal is then scaled back to its input intensity range:

\par \f$ I = \frac{N(I_0 \times s)}{s} \f$

\par where \f$ s \f$ is the scale factor.

The Poisson-distributed variable \f$ \lambda \f$ is computed by using the algorithm:

\par \f$ \begin{array}{l} k \leftarrow 0 \ p \leftarrow 1 \ \textbf{repeat} \ \left{ \begin{array}{l} k \leftarrow k+1 \ p \leftarrow p \ast U() \end{array} \right. \ \textbf{until } p > e^{\lambda} \ \textbf{return} (k) \end{array} \f$

\par where \f$ U() \f$ provides a uniformly distributed random variable in the interval \f$ [0,1] \f$ .

This algorithm is very inefficient for large values of \f$ \lambda \f$ , though. Fortunately, the Poisson distribution can be accurately approximated by a Gaussian distribution of mean and variance \f$ \lambda \f$ when \f$ \lambda \f$ is large enough. In this implementation, this value is considered to be 50. This leads to the faster algorithm:

\par \f$ \lambda + \sqrt{\lambda} \times N()\f$

\par where \f$ N() \f$ is a normally distributed random variable of mean 0 and variance 1.

\author Gaetan Lehmann

This code was contributed in the Insight Journal paper "Noise Simulation". https://hdl.handle.net/10380/3158

Parameters

Name Type Description
Scale double Set/Get the value to map the pixel value to the actual particle counting. The scaling can be seen as the inverse of the gain used during the acquisition. The noisy signal is then scaled back to its input intensity range. Defaults to 1.0.
Seed double N/A

Required Geometry

Image

Required Objects

Kind Default Name Type Component Dimensions Description
Cell Attribute Array None N/A (1) Array containing input image

Created Objects

Kind Default Name Type Component Dimensions Description
Cell Attribute Array None (1) Array containing filtered image

References

[1] T.S. Yoo, M. J. Ackerman, W. E. Lorensen, W. Schroeder, V. Chalana, S. Aylward, D. Metaxas, R. Whitaker. Engineering and Algorithm Design for an Image Processing API: A Technical Report on ITK - The Insight Toolkit. In Proc. of Medicine Meets Virtual Reality, J. Westwood, ed., IOS Press Amsterdam pp 586-592 (2002). [2] H. Johnson, M. McCormick, L. Ibanez. The ITK Software Guide: Design and Functionality. Fourth Edition. Published by Kitware Inc. 2015 ISBN: 9781-930934-28-3 [3] H. Johnson, M. McCormick, L. Ibanez. The ITK Software Guide: Introduction and Development Guidelines. Fourth Edition. Published by Kitware Inc. 2015 ISBN: 9781-930934-27-6

Example Pipelines

Please see the description file distributed with this plugin.

DREAM3D Mailing Lists

If you need more help with a filter, please consider asking your question on the DREAM3D Users mailing list: https://groups.google.com/forum/?hl=en#!forum/dream3d-users