Speckle (interference)

Speckle, speckle pattern, or speckle noise designates the granular structure observed in coherent light, resulting from random interference. Speckle patterns are used in a wide range of metrology techniques, as they generally allow high sensitivity and simple setups. They can also be a limiting factor in imaging systems, such as radar, synthetic aperture radar (SAR), medical ultrasound and optical coherence tomography.^[1][2][3][4] Speckle is not external noise; rather, it is an inherent fluctuation in diffuse reflections, because the scatterers are not identical for each cell, and the coherent illumination wave is highly sensitive to small variations in phase changes.^[5]

Speckle patterns arise when coherent light is randomised. The simplest case of such randomisation is when light reflects off an optically rough surface. Optically rough means that the surface profile contains fluctuations larger than the wavelength. Most common surfaces are rough to visible light, such as paper, wood, or paint.

The vast majority of surfaces, synthetic or natural, are extremely rough on the scale of the wavelength. We see the origin of this phenomenon if we model our reflectivity function as an array of scatterers. Because of the finite resolution, at any time we are receiving from a distribution of scatterers within the resolution cell. These scattered signals add coherently; that is, they add constructively and destructively depending on the relative phases of each scattered waveform. Speckle results from these patterns of constructive and destructive interference shown as bright and dark dots in the image.^[6]

Speckle in conventional radar increases the mean grey level of a local area.^[7] Speckle in SAR is generally serious, causing difficulties for image interpretation.^[7][8] It is caused by coherent processing of backscattered signals from multiple distributed targets. In SAR oceanography, for example, speckle is caused by signals from elementary scatterers, the gravity-capillary ripples, and manifests as a pedestal image, beneath the image of the sea waves.^[9][10]

The speckle can also represent some useful information, particularly when it is linked to the laser speckle and to the dynamic speckle phenomenon, where the changes of the spatial speckle pattern over time can be used as a measurement of the surface's activity, such as which is useful for measuring displacement fields via digital image correlation.

Although scientists have investigated this phenomenon since the time of Newton,^[11] speckles have come into prominence since the invention of the laser.^[12]

The speckle effect is a result of the interference of many waves of the same frequency, having different phases and amplitudes, which add together to give a resultant wave whose amplitude, and therefore intensity, varies randomly. If we model each wave by a vector, we can then see that if we add a number of vectors with random angles together, the length of the resulting vector can be anything from zero to the sum of the individual vector lengths—a 2-dimensional random walk, sometimes known as a drunkard's walk. In the limit of many interfering waves, and for polarised waves, the distribution of intensities (which go as the square of the vector's length) becomes exponential ${\textstyle P(I)={\frac {1}{\langle I\rangle }}\exp \left({\frac {-I}{\langle I\rangle }}\right)}$, where ${\displaystyle \langle I\rangle }$ is the mean intensity.^{[1][2][13][14]}

When a surface is illuminated by a light wave, according to diffraction theory, each point on an illuminated surface acts as a source of secondary spherical waves. The light at any point in the scattered light field is made up of waves which have been scattered from each point on the illuminated surface. If the surface is rough enough to create path-length differences exceeding one wavelength, giving rise to phase changes greater than 2π, the amplitude, and hence the intensity, of the resultant light varies randomly.

If light of low coherence (i.e., made up of many wavelengths) is used, a speckle pattern will not normally be observed, because the speckle patterns produced by individual wavelengths have different dimensions and will normally average one another out. However, we can observe speckle patterns in polychromatic light in some conditions, such as with partially coherent light,^[15] or spatially coherent light.^[16]

When a rough surface which is illuminated by a coherent light (e.g. a laser beam) is imaged, a speckle pattern is observed in the image plane; this is called a "subjective speckle pattern" – see image above. It is called "subjective" because the detailed structure of the speckle pattern depends on the viewing system parameters; for instance, if the size of the lens aperture changes, the size of the speckles change. If the position of the imaging system is altered, the pattern will gradually change and will eventually be unrelated to the original speckle pattern.

We can explain this as follows. We can consider each point in the image to be illuminated by a finite area in the object.^{[clarification needed]} We determine the size of this area by the diffraction-limited resolution of the lens which is given by the Airy disk whose diameter is 2.4λu/D, where λ is the wavelength of the light, u is the distance between the object and the lens, and D is the diameter of the lens aperture. (This is a simplified model of diffraction-limited imaging.)

