Skip to main content

mcpc3ds_single_coil

qsm_core::utils::multi_echo

Function mcpc3ds_single_coil 

Source 
pub fn mcpc3ds_single_coil(
    phases: &[impl AsRef<[f64]>],
    mags: &[impl AsRef<[f64]>],
    tes: &[f64],
    mask: &[u8],
    sigma: [f64; 3],
    echoes: [usize; 2],
    nx: usize,
    ny: usize,
    nz: usize,
) -> (Vec<Vec<f64>>, Vec<f64>)
Expand description

MCPC-3D-S phase offset estimation for single-coil multi-echo data

Implements the MCPC-3D-S algorithm from MriResearchTools.jl for single-coil data. This estimates and removes the phase offset (φ₀) from each echo.

§Arguments

  • phases - Phase data for all echoes, shape [n_echoes][nxnynz]
  • mags - Magnitude data for all echoes, shape [n_echoes][nxnynz]
  • tes - Echo times in ms
  • mask - Binary mask
  • sigma - Smoothing sigma in voxels [sx, sy, sz], default [10, 10, 5]
  • echoes - Which echoes to use for HIP calculation, default [0, 1] (first two)
  • nx, ny, nz - Dimensions

§Returns

(corrected_phases, phase_offset) where:

  • corrected_phases: phases with offset removed
  • phase_offset: estimated phase offset