jiaodian.asv

人脸定位的MATLAB源程序

clear all
p=imread('pic.bmp');
a=edge(p);
[px,py]=size(a);
for i=1:px
    for j=1:py
        if j==1
            b(i,j)=0;
        elseif j==py
            b(i,j)=0;
        else
            b(i,j)=((a(i,j)-a(i,j-1))^2+(a(i,j)-a(i,j+1))^2)^3;
        end
    end
end
for i=1:py
    for j=1:px
        if j==1
           c(j,i)=0;
        elseif j==px
           c(j,i)=0;
        else
            c(j,i)=((a(j,i)-a(j-1,i))^2+(a(j,i)-a(j+1,i))^2)^3;
        end
    end
end
d=b+c;
for i=1:px
    for j=1:py
        if d(i,j)>=9
            e(i,j)=1;
        else
            e(i,j)=0;
        end
    end
end
figure,imshow(e);