rnd = mmbinary(rand(shape(C)[0],shape(C)[1]) > 0.996)
Ba=mmintersec(mmgradm(C), rnd)
Bg=mmdil(Ba,mmimg2se(mmgradm(mmpad(mmseshow(B),mmsebox(2)))))
Bg=mmunion(mmdil(Ba,mmsecross(2)),Bg)
AC=mmsymdif(mmunion(mmgray(A,'uint8',1),mmgray(C,'uint8',2)),mmgray(Bg,'uint8',1))
mmshow(AC)