function [bw1,bw2] = getbinaryrgb(im,ch1,ch2,chempty)

bw1 = logical(im(:,:,:,ch1) - im(:,:,:,chempty));
bw2 = logical(im(:,:,:,ch2) - im(:,:,:,chempty));
