Function iharperella 

pub fn iharperella(
    phase: &[f64],
    mask: &[u8],
    grid: &Grid,
    params: &HarperellaParams,
    progress: impl FnMut(usize, usize),
) -> (Vec<f64>, Vec<u8>)

iHARPERELLA — improved integrated phase unwrapping and background removal

More robust low-frequency background suppression than HARPERELLA. Estimates exterior Laplacian by directly minimizing the weighted resulting phase rather than enforcing SMV uniformity.

Arguments

• phase - Wrapped phase in radians (nx * ny * nz)
• mask - Binary brain mask (nx * ny * nz), 1 = inside
• grid - Volume dimensions and voxel sizes
• params - Algorithm parameters (radius is unused in iHARPERELLA)
• progress - Progress callback (iteration, max_iter)

Returns

(tissue_phase, mask)