matlab畢業(yè)設(shè)計(jì)外文翻譯--復(fù)雜脊波圖像去噪

1、<p><b>　　中文2990字</b></p><p><b>　　復(fù)雜脊波圖像去噪</b></p><p>　　作者：G. Y. Chen and B. Kegl  刊名：Pattern Recognition; 出版日期：2007</p><p><b>　　1.介

2、紹</b></p><p>　　小波變換已成功地應(yīng)用于many scientific fields such as image compression, i許多科學(xué)領(lǐng)域，如圖像壓縮，圖像age denoising, signal processing, com去噪，信號(hào)處理，計(jì)算機(jī)圖形，ics, and pattern recognition, to name only a feIC和模式識(shí)別，僅舉幾

3、例。Donoho和他的同事們提出了小波閾值去噪通過軟閾值和閾值.這種方法的出現(xiàn)對(duì)于大量的應(yīng)用程序是一個(gè)好的選擇。這是因?yàn)橐粋€(gè)小波變換能結(jié)合的能量,在一小部分的大型系數(shù)和大多數(shù)的小波系數(shù)中非常小,這樣他們可以設(shè)置為零。這個(gè)閾值的小波系數(shù)是可以做到的只有細(xì)節(jié)的小波分解子帶。我們有一些低頻波子帶不能碰觸，讓他們不閾值。眾所周知，Donoho提出的方法的優(yōu)勢(shì)是光滑和自適應(yīng)。然而，Coifman和Donoho指出，這種算法展示出一個(gè)視覺產(chǎn)出：吉布

4、斯現(xiàn)象在鄰近的間斷。因此，他們提出對(duì)這些產(chǎn)出去噪通過平均抑制所有循環(huán)信號(hào)。實(shí)驗(yàn)結(jié)果證實(shí)單目標(biāo)識(shí)別小波消噪優(yōu)于沒有目標(biāo)識(shí)別的情況。Bui和Chen擴(kuò)展了這個(gè)目標(biāo)識(shí)別計(jì)劃，他們發(fā)現(xiàn)多小波的目標(biāo)</p><p>　　研究脊波變換的數(shù)多年來打破了小波變換的局限性。將小波變換產(chǎn)生的二維圖像在每個(gè)規(guī)模大的小波系數(shù)的分解。有這么多的大系數(shù)，對(duì)于圖像去噪有很多困難。我們知道脊波變換已經(jīng)成功用于分析數(shù)字圖像。不像小波變換,脊波變換

5、過程首先計(jì)算積分在不同的方向和位置的數(shù)據(jù)。沿著“x1cos_ + x2sin_ = 常數(shù)” 一條線的脊波是不變的。在這些脊的方向正交小波變換是一。最近脊波已成功應(yīng)用于圖像去噪。在本文中，我們結(jié)合dual-tree complex wavelet in the ridgelet transfo二元樹復(fù)小波的脊波變換中并將其應(yīng)用到圖像降噪。這種近似二元樹性能的復(fù)雜變性小波和良好性能的脊波使我們有更好的方法去圖像去噪。實(shí)驗(yàn)結(jié)果表明，采用二元樹

6、復(fù)雜脊波在所有去噪圖像和許多不同噪音中我們的算法獲得較高的峰值信噪比（PSNR）。</p><p>　　這篇文章大體是這樣的。在第二部分，我們將解釋如何將二元樹復(fù)雜的波變換成脊波去圖像去噪。實(shí)驗(yàn)結(jié)果在第3節(jié)。第4節(jié)是最后得出的結(jié)論和未來需要做的工作。</p><p>　　2.用復(fù)雜脊波圖像去噪</p><p>　　離散脊波變換提供接近理想的稀松代表光滑的物體邊緣。高

7、斯去噪是一個(gè)接近最優(yōu)的方法。脊波變換能夠壓縮圖像能量成為少量的脊波系數(shù)。在另一方面,利用小波變換產(chǎn)生的多大的小波系數(shù)對(duì)每個(gè)尺度邊緣二維小波分解。這句話的意思是說許多小波系數(shù)進(jìn)行重構(gòu)在圖像的邊緣。我們知道近似氡轉(zhuǎn)化為數(shù)字?jǐn)?shù)據(jù)可以基于離散傅立葉變換。普通的脊波變換即可達(dá)到如下：</p><p>　　計(jì)算出二維FFT的圖像。</p><p>　　替補(bǔ)的采樣傅里葉廣域上變換得到晶格和極性格的采樣值

