Fast realistic modeling in bioelectromagnetism using lead-field interpolation

The practical use of realistic models in bioelectromagnetism is limited by the time-consuming amount of numerical calculations. We propose a method leading to much higher speed than currently available, and compatible with any kind of numerical methods (boundary elements (BEM), finite elements, finite differences). Illustrated with the BEM for EEG and MEG, it applies to ECG and MCG as well. The principle is two-fold. First, a Lead-Field matrix is calculated (once for all) for a grid of dipoles covering the brain volume. Second, any forward solution is interpolated from the pre-calculated Lead-Fields corresponding to grid dipoles near the source. Extrapolation is used for shallow sources falling outside the grid. Three interpolation techniques were tested: trilinear, second-order Bézier (Bernstein polynomials), and 3D spline. The trilinear interpolation yielded the highest speed gain, with factors better than x10,000 for a 9,000-triangle BEM model. More accurate results could be obtained with the Bézier interpolation (speed gain approximately 1,000), which, combined with a 8-mm step grid, lead to intrinsic localization and orientation errors of only 0.2 mm and 0.2 degrees. Further improvements in MEG could be obtained by interpolating only the contribution of secondary currents. Cropping grids by removing shallow points lead to a much better estimation of the dipole orientation in EEG than when solving the forward problem classically, providing an efficient alternative to locally refined models. This method would show special usefulness when combining realistic models with stochastic inverse procedures (simulated annealing, genetic algorithms) requiring many forward calculations.