def pbHeadModel_Callback(self):
options = QtGui.QFileDialog.DontResolveSymlinks | QtGui.QFileDialog.ShowDirsOnly
HMdir = QtGui.QFileDialog.getExistingDirectory(
self.pbHeadModel,
"Select a Head Model Directory.",
self.pbHeadModel.text(),
options=options)
import os
from braink import read_lfm
import config
import glob
lfm_fname = glob.glob(HMdir + '/*.lfm')
if os.path.exists(HMdir + '/fdmForwardMatrixOriented'):
K = read_lfm.forward(HMdir + '/fdmForwardMatrixOriented',
config.nE)
print("fdmForwardMatrixOriented is loaded. ")
elif len(
lfm_fname) > 0: # elif os.path.exists(HMdir+'/Leadfield.lfm'):
K = read_lfm.lfm(lfm_fname[0], config.nE)
print("Leadfield.lfm is loaded. ")
else:
print("Either fdmForwardMatrixOriented or *.lfm are not loaded. ")
# nV = K.shape[1]
config.K = K
config.nV = K.shape[1]
return K
def pbHeadModel_Callback(self):
options = QtGui.QFileDialog.DontResolveSymlinks | QtGui.QFileDialog.ShowDirsOnly
HMdir = QtGui.QFileDialog.getExistingDirectory(self.pbHeadModel,
"Select a Head Model Directory.",
self.pbHeadModel.text(), options=options)
import os
from braink import read_lfm
import config
import glob
lfm_fname = glob.glob(HMdir+'/*.lfm')
if os.path.exists(HMdir+'/fdmForwardMatrixOriented'):
K = read_lfm.forward(HMdir+'/fdmForwardMatrixOriented', config.nE)
print("fdmForwardMatrixOriented is loaded. ")
elif len(lfm_fname) > 0: # elif os.path.exists(HMdir+'/Leadfield.lfm'):
K = read_lfm.lfm(lfm_fname[0], config.nE)
print("Leadfield.lfm is loaded. ")
else:
print("Either fdmForwardMatrixOriented or *.lfm are not loaded. ")
# nV = K.shape[1]
config.K = K
config.nV = K.shape[1]
return K
def pbHM_Callback(self):
options = QtWidgets.QFileDialog.DontResolveSymlinks | QtWidgets.QFileDialog.ShowDirsOnly
HMdir = QtWidgets.QFileDialog.getExistingDirectory(self.pbHM,
"Select a Head Model Directory.",
self.pbHM.text(), options=options)
import os
from braink import read_lfm
import config
nE = config.nE
nC = nE + 1
if os.path.exists(HMdir+'/fdmForwardMatrixOriented'):
K = read_lfm.forward(HMdir+'/fdmForwardMatrixOriented', config.nE)
print K
print("fdmForwardMatrixOriented is loaded. ")
elif os.path.exists(HMdir+'/Leadfield.lfm'):
K = read_lfm.lfm(HMdir+'/Leadfield.lfm', config.nE)
print("Leadfield.lfm is loaded. ")
print K
else:
print("LFM is not loaded. ")
nV = K.shape[1]
config.K = K
config.nV = nV
config.nC = nC
os.system("open /Applications/EAV/EAV.app")
return
def pbHM_Callback(self):
options = QtWidgets.QFileDialog.DontResolveSymlinks | QtWidgets.QFileDialog.ShowDirsOnly
HMdir = QtWidgets.QFileDialog.getExistingDirectory(
self.pbHM,
"Select a Head Model Directory.",
self.pbHM.text(),
options=options)
import os
from braink import read_lfm
import config
nE = config.nE
nC = nE + 1
if os.path.exists(HMdir + '/fdmForwardMatrixOriented'):
K = read_lfm.forward(HMdir + '/fdmForwardMatrixOriented',
config.nE)
print K
print("fdmForwardMatrixOriented is loaded. ")
elif os.path.exists(HMdir + '/Leadfield.lfm'):
K = read_lfm.lfm(HMdir + '/Leadfield.lfm', config.nE)
print("Leadfield.lfm is loaded. ")
print K
else:
print("LFM is not loaded. ")
nV = K.shape[1]
config.K = K
config.nV = nV
config.nC = nC
os.system("open /Applications/EAV/EAV.app")
return
# -*- coding: utf-8 -*-
"""
Created on Wed Jun 3 13:24:49 2015
@author: jesong1126
"""
import numpy as np
import matplotlib.pyplot as plt
from braink import read_lfm, read_bk
nE = 256
#K = read_lfm.lfm('/Users/jesong1126/Python27/GeoPy/data/108_HM/Leadfield.lfm', nE) # OK!!! left=1 and right=2
K = read_lfm.lfm(
'/Users/jesong1126/Work/EAV_MIE_BrainK/r842_ihm_test1/r842_108_LFM/r842_108_output/r842_108_ori.lfm',
nE)
#K = read_lfm.lfm('/Users/jesong1126/Work/Data/EAV/2941_HM/Leadfield.lfm', nE)
#/Users/jesong1126/Work/EAV_MIE_BrainK/r842_ihm_test1/r842_108_ori.lfm/
nD = K.shape[1] # 2383
import os
os.system("open /Applications/EAV/Mimir.app")
from RabbitMQ import Connection
i = 83
DipI = np.zeros(nD)
DipI[i] = 1
plt.plot(K[:, i])
plt.show()
@author: jesong1126
"""
import numpy as np
import os
os.system("open /Applications/EAV/EAV.app")
from RabbitMQ import Connection
from braink import read_lfm
#from swcm import Imatrix
nC = 257
nE = nC -1
K = read_lfm.lfm('/Users/jesong1126/Python27/data_GeoPy/108_HM/Leadfield.lfm', nE)
nV = K.shape[1]
nV = 2383
Dip1 = np.zeros(nV)
Dip1[0] = 1
c = Connection.Connection('localhost')
c.connect()
#c.sendEEGData(Tmaps[:, i], "Tmap", "")
c.sendOrientedData(Dip1, "Oriented", "4:00pm")
#c.sendTriplesData(bbb, "Triples", "4:00pm")
#c.sendOrientedData(abs(sdenTmaps[:, i]), "sdenTmap", "")
#c.sendTriplesData(sdenTmaps[:,i] , "Triples", "4:00pm")
c.disconnect()
"""
import os
from braink import read_lfm
import glob
HMdir = '/Users/jesong1126/Python27/GeoPy/data/108_HM'
nE = 256
lfm_fname = glob.glob(HMdir+'/*.lfm')
if os.path.exists(HMdir+'/fdmForwardMatrixOriented'):
K = read_lfm.forward(HMdir+'/fdmForwardMatrixOriented', nE)
print("fdmForwardMatrixOriented is loaded. ")
elif len(lfm_fname) > 0: # elif os.path.exists(HMdir+'/Leadfield.lfm'):
K = read_lfm.lfm(lfm_fname[0], nE)
print("Leadfield.lfm is loaded. ")
else:
print("Either fdmForwardMatrixOriented or *.lfm are not loaded. ")
nV = K.shape[1]
os.system("open /Applications/EAV/EAV.app")
#os.system("open /Users/jesong1126/Desktop/EAV.app")
from RabbitMQ import Connection
SampleTime = 0
c = Connection.Connection('localhost')
c.connect()
c.sendEEGData(K[:,SampleTime], "StatEEG", "")
c=Tmaps[:, ii],
s=30,
cmap=plt.get_cmap('seismic'),
alpha=.5)
axes[ii].set_alpha(0.75)
axes[ii].set_title(('MicroState %d' % ii))
axes[ii].set_xticks([])
axes[ii].set_yticks([])
plt.savefig('VGT130MicroStateTMaps.png')
##-----------------------------------------------------------------------------
from braink import read_lfm
from swcm import Imatrix
nE = nC - 1
K = read_lfm.lfm('/Users/jesong1126/Python27/data_GeoPy/108_HM/Leadfield.lfm',
nE)
nV = K.shape[1]
alpha = 0.001
Imat = Imatrix.MN(alpha, K)
sdenTmaps = np.zeros((nV, Tmaps.shape[1]))
for i in range(Tmaps.shape[1]):
sdenTmaps[:, i] = np.dot(Imat, Tmaps[:, i])
import os
os.system("open /Applications/EAV/EAV.app")
from RabbitMQ import Connection
# -*- coding: utf-8 -*-
"""
Created on Wed Jun 3 13:24:49 2015
@author: jesong1126
"""
import numpy as np
import matplotlib.pyplot as plt
from braink import read_lfm, read_bk
nE = 256
#K = read_lfm.lfm('/Users/jesong1126/Python27/GeoPy/data/108_HM/Leadfield.lfm', nE) # OK!!! left=1 and right=2
K = read_lfm.lfm('/Users/jesong1126/Work/EAV_MIE_BrainK/r842_ihm_test1/r842_108_LFM/r842_108_output/r842_108_ori.lfm', nE)
#K = read_lfm.lfm('/Users/jesong1126/Work/Data/EAV/2941_HM/Leadfield.lfm', nE)
#/Users/jesong1126/Work/EAV_MIE_BrainK/r842_ihm_test1/r842_108_ori.lfm/
nD = K.shape[1] # 2383
import os
os.system("open /Applications/EAV/Mimir.app")
from RabbitMQ import Connection
i=83; DipI= np.zeros(nD); DipI[i] = 1
plt.plot(K[:,i]); plt.show()
c = Connection.Connection('localhost')
c.connect()
c.sendOrientedData(DipI, "Ori", "4:00pm")
c.sendEEGData(K[:, i], "EEG", "")