8、。</p><p>　　計(jì)算一維逆FFT每一個(gè)角的線。</p><p>　　執(zhí)行一維標(biāo)量小波對(duì)角線結(jié)果，獲取脊波系數(shù)。</p><p>　　眾所周知，普通的離散小波變換在變換期間是不移位和不轉(zhuǎn)變的。輸入信號(hào)的一個(gè)小小的改變能夠引起輸出小波系數(shù)很大的變化。為了克服這個(gè)問題，Kingsbury發(fā)明了一種新型的小波變換，叫做二元樹復(fù)雜小波變換，它能夠轉(zhuǎn)移性能和提高近似角分

9、辨率不變。由于標(biāo)量波不是轉(zhuǎn)移不變的，在脊波變換中就更好的應(yīng)用二元樹復(fù)雜小波變換這樣我們就可以叫我們的復(fù)雜脊波。這樣可以通過取代一維標(biāo)量小波的一維二元樹復(fù)雜小波在最后一步進(jìn)行脊波變換。用這種方法我們可以優(yōu)秀品質(zhì)的脊波變換用來替換二元樹發(fā)雜脊波。</p><p>　　這個(gè)復(fù)雜的脊波變換可以應(yīng)用到整體圖像，或者我們可以應(yīng)用到分割圖像大量重疊的平方或者在每一平方上運(yùn)用脊波變換。我們分解一組n*n的影像重疊順利進(jìn)入邊長R的

10、象素是重疊的是兩個(gè)相鄰長方形的數(shù)組大小為R/2*R兩者之間重疊的相鄰區(qū)域就是一個(gè)長方形的大小R*R/2。對(duì)于一個(gè)n*n的圖像，我們能夠計(jì)數(shù)2n=R對(duì)于不同方向的模塊，這個(gè)分區(qū)就引入了4倍的冗余。為了得到降噪的復(fù)雜脊波系數(shù)我們通常在當(dāng)前象素地位對(duì)降噪的復(fù)雜脊波系數(shù)進(jìn)行平均4份。復(fù)雜的脊波變換閾值類似于曲波閾值。當(dāng)我們求閾值時(shí)一個(gè)不同是我們采取的是復(fù)雜的脊波系數(shù)。當(dāng)yλ是帶噪的脊波系數(shù)。我們使用下列硬閾值規(guī)則估算未知的脊波系數(shù)。當(dāng)│yλ┃&

11、gt; kσ?, 我們令λ= ?λλ.否則, ^y_ = 0.在這里，?σ是通過用蒙特卡羅模擬接近。采用的系數(shù)k是依賴于噪聲系數(shù)。當(dāng)這個(gè)小于30時(shí)，我們用k=5首先分解尺度和k=4分解其他尺度。當(dāng)這個(gè)噪音系數(shù)大于30時(shí)，我們用k=6首次分解尺度和k=5分解其他尺度。這個(gè)復(fù)雜的脊波去噪算法能夠被描述如下：</p><p>　　圖像分割成R*R塊，兩個(gè)垂直相鄰的R /2*R重疊，兩個(gè)水平象素塊R*R/2重疊。<

12、/p><p>　　對(duì)于每一塊，應(yīng)用所提出的復(fù)雜脊波，復(fù)雜脊波系數(shù)的閾值，復(fù)雜脊波的逆換算。</p><p>　　在同一位置以平均象素對(duì)圖像去噪。</p><p>　　我們稱這種算法叫，同時(shí)我們使用普通的脊波。這個(gè)計(jì)算復(fù)雜度的ComRidgeletShrink是和小波RidgeletShrink的標(biāo)量相似。唯一的區(qū)別是我們?nèi)〈艘痪S小波變換與一維二元樹發(fā)雜小波變換。這個(gè)數(shù)

13、量的計(jì)算是一維二元樹復(fù)數(shù)小波的變換是一維小波變換的兩倍。該算法的其他計(jì)算步驟保持相同。我們的實(shí)驗(yàn)結(jié)果顯示ComRidgeletShrink優(yōu)于V isuShrink, RidgeletShink, and 過濾器wiener2等所有測(cè)試案例。在某些情況下，我們?cè)赗idgeletShink中能夠提高0.8db的信噪比。通過V isuShrink,能夠改善更大的去噪圖像。這表明ComRidgeletSrink對(duì)于自然圖像去噪是一個(gè)很好的選擇

