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

Timestamp:

Feb 26, 2010, 2:27:07 PM (16 years ago)

Author:

Paul Price

Message:

Make 1/noise^{2 weighting optional. This is an attempt to reduce the massive print-through with PR images of M42, where there is no background subtraction.}

File:

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

trunk/ppStack/src/ppStackMatch.c

@@ -628 +628 @@
 
     // Measure the variance level for the weighting
-    if (!psImageBackground(bg, NULL, readout->variance, readout->mask, maskVal | maskBad, rng)) {
-        psError(PPSTACK_ERR_DATA, false, "Can't measure mean variance for image.");
-        psFree(rng);
-        psFree(bg);
-        return false;
-    }
-    options->weightings->data.F32[index] = 1.0 / (psStatsGetValue(bg, PS_STAT_ROBUST_MEDIAN) *
-                                                  psImageCovarianceFactor(readout->covariance));
+    if (psMetadataLookupBool(NULL, recipe, "WEIGHTS")) {
+        if (!psImageBackground(bg, NULL, readout->variance, readout->mask, maskVal | maskBad, rng)) {
+            psError(PPSTACK_ERR_DATA, false, "Can't measure mean variance for image.");
+            psFree(rng);
+            psFree(bg);
+            return false;
+        }
+        options->weightings->data.F32[index] = 1.0 / (psStatsGetValue(bg, PS_STAT_ROBUST_MEDIAN) *
+                                                      psImageCovarianceFactor(readout->covariance));
+    } else {
+        options->weightings->data.F32[index] = 1.0;
+    }
     psLogMsg("ppStack", PS_LOG_INFO, "Weighting for image %d is %f\n",
              index, options->weightings->data.F32[index]);
-    psMetadataAddF32(readout->analysis, PS_LIST_TAIL, "PPSTACK.WEIGHTING", 0,
-                     "Weighting by 1/noise^2 for stack", options->weightings->data.F32[index]);
 
     psFree(rng);