High-order Taylor series expansion methods for error propagation in geographic information systems

Abstract

The quality of modeling results in GIS operations depends on how well we can track error propagating from inputs to outputs. Monte Carlo simulation, moment design and Taylor series expansion have been employed to study error propagation over the years. Among them, first-order Taylor series expansion is popular because error propagation can be analytically studied. Because most operations in GIS are nonlinear, first-order Taylor series expansion generally cannot meet practical needs, and higher-order approximation is thus necessary. In this paper, we employ Taylor series expansion methods of different orders to investigate error propagation when the random error vectors are normally and independently or dependently distributed. We also extend these methods to situations involving multi-dimensional output vectors. We employ these methods to examine length measurement of linear segments, perimeter of polygons and intersections of two line segments basic in GIS operations. Simulation experiments indicate that the fifth-order Taylor series expansion method is most accurate compared with the first-order and third-order method. Compared with the third-order expansion; however, it can only slightly improve the accuracy, but on the expense of substantially increasing the number of partial derivatives that need to be calculated. Striking a balance between accuracy and complexity, the third-order Taylor series expansion method appears to be a more appropriate choice for practical applications.