[發(fā)明專利]一種三維超聲圖像非局部均值濾波方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310428718.6
• 申請(qǐng)日: 2013-09-17
• 公開（公告）號(hào): CN103544682A
• 公開（公告）日: 2014-01-29
• 發(fā)明（設(shè)計(jì)）人: 侯文廣;李柳;王學(xué)文;翁飛;張旭明;丁明躍
• 申請(qǐng)（專利權(quán)）人: 華中科技大學(xué)
• 主分類號(hào): G06T5/00
• 分類號(hào): G06T5/00;A61B8/00
• 代理公司: 華中科技大學(xué)專利中心 42201
• 代理人: 朱仁玲

1.一種三維超聲圖像非局部均值濾波方法，其特征在于，包括以下步驟：

（1）載入待濾波的目標(biāo)三維超聲圖像信息；

（2）將目標(biāo)三維超聲圖像拆分為多個(gè)基準(zhǔn)塊；

（3）對(duì)于每一個(gè)基準(zhǔn)塊，分別在目標(biāo)三維超聲圖像中確定其相應(yīng)的搜索區(qū)域；

（4）對(duì)于每一個(gè)基準(zhǔn)塊，逐體素遍歷該基準(zhǔn)塊相應(yīng)搜索區(qū)域中與基準(zhǔn)塊大小相同的圖像塊，將這些圖像塊作為該基準(zhǔn)塊的相似塊，計(jì)算這些圖像塊與所述基準(zhǔn)塊之間的相似度；

（5）對(duì)于每一個(gè)基準(zhǔn)塊，將其所有相似塊的加權(quán)平均值作為該基準(zhǔn)塊的濾波結(jié)果，其中相似塊的權(quán)重為其與基準(zhǔn)塊的相似度；

（6）將所有基準(zhǔn)塊濾波結(jié)果整合，得到目標(biāo)三維超聲圖像最終濾波結(jié)果并輸出。

2.如權(quán)利要求1所述的三維超聲圖像非局部均值濾波方法，其特征在于，所述相似塊與基準(zhǔn)塊間相似度根據(jù)伽馬模型模擬貝葉斯分布計(jì)算，相似度計(jì)算公式為：

w(Bi,Bj)=1ziexp(-1h2Σp=1P(u(Bi)(p)-u(Bj)(p))2u(Bj)(p)),]]>

其中，B_i為中心是第i體素的基準(zhǔn)塊，B_j為中心是第j體素的相似塊，w(B_i,B_j)為相似度，Z_i為歸一化系數(shù),h為濾波因子,P為圖像塊B_i或B_j的大小，為基準(zhǔn)塊內(nèi)第p個(gè)體素點(diǎn)的圖像灰度值，為相似塊內(nèi)第p個(gè)體素點(diǎn)的圖像灰度值。

3.如權(quán)利要求1或2所述的三維超聲圖像非局部均值濾波方法，其特征在于，所述基準(zhǔn)塊為邊長(zhǎng)在5到11體素的三維立方體圖像。

4.如權(quán)利要求1至3中任意一項(xiàng)所述的三維超聲圖像非局部均值濾波方法，其特征在于，所述相鄰基準(zhǔn)塊中心之間的距離為3到6體素。

5.如權(quán)利要求1至4中任意一項(xiàng)所述的三維超聲圖像非局部均值濾波方法，其特征在于，所述搜索區(qū)域?yàn)橐运龌鶞?zhǔn)塊中心為中心的三維立方體圖像區(qū)域。

6.如權(quán)利要求1至5中任意一項(xiàng)所述的三維超聲圖像非局部均值濾波方法，其特征在于，所述搜索區(qū)域邊長(zhǎng)為11體素。

7.如權(quán)利要求1到6中任意一項(xiàng)所述的三維超聲圖像非局部均值濾波方法，其特征在于，所述相鄰基準(zhǔn)塊其中心之間的距離小于等于基準(zhǔn)塊邊長(zhǎng)。