[發明專利]一種針對高通量基因測序圖像的去噪方法

有效

• 申請號: 201510433962.0
• 申請日: 2015-07-22
• 公開（公告）號: CN105205788B
• 公開（公告）日: 2018-06-01
• 發明（設計）人: 顏珂;徐勇
• 申請（專利權）人: 哈爾濱工業大學深圳研究生院
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 深圳市科吉華烽知識產權事務所(普通合伙) 44248
• 代理人: 張英;韓云涵

1.一種針對高通量基因測序圖像的去噪方法，其包括以下步驟：

(1)對測序圖像使用小波函數進行小波分解，得到每一層小波系數ωi，具體為采用àtrous小波和它的β-3spline version對圖像進行小波分解，包括如下步驟：

a.初始化i＝0，則第0層的圖像為原圖M0；

b.變量i自加，圖像Mi-1每行每列都與一個一維的核h進行卷積，卷積后圖像表示為Mi，核h表示為矩陣并且在矩陣的元素之間插入(2i-1-1)個零；

c.計算每一層小波系數：ωi(k)＝Mi-1(k)-Mi(k)；

(2)針對每一層小波系數，計算出當前小波系數對應的全局閾值λi：

λ i = k m 1 * m 2 | | ω i - m e d i a n ( ω i ) | | 1 , ]]>

其中median(ωi)代表小波系數ωi對應的中值，m1,m2代表圖像的行和列，k代表系數；

(3)通過每一層小波系數ωi和對應的全局閾值λi，求出每一層小波系數對應的估計小波系數；

ω ~ i = 0 , | ω i | ≤ λ i sgn ( ω i ) ( | ω i | - r * λ i * α 1 | | ω i - λ i | | 2 ) , | ω i | > λ i , ]]>

其中，α＞1,r是調整因子，sgn(x)代表信號函數，當ωi大于0時，信號函數值為1；當ωi小于0，信號函數值為-1；

(4)使用小波重構算法得到去噪后的圖像：針對每一層估計小波系數得到去噪后的圖像，

M ( x , y ) = M N ( x , y ) + Σ i = 1 N ω ~ i ( x , y ) , ]]>

其中N代表最大小波分解層數，MN(x,y)代表經過小波分解算法后得到的低頻分量，代表每一層估計小波系數。