Changeset 4977 for trunk/psphot/src/pmSourceUtils.c

Timestamp:

Sep 7, 2005, 6:32:40 PM (21 years ago)

Author:

eugene

Message:

cleanup, organization

File:

• trunk/psphot/src/pmSourceUtils.c (modified) (1 diff)

trunk/psphot/src/pmSourceUtils.c

@@ -1 +1 @@
 # include "psphot.h"
-
-bool pmSourceDefinePixels(pmSource *mySource, 
-                          const psImageData *imdata,
-                          psF32 x, 
-                          psF32 y,
-                          psF32 Radius)
-{
-    psRegion srcRegion;
-
-    // Grab a subimage of the original image of size (2 * outerRadius).
-    srcRegion = psRegionForSquare (x, y, Radius);
-    srcRegion = psRegionForImage (imdata->image, &srcRegion);
-
-    mySource->pixels = psImageSubset(imdata->image, srcRegion);
-    mySource->noise  = psImageSubset(imdata->noise, srcRegion);
-    mySource->mask   = psImageSubset(imdata->mask,  srcRegion);
-
-    return(mySource);
-}
-
-bool pmSourcePhotometry (float *fitMag, float *obsMag, pmModel *model, psImage *image, psImage *mask) {
-
-    float obsSum = 0;
-    float fitSum = 0;
-    float sky = model->params->data.F32[0];
-
-    pmModelFlux modelFluxFunc = pmModelFlux_GetFunction (model->type);
-    fitSum = modelFluxFunc (model->params);
-
-    for (int ix = 0; ix < image->numCols; ix++) {
-        for (int iy = 0; iy < image->numRows; iy++) {
-            if (mask->data.U8[iy][ix]) continue;
-            obsSum += image->data.F32[iy][ix] - sky;
-        }
-    }
-    if (obsSum <= 0) return false;
-    if (fitSum <= 0) return false;
-
-    *fitMag = -2.5*log10(fitSum);
-    *obsMag = -2.5*log10(obsSum);
-    return (true);
-} 
-