function C = setcolor(C,col,bw)

ind = find(bw);
C(ind) = col(1);
C(ind+size(C,1)*size(C,2)) = col(2);
C(ind+2*size(C,1)*size(C,2)) = col(3);