Background field removal by solving the Laplacian boundary value problem

The removal of the background magnetic field is a critical step in generating phase images and quantitative susceptibility maps, which have recently been receiving increasing attention. Although it is known that the background field satisfies Laplace's equation, the boundary values of the background field for the region of interest have not been explicitly addressed in the existing methods, and they are not directly available from MRI measurements. In this paper, we assume simple boundary conditions and remove the background field by explicitly solving the boundary value problems of Laplace's or Poisson's equation. The proposed Laplacian boundary value (LBV) method for background field removal retains data near the boundary and is computationally efficient. Tests on a numerical phantom and an experimental phantom showed that LBV was more accurate than two existing methods.