next up previous
Next: Modelling the spatial distribution Up: Two New Thresholding Methods Previous: Two New Thresholding Methods

Modelling the spatial distribution of the noise

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 tex2html_wrap_inline390 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, tex2html_wrap_inline394 , from which the mean, tex2html_wrap_inline396 , and variance, tex2html_wrap_inline398 , of the tex2html_wrap_inline400 can be found. Although the test is sensitive to the window size and point density it works adequately as long as tex2html_wrap_inline396 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).



Paul L Rosin
Tue Aug 25 16:29:46 BST 1998