function [px1,px2,py1,py2]=location(d)
%ͼƶλ
[y,x]=size(d);
myI=double(d);
a=1;
for i=1:x
    for j=1:y
        if myI(j,i,1)==1
            weiz1(1,a)=i;
            a=a+1;
            break;
        end
    end
end
px1=min(weiz1)+3;
px2=max(weiz1)-3;
a=1;
for i=1:y
    for j=1:x
        if myI(i,j,1)==1
            weiz2(1,a)=i;
            a=a+1;
            break;
        end
    end
end
py1=min(weiz2)+3;
py2=max(weiz2)-3;