Changeset 30096

Timestamp:

Dec 17, 2010, 9:20:43 AM (15 years ago)

Author:

eugene

Message:

new peaks for new sources need a footprint, at least to determine extended source analysis area

File:

• branches/eam_branches/ipp-20101205/psphot/src/psphotSourceMatch.c (modified) (4 diffs)

branches/eam_branches/ipp-20101205/psphot/src/psphotSourceMatch.c

@@ -187 +187 @@
         pmPhotObj *obj = objects->data[i]; 
 
+        // we will find the input source with the max number of spans and reproduce that footprint
+        int nSpansMax = 0;
+        int iSpansMax = -1;
+
         // mark the images for which sources have been found
         psVectorInit (found, 0);
@@ -196 +200 @@
             psAssert (index < found->n, "invalid index");
 
+            if (src->peak && src->peak->footprint && src->peak->footprint->nspans > nSpansMax) {
+                nSpansMax = src->peak->footprint->nspans;
+                iSpansMax = j;
+            }
+
             found->data.U8[index] = 1;
+        }
+
+        // we make a copy of the largest footprint; this will be used for all new sources associated with this object
+        pmFootprint *footprint = NULL;
+        if (iSpansMax != -1) { // copy the footprint info
+            pmSource *src = obj->sources->data[iSpansMax]; 
+            psAssert(src->peak, "source does not exist?");
+            psAssert(src->peak->footprint, "footprint does not exist");
+            psAssert(src->peak->footprint->nspans == nSpansMax, "wrong footprint?");
+            
+            // we only care about the spans, do not worry about the image of this footprint
+            footprint = pmFootprintCopyData(src->peak->footprint, NULL);
         }
 
@@ -219 +240 @@
             peak->dy = NAN;
             
-            // XXX assign to a footprint?
+            // assign to a footprint on this readout->image
+            peak->footprint = pmFootprintCopyData(footprint, readout->image);
 
             // create a new source
@@ -241 +263 @@
             psFree (peak);
         }
+        psFree (footprint);
     }