Changeset 30871

Timestamp:

Mar 10, 2011, 5:31:06 PM (15 years ago)

Author:

eugene

Message:

propagate the header information to the output files in the psphotStack context

Location:

branches/eam_branches/ipp-20110213/psphot/src

Files:

• psphotDeblendSatstars.c (modified) (1 diff)
• psphotStackImageLoop.c (modified) (8 diffs)
• psphotStackMatchPSFs.c (modified) (1 diff)

branches/eam_branches/ipp-20110213/psphot/src/psphotDeblendSatstars.c

@@ -86 +86 @@
         // XXX filter? if (source->mode & PM_SOURCE_MODE_SATSTAR) continue;
         if (source->mode & PM_SOURCE_MODE_BLEND) continue;
-        if (source->peak->flux < SATURATION) continue;
+        if (source->peak->flux < SAT_TEST_LEVEL) continue;
 
         // save these for reference below

branches/eam_branches/ipp-20110213/psphot/src/psphotStackImageLoop.c

@@ -1 +1 @@
 # include "psphotStandAlone.h"
-#define WCS_NONLIN_TOL 0.001            // Non-linear tolerance for header WCS
 
 # define ESCAPE(MESSAGE) {                              \
@@ -8 +7 @@
     }
 
-bool GetAstrometryFPA (pmConfig *config, pmFPAview *view);
-bool SetAstrometryFPA (pmConfig *config, pmFPAview *view, bool bilevelAstrometry);
-bool GetAstrometryChip (pmConfig *config, pmFPAview *view, bool bilevelAstrometry);
+// XXX this implementation is not smart about multi-level astrometry headers
+bool UpdateHeadersForFPA (pmConfig *config, pmFPAview *view);
+bool UpdateHeadersForChip (pmConfig *config, pmFPAview *view);
 
 bool psphotStackImageLoop (pmConfig *config) {
@@ -34 +33 @@
     // XXX for now, just load the full set of images up front
     if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE ("failed input for fpa in psphot.");
-
-    bool bilevelAstrometry = GetAstrometryFPA(config, view);
 
     // for psphot, we force data to be read at the chip level
@@ -83 +80 @@
         }
 
-        GetAstrometryChip(config, view, bilevelAstrometry);
+        UpdateHeadersForChip(config, view);
 
         // save output which is saved at the chip level
@@ -90 +87 @@
     psMemDump("doneloop");
 
-    SetAstrometryFPA(config, view, bilevelAstrometry);
-    
+    UpdateHeadersForFPA(config, view);
+
     // save output which is saved at the fpa level
     if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE ("failed ouput for FPA in psphot.");
@@ -104 +101 @@
 }
 
