Changeset 27924

Timestamp:

May 11, 2010, 5:39:56 PM (16 years ago)

Author:

eugene

Message:

adding new structure and IO function to support simple circular radial apertures

Location:

branches/eam_branches/psModules.20100506/src/objects

Files:

• Makefile.am (modified) (1 diff)
• pmSource.c (modified) (2 diffs)
• pmSource.h (modified) (1 diff)
• pmSourceExtendedPars.c (modified) (1 diff)
• pmSourceExtendedPars.h (modified) (2 diffs)
• pmSourceIO.c (modified) (5 diffs)
• pmSourceIO.h (modified) (2 diffs)
• pmSourceIO_CMF_PS1_SV1.c (added)

branches/eam_branches/psModules.20100506/src/objects/Makefile.am

@@ -43 +43 @@
         pmSourceIO_CMF_PS1_V1.c \
         pmSourceIO_CMF_PS1_V2.c \
+        pmSourceIO_CMF_PS1_SV1.c \
         pmSourceIO_CMF_PS1_DV1.c \
         pmSourceIO_MatchedRefs.c \

branches/eam_branches/psModules.20100506/src/objects/pmSource.c

@@ -54 +54 @@
     psFree(tmp->moments);
     psFree(tmp->diffStats);
+    psFree(tmp->radial);
     psFree(tmp->blends);
     psTrace("psModules.objects", 10, "---- end ----\n");
@@ -116 +117 @@
     source->extpars = NULL;
     source->diffStats = NULL;
+    source->radial = NULL;
 
     source->region = psRegionSet(NAN, NAN, NAN, NAN);

branches/eam_branches/psModules.20100506/src/objects/pmSource.h

@@ -91 +91 @@
     pmSourceExtendedPars *extpars;      ///< extended source parameters
     pmSourceDiffStats *diffStats;       ///< extra parameters for difference detections
+    pmSourceRadialApertures *radial;    ///< radial flux in circular apertures
     int imageID;
 };

branches/eam_branches/psModules.20100506/src/objects/pmSourceExtendedPars.c

@@ -65 +65 @@
 }
 
+// pmSourceRadialApertures carries the raw radial flux information, including angular bins
+static void pmSourceRadialAperturesFree(pmSourceRadialApertures *radial)
+{
+    if (!radial) return;
+    psFree(radial->flux);
+    psFree(radial->fill);
+}
+
+pmSourceRadialApertures *pmSourceRadialAperturesAlloc()
+{
+    pmSourceRadialApertures *radial = (pmSourceRadialApertures *)psAlloc(sizeof(pmSourceRadialApertures));
+    psMemSetDeallocator(radial, (psFreeFunc) pmSourceRadialAperturesFree);
+
+    radial->flux = NULL;
+    radial->fill = NULL;
+    return radial;
+}
+
+bool psMemCheckSourceRadialApertures(psPtr ptr)
+{
+    PS_ASSERT_PTR(ptr, false);
+    return ( psMemGetDeallocator(ptr) == (psFreeFunc) pmSourceRadialAperturesFree);
+}
+
 // pmSourceEllipticalFlux carries the elliptical renormalized radial flux info
 static void pmSourceEllipticalFluxFree(pmSourceEllipticalFlux *flux)

branches/eam_branches/psModules.20100506/src/objects/pmSourceExtendedPars.h

@@ -20 +20 @@
     psVector *isophotalRadii;           // isophotal radius for the above angles
 } pmSourceRadialFlux;
+
+typedef struct {
+    psVector *flux;                     // fluxes measured at above radii
+    psVector *fill;                     // angles corresponding to above radial profiles
+} pmSourceRadialApertures;
 
 typedef struct {
@@ -62 +67 @@
 bool psMemCheckSourceRadialFlux(psPtr ptr);
 
+pmSourceRadialApertures *pmSourceRadialAperturesAlloc();
+bool psMemCheckSourceRadialApertures(psPtr ptr);
+
 pmSourceEllipticalFlux *pmSourceEllipticalFluxAlloc();
 bool psMemCheckSourceEllipticalFlux(psPtr ptr);

branches/eam_branches/psModules.20100506/src/objects/pmSourceIO.c

@@ -539 +539 @@
                 status &= pmSourcesWrite_CMF_PS1_V2 (file->fits, readout, sources, file->header, outhead, dataname);
             }
