pub fn compute_removed_voxels(
mask: &[f64],
r_0: &[f64],
vasc_only: &[f64],
) -> Vec<f64>Expand description
Compute removed voxels mask for offset adjustment
removed_voxels = (mask * R_0) - vasc_only This represents voxels processed in stage 1 but not in stage 2