Skip to main content

dipole_kernel

qsm_core::kernels::dipole

Function dipole_kernel 

Source 
pub fn dipole_kernel(
    nx: usize,
    ny: usize,
    nz: usize,
    vsx: f64,
    vsy: f64,
    vsz: f64,
    bdir: (f64, f64, f64),
) -> Vec<f64>
Expand description

Generate dipole kernel in k-space

Creates the dipole kernel D(k) = 1/3 - (kz)²/|k|² (for B = [0,0,1]) centered at index (0, 0, 0) (not shifted).

§Arguments

  • nx, ny, nz - Array dimensions
  • vsx, vsy, vsz - Voxel sizes in mm
  • bdir - B0 field direction as (bx, by, bz), default (0, 0, 1)

§Returns

Flattened dipole kernel array of size nxnynz in C order