Changeset 18435

Timestamp:

Jul 8, 2008, 9:51:39 AM (18 years ago)

Author:

Paul Price

Message:

Propagate WCS.

Location:

trunk/ppStack/src

Files:

• ppStackLoop.c (modified) (2 diffs)
• ppStackReadout.c (modified) (2 diffs)

trunk/ppStack/src/ppStackLoop.c

@@ -13 +13 @@
 
 //#define TESTING
+
+#define WCS_TOLERANCE 0.001             // Tolerance for WCS
 
 // Here follows lists of files for activation/deactivation at various stages.  Each must be NULL-terminated.
@@ -702 +704 @@
 
         // Close up
+        bool wcsDone = false;           // Have we done the WCS?
         for (int i = 0; i < num; i++) {
             if (!readouts->data[i]) {
                 continue;
             }
+            if (!wcsDone) {
+                // Copy astrometry over
+                wcsDone = true;
+                pmReadout *inRO = readouts->data[i]; // Template readout
+                pmHDU *inHDU = pmHDUFromCell(inRO->parent); // Template HDU
+                pmHDU *outHDU = pmHDUFromCell(outRO->parent); // Output HDU
+                pmChip *outChip = outRO->parent->parent; // Output chip
+                pmFPA *outFPA = outChip->parent; // Output FPA
+                if (!outHDU || !inHDU) {
+                    psWarning("Unable to find HDU at FPA level to copy astrometry.");
+                } else {
+                    if (!pmAstromReadWCS(outFPA, outChip, inHDU->header, 1.0)) {
+                        psErrorClear();
+                        psWarning("Unable to read WCS astrometry from input FPA.");
+                        wcsDone = false;
+                    } else {
+                        if (!outHDU->header) {
+                            outHDU->header = psMetadataAlloc();
+                        }
+                        if (!pmAstromWriteWCS(outHDU->header, outFPA, outChip, WCS_TOLERANCE)) {
+                            psErrorClear();
+                            psWarning("Unable to write WCS astrometry to output FPA.");
+                            wcsDone = false;
+                        }
+                    }
+                }
+            }
+
             psFitsClose(imageFits->data[i]);
             psFitsClose(maskFits->data[i]);

trunk/ppStack/src/ppStackReadout.c

@@ -9 +9 @@
 
 #include "ppStack.h"
-
-#define WCS_TOLERANCE 0.001             // Tolerance for WCS
 
 //#define REJECTION_FILES                 // Write rejection mask?
@@ -185 +183 @@
         }
         totExposure += exposure;        // Total exposure time
-
-#if 0
-        if (i == 0) {
-            // Copy astrometry over
-            pmHDU *hdu = fpa->hdu; // Template HDU
-            pmHDU *outHDU = outFPA->hdu; // Output HDU
-            if (!outHDU || !hdu) {
-                psWarning("Unable to find HDU at FPA level to copy astrometry.");
-            } else {
-                if (!pmAstromReadWCS(outFPA, outChip, hdu->header, 1.0)) {
-                    psErrorClear();
-                    psWarning("Unable to read WCS astrometry from input FPA.");
-                } else {
-                    if (!outHDU->header) {
-                        outHDU->header = psMetadataAlloc();
-                    }
-                    if (!pmAstromWriteWCS(outHDU->header, outFPA, outChip, WCS_TOLERANCE)) {
-                        psErrorClear();
-                        psWarning("Unable to write WCS astrometry to output FPA.");
-                    }
-                }
-            }
-        }
-#endif
 
         // Ensure there is a mask, or pmStackCombine will complain