-/* 
-   the easiest way to implement this is to assume we can pre-load the full set of images up front.
-   with 5 filters and 6000^2 (image, mask, var = 10 byte per pixel), we need 1.8GB, which is not too bad.
-*/
-
-# define UPDATE_HEADER 1
-
-bool GetAstrometryFPA (pmConfig *config, pmFPAview *view) {
-
-    bool foundAstrometry = false;
-    bool bilevelAstrometry = false;
+bool UpdateHeadersForChip (pmConfig *config, pmFPAview *view) {
 
     int num = psphotFileruleCount(config, "PSPHOT.INPUT");
@@ -130 +117 @@
         psAssert (input, "missing input file");
 
-        // find the FPA phu
-        pmHDU *phu = pmFPAviewThisPHU(view, input->fpa);
-        if (!phu) {
-            psWarning("Unable to read bilevel mosaic astrometry for input FPA entry %d", i);
-            continue;
+        // just copy the input headers to the output headers, then update version info
+        pmChip *inChip = pmFPAviewThisChip(view, input->fpa); ///< Chip in the input
+        pmHDU *inHDU = pmFPAviewThisHDU (view, input->fpa);
+
+        pmChip *outChip = pmFPAviewThisChip(view, output->fpa); ///< Chip in the output
+        pmHDU *outHDU = pmFPAviewThisHDU (view, output->fpa);
+        if (!outHDU) {
+            pmFPAAddSourceFromView(output->fpa, view, output->format);
+            outHDU = pmFPAviewThisHDU (view, output->fpa);
+            psAssert (outHDU, "failed to make HDU");
         }
-
-        char *ctype = psMetadataLookupStr(NULL, phu->header, "CTYPE1");
-        if (!ctype) {
-            psWarning("Error in WCS keywords for input FPA entry %d", i);
-            continue;
+        if (!outHDU->header) {
+            outHDU->header = psMetadataCopy(NULL, inHDU->header);
         }
-
-        if (!foundAstrometry) {
-            bilevelAstrometry = !strcmp (&ctype[4], "-DIS");
-        } else {
-            if (bilevelAstrometry != !strcmp (&ctype[4], "-DIS")) {
-                psAbort("astrometry type mis-match");
-            }
-        }
-
-        if (bilevelAstrometry) {
-            // update the output structures
-            if (!pmAstromReadBilevelMosaic(output->fpa, phu->header)) {
-                psWarning("Unable to read bilevel mosaic astrometry for input FPA.");
-            }
-        }
+        outChip->toFPA = psMemIncrRefCounter(inChip->toFPA);
+        outChip->fromFPA = psMemIncrRefCounter(inChip->fromFPA);
     }
-    return bilevelAstrometry;
+    return true;
 }
 
-bool GetAstrometryChip (pmConfig *config, pmFPAview *view, bool bilevelAstrometry) {
+bool UpdateHeadersForFPA (pmConfig *config, pmFPAview *view) {
 
     int num = psphotFileruleCount(config, "PSPHOT.INPUT");
@@ -177 +153 @@
         psAssert (input, "missing input file");
 
-        // Need to update the output for astrometry.  Read WCS data from the corresponding
-        // header and save in the output fpa
-        pmHDU *inHDU = pmFPAviewThisHDU (view, input->fpa);
-        pmChip *outChip = pmFPAviewThisChip(view, output->fpa); ///< Chip in the output
+        output->fpa->toTPA = psMemIncrRefCounter(input->fpa->toTPA);
+        output->fpa->fromTPA = psMemIncrRefCounter(input->fpa->fromTPA);
+        output->fpa->toSky = psMemIncrRefCounter(input->fpa->toSky);
 
-        pmHDU *outHDU = pmFPAviewThisHDU (view, output->fpa);
-        if (!outHDU) {
-            pmFPAAddSourceFromView(output->fpa, view, output->format);
-            outHDU = pmFPAviewThisHDU (view, output->fpa);
-            psAssert (outHDU, "failed to make HDU");
-        }
-        if (!outHDU->header) {
-          outHDU->header = psMetadataAlloc();
-        }
+        pmConceptsCopyFPA(output->fpa, input->fpa, true, true);
 
-        if (bilevelAstrometry) {
-            if (!pmAstromReadBilevelChip (outChip, inHDU->header)) {
-                psWarning("Unable to read bilevel chip astrometry for input FPA.");
-                continue;
-            }
-            if (!pmAstromWriteBilevelChip(outHDU->header, outChip, WCS_NONLIN_TOL)) {
-                psWarning("Unable to generate WCS header.");
-                continue;
-            }
-        } else {
-            // we use a default FPA pixel scale of 1.0
-            if (!pmAstromReadWCS (output->fpa, outChip, inHDU->header, 1.0)) {
-                psWarning("Unable to read WCS astrometry for input FPA.");
-                continue;
-            }
-            if (!pmAstromWriteWCS(outHDU->header, output->fpa, outChip, WCS_NONLIN_TOL)) {
-                psWarning("Unable to generate WCS header.");
-                continue;
-            }
-        }
+        // XXX TEST
+        // pmFPASetFileStatus(output->fpa, true);
+        // pmFPASetDataStatus(output->fpa, true);
+        // pmChip *chip = output->fpa->chips->data[0];
+        // pmCell *cell = chip->cells->data[0];
+        // pmReadout *readout = cell->readouts->data[0];
     }
-
     return true;
 }
 
-bool SetAstrometryFPA (pmConfig *config, pmFPAview *view, bool bilevelAstrometry) {
+/* 
+   the easiest way to implement this is to assume we can pre-load the full set of images up front.
+   with 5 filters and 6000^2 (image, mask, var = 10 byte per pixel), we need 1.8GB, which is not too bad.
+*/
 
-    if (!bilevelAstrometry) return true;
-
-    int num = psphotFileruleCount(config, "PSPHOT.INPUT");
-
-    // loop over the available readouts
-    for (int i = 0; i < num; i++) {
-
-        // find the currently selected readout
-        pmFPAfile *output = pmFPAfileSelectSingle(config->files, "PSPHOT.STACK.OUTPUT.IMAGE", i); // File of interest
-        psAssert (output, "missing file?");
-
-        pmHDU *PHU = pmFPAviewThisPHU(view, output->fpa);
-        if (!PHU) {
-            pmFPAAddSourceFromView(output->fpa, view, output->format);
-            PHU = pmFPAviewThisPHU (view, output->fpa);
-            psAssert (PHU, "failed to make PHU");
-        }
-        if (!PHU->header) {
-          PHU->header = psMetadataAlloc();
-        }
-
-        if (!pmAstromWriteBilevelMosaic(PHU->header, output->fpa, WCS_NONLIN_TOL)) {
-            psWarning("Unable to generate WCS header.");
-        }
-    }
-
-    return true;
-}

branches/eam_branches/ipp-20110213/psphot/src/psphotStackMatchPSFs.c

@@ -92 +92 @@
     }
 
+    // copy the header data from Src to Out
+    // pmHDU *hduSrc = pmHDUFromReadout(readoutSrc);
+    // psAssert (hduSrc, "input missing hdu?");
+
+
     // set NAN pixels to 'SAT'
     psImageMaskType maskSat = pmConfigMaskGet("SAT", config);