人脸定位代码

                                                            %-------------------------------------------------------------------------

clear all

x=imread('C:\wuzun.jpg');  %读取原始图像

figure(1);subplot(1,4,1);imshow(x);title('原图像1');

y=rgb2gray(x);

subplot(1,4,2);imshow(y);title('图像1的灰度图');           % 图1灰度图

u1=imnoise(y,'salt & pepper',0.07);

subplot(1,4,3);imshow(u1);title('图像1加噪声图');     %给图1加椒盐噪声

zz=medfilt2(u1,[3 3]);z=medfilt2(zz,[5 5]);%(2次中值滤波)

subplot(1,4,4);imshow(z);title('图像1中值滤波');  %图1中值滤波（3*3窗口）

% figure(6);subplot(1,2,1);imshow(zz);title('一次滤波');

% subplot(1,2,2);imshow(z);title('二次滤波');

x1=imread('C:\wuzun1.jpg');

figure(2);subplot(1,4,1);imshow(x1);title('原图像2');

y1=rgb2gray(x1);

subplot(1,4,2);imshow(y1);title('图像2灰度图');            %图2灰度图

u2=imnoise(y1,'salt & pepper',0.13);

subplot(1,4,3);imshow(u2);title('图像2叫噪声图');    %图2加椒盐噪声

zz1=medfilt2(u2,[3 3]);z1=medfilt2(zz1,[5 5]);%(2次中值滤波)

subplot(1,4,4);imshow(z1);title('图像2滤波图');  %图2中值滤波

% figure(5);subplot(1,2,1);imshow(zz1),title('一次滤波');

% subplot(1,2,2);imshow(z1);title('二次滤波');

% figure(3);subplot(1,2,1);imshow(z);title('处理后图1');

% figure(3);subplot(1,2,2);imshow(z1);title('处理后图2');

f=imsubtract(z1,z);                                               %图1图2差分

figure(4);subplot(1,1,1);imshow(f);title('差分后图像');

% agin=medfilt2(f,[3 3]);figure(10);imshow(agin);f=agin;   %差分图中值滤波

%-------------------------------------------------------------------------