clc;
clear all; %清除所有数据
close all; %关闭所有程序
image=imread('D:\buildingGray2.bmp'); %读取实验图像
image=rgb2gray(image); %将图像转化为2维图像
figure,imshow(image),title('原来图像');
w1=1/3.*[1 1 1]; %创建水平方向上滤波的模板
w2=1/3.*[1;1;1]; %创建垂直方向上的滤波模板
w3=[0.2 0 0.2;0 0.2 0;0.2 0 0.2]; %创建对角线上的滤波模板
w4=[0 0.2 0;0.2 0.2 0.2; 0 0.2 0];
[a1,h1,v1,d1]=dwt2(image,'bior3.7'); %对图像进行单尺度分解
[a2,h2,v2,d2]=dwt2(a1,'bior3.7'); %对a1进行单尺度分解
[a3,h3,v3,d3]=dwt2(a2,'bior3.7'); %对a2进行单尺度分解
h33=imfilter(h3,w1); %对h向量进行水平方向模板进行滤波
v33=imfilter(v3,w2); %对v向量进行垂直方向模板进行滤波
d33=imfilter(d3,w4); %对d向量进行X方向进行滤波
x=idwt2(a3,h33,v33,d33,'bior3.7'); %重构出a2
x=x(1:131,1:30); %截取a2的一部分
h22=imfilter(h2,w1); %对h向量进行水平方向模板进行滤波
v22=imfilter(v2,w2); %对v向量进行垂直方向模板进行滤波
d22=imfilter(d2,w4); %对d向量进行X方向进行滤波
y=idwt2(x,h22,v22,d22,'bior3.7'); %重构出a1
y=y(1:247,1:349); %截取y的一部分
h11=imfilter(h1,w1); %对h向量进行水平方向模板进行滤波
v11=imfilter(v1,w2); %对v向量进行垂直方向模板进行滤波
d11=imfilter(d1,w4); %对d向量进行X方向进行滤波
z=idwt2(y,h11,v11,d11,'bior3.7'); %重构出图像
figure,imshow(z,[]); %将z归一化之后显示出来
评论0