Steepest Descent Method for Representing Spatially Correlated Uncertainty in GIS

Abstract

All spatial data in a geographic information system (GIS) intrinsically contain uncertainty. Simulations could be used for many GIS applications in order to estimate confidence ranges of certain analyses and project worst-case scenarios. For those applications, generation of Gaussian random fields is essential to simulate the uncertainty effects, because errors in spatial data are assumed dependent upon the Gaussian distribution. Gaussian fields with no spatial dependency could be assumed because of their simple concept and easy computation, but the reality is that spatial errors have a spatially correlated nature. For this reason, the intensive matrix computation for generating spatially autocorrelated Gaussian random fields requires the solution of a large, sparse linear system: