A robust method of thin plate spline and its application to DEM construction

  • Authors:

  • Chuanfa Chen;Yanyan Li

  • Affiliations:

  • Geomatics College, Shandong University of Science and Technology, 266590 Qingdao, China and Key Laboratory of Surveying and Mapping Technology on Island and Reef, National Administration of Survey ...;Shanghai Astronomical Observatory, 200030 Shanghai, China

  • Venue:
  • Computers & Geosciences
  • Year:
  • 2012

Quantified Score

Hi-index 0.00

Visualization

Abstract

In order to avoid the ill-conditioning problem of thin plate spline (TPS), the orthogonal least squares (OLS) method was introduced, and a modified OLS (MOLS) was developed. The MOLS of TPS (TPS-M) can not only select significant points, termed knots, from large and dense sampling data sets, but also easily compute the weights of the knots in terms of back-substitution. For interpolating large sampling points, we developed a local TPS-M, where some neighbor sampling points around the point being estimated are selected for computation. Numerical tests indicate that irrespective of sampling noise level, the average performance of TPS-M can advantage with smoothing TPS. Under the same simulation accuracy, the computational time of TPS-M decreases with the increase of the number of sampling points. The smooth fitting results on lidar-derived noise data indicate that TPS-M has an obvious smoothing effect, which is on par with smoothing TPS. The example of constructing a series of large scale DEMs, located in Shandong province, China, was employed to comparatively analyze the estimation accuracies of the two versions of TPS and the classical interpolation methods including inverse distance weighting (IDW), ordinary kriging (OK) and universal kriging with the second-order drift function (UK). Results show that regardless of sampling interval and spatial resolution, TPS-M is more accurate than the classical interpolation methods, except for the smoothing TPS at the finest sampling interval of 20m, and the two versions of kriging at the spatial resolution of 15m. In conclusion, TPS-M, which avoids the ill-conditioning problem, is considered as a robust method for DEM construction.