pub fn nltv(
local_field: &[f64],
mask: &[u8],
grid: &Grid,
bdir: (f64, f64, f64),
params: &NltvParams,
progress: impl FnMut(usize, usize),
) -> Vec<f64>Expand description
NLTV dipole inversion using iteratively reweighted ADMM
§Arguments
local_field- Local field values (nx * ny * nz)mask- Binary mask (nx * ny * nz), 1 = inside ROIgrid- Volume grid (dimensions and voxel sizes)bdir- B0 field directionparams- NLTV parametersprogress- Progress callback(iteration, total_iterations)
§Returns
Susceptibility map