pub fn hdqsm(
local_field: &[f64],
mask: &[u8],
grid: &Grid,
bdir: (f64, f64, f64),
params: &HdQsmParams,
progress: impl FnMut(usize, usize),
) -> Vec<f64>Expand description
HD-QSM dipole inversion.
§Arguments
local_field- Local field values in ppm (nx * ny * nz).mask- Binary ROI mask (nx * ny * nz), 1 = inside.grid- Volume grid (dimensions and voxel sizes).bdir- B0 field direction.params- HD-QSM parameters.progress- Progress callback(global_iter, max_iter_l1 + max_iter_l2).
§Returns
Susceptibility map in ppm.