Meg Source Space Projection

---

Sourcespace projection code

show code

collapse code

""" Source Space Projection
from meg import sourcespaceprojection
 
qn=array([.01, .2, -1])
grid=array([1, 1, 3])
L=leadfield.calclf(grid)
slices=pdf.GetSliceRangeMEGChannels(0,678)
 
SSP=sourcespaceprojection.calc(data, L, qn)
"""
 
from numpy import invert,linalg, shape,array, dot, zeros
from msiread import getposted
 
def calc(data, L, qn):
    Li=linalg.pinv(L); #3x248 inverted leadfield
    #s=(qn.conj().transpose())*(Li) #s=(qn')*(Li); 
    s=dot(qn,Li); #should be NumChan X 3
    SSPvirch=zeros((1,shape(data)[1]),float) #make an empty virchan array. Should by 1 X number of timepoints.
    datadim = shape(data) #data dimensions. Assuming 2D for now. Maybe later can accomadate 3D.
    nump = datadim[1] #number of points
 
    for eachpoint in range(0,nump):
        SSPvirch[0,eachpoint] = dot(data[:,eachpoint],s)
 
    return SSPvirch

show howto

collapse howto

from meg import sourcespaceprojection

qn=array([.01, .2, -1])
grid=array([1, 1, 3])
L=leadfield.calclf(grid)
slices=pdf.GetSliceRangeMEGChannels(0,678)

SSP=sourcespaceprojection.calc(data, L, qn)

---