+            if (!strcmp (exttype, "PS1_SV1")) {
+                status &= pmSourcesWrite_CMF_PS1_SV1 (file->fits, readout, sources, file->header, outhead, dataname, recipe);
+            }
             if (!strcmp (exttype, "PS1_DV1")) {
                 status &= pmSourcesWrite_CMF_PS1_DV1 (file->fits, readout, sources, file->header, outhead, dataname);
@@ -556 +559 @@
                     status &= pmSourcesWrite_CMF_PS1_V2_XSRC (file->fits, readout, sources, file->header, xsrcname, recipe);
                 }
+                if (!strcmp (exttype, "PS1_SV1")) {
+                    status &= pmSourcesWrite_CMF_PS1_SV1_XSRC (file->fits, readout, sources, file->header, xsrcname, recipe);
+                }
                 if (!strcmp (exttype, "PS1_DV1")) {
                     status &= pmSourcesWrite_CMF_PS1_DV1_XSRC (file->fits, readout, sources, file->header, xsrcname, recipe);
@@ -573 +579 @@
                     status &= pmSourcesWrite_CMF_PS1_V2_XFIT (file->fits, readout, sources, xfitname);
                 }
+                if (!strcmp (exttype, "PS1_SV1")) {
+                    status &= pmSourcesWrite_CMF_PS1_SV1_XFIT (file->fits, readout, sources, xfitname);
+                }
                 if (!strcmp (exttype, "PS1_DV1")) {
                     status &= pmSourcesWrite_CMF_PS1_DV1_XFIT (file->fits, readout, sources, xfitname);
@@ -981 +990 @@
         // we need to find the corresponding table EXTNAME.
         // first check the header
-        char *extdata = psMetadataLookupStr (NULL, hdu->header, "EXTDATA");
+        char *extdata = psMetadataLookupStr (&status, hdu->header, "EXTDATA");
         if (extdata) {
             // if EXTDATA is defined in the header, use that value for 'dataname'
@@ -1023 +1032 @@
             if (!strcmp (exttype, "PS1_V2")) {
                 sources = pmSourcesRead_CMF_PS1_V2 (file->fits, hdu->header);
+            }
+            if (!strcmp (exttype, "PS1_SV1")) {
+                sources = pmSourcesRead_CMF_PS1_SV1 (file->fits, hdu->header);
             }
             if (!strcmp (exttype, "PS1_DV1")) {

branches/eam_branches/psModules.20100506/src/objects/pmSourceIO.h

@@ -43 +43 @@
 bool pmSourcesWrite_CMF_PS1_V2_XFIT (psFits *fits, pmReadout *readout, psArray *sources, char *extname);
 
+bool pmSourcesWrite_CMF_PS1_SV1 (psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, psMetadata *tableHeader, char *extname, psMetadata *recipe);
+bool pmSourcesWrite_CMF_PS1_SV1_XSRC (psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, char *extname, psMetadata *recipe);
+bool pmSourcesWrite_CMF_PS1_SV1_XFIT (psFits *fits, pmReadout *readout, psArray *sources, char *extname);
+
 bool pmSourcesWrite_CMF_PS1_DV1 (psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, psMetadata *tableHeader, char *extname);
 bool pmSourcesWrite_CMF_PS1_DV1_XSRC (psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, char *extname, psMetadata *recipe);
@@ -57 +61 @@
 psArray *pmSourcesRead_CMF_PS1_V1 (psFits *fits, psMetadata *header);
 psArray *pmSourcesRead_CMF_PS1_V2 (psFits *fits, psMetadata *header);
+psArray *pmSourcesRead_CMF_PS1_SV1 (psFits *fits, psMetadata *header);
 psArray *pmSourcesRead_CMF_PS1_DV1 (psFits *fits, psMetadata *header);