pub fn tkd(
local_field: &[f64],
mask: &[u8],
grid: &Grid,
bdir: (f64, f64, f64),
params: &TkdParams,
) -> Vec<f64>Expand description
Truncated k-space division (TKD) for dipole inversion
Computes susceptibility map from local field using direct k-space division with threshold-based truncation.
§Arguments
local_field- Local field values (unwrapped phase / TE / gamma / B0)mask- Binary mask (1 = inside ROI, 0 = outside)grid- Volume grid (dimensions and voxel sizes)bdir- B0 field direction (bx, by, bz)params- TKD parameters (truncation threshold)
§Returns
Susceptibility map (same size as input)