Changeset 6080 for branches/eam_rel9_p0/psModules/src/astrom/pmFPA.c

Timestamp:

Jan 19, 2006, 11:47:06 PM (21 years ago)

Author:

Paul Price

Message:

Readouts now contain a subimage

File:

• branches/eam_rel9_p0/psModules/src/astrom/pmFPA.c (modified) (8 diffs)

branches/eam_rel9_p0/psModules/src/astrom/pmFPA.c

@@ -12 +12 @@
 * XXX: Should we implement non-linear cell->chip transforms?
 *
-*  @version $Revision: 1.1.2.1 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2006-01-20 02:36:41 $
+*  @version $Revision: 1.1.2.2 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2006-01-20 09:47:06 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -23 +23 @@
 #include <string.h>
 #include <math.h>
+#include <assert.h>
 #include "pslib.h"
 
 #include "pmFPA.h"
 #include "pmConcepts.h"
+#include "pmMaskBadPixels.h"
 
 /******************************************************************************
@@ -163 +165 @@
     tmpReadout->mask = NULL;
     tmpReadout->weight = NULL;
+    tmpReadout->bias = psListAlloc(NULL);
     tmpReadout->analysis = psMetadataAlloc();
     tmpReadout->parent = cell;
@@ -390 +393 @@
 }
 
+bool pmChipSelectCell(pmChip *chip,
+                      int cellNum
+                     )
+{
+    assert(chip);
+
+    psArray *cells = chip->cells;       // Component cells
+    if (!cells || cellNum > cells->n) {
+        return false;
+    }
+
+    for (int i = 0; i < cells->n; i++) {
+        pmCell *cell = cells->data[i];
+        if (!cell) {
+            continue;
+        }
+        cell->process = (i == cellNum);
+    }
+
+    return true;
+}
 
 /*****************************************************************************
@@ -437 +461 @@
 }
 
+int pmChipExcludeCell(pmChip *chip,
+                      int cellNum
+                     )
+{
+    assert(chip);
+
+    psArray *cells = chip->cells;       // The component cells
+    if (!cells || cellNum > cells->n) {
+        return 0;
+    }
+
+    int numCells = 0;                   // Number of cells to be processed
+    for (int i = 0; i < cells->n; i++) {
+        pmCell *cell = cells->data[i];
+        if (!cell) {
+            continue;
+        }
+        if (i == cellNum) {
+            cell->process = false;
+        } else {
+            numCells++;
+        }
+    }
+
+    return numCells;
+}
+
 
 bool pmCellSetWeights(pmCell *cell // Cell for which to set weights
@@ -444 +495 @@
     float readnoise = psMetadataLookupF32(NULL, cell->concepts, "CELL.READNOISE"); // Cell read noise
     psRegion *trimsec = psMetadataLookupPtr(NULL, cell->concepts, "CELL.TRIMSEC"); // Trim section
+    float saturation = psMetadataLookupF32(NULL, cell->concepts, "CELL.SATURATION"); // Saturation level
+    float bad = psMetadataLookupF32(NULL, cell->concepts, "CELL.BAD"); // Bad level
 
     p_pmHDU *hdu = cell->hdu;           // The data unit, containing the weight and mask originals
@@ -478 +531 @@
             psImage *image = pixels->data[i];
             masks->data[i] = psImageAlloc(image->numCols, image->numRows, PS_TYPE_U8);
-            psImageInit(masks->data[i], 0);
         }
         hdu->masks = masks;
@@ -489 +541 @@
 
         if (! readout->weight) {
-            readout->weight = weights->data[i];
+            readout->weight = psImageSubset(weights->data[i], *trimsec);
         }
         if (! readout->mask) {
-            readout->mask = masks->data[i];
-        }
-
-        // Mask is already set to 0
+            readout->mask = psImageSubset(masks->data[i], *trimsec);
+        }
+
+        // Set up the mask
+        psImage *image = readout->image;// Pixels
+        psImage *mask = readout->mask;  // Mask image
+        for (int i = 0; i < image->numRows; i++) {
+            for (int j = 0; j < image->numCols; j++) {
+                if (image->data.F32[i][j] > saturation) {
+                    mask->data.F32[i][j] = PM_MASK_SAT;
+                }
+                if (image->data.F32[i][j] < bad) {
+                    mask->data.F32[i][j] = PM_MASK_BAD;
+                }
+            }
+        }
 
         // Set weight image to the variance = g*f + rn^2
-        psImage *image = psImageSubset(readout->image, *trimsec); // The pixels
-        psImage *weight = psImageSubset(readout->weight, *trimsec); // The weight map
-        psBinaryOp(weight, image, "*", psScalarAlloc(gain, PS_TYPE_F32));
-        psBinaryOp(weight, weight, "+", psScalarAlloc(readnoise*readnoise, PS_TYPE_F32));
+        psBinaryOp(readout->weight, image, "*", psScalarAlloc(gain, PS_TYPE_F32));
+        psBinaryOp(readout->weight, readout->weight, "+", psScalarAlloc(readnoise*readnoise, PS_TYPE_F32));
     }