Changeset 17730

Timestamp:

May 16, 2008, 4:38:46 PM (18 years ago)

Author:

Paul Price

Message:

Adding function to set image,weight pixels to NAN where they are masked.

Location:

trunk/psModules/src/camera

Files:

• pmFPAMaskWeight.c (modified) (1 diff)
• pmFPAMaskWeight.h (modified) (2 diffs)

trunk/psModules/src/camera/pmFPAMaskWeight.c

@@ -404 +404 @@
     return true;
 }
+
+
+
+bool pmSubtractionMaskApply(pmReadout *readout, psMaskType maskVal)
+{
+    PM_ASSERT_READOUT_NON_NULL(readout, false);
+    PM_ASSERT_READOUT_IMAGE(readout, false);
+    PM_ASSERT_READOUT_MASK(readout, false);
+
+    int numCols = image->numCols, numRows = image->numRows; // Size of image
+    psMaskType **maskData = mask->data.PS_TYPE_MASK_DATA; // Dereference mask
+    psF32 **imageData = image->data.F32;// Dereference image
+    psF32 **weightData = weight->data.F32; // Dereference weight map
+
+    for (int y = 0; y < numRows; y++) {
+        for (int x = 0; x < numCols; x++) {
+            if (maskData[y][x] & maskVal) {
+                imageData[y][x] = NAN;
+                if (weight) {
+                    weightData[y][x] = NAN;
+                }
+            }
+        }
+    }
+
+    return true;
+}

trunk/psModules/src/camera/pmFPAMaskWeight.h

@@ -5 +5 @@
  * @author Eugene Magnier, IfA
  *
- * @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
- * @date $Date: 2008-05-15 03:25:10 $
+ * @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
+ * @date $Date: 2008-05-17 02:38:46 $
  * Copyright 2005-2006 Institute for Astronomy, University of Hawaii
  */
@@ -115 +115 @@
     );
 
+/// Apply a mask to the image and weight map
+///
+/// Unfortunately, image subtraction may result in a bi-modal image in masked areas, which can upset image
+/// statistics (very important for quantising images so that a product can be written out!).  This function
+/// sets masked areas to NAN in the image and weight.
+bool pmSubtractionMaskApply(pmReadout *readout, ///< Readout to mask
+                            psMaskType maskVal ///< Mask value for which to apply mask
+    );
+
+
 /// @}
 #endif