Changeset 10954 for trunk/pswarp/src/pswarpDataLoad.c

Timestamp:

Jan 7, 2007, 6:48:50 AM (19 years ago)

Author:

eugene

Message:

functioning pswarp, needs opt code

File:

• trunk/pswarp/src/pswarpDataLoad.c (modified) (3 diffs)

trunk/pswarp/src/pswarpDataLoad.c

@@ -32 +32 @@
 
     // de-activate PSWARP.SKYCELL and PSWARP.OUTPUT
-    // 
+    pmFPAfileActivate (config->files, false, "PSWARP.SKYCELL");
     pmFPAfileActivate (config->files, false, "PSWARP.OUTPUT");
-    pmFPAfileActivate (config->files, false, "PSWARP.SKYCELL");
 
     pmFPAview *view = pmFPAviewAlloc (0);
+
+    // find the FPA phu
+    // XXX wrap the test below into a function
+    // XXX need to read only the headers for the skycell
+    // XXX need to optionally load the astrometry datafile
+    bool bilevelAstrometry = false;
+    pmHDU *phu = pmFPAviewThisPHU (view, input->fpa);
+    if (phu) {
+      char *ctype = psMetadataLookupStr (NULL, phu->header, "CTYPE1");
+      if (ctype) {
+        bilevelAstrometry = !strcmp (&ctype[4], "-DIS");
+      }
+    }
+    if (bilevelAstrometry) {
+      pmAstromReadBilevelMosaic (input->fpa, phu->header);
+    } 
 
     // files associated with the science image
@@ -45 +60 @@
         if (!chip->process || !chip->file_exists) { continue; }
         pmFPAfileIOChecks (config, view, PM_FPA_BEFORE);
+
+        // read WCS data from the corresponding header
+        pmHDU *hdu = pmFPAviewThisHDU (view, input->fpa);
+        if (bilevelAstrometry) {
+          pmAstromReadBilevelChip (chip, hdu->header); 
+        } else {
+          // XXX get pixelScale from recipes.  does it matter?
+          float pixelScale = 13.5;
+          pmAstromReadWCS (input->fpa, chip, hdu->header, pixelScale);
+        }
 
         while ((cell = pmFPAviewNextCell (view, input->fpa, 1)) != NULL) {
@@ -56 +81 @@
                 if (! readout->data_exists) { continue; }
 
-                pswarpConvertReadout (output, readout, config);
+                pswarpTransformReadout (output, readout, config);
 
                 pmFPAfileIOChecks (config, view, PM_FPA_AFTER);