FELEE -- Enhanced Lee SAR Speckle Filter

Applies an Enhanced Lee speckle filter on a SAR image. This filter is primarily used on radar data to remove high frequency noise (speckle), while preserving high frequency features (edges).

See Also: FSPEC, FAV, FED, FGA, FME, FMO, FSHARP, FSOBEL, FPR, FPRE

INPUT PORT(S)

Input

Type:Raster
Connection:Mandatory
Minimum Layers:1
Maximum Layers:1024

Contains the input layer(s) to be filtered.

Mask

Type:Bitmap
Connection:Optional
Minimum Layers:0
Maximum Layers:1

Contains the input bitmap layer that indicates the area under which the input raster should be processed.

OUTPUT PORT(S)

Output

Type:Raster
Connection:Optional
Minimum Layers:0
Maximum Layers:1024

Contains the output layer(s) that will receive the filtered image data. Only the area under "Mask" is written to "Output".

INPUT PARAMETERS

FELEE is controlled by the following parameters:

Image Type

Name:ImageType
Type:Text
Valid Values:AMP, POW
Default:AMP
Requirement:Mandatory

Specifies the image format of the radar image, which can be supplied in either the amplitude or power (intensity) format.

Amplitude is the square root of power. Most radar images are supplied in the amplitude format to preserve the values. 

       AMP     | image is in amplitude format
       POW     | image is in power format

Filter X Size (Pixels)

Name:FiltPixels
Type:Integer
Valid Values:3, 5, 7, 9, 11
Default:3
Requirement:Mandatory

Specifies the filter's X size in units of pixels. The filter must be square (i.e., the X size must equal the Y size).

Filter Y Size (Lines)

Name:FiltLines
Type:Integer
Valid Values:3, 5, 7, 9, 11
Default:3
Requirement:Mandatory

Specifies the filter's Y size in units of lines. The filter must be square (i.e., the X size must equal the Y size).

Number of Looks

Name:NumberLooks
Type:Integer
Valid Values:0 <= x <= 100
Default:1
Requirement:Mandatory

Specifies the (effective) number of looks of the image. This parameter is used to derive noise variance. By adjusting NumberLooks, the user can control the amount of smoothing applied to the image.

Damping Factor

Name:DampFactor
Type:Real
Valid Values:0.0 <= x <= 10.0
Default:1.0
Requirement:Mandatory

Specifies the damping constant for the filter. This constant specifies the extent of the damping effect of the filtering. The default value of 1.0 is sufficient for most SAR images.

The use of large values for DampFactor allows for better preservation of sharp edges, but reduces the smoothing effect. The use of small values for DampFactor increases the smoothing effect, but does not preserve sharp edges well.

DETAILS

The Enhanced Lee filter is used primarily to filter speckled radar data. It is designed to smooth out noise while retaining edges or shape features in the image.

The filter size can be specified through the FilterSize parameters. Different filter sizes will greatly affect the quality of processed images. If the filter is too small, the noise filtering algorithm is not effective. If the filter is too large, subtle details of the image will be lost in the filtering process. A 7x7 filter usually gives the best results.

The NumberLooks parameter is used to estimate noise variance and it effectively controls the amount of smoothing applied to the image by the filter. Theoretically, the correct value for NumberLooks should be the effective number of looks of the radar image. It should be close to the actual number of looks, but may be different if the image has undergone resampling. The user may experimentally adjust the NumberLooks value so as to control the effect of the filter. A smaller NumberLooks value leads to more smoothing; a larger NumberLooks value preserves more image features.

A damping factor (DampFactor) is required by the Enhanced Lee filter. The value of DampFactor defines the extent of exponential damping (the smaller the value, the smaller the damping effect). It depends on the non-filtered image and may require trial-and-error experiments to determine the best value. The default value for DampFactor is 1.

FELEE performs spatial filtering on each individual pixel in an image using the grey level values in a square window surrounding each pixel. The dimensions of the filter must be odd, and must be at least 3x3.

All pixels are filtered. In order to filter pixels located near the edges of the image, edge-pixel values are replicated to give sufficient data.

A bitmap specifies the area within the input layer which will be processed. Only this area will be filtered and the rest of the image will be unchanged. If no bitmap is connected, the entire database is processed.

The Enhanced Lee filter model requires that the signal represents power. If the input image is in amplitude format, each grey level will be squared to derive power and finally square root will be applied to the filtered result.

ALGORITHM

Implementation of the speckle filters are based on the following papers, and especially the review paper by Shi and Fung. 

 Jong-Sen Lee, "Digital Image Enhancement and Noise Filtering
 by Use of Local Statistics", IEEE Transactions on Pattern
 Analysis and Machine Intelligence, Vol. PAM1-2, No. 2, March, 1980.

 J.S.Lee, "Refined Filtering of Image Noise Using Local Statistics"
 Computer Graphic and Image Processing 15, 380-389 (1981)

 D.T. Kuan, A.A. Sawchuk, T.C. Strand, and P. Chavel, 
 "Adaptive restoration of images with speckle," IEEE Trans. ASSP.,
 Vol. 35, no. 3, pp. 373-383, March 1987.

 A. Lopes, R. Touzi and E. Nezry, "Adaptive speckle filters and
 Scene heterogeneity", IEEE Transaction on Geoscience and Remote
 Sensing, Vol. 28, No. 6, pp. 992-1000, Nov. 1990.

 V.S. Frost, J.A. Stiles, K.S. Shanmugan, and J.C. Holtzman,
 "A model for radar images and its application to adaptive digital 
 filtering of multiplicative noise," IEEE Trans. Pattern Analysis
 and Machine Intelligence, vol. 4, no. 2, pp. 157-166, March 1982. 

 A. Lopes, E. Nezry, R. Touzi, and H. Laur, "Structure detection
 and statistical adaptive speckle filtering in SAR images",
 International Journal of Remote Sensing, Vol. 14, No. 9, 
 pp. 1735-1758, 1993.

 A. Lopes, R. Touzi and E. Nezry, "Adaptive speckle filters and
 Scene heterogeneity", IEEE Transaction on Geoscience and Remote
 Sensing, Vol. 28, No. 6, pp. 992-1000, Nov. 1990.

 Zhenghao Shi and Ko B. Fung, 1994, "A Comparison of Digital
 Speckle Filters", Proceedings of IGRASS 94, August 8-12, 1994.

ENHANCED LEE FILTER

Enhanced Lee filter was modified by Lopes. He propose to divide an image into areas of three classes. The first class corresponds to the homogeneous areas in which the speckles may be eliminated simply by applying an LP filter (or equivalently, averaging, multi-look processing). The second class corresponds to the heterogeneous areas in which the speckles are to be reduced while preserving texture. And the third class are areas containing isolated point targets, which, in this case, the filter should preserve the observed value.

The resulting grey-level value R for the smoothed pixel is: 

       R = Im                  for Ci <= Cu
       R = Im * W + Ic * (1-W) for Cu < Ci < Cmax
       R = Ic                  for Ci >= Cmax
where: 

       W   =   exp (-DampFactor (Ci-Cu)/(Cmax - Ci))
       Cu  =   SQRT(1/NumberLooks)
       Ci  =   S / Im
       Cmax =  SQRT(1+2/NumberLooks)           
       Ic  =   center pixel in filter window
       Im  =   mean value of intensity within window
       S   =   standard deviation of intensity within window
About PCI Help Gateway