能够对图像进行小波模极大值边缘检测,边缘跟踪,以及将图像分割出来

%小波模极大值边缘检测
clc;clear

% 下载图像
%load woman
file='E:\学习\论文\中期准备资料\region_seg\滤波+叠加.bmp';
I=imread(file);
%X=rgb2gray(I);
%image(X);%
%colormap(map);%
X=double(I);
a=min(size(X));
b=max(size(X));
%X=X(1:a,1:a);
X(a+1:b,1:b)=1;
SIZE=length(X); % 图像尺寸



% 多尺度
m=2.0;
delta=2^m;


% 构造高斯函数的偏导
N=20; % 滤波器长度（需要调整，必须是偶数）
A=-1/sqrt(2*pi); % 幅度

for index_x=1:N;
for index_y=1:N;
x=index_x-(N+1)/2;
y=index_y-(N+1)/2;
phi_x(index_x,index_y)=A*(x/delta^2).*exp(-(x.*x+y.*y)/(2*delta^2));
phi_y(index_x,index_y)=A*(y/delta^2).*exp(-(x.*x+y.*y)/(2*delta^2));
end
end;

phi_x=phi_x/norm(phi_x); % 能量归一化
phi_y=phi_y/norm(phi_y); % 能量归一化

% 对图象做行列卷积，求小波系数
Gx=conv2(X,phi_x,'same');
Gy=conv2(X,phi_y,'same');

% 求梯度
Grads=sqrt((Gx.*Gx)+(Gy.*Gy));

% 求幅角（梯度方向）
angle_array=zeros(SIZE,SIZE); % 角度

% 遍历
for i=1:SIZE
for j=1:SIZE
if (abs(Gx(i,j))>eps*100) % x的绝对值足够大
p=atan(Gy(i,j)/Gx(i,j))*180/pi; % 反正切求角度值(1,4象限)
if (p<0) % 负的幅角（4象限）
p=p+360;
end;
if (Gx(i,j)<0 && p>180) % 2象限的特殊处理
p=p-180;
elseif (Gx(i,j)<0 && p<180) % 3象限的特殊处理
p=p+180;
end
else % 90或270度
p=90;
end
angle_array(i,j)=p; % 幅角赋值
end
end;

% 找边缘
edge_array=zeros(SIZE,SIZE);

% 遍历
for i=2:SIZE-1
for j=2:SIZE-1
if ((angle_array(i,j)>=(-22.5) && angle_array(i,j)<=22.5) || ...
(angle_array(i,j)>=(180-22.5) && angle_array(i,j)<=(180+22.5))) % 0/180
if (Grads(i,j)>Grads(i+1,j) && Grads(i,j)>Grads(i-1,j))
edge_array(i,j)=Grads(i,j);
end
elseif ((angle_array(i,j)>=(90-22.5) && angle_array(i,j)<=(90+22.5)) || ...
(angle_array(i,j)>=(270-22.5) && angle_array(i,j)<=(270+22.5))) % 90/270
if (Grads(i,j)>Grads(i,j+1) && Grads(i,j)>Grads(i,j-1))
edge_array(i,j)=Grads(i,j);
end
elseif ((angle_array(i,j)>=(45-22.5) && angle_array(i,j)<=(45+22.5)) || ...
(angle_array(i,j)>=(225-22.5) && angle_array(i,j)<=(225+22.5))) % 45/225
if (Grads(i,j)>Grads(i+1,j+1) && Grads(i,j)>Grads(i-1,j-1))
edge_array(i,j)=Grads(i,j);
end
else % 135/215
if (Grads(i,j)>Grads(i+1,j-1) && Grads(i,j)>Grads(i-1,j+1))
edge_array(i,j)=Grads(i,j);
end
end
end
end

% 去除伪边缘
MAX_E=max(max(edge_array).'); % 最大幅度值
edge_array=edge_array/M