Changeset 4215 for trunk/psphot/src/source_moments.c

Timestamp:

Jun 11, 2005, 5:43:47 AM (21 years ago)

Author:

eugene

Message:

adding mask & noise images

File:

• trunk/psphot/src/source_moments.c (modified) (3 diffs)

trunk/psphot/src/source_moments.c

@@ -1 +1 @@
 # include "psphot.h"
 
-psArray *source_moments (psImage *image, psMetadata *config, psArray *allpeaks) 
+psArray *source_moments (psImageData *imdata, psMetadata *config, psArray *allpeaks) 
 {
     bool     status  = false;
@@ -8 +8 @@
 
     psTimerStart ("psphot");
-
-    // reduce the initial peaks to a subset of the better peaks
-    // is this the appropriate place for this culling?
-
-    float XBORDER  = psMetadataLookupF32 (&status, config, "XBORDER");
-    float YBORDER  = psMetadataLookupF32 (&status, config, "YBORDER");
-    float SATURATE = psMetadataLookupF32 (&status, config, "SATURATE");
-
-    psRegion *keep = psRegionAlloc (XBORDER, -XBORDER, YBORDER, -YBORDER);
-    keep           = psRegionForImage (keep, image, keep);
-    // I need to exclude the saturated objects here, but include them again later...?
-    peaks          = pmPeaksSubset (allpeaks, 100000, keep);
 
     // determine properties (sky, moments) of initial sources
@@ -32 +20 @@
     for (int i = 0; i < peaks->n; i++) {
         // provide the user arguments as metadata?
-        psSource *source  = pmSourceLocalSky (image, (psPeak *)peaks->data[i], PS_STAT_SAMPLE_MEDIAN, INNER, OUTER);
+        psSource *source  = pmSourceDefinePixels (imdata, (psPeak *)peaks->data[i], OUTER);
         if (source == NULL) continue;
-        pmSourceMaskRegion (source, keep);
+
+        pmSourceLocalSky_EAM (source, PS_STAT_SAMPLE_MEDIAN, INNER);
         pmSourceMoments (source, RADIUS);
         psArrayAdd (sources, 100, source);