Changeset 7436

Timestamp:

Jun 8, 2006, 10:55:52 AM (20 years ago)

Author:

eugene

Message:

account for masked pixels when generating smoothed image

File:

• trunk/psphot/src/psphotFindPeaks.c (modified) (4 diffs)

trunk/psphot/src/psphotFindPeaks.c

@@ -8 +8 @@
     float SIGMA;
     float threshold;
+    float value;
 
     // smooth the image and weight map
@@ -15 +16 @@
     NSIGMA = psMetadataLookupF32 (&status, recipe, "PEAKS_SMOOTH_NSIGMA");
 
+    // XXX we are somewhat doing the wrong thing here with the mask.
+    // it would be better to smooth the image using the mask information
     psImage *smooth_im = psImageCopy (NULL, readout->image, PS_TYPE_F32);
     psImageSmooth (smooth_im, SIGMA, NSIGMA);
@@ -23 +26 @@
     psLogMsg ("psphot", 4, "smooth weight: %f sec\n", psTimerMark ("psphot"));
 
+    psImage *mask = readout->mask;
+
     // psphotSaveImage (NULL, smooth_im, "imsmooth.fits");
     // psphotSaveImage (NULL, smooth_wt, "wtsmooth.fits");
@@ -29 +34 @@
     for (int j = 0; j < smooth_im->numRows; j++) {
         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(smooth_im->data.F32[j][i]) / smooth_wt->data.F32[j][i];
+            smooth_im->data.F32[j][i] = PS_SQR(value) / smooth_wt->data.F32[j][i];
         }
     }