Changeset 6556 for branches/rel10_ifa/psModules/src/objects/pmSourceIO_CMF.c

Timestamp:

Mar 8, 2006, 5:14:23 PM (20 years ago)

Author:

magnier

Message:

major rework of objects code

File:

• branches/rel10_ifa/psModules/src/objects/pmSourceIO_CMF.c (modified) (6 diffs)

branches/rel10_ifa/psModules/src/objects/pmSourceIO_CMF.c

@@ -1 @@
-# include "pmSource.h"
+/** @file  pmSourceIO.c
+ *
+ *  @author EAM, IfA
+ *
+ *  @version $Revision: 1.1.2.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-09 03:14:23 $
+ *
+ *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
+ *
+ */
+
+#include <stdio.h>
+#include <math.h>
+#include <string.h>
+#include "pslib.h"
+#include "psEllipse.h"
+#include "pmHDU.h"
+#include "pmFPA.h"
+#include "pmPeaks.h"
+#include "pmMoments.h"
+#include "pmGrowthCurve.h"
+#include "pmModel.h"
+#include "pmPSF.h"
+#include "pmSource.h"
+#include "pmModelGroup.h"
+#include "pmSourceIO.h"
 
 // elixir-style FITS table output (header + table in 1st extension)
@@ -12 +37 @@
     int i, type;
     psF32 *PAR, *dPAR;
-    float dmag, lsky;
+    float lsky;
     bool status;
+    EllipseShape shape;
+    EllipseAxes axes;
 
     // find config information for output header
