Skip to main content

Table 1 OpenMEEG demo script in Python

From: OpenMEEG: opensource software for quasistatic bioelectromagnetics

Python code

Import openmeeg as om

# Load data

condFile = 'om_demo.cond'

geomFile = 'om_demo.geom'

dipoleFile = 'cortex.dip'

squidsFile = 'meg_squids.txt'

electrodesFile = 'eeg_electrodes.txt'

geom = om.Geometry()

geom.read(geomFile,condFile)

dipoles = om.Matrix()

dipoles.load(dipoleFile)

squids = om.Sensors()

squids.load(squidsFile)

electrodes = om.Matrix()

electrodes.load(electrodesFile)

# Compute forward problem

gaussOrder = 3; # Integration order

hm = om.HeadMat(geom,gaussOrder)

hminv = hm.inverse()

dsm = om.DipSourceMat (geom, dipoles, gaussOrder)

ds2 mm = om.DipSource2MEGMat (dipoles, squids).

h2 mm = om.Head2MEGMat (geom, squids)

h2em = om.Head2EEGMat (geom, electrodes).

gain_meg = om.GainMEG (hminv, dsm, h2 mm, ds2mm)

gain_eeg = om.GainEEG (hminv, dsm, h2em)

  1. Demo script for computing MEG and EEG forward problems with OpenMEEG in Python.