Changeset 10295

Timestamp:

Nov 29, 2006, 1:32:49 PM (19 years ago)

Author:

rhl

Message:

Don't find negative peaks

File:

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

trunk/psphot/src/psphotFindPeaks.c

@@ -40 +40 @@
         for (int i = 0; i < smooth_im->numCols; i++) {
             value = smooth_im->data.F32[j][i];
-            smooth_im->data.F32[j][i] = 0.0;
-            if (mask->data.U8[j][i]) continue;
-            if (smooth_wt->data.F32[j][i] <= 0) continue;
-            smooth_im->data.F32[j][i] = PS_SQR(value) / smooth_wt->data.F32[j][i];
+            if (value < 0 || smooth_wt->data.F32[j][i] <= 0 || mask->data.U8[j][i]) {
+                smooth_im->data.F32[j][i] = 0.0;
+            } else {
+                smooth_im->data.F32[j][i] = PS_SQR(value) / smooth_wt->data.F32[j][i];
+            }
         }
     }