@@ -23 +50 @@
     for (i = 0; i < sources->n; i++) {
         pmSource *source = (pmSource *) sources->data[i];
-        pmModel *model = pmModelSelect (source);
+        pmModel *model = pmSourceSelectModel (source);
         if (model == NULL)
             continue;
@@ -30 +57 @@
         dPAR = model->dparams->data.F32;
 
-        dmag = dPAR[1] / PAR[1];
         type = pmSourceDophotType (source);
         lsky = (PAR[0] < 1.0) ? 0.0 : log10(PAR[0]);
+
+        shape.sx  = PAR[4];
+        shape.sy  = PAR[5];
+        shape.sxy = PAR[6];
+        axes = EllipseShapeToAxes (shape);
 
         row = psMetadataAlloc ();
         psMetadataAdd (row, PS_LIST_TAIL, "X_PIX",   PS_DATA_F32, "", PAR[2]);
         psMetadataAdd (row, PS_LIST_TAIL, "Y_PIX",   PS_DATA_F32, "", PAR[3]);
-        psMetadataAdd (row, PS_LIST_TAIL, "MAG_RAW", PS_DATA_F32, "", source->fitMag + ZERO_POINT);
-        psMetadataAdd (row, PS_LIST_TAIL, "MAG_ERR", PS_DATA_F32, "", (int)(1000*dmag));
-        psMetadataAdd (row, PS_LIST_TAIL, "MAG_GAL", PS_DATA_F32, "", 32.0);
-        psMetadataAdd (row, PS_LIST_TAIL, "MAG_AP",  PS_DATA_F32, "", source->apMag + ZERO_POINT);
+        psMetadataAdd (row, PS_LIST_TAIL, "MAG_RAW", PS_DATA_F32, "", PS_MIN (99.0, source->psfMag + ZERO_POINT));
+        psMetadataAdd (row, PS_LIST_TAIL, "MAG_GAL", PS_DATA_F32, "", PS_MIN (99.0, source->extMag + ZERO_POINT));
+        psMetadataAdd (row, PS_LIST_TAIL, "MAG_AP",  PS_DATA_F32, "", PS_MIN (99.0, source->apMag + ZERO_POINT));
+        psMetadataAdd (row, PS_LIST_TAIL, "MAG_ERR", PS_DATA_F32, "", PS_MIN (999, 1000*source->errMag));
         psMetadataAdd (row, PS_LIST_TAIL, "LOG_SKY", PS_DATA_F32, "", lsky);
-        psMetadataAdd (row, PS_LIST_TAIL, "FWHM_X",  PS_DATA_F32, "", type);
-        psMetadataAdd (row, PS_LIST_TAIL, "FWHM_Y",  PS_DATA_F32, "", PAR[4]);
-        psMetadataAdd (row, PS_LIST_TAIL, "THETA",   PS_DATA_F32, "", PAR[5]);
-        psMetadataAdd (row, PS_LIST_TAIL, "DOPHOT",  PS_DATA_STRING, "", "0");
+        psMetadataAdd (row, PS_LIST_TAIL, "FWHM_X",  PS_DATA_F32, "", axes.major);
+        psMetadataAdd (row, PS_LIST_TAIL, "FWHM_Y",  PS_DATA_F32, "", axes.minor);
+        psMetadataAdd (row, PS_LIST_TAIL, "THETA",   PS_DATA_F32, "", axes.theta);
+        psMetadataAdd (row, PS_LIST_TAIL, "DOPHOT",  PS_DATA_U8,  "", type);
         psMetadataAdd (row, PS_LIST_TAIL, "DUMMY",   PS_DATA_STRING, "", "123");
 
@@ -62 +93 @@
 
 // read in a readout from the fits file
-psArray *pmReadoutReadCMF (psFits *fits, psMetadata *header)
+psArray *pmSourcesReadCMF (psFits *fits, psMetadata *header)
 {
 
     bool status;
+    psF32 *PAR, *dPAR;
+    EllipseShape shape;
+    EllipseAxes axes;
+    float lsky;
+
+    // define PSF model type
+    int modelType = pmModelSetType ("PS_MODEL_GAUSS");
+
+    char *PSF_NAME = psMetadataLookupStr (&status, header, "PSF_NAME");
+    if (PSF_NAME != NULL) {
+        modelType = pmModelSetType (PSF_NAME);
+    }
 
     // find config information for output header
     float ZERO_POINT = psMetadataLookupF32 (&status, header, "ZERO_PT");
+    if (!status)
+        ZERO_POINT = 25.0;
 
     psArray *table = psFitsReadTable (fits);
@@ -76 +121 @@
 
     // convert the table to the pmSource entries
+    // XXX need to chooose PSF vs EXT, based on type?
     for (int i = 0; i < table->n; i++) {
-        pmStar *source = pmSourceAlloc ();
+        pmSource *source = pmSourceAlloc ();
+        pmModel *model = pmModelAlloc (modelType);
+        source->modelPSF  = model;
+
+        PAR = model->params->data.F32;
+        dPAR = model->dparams->data.F32;
 
         psMetadata *row = table->data[i];
 
-        source->x      = psMetadataLookupF32 (&status, row, "X_PIX");
-        source->y      = psMetadataLookupF32 (&status, row, "Y_PIX");
-        source->fitMag = psMetadataLookupF32 (&status, row, "MAG_RAW") - ZERO_POINT;
-        source->galMag = psMetadataLookupF32 (&status, row, "MAG_GAL") - ZERO_POINT;
+        lsky           = psMetadataLookupF32 (&status, row, "LOG_SKY");
+        PAR[0]         = pow(lsky, 10.0);
+        PAR[2]         = psMetadataLookupF32 (&status, row, "X_PIX");
+        PAR[3]         = psMetadataLookupF32 (&status, row, "Y_PIX");
+        axes.major     = psMetadataLookupF32 (&status, row, "FWHM_X");
+        axes.minor     = psMetadataLookupF32 (&status, row, "FWHM_Y");
+        axes.theta     = psMetadataLookupF32 (&status, row, "THETA");
+
+        shape = EllipseAxesToShape (axes);
+
+        PAR[4] = shape.sx;
+        PAR[5] = shape.sy;
+        PAR[6] = shape.sxy;
+
+        source->psfMag = psMetadataLookupF32 (&status, row, "MAG_RAW") - ZERO_POINT;
+        source->extMag = psMetadataLookupF32 (&status, row, "MAG_GAL") - ZERO_POINT;
+        source->errMag = psMetadataLookupF32 (&status, row, "MAG_ERR");
         source->apMag  = psMetadataLookupF32 (&status, row, "MAG_AP")  - ZERO_POINT;
-        source->dMag   = psMetadataLookupF32 (&status, row, "MAG_ERR");
-        source->sky    = psMetadataLookupF32 (&status, row, "LOG_SKY");
-        source->sx     = psMetadataLookupF32 (&status, row, "FWHM_X");
-        source->sy     = psMetadataLookupF32 (&status, row, "FWHM_Y");
-        source->theta  = psMetadataLookupF32 (&status, row, "THETA");
 
         sources->data[i] = source;