For 4 binary images, you have totally 2^4 = 16 possible "colors" per pixel according to it being "water"/"soil" or a combination of these classes (if it is possible in your system).
Thus you can convert b1
...b4
to a single image:
>> l = b1 + 2*b2 + 4*b3 + 8*b4;
>> imshow( l ); colormap( rand(16,3) ); % random color map