function [autocorr, startx, starty, mfit, bwmodel_rot] = isThisAFace(mult, model,ly,wx,cx, ...                                                  cy, angle)  model_rot = imrotate(model,-angle,'bilinear');  % Use the frontal model to do the comparison  % model_rot = imresize(model_rot,[ly wx],'bilinear');  [l,r,u,d] = recsize(model_rot);  model_rot=imcrop(model_rot,[l u (r-l) (d-u)]); % figure; % imshow(model_rot);  bwmodel_rot = zeros(size(model_rot));  bwmodel_rot(find(model_rot > 1))=1;  % Compute the center of the rotated face model  [modx,mody] =center(bwmodel_rot);  % Get the size of the rotated face model  [morig, norig] = size(bwmodel_rot);  startx = cx-modx;  starty = cy-mody;  endx = startx + norig-1;  endy = starty + morig-1;  mfit = zeros(size(mult));  mfitbw = mfit & 0;  [limy, limx] = size(mfit);  startx = checklimit(startx,limx);  starty = checklimit(starty,limy);  endx = checklimit(endx,limx);  endy = checklimit(endy,limy);  % -------------------------------------------------------------  % The following is to generate a new image having the same  % size as the original one, but with the face of the model on it  % rotated accordingly.  for i=starty:endy,   for j=startx:endx,     mfit(i,j) = model_rot(i-starty+1,j-startx+1);     mfitbw(i,j) = mfit(i,j);    end;  end; figure;  imshow(mfit,[0 255]);  % Get the autocorrelation value  autocorr = corr2(mfit,mult)% Verifies that the coordinate is between the image region.function newcoord = checklimit(coord,maxval)  newcoord = coord;  if (newcoord<1) newcoord=1; end;   if (newcoord>maxval) newcoord=maxval; end;% EOF.