Changeset 42710 for branches/eam_branches/ipp-20230313/psModules/src/objects/pmSource.c

Timestamp:

Jul 5, 2024, 12:29:41 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; reduce chance of choosing galaxies for PSF stars

File:

• branches/eam_branches/ipp-20230313/psModules/src/objects/pmSource.c (modified) (4 diffs)

branches/eam_branches/ipp-20230313/psModules/src/objects/pmSource.c

@@ -407 +407 @@
 *****************************************************************************/
 
-pmPSFClump pmSourcePSFClump(psImage **savedImage, psRegion *region, psArray *sources, float PSF_SN_LIM, float PSF_CLUMP_GRID_SCALE, psF32 SX_MAX, psF32 SY_MAX, psF32 SX_MIN, psF32 SY_MIN, psF32 AR_MAX)
+pmPSFClump pmSourcePSFClump(psImage **savedImage, psRegion *region, psArray *sources, float PSF_SN_LIM, float PSF_CLUMP_GRID_SCALE, psF32 SX_MAX, psF32 SY_MAX, psF32 SX_MIN, psF32 SY_MIN, psF32 AR_MAX, bool useClippedMean)
 {
     psTrace("psModules.objects", 10, "---- begin ----\n");
@@ -601 +601 @@
         }
 
-        // measures stats of Sx, Sy
-        stats = psStatsAlloc (PS_STAT_CLIPPED_MEAN | PS_STAT_CLIPPED_STDEV);
-
-        if (!psVectorStats (stats, tmpSx, NULL, NULL, 0)) {
-            psError(PS_ERR_UNKNOWN, false, "failed to measure Sx stats");
-            return (emptyClump);
-        }
-        psfClump.X  = stats->clippedMean;
-        psfClump.dX = hypot(stats->clippedStdev, PSF_CLUMP_GRID_SCALE);
-
-        if (!psVectorStats (stats, tmpSy, NULL, NULL, 0)) {
-            psError(PS_ERR_UNKNOWN, false, "failed to measure Sy stats");
-            return (emptyClump);
-        }
-        psfClump.Y  = stats->clippedMean;
-        psfClump.dY = hypot(stats->clippedStdev, PSF_CLUMP_GRID_SCALE);
+        // measures stats of Sx, Sy
+        if (useClippedMean) {
+          stats = psStatsAlloc (PS_STAT_CLIPPED_MEAN | PS_STAT_CLIPPED_STDEV);
+        } else {
+          stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
+        }
+
+        if (!psVectorStats (stats, tmpSx, NULL, NULL, 0)) {
+          psError(PS_ERR_UNKNOWN, false, "failed to measure Sx stats");
+          return (emptyClump);
+        }
+        psfClump.X  = psStatsGetValue (stats, psStatsMeanOption (stats->options)); 
+        psfClump.dX = psStatsGetValue (stats, psStatsStdevOption (stats->options));
+        psfClump.dX = hypot(psfClump.dX, PSF_CLUMP_GRID_SCALE);
+
+        if (!psVectorStats (stats, tmpSy, NULL, NULL, 0)) {
+          psError(PS_ERR_UNKNOWN, false, "failed to measure Sy stats");
+          return (emptyClump);
+        }
+        psfClump.Y  = psStatsGetValue (stats, psStatsMeanOption (stats->options)); 
+        psfClump.dY = psStatsGetValue (stats, psStatsStdevOption (stats->options));
+        psfClump.dY = hypot(psfClump.dX, PSF_CLUMP_GRID_SCALE);              
 
         psTrace ("psModules.objects", 2, "clump  X,  Y: %f, %f\n", psfClump.X, psfClump.Y);
@@ -749 +755 @@
             }
 
+            psF32 radius = hypot ((sigX-clump.X)/clump.dX, (sigY-clump.Y)/clump.dY);
+
+            // XXX this is crude cut to avoid using extended sources in the PSF model
+            // psphot will update this analysis based on the PSF - Kron mags
             // likely unsaturated extended source (too large to be stellar)
-            if (sigX > clump.X + 3*clump.dX || sigY > clump.Y + 3*clump.dY) {
+            // XXX old test: if (sigX > clump.X + 3*clump.dX || sigY > clump.Y + 3*clump.dY) {
+
+            if (radius > 2.0) {
                 source->type = PM_SOURCE_TYPE_EXTENDED;
                 Next ++;
@@ -768 +780 @@
 
             // PSF star (within 1.5 sigma of clump center, S/N > limit)
-            psF32 radius = hypot ((sigX-clump.X)/clump.dX, (sigY-clump.Y)/clump.dY);
             if ((source->moments->SN > PSF_SN_LIM) && (radius < PSF_CLUMP_NSIGMA)) {
                 source->type = PM_SOURCE_TYPE_STAR;