CImCombine:MtmSigmaClip CImCombine:SigmaClip

CImCombine:RankClip


The RankClip method merges the pixel values after clipping outliers in the upper and lower rank percentiles of the distribution at each point.

Syntax

CImage = CImCombine:RankClip( CImageSet, pHigh=2.5, pLow=2.5 )

Remarks

Each pixel in the output image is the mean value of all images after rejection of images having values inside the upper nHigh percentile and images having pixels in the lower nLow percentile of values at each location.

Example

Suppose a CImageSet named S exists. The following script fragment combines the images by the RankClip method by excluding points above the 95th percentile and below the 2nd percentile. The images are normalized by the Median statistic computed inside the central 10% of the image:

 

C = CImCombine:new()

-- create a CImCombine object

R = CRect:new(10,50,10,50)

-- central 10% of the image

C:SetRelRect( R )

-- specify CImage and CRect to measure

C:SetNormStat("median")

-- use the RankClip as the normalization statistic

C:SetNormMethod("scale")

-- scale to normalize the image set

I = C:RankClip( S, 5, 2 )

-- combine the image set and return a new image

I:Display()

-- show the new image

C:delete()

-- when done with S, remove it from memory

Related Topics

CImCombine, CImageSet, CImage, CRect