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) |