小波去噪主程序.doc

小波去噪主程序 %xiaoboquzao.m %对彩色图像进行去噪 I = imread(菊花.png,png); % 读入图像 X = im2double(I); % 转换成双精度类型 x_noise = imnoise(X, gaussian, 0.01); % 加入高斯噪声 %提取三个通道信息 xr = x_noise(:, :, 1); % R通道 xg = x_noise(:, :, 2); % G通道 xb = x_noise(:, :, 3); % B通道 %估计三个通道的阈值 [Cr, Sr] = wavedec2(xr, 2, sym4); [Cg, Sg] = wavedec2(xg, 2, sym4); [Cb, Sb] = wavedec2(xb, 2, sym4); thr_r = Donoho(xr); % R通道全局阈值 thr_g = Donoho(xg); % G通道全局阈值 thr_b = Donoho(xb); % B通道全局阈值 %对三个通道分别进行去噪 % Donoho全局阈值 软阈值公式 x_soft_r = wdenoise(xr, gbl, s, thr_r, sym4, 2); x_soft_g = wdenoise(xg, gbl, s, thr_g, sym4, 2); x_soft_b = wdenoise(xb, gbl, s, thr_b, sym4, 2); % Donoho全局阈值 硬阈值公式---------------------------------------------- x_hard_r = wdenoise(xr, gbl, h, thr_r, sym4, 2); x_hard_g = wdenoise(xg, gbl, h, thr_g, sym4, 2); x_hard_b = wdenoise(xb, gbl, h, thr_b, sym4, 2); % Birge-Massart策略 软阈值公式---------------------------------------------- thr_lvd_r=thr_r x_soft_lvd_r = wdenoise(xr, lvd, s, thr_lvd_r, sym4, 1); thr_lvd_g=thr_g x_soft_lvd_g = wdenoise(xg, lvd, s, thr_lvd_g, sym4, 1); thr_lvd_b=thr_b x_soft_lvd_b = wdenoise(xb, lvd, s, thr_lvd_b, sym4, 1); %半软阈值--------------------------------------------------------------- x1_r = den1(xr, sym4, 2, thr_r); x1_g = den1(xg, sym4, 2, thr_g); x1_b = den1(xb, sym4, 2, thr_b); %半软阈值 + 均值滤波---------------------------------------------------- x1_5_r = den1_5_1(xr, sym4, 2, thr_r, 0.5*thr_r); x1_5_g = den1_5_1(xg, sym4, 2, thr_g, 0.5*thr_g); x1_5_b = den1_5_1(xb, sym4, 2, thr_b, 0.5*thr_b); %自适应阈值------------------------------------------------------------- x4_r = den4(xr, sym4, 2); x4_g = den4(xg, sym4, 2); x4_b = den4(xb, sym4, 2); % 恢复去噪后的图像 x_soft = cat(3, x_soft_r, x_soft_g, x_soft_b);