Skip to content

Commit 916b305

Browse files
committed
fixing noise_mask in Denoise
1 parent eca0251 commit 916b305

File tree

1 file changed

+6
-5
lines changed

1 file changed

+6
-5
lines changed

nipype/interfaces/dipy/preprocess.py

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -230,12 +230,13 @@ def nlmeans_proxy(in_file, settings,
230230

231231
nmask = data[..., 0] > 80
232232
if noise_mask is not None:
233-
if noise_mask.ndim != data.ndim:
234-
nmask = np.array([noise_mask] * data.shape[-1])
233+
noise_mask = np.squeeze(noise_mask)
234+
nmask = np.zeros_like(noise_mask)
235+
nmask[noise_mask > 0] = 1
236+
if nmask.ndim != data.ndim:
237+
nmask = np.array([nmask] * data.shape[-1])
235238

236-
nmask[nmask > 0] = 1
237-
238-
sigma = np.std(data[nmask == 1])
239+
sigma = np.std(data[nmask > 0])
239240
den = nlmeans(data, sigma, **settings)
240241

241242
nb.Nifti1Image(den.astype(hdr.get_data_dtype()), aff,

0 commit comments

Comments
 (0)