Changeset 6556 for branches/rel10_ifa/psModules/src/objects/pmSourceIO_CMP.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_CMP.c (modified) (8 diffs)

branches/rel10_ifa/psModules/src/objects/pmSourceIO_CMP.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 "psLine.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 pseudo FITS table (header + ascii list)
@@ -10 +36 @@
     float dmag, lsky;
     bool status;
+    EllipseShape shape;
+    EllipseAxes axes;
 
     // find config information for output header
@@ -19 +47 @@
     for (i = nSrc = 0; i < sources->n; i++) {
         pmSource *source = (pmSource *) sources->data[i];
-        pmModel *model = pmModelSelect (source);
+        pmModel *model = pmSourceSelectModel (source);
         if (model == NULL)
             continue;
@@ -49 +77 @@
     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;
@@ -56 +84 @@
         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);
+
         psLineInit (line);
-        psLineAdd (line, "%6.1f ", PAR[2]);
-        psLineAdd (line, "%6.1f ", PAR[3]);
-        psLineAdd (line, "%6.3f ", PS_MIN (99.999, source->fitMag + ZERO_POINT));
-        psLineAdd (line, "%03d ",  PS_MIN (999, (int)(1000*dmag)));
-        psLineAdd (line, "%2d ",   type);
-        psLineAdd (line, "%3.1f ", lsky);
-        psLineAdd (line, "%6.3f ", 99.999); // should be 'Mgal
-        psLineAdd (line, "%6.3f ", PS_MIN (99.999, source->apMag + ZERO_POINT));
-        psLineAdd (line, "%6.2f ", PAR[4]); // should be 'FHWM x'
-        psLineAdd (line, "%6.2f ", PAR[5]); // should be 'FHWM y'
-        psLineAdd (line, "%5.1f\n", 0); // should be theta
+        psLineAdd (line, "%6.1f ",  PAR[2]);
+        psLineAdd (line, "%6.1f ",  PAR[3]);
+        psLineAdd (line, "%6.3f ",  PS_MIN (99.0, source->psfMag + ZERO_POINT));
+        psLineAdd (line, "%03d ",   PS_MIN (999, (int)(1000*dmag)));
+        psLineAdd (line, "%2d ",    type);
+        psLineAdd (line, "%3.1f ",  lsky);
+        psLineAdd (line, "%6.3f ",  PS_MIN (99.0, source->extMag + ZERO_POINT));
+        psLineAdd (line, "%6.3f ",  PS_MIN (99.0, source->apMag  + ZERO_POINT));
+        psLineAdd (line, "%6.2f ",  axes.major);
+        psLineAdd (line, "%6.2f ",  axes.minor);
+        psLineAdd (line, "%5.1f\n", axes.theta);
         fwrite (line->line, 1, line->Nline, f);
     }
@@ -82 +114 @@
 
 // elixir-style pseudo FITS table (header + ascii list)
-psArray *sources pmSourcesReadCMP (char *filename, psMetadata *header)
+psArray *pmSourcesReadCMP (char *filename, psMetadata *header)
 {
 
-    int i, type, nSrc;
-    psMetadataItem *mdi;
+    int Ninstar;
     psF32 *PAR, *dPAR;
-    float dmag, lsky;
     bool status;
+    EllipseShape shape;
+    EllipseAxes axes;
+
+    // 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
@@ -113 +153 @@
 
     // prepare array to store data
-    int nStars = psMetadataLookupS32 (&status, myHeader, "NSTARS");
+    int nStars = psMetadataLookupS32 (&status, header, "NSTARS");
     psArray *sources = psArrayAlloc (nStars);
     sources->n = 0;
 
     // we have fixed bytes / line : use that info
-    char *buffer = psAlloc (NBYTES_STAR*nStars);
+    char *buffer = psAlloc (BYTES_STAR*nStars);
 
     int Nextra = 0;
@@ -161 +201 @@
 
         /* extract data for stars */
-        int Ninstar = nbytes / BYTES_STAR;
+        Ninstar = nbytes / BYTES_STAR;
         Nextra = nbytes % BYTES_STAR;
-        for (j = 0; j < Ninstar; j++, N++) {
-            psString line = psStringNCopy (&buffer[j*BYTES_STAR]);
+        for (int j = 0; j < Ninstar; j++) {
+            psString line = psStringNCopy (&buffer[j*BYTES_STAR], BYTES_STAR);
             psList *list = psStringSplit (line, " ");
             psArray *array = psListToArray (list);
 
             pmSource *source = pmSourceAlloc ();
-
-            source->x      = atof (array->data[0]);
-            source->y      = atof (array->data[1]);
-            source->fitMag = atof (array->data[2]);
-            source->dMag   = 0.001*atof (array->data[3]);
-            source->type   = atof (array->data[4]);
-            source->sky    = pow (atof (array->data[5]), 10.0);
-            source->galMag = atof (array->data[6]);
+            pmModel *model = pmModelAlloc (modelType);
+            source->modelPSF  = model;
+
+            PAR = model->params->data.F32;
+            dPAR = model->dparams->data.F32;
+
+            PAR[0]         = pow (atof (array->data[5]), 10.0);
+            PAR[2]         = atof (array->data[0]);
+            PAR[3]         = atof (array->data[1]);
+            source->psfMag = atof (array->data[2]);
+            source->extMag = atof (array->data[6]);
+            source->errMag = atof (array->data[3]) / 1000.0;
             source->apMag  = atof (array->data[7]);
-            source->sx     = atof (array->data[8]);
-            source->sy     = atof (array->data[9]);
-            source->theta  = atof (array->data[10]);
+            axes.major     = atof (array->data[8]);
+            axes.minor     = atof (array->data[9]);
+            axes.theta  = atof (array->data[10]);
+
+            shape = EllipseAxesToShape (axes);
+
+            PAR[4] = shape.sx;
+            PAR[5] = shape.sy;
+            PAR[6] = shape.sxy;
+
+            // source->type   = atof (array->data[4]);
 
             psArrayAdd (sources, 100, source);