Changeset 41851 for branches/eam_branches/ipp-dev-20210817

Timestamp:

Oct 19, 2021, 4:53:35 PM (5 years ago)

Author:

eugene

Message:

adding pmPattern amplitude I/O

Location:

branches/eam_branches/ipp-dev-20210817/psModules/src/detrend

Files:

• pmPattern.c (modified) (1 diff)
• pmPatternIO.c (modified) (1 diff)

branches/eam_branches/ipp-dev-20210817/psModules/src/detrend/pmPattern.c

@@ -476 +476 @@
     psImage *mask = ro->mask;           // Mask for image
     int numCols = image->numCols, numRows = image->numRows; // Size of image
+
+    // use the background and typical amplitude value to choose to correct or not
+    // pmCell *cell = ro->parent;
 
     psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);

branches/eam_branches/ipp-dev-20210817/psModules/src/detrend/pmPatternIO.c

@@ -452 +452 @@
     return pmReadoutReadPattern(readout, file->fits);
 }
+
+/**************** PatternRowAmp(litude) I/O *************************/
+
+# if (0)
+/********************* Read Data functions *****************************/
+
+bool pmPatternRowAmpReadForView (const pmFPAview *view, pmFPAfile *file, const pmConfig *config)
+    {
+        // read the full model in one pass: require the level to be FPA
+        if (view->chip != -1) {
+            psError(PS_ERR_IO, false, "Pattern Row Amplitude must be read at the FPA level");
+            return false;
+        }
+
+        if (!pmPatternRowAmpReadFPA (file)) {
+            psError(PS_ERR_IO, false, "Failed to read Pattern Row Amplitude for fpa");
+            return false;
+        }
+        return true;
+    }
+
+// read in all chip-level Pattern Row Amplitude data for this FPA
+bool pmPatternRowAmpReadFPA (pmFPAfile *file) {
+
+    if (!pmPatternRowAmpReadChips (file)) {
+        psError(PS_ERR_IO, false, "Failed to read Pattern Row Amplitude for chips");
+        return false;
+    }
+
+    return true;
+}
+
+// read the set of tables, one for each chip
+bool pmPatternRowAmpReadChips (pmFPAfile *file) {
+
+    bool haveData, status;
+
+    // loop over the extensions
+    // for each extension, use the extname (eg, XY01.ptn) to assign to a chip
+
+    // move to the start of the file
+    haveData = psFitsMoveExtNum (file->fits, 1, false);
+    if (!haveData) {
+        psError(PS_ERR_IO, false, "Failed to read even the first extension?");
+        return false;
+    }
+
+    while (haveData) {
+
+        // load the header
+        psMetadata *header = psFitsReadHeader(NULL, file->fits); // The FITS header
+        if (!header) psAbort("cannot read model header");
+
+        // load the full model in one shot
+        psArray *model = psFitsReadTable (file->fits);
+        if (!model) psAbort("cannot read model");
+        
+        // determine the chip:
+        char *extname = psMetadataLookupStr (&status, header, "EXTNAME");
+        psLogMsg ("psModules.detrend", 4, "read %ld rows from Pattern Row Amplitude file, extname %s\n", model->n, extname);
+
+        // I expect to find a name of the form: chipName.ptn (eg, XY01.ptn)
+        // where chipName like 'XY01'
+        psAssert (strlen(extname) == 10, "invalid extension %s", extname);
+        psAssert (extname[5] == 'p', "invalid extension %s", extname);
+        psAssert (extname[6] == 't', "invalid extension %s", extname);
+        psAssert (extname[7] == 'n', "invalid extension %s", extname);
+
+        char chipName[5];
+        strncpy (chipName, extname, 4);
+        chipName[4] = 0;
+
+        pmChip *chip = pmConceptsChipFromName (file->fpa, chipName);
+        if (!chip) psAbort ("invalid chip?");
+
+        // parse the model entries
+        for (int i = 0; i < model->n; i++) {
+            psMetadata *row = model->data[i];
+            char *cellName = psMetadataLookupStr(&status, row, "cellname");
+            float *amplitude = psMetadataLookupF32(&status, row, "medianValue");
+
+            int cellNumber = pmChipFindCell (chip, cellName);
+            pmCell *cell = chip->data[cellNumber];
+
+            psMetadataAddF32 (cell->analysis, PS_LIST_TAIL, "PTN.ROW.AMP", PS_META_REPLACE, "", amplitude);
+        }
+        psFree (model);
+        psFree (header);
+
+        // move to the next extension
+        haveData = psFitsMoveExtNum (file->fits, 1, true);
+    }
+
+    return true;
+}
+# endif