Skip to main content

whqsm

qsm_core::inversion::whqsm

Function whqsm 

Source 
pub fn whqsm(
    local_field: &[f64],
    mask: &[u8],
    grid: &Grid,
    bdir: (f64, f64, f64),
    params: &WhQsmParams,
    progress: impl FnMut(usize, usize),
) -> Vec<f64>
Expand description

Weak-Harmonic QSM dipole inversion.

Jointly estimates the susceptibility map and a residual harmonic background field, returning the susceptibility map (out.x).

§Arguments

  • local_field - Local field values (nx * ny * nz).
  • mask - Binary mask (nx * ny * nz), 1 = inside ROI.
  • grid - Volume grid (dimensions and voxel sizes).
  • bdir - B0 field direction.
  • params - WH-QSM parameters.
  • progress - Progress callback (iteration, max_iter).

§Returns

Susceptibility map.