Changeset 30624 for trunk/psphot/src/psphotStackImageLoop.c

Timestamp:

Feb 13, 2011, 12:33:05 PM (15 years ago)

Author:

eugene

Message:

adding some of the metadata needed by PSPS to output headers; skip models with few valid pixels (eg, only the edge is showing); only do the linear fit on pixels within the fit radius; modify psf-match auto-scaling process; enable multiple target psfs for matched-psf aperture photometry; speed up analysis of the radial apertures

Location:

trunk/psphot

Files:

• . (modified) (1 prop)
• src/psphotStackImageLoop.c (modified) (5 diffs)

trunk/psphot/src/psphotStackImageLoop.c

@@ -101 +101 @@
 */
 
-# define UPDATE_HEADER 0
+# define UPDATE_HEADER 1
 
 bool GetAstrometryFPA (pmConfig *config, pmFPAview *view) {
@@ -174 +174 @@
         pmChip *outChip = pmFPAviewThisChip(view, output->fpa); ///< Chip in the output
 
-# if (UPDATE_HEADER)
         pmHDU *outHDU = pmFPAviewThisHDU (view, output->fpa);
         if (!outHDU) {
-            pmFPAAddSourceFromView(output->fpa, "name", view, output->format);
+            pmFPAAddSourceFromView(output->fpa, view, output->format);
             outHDU = pmFPAviewThisHDU (view, output->fpa);
             psAssert (outHDU, "failed to make HDU");
         }
-# endif
+        if (!outHDU->header) {
+          outHDU->header = psMetadataAlloc();
+        }
 
         if (bilevelAstrometry) {
@@ -188 +189 @@
                 continue;
             }
-# if (UPDATE_HEADER)
             if (!pmAstromWriteBilevelChip(outHDU->header, outChip, WCS_NONLIN_TOL)) {
                 psWarning("Unable to generate WCS header.");
                 continue;
             }
-# endif
         } else {
             // we use a default FPA pixel scale of 1.0
@@ -200 +199 @@
                 continue;
             }
-# if (UPDATE_HEADER)
-            if (UPDATE_HEADER && !pmAstromWriteWCS(outHDU->header, output->fpa, outChip, WCS_NONLIN_TOL)) {
+            if (!pmAstromWriteWCS(outHDU->header, output->fpa, outChip, WCS_NONLIN_TOL)) {
                 psWarning("Unable to generate WCS header.");
                 continue;
             }
-# endif
         }
     }
@@ -225 +222 @@
         psAssert (output, "missing file?");
 
-# if (UPDATE_HEADER)
         pmHDU *PHU = pmFPAviewThisPHU(view, output->fpa);
         if (!PHU) {
-            pmFPAAddSourceFromView(output->fpa, "name", view, output->format);
+            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.");
         }
-# endif
     }