matlab x y矢量大小,请教 错误使用 zeros 大小矢量必须是包含实数元素的行矢量

clear;

[fname,pname]=uigetfile({'*.jpg';'*.bmp';'*.tif';'*.gif'},'Please choose a color picture...');% 返回打开的图片名与图片路径名

[u,v]=size(fname);

y=fname(v);% 图片格式代表值

switch y

case 0

errordlg('You Should Load Image File First...','Warning...');

case{'g';'G';'p';'P';'f';'F'};% 图片格式若是JPG/jpg、BMP/bmp、TIF/tif或者GIF/gif，才打开

I=cat(2,pname,fname);

Ori_Face=imread(I);

subplot(2,3,1),imshow(Ori_Face);

otherwise

errordlg('You Should Load Image File First...','Warning...');

end

p=rgb_RGB(Ori_Face);          %肤色分割

subplot(2,3,2);imshow(p);     % 显示皮肤灰度图像

low_pass=1/9*ones(3);

image_low=filter2(low_pass,p);% 低通滤波去噪声

subplot(2,3,3);imshow(image_low);

% 自适应阀值程序

previousSkin2=zeros(i,j);

changelist=[];

for threshold=0.55:-0.1:0.05

two_value=zeros(i,j);

two_value(find(image_low>threshold))=1;

change=sum(sum(two_value-previousSkin2));

changelist=[changelist change];

previousSkin2=two_value;

end

[C,I]=min(changelist);

optimalThreshold=(7-I)*0.1;

two_value=zeros(i,j);

two_value(find(image_low>optimalThreshold))=1;% 二值化

subplot(2,3,4);imshow(two_value);             % 显示二值图像

frontalmodel=imread('templet.jpg');           % 读入人脸模板照片

FaceCoord=[];

imsourcegray=rgb2gray(Ori_Face);             % 将原照片转换为灰度图像

[L,N]=bwlabel(two_value,8);% 标注二值图像中连接的部分,L为数据矩阵，N为颗粒的个数

for i=1:N,

[x,y]=find(bwlabel(two_value)==i);% 寻找矩阵中标号为i的行和列的下标

bwsegment=bwselect(two_value,y,x,8);% 选择出第i个颗粒

numholes=1-bweuler(bwsegment,4);% 计算此区域的空洞数

if(numholes>=1)% 若此区域至少包含一个洞，则将其选出进行下一步运算

RectCoord=-1;

[m,n]=size(bwsegment);

[cx,cy]=center(bwsegment);% 求此区域的质心

bwnohole=bwfill(bwsegment,'holes');% 将洞封住(将灰度值赋为1)

justface=uint8(double(bwnohole).*double(imsourcegray));% 只在原照片的灰度图像中保留该待识别区域

angle=orient(bwsegment,cx,cy);% 求此区域的偏转角度

bw=imrotate(bwsegment,angle,'bilinear');

bw=bwfill(bw,'holes');

[l,r,u,d]=bianjie(bw);

wx=(r-l+1);% 宽度

ly=(d-u+1);% 高度

wratio=ly/wx;% 高宽比

if((0.8<=wratio)&&(wratio<=2))% 如果目标区域的高度/宽度比例大于0.8且小于2.0，则将其选出进行下一步运算

S=ly*wx;% 计算包含此区域矩形的面积

A=bwarea(bwsegment);% 计算此区域面积

if(A/S>0.35)

[ccorr,mfit,RectCoord]=mobanpipei(justface,frontalmodel,ly,wx,cx,cy,angle);

end

if(ccorr>=0.6)

mfitbw=(mfit>=1);

invbw=xor(mfitbw,ones(size(mfitbw)));

source_with_hole=uint8(double(invbw).*double(imsourcegray));

final_image=uint8(double(source_with_hole)+double(mfit));

subplot(2,3,5);imshow(final_image);% 显示覆盖了模板脸的灰度图像

imsourcegray=final_image;

subplot(2,3,6);imshow(Ori_Face);% 显示检测效果图

end;

if(RectCoord~=-1)

FaceCoord=[FaceCoord;RectCoord];

end

end

end

end

% 在认为是人脸的区域画矩形

[numfaces x]=size(FaceCoord);

for i=1:numfaces,

hd=rectangle('Position',FaceCoord(i,:));

set(hd,'edgecolor','y');

end