[發(fā)明專利]光聲圖像去噪方法及裝置

在審

• 申請(qǐng)?zhí)枺?/td>: 201610130722.8
• 申請(qǐng)日: 2016-03-08
• 公開（公告）號(hào): CN105631830A
• 公開（公告）日: 2016-06-01
• 發(fā)明（設(shè)計(jì)）人: 孫明健;劉婷;王明華;李百祺;史雅慧;付穎;冷官冀
• 申請(qǐng)（專利權(quán)）人: 哈爾濱工業(yè)大學(xué)（威海）
• 主分類號(hào): G06T5/00
• 分類號(hào): G06T5/00
• 代理公司: 北京路浩知識(shí)產(chǎn)權(quán)代理有限公司 11002
• 代理人: 李相雨

1.一種光聲圖像去噪方法，其特征在于，所述方法包括以下步驟：

S1、輸入源光聲圖像矩陣，其中所述源光聲圖像矩陣中含有噪聲；

S2、利用低秩矩陣近似的方法對(duì)所述源光聲圖像矩陣進(jìn)行去噪處 理，得到目標(biāo)光聲圖像矩陣。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟S2中利 用低秩矩陣近似的方法去除所述源光聲圖像矩陣中的稀疏脈沖噪聲和 高斯噪聲。

3.根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述步驟S2 包括以下步驟：

S21、根據(jù)所述源光聲圖像矩陣建立去噪模型；

S22、利用所述低秩矩陣近似的方法求解所述去噪模型，得到去噪 后的所述目標(biāo)光聲圖像矩陣。

4.根據(jù)權(quán)利要求3所述的方法，其特征在于，所述去噪模型為：

min||X-L-S||F2]]>

s.t.rank(L)≤r,car(S)≤k

式中，X為所述源光聲圖像矩陣，L為所述目標(biāo)光聲圖像矩陣，S 為稀疏脈沖噪聲矩陣，rank(L)為所述目標(biāo)光聲圖像矩陣L的秩，car(S) 為所述稀疏脈沖噪聲矩陣包含的元素的個(gè)數(shù)，r為預(yù)定秩值，k為所述 稀疏脈沖噪聲矩陣包含的元素個(gè)數(shù)的預(yù)定值。

5.根據(jù)權(quán)利要求4所述的方法，其特征在于，所述步驟S22包括 以下步驟：

S221、設(shè)定所述預(yù)定秩值r、所述稀疏脈沖噪聲矩陣包含的元素個(gè) 數(shù)的預(yù)定值k、最大迭代次數(shù)t_max、迭代停止判別值ε；同時(shí)初始化迭 代次數(shù)t為零，初始化所述目標(biāo)光聲圖像矩陣L₀等于所述源光聲圖像 矩陣，初始化所述稀疏脈沖噪聲矩陣S₀為零；

S222、所述迭代次數(shù)t加1；

S223、構(gòu)建兩個(gè)隨機(jī)高斯矩陣A₁和A₂：

其中，m為所述源光聲圖像矩陣的行數(shù)，n為所述源光聲圖像矩陣 的列數(shù)；

S224、求取矩陣X-S_t-1的雙邊隨機(jī)映射值Y₁和Y₂：

Y₁＝(X-S_t-1)A₁

A₂＝Y(jié)₁

Y₂＝(X-S_t-1)^TA₂

S225、比較所述預(yù)定秩值r與的大小，若則進(jìn)行S226；否則回到S222；其中為矩陣的秩；

S226、更新所述目標(biāo)光聲圖像矩陣為L(zhǎng)t：

S227、更新所述稀疏脈沖噪聲矩陣為S_t：S_t←P_Ω(X-L_t)，其中P_Ω為 |X-L_t|前k個(gè)最大的不為0的元素；

S228、計(jì)算的值，判斷該值是否小于ε或者t≥t_max是否成立，若該值小于ε或者t≥t_max，則停止迭代，該時(shí)刻的目標(biāo)光聲 圖像矩陣L_t為最終的目標(biāo)光聲圖像矩陣，否則返回所述步驟S222。

6.一種光聲圖像去噪裝置，其特征在于，所述裝置包括源光聲圖 像矩陣獲取單元以及去噪處理單元；

所述源光聲圖像矩陣獲取單元用于輸入源光聲圖像矩陣，其中所 述源光聲圖像矩陣中含有噪聲；

所述去噪處理單元用于利用低秩矩陣近似的方法對(duì)所述源光聲圖 像矩陣進(jìn)行去噪處理，得到目標(biāo)光聲圖像矩陣。