Changeset 29606 for trunk/psphot/src/psphotDeblendSatstars.c

Timestamp:

Oct 29, 2010, 10:59:01 AM (16 years ago)

Author:

eugene

Message:

do not completely trust the header/concept value for saturation: minimum allowed set in recipe (40k default)

File:

• trunk/psphot/src/psphotDeblendSatstars.c (modified) (1 diff)

trunk/psphot/src/psphotDeblendSatstars.c

@@ -53 +53 @@
 
     pmCell *cell = readout->parent;
-    float SATURATION = 0.75*psMetadataLookupF32 (&status, cell->concepts, "CELL.SATURATION");
+
+    float SATURATION = NAN;
+
+    // do not completely trust the values in the header...
+    float CELL_SATURATION = psMetadataLookupF32 (&status, cell->concepts, "CELL.SATURATION");
+    float MIN_SATURATION = psMetadataLookupF32 (&status, recipe, "DEBLEND_MIN_SATURATION");
+    if (!status || !isfinite(MIN_SATURATION)) {
+        MIN_SATURATION = 40000.0;
+    }
+    if (!isfinite(CELL_SATURATION)) {
+        SATURATION = MIN_SATURATION;
+    } else {
+        SATURATION = PS_MAX(MIN_SATURATION, CELL_SATURATION);
+    }
     float SAT_TEST_LEVEL = 0.5*SATURATION;