Changeset 42715 for branches/2dbias/psModules/src/detrend/pmOverscan.c

Timestamp:

Aug 6, 2024, 5:35:16 PM (2 years ago)

Author:

hgao

Message:

enable clipped mean in psVectorBoxcar

File:

• branches/2dbias/psModules/src/detrend/pmOverscan.c (modified) (6 diffs)

branches/2dbias/psModules/src/detrend/pmOverscan.c

@@ -209 +209 @@
 
                 // Reduce the overscans
-                psVector *reduced = pmOverscanVector(&chi2, overscanOpts->primary, pixels);
+                psVector *reduced = pmOverscanVector(&chi2, overscanOpts->primary, pixels, false);
                 psFree(pixels);
                 if (!reduced)
@@ -293 +293 @@
 
                 // Reduce the overscans
-                psVector *reduced = pmOverscanVector(&chi2, overscanOpts->primary, pixels);
+                psVector *reduced = pmOverscanVector(&chi2, overscanOpts->primary, pixels, false);
                 psFree(pixels);
                 if (!reduced)
@@ -415 +415 @@
                 // Reduce the overscans
                 // XXX need to save 2 different chi-square values
-                psVector *yReduced = pmOverscanVector(&chi2, overscanOpts->primary, yscanPixels);
+                psVector *yReduced = pmOverscanVector(&chi2, overscanOpts->primary, yscanPixels, false);
                 psFree(yscanPixels);
                 if (!yReduced)
@@ -430 +430 @@
                 // Reduce the overscans
                 // XXX need to save 2 different chi-square values
-                psVector *xReduced = pmOverscanVector(&chi2, overscanOpts->secondary, xscanPixels);
+                psVector *xReduced = pmOverscanVector(&chi2, overscanOpts->secondary, xscanPixels, true);
                 psFree(xscanPixels);
                 if (!xReduced)
@@ -560 +560 @@
 psVector *pmOverscanVector(float *chi2,                                                 // chi^2 from fit
                                                    pmOverscanStatOptions *overscanOpts, // Overscan statistic options
-                                                   const psArray *pixels                                // Array of vectors containing the pixel values
-)
+                                                   const psArray *pixels,                               // Array of vectors containing the pixel values
+                                                   bool robust)                                                 // Use robust statistics for boxcar smoothing
 {
         assert(overscanOpts);
@@ -611 +611 @@
         if (overscanOpts->boxcar > 0)
         {
-                psVector *smoothed = psVectorBoxcar(NULL, reduced, overscanOpts->boxcar); // Smoothed vector
+                psVector *smoothed = psVectorBoxcar(NULL, reduced, overscanOpts->boxcar, robust); // Smoothed vector
                 psFree(reduced);
                 reduced = smoothed;