Changeset 13189 for trunk/psModules/src/camera/pmFPAMosaic.c

Timestamp:

May 2, 2007, 8:07:43 PM (19 years ago)

Author:

magnier

Message:

added ifdef-ed out code to set the WCS for a chip in pmChipMosaic

File:

• trunk/psModules/src/camera/pmFPAMosaic.c (modified) (2 diffs)

trunk/psModules/src/camera/pmFPAMosaic.c

@@ -1072 +1072 @@
         item->data.S32 = yBin;
     }
-
     psFree(sourceCells);
     psFree(chipRegion);
@@ -1117 +1116 @@
         return false;
     }
+    
+// this code extracts the WCS information from the header of one of the cells and applies it to
+// the mosaicked chip.  XXX EAM : I am leaving it at this point for PAP to get the binning, etc
+// from the concepts
+# if (0)
+    // if the cells contain the headers, we need to apply the WCS terms from (one of?) the cells
+    for (int i = 0; i < source->cells->n; i++) {
+        pmCell *cell = source->cells->data[i];
+        if (cell == NULL) continue;
+        if (cell->hdu == NULL) continue;
+
+        // extract the WCS terms from this cell
+        pmAstromWCS *wcs = pmAstromWCSfromHeader (cell->hdu->header);
+
+        // modify the wcs terms for the cell offset, binning, and parity
+        xBinRatio = xBinCell / xBinChip;
+        if (xFlip) {
+            wcs->crpix1 = xCell - wcs->crpix1*xBinRatio;
+            wcs->cdelt1 *= -1;
+        } else {
+            wcs->crpix1 = xCell + wcs->crpix1*xBinRatio;
+        }
+        wcs->cdelt1 *= xBinRatio;
+
+        yBinRatio = yBinCell / yBinChip;
+        if (yFlip) {
+            wcs->crpix2 = yCell - wcs->crpix2*yBinRatio;
+            wcs->cdelt2 *= -1;
+        } else {
+            wcs->crpix2 = xCell + wcs->crpix2*yBinRatio;
+        }
+        wcs->cdelt2 *= yBinRatio;
+
+        pmAstromWCStoHeader (targetHDU->header);
+        break;
+        // XXX rather than quitting at this point, we could save this wcs structure and compare
+        // its values to the equivalent version from one of the other cells.
+    }
+# endif
 
     return true;