diff --git a/etspy/utils.py b/etspy/utils.py index db4a76b..3b5a8ba 100644 --- a/etspy/utils.py +++ b/etspy/utils.py @@ -231,7 +231,6 @@ def weight_stack( weighted_stack.data * window[np.newaxis, :, :], axis=(1, 2), ) - # Update the status and adjustments for slices with weighted sums below min_mass update_mask = (status != 0) & (weighted_mass < min_mass) status[update_mask] = 0 @@ -249,6 +248,7 @@ def weight_stack( # In this step the severity of the window is calculated again using the value # calculated in the coarse step and the window is made more restrictive in 1 # pixel increments. + status = np.ones(ntilts) status[min_slice] = 0