Changeset 35868 for trunk/ppStack/src/ppStackMatch.c

Timestamp:

Jul 31, 2013, 12:53:02 PM (13 years ago)

Author:

watersc1

Message:

Fix SCL_ values to now have scale = 0.0 when an image is rejected. Fix covariance factor that was causing over-rejection on simple 1DG kernels.

File:

• trunk/ppStack/src/ppStackMatch.c (modified) (1 diff)

trunk/ppStack/src/ppStackMatch.c

@@ -467 +467 @@
                 sum += kernels->rms;
                 num++;
-                //              psLogMsg("ppStack", PS_LOG_INFO, "KERNMATCHDUMP %d %d %g %g %g\n",
-                //                       index,num,kernels->mean,kernels->rms,sum);
+                //                              psLogMsg("ppStack", PS_LOG_INFO, "KERNMATCHDUMP %d %d %g %g %g %g\n",
+                //                                       index,num,kernels->mean,kernels->rms,sum,psImageCovarianceFactor(readout->covariance));
             }
             psFree(iter);
-            options->matchChi2->data.F32[index] = sum / (psImageCovarianceFactor(readout->covariance) * num);
+
+            // Apply the covariance factor to the chi^2 value
+            double covar_factor = 1.0;
+            if (type != PM_SUBTRACTION_KERNEL_SIMPLE) { // Except if we're using simple kernels.
+              covar_factor = psImageCovarianceFactor(readout->covariance);
+            }
+            options->matchChi2->data.F32[index] = sum / (covar_factor * num);
 
         }