Merge pull request #2883 from yarikoptic/enh-stddev-once

ENH: Threshold stddev once only in TSNR

Author: effigies

nipype/algorithms/confounds.py

@@ -799,8 +799,9 @@ def _run_interface(self, runtime):
         meanimg = np.mean(data, axis=3)
         stddevimg = np.std(data, axis=3)
         tsnr = np.zeros_like(meanimg)
-        tsnr[stddevimg > 1.e-3] = meanimg[stddevimg > 1.e-3] / stddevimg[
-            stddevimg > 1.e-3]
+        stddevimg_nonzero = stddevimg > 1.e-3
+        tsnr[stddevimg_nonzero] = meanimg[stddevimg_nonzero] / stddevimg[
+            stddevimg_nonzero]
         img = nb.Nifti1Image(tsnr, img.affine, header)
         nb.save(img, op.abspath(self.inputs.tsnr_file))
         img = nb.Nifti1Image(meanimg, img.affine, header)