Robustness of noise filtering by kriging analysis

In geostatistics, factorial kriging is often proposed to filter noise. This filter is built from a linear model which is ideally suited to a Gaussian signal with additive independent noise. Robustness of the performance of factorial kriging is evaluated in less congenial situations. Three different types of noise are considered all perturbing a lognormally distributed signal. The first noise model is independent of the signal. The second noise model is heteroscedastic; its variance depends on the signal, yet noise and signal are uncorrelated. The third noise model is both heteroscedastic and linearly correlated with the signal. In ideal conditions, exhaustive sampling and additive independent noise, factorial kriging succeeds to reproduce the spatial patterns of high signal values. This score remains good in presence of heteroscedastic noise variance but falls quickly in presence of noise-to-signal correlation as soon as the sample becomes sparser.