小波边缘检测matlab程序,一个小波边缘检测的MATLAB程序

一个小波边缘检测的MATLAB程序

%load woman;

X=imread('Miss.bmp');

clear coef1;

clear coef2;

clear coef3;

clear form;

x=X;

x=double(x);%200 line,300 column;

line=250;

column=250;

J=6;

Tu=7;

smooth_area_num=20;

edge_area_num=40;

for i=1:line

coef1(i, :)=cwt(x(i,:),1,'bior4.4');

coef2(i,:)=cwt(x(i,:),2,'bior4.4');

coef3(i,:)=cwt(x(i,:),J,'bior4.4');

end

%正值1，负值0的矩阵temp表示(突变点左位置)

max=0;

max_loc=0;

max_num=15;

for i=1:line

j=1;

while j

if (coef3(i,j)>0 & coef3(i,j+1)<0) | (coef3(i,j)<0 & coef3(i,j+1)>0 )%找到左位置j

for k=j+1:column-1

if coef3(i,k)*coef3(i,j)>=0%找到右位置k

for m=j:k

if abs(coef3(i,m))>max

max=abs(coef3(i,m));

max_loc=m;

end

end

if max>max_num

form(i,max_loc)=1;

end

j=k;

max=0;

break;

end

end%for k

if k==column-1

j=column-1;

end

else

j=j+1;

end%if

end%while

end

%for i=1:line

%    for j=1:column

%        if abs(coef1(i,j))>=2*abs(coef2(i,j)) %是否加绝对值？

%           form(i,j)=0;%第一种方法去噪声

%        end

%    end

%end

%对列作小波变换

x=x';

for i=1:line

coef1(i,:)=cwt(x(i,:),1,'bior4.4');

coef2(i,:)=cwt(x(i,:),2,'bior4.4');

coef3(i,:)=cwt(x(i,:),J,'bior4.4');

end

%正值1，负值0的矩阵temp表示(突变点左位置)

max=0;

max_loc=0;

for i=1:line

j=1;

while j

if (coef3(i,j)>0 & coef3(i,j+1)<0) | (coef3(i,j)<0 & coef3(i,j+1)>0 )%找到左位置j

for k=j+1:column-1

if coef3(i,k)*coef3(i,j)>=0%找到右位置k

for m=j:k

if abs(coef3(i,m))>max

max=abs(coef3(i,m));

max_loc=m;

end

end

if max>max_num

form(max_loc,i)=1;

end

j=k;

max=0;

break;

end

end%for k

if k==column-1

j=column-1;

end

else

j=j+1;

end%if

end%while

end

%包括噪声在内的所有极值点①

%for i=1:line

%    for j=1:column

%        if abs(coef1(i,j))>=2*abs(coef2(i,j)) %是否加绝对值？!

%           form(i,j)=0;%第一种方法去噪声

%        end

%    end

%end

%此处添加用第二种方法去噪!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!② 结束

for i=2:line-2

for j=2:column-2

%链长为1

if form(i,j-1)==0 & form(i,j+1)==0 & form(i-1,j-1)==0 & form(i-1,j)==0 & form(i-1,j+1)==0 & form(i+1,j-1)==0 & form(i+1,j)==0 & form(i+1,j+1)==0

form(i,j)=0;%

end

%****************************************

%链长为2

%0000

%0110

%0000

if form(i,j)==1 & form(i,j+1)==1 & sum(sum(form(i-1:i+1,j-1:j+2)))==2

form(i,j)=0;

form(i,j+1)=0;

end

%000

%010

%010

%000

if form(i,j)==1 & form(i+1,j)==1 & sum(sum(form(i-1:i+2,j-1:j+1)))==2

form(i,j)=0;

form(i+1,j)=0;

end

%0000

%0100

%0010

%0000

if form(i,j)==1 & form(i+1,j+1)==1 & sum(sum(form(i:i+1,j:j+1)))==2

if sum(sum(form(i-1:i+2,j-1:j+2)))==3 & ( form(i-1,j+2)==1 | form(i+2,j-1)==1)

form(i,j)=0;

form(i+1,j+1)=0;

end

if sum(sum(form(i-1:i+2,j-1:j+2)))==4 &  form(i-1,j+2)==1 & form(i+2,j-1)==1

form(i,j)=0;

form(i+1,j+1)=0;

end

end

%0000

%0010

%0100

%0000

if form(i,j+1)==1 & form(i+1,j)==1 & sum(sum(form(i:i+1,j:j+1)))==2

if sum(sum(form(i-1:i+2,j-1:j+2)))==3 & ( form(i-1,j-1)==1 | form(i+2,j+2)==1)

form(i,j+1)=0;

form(i+1,j)=0;

end

if sum(sum(form(i-1:i+2,j-1:j+2)))==4 &  form(i-1,j-1)==1 & form(i+2,j+2)==1

form(i,j+1)=0;

form(i+1,j)=0;

end

end

end

end

%!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1!!

max=0;

for i=3:line-3

for j=3:column-3

s=sum(sum(form(i-2:i+2,j-2:j+2)));

for m=i-2:i+2

for n=j-2:j+2

if s>=7 & abs(coef3(m,n))<=s*3

form(m,n)=0;

end

end

end

end

end

end

%显示form即为边缘检测图像

[本帖最后由 mooni 于 2009-4-28 19:39 编辑]