Changeset 42842 for trunk/psphot/src/psphotMaskReadout.c

Timestamp:

May 9, 2025, 11:21:40 AM (14 months ago)

Author:

eugene

Message:

add psphotFindFeatures; enable forced-photometry between cameras; do NOT skip negative-flux sources in full-force photometry; make PSF residual component generation more robust: IRLS vs median; exclude pixels with poor values (scaled value > 1.5, relative to normalized PSF); exclude pixels with too few inputs (<5); allow PSF clump to use clipped stats or robust stats

Location:

trunk/psphot

Files:

• . (modified) (1 prop)
• src/psphotMaskReadout.c (modified) (1 diff)

trunk/psphot/src/psphotMaskReadout.c

@@ -58 +58 @@
             return false;
         }
+    }
+
+    // XXX test to repair variance
+    if (true) {
+      psImage *im = readout->image;
+      psImage *wt = readout->variance;
+      float maxVar = 0;
+      // find the max value
+      for (int j = 0; j < im->numRows; j++) {
+        for (int i = 0; i < im->numCols; i++) {
+            if (!isfinite(wt->data.F32[j][i])) continue;
+            maxVar = PS_MAX (wt->data.F32[j][i], maxVar);
+        }
+      }
+      // set pixels with nan variance and non-nan image to maxVar
+      for (int j = 0; j < im->numRows; j++) {
+        for (int i = 0; i < im->numCols; i++) {
+            if (!isfinite(im->data.F32[j][i])) continue;
+            if (isfinite(wt->data.F32[j][i])) continue;
+            wt->data.F32[j][i] = maxVar;
+        }
+      }
     }