imshow('D:\MATLAB7\work\z7.bmp');
bw1=imread('D:\MATLAB7\work\z7.bmp');
bw2=imread('D:\MATLAB7\work\a9.bmp');
[height,width]=size(bw1);
[theight,twidth]=size(bw2);
wcyz=17;			%誤差閾值,9、8、2的閾值取60；7、4、3、0取40；1取30；
count=0;
w=zeros(1,30);
h=zeros(1,30);
for i=1:2:height-theight	%隔行隔列掃描匹配
    for j=1:2:width-twidth
	nerror=0;
	for m=1:2:theight
	    for n=1:2:twidth
		if(bw1(i+m-1,j+n-1)~=bw2(m,n))
		   nerror=nerror+1;
		end
	    end 
	       if(nerror>wcyz),break;
	       end
	end
		if(nerror<=wcyz)
		  count=count+1;
		  w(count)=j;
		  h(count)=i;
		end
    end
end
		bw3=imread('D:\MATLAB7\work\z7.bmp');
	    for j1=1:count
		for i1=0:twidth,bw3(h(j1),w(j1)+i1)=0;
			bw3(h(j1)+theight,w(j1)+i1)=0;

		end
	    for i1=0:theight,bw3(h(j1)+i1,w(j1))=0;
			bw3(h(j1)+i1,w(j1)+twidth)=0;
    end
end
    figure,imshow(bw3);imwrite(bw3,'D:\MATLAB7\work\lw2.bmp','bmp');