Changeset 18724 for trunk/ppImage/src/ppImageRebinReadout.c

Timestamp:

Jul 24, 2008, 3:58:31 PM (18 years ago)

Author:

eugene

Message:

set masked pixels in binned images to maskedValue if requested (needs modified version of API)

File:

• trunk/ppImage/src/ppImageRebinReadout.c (modified) (4 diffs)

trunk/ppImage/src/ppImageRebinReadout.c

@@ -5 +5 @@
 #include "ppImage.h"
 
-bool ppImageRebinChip (pmConfig *config, pmFPAview *view, char *outName) {
+bool ppImageRebinChip (pmConfig *config, pmFPAview *view, ppImageOptions *options, char *outName) {
 
     pmCell *cell;
@@ -33 +33 @@
 
             // run the rebin code
-            ppImageRebinReadout (outReadout, inReadout, outFile);
+            ppImageRebinReadout (outReadout, inReadout, outFile, options);
         }
 
@@ -61 +61 @@
 
 // XXX this should be made consistent with psImageBinning
-bool ppImageRebinReadout (pmReadout *output, pmReadout *input, pmFPAfile *outFile)
+bool ppImageRebinReadout (pmReadout *output, pmReadout *input, pmFPAfile *outFile, ppImageOptions *options)
 {
     PS_ASSERT_PTR_NON_NULL(output, false);
@@ -80 +80 @@
     int nY = input->image->numRows;
 
+    // we should *either* skip masked pixels or replace masked pixels with the specified value 
+    
+
     // do the rebinning by hand, mean only for test
     for (int yOut = 0; yOut < output->image->numRows; yOut++) {
         for (int xOut = 0; xOut < output->image->numCols; xOut++) {
+            float maskedValue = 0.0;
             float value = 0;
             int nPix = 0;
             for (int yIn = yOut * dY; (yIn < yOut * dY + dY) && (yIn < nY); yIn ++) {
                 for (int xIn = xOut * dX; (xIn < xOut * dX + dX) && (xIn < nX); xIn ++) {
-                    if (input->mask && input->mask->data.U8[yIn][xIn]) continue;
+                    if (input->mask && input->mask->data.U8[yIn][xIn]) {
+                      maskedValue = input->image->data.F32[yIn][xIn];
+                      continue;
+                    }
                     value += input->image->data.F32[yIn][xIn];
                     nPix ++;
                 }
             }
-            output->image->data.F32[yOut][xOut] = value / nPix;
+            if (nPix > 0) {
+              output->image->data.F32[yOut][xOut] = value / nPix;
+            } else {
+              if (options->replaceMasked) {
+                output->image->data.F32[yOut][xOut] = maskedValue;
+              } else {
+                output->image->data.F32[yOut][xOut] = 0.0;
+              }
+            }
         }
     }