Changeset 42711 for branches/eam_branches/ipp-20230313/psphot/src/psphotMaskReadout.c

Timestamp:

Jul 5, 2024, 12:30:50 PM (2 years ago)

Author:

eugene

Message:

add option in psphot to fit the PSF star positions as well as shapes (default is false); add option to use robustMedian vs clippedMean for moments-based source sizes and PSF star selection; saturate variance to max value if signal is non NAN

File:

• branches/eam_branches/ipp-20230313/psphot/src/psphotMaskReadout.c (modified) (1 diff)

branches/eam_branches/ipp-20230313/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;
+        }
+      }
     }