Changeset 21519 for trunk/psphot/src/psphotFindPeaks.c

Timestamp:

Feb 16, 2009, 12:34:24 PM (17 years ago)

Author:

eugene

Message:

added additional flag values, threaded functions always use the threaded version

File:

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

trunk/psphot/src/psphotFindPeaks.c

@@ -21 +21 @@
     }
 
-    // correct the peak values to S/N = sqrt(significance)
-    // get the peak flux from the unsmoothed image
-    // the peak pixel coords are guaranteed to be on the image
+    // Convert the peak values to S/N = sqrt(significance).
+    // Get the peak flux from the unsmoothed image.
+    // Rescale the peak position errors using the peak variance
+    // The peak pixel coords are guaranteed to be on the image
     int row0 = readout->image->row0;
     int col0 = readout->image->col0;
@@ -30 +31 @@
         peak->SN = sqrt(peak->value);
         peak->flux = readout->image->data.F32[peak->y-row0][peak->x-col0];
+        if (readout->variance && isfinite (peak->dx)) {
+            peak->dx *= sqrt(readout->variance->data.F32[peak->y-row0][peak->x-col0]);
+        }
+        if (readout->variance && isfinite (peak->dy)) {
+            peak->dy *= sqrt(readout->variance->data.F32[peak->y-row0][peak->x-col0]);
+        }
     }