Skip to main content

rts

qsm_core::inversion::rts

Function rts 

Source 
pub fn rts(
    local_field: &[f64],
    mask: &[u8],
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f64,
    vsy: f64,
    vsz: f64,
    bdir: (f64, f64, f64),
    delta: f64,
    mu: f64,
    rho: f64,
    tol: f64,
    max_iter: usize,
    lsmr_iter: usize,
) -> Vec<f64>
Expand description

RTS dipole inversion (optimized)

§Arguments

  • local_field - Local field values (nx * ny * nz)
  • mask - Binary mask (nx * ny * nz), 1 = inside ROI
  • nx, ny, nz - Array dimensions
  • vsx, vsy, vsz - Voxel sizes in mm
  • bdir - B0 field direction
  • delta - Threshold for ill-conditioned region (typically 0.15)
  • mu - Regularization parameter for well-conditioned region (typically 1e5)
  • rho - ADMM penalty parameter (typically 10)
  • tol - Convergence tolerance
  • max_iter - Maximum ADMM iterations
  • lsmr_iter - LSMR iterations for step 1 (typically 4)

§Returns

Susceptibility map