pub fn l1qsm(
local_field: &[f64],
mask: &[u8],
grid: &Grid,
bdir: (f64, f64, f64),
params: &L1QsmParams,
progress: impl FnMut(usize, usize),
) -> Vec<f64>Expand description
L1-QSM nonlinear dipole inversion.
§Arguments
local_field- Local field values, ppm-scale (nx * ny * nz).mask- Binary mask (nx * ny * nz), non-zero = inside ROI.grid- Volume grid (dimensions and voxel sizes).bdir- B0 field direction.params- L1-QSM parameters.progress- Progress callback(iteration, max_iter).
§Returns
Estimated susceptibility map (ppm-scale), masked to the ROI.