Changeset 14520 for trunk/ppStack/src/ppStackReadout.c

Timestamp:

Aug 15, 2007, 4:18:51 PM (19 years ago)

Author:

Paul Price

Message:

Ensuring seeing, weighting, scale are set to reasonable values.

File:

• trunk/ppStack/src/ppStackReadout.c (modified) (2 diffs)

trunk/ppStack/src/ppStackReadout.c

@@ -52 +52 @@
         psListAdd(cellList, PS_LIST_TAIL, ro->parent);
 
-        float seeing = psMetadataLookupF32(NULL, inputFile->fpa->analysis, "PPSTACK.SEEING"); // Seeing FWHM
-        float weighting = psMetadataLookupF32(NULL, inputFile->fpa->analysis,
+        bool mdok;                      // Status of MD lookup
+        float seeing = psMetadataLookupF32(&mdok, inputFile->fpa->analysis, "PPSTACK.SEEING"); // Seeing FWHM
+        if (!mdok || !isfinite(seeing)) {
+            psWarning("No SEEING supplied for image %d --- set to unity.", fileNum);
+            seeing = 1.0;
+        }
+        float weighting = psMetadataLookupF32(&mdok, inputFile->fpa->analysis,
                                               "PPSTACK.WEIGHTING"); // Relative weighting
-        float scale = psMetadataLookupF32(NULL, inputFile->fpa->analysis, "PPSTACK.SCALE"); // Rel. scale
+        if (!mdok || !isfinite(weighting)) {
+            psWarning("No WEIGHTING supplied for image %d --- set to unity.", fileNum);
+            weighting = 1.0;
+        }
+        float scale = psMetadataLookupF32(&mdok, inputFile->fpa->analysis, "PPSTACK.SCALE"); // Rel. scale
+        if (!mdok || !isfinite(scale)) {
+            psWarning("No SCALE supplied for image %d --- set to unity.", fileNum);
+            scale = 1.0;
+        }
 
         // Brain-dead background subtraction
@@ -69 +82 @@
             return false;
         }
+        psTrace("ppStack", 3, "Background for image %d is %f\n", fileNum, stats->robustMedian);
         (void)psBinaryOp(ro->image, ro->image, "-", psScalarAlloc(stats->robustMedian, PS_TYPE_F32));