于健.txt

自编程序

clear;
close all;
I=imread('D:\于健-驾驶中的眼睛定位方法研究\WL001.BMP');
figure;
imshow(I);
I2=rgb2ycbcr(I);%
I3=rgb2gray(I);
I4=I3;
figure;
imshow(I3);
Cb=I2(:,:,2);
Cr=I2(:,:,3);
for i=1:288
    for j=1:352
          if Cb(i,j)<110.08|Cb(i,j)>125.44
             I3(i,j)=0;
          end 
     end
end
for i=1:288 
    for j=1:352
        if Cr(i,j)<135.68|Cr(i,j)>166.4
          I3(i,j)=0; 
          end 
     end
end
figure;
imshow(I3);
A=imcrop(I4,[125 130 80 60]);
figure,imshow(A);
B=zeros(61,81);
for i=1:61
    for  j=1:81
        if i<61
            B(i,j)=double(A(i+1,j))-double(A(i,j));
        end
    end
end
C=zeros(61,1);
C=sum(B,2);
figure,plot(C);
[max1,n1]=max(C)
[max2,n2]=min(C)
D=imcrop(A,[1 22 80 8]);
figure,imshow(D);
B=zeros(9,81);
for i=1:9
    for  j=1:81
        if i<9
            B(i,j)=double(D(i+1,j))-double(D(i,j));
        end
    end
end
C=zeros(81,1);
C=sum(D,1);
figure,plot(C);
[min1,n1]=min(C)
E=imcrop(D,[16 1 6 8]);
figure,imshow(E);