for i=1:s1
    for j=1:s2
        if outimgindx(i,j)== 1
            outimg(i,j)= 0;
        elseif outimgindx(i,j)== 2
            outimg(i,j)= 0;
        elseif outimgindx(i,j)== 3
            outimg(i,j)= 0;
        elseif outimgindx(i,j)== 4
            outimg(i,j)= 1;
        end
    end
end
img(repmat(outimg,[1,1,3])==0)=0;