脑电地形图绘制

（ⅱ）Time scheme
Ⅳ-ⅱ 特征提取——小波包熵
1.特征提取依据——大脑ERD 、ERS 特性
当人们实际做或仅想象单侧肢体运动时，大脑的感觉运动皮层会出现 EEG的节律性活动，在9HZ~13HZ（u节律，主要来自中枢后躯体感觉皮层） 和18HZ~22HZ（β节律，主要来自中央前运动皮层）两个频段的EEG信号幅 度将发生相应的改变。 ↓ ERD (event-related de-synchronization, ERD )事件相关去同步 ↑ ERS (event-related synchronization, ERS )事件相关同步
医学
石油勘探
应用领域
航空航天
气象
金融
Ⅱ 脑电地形图的意义
脑电地形图（BEAM）是一种集中表达大脑电生理信息的图形技术，能比较直观地 反映大脑神经活动的图形系统。
为了减轻临床医生的头脑中建立综合图像的 困难，工程人员致力于把现代计算机技 术和信号处理技术引入神经医学领域，把计算机的高速计算、高质量彩色图、易于操 作的屏幕控制结合起来构成了各种节律的脑电地形图。用直观的彩色图像（或灰度差 图像）取代了对多道原始EEG数据所包含的空间信息的表达。
pn En E
小波包熵，即是从小波包分解后的信号序列计算的一种熵值。根据香农熵的定义， 定义小波包熵为：
HIS色彩空间：从人的视觉系统出发，用色调（Hue）、饱和度 （Saturation）和亮度（Intensity）描述色彩。
优点：1.亮度分量和色度分量是分开的。 2.色调H和饱和度S的概念互相独立并与人的感知紧密相连。
解决方法：HIS转RGB模 型 1）当H在[0,120]之间：
B = I(1-S) R = I(1+ScosH/cos(60-H)) G = 3I-(B+R)
δ（1－3Hz）、θ（4－7Hz）、α（8－13Hz）、β（14－30Hz）
频段
频率下限
频率上限
δ
0
4.0
θ
4.0
8.0
α1
8.0
12.0
α1
12.0
16.0
β1
16.0
20.0
β2
20.0
30.0
Ⅲ-ⅰ 空间插值
由于采集的脑电信号仅来自于被测者头上有限固定的几个采集点，各点之间的空白 处需要采用一定的插值公式进行插值填充，这些插值的数据是依靠各采集点的功率大 小以及各采集点到需要进行插值计算点的距离而定，一般采用的插值公式为：
Ⅳ-ⅰ 离线分析数据
BCI competition 2003，University of Technology Graz, Austria
C3 CZ C4
+++
（ⅰ）Electrode locations
0
1
23
4
56
7
8
9 sec
Trigger Beep
Feedback period with cue
功率值 max
计算H值大小
（H，功率）
240
H（角度）
根据H值大小和公式计算相应 的RGB分量
RGB（R,G,B）
HIS转RGB模型测试
取功率为0-50线性增长，步长为一，按照上面的流程将其映射为蓝色到红色。
Ⅳ 离线脑电地形图
开发环境：Matlab 频率分段及功率表现方式：小波包+熵 插值方法：一般方式 显示方法：动画
2）当H在[120,240]之间：
R = I(1-S) G = I(1+Scos(H-120) /cos(180-H)) B = 3I-(G+R)
3）当H在[240,360]之间：
G = I(1-S) B = I(1+Scos(H-240) /cos(300-H)) R = 3I-(B+G)
获得头皮某点 的功率
脑电地形图绘制
报告内容
数据可视化 脑电地形图的意义 脑电地形图绘制的关键 离线地形图 实时地形图
Ⅰ 数据可视化
概念：借助图形化手段，将数据以图形图像形式表示，并利用 数据分析和开发工具发现其中未知信息的处理过程。
意义：提供象人眼一样的直觉的、交互的和反应灵敏的可视化 环境。
手段：图象、曲线、二维图形、三维体和动画等。
Electroencephalogram
BEAM
Ⅲ 脑电地形图绘制的关键
➢ 脑电信号采集（脑电仪）和预处理（滤波等）。 ➢ 计算各采集点各个频段（FFT、小波等）处功率谱强度。 ➢ 空间插值。（最近邻点、克里金、径向基函数插值法等） ➢ 功率值-彩色映射。（HIS模型转RGB）
Ⅲ-ⅰ 计算采集点处功ຫໍສະໝຸດ Baidu谱
Ⅲ-ⅰ 彩色映射
问题：
功率强度 0~max
彩色模式 蓝~红
RGB色彩模式：为图像中的每一个RGB分量分配一个0~255范围内的强度值。目前 的显示器大都采用这种模式。
VC中的像素点赋值函数：COLORREF SetPixel( int x, int y, COLORREF crColor ); crColor：RGB (red, green, blue)
8~12
20~24 16~20
mu rhythm
(8–13Hz)
β rhythm
(18–22Hz)
小波包节点能量可以有效表示信号的能量。子信号的能量可由该空间小波包系数的平方和计
算，如式：
En
D j,n 2 k
k
信号总能量即为每个子信号能量之和：
E En
n
则可以定义相对小波包能量为每个子空间信号能量和信号总能量之比：
f(X )a 1 //X X A A 2 2 b 1//X X B B 2 2 ...... 1 p //X X P P 22
式 中 X 为 所 需 进 行 插 值 计 算 的 点 的 位 置 ， a 、 b…p 代 表 各 采 集 点 的 功 率 值 ， XA 、 XB…XP为所需进行插值计算的点到各采集点的距离。
小波包分解结构
信号F(t)
S(2,1) 0~16HZ
S(1,1) 0~32HZ
S(2,2) 16~32HZ
S(2,4)
S(1,2)
S(2,3)
S(3,1) 0~8HZ
S(3,2) 8~16HZ
S(3,4) 24~32
S(3,3) 16~24
S(3,7)
S(3,8)
S(3,6)
S(3,5)
S(4,1) S(4,2) S(4,4) S(4,3) S(4,7) S(4,8) S(4,6) S(4,5) S(4,13) S(4,14) S(4,16) S(4,15) S(4,11) S(4,12) S(4,10) S(4,9)