function b = bmaskc(bands,f)
% b = bmaskc(bands,f)
[m n] = size(f);
b = zeros(m,n);
[bm bn] = size(bands);
for i1 = 1:bm
  if bands(i1,15)~=0   % it's a rep
      b(bands(i1,15),bands(i1,16)) = 255;
    end
  end
end
