Changeset 10240

Timestamp:

Nov 28, 2006, 12:33:43 PM (20 years ago)

Author:

Paul Price

Message:

For non-positive gain, use read noise only.

File:

• trunk/psModules/src/camera/pmFPAMaskWeight.c (modified) (1 diff)

trunk/psModules/src/camera/pmFPAMaskWeight.c

@@ -207 +207 @@
     }
 
-    // Set weight image to the variance in ADU = f/g + rn^2
-    psImage *image = readout->image;    // The image pixels
-    readout->weight = (psImage*)psBinaryOp(readout->weight, image, "/", psScalarAlloc(gain, PS_TYPE_F32));
-
-    // a negative weight is non-sensical. if the image value drops below 1, the weight must be 1.
-    readout->weight = (psImage*)psUnaryOp(readout->weight, readout->weight, "abs");
-    readout->weight = (psImage*)psBinaryOp(readout->weight, readout->weight, "max", psScalarAlloc(1, PS_TYPE_F32));
+    if (gain > 0.0) {
+        // Set weight image to the variance in ADU = f/g + rn^2
+        psImage *image = readout->image;    // The image pixels
+        readout->weight = (psImage*)psBinaryOp(readout->weight, image, "/", psScalarAlloc(gain, PS_TYPE_F32));
+
+        // a negative weight is non-sensical. if the image value drops below 1, the weight must be 1.
+        readout->weight = (psImage*)psUnaryOp(readout->weight, readout->weight, "abs");
+        readout->weight = (psImage*)psBinaryOp(readout->weight, readout->weight, "max",
+                                               psScalarAlloc(1, PS_TYPE_F32));
+    } else {
+        // Gain is negative or zero --- just use the read noise
+        psImageInit(readout->weight, 0.0);
+    }
 
     readout->weight = (psImage*)psBinaryOp(readout->weight, readout->weight, "+",