Recent advances in the analysis of biomagnetic signals
|
|
- Vanessa Logan
- 5 years ago
- Views:
Transcription
1 Recent advances in the analysis of biomagnetic signals Kensuke Sekihara Mind Articulation Project, Japan Science and echnology Corporation okyo Metropolitan Institute of echnology
2 Application of spatial filter techniques to biomagnetic source localization
3 Definitions b () t 1 b () t 2 data vector: b(t) = b () t M b () t 1 b () t 2 b () t M b (): t the jth sensor recording at t j data covariance matrix: D = b() tb () t represents time average
4 Source moment magnitude at r = [ xyz,, ] orientation: and at t: s( r,) t h( r,t ) = [ η x ( r,t ), η y ( r,t ), η z ( r,t)] source moment vector: ηx( r,t) s x( r,t) sr (,) t = s ( r,t ) ηy( r,t) = s y( r,t) ηz( r,t) s z( r,t) z θ y ϕ x η = sinθ cosϕ x η = sinθ sinϕ η y x = cosθ
5 Lead field vector for the j th sensor l ( r) = [ l x ( r), l y ( r), l z ( r)] j j j j Lead field matrix for the whole sensor array x y z l ( r) ( ) ( ) ( ) 1 l r l r l r x y z l ( r) ( ) ( ) ( ) 2 l r l r l r Lr ( ) = = x y z l ( r) l ( ) l ( ) l ( ) M r r r M M M
6 Basic relationship b () t = ( ) ( ) j l rsr,t dr j or b() t = Lrsr ( ) (,t) dr Problem of source localization: Estimate sr (,) t from the measurement b() t Estimate
7 Spatial filter bt () 1 sr ˆ(,) t = w ()() rb t = [ w (), r, w ()] r 1 M b () t M estimate of sr (,) t weight vector (neglecting the explicit time notation) b = Lrsr ( ) ( ) dr } ˆ sr ( ) = ( ) ( ') ( ') d ' w rlr sr r ˆs( r)= w ( rb ) R( r, r' ) Resolution kernel
8 Resolution kernel: sr ˆ( ) = R(, rrsr ') ( ') dr' he weight w(r) must be chosen so that the resolution kernel has a reasonable shape. What is reasonable shape? peak at r small mail-lobe width low side-lobe amplitude r main-lobe width side-lobe amplitude
9 Non-adaptive weight wr () is data independent Adaptive weight wr () is data dependent
10 Data-independent (non-adaptive) weight Spatial resolution is limited by sensor-array configuration. Final results are not influenced by source temporal correlation. Data-dependent (adaptive) weight Spatial resolution can exceed the limit imposed by the sensor-array configuration. Strong temporal correlation among source activities severely degrade the quality of final estimation results
11 Data-independent (non-adaptive) weight minimum-norm estimate (Hamalainen et al.) he weight wr ( ) is obtained by [ R ] 2 min (, rr') δ ( r r ) dr 1 w () r = L ( rg ), where G = lrl ( ) ( r) dr ij, i j 1 Inverse solution: sr ˆ( ) = L ( rgb )
12 Property of G matrix G = l ( rl ) ( r) dr i,j i j Overlaps of sensor lead fields is large Overlaps of sensor lead fields is small G is poorly conditioned G has a small condition number
13 G is usually calculated by introducing pixel grid r j b = Lrsr ( ) ( ) dr = Lrsr ( ) ( ) j = 1 sr ( ) 1 = Lr ( ),, Lr ( ) = Ls 1 L ( ) N sr N sn herefore G = LL and w () r = L ()( r LL) or N N N N N j j N N N 1 w () r = L ()( r LL + I) N N γ 1 [regularized versi on]
14 Minimum-norm weight with normalized lead field wr LL I Lr Lr 1 () = ( + γ ) ()/ () N N where L N Lr ( ) Lr ( ) 1 N =,, Lr ( ) Lr ( ) 1 N Minimum-norm estimate with normalized weight Dale et al., Valdes et al. Calculate 2 2 q ˆ( r) = w () rb / wr ()
15 no normalization z (cm) normalized lead field used normalized weight used
16 Source imaging experiments Auditory-evoked field were measured using 148-channel whole-head sensor array (Magnes 2500). Stimulus: 1-kHz pure tone applied to subject s left ear Data acquisition: 1 khz sampling frequency, Hz bandpass filtering, 100 epochs averaged
17
18
19
20 Linear-estimation-based methods LOREA: impose the maximum-smoothness constraint. (Pascual-Marqui et al.) fmri constraint: constrain solution based on fmri results. (Dale et al.) FOCUSS: obtain a focal solution iteratively. (Gorodnitsky et al.) Bayesian approach: impose prior assumptions. (Schmidt et al.) l 1 -norm approach: use the l 1 -norm, instead of using the l 2 -norm. (Matuura et al., Uutela et al., Beucker et al.)
21 Data dependent (adaptive) weight minimum-variance beamformer min w wdw subject to wlrh (, ) = l (, r h) D l (, r h) Db w (, r h) = and sr ˆ(, h) = 1 1 l (, r h) Dlr (, h) l (, r h) Dlr (, h) Lrh () beamformer pointing orientation
22 Generalized Wiener estimate: ˆ ( ) 1 s = RL LRL + C b N N N N source covariance matrix noise covariance matrix Use R I, and C γi (assume no prior-information) ˆ ( ) 1 s = L LL + γi b N N N N minimum-norm estimate
23 Generalized Wiener estimate: ˆ ( ) 1 s = RL LRL + C b N N N N iuse DLRL = + C, and N N iassume off-diagonals of R are all zero Relashonship, sˆ p 1 = Rl pp Db P 2 1 sˆ = R leads to R = 1/( l D l ) p pp pp P p s 1 1 ˆ ( l p )/( l = Db D l p) P P minimum-variance estimate
24 Problems when applying MV beamformer to MEG source localization (1) How to determine beamformer orientation h. (2) Lr () norm artifact. (3) Output SNR degradation due to the matrix inversion.
25 How to determine beamformer orientation h? he weight wr (, h) is calculated for r and h. search all directions (Robinson and Vrba). use the MUSIC algorithm to determine h (Sekihara et al.). use vector beamformer formulation. (Spencer et al.van Veen et al.) What happens if the beamformer orientation is different from the source orientation? Severe signal-intensity loss
26 Computer simulation for calculating beamformer angular response 37-channel sensor array A single source exists 6-cm below the sensor array 6cm δθ assumed source waveform beamformer orientation different from the source orientation with δθ normalized value δθ degree Calculated angular response
27 A single source exist at r with orientation h and time course st () measured data: b() t = [ ηl () r + ηl () r + ηlr ()] st () xx yy zz When beamformer weight has a wrong orientation h ( h) w (, r h )() b t = w (, r h )[ ηxx l ()() r st + ηyy l ()() r st + ηzlr z()() st] virtual three coherent sources 0
28 Vector beamformer formulation hree weight vectors [ wx(), r wy(), r wr z()] detects the x, y, and z component of the source moment, separately. min wdw subject to wlrf (, ) = 1, wlrf (, ) = 0, wlrf (, ) = 0 x x x x x y x z min wdw subject to wlrf (, ) = 0, wlrf (, ) = 1, wlrf (, ) = 0 y y y x y y y z min wdw subject to wlrf (, ) = 0, wlrf (, ) = 0, wlrf (, ) = 1 z z z x z y z z [ w (), r w (), r w ()] r = DLr ()[ L () rdlr ()] x y z and ˆ( t) = [ () ()] () () t s L rdlr L rdb f = [1,0,0], f = [0,1,0], f = [0,0,1] x y z
29 Calculated beamformer angular response cosine curve normalized value normalized value degree scalar beamformer degree vector beamformer
30 Lr () norm artifact Severe artifacts arise near the center of the sphere, when the spherical conductor model is used. o avoid these artifacts use normalized lead field (Van Veen et al.) use normalized weight (Robinson et al., Sekihara et al.)
31 assumed source waveform 37-channel sensor array latency (ms) z (cm) } three sources -15 reconstruction region relative magnitude latency (ms) y (cm) sphere for the forward modeling generated magnetic field
32 ime-averaged reconstruction sr ˆ(,) t 2 z (cm) z (cm) z (cm) y (cm) no normalization y (cm) normalized lead field used
33 Borgiotti-Kaplan beamformer min w wdw subject to ww = 1 w (, r h) = l (, r h) D and 1 2 (, h) (, h) l r Dlr sr ˆ(, h) 2 = 1 (, h) Dl(, r h) 2 (, h) Dlr (, h) l r l r
34 Vector extension of Borgiotti-Kaplan beamformer min wdw subject to ww = 1, wlrf (, ) = 0, wlrf (, ) = 0 x x x x x y x z min wdw subject to wlrf (, ) = 0, ww = 1, wlrf (, ) = 0 y y y x y y y z min wdw subject to wlrf (, ) = 0, wlrf (, ) = 0, ww = 1 z z z x z y z z let µ = x, y or z w () r µ = DLr ()[ L () rdlr ()] fµ f Wf µ µ Ω = [ L () rdlr ()] L () rd Lr ()[ L () rdlr ()]
35 ime-averaged reconstruction sr ˆ(,) t 2 z (cm) z (cm) y (cm) normalized lead field used y (cm) normalized weight used (B-K beamformer results)
36 Resolution kernel for BK beamformer
37 Output SNR degradation for spatio-temporal reconstruction Signal-to-noise ratio of the beamformer output is severely degraded even by a small error in the estimated lead field his is caused by the use of direct matrix inversion o avoid this, use regularized inverse (Robinson et al.) use eigenspace projection (Sekihara et al.)
38 Eigendecomposition of D Eigenspace projection λ L 0 S D = U λ P U = U U 0 LN λ M U = [ e,, e e,, e ] = [ E E ] 1 P P+ 1 M S N E E S N Extension to eigenspace projection beamformer w = EEw, where µ= xy, or z S S µ µ
39 SNR consideration Output of eigenspace-projected BK beamformer [ l () r G lr ()] 2 S 1-1 where G = EL E, G = E L E 2 S S S S N N N N G lr S [ l () r ()] Output of BK beamformer [ l () r G lr ()] 2 S 2 2 [ l () r G lr () + ] S e G e N overall error in estimating lr ( ) Even when e is small, e 2 G N e may not be small
40 assumed source waveform 37-channel sensor array latency (ms) z (cm) } three sources -15 reconstruction region relative magnitude latency (ms) y (cm) sphere for the forward modeling generated magnetic field
41 Spatio-tempotal reconstruction by vector-extended BK beamformer
42 Spatio-tempotal reconstruction by vector-extended BK beamformer 1 with regularized inverse, ( D + γi)
43 Spatio-tempotal reconstruction by vector-extended BK beamformer with eigen-space projection
44 Eigen-space projection does not preserve the null constraints hat is, EEwx ly() r 0, EEwx lz() r 0, S S S S EEwy lx() r 0, EEwy lz() r 0, S S S S EEwz lx() r 0, EEwz ly() r 0. S S S S his fact does not cause a problem. (Poster: 167b)
45
46
47
48 7 mm Poster: 92a
49
50 right auditory cortex activation left auditory cortex activation correlation coefficient: 0.97
51 Summary Adaptive spatial filter techniques can provide a spatial resolution higher than that of non adaptive techniques. his is because the spatial resolution for non-adaptive techniques is limited by the sensor configuration but adaptive techniques can exceed this limit. Correlated source activities, however, affect the quality of the results obtained by adaptive techniques. herefore Adaptive techniques may be suited to observe relatively small cortical regions with high spatial resolution, and non-adaptive techniques may be suited to observe whole-brain activities.
52 Collaborators Kanazawa Institute of echnology Dr. Isao Hashimoto Laboratory of Physiological Sciences Dr. Keiji Sakuma Communications Research Lab. Dr. ouru Miyauchi Communications Research Lab. Dr. Norio Fujimaki Siraume University Dr. Ryosuke akino University of Utah Dr. Srikantan S. Nagarajan M. I.. Dr. Alec Marantz University of Maryland Dr. David Poeppel University of California Dr. im Roberts University of okyo Dr. Yasushi Miyashita
FORSCHUNGSZENTRUM JÜLICH GmbH Zentralinstitut für Angewandte Mathematik D Jülich, Tel. (02461)
FORSCHUNGSZENTRUM JÜLICH GmbH Zentralinstitut für Angewandte Mathematik D-52425 Jülich, Tel. (2461) 61-642 Interner Bericht Temporal and Spatial Prewhitening of Multi-Channel MEG Data Roland Beucker, Heidi
More informationEstimating Neural Sources from Each Time-Frequency Component of Magnetoencephalographic Data
642 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 47, NO. 5, MAY 2000 Estimating Neural Sources from Each Time-Frequency Component of Magnetoencephalographic Data Kensuke Sekihara*, Member, IEEE, Srikantan
More informationALTHOUGH the single equivalent current dipole (ECD)
IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL 44, NO 9, SEPTEMBER 1997 839 Noise Covariance Incorporated MEG-MUSIC Algorithm: A Method for Multiple-Dipole Estimation Tolerant of the Influence of Background
More informationPARTICLE FILTERING, BEAMFORMING AND MULTIPLE SIGNAL CLASSIFICATION FOR THE ANALYSIS OF MAGNETOENCEPHALOGRAPHY TIME SERIES: A COMPARISON OF ALGORITHMS
Volume X, No. X, X, X XX Web site: http://www.aimsciences.org PARTICLE FILTERING, BEAMFORMING AND MULTIPLE SIGNAL CLASSIFICATION FOR THE ANALYSIS OF MAGNETOENCEPHALOGRAPHY TIME SERIES: A COMPARISON OF
More informationEEG/MEG Inverse Solution Driven by fmri
EEG/MEG Inverse Solution Driven by fmri Yaroslav Halchenko CS @ NJIT 1 Functional Brain Imaging EEG ElectroEncephaloGram MEG MagnetoEncephaloGram fmri Functional Magnetic Resonance Imaging others 2 Functional
More informationMaximum-Likelihood Estimation of Low-Rank Signals for Multiepoch MEG/EEG Analysis
IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 51, NO. 11, NOVEMBER 2004 1981 Maximum-Likelihood Estimation of Low-Rank Signals for Multiepoch MEG/EEG Analysis Boris V. Baryshnikov*, Barry D. Van Veen,
More informationBEAMFORMING DETECTORS WITH SUBSPACE SIDE INFORMATION. Andrew Bolstad, Barry Van Veen, Rob Nowak
BEAMFORMING DETECTORS WITH SUBSPACE SIDE INFORMATION Andrew Bolstad, Barry Van Veen, Rob Nowak University of Wisconsin - Madison 45 Engineering Drive Madison, WI 5376-69 akbolstad@wisc.edu, vanveen@engr.wisc.edu,
More informationRecipes for the Linear Analysis of EEG and applications
Recipes for the Linear Analysis of EEG and applications Paul Sajda Department of Biomedical Engineering Columbia University Can we read the brain non-invasively and in real-time? decoder 1001110 if YES
More informationSpatial Source Filtering. Outline EEG/ERP. ERPs) Event-related Potentials (ERPs( EEG
Integration of /MEG/fMRI Vince D. Calhoun, Ph.D. Director, Image Analysis & MR Research The MIND Institute Outline fmri/ data Three Approaches to integration/fusion Prediction Constraints 2 nd Level Fusion
More informationM/EEG source analysis
Jérémie Mattout Lyon Neuroscience Research Center Will it ever happen that mathematicians will know enough about the physiology of the brain, and neurophysiologists enough of mathematical discovery, for
More informationA Rao-Blackwellized particle filter for magnetoencephalography
A Rao-Blackwellized particle filter for magnetoencephalography C Campi Dipartimento di Matematica, Università di Genova, via Dodecaneso 35, 16146 Genova, Italy A Pascarella Dipartimento di Matematica,
More informationThe General Linear Model (GLM)
he General Linear Model (GLM) Klaas Enno Stephan ranslational Neuromodeling Unit (NU) Institute for Biomedical Engineering University of Zurich & EH Zurich Wellcome rust Centre for Neuroimaging Institute
More informationA Quantum Mechanical Model for the Vibration and Rotation of Molecules. Rigid Rotor
A Quantum Mechanical Model for the Vibration and Rotation of Molecules Harmonic Oscillator Rigid Rotor Degrees of Freedom Translation: quantum mechanical model is particle in box or free particle. A molecule
More information840 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 53, NO. 5, MAY 2006
840 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 53, NO. 5, MAY 2006 Performance Analysis of Reduced-Rank Beamformers for Estimating Dipole Source Signals Using EEG/MEG David Gutiérrez, Member, IEEE,
More informationPaired MEG Data Set Source Localization Using Recursively Applied and Projected (RAP) MUSIC
1248 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 47, NO. 9, SEPTEMBER 2000 Paired MEG Data Set Source Localization Using Recursively Applied and Projected (RAP) MUSIC John J. Ermer, Member, IEEE,
More informationAzimuth-elevation direction finding, using one fourcomponent acoustic vector-sensor spread spatially along a straight line
Volume 23 http://acousticalsociety.org/ 169th Meeting of the Acoustical Society of America Pittsburgh, Pennsylvania 18-22 May 2015 Signal Processing in Acoustics: Paper 4aSP4 Azimuth-elevation direction
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More informationBeamforming Techniques Applied in EEG Source Analysis
Beamforming Techniques Applied in EEG Source Analysis G. Van Hoey 1,, R. Van de Walle 1, B. Vanrumste 1,, M. D Havé,I.Lemahieu 1 and P. Boon 1 Department of Electronics and Information Systems, University
More informationNeuroscience Introduction
Neuroscience Introduction The brain As humans, we can identify galaxies light years away, we can study particles smaller than an atom. But we still haven t unlocked the mystery of the three pounds of matter
More informationA hierarchical Krylov-Bayes iterative inverse solver for MEG with physiological preconditioning
A hierarchical Krylov-Bayes iterative inverse solver for MEG with physiological preconditioning D Calvetti 1 A Pascarella 3 F Pitolli 2 E Somersalo 1 B Vantaggi 2 1 Case Western Reserve University Department
More informationTransformation of Whole-Head MEG Recordings Between Different Sensor Positions
Transformation of Whole-Head MEG Recordings Between Different Sensor Positions Transformation von Ganzkopf-MEG-Messungen zwischen verschiedenen Sensorpositionen Thomas R. Knösche Max Planck Institute of
More informationM. Holschneider 1 A. Eicker 2 R. Schachtschneider 1 T. Mayer-Guerr 2 K. Ilk 2
SPP project TREGMAT: Tailored Gravity Field Models for Mass Distributions and Mass Transport Phenomena in the Earth System by by M. 1 A. Eicker 2 R. Schachtschneider 1 T. Mayer-Guerr 2 K. Ilk 2 1 University
More informationPart 2: Multivariate fmri analysis using a sparsifying spatio-temporal prior
Chalmers Machine Learning Summer School Approximate message passing and biomedicine Part 2: Multivariate fmri analysis using a sparsifying spatio-temporal prior Tom Heskes joint work with Marcel van Gerven
More informationSolutions: Homework 5
Ex. 5.1: Capacitor Solutions: Homework 5 (a) Consider a parallel plate capacitor with large circular plates, radius a, a distance d apart, with a d. Choose cylindrical coordinates (r,φ,z) and let the z
More informationSpatial Harmonic Analysis of EEG Data
Spatial Harmonic Analysis of EEG Data Uwe Graichen Institute of Biomedical Engineering and Informatics Ilmenau University of Technology Singapore, 8/11/2012 Outline 1 Motivation 2 Introduction 3 Material
More informationHST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS
HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS EMERY N. BROWN AND CHRIS LONG NEUROSCIENCE STATISTICS RESEARCH LABORATORY DEPARTMENT
More informationNew Machine Learning Methods for Neuroimaging
New Machine Learning Methods for Neuroimaging Gatsby Computational Neuroscience Unit University College London, UK Dept of Computer Science University of Helsinki, Finland Outline Resting-state networks
More informationNeural mass model parameter identification for MEG/EEG
Neural mass model parameter identification for MEG/EEG Jan Kybic a, Olivier Faugeras b, Maureen Clerc b, Théo Papadopoulo b a Center for Machine Perception, Faculty of Electrical Engineering, Czech Technical
More informationTutorial on Blind Source Separation and Independent Component Analysis
Tutorial on Blind Source Separation and Independent Component Analysis Lucas Parra Adaptive Image & Signal Processing Group Sarnoff Corporation February 09, 2002 Linear Mixtures... problem statement...
More informationRobust regression and non-linear kernel methods for characterization of neuronal response functions from limited data
Robust regression and non-linear kernel methods for characterization of neuronal response functions from limited data Maneesh Sahani Gatsby Computational Neuroscience Unit University College, London Jennifer
More informationSpatial sound. Lecture 8: EE E6820: Speech & Audio Processing & Recognition. Columbia University Dept. of Electrical Engineering
EE E6820: Speech & Audio Processing & Recognition Lecture 8: Spatial sound 1 Spatial acoustics 2 Binaural perception 3 Synthesizing spatial audio 4 Extracting spatial sounds Dan Ellis
More informationTHE PROBLEM of localizing the sources of event related. Recursive MUSIC: A Framework for EEG and MEG Source Localization
1342 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 45, NO. 11, NOVEMBER 1998 Recursive MUSIC: A Framework for EEG and MEG Source Localization John C. Mosher,* Member, IEEE, and Richard M. Leahy, Member,
More informationWhat is NIRS? First-Level Statistical Models 5/18/18
First-Level Statistical Models Theodore Huppert, PhD (huppertt@upmc.edu) University of Pittsburgh Departments of Radiology and Bioengineering What is NIRS? Light Intensity SO 2 and Heart Rate 2 1 5/18/18
More informationREVIEW OF ORIGINAL SEMBLANCE CRITERION SUMMARY
Semblance Criterion Modification to Incorporate Signal Energy Threshold Sandip Bose*, Henri-Pierre Valero and Alain Dumont, Schlumberger Oilfield Services SUMMARY The semblance criterion widely used for
More informationFirst Technical Course, European Centre for Soft Computing, Mieres, Spain. 4th July 2011
First Technical Course, European Centre for Soft Computing, Mieres, Spain. 4th July 2011 Linear Given probabilities p(a), p(b), and the joint probability p(a, B), we can write the conditional probabilities
More information11.1 Three-Dimensional Coordinate System
11.1 Three-Dimensional Coordinate System In three dimensions, a point has three coordinates: (x,y,z). The normal orientation of the x, y, and z-axes is shown below. The three axes divide the region into
More informationOn Linearly Constrained Minimum Variance Beamforming
Journal of Machine Learning Research 99- Submitted /; Published / On Linearly Constrained Minimum Variance Beamforming Jian Zhang Chao Liu School of Mathematics, Statistics and Actuarial Science University
More informationSimultaneous Multi-frame MAP Super-Resolution Video Enhancement using Spatio-temporal Priors
Simultaneous Multi-frame MAP Super-Resolution Video Enhancement using Spatio-temporal Priors Sean Borman and Robert L. Stevenson Department of Electrical Engineering, University of Notre Dame Notre Dame,
More informationDirectional Sources and Beamforming
156 th meeting of the Acoustical Society of America Miami, Florida, November 13, 2008 Directional Sources and Beamforming Christian Bouchard Christian.Bouchard@nrc-cnrc.gc.ca David I. Havelock Institute
More informationIn: W. von der Linden, V. Dose, R. Fischer and R. Preuss (eds.), Maximum Entropy and Bayesian Methods, Munich 1998, Dordrecht. Kluwer, pp
In: W. von der Linden, V. Dose, R. Fischer and R. Preuss (eds.), Maximum Entropy and Bayesian Methods, Munich 1998, Dordrecht. Kluwer, pp. 17-6. CONVERGENT BAYESIAN FORMULATIONS OF BLIND SOURCE SEPARATION
More informationPerformance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna
Performance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna Kwan Hyeong Lee Dept. Electriacal Electronic & Communicaton, Daejin University, 1007 Ho Guk ro, Pochen,Gyeonggi,
More informationRobust Capon Beamforming
Robust Capon Beamforming Yi Jiang Petre Stoica Zhisong Wang Jian Li University of Florida Uppsala University University of Florida University of Florida March 11, 2003 ASAP Workshop 2003 1 Outline Standard
More informationLight Propagation in Free Space
Intro Light Propagation in Free Space Helmholtz Equation 1-D Propagation Plane waves Plane wave propagation Light Propagation in Free Space 3-D Propagation Spherical Waves Huygen s Principle Each point
More informationAdaptive Array Detection, Estimation and Beamforming
Adaptive Array Detection, Estimation and Beamforming Christ D. Richmond Workshop on Stochastic Eigen-Analysis and its Applications 3:30pm, Monday, July 10th 2006 C. D. Richmond-1 *This work was sponsored
More informationCOMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION. Aditya K. Jagannatham and Bhaskar D.
COMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION Aditya K Jagannatham and Bhaskar D Rao University of California, SanDiego 9500 Gilman Drive, La Jolla, CA 92093-0407
More informationADAPTIVE ANTENNAS. SPATIAL BF
ADAPTIVE ANTENNAS SPATIAL BF 1 1-Spatial reference BF -Spatial reference beamforming may not use of embedded training sequences. Instead, the directions of arrival (DoA) of the impinging waves are used
More informationSingular value decomposition. If only the first p singular values are nonzero we write. U T o U p =0
Singular value decomposition If only the first p singular values are nonzero we write G =[U p U o ] " Sp 0 0 0 # [V p V o ] T U p represents the first p columns of U U o represents the last N-p columns
More informationExercises. Chapter 1. of τ approx that produces the most accurate estimate for this firing pattern.
1 Exercises Chapter 1 1. Generate spike sequences with a constant firing rate r 0 using a Poisson spike generator. Then, add a refractory period to the model by allowing the firing rate r(t) to depend
More informationSignal Processing for Functional Brain Imaging: General Linear Model (2)
Signal Processing for Functional Brain Imaging: General Linear Model (2) Maria Giulia Preti, Dimitri Van De Ville Medical Image Processing Lab, EPFL/UniGE http://miplab.epfl.ch/teaching/micro-513/ March
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationIndependent Component Analysis (ICA) Bhaskar D Rao University of California, San Diego
Independent Component Analysis (ICA) Bhaskar D Rao University of California, San Diego Email: brao@ucsdedu References 1 Hyvarinen, A, Karhunen, J, & Oja, E (2004) Independent component analysis (Vol 46)
More informationwissen leben WWU Münster
MÜNSTER Sparsity Constraints in Bayesian Inversion Inverse Days conference in Jyväskylä, Finland. Felix Lucka 19.12.2012 MÜNSTER 2 /41 Sparsity Constraints in Inverse Problems Current trend in high dimensional
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More information32 +( 2) ( 4) ( 2)
Math 241 Exam 1 Sample 2 Solutions 1. (a) If ā = 3î 2ĵ+1ˆk and b = 4î+0ĵ 2ˆk, find the sine and cosine of the angle θ between [10 pts] ā and b. We know that ā b = ā b cosθ and so cosθ = ā b ā b = (3)(
More informationReal Plasma with n, T ~ p Equilibrium: p = j B
Real Plasma with n, T ~ p Equilibrium: p = j B B lines must lie in isobaric surfaces. Since B = 0, only possible if isobaric surfaces are topological tori. Magnetic field lines must form nested tori. Equilibrium
More informationCovariance smoothing and consistent Wiener filtering for artifact reduction in audio source separation
Covariance smoothing and consistent Wiener filtering for artifact reduction in audio source separation Emmanuel Vincent METISS Team Inria Rennes - Bretagne Atlantique E. Vincent (Inria) Artifact reduction
More informationStatistical Approaches to the Inverse Problem
50 Statistical Approaches to the Inverse Problem A. Pascarella 1 and A. Sorrentino 2 1 Dipartimento di Neuroscienze, Universitá di Parma 2 Department of Statistics, University of Warwick 1 Italy 2 United
More informationSOUND SOURCE LOCALIZATION VIA ELASTIC NET REGULARIZATION
BeBeC-2014-02 SOUND SOURCE LOCALIZATION VIA ELASTIC NET REGULARIZATION Xiaodong Li, Weiming Tong and Min Jiang School of Energy and Power engineering, Beihang University 100191, Beijing, China ABSTRACT
More informationApplications of Robust Optimization in Signal Processing: Beamforming and Power Control Fall 2012
Applications of Robust Optimization in Signal Processing: Beamforg and Power Control Fall 2012 Instructor: Farid Alizadeh Scribe: Shunqiao Sun 12/09/2012 1 Overview In this presentation, we study the applications
More informationROBUST ADAPTIVE BEAMFORMING BASED ON CO- VARIANCE MATRIX RECONSTRUCTION FOR LOOK DIRECTION MISMATCH
Progress In Electromagnetics Research Letters, Vol. 25, 37 46, 2011 ROBUST ADAPTIVE BEAMFORMING BASED ON CO- VARIANCE MATRIX RECONSTRUCTION FOR LOOK DIRECTION MISMATCH R. Mallipeddi 1, J. P. Lie 2, S.
More informationClassical Description of NMR Parameters: The Bloch Equations
Classical Description of NMR Parameters: The Bloch Equations Pascale Legault Département de Biochimie Université de Montréal 1 Outline 1) Classical Behavior of Magnetic Nuclei: The Bloch Equation 2) Precession
More informationOutlines: (June 11, 1996) Instructor:
Magnetic Resonance Imaging (June 11, 1996) Instructor: Tai-huang Huang Institute of Biomedical Sciences Academia Sinica Tel. (02) 2652-3036; Fax. (02) 2788-7641 E. mail: bmthh@ibms.sinica.edu.tw Reference:
More informationINTRODUCTION Noise is present in many situations of daily life for ex: Microphones will record noise and speech. Goal: Reconstruct original signal Wie
WIENER FILTERING Presented by N.Srikanth(Y8104060), M.Manikanta PhaniKumar(Y8104031). INDIAN INSTITUTE OF TECHNOLOGY KANPUR Electrical Engineering dept. INTRODUCTION Noise is present in many situations
More information7a3 2. (c) πa 3 (d) πa 3 (e) πa3
1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin
More informationArray signal processing for radio-astronomy imaging and future radio telescopes
Array signal processing for radio-astronomy imaging and future radio telescopes School of Engineering Bar-Ilan University Amir Leshem Faculty of EEMCS Delft University of Technology Joint wor with Alle-Jan
More informationSPATIAL FILTERING OF MAGNETOENCEPHALOGRAPHIC DATA IN SPHERICAL HARMONICS DOMAIN
SPATIAL FILTERING OF MAGNETOENCEPHALOGRAPHIC DATA IN SPHERICAL HARMONICS DOMAIN by Tolga Esat Özkurt B.Sc in E.C.E., Istanbul Technical University, 2002 M.Sc in C.S., Istanbul Technical University, 2004
More informationAnalysis of Empirical Bayesian Methods for Neuroelectromagnetic Source Localization
Analysis of Empirical Bayesian Methods for Neuroelectromagnetic Source Localization David Wipf 1, Rey Ramírez 2, Jason Palmer 1,2, Scott Makeig 2, & Bhaskar Rao 1 1 Signal Processing and Intelligent Systems
More informationLinear Optimum Filtering: Statement
Ch2: Wiener Filters Optimal filters for stationary stochastic models are reviewed and derived in this presentation. Contents: Linear optimal filtering Principle of orthogonality Minimum mean squared error
More informationErrata for Robot Vision
Errata for Robot Vision This is a list of known nontrivial bugs in Robot Vision 1986) by B.K.P. Horn, MIT Press, Cambridge, MA ISBN 0-6-08159-8 and McGraw-Hill, New York, NY ISBN 0-07-030349-5. Thanks
More information[POLS 8500] Review of Linear Algebra, Probability and Information Theory
[POLS 8500] Review of Linear Algebra, Probability and Information Theory Professor Jason Anastasopoulos ljanastas@uga.edu January 12, 2017 For today... Basic linear algebra. Basic probability. Programming
More informationChapter 1. Vector Analysis
Chapter 1. Vector Analysis Hayt; 8/31/2009; 1-1 1.1 Scalars and Vectors Scalar : Vector: A quantity represented by a single real number Distance, time, temperature, voltage, etc Magnitude and direction
More informationAN ITERATIVE DISTRIBUTED SOURCE METHOD FOR THE DIVERGENCE OF SOURCE CURRENT IN EEG INVERSE PROBLEM
J. KSIAM Vol.12, No.3, 191 199, 2008 AN ITERATIVE DISTRIBUTED SOURCE METHOD FOR THE DIVERGENCE OF SOURCE CURRENT IN EEG INVERSE PROBLEM JONGHO CHOI 1, CHNAG-OCK LEE 2, AND HYUN-KYO JUNG 1 1 SCHOOL OF ELECTRICAL
More informationThe Hebb rule Neurons that fire together wire together.
Unsupervised learning The Hebb rule Neurons that fire together wire together. PCA RF development with PCA Classical Conditioning and Hebbʼs rule Ear A Nose B Tongue When an axon in cell A is near enough
More informationNumerical Methods. Rafał Zdunek Underdetermined problems (2h.) Applications) (FOCUSS, M-FOCUSS,
Numerical Methods Rafał Zdunek Underdetermined problems (h.) (FOCUSS, M-FOCUSS, M Applications) Introduction Solutions to underdetermined linear systems, Morphological constraints, FOCUSS algorithm, M-FOCUSS
More informationDynamic Causal Modelling for EEG and MEG. Stefan Kiebel
Dynamic Causal Modelling for EEG and MEG Stefan Kiebel Overview 1 M/EEG analysis 2 Dynamic Causal Modelling Motivation 3 Dynamic Causal Modelling Generative model 4 Bayesian inference 5 Applications Overview
More informationMATH 0350 PRACTICE FINAL FALL 2017 SAMUEL S. WATSON. a c. b c.
MATH 35 PRACTICE FINAL FALL 17 SAMUEL S. WATSON Problem 1 Verify that if a and b are nonzero vectors, the vector c = a b + b a bisects the angle between a and b. The cosine of the angle between a and c
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationDynamical MEG Source Modeling with Multi-Target Bayesian Filtering
r Human Brain Mapping 000:000 000 (009) r Dynamical MEG Source Modeling with Multi-Target Bayesian Filtering Alberto Sorrentino, 1 * Lauri Parkkonen,,3 Annalisa Pascarella, 1,4 Cristina Campi, 1,5 and
More informationBayesian analysis of the neuromagnetic inverse
Bayesian analysis of the neuromagnetic inverse problem with l p -norm priors Toni Auranen, a, Aapo Nummenmaa, a Matti S. Hämäläinen, b Iiro P. Jääskeläinen, a,b Jouko Lampinen, a Aki Vehtari, a and Mikko
More informationA spatiotemporal dynamic distributed solution to the MEG inverse problem
A spatiotemporal dynamic distributed solution to the MEG inverse problem The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationRigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY)
Rigid Body Rotation Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content Introduction Angular Velocity Angular Momentum
More informationTechniques to Estimate Brain Connectivity from Measurements with Low Spatial Resolution
Techniques to Estimate Brain Connectivity from Measurements with Low Spatial Resolution 1. What is coherence? 2. The problem of volume conduction 3. Recent developments G. Nolte Dept. of Neurophysiology
More informationRobust Range-rate Estimation of Passive Narrowband Sources in Shallow Water
Robust Range-rate Estimation of Passive Narrowband Sources in Shallow Water p. 1/23 Robust Range-rate Estimation of Passive Narrowband Sources in Shallow Water Hailiang Tao and Jeffrey Krolik Department
More informationDifferential Equations
Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed
More informationPassive Sonar Detection Performance Prediction of a Moving Source in an Uncertain Environment
Acoustical Society of America Meeting Fall 2005 Passive Sonar Detection Performance Prediction of a Moving Source in an Uncertain Environment Vivek Varadarajan and Jeffrey Krolik Duke University Department
More informationContents. Introduction The General Linear Model. General Linear Linear Model Model. The General Linear Model, Part I. «Take home» message
DISCOS SPM course, CRC, Liège, 2009 Contents The General Linear Model, Part I Introduction The General Linear Model Data & model Design matrix Parameter estimates & interpretation Simple contrast «Take
More informationBeamforming. A brief introduction. Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University
Beamforming A brief introduction Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University March 2008 References Barry D. Van Veen and Kevin Buckley, Beamforming:
More informationExtraction of Sleep-Spindles from the Electroencephalogram (EEG)
Extraction of Sleep-Spindles from the Electroencephalogram (EEG) Allan Kardec Barros Bio-mimetic Control Research Center, RIKEN, 2271-13 Anagahora, Shimoshidami, Moriyama-ku, Nagoya 463, Japan Roman Rosipal,
More informationHST.582J/6.555J/16.456J
Blind Source Separation: PCA & ICA HST.582J/6.555J/16.456J Gari D. Clifford gari [at] mit. edu http://www.mit.edu/~gari G. D. Clifford 2005-2009 What is BSS? Assume an observation (signal) is a linear
More informationErrata for Robot Vision
Errata for Robot Vision This is a list of known nontrivial bugs in Robot Vision 1986) by B.K.P. Horn, MIT Press, Cambridge, MA ISBN 0-6-08159-8 and McGraw-Hill, New York, NY ISBN 0-07-030349-5. If you
More informationLarge brain effective network from EEG/MEG data and dmr information
Large brain effective network from EEG/MEG data and dmr information Brahim Belaoucha, Théodore Papadopoulo To cite this version: Brahim Belaoucha, Théodore Papadopoulo Large brain effective network from
More informationHomogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different Galilean Reference Systems
Progress In Electromagnetics Research Symposium Proceedings, Cambridge, USA, July 5 8, 2010 489 Homogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different
More informationMRI beyond Fourier Encoding: From array detection to higher-order field dynamics
MRI beyond Fourier Encoding: From array detection to higher-order field dynamics K. Pruessmann Institute for Biomedical Engineering ETH Zurich and University of Zurich Parallel MRI Signal sample: m γκ,
More informationConsider the following spike trains from two different neurons N1 and N2:
About synchrony and oscillations So far, our discussions have assumed that we are either observing a single neuron at a, or that neurons fire independent of each other. This assumption may be correct in
More informationarxiv: v1 [stat.ap] 29 Sep 2014
Sequential Monte Carlo Samplers for Semi Linear Inverse Problems and Application to Magnetoencephalography arxiv:1409.8109v1 [stat.ap] 29 Sep 2014 1. Introduction Sara Sommariva and Alberto Sorrentino
More informationClassical Description of NMR Parameters: The Bloch Equations
Classical Description of NMR Parameters: The Bloch Equations Pascale Legault Département de Biochimie Université de Montréal 1 Outline 1) Classical Behavior of Magnetic Nuclei: The Bloch Equation 2) Precession
More informationA Statistical Kirchhoff Model for EM Scattering from Gaussian Rough Surface
Progress In Electromagnetics Research Symposium 2005, Hangzhou, China, August 22-26 187 A Statistical Kirchhoff Model for EM Scattering from Gaussian Rough Surface Yang Du 1, Tao Xu 1, Yingliang Luo 1,
More informationPolarization analysis on three-component seismic data
Bollettino di Geofisica Teorica ed Applicata Vol., n. -, pp. 9-; Sep.-Dec. Polarization analysis on three-component seismic data P. Capizzi, L. De Luca and M. Vitale Department of Chemistry and Physics
More informationMulti-Baseline SAR interferometry
ulti-baseline SAR interferometry A. onti Guarnieri, S. Tebaldini Dipartimento di Elettronica e Informazione - Politecnico di ilano G. Fornaro, A. Pauciullo IREA-CNR ESA cat-1 project, 3173. Page 1 Classical
More information