for
each pixel [r,c]
if
(|I1[r,c] - I2[r,c]| >
threshold) then Iout[r,c] = 1 else Iout[r,c] = 0
Perform
connected components on Iout.
Remove small regions.
Perform
a closing with a small disk for merging close neighbors.
Compute
and return the bounding boxes B
of each remaining region.