% Using each image and its mask,
% combine all the images. In areas
% where the mask indicates a difference
% do not add the part of the image
% to the average (ignore it).

mclean = double(m1);
bclean = double(b1n) .* double(m1);

mclean = mclean + double(m2);
bclean = bclean + double(b2n) .* double(m2);

mclean = mclean + double(m3);
bclean = bclean + double(b3n) .* double(m3);

mclean = mclean + double(m4);
bclean = bclean + double(b4n) .* double(m4);

mclean = mclean + double(m5);
bclean = bclean + double(b5n) .* double(m5);

mclean = mclean + double(m6);
bclean = bclean + double(b6n) .* double(m6);

mclean = mclean + double(m7);
bclean = bclean + double(b7n) .* double(m7);

mclean = mclean + double(m8);
bclean = bclean + double(b8n) .* double(m8);

mclean = mclean + double(m9);
bclean = bclean + double(b9n) .* double(m9);

mclean = mclean + double(ma);
bclean = bclean + double(ban) .* double(ma);

bclean = bclean ./ mclean;

image(abs(bclean - bt)>2);