小波去噪主程序%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); x...
