Changeset 35149 for branches/czw_branch/20130211/ppImage/src/ppImageDetrendPattern.c

Timestamp:

Feb 11, 2013, 4:50:12 PM (13 years ago)

Author:

watersc1

Message:

Changes that should not conflict with the trunk. This is the camera stage background continuity and the pattern.median example code.

File:

• branches/czw_branch/20130211/ppImage/src/ppImageDetrendPattern.c (modified) (2 diffs)

branches/czw_branch/20130211/ppImage/src/ppImageDetrendPattern.c

@@ -18 +18 @@
     pmCell *cell = NULL;
 
-    assert(options->doPatternRow || options->doPatternCell || options->doPatternContinuity); // do not call if not needed
+    assert(options->doPatternMedian || options->doPatternRow ||
+           options->doPatternCell || options->doPatternContinuity); // do not call if not needed
     assert(inputView->chip != -1);
     assert(inputView->cell == -1);
@@ -35 +36 @@
     // chip.
 
+    if (options->doPatternMedian) {
+      //      bool status;
+      int numCells = chip->cells->n;       // Number of cells
+      psVector *tweak = psVectorAlloc(numCells, PS_TYPE_U8); // Tweak cell?
+      //      pmFPAfile *input = psMetadataLookupPtr(&status, config->files, "PPIMAGE.INPUT");
+      pmFPAview *view = pmFPAviewAlloc(0); // View for local processing
+      
+
+      pmHDU *hdu = pmHDUFromChip(chip);
+      if (psMetadataLookupBool(NULL,hdu->header, "PTRN_MED")) {
+        psLogMsg("ppImage", PS_LOG_INFO, "Not performing median pattern correction as it has already been done.");
+        goto pattern_row;
+      }
+      for (int i = 0; i < chip->cells->n; i++) {
+        view->cell = i;
+        
+        pmCell *cell = chip->cells->data[i]; // Cell of interest
+        
+        if (cell->readouts->n > 1) {
+          psLogMsg("ppImage", PS_LOG_INFO, "Not performing cell continuity correction on video cell.");
+          continue;
+        }
+        
+        bool doPattern = false;
+        if (!doPatternForView(&doPattern, config, chip, view, RECIPE_NAME, "PATTERN.CONTINUITY.SUBSET")) {
+          ESCAPE(false, "Unable to determine whether row pattern matching should be applied.");
+        }
+        if (doPattern) {
+          tweak->data.U8[i] = 0xFF;
+        }
+      }
+      
+      // Tweak the cells
+      if (!pmPatternMedian(chip, tweak, options->patternCellBG, options->patternCellMean,
+                           options->maskValue, options->darkMask,true,true)) {
+        psFree(tweak);
+        psFree(view);
+        return false;
+      }
+      psFree(tweak);
+      psFree(view);
+      
+      psMetadataAddBool(hdu->header, PS_LIST_TAIL, "PTRN_MED",PS_META_REPLACE,"PATTERN.MEDIAN correction applied",true);
+      
+    }
+ pattern_row:
     if (options->doPatternRow) {
         bool status;