[發(fā)明專利]血樣水平依賴性功能磁共振信號震蕩頻率聚類分析方法

無效

• 申請?zhí)枺?/td>: 201410053308.2
• 申請日: 2014-02-17
• 公開（公告）號: CN103767707A
• 公開（公告）日: 2014-05-07
• 發(fā)明（設(shè)計）人: 宋瀟鵬;張毅;劉一軍
• 申請（專利權(quán)）人: 南通南大四維影通科技有限公司
• 主分類號: A61B5/055
• 分類號: A61B5/055;A61B5/145
• 代理公司: 北京世譽(yù)鑫誠專利代理事務(wù)所(普通合伙) 11368
• 代理人: 郭官厚

1.血樣水平依賴性功能磁共振信號震蕩頻率聚類分析方法，其特征在于:包括以下步驟：

步驟一：對從磁共振機(jī)器采集到的靜息態(tài)數(shù)據(jù)進(jìn)行預(yù)處理分析；

步驟二：采用希爾伯特黃變換對預(yù)處理后的信號進(jìn)行分解，得到不同頻率的固有模態(tài)函數(shù)；希爾伯特黃變換包括經(jīng)驗(yàn)?zāi)Ｊ椒纸夂拖柌刈儞Q兩步；

步驟三：計算不同頻率成分的量化指標(biāo)；

步驟四：以步驟三的量化指標(biāo)為分類依據(jù)，將腦內(nèi)不同體素聚類，統(tǒng)計分析以上不同指標(biāo)。

2.根據(jù)權(quán)利要求1所述的血樣水平依賴性功能磁共振信號震蕩頻率聚類分析方法，其特征在于：所述步驟一中，對從磁共振機(jī)器采集到的靜息態(tài)數(shù)據(jù)為BOLD-fMRI信號數(shù)據(jù)，對BOLD-fMRI信號數(shù)據(jù)進(jìn)行預(yù)處理分析的具體實(shí)現(xiàn)步驟為：

1）、利用基于MATLAB平臺的SPM軟件，對采集到的靜息態(tài)數(shù)據(jù)進(jìn)行時間差異矯正；

2）、對時間矯正過后的數(shù)據(jù)進(jìn)行頭動對齊；

3）、對頭動對齊后的數(shù)據(jù)進(jìn)行空間標(biāo)準(zhǔn)化；

4）、對空間標(biāo)準(zhǔn)化后的數(shù)據(jù)進(jìn)行去線性漂移和時域標(biāo)準(zhǔn)化。

3.根據(jù)權(quán)利要求1所述的血樣水平依賴性功能磁共振信號震蕩頻率聚類分析方法，其特征在于：所述步驟二中，采用希爾伯特黃變換對預(yù)處理后的信號進(jìn)行分解，得到不同頻率的成分，即固有模態(tài)函數(shù)，具體實(shí)現(xiàn)步驟為：

1）、采用經(jīng)驗(yàn)?zāi)Ｊ椒纸鈱π盘栠M(jìn)行分解，得到若干固有模態(tài)函數(shù)；

2）、對經(jīng)驗(yàn)?zāi)Ｊ椒纸夂蟮玫降墓逃心B(tài)函數(shù)進(jìn)行希爾伯特變換。

4.根據(jù)權(quán)利要求3所述的血樣水平依賴性功能磁共振信號震蕩頻率聚類分析方法，其特征在于：采用希爾伯特黃變換對預(yù)處理后的信號進(jìn)行分解，其中步驟1）中采用經(jīng)驗(yàn)?zāi)Ｊ椒纸鈱π盘栠M(jìn)行分解，得到若干固有模態(tài)函數(shù)的具體實(shí)現(xiàn)步驟為：

步驟1.1）、找到原始信號x(t)的所有極值點(diǎn)；

步驟1.2）、在所有極大值點(diǎn)之間插值得到x(t)的上包絡(luò)eup(t)，在所有極小值點(diǎn)之間插值得到x(t)的下包絡(luò)elow(t);

步驟1.3）、計算eup(t)和elow(t)的均值時間序列emean(t)=[elow(t)+eup(t)]/2；

步驟1.4）、x(t)減去均值時間序列得到一個震蕩模式：r(t)=x(t)-emean(t);

步驟1.5）、判斷斷震蕩模式r(t)是否符合固有模態(tài)函數(shù)的特征或是單調(diào)曲線，若符合IMF特征，則該“斷震蕩模式”成為一個固有模態(tài)函數(shù)，IMF_i(t)=r(t)；單調(diào)曲線則否r(t)成為殘差，EMD過程結(jié)束，否則設(shè)x(t)=r(t)重復(fù)以上步驟直到得到一個IMF；IMF特征為上下包絡(luò)對稱，IMF曲線過零點(diǎn)的次數(shù)與極值點(diǎn)個數(shù)相等或最多相差1個。