课程设计报告基于肤色模型的人脸分割定位

1、目 录1. 前 言12. 设计目标12.1基于肤色的人脸标定系统设计目标12.2基于肤色的人脸标定系统设计原理13. 基于肤色的人脸标定系统实现算法13.1基于RGB彩色空间的肤色模型13.2基于YCbCr彩色空间的肤色模型14. 基于肤色的人脸标定系统实验结果14.1肤色提取算法实验14.2人脸标定算法实验14.2人眼标定算法实验15. 总结1参考文献24 / 5文档可自由编辑打印1. 前 言2. 设计目标2.1基于肤色的人脸标定系统设计目标2.2基于肤色的人脸标定系统设计原理3. 基于肤色的人脸标定系统实现算法3.1基于RGB彩色空间的肤色模型3.2基于YCbCr彩色空间的肤色模型4. 基

2、于肤色的人脸检测4.1肤色相似度计算4.2动态阈值确定4.2人脸区域边框生成5. 总结参考文献1 阮秋琦 著.数字图像处理学.北京：电子工业出版社，20014 沈凤麟 陈和晏生物医学随机信号处理合肥：中国科学技术大学出版社1999 6 Dimitris G. Manokakis,Vinay K. Ingle,Stephen M.Kogon著, 周正 等译.统计与自适应滤波.北京：电子工业出版社，20038 董长虹主编，赖志国，余啸海编著，MATLAB 图像处理与应用。国防工业出版社，2004 10 Miroslaw D. Lutovac , Dejan V . Tosic , Brian L

3、. Evans 著，朱义胜，董辉等译，信号处理滤波器设计基于MATLAB和Mathematica 设计方法。电子工业出版社，2004 13 胡广书编著，数字信号处理理论、算法与实现（第二版）。北京：清华大学出版社，2003 %1.选择图像dir *.jpg;str(Please input the :,s);str(strFile,.jpg);imgrgb=imread(strFile1);m,n,c=size(imgrgb);% delete strFile,c;%2.将RGB色彩空间转换为Ycbcr色彩空间imgrgb3=rgb2ycbcr(imgrgb);y=imgrgb3(:,:,1)

4、;cb=imgrgb3(:,:,2);cr=imgrgb3(:,:,3);%3.在Ycbcr色彩空间中分割肤色区域。cb=double(cb);cr=double(cr);y=double(y);bwycbcr=zeros(m,n);for i=1:m for j=1:n if y(i,j)188 cb1=108+(y(i,j)-188)*10/47; cr1=154+(y(i,j)-188)*22/47; wcb=14+(235-y(i,j)*32.97/47; wcr=10+(235-y(i,j)*28.76/47; cb1=(cb(i,j)-cb1)*46.97/wcb+cb1; cr1

5、=(cr(i,j)-cr1)*38.76/wcr+cr1; else cb1=cb(i,j); cr1=cr(i,j); end x1=-0.819 0.574*cb1-109.38;cr1-152.02; y1=-0.574 -0.819*cb1-109.38;cr1-152.02; if (x1-1.60).2/644.6521+(y1-2.41).2/196.8409=1 bwycbcr(i,j)=1; imgrgb4(i,j,1)=y(i,j); imgrgb4(i,j,2)=cb(i,j); imgrgb4(i,j,3)=cr(i,j); else bwycbcr(i,j)=0; i

6、mgrgb4(i,j,1)=0; imgrgb4(i,j,2)=0; imgrgb4(i,j,3)=0; end end end%strWrite=strcat(a,strFile1);%imwrite(bwycbcr,strWrite,jpg);% *%标定人脸%获取图像总尺寸m n = size(y);isTop = 0;Left_min = 1 n;Right_max=1 1;length_R=0;%逐行检查for i=1:m for j=1:n%逐列检查 %查到第一个最上点 if bwycbcr(i,j)=1 & isTop =0 top=i j; isTop =1; end %最左边

7、的点 if bwycbcr(i,j)=1 & j= Right_max(2) Right_max=i j; end %最低点 %elseif bwycbcr(i,j)=1 & endend%画矩形imgrgb5=imgrgb;%画上方横线length_R=Right_max(2)-Left_min(2)+1;imgrgb5(top(1),Left_min(2):Right_max(2),1)=255*ones(1,length_R);imgrgb5(top(1),Left_min(2):Right_max(2),2)=255*ones(1,length_R);imgrgb5(top(1),Left_min(2):Right_max(2),3)=255*ones(1,length_R);%画下方横线bottom=160;imgrgb5(90,Left_min(2):Right_max(2),1)=255*ones(1,length_R);imgrgb5(90,Left_min(2):Right_max(2),2)=255*ones(1,length_R);imgrgb5(90,Left_min(2):Right_max(2),3)=255