[FieldTrip] Converting units of hdr.sens.tra matrix when computing leadfields

[Knipscheer, B. (Bram)]

Hello,

I am trying to compute leadfields for MEG data with magnetometers and planar gradiometers (neuromag306).
After loading the header I convert the units to meters (hdr.grad  = ft_convert_units(hdr.grad,'m')) to keep units everywhere the same (my headmodel and sourcemodel are also in m).
While this converts variables such as chanpos and coilpos correctly, it also converts the values corresponding to the gradiometers in the hdr.grad.tra matrix.
These values are e.g. -0.5952 and 0.5952 for the gradiometers and are 1 for magnetometers and are converted to -59.52 and 59.92 and 1 respectively.
The next step is using ft_prepare_leadfield and therefore ft_compute_leadfield, which computes the lf for every coil and then transforms this to lf for every sensor with the line lf=sens.tra*lf.
The results is that the lf matrix gets scaled because the units of sens.tra were converted to m. This seems incorrect, as this step should only transform lf per coil to lf per sensor.
Is this indeed incorrect? The question is if the values in hdr.grad.tra for gradiometers have a unit (T/cm or T/m) or that this should be a unitless scaling factor? Additionally, what is the unit of the leadfields?

For clarity, this is the code I'm using:
hdr                                       = ft_read_header(fif-file);
hdr.grad                             = ft_convert_units(hdr.grad,'m');
cfg                                       = [];
cfg.grad                             = hdr.grad;
cfg.headmodel                 = headmodel; %The headmodel is based on OpenMEEG
cfg.sourcemodel              = sourcemodel;
lf                                          = ft_prepare_leadfield(cfg);

Kind regards,

Bram Knipscheer

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

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/20210316/2ea7556f/attachment.htm>