clear all
q=imread('pic.bmp'); %導入圖像
q1=imread('pic1.bmp');
for i=1:300
    for j=1:250
         rect=[i j 99 39];
         q2=imcrop(q,rect);
         r(i,j)=corr2(q1,q2);
         objvalue(i,j)=[(r(i,j)+1)/2]^2;
     end
 end
 z=objvalue(1,1);
 for i=1:300
     for j=1:250
         if objvalue(i,j)>z;
             z=objvalue(i,j);
             x=i;
             y=j;
         end
     end
 end
bx=[x x+99 x+99 x];
by=[y y y+109 y+109];
BW=roipoly(q,bx,by);
bx1=[x-1 x+100 x+100 x-1];
by1=[y-1 y-1 y+110 y+110];
BW1=roipoly(q,bx1,by1);
BW2=BW1-BW;
q3=double(q)/255;
Q=q3+BW2;
figure,imshow(Q);