CImCombine:SigmaClip
The SigmaClip method merges the pixel values using iterative sigma-clipping to remove high and lower value outliers more than some number of sigmas from the mean at each location.
CImage = CImCombine:SigmaClip( CImageSet, nHigh=2.5, nLow=2.5, nIter=5 ) |
Each pixel in the output image is the mean value of all images except those rejected by the sigma clipping criteria. The value of sigma is computed independently at each location. In order for the SigmaClip method to work properly, normalization must be done as part of the combining.
Suppose a CImageSet named S exists. The following script fragment combines the images by the SigmaClip method by excluding points more then 2.5 sigmas above the mean and 5 sigmas below the mean. The images are normalized by the Median statistic computed inside the central 10% of the image:
|
-- create a CImCombine object |
|
-- central 10% of the image |
|
-- specify CImage and CRect to measure |
|
-- use the SigmaClip as the statistic |
|
-- scale to normalize the image set |
|
-- combine image set and return new image |
|
-- show the new image |
|
-- when done with S, remove it from memory |