14、。</p><p><b>　　3.實(shí)驗(yàn)結(jié)果</b></p><p>　　我們通過對(duì)眾所周知的蕾娜進(jìn)行處理，通過Donoho等人我們得到了這種圖片的自由軟體包WaveLab。帶有不同噪音的噪音圖像時(shí)通過對(duì)原無噪音圖像添加高斯白噪音得到的。與之相比，我們實(shí)行VisuShrink, RidgeletShrink, ComRidgeletShrink and wiener2

15、。VisuShrink是通用軟閾值去噪技術(shù)。這個(gè)wiener2函數(shù)是可以從MatLab圖像工具箱得到，我們用一個(gè)5*5的相鄰圖像在每個(gè)象素中。該wiener2適用于一個(gè)維納濾波器（一種線性的濾波器）圖形自適應(yīng)。剪裁自己的圖像局部方差。峰值信噪比的實(shí)驗(yàn)結(jié)果顯示的表1.我們發(fā)現(xiàn)對(duì)于分區(qū)塊的大小32*32或者64*64是最好的選擇。表1是對(duì)蕾娜圖像進(jìn)行去噪，根據(jù)不同的噪聲水平固定分區(qū)和一素塊為32*32。表格中的第一欄是原來帶噪圖片的信噪比，

16、其他列是通過不同去噪方法得到的去噪圖像信噪比。這個(gè)信噪比被定義PSNR = 􀀀10 log10Pi;j (B(i; j) 􀀀 A(j))2n22552 :; 其中B是去噪圖像A是無噪音圖像。從表1.我們可以看出VisuShrink ,ComRidgeletSh</p><p><b>　　4．結(jié)論和未來工作</b></p><p>

17、　　在這篇文章中我們研究了用復(fù)雜脊波對(duì)圖像去噪。復(fù)雜脊波變換是通過執(zhí)行一維二元樹復(fù)雜小波在空氣中氡的變換系數(shù)獲得的。氡變換是通過投影片定理得到的。對(duì)于圖像去噪近似轉(zhuǎn)換不變性質(zhì)的二元樹復(fù)數(shù)小波變換對(duì)于復(fù)雜小波變換是一個(gè)很好的選擇。復(fù)雜脊波變換提供了近乎完美的對(duì)于光滑的物體和表現(xiàn)對(duì)象與邊緣。這使噪音閾值的脊波系數(shù)更接近高斯白噪音的消噪。我們測(cè)試了我們新的去噪方法和幾個(gè)標(biāo)準(zhǔn)圖像和加入高斯白噪音的圖像。用一個(gè)非常簡單的硬閾值復(fù)雜脊波系數(shù)。在這些

18、脊波實(shí)驗(yàn)中，實(shí)驗(yàn)結(jié)果表明復(fù)雜的脊波有更好的去噪能力比起VisuShrink和普通的wiener2。我們建議ComRidgeletShrink用于實(shí)際的圖像去噪中。未來工作主要是考慮在復(fù)雜圖像應(yīng)用曲波復(fù)雜脊波。同樣，復(fù)雜脊波還可以應(yīng)用的不變特征提取模式識(shí)別方法。</p><p>　　Complex Ridgelets for Image Denoising</p><p>　　G. Y. C

19、hen and B. Kegl</p><p>　　1 Introduction</p><p>　　Wavelet transforms have been successfully used in many scientific fields such as image compression, image denoising, signal processing, computer

20、graphics,and pattern recognition, to name only a few.Donoho and his coworkers pioneered a wavelet denoising scheme by using soft thresholding and hard thresholding. This approach appears to be a good choice for a number

21、of applications. This is because a wavelet transform can compact the energy of the image to only a small number of large coefficients a</p><p>　　and Donoho pointed out, this algorithm exhibits visual artifac

22、ts: Gibbs phenomena in the neighbourhood of discontinuities. Therefore, they propose in a translation invariant (TI) denoising scheme to suppress such artifacts by averaging over the denoised signals of all circular shif

23、ts. The experimental results in confirm that single TI wavelet denoising performs better than the non-TI case. Bui and Chen extended this TI scheme to the multiwavelet case and they found that TI multiwavelet denoising&l

24、t;/p><p>　　along the lines x1cos_ + x2sin_ = constant. In the direction orthogonal to these ridges it is a wavelet.Ridgelets have been successfully applied in image denoising recently. In this paper, we combine

25、 the dual-tree complex wavelet in the ridgelet transform and apply it to image denoising. The approximate shift invariance property of the dual-tree complex wavelet and the good property of the ridgelet make our method a

