Changeset 24154 for branches/pap_magic/ppSub/src/ppSubReadoutInverse.c

Timestamp:

May 12, 2009, 6:15:31 PM (17 years ago)

Author:

Paul Price

Message:

Generate WCS for inverse images too!

File:

• branches/pap_magic/ppSub/src/ppSubReadoutInverse.c (modified) (2 diffs)

branches/pap_magic/ppSub/src/ppSubReadoutInverse.c

@@ -7 +7 @@
 
 #include "ppSub.h"
-
 
 bool ppSubReadoutInverse(pmConfig *config)
@@ -24 +23 @@
     invRO->data_exists = invRO->parent->data_exists = invRO->parent->parent->data_exists = true;
 
-    pmConceptsCopyFPA(invRO->parent->parent->parent, outRO->parent->parent->parent, true, true);
+    pmChip *outChip = outRO->parent->parent;       // Output chip
+    pmFPA *outFPA = outChip->parent;               // Output FPA
+    pmFPA *invFPA = invRO->parent->parent->parent; // Inverse FPA
+
+    pmConceptsCopyFPA(invFPA, outFPA, true, true);
+
+#if 0
+    // Copy astrometry over
+    pmHDU *outHDU = outFPA->hdu;          // Output HDU
+    pmChip *invChip = invRO->parent->parent; // Inverse chip
+    psFree(view);
+    if (!outHDU || !inHDU) {
+        psError(PS_ERR_UNKNOWN, false, "Unable to find HDU at FPA level to copy astrometry.");
+        return false;
+    }
+    if (!pmAstromReadWCS(invFPA, invChip, outHDU->header, 1.0)) {
+        psError(PS_ERR_UNKNOWN, false, "Unable to read WCS astrometry from PPSUB.OUTPUT.");
+        return false;
+    }
+#endif
+    pmHDU *invHDU = invFPA->hdu;          // Inverse HDU
+    if (!pmAstromWriteWCS(invHDU->header, outFPA, outChip, WCS_TOLERANCE)) {
+        psError(PS_ERR_UNKNOWN, false, "Unable to write WCS astrometry to PPSUB.INVERSE.");
+        return false;
+    }
 
     return true;