Changeset 4694 for trunk/archive/scripts/src/test_pmFPARead.c

Timestamp:

Aug 1, 2005, 2:21:14 PM (21 years ago)

Author:

Paul Price

Message:

Rough draft of pmFPAMorph working

File:

• trunk/archive/scripts/src/test_pmFPARead.c (modified) (2 diffs)

trunk/archive/scripts/src/test_pmFPARead.c

@@ -26 +26 @@
 
     (void)psTraceSetLevel(".", 0);
+    (void)psTraceSetLevel("readExtension", 10);
     (void)psTraceSetLevel("readMultipleRegions", 0);
     (void)psTraceSetLevel("portions", 0);
-    (void)psTraceSetLevel("pmFPARead", 10);
+    (void)psTraceSetLevel("pmFPAPrint", 10);
     (void)psTraceSetLevel("pmFPAfromHeader", 10);
     (void)psTraceSetLevel("pmCameraFromHeader", 10);
-
+    (void)psTraceSetLevel("pmFPAMorph", 10);
+    (void)psTraceSetLevel("spliceCells", 10);
 
 
@@ -53 +55 @@
     psMetadata *camera = pmCameraFromHeader(header, ipprc);
 
-    papFPA *fpa = pmFPARead(fits, camera, NULL);
+    papFPA *fpa = pmFPAConstruct(camera, NULL);
+
+    // Cut off a lot of the chips so I can fit things in memory
+    {
+        psArray *chips = fpa->chips;
+        papChip *chip = chips->data[0];
+        psFree(chips);
+        chips = psArrayAlloc(1);
+        chips->data[0] = chip;
+    }
+
+    (void)pmFPARead(fpa, fits);
     (void)pmFPAPrint(fpa);
 
-    printf("Test: %d %d\n", psCellGetValueS32(fpa, "ccd15", "left", "CELL.XPARITY"),
-           psCellGetValueS32(fpa, "ccd15", "right", "CELL.XPARITY"));
-    printf("Test: %d %d\n", psCellGetValueS32(fpa, "ccd15", "left", "CELL.YPARITY"),
-           psCellGetValueS32(fpa, "ccd15", "right", "CELL.YPARITY"));
-    printf("Test: %d %d\n", psCellGetValueS32(fpa, "ccd15", "left", "CELL.YPARITY"),
-           psCellGetValueS32(fpa, "ccd15", "right", "CELL.YPARITY"));
+    // Morph to a splice image
+    psMetadata *newCamera = psMetadataConfigParse(NULL, &badLines, "megacam_splice.config", true);
+    papFPA *newfpa = pmFPAConstruct(newCamera, NULL);
+
+    // Cut off a lot of the chips so I can fit things in memory
+    {
+        psArray *chips = newfpa->chips;
+        papChip *chip = chips->data[0];
+        psFree(chips);
+        chips = psArrayAlloc(1);
+        chips->data[0] = chip;
+    }
+
+    pmFPAPrint(newfpa);
+    pmFPAMorph(newfpa, fpa, 0, 0);
+    pmFPAPrint(newfpa);
+
+#if 0
+    // Write out
+    psFits *newfits = psFitsAlloc("test.fits");
+    pmFPAWrite(fits, newfpa);
+    psFree(newfits);
+#endif
+
+#if 1
+    psArray *chips = newfpa->chips;
+    for (int i = 0; i < chips->n; i++) {
+        papChip *chip = chips->data[i]; // Chip of interest
+        psArray *cells = chip->cells;   // Array of cells
+
+        if (chip->extname) {
+            psFits *fp = psFitsAlloc(chip->extname);
+            psFitsWriteImage(fp, chip->header, chip->pixels->data[i], i);
+            psFree(fp);
+        }
+
+        for (int j = 0; j < cells->n; j++) {
+            papCell *cell = cells->data[j]; // Cell of interest
+
+            if (cell->extname) {
+                psFits *fp = psFitsAlloc(cell->extname);
+                psFitsWriteImage(fp, cell->header, cell->pixels->data[i], i);
+                psFree(fp);
+            }
+        }
+    }
+#endif
+
+    psFree(newfpa);
+    psFree(newCamera);
 
     // Tidy up