Index: trunk/psModules/src/detrend/pmDark.c
===================================================================
--- trunk/psModules/src/detrend/pmDark.c	(revision 17058)
+++ trunk/psModules/src/detrend/pmDark.c	(revision 17228)
@@ -227,4 +227,16 @@
     }
 
+    psImage *counts = pmReadoutAnalysisImage(output->readouts->data[0], PM_READOUT_STACK_ANALYSIS_COUNT,
+                                             xSize, ySize, PS_TYPE_U16, 0);
+    if (!counts) {
+        return false;
+    }
+    psImage *sigma = pmReadoutAnalysisImage(output->readouts->data[0], PM_READOUT_STACK_ANALYSIS_SIGMA,
+                                            xSize, ySize, PS_TYPE_F32, NAN);
+    if (!sigma) {
+        psFree(counts);
+        return false;
+    }
+
     // Iterate over pixels, fitting polynomial
     psVector *pixels = psVectorAlloc(numInputs, PS_TYPE_F32); // Stack of pixels
@@ -269,4 +281,6 @@
                 ro->image->data.F32[yOut][xOut] = poly->coeff->data.F64[k];
             }
+            counts->data.U16[yOut][xOut] = poly->numFit;
+            sigma->data.F32[yOut][xOut] = poly->stdevFit;
         }
     }