Module multi_echo 

Multi-echo phase combination utilities

Implements MCPC-3D-S (Multi-Channel Phase Combination - 3D - Smoothed) algorithm and weighted B0 calculation.

Reference: Eckstein, K., Dymerska, B., Bachrata, B., Bogner, W., Poljanc, K., Trattnig, S., Robinson, S.D. (2018). “Computationally Efficient Combination of Multi-channel Phase Data From Multi-echo Acquisitions (ASPIRE).” Magnetic Resonance in Medicine, 79:2996-3006. https://doi.org/10.1002/mrm.26963

Reference implementation: https://github.com/korbinian90/MriResearchTools.jl

Structs

LinearFitParams

Parameters for multi-echo linear fit.

LinearFitResult

Result of multi-echo linear fit

Mcpc3dsParams

Parameters for MCPC-3D-S phase combination.

Enums

B0WeightType

B0 weighting types matching MriResearchTools.jl

Constants

TWO_PI 🔒

Functions

calculate_b0_weighted

Calculate B0 field from unwrapped phase using weighted averaging

compute_reliability_mask 🔒

Compute reliability mask by thresholding fit residuals

field_to_hz

Convert field from rad/s to Hz

find_seed_point 🔒

Find a good seed point (center of mass of the mask)

gaussian_smooth_3d_phase

3D Gaussian smoothing for phase data (handles phase wrapping)

gaussian_smooth_3d_separable 🔒

Separable 3D Gaussian smoothing

hermitian_inner_product

Compute Hermitian Inner Product (HIP) between two echoes

idx3d 🔒

Index into 3D array (Fortran/column-major order)

make_gaussian_kernel 🔒

Create 1D Gaussian kernel

mcpc3ds_b0_pipeline

Full MCPC-3D-S + B0 calculation pipeline

mcpc3ds_single_coil

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

multi_echo_linear_fit

Multi-echo linear fit with magnitude weighting

unwrap_with_romeo 🔒

Unwrap phase using ROMEO algorithm

wrap_to_pi 🔒

Wrap angle to [-π, π]