laura: Local Auto-Regressive Average

This repository contains the source code for the Local AUto-Regressive Average (LAURA) inverse solution as described by de Peralta Menendez et al. (2001, 2004). The code is based on mne-python, a powerful EEG library for python.

Personally, I think this linear inverse solution finds the neural sources underlying M/EEG measurements with great success and is a valuable option among other inverse solutions such as minimum norm estimates and (e)LORETA.

Dependencies

• mne
• scipy
• numpy

That's it!

Installation from PyPi

Use pip to install laura and all its dependencies from PyPi:

pip install laura

Quick Start

The following code demonstrates how to use this package:

from laura import compute_laura

stc = compute_laura(evoked, forward)
stc.plot()

, where evoked is an instance of mne.Evoked and forward is an instance of mne.Forward. For further explanation on mne and its objects please refer to the mne website.

For a more comprehensive tutorial hop over to this notebook!

Feedback

Please leave your feedback and bug reports at lukas_hecker@web.de.

References

Please cite the authors of the LAURA inverse solution appropriately:

[1] Menendez, R. G. D. P., Andino, S. G., Lantz, G., Michel, C. M., & Landis, T. (2001). Noninvasive localization of electromagnetic epileptic activity. I. Method descriptions and simulations. Brain topography, 14(2), 131-137.

[2] de Peralta Menendez, R. G., Murray, M. M., Michel, C. M., Martuzzi, R., & Andino, S. L. G. (2004). Electrical neuroimaging based on biophysical constraints. Neuroimage, 21(2), 527-539.

I would be happy if you would cite this package, too:

LAURA was calculated using the laura python package available at https://github.com/LukeTheHecker/laura.

Limitations

The current implementation is limited to:

• fixed dipole orientations
• time-domain EEG data

Feel free to modify the code and start a pull request!

Troubleshooting

• Having problems with the installation? Check the package requirements.