面部和手势提取

一时失言乱红尘 2022-06-06 10:55 85阅读 0赞

clc,clear,close all;
    f=imread('hhh.jpg');
    %rgb转换为yCbCr
    y=rgb2ycbcr(f);
    %提取Cr
    Cr=y(:,:,3);
    figure,imshow(Cr);
    %肤色范围140-160 肤色范围赋值为255（白色） 其他赋值为0（黑色）
    Cr(Cr>140&Cr<160)=255;
    Cr(Cr~=255)=0;
    figure,imshow(Cr);
    Cr=im2bw(Cr);
    % bwareaopen(BW,P)
    % 删除二值图像BW中面积小于P的对象并返回修改后的二值图
    Cr=bwareaopen(Cr,200);
    %填充二值图中的坑洞并返回修改后的二值图
    Cr=imfill(Cr,'holes');
    figure,imshow(Cr);
    [h,w]=size(Cr);
    a=1;b=1;
    hs=sum(Cr);
    figure,plot(hs);
    %循环条件 小于图片宽度
    while(a<w)
    %寻找临界条件 垂直投影的列的值大于0
    while(hs(a)==0&&a<w)
        a=a+1;
    end;
    b=a;
    while(hs(b)>0)
        b=b+1;
    end;
    %排除最后一次循环a=b输出的一条线
    if(a<b)
    C1=Cr(:,a:b);
    %[hss,wss]=find(C1)返回矩阵C1中非零元素的行和列的索引值。
    [hss,wss]=find(C1);
    c=min(hss);
    d=max(hss);
    %f(c:d,a:b,)二维图 从上到下，从左到右
    C1=C1(c:d,:);
    figure,imshow(C1);
    C2=cat(3,C1,C1,C1);
    %result=f(c:d,a:b,:)三维图 从上到下，从左到右
    result=f(c:d,a:b,:);
    %同类型才能执行运算 类型强转
    C2=uint8(C2).*result;
    figure,imshow(result);
    figure,imshow(C2);
    end;
    a=b;
    end;

执行y=rgb2ycbcr(f);Cr=y(:,:,3);显示的结果

![这里写图片描述][SouthEast]  
转换为二值图并把肤色转换为白色 其他变为黑色  
![这里写图片描述][SouthEast 1]  
去除小面积的白圈Cr=bwareaopen(Cr,200);  
![这里写图片描述][SouthEast 2]  
填充脸部的坑洞Cr=imfill(Cr,’holes’);  
![这里写图片描述][SouthEast 3]  
显示垂直投影直方图hs=sum(Cr);figure,plot(hs);  
![这里写图片描述][SouthEast 4]  
经过while循环提取面部和手势  
![这里写图片描述][SouthEast 5]

[SouthEast]: /images/20220606/2ad4a3f17d534f50bb9c42e586149782.png
[SouthEast 1]: /images/20220606/a7b091c437c14cedb771c64ba06138e7.png
[SouthEast 2]: /images/20220606/a4169d5906cf4219a2af94f6feb49181.png
[SouthEast 3]: /images/20220606/939aa3488da3478596d1356150448846.png
[SouthEast 4]: /images/20220606/376e861daefe46d8952096361b12f87f.png
[SouthEast 5]: /images/20220606/edff75576fef47febbc0a23fca523177.png