[發(fā)明專利]一種基于結(jié)構(gòu)相似性的腦磁MEG源定位方法

有效

• 申請(qǐng)?zhí)枺?/td>: 202210048425.4
• 申請(qǐng)日: 2022-01-17
• 公開（公告）號(hào): CN114065825B
• 公開（公告）日: 2022-04-19
• 發(fā)明（設(shè)計(jì)）人: 高陽;李文;寧曉琳;安楠
• 申請(qǐng)（專利權(quán)）人: 北京航空航天大學(xué)杭州創(chuàng)新研究院
• 主分類號(hào): G06K9/00
• 分類號(hào): G06K9/00;G06K9/62;G06V10/74;A61B5/245;A61B5/00
• 代理公司: 北京易捷勝知識(shí)產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11613
• 代理人: 齊勝杰

1.一種基于結(jié)構(gòu)相似性的腦磁MEG源定位方法，其特征在于，包括：

S1、針對(duì)預(yù)先獲取的腦磁MEG信號(hào)B進(jìn)行預(yù)處理，得到預(yù)處理后的腦磁MEG信號(hào)B；

所述預(yù)先獲取的腦磁MEG信號(hào)B為在磁屏蔽的環(huán)境下，采用測量腦磁信號(hào)的傳感器所采集的；

S2、基于預(yù)先采集的被試者的MRI數(shù)據(jù)構(gòu)建真實(shí)的頭部模型，并得出引導(dǎo)場矩陣L；

S3、針對(duì)所述預(yù)處理后的腦磁MEG信號(hào)B和所述引導(dǎo)場矩陣L，獲取多個(gè)初始源定位結(jié)果*S以及分別與多個(gè)初始源定位結(jié)果*S一一對(duì)應(yīng)的估計(jì)腦磁信號(hào)*B；

S4、基于多個(gè)初始源定位結(jié)果*S一一對(duì)應(yīng)的估計(jì)腦磁信號(hào)*B和所述預(yù)處理后的腦磁MEG信號(hào)B，確定最終源估計(jì)的值，并成像顯示；

所述S2具體包括：

S21、采用Freesurfer軟件對(duì)所述預(yù)先采集的受試者的MRI數(shù)據(jù)進(jìn)行分割，把受試者的MRI數(shù)據(jù)的頭部部分分割為：頭皮組織、顱骨、灰質(zhì)和白質(zhì)；

S22、采用openMEEG軟件設(shè)置所述灰質(zhì)和白質(zhì)、頭皮組織和顱骨的電導(dǎo)率，構(gòu)建受試者的頭部模型；

S23、采用高精度光學(xué)掃描儀將所述頭部模型所在的坐標(biāo)系和測量腦磁信號(hào)的傳感器所在的坐標(biāo)系進(jìn)行配準(zhǔn)得到配準(zhǔn)矩陣，并進(jìn)一步基于所述配準(zhǔn)矩陣和所述頭部模型采用matlab中fieldtrip工具包計(jì)算得出引導(dǎo)場矩陣L；

所述S3具體包括：

S31、分別采用預(yù)先設(shè)定的不同參數(shù)值下的偶極子擬合法針對(duì)所述預(yù)處理后的腦磁MEG信號(hào)B和所述引導(dǎo)場矩陣L進(jìn)行處理，分別獲取不同參數(shù)值所對(duì)應(yīng)的初始源定位結(jié)果*S；

所述參數(shù)值包括偶極子擬合法中的源的位置、方向和大小參數(shù)的值；

S32、基于所述初始源定位結(jié)果*S和所述引導(dǎo)場矩陣L，采用公式(1)，獲取與該初始源定位結(jié)果*S所對(duì)應(yīng)的估計(jì)腦磁信號(hào)*B；

所述公式(1)為：

*B＝L×*S；

所述S4包括：

S41、將采用公式(2)分別將所述估計(jì)腦磁信號(hào)*B和預(yù)處理后的腦磁MEG信號(hào)B中的任一值歸一化成0至255之間，獲取歸一化后的估計(jì)腦磁信號(hào)*B和預(yù)處理后的腦磁MEG信號(hào)B，并基于歸一化后的估計(jì)腦磁信號(hào)*B和預(yù)處理后的腦磁MEG信號(hào)B，獲取所述估計(jì)腦磁信號(hào)*B在灰度圖片格式下的矩陣數(shù)據(jù)和所述預(yù)處理后的腦磁MEG信號(hào)B在灰度圖片格式下的矩陣數(shù)據(jù)；

所述公式(2)為：

*A_i＝255×(A_i-a)/(b-a)；

A_i為所述估計(jì)腦磁信號(hào)*B或預(yù)處理后的腦磁MEG信號(hào)B中的任一值；

*A_i為所述估計(jì)腦磁信號(hào)*B或預(yù)處理后的腦磁MEG信號(hào)B中的任一值歸一化后所對(duì)應(yīng)的像素值；

a為所述估計(jì)腦磁信號(hào)*B或預(yù)處理后的腦磁MEG信號(hào)B中的最小值；

b為所述估計(jì)腦磁信號(hào)*B或預(yù)處理后的腦磁MEG信號(hào)B中的最大值；

S42、基于所述估計(jì)腦磁信號(hào)*B在灰度圖片格式下的矩陣數(shù)據(jù)和所述預(yù)處理后的腦磁MEG信號(hào)B在灰度圖片格式下的矩陣數(shù)據(jù)，獲取所述估計(jì)腦磁信號(hào)*B在灰度圖片格式下的矩陣數(shù)據(jù)和所述預(yù)處理后的腦磁MEG信號(hào)B在灰度圖片格式下的矩陣數(shù)據(jù)的相似度的值；

S43、基于每一所述估計(jì)腦磁信號(hào)*B在灰度圖片格式下的矩陣數(shù)據(jù)和所述預(yù)處理后的腦磁MEG信號(hào)B在灰度圖片格式下的矩陣數(shù)據(jù)的相似度的值，確定最終源估計(jì)的值，并成像顯示。