pub fn l1qsm(
local_field: &[f64],
mask: &[u8],
magnitude: &[f64],
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 in ppm (nx * ny * nz).mask- Binary mask (nx * ny * nz), non-zero = inside ROI.magnitude- Magnitude image, the data weight (see module docs).grid- Volume grid (dimensions and voxel sizes).bdir- B0 field direction.params- L1-QSM parameters.progress- Progress callback(iteration, max_iter).
§Returns
Susceptibility map in ppm, masked to the ROI.