26、 very good method for image denoising.Experimental results show that by using dua</p><p>　　complex wavelets into the ridgelet transform for image denoising. Experimental results are conducted in Section 3. F

27、inally we give the conclusion and future work to be done in section 4.</p><p>　　2 Image Denoising by using Complex</p><p>　　Ridgelets Discrete ridgelet transform provides near-ideal sparsity of

28、representation of both smooth objects and of objects with edges. It is a near-optimal method of denoising for Gaussian noise. The ridgelet transform can compress the energy of the image into a smaller number of ridgelet

29、coe_cients. On the other hand, the wavelet transform produces many large wavelet coe_cients on the edges on every scale of the 2D wavelet decomposition. This means that many wavelet coe_cients are needed in ord</p>

30、<p>　　1. Compute the 2D FFT of the image.</p><p>　　2. Substitute the sampled values of the Fourier transform obtained on the square lattice with sampled values on a polar lattice.</p><p>

31、　　3. Compute the 1D inverse FFT on each angular line.</p><p>　　4. Perform the 1D scalar wavelet transform on the resulting angular lines in order to obtain the ridgelet coe_cients.</p><p>　　It i

32、s well known that the ordinary discrete wavelet transform is not shift invariant because of the decimation operation during the transform. A small shift in the input signal can cause very di_erent output wavelet coe_cien

33、ts. In order to overcome this problem, Kingsbury introduced a new kind of wavelet transform, called the dual-tree complex wavelet transform, that exhibits approximate shift invariant property and improved angular resolut

34、ion. Since the scalar wavelet is not shift invariant, it</p><p>　　The complex ridgelet transform can be applied to the entire image or we can partition the image into a number of overlapping squares and we a

35、pply the ridgelet transform to each square. We decompose the original n _ n image into smoothly overlapping blocks of sidelength R pixels so that the overlap between two vertically adjacent blocks is a rectangular array

36、of size R=2 _ R and the overlap between two horizontally adjacent blocks is a rectangular array of size R _ R=2 . For an n _ n image, we co</p><p>　　The thresholding for the complex ridgelet transform is sim

37、ilar to the curvelet thresholding [10]. One difference is that we take the magnitude of the complex ridgelet coe_cients when we do the thresholding. Let y_ be the noisy ridgelet coe_cients. We use the following hard thre

38、sholding rule for estimating the unknown ridgelet coe_cients. When jy_j > k_~_, we let ^y_ = y_. Otherwise, ^y_ = 0. Here, ~It is approximated by using Monte-Carlo simulations. The constant k used is dependent on the

39、noise</p><p>　　The complex ridgelet image denoising algorithm can be described as follows:</p><p>　　1. Partition the image into R*R blocks with two vertically adjacent blocks overlapping R=2*R p

40、ixels and two horizontally adjacent blocks overlapping R _ R=2 pixels</p><p>　　2. For each block, Apply the proposed complex ridgelets, threshold the complex ridgelet coefficients, and perform inverse comple

41、x ridgelet transform.</p><p>　　3. Take the average of the denoising image pixel values at the same location.</p><p>　　We call this algorithm ComRidgeletShrink,while the algorithm using the ordin

42、ary ridgelets RidgeletShrink. The computational complexity of ComRidgeletShrink is similar to that of RidgeletShrink by using the scalar wavelets. The only di_erence is that we replaced the 1D wavelet transform with the

43、1D dual-tree complex wavelet transform. The amount of computation for the 1D dual-tree complex wavelet is twice that of the 1D scalar wavelet transform. However, other steps of the algorithm keep the same</p><

44、p>　　3 Experimental Results</p><p>　　We perform our experiments on the well-known image Lena. We get this image from the free software package WaveLab developed by Donoho et al. at Stanford University. Noi

45、sy images with di_erent noise levels are generated by adding Gaussian white noise to the original noise-free images. For comparison, we implement VisuShrink, RidgeletShrink, ComRidgeletShrink and wiener2. VisuShrink is t

46、he universal soft-thresholding denoising technique. The wiener2 function is available in the MATLAB Image Proces</p><p>　　4 Conclusions and Future Work</p><p>　　In this paper, we study image den

47、oising by using complex ridgelets. Our complex ridgelet transform is obtained by performing 1D dual-tree complex wavelet onto the Radon transform coe_cients. The Radon transform is done by means of the projection-slice t

48、heorem. The approximate shift invariant property of the dual-tree complex wavelet transform makes the complex ridgelet transform an excellent choice for image denoising. The complex ridgelet transform provides near-ideal