 
  
  
   
If we assume that the noise is white then its spatial distribution over the 
image will be random.
For the analysis of spatial data there are many measures of
randomness [18], often based on the assumption that the
observations follow a Poisson distribution.
This can be tested using the relative variance
  It is calculated by first counting the number of observations (in our case
the number of above threshold pixels in the difference map) in n windows,
 
It is calculated by first counting the number of observations (in our case
the number of above threshold pixels in the difference map) in n windows,
  , from which the mean,
 , from which the mean,   , and variance,
 , and variance,   , of
the
 , of
the   can be found.
Although the test is sensitive to the window size and point density it works
adequately as long as
  can be found.
Although the test is sensitive to the window size and point density it works
adequately as long as   is sufficiently large.
  is sufficiently large.
For our purposes we do not wish to detect the spatially random noise, but rather to avoid it in our thresholded image. We therefore select the threshold which maximises the relative variance, thereby maximising ``clumpiness'' (regions of change) and minimising the Poisson distribution (noise).