Generalized susceptibility voxel convolution (gSVC) code for susceptibility-induced B0 calculation

Artifact-free method to compute B0 maps from voxelated susceptibility distribution.

This is an efficient method to calculate magnetic field (B0) induced by
susceptibility voxels (chi). Our method computes B0 by circular convolution
between a zero-filled susceptibility matrix and a shifted,
voxel-integrated dipolar field kernel on a grid of size NS + NT - 1 in each
dimension, where NS and NT are the sizes of the susceptibility source and
B0 target grids, respectively. The method is called generalized
susceptibility voxel convolution (gSVC). gSVC can be a fast and reliable
way to compute susceptibility-induced B0 when the susceptibility source
is not co-located with the B0 target volume of interest.