pub fn nltv(
local_field: &[f64],
mask: &[u8],
magnitude: &[f64],
grid: &Grid,
bdir: (f64, f64, f64),
params: &NltvParams,
progress: impl FnMut(usize, usize),
) -> Vec<f64>Expand description
NLTV dipole inversion (QSM.jl nltv)
§Arguments
local_field- Local field in ppm (nx * ny * nz)mask- Binary mask (nx * ny * nz); the result is zeroed outside itmagnitude- Magnitude image (required), made into the data weightW(see module docs; zero in the padding)grid- Volume grid (dimensions and voxel sizes)bdir- B0 field directionparams- NLTV parametersprogress- Progress callback(iteration, max_iter)
§Returns
Susceptibility map in ppm, masked to the ROI.