The light at neighboring points in the image has been scattered from areas which have many points in common and the intensity of two such points will not differ much. However, two points in the image which are illuminated by areas in the object which are separated by the diameter of the Airy disk, have light intensities which are unrelated. This corresponds to a distance in the image of 2.4λv/D where v is the distance between the lens and the image. Thus, the "size" of the speckles in the image is of this order.

We can observe the change in speckle size with lens aperture by looking at a laser spot on a wall directly, and then through a very small hole. The speckles will be seen to increase significantly in size. Also, the speckle pattern itself will change when moving the position of the eye while keeping the laser pointer steady. A further proof that the speckle pattern is formed only in the image plane (in the specific case the eye's retina) is that the speckles will stay visible if the eye's focus is shifted away from the wall (this is different for an objective speckle pattern, where the speckle visibility is lost under defocusing).

When laser light which has been scattered off a rough surface falls on another surface, it forms an "objective speckle pattern". If a photographic plate or another 2-D optical sensor is located within the scattered light field without a lens, a speckle pattern is obtained whose characteristics depend on the geometry of the system and the wavelength of the laser. The speckle pattern in the figure was obtained by pointing a laser beam at the surface of a mobile phone so that the scattered light fell onto an adjacent wall. A photograph was then taken of the speckle pattern formed on the wall. Strictly speaking, this also has a second subjective speckle pattern but its dimensions are much smaller than the objective pattern so it is not seen in the image.

Contributions from the whole of the scattering surface make up the light at a given point in the speckle pattern. The relative phases of these scattered waves vary across the scattering surface, so that the resulting phase on each point of the second surface varies randomly. The pattern is the same regardless of how it is imaged, just as if it were a painted pattern.

The "size" of the speckles is a function of the wavelength of the light, the size of the laser beam which illuminates the first surface, and the distance between this surface and the surface where the speckle pattern is formed. This is the case because when the angle of scattering changes such that the relative path difference between light scattered from the centre of the illuminated area compared with light scattered from the edge of the illuminated area changes by λ, the intensity becomes uncorrelated. Dainty^[1] derives an expression for the mean speckle size as λz/L where L is the width of the illuminated area and z is the distance between the object and the location of the speckle pattern.

Objective speckles are usually obtained in the far field (also called Fraunhofer region, that is the zone where Fraunhofer diffraction happens). This means that they are generated "far" from the object that emits or scatters light. We can also observe speckles close to the scattering object, in the near field (also called Fresnel region, that is, the region where Fresnel diffraction happens). We call this kind of speckles near-field speckles. See near and far field for a more rigorous definition of "near" and "far".

The statistical properties of a far-field speckle pattern (i.e., the speckle form and dimension) depend on the form and dimension of the region hit by laser light. By contrast, a very interesting feature of near field speckles is that their statistical properties are closely related to the form and structure of the scattering object: objects that scatter at high angles generate small near field speckles, and vice versa. Under Rayleigh–Gans condition, in particular, speckle dimension mirrors the average dimension of the scattering objects, while, in general, the statistical properties of near field speckles generated by a sample depend on the light scattering distribution.^[17][18]

Actually, the condition under which the near field speckles appear has been described as more strict than the usual Fresnel condition.^[19]

Speckle patterns have been used in a variety of applications in microscopy,^[20][21] imaging,^[22][23] and optical manipulation.^[24][25][26]

When lasers were first invented, the speckle effect was considered to be a severe drawback in using lasers to illuminate objects, particularly in holographic imaging because of the grainy image produced. Researchers later realized that speckle patterns could carry information about the object's surface deformations, and exploited this effect in holographic interferometry and electronic speckle pattern interferometry.^[27] Speckle imaging and eye testing using speckle also use the speckle effect.

Speckle is the chief limitation of coherent lidar and coherent imaging in optical heterodyne detection.

In the case of near field speckles, the statistical properties depend on the light scattering distribution of a given sample. This allows the use of near field speckle analysis to detect the scattering distribution; this is the so-called near-field scattering technique.^[28]

When the speckle pattern changes in time, due to changes in the illuminated surface, the phenomenon is known as dynamic speckle, and it can be used to measure activity, by means of, for example, an optical flow sensor (optical computer mouse). In biological materials, the phenomenon is known as biospeckle.

In a static environment, changes in speckle can also be used as a sensitive probe of the light source. This can be used in a wavemeter configuration, with a resolution around 1 attometre,^[29] (equivalent to 1 part in 10¹² of the wavelength, equivalent to measuring the length of a football field at the resolution of a single atom^[30]) and can also stabilise the wavelength of lasers^[31] or measure polarization.^[32]

The disordered pattern produced by speckle has been used in quantum simulations with cold atoms. The randomly-distributed regions of bright and dark light act as an analog of disorder in solid-state systems, and are used to investigate localization phenomena.^[33]

In fluorescence microscopy, a sub-diffraction-limited resolution can be obtained in 2D from saturable/photoconvertible pattern illumination techniques like stimulated emission depletion (STED) microscopy, ground state depletion (GSD) microscopy, and reversible saturable optical fluorescence transitions (RESOLFT). Adapting speckle patterns for use in these applications enables parallel 3D super-resolution imaging.^[34]

Speckle is considered to be a problem in laser based display systems like the Laser TV. Speckle is usually quantified by the speckle contrast. Speckle contrast reduction is essentially the creation of many independent speckle patterns, so that they average out on the retina/detector. This can be achieved by,^[35]

• Angle diversity: illumination from different angles
• Polarization diversity: use of different polarization states
• Wavelength diversity: use of laser sources which differ in wavelength by a small amount

Rotating diffusers—which destroys the spatial coherence of the laser light—can also be used to reduce the speckle. Moving/vibrating screens or fibers may also be solutions.^[36] The Mitsubishi Laser TV appears to use such a screen which requires special care according to their product manual. A more detailed discussion on laser speckle reduction can be found here.^[37]

Synthetic array heterodyne detection was developed to reduce speckle noise in coherent optical imaging and coherent differential absorption LIDAR.

In scientific applications, a spatial filter can be used to reduce speckle.

Several different methods are used to eliminate speckle, based upon different mathematical models of the phenomenon.^[9] One method, for example, employs multiple-look processing (a.k.a. multi-look processing), averaging out the speckle by taking several "looks" at a target in a single radar sweep.^[7][8] The average is the incoherent average of the looks.^[8]

A second method involves using adaptive and non-adaptive filters on the signal processing (where adaptive filters adapt their weightings across the image to the speckle level, and non-adaptive filters apply the same weightings uniformly across the entire image). Such filtering also eliminates actual image information as well, in particular high-frequency information, and the applicability of filtering and the choice of filter type involves tradeoffs. Adaptive speckle filtering is better at preserving edges and detail in high-texture areas (such as forests or urban areas). Non-adaptive filtering is simpler to implement, and requires less computational power, however.^[7][8]

There are two forms of non-adaptive speckle filtering: one based on the mean and one based upon the median (within a given rectangular area of pixels in the image). The latter is better at preserving edges whilst eliminating spikes, than the former is. There are many forms of adaptive speckle filtering,^[38] including the Lee filter, the Frost filter, and the refined gamma maximum-A-posteriori (RGMAP) filter. They all rely upon three fundamental assumptions in their mathematical models, however:^[7]

• Speckle in SAR is a multiplicative, i.e. it is in direct proportion to the local grey level in any area.^[7]
• The signal and the speckle are statistically independent of each other.^[7]
• The sample mean and variance of a single pixel are equal to the mean and variance of the local area that is centred on that pixel.^[7]

The Lee filter converts the multiplicative model into an additive one, thereby reducing the problem of dealing with speckle to a known tractable case.^[39]

Recently, the use of wavelet transform has led to significant advances in image analysis. The main reason for the use of multiscale processing is the fact that many natural signals, when decomposed into wavelet bases are significantly simplified and can be modeled by known distributions. Besides, wavelet decomposition is able to separate signals at different scales and orientations. Therefore, the original signal at any scale and direction can be recovered and useful details are not lost.^[40]

The first multiscale speckle reduction methods were based on the thresholding of detail subband coefficients.^[41] Wavelet thresholding methods have some drawbacks: (i) the choice of threshold is made in an ad hoc manner, supposing that wanted and unwanted components of the signal obey their known distributions, irrespective of their scale and orientations; and (ii) the thresholding procedure generally results in some artifacts in the denoised image. To address these disadvantages, non-linear estimators, based on Bayes' theory were developed.^[40][42]

Speckle patterns can also be observed over time instead of space. This is the case of phase sensitive optical time-domain reflectometry, where multiple reflections of a coherent pulse generated at different instants interfere to produce a pseudorandom time-domain signal.^[43]

Speckle interference pattern may be decomposed in the sum of plane waves. There exist a set of points where amplitude of electromagnetic field is exactly zero. Researchers had recognized these points as dislocations of wave trains.^[44] We know these phase dislocations of electromagnetic fields as optical vortices.

There is a circular energy flow around each vortex core. Thus each vortex in the speckle pattern carries optical angular momentum. The angular momentum density is given by:^[45]

${\displaystyle {\begin{aligned}{\vec {\mathbf {L} }}\left({\vec {\mathbf {r} }},t\right)&={\vec {\mathbf {r} }}\times {\vec {\mathbf {S} }}\left({\vec {\mathbf {r} }},t\right)\\{\vec {\mathbf {S} }}\left({\vec {\mathbf {r} }},t\right)&=\epsilon _{0}c^{2}{\vec {\mathbf {E} }}\left({\vec {\mathbf {r} }},t\right)\times {\vec {\mathbf {B} }}\left({\vec {\mathbf {r} }},t\right).\end{aligned}}}$

Typically vortices appear in speckle pattern in pairs. These vortex - antivortex pairs are placed randomly in space. One may show that electromagnetic angular momentum of each vortex pair is close to zero.^[46] In phase conjugating mirrors based on stimulated Brillouin scattering optical vortices excite acoustical vortices.^[47]

Apart from formal decomposition in Fourier series the speckle pattern may be composed for plane waves emitted by tilted regions of the phase plate. This approach significantly simplifies numerical modelling. 3D numerical emulation demonstrates the intertwining of vortices which leads to formation of ropes in optical speckle.^[48]

• Diffusing-wave spectroscopy
• Gaussian noise
• Laser speckle contrast imaging
• Salt and pepper noise

• Cheng Hua & Tian Jinwen (2009). "Speckle Reduction of Synthetic Aperture Radar Images Based on Fuzzy Logic". First International Workshop on Education Technology and Computer Science, Wuhan, Hubei, China, March 07–08 2009. Vol. 1. pp. 933–937. doi:10.1109/ETCS.2009.212.
• Forouzanfar, M., Abrishami-Moghaddam, H., and Dehghani, M., (2007) "Speckle reduction in medical ultrasound images using a new multiscale bivariate Bayesian MMSE-based method," IEEE 15th Signal Processing and Communication Applications Conf. (SIU'07), Turkey, June 2007, pp. 1–4.
• Sedef Kent; Osman Nuri Oçan & Tolga Ensari (2004). "Speckle Reduction of Synthetic Aperture Radar Images Using Wavelet Filtering". In ITG; VDE; FGAN; DLR; EADS & astrium (eds.). EUSAR 2004 — Proceedings — 5th European Conference on Synthetic Aperture Radar, May 25–27, 2004, Ulm, Germany. Margret Schneider. pp. 1001–1003. ISBN 9783800728282.
• Andrew K. Chan & Cheng Peng (2003). "Wavelet applications to the processing of SAR images". Wavelets for sensing technologies. Artech House remote sensing library. Artech House. ISBN 9781580533171.
• Jong-Sen Lee & Eric Pottier (2009). "Polarimetric SAR speckle filtering". Polarimetric Radar Imaging: From Basics to Applications. Optical science and engineering series. Vol. 142. CRC Press. ISBN 9781420054972.

• Seeing speckle in your fingernail
• Research group on light scattering and photonic materials
• Brogioli, Doriano; Vailati, Alberto; Giglio, Marzio (2009). "Near Field Speckles". arXiv:0907.3376 [physics.optics].