# Use your own forward leadfield model in an inverse beamformer computation

You can use externally precomputed leadfield matrices for each dipole location on which you want to compute the beamformer solution. That is usually done using the **ft_prepare_leadfield** function, and the main purpose is to speed up the source reconstruction. The output of **ft_prepare_leadfield** can be used as input in **ft_sourceanalysis**, and you should construct a matlab structure that resembles the output of **ft_prepare_leadfield**. The best way of getting to know the format of the precomputed leadfields is by computing the leadfield for a spherical model and look at it.

That can be demonstrated using the following code

% construct some random EEG electrodes on a sphere elec = []; elec.pnt = randn(128,3); dum = sqrt(sum(elec.pnt.^2,2)); elec.pnt = 12 * elec.pnt ./ [dum dum dum]; % scale them to a 12 cm sphere for i=1:128 elec.label{i} = sprintf('%03d', i); end

I suggest that you do not use these electrodes, but instead that you get more meaningfull electrode locations. See for example http://robertoostenveld.nl/?p=5

% create a concentric 3-sphere volume conductor, the radius is the same as for the electrodes vol = []; vol.r = 12 * [0.88 0.92 1.00]; % radii of spheres, the head radius is 12 cm vol.c = [1 1/80 1]; % conductivity vol.o = [0 0 0]; % center of sphere

Using the geometrical description of the sensor locations and the geometrical description of the volume conductor model, you now can compute the leadfield for the locations of interest. By default we scan with the beamformer on a regularly spaced 3D grid of approximately 1 cm resolution. Here is an example for computing the leadfields.

% compute the leadfields that can be used for the beamforming cfg = []; cfg.elec = elec; cfg.vol = vol; cfg.grid.resolution = 2; cfg.grid.unit = 'cm';% same unit as above, i.e. in cm grid = ft_prepare_leadfield(cfg);

The idea now is that, using your own forward model, you construct a “grid” structure with the same elements. Note that the channel numbering should correspond with the channels in the data, and that the EEG leadfields should be averaged referenced.