Speckle function

Speckle is the high-frequency noise on a radar data. The images generated by synthetic aperture radar (SAR) systems are highly subject to speckling effects due to the processing of scattered signals and interference of electromagnetic waves scattered from surfaces or objects. This Speckle function filters the speckled radar dataset and smooths out the noise while retaining the edges or sharp features in the image.

Smoothing algorithms in this function effectively remove speckles on the basis of a noise model. The speckle reduction filtering techniques used in this function are as follows:

These filters in the Speckle function preserve the edges and details while efficiently reducing noise speckle in the image. All these filters have parameters that you can adjust to optimize the results.

The main inputs for the Speckle function are the following:

Lee

The Lee filter reduces the speckle noise by applying a spatial filter to each pixel in an image, which filters the data based on local statistics calculated within a square window. The value of the center pixel is replaced by a value calculated using the neighboring pixels.

Inputs

The Lee filter is applied depending on the following three noise models that are exhibited by the input raster:

  • Additive
  • Multiplicative
  • Additive and Multiplicative

Basically, speckles have the nature of multiplicative noise.

The following inputs depend on the noise model selected:

Noise model

Parameter

Description

Default value

Additive

Noise Variance

This is the noise variance of the image in case of additive noise as well as additive and multiplicative noise.

0.25

Multiplicative (default)

Multiplicative Noise Mean

Mean value of multiplicative noise.

1

Number of Looks

Specifies the number of looks of the image. This is used to calculate the noise variance in case of multiplicative noise.

1

Additive and Multiplicative

Noise Variance

This is the noise variance of the image in case of additive noise as well as additive and multiplicative noise.

0.25

Additive Noise Mean

Mean value of additive noise.

0

Multiplicative Noise Mean

Mean value of multiplicative noise.

1

Lee filter noise models
NoteNote:

Additive Noise Mean is usually 0. Multiplicative Noise Mean is usually 1.

Algorithms

The algorithms used in the implementation of the Lee filter are as follows:

Noise model

Algorithm

Additive

Value of filtered pixel = LM + K * (PC - LM)

where

K (weight function) = LV / (LV + AV)

Multiplicative

Value of filtered pixel = LM + K * (PC - M * LM)

where

K (weight function) = M * LV / ((LM * LM * MV) + (M * M * LV))

where

MV = 1 / NLooks

Additive and Multiplicative

Value of filtered pixel = LM + K * (PC - M * LM - A)

where

K (weight function) = M * LV / ((LM * LM * MV) + (M * M * LV) + AV)

where

MV = (SD / LM)2

Lee filter noise model algorithms

where

PC—Center pixel value of window

LM—Local mean of filter window

LV—Local variance of filter window

M—Multiplicative noise mean (input parameter)

A—Additive noise mean (input parameter)

AV—Additive noise variance (input parameter)

MV—Multiplicative noise variance (input parameter)

SD—Standard deviation of filter window

NLooks—Number of looks (input parameter)

Enhanced Lee

The enhanced Lee filter is an altered version of the Lee filter reducing the speckle noise effectively by preserving image sharpness and detail. It requires a damping factor and number of looks.

The Number of Looks parameter controls image smoothing and estimates noise variance. The smaller the value is, the better the smoothing effect and filter performance. A larger value retains more image features.

The Damping Factor value defines the extent of exponential damping. A larger value leads to increased smoothing ability.

Inputs

The inputs for the enhanced Lee filter are as follows:

Parameter

Description

Default value

Number of Looks

Specifies the number of looks of the image

1

Damping Factor

Specifies the damping factor to define the extent of smoothing

1.0

Enhanced Lee filter inputs

Algorithm

The algorithm used in the implementation of the enhanced Lee filter is as follows:

Value of smoothed center pixel:

LM for CI <= CU

LM * K + PC * (1 - K) for CU < CI < Cmax

PC for CI >= Cmax

where

PC—Center pixel value of window

LM—Local mean of filter window

SD—Standard deviation in filter window

NLooks—Number of looks (input parameter)

D—Damping factor (input parameter)

CU = 1 / square root (NLooks) (Noise variation coefficient)

Cmax = square root (1 + 2 / NLooks) (Maximum noise variation coefficient)

CI = SD / LM(Image variation coefficient)

K = e(- D (Ci - CU) / (Cmax - CI))

Frost

The Frost filter reduces speckle noise and preserves important image features at the edges with an exponentially damped circularly symmetric filter that uses local statistics within individual filter windows.

Scene reflectivity is an important factor that differentiates the Frost filter from the Lee and Kuan filters, which is calculated by combining the observed image with the impulse response of the SAR system.

The Frost filter requires a damping factor. The Damping Factor value defines the extent of exponential damping. The smaller the value is, the better the smoothing ability and filter performance.

After application of the Frost filter, the denoised images show better sharpness at the edges.

Input

The input for the Frost filter is as follows:

Parameter

Description

Default value

Damping Factor

Specifies the damping factor to define the extent of smoothing

1.0

Frost filter input

Algorithm

The implementation of this filter consists of defining a circularly symmetric filter with a set of weighting values M for each pixel. The algorithm used in the implementation of the Frost filter is as follows:

K = e (- B * S)

where

B = D * (LV / LM * LM)

S—Absolute value of the pixel distance from the center pixel to its neighbors in the filter window

D—Exponential damping factor (input parameter)

LM—Local mean of filter window

LV—Local variance of filter window

The resulting gray-level value of the filtered pixel is

R = (P1 * K1 + P2 * K2 + ... + Pn * Kn) / (K1 + K2 + ... + Kn)

where

P1,P2,...Pn are gray levels of each pixel in the filter window

K1,K2,...Kn are weights (as defined above) for each pixel

Kuan

The Kuan filter follows a similar filtering process to the Lee filter in reducing speckle noise. This filter also applies a spatial filter to each pixel in an image, filtering the data based on local statistics of the centered pixel value that is calculated using the neighboring pixels.

The Number of Looks parameter controls image smoothing and estimates noise variance; these estimates are used in various ways to control the filter process. The smaller the value is, the better the smoothing effect and filter performance. A larger Number of Looks value retains more image features.

Inputs

The input for the Kuan filter is as follows:

Parameter

Description

Default value

Number of Looks

Specifies the number of looks of the image

1

Kuan filter input

Algorithm

The algorithm for implementation of the Kuan filter is as follows:

The resulting filtered pixel value is:

R = PC * K + LM * (1 - K)

where

CU = 1 / sqrt (NLooks)—Noise variation coefficient

CI = sqrt (LV) / LM—Image variation coefficient

K = (1 - ((CU * CU) / (CI * CI))) / (1 + (CU * CU))

PC—Center pixel value of window

LM—Local mean of filter window

LV—Local variance of filter window

NLooks—Number of looks

For optimal speckle reduction, you can try the following: Different filter sizes greatly affect the quality of processed images. A 7 x 7 filter usually gives the best results. The number of looks is used to estimate noise variance, and it effectively controls the amount of smoothing applied to the image by the filter. A smaller Number of Looks value leads to more smoothing; a larger Number of Looks value preserves more image features.

It is recommended that you apply a histogram stretch to adjust the image's contrast or brightness to help pull out features for optimal display results.

Related Topics

5/18/2014