问题描述

MATLAB中如何沿边缘法线方向在边缘左右各取3个点并求得这些点的梯度。（已求得边缘法线方向）

clear all;
clc;

%%%%%%%%%%%%%%%%%%
%%%%读入图像%%%%
load('Tar0');
imshow(Im00);
title('origin image');
%I=imread('噪声原型.jpg');

%%%%%%%%%%%%%%%%%%
%%%%初始化%%%%%%%%
i=0;j=0;
num1=1;num2=5;
dir=0;

%%%%%%%%%%%%%%%%%%
%%%%求d(xy)值%%%%
f=Im00;
[m1n1]=size(f);
f1=[zeros(m11) f zeros(m11)];
f1=[zeros(1n1+2);f1;zeros(1n1+2)];
for i=1:m1
for j=1:n1
d(ij)=sqrt([f1(ij)-f1(i+2j)]^2+[f1(ij)-f1(ij+2)]^2);
end
end

%%%%%%%%%%%%%%%%%%
%%%%图像边缘%%%%%%%
I=edge(f'canny'[0.050.1]2.0);
figure(2);
imshow(I);
title('edge image');

%%%%%%%%%%%%%%%%%%
%%%%求四个方向值%%%
[m2n2]=size(I);
I=[zeros(m21) I zeros(m21)];
I=[zeros(1n2+2);I;zeros(1n2+2)];

for i=3:m2-2
for j=3:n2-2

    dir0=d(i-2j)+d(i-1j)+d(ij)+d(i+1j)+d(i+2j);    dir90=d(ij-2)+d(ij-1)+d(ij)+d(ij+1)+d(ij+2);    dir45=d(i-2j+2)+d(i-1j+1)+d(ij)+d(i+1j-1)+d(i+2j-2);    dir135=d(i-2j-2)+d(i-1j-1)+d(ij)+d(i+1j+1)+d(i+2j+2);    if(dir0>=dir90&&dir0>=dir45&&dir0>=dir135)        dir=dir0;    elseif (dir90>=dir0&&dir90>=dir45&&dir90>=dir135)        dir=dir90;    elseif (dir45>=dir0&&dir45>=dir90&&dir45>=dir135)        dir=dir45;    elseif (dir135>=dir0&&dir135>=dir90&&dir135>=dir45)        dir=dir135;    end        end

end

dir1=-1/dir;