Skip to content

Commit

Permalink
Fix inconsistent indentation
Browse files Browse the repository at this point in the history
  • Loading branch information
mcraig-ibme committed Jun 17, 2021
1 parent 996be53 commit f08aabc
Showing 1 changed file with 15 additions and 16 deletions.
31 changes: 15 additions & 16 deletions quantiphyse/packages/core/smoothing/process.py
Original file line number Diff line number Diff line change
Expand Up @@ -50,23 +50,22 @@ def run(self, options):
if data.nvols > 1:
sigmas += [0, ]

#output = scipy.ndimage.filters.gaussian_filter(data.raw(), sigmas, order=order, mode=mode)
output = self._norm_conv(data.raw(), sigmas, order=order, mode=mode)
self.ivm.add(NumpyData(output, grid=data.grid, name=output_name), make_current=True)

def _norm_conv(self, data, sigma, **kwargs):
"""
Normalized convolution
This is a way to compensate for data having nan/infinite values.
Taken from stackoverflow.com/questions/18697532/gaussian-filtering-a-image-with-nan-in-python
"""
v = data.copy()
v[~np.isfinite(data)] = 0
vv = scipy.ndimage.filters.gaussian_filter(v, sigma, **kwargs)

w = 0*data.copy()+1
w[~np.isfinite(data)] = 0
ww = scipy.ndimage.filters.gaussian_filter(w, sigma, **kwargs)
return vv/ww
"""
Normalized convolution
This is a way to compensate for data having nan/infinite values.
Taken from stackoverflow.com/questions/18697532/gaussian-filtering-a-image-with-nan-in-python
"""
v = data.copy()
v[~np.isfinite(data)] = 0
vv = scipy.ndimage.filters.gaussian_filter(v, sigma, **kwargs)

w = 0*data.copy()+1
w[~np.isfinite(data)] = 0
ww = scipy.ndimage.filters.gaussian_filter(w, sigma, **kwargs)

return vv/ww

0 comments on commit f08aabc

Please sign in to comment.