Changeset 5084 for trunk/psphot/src/psPolynomials.c

Timestamp:

Sep 20, 2005, 6:48:16 PM (21 years ago)

Author:

eugene

Message:

alternative clipping for pmPSFtryMetric

File:

• trunk/psphot/src/psPolynomials.c (modified) (4 diffs)

trunk/psphot/src/psPolynomials.c

@@ -507 +507 @@
     psVector *zFit   = NULL;
     psVector *zResid = psVectorAlloc (x->n, PS_TYPE_F64);
-
-    // XXX EAM : use SAMPLE_MEAN and SAMPLE_STDEV for stats:
+    float minClipSigma = stats->min;
+    float maxClipSigma = stats->max;
+    float minClipValue;
+    float maxClipValue;
+
+    // XXX EAM : use SAMPLE_MEDIAN and SAMPLE_STDEV for stats:
     stats->options |= (PS_STAT_SAMPLE_MEDIAN | PS_STAT_SAMPLE_STDEV);
 
-    for (int N = 0; N < 3; N++) {
+    for (int N = 0; N < stats->clipIter; N++) {
         int Nkeep = 0;
 
@@ -519 +523 @@
 
         stats  = psVectorStats (stats, zResid, NULL, mask, maskValue);
-        psTrace (".psphot.RobustFit", 4, "residual stats for robust fit:  %g +/- %g\n", 
+        minClipValue = -1*fabs(minClipSigma)*stats->sampleStdev;
+        maxClipValue = +1*fabs(maxClipSigma)*stats->sampleStdev;
+
+        psTrace (".psphot.VectorClipFit", 4, "residual stats for robust fit:  %g +/- %g\n", 
                  stats->sampleMedian, stats->sampleStdev);
 
@@ -527 +534 @@
         for (int i = 0; i < zResid->n; i++) {
             if (mask->data.U8[i]) continue;
-            if (fabs(zResid->data.F64[i] - stats->sampleMedian) > 2*stats->sampleStdev) {
+            if (zResid->data.F64[i] - stats->sampleMedian > maxClipValue) {
+                mask->data.U8[i] |= 0x01;
+                continue;
+            }       
+            if (zResid->data.F64[i] - stats->sampleMedian < minClipValue) {
                 mask->data.U8[i] |= 0x01;
                 continue;
@@ -534 +545 @@
         }
 
-        psTrace (".psphot.RobustFit", 4, "keeping %d of %d pts for fit\n", 
+        psTrace (".psphot.VectorClipFit", 4, "keeping %d of %d pts for fit\n", 
                  Nkeep, x->n);