pub fn ndi(
local_field: &[f64],
mask: &[u8],
grid: &Grid,
bdir: (f64, f64, f64),
params: &NdiParams,
progress: impl FnMut(usize, usize),
) -> Vec<f64>Expand description
Nonlinear Dipole Inversion (NDI)
§Arguments
local_field- Local field values, ppm-scale (nx * ny * nz)mask- Binary mask (nx * ny * nz), non-zero = inside ROIgrid- Volume grid (dimensions and voxel sizes)bdir- B0 field directionparams- NDI parametersprogress- Progress callback(iteration, max_iter)
§Returns
Estimated susceptibility map (ppm-scale), masked to the ROI.