发明名称 基于灰建模的EEG-fMRI融合方法
摘要 本发明提出了一种基于灰建模的EEG-fMRI的融合方法,其步骤包括:(1)对原始EEG进行降采样及溯源,以便于对EEG和fMRI数据在空域和时域进行配准;(2)对EEG和fMRI数据进行预处理;(3)对数据进行累加生成,为建立灰色模型做必要的准备;(4)利用GM(1,1)模型提取特征,求取相应的参数;(5)对EEG和fMRI的特征参数进行融合;(6)将融合后的特征反投影回观测域得到最终的结果。本发明通过灰建模方法可以有效地融合EEG和fMRI的特点,得到同时具有高的时间分辨率和高空间分辨率的图像。
申请公布号 CN106355189A 申请公布日期 2017.01.25
申请号 CN201510420802.2 申请日期 2015.07.13
申请人 西北工业大学 发明人 谢松云;李亚兵;张筠鹏;张娟丽;侯文刚
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 代理人
主权项 基于灰建模的EEG‑fMRI融合方法,包括如下步骤:(1)对原始EEG进行降采样及溯源,对EEG和fMRI数据进行预处理,为融合做准备。EEG的预处理是指将EEG源的数值标准化;fMRI的预处理包含时间重排列及与sMRI的配准;(2)将原始数据进行累加生成,为下一步建立灰色模型做预处理;(3)利用灰色理论的GM(1,1)模型提取特征,求取相应的参数:发展系数a和灰作用量b<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><msup><mi>dx</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mi>dt</mi></mfrac><mo>+</mo><msup><mi>ax</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>=</mo><mi>b</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000119229490000011.GIF" wi="1070" he="139" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mfrac><mi>b</mi><mi>a</mi></mfrac><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>ak</mi></mrow></msup><mo>+</mo><mfrac><mi>b</mi><mi>a</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000119229490000012.GIF" wi="1233" he="126" /></maths>(4)分别将提取出来的EEG和fMRI特征进行融合;由于EEG和fMRI根据GM(1,1)提取出的特征进行融合,而非对原始图像进行融合,且经过溯源的EEG和fMRI数据量纲并不统一,无法利用常用的图像融合方法进行融合。本技术采用平均值的数学运算进行特征叠加。发展系数a反映了数据的发展趋势,灰作用量b反映了数据的变化关系。根据公式(1)和(2)可以看出,指数部分的解影响较大,为了减小解的震荡,故对发展系数a求均方根;b代表数据的整体偏离程度,故对灰作用量b求算术平均,从而获得融合结果的参数a和b:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>a</mi><mi>fusion</mi></msub><mo>=</mo><mo>-</mo><msqrt><mfrac><mrow><msup><msub><mi>a</mi><mi>EEG</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>a</mi><mi>fMRI</mi></msub><mn>2</mn></msup></mrow><mn>2</mn></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000119229490000013.GIF" wi="1154" he="178" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>b</mi><mi>fusion</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mi>EEG</mi></msub><mo>+</mo><msub><mi>b</mi><mi>fMRI</mi></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000119229490000014.GIF" wi="1095" he="137" /></maths>(5)将融合后的特征反投影回观测域以得到最终的融合结果。
地址 710072 陕西省西安市友谊西路127号