[發(fā)明專利]一種基于噪聲相關和小波閾值法的磁共振信號降噪方法

在審

• 申請?zhí)枺?/td>: 202010400933.5
• 申請日: 2020-05-13
• 公開（公告）號: CN111650653A
• 公開（公告）日: 2020-09-11
• 發(fā)明（設計）人: 林婷婷;于思佳;張揚;萬玲
• 申請（專利權）人: 吉林大學
• 主分類號: G01V3/38
• 分類號: G01V3/38;G01V3/14
• 代理公司: 沈陽銘揚聯(lián)創(chuàng)知識產權代理事務所(普通合伙) 21241
• 代理人: 屈芳

1.一種基于噪聲相關和小波閾值法的磁共振信號降噪方法，其特征在于，包括以下步驟：

a、使用核磁共振地下水探測儀同時收錄探測線圈和參考線圈采集的數(shù)據(jù)，包括探測線圈采集的空采噪聲s(t)和含噪磁共振信號s_x(t)，參考線圈同步采集的參考噪聲數(shù)據(jù)r(t)和參考噪聲數(shù)據(jù)r_x(t)，t為采樣時間；

b、根據(jù)小波閾值法的原理，對參考線圈采集的數(shù)據(jù)r_x(t)進行處理，濾除其中含有的有效磁共振信號r_s(t)，得到純凈的噪聲數(shù)據(jù)r_N(t)＝r_x(t)-r_s(t)；

c、依據(jù)空采噪聲在探測線圈和參考線圈中的相關關系h，計算探測線圈中的噪聲數(shù)據(jù)s_N(t)＝h·r_N(t)，并從含噪信號中去除，得到消噪后的磁共振信號s_c(t)，s_c(t)＝s_x(t)-s_N(t)。

2.按照權利要求1所述的方法，其特征在于，所述步驟b中的小波閾值法的具體步驟為：

1)、對參考線圈采集的數(shù)據(jù)r_x(t)進行小波分解，提取出各層細節(jié)系數(shù)d₁,…,d_k以及最后一層近似系數(shù)a_k，k為小波分解層數(shù)；

2)、對各層細節(jié)系數(shù)進行硬閾值處理，將大于閾值的細節(jié)系數(shù)置零，小于閾值的細節(jié)系數(shù)保持不變；

3)、利用最后一層小波近似系數(shù)和閾值處理后的各層細節(jié)系數(shù)重構參考通道中采集到的有效磁共振信號r_s(t)。

3.按照權利要求1所述的方法，其特征在于，所述步驟c中的相關關系h的具體計算步驟為：

探測線圈和參考線圈中空采噪聲的相關關系為rs為r(t)和s(t)的互相關系數(shù)，rr為r(t)的自相關系數(shù)，h為傳遞系數(shù)，通過矩陣形式相關關系h表示為h＝(r^Tr)^-1·r^Ts，T為矩陣轉置符號。