fn compute_reliability(
phase: &[f64],
seed: u64,
nx: usize,
ny: usize,
nz: usize,
) -> Vec<f64>Expand description
Per-voxel reliability: the summed squared second difference of the wrapped phase over the 13 directions of the 26-neighbourhood. Lower is better.
Every voxel with a complete neighbourhood — that is, every voxel off the volume
border — gets a real value, whether or not it or its neighbours are in the mask.
Border voxels keep the pseudo-random stand-in, which sorts them last. See the
module docs on extend_mask for why the mask is not consulted here.