[FieldTrip] Normalization of beamformer leadfields

[Klink-3, N.E.C. van]

Dear all,

I want to do SAM beamformer source localization on single trial EEG data. I would like to normalize the leadfields to correct for depth, like mentioned in the lmcv beamformer tutorial: (http://www.fieldtriptoolbox.org/tutorial/beamformer_lcmv)

cfg                  = [];
cfg.elec             = hdr.elec;  % electrode distances
cfg.headmodel        = vol;   % volume conduction headmodel
cfg.grid             = grid;  % normalized grid positions
cfg.channel          = {'EEG'};
cfg.normalize        = 'yes'; % to remove depth bias (Q in eq. 27 of van Veen et al, 1997).
lf                   = ft_prepare_leadfield(cfg);

However when I look what happens with cfg.normalize='yes', the following script is used in ft_compute_leadfield, from line 570:

  case 'yes'
    for ii=1:Ndipoles
      tmplf = lf(:, (3*ii-2):(3*ii));
      if normalizeparam==0.5
        % normalize the leadfield by the Frobenius norm of the matrix
        % this is the same as below in case normalizeparam is 0.5
        nrm = norm(tmplf, 'fro');
      else
        % normalize the leadfield by sum of squares of the elements of the leadfield matrix to the power "normalizeparam"
        % this is the same as the Frobenius norm if normalizeparam is 0.5
        nrm = sum(tmplf(:).^2)^normalizeparam;
      end
      if nrm>0
        tmplf = tmplf ./ nrm;
      end
      lf(:, (3*ii-2):(3*ii)) = tmplf;
    end

This seems to me as independent of the dipole location, and does not use an estimate of the noise spectrum as in Eq 27 of van Veen et al 1997.
DICS beamformer has the option to estimate the noise spectrum with 'projectnoise', but SAM beamformer does not have that option. SAM does something with noise and a lambda, which is noise regularization I guess (beamformer_sam from line 102). I use Fieldtrip 20170212.

My main question: how do I correct the leadfields for depth bias?

Thanks in advance,
Nicole

------------------------------------------------------------------------------

De informatie opgenomen in dit bericht kan vertrouwelijk zijn en is
uitsluitend bestemd voor de geadresseerde. Indien u dit bericht onterecht
ontvangt, wordt u verzocht de inhoud niet te gebruiken en de afzender direct
te informeren door het bericht te retourneren. Het Universitair Medisch
Centrum Utrecht is een publiekrechtelijke rechtspersoon in de zin van de W.H.W.
(Wet Hoger Onderwijs en Wetenschappelijk Onderzoek) en staat geregistreerd bij
de Kamer van Koophandel voor Midden-Nederland onder nr. 30244197.

Denk s.v.p aan het milieu voor u deze e-mail afdrukt.

------------------------------------------------------------------------------

This message may contain confidential information and is intended exclusively
for the addressee. If you receive this message unintentionally, please do not
use the contents but notify the sender immediately by return e-mail. University
Medical Center Utrecht is a legal person by public law and is registered at
the Chamber of Commerce for Midden-Nederland under no. 30244197.

Please consider the environment before printing this e-mail.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20170301/d0730255/attachment-0001.html>