【正文】
。,12)。,39。title(39。image(a2)。axis square。fontsize39。第一次消噪圖像 39。colormap(map)。 ? subplot(2,2,3)。,12)。,39。title(39。 ? subplot(2,2,2)。,12)。,39。title(39。image(X)。,2)。,c,s,39。%利用小波分解中的第一層低頻系數(shù)進(jìn)行重構(gòu),即實(shí)現(xiàn)低通濾波消噪 ? a2 = wrcoef2(39。 ? [c,s] = wavedec2(X,2,‘db1’)。gaussian39。)。 圖象去噪的關(guān)鍵問(wèn)題是既要去除(或減弱)噪聲所對(duì)應(yīng)的 高頻成分,又要保留(或增強(qiáng))邊緣所對(duì)應(yīng)的高頻成分 小波變換應(yīng)用于圖像去噪 Step1 對(duì)含噪信號(hào)進(jìn)行小波分解,獲取近似序列&細(xì)節(jié)序列 Step2 對(duì)細(xì)節(jié)系數(shù)進(jìn)行處理 Step3 利用 step1中獲取的近似序列& step2中獲取的處理后的細(xì)節(jié)序列進(jìn)行重構(gòu),獲取去噪后的有用信號(hào) ? clc ? clear all ? [X,map]=imread(39。如果 是 的MRA,則 是 的 MRA 2()LR 2 ()nLR22()LR 22()LR( , )xy?2jV2j j jV V V?? jV2jV2()LR22()LR分離變量方法 ???????????????)()(),()()(),()()(),()()(),()3()2()1(yxyxyxyxyxyxyxyx????????二維小波分解算法 ? 行處理:將 的每一行( n 取定值)看成是一個(gè)一維信號(hào),分別通過(guò)低通濾波器 和高通濾波器 ? 列處理:將上述結(jié)果的每一列當(dāng)成一維信號(hào),再次通過(guò)低通濾波器 和高通濾波器 ( , )ja m n{}g?{}h?{}g?{}h?二維小波分解數(shù)據(jù)總量保持不變 是 、 均為低頻的信號(hào)分量 x? y?是 為低頻、 為高頻的信號(hào)分量 是 、 均為高頻的信號(hào)分量 是 為高頻、 為低頻的信號(hào)分量 x?x?x?y?y?y?2D- DWT 二維小波分解算法 ? 在可分離變量的情況下,二維重構(gòu)算法也可通過(guò)行處理和列處理的兩個(gè)步驟進(jìn)行 小波變換在圖像去噪中的應(yīng)用 ? 數(shù)據(jù)在采集&傳輸?shù)倪^(guò)程中可能受到噪聲污染 ? 噪聲可用平穩(wěn) Gaussian隨機(jī)過(guò)程 來(lái)描述 ? 當(dāng)噪聲功率譜為常數(shù)時(shí),稱(chēng)為 Gaussian白噪聲。) 0 20 40 60 80 100 12021012S i g n a l0 20 40 60 80 100 120020406080A m p l i t u d e S p e c t r u mfigure1 figure2 0 20 40 60 80 100 12000 . 511 . 5L o w p a s s F i l t e r ( V0)0 20 40 60 80 100 12000 . 511 . 5H i g h p a s s F i l t e r ( W0)0 20 40 60 80 100 12021012L o w f r e q u e n c y C o m p o n e n t0 20 40 60 80 100 1201012H i g h f r e q u e n c y C o m p o n e n tfigure3 figure4 0 20 40 60 80 100 120020406080100A m p l i t u d e S p e c t r u m o f L o w f r e q u e n c y C o m p o n e n t0 20 40 60 80 100 120050100A m p l i t u d e S p e c t r u m o f H i g h f r e q u e n c y C o m p o n e n tfigure5 0 20 40 60 80 100 1201 0 . 500 . 51R e c o n s t r u c t e d L o w f r e q u e n c y S i g n a l0 20 40 60 80