Changeset 26612

Timestamp:

Jan 14, 2010, 4:49:42 PM (16 years ago)

Author:

eugene

Message:

limit the number of sources used for the PSF model measurement

File:

• branches/eam_branches/20091201/ppSub/src/ppSubMakePSF.c (modified) (5 diffs)

branches/eam_branches/20091201/ppSub/src/ppSubMakePSF.c

@@ -21 +21 @@
 
 #include "ppSub.h"
+
+psArray *ppSubSelectPSFSources(psArray *sources);
 
 bool ppSubMakePSF(ppSubData *data)
@@ -79 +81 @@
     // Here, we assume the image is background-subtracted
     psArray *sources = psMetadataLookupPtr(&mdok, minuend->analysis, "PSPHOT.SOURCES");
-    if (!psphotReadoutFindPSF(config, view, sources)) {
+    
+    // XXX filter sources?  limit the total number and return only brighter objects?
+    // use flags to toss totally bogus entries?
+    psArray *goodSources = ppSubSelectPSFSources (sources);
+    if (!psphotReadoutFindPSF(config, view, goodSources)) {
         // This is likely a data quality issue
         // XXX Split into multiple cases using error codes?
@@ -86 +92 @@
         ppSubDataQuality(data, PSPHOT_ERR_PSF, PPSUB_FILES_PHOT_SUB | PPSUB_FILES_PHOT_INV);
         psFree(view);
+        psFree(goodSources);
         return true;
     }
@@ -100 +107 @@
     psMetadataRemoveKey(photRO->analysis, "PSPHOT.HEADER");
 
+    psFree(goodSources);
     psFree(view);
 
@@ -156 +164 @@
     return true;
 }
+
+
+// XXX hardwired MAX for now
+# define MAX_NPSF 500
+
+psArray *ppSubSelectPSFSources(psArray *sources){ 
+
+    sources = psArraySort (sources, pmSourceSortBySN);
+
+    psArray *subset = psArrayAllocEmpty(MAX_NPSF);
+
+    int nPSF = 0;
+    for (int i = 0; (nPSF < MAX_NPSF) && (i < sources->n); i++) {
+
+        pmSource *source = sources->data[i];
+        if (!source) continue;
+
+        // skip non-astronomical objects (very likely defects)
+        if (source->type == PM_SOURCE_TYPE_DEFECT) continue;
+        if (source->type == PM_SOURCE_TYPE_SATURATED) continue;
+        if (source->mode & PM_SOURCE_MODE_SATSTAR) continue;
+        
+        psArrayAdd (subset, 100, source);
+        nPSF++;
+    }
+    
+    return subset;
+}