A Robust Multiquadratic Method and Its Application to DEM Construction

  • Geomatics College, Shandong University of Science and Technology, Qingdao 266590, China

Received date: 2013-11-08

  Revised date: 2013-12-01

  Online published: 2013-12-25


In order to resist the effect of outliers on DEM construction, a robust multiquadric method (MQ-R) has been developed. MQ-R firstly takes the estimation of the classical MQ as the initial values to compute the residuals of all sampling points, and then a weighted function has been constructed to determine the weights of sampling points based on the above residuals. Finally, a iteratively re-weighted MQ is formed to decrease the effect of outliers on DEM construction. At the same time, the smoothing parameter of MQ and MQ-R is determined based on a k-fold cross-validation. A synthetic surface was employed to comparatively analyze the estimation accuracies of MQ-R and the classcial MQ, where the sampling points are contaminated by three groups of errors with different distributions. These include the standard normal distribution, contaminated normal distribution with the contaminating proportion of 10%, 20% and 30%, and Cauchy distribution. Numerical tests indicate that when sampling errors are from the standard normal distribution, the accuracy of MQ-R is comparative to that of MQ. As the contaminating proportion increases, the accuracy of MQ becomes lower, whereas MQ-R can resist oultiers very well. When sampling errors are from Cauchy distribution, the results of MQ are completely destroyed, but those of MQ-R are still satisfactory. In conclusions, MQ-R with a high efficiency and a high robustness can be used to resist outliers in DEM construction.

Cite this article

CHEN Chuan-Fa, LI Wei, LI Meng-Fei, DAI Hong-Lei . A Robust Multiquadratic Method and Its Application to DEM Construction[J]. Journal of Geo-information Science, 2013 , 15(6) : 840 -845 . DOI: 10.3724/SP.J.1047.2013.00840


[1] Chen C F, Yue T X, Li Y Y. A high speed method of SMTS[J]. Computers & Geosciences, 2012(41):64-71.

[2] Chen C F, Fan Z M, Yue T X, et al. A robust estimator for the accuracy assessment of remote-sensing-derived DEMs[J]. International Journal of Remote Sensing, 2012, 33(8):2482-2497.

[3] Gong J, Li Z, Zhu Q, et al. Effects of various factors on the accuracy of DEMs: an intensive experimental investigation[J]. Photogrammetric Engineering and Remote Sensing, 2000, 66(9):1113-1117.

[4] Aguilar F J, Aguera F, Aguilar M A, et al. Effects of terrain morphology, sampling density, and interpolation methods on grid DEM accuracy[J]. Photogrammetric Engineering and Remote Sensing, 2005, 71(7):805-816.

[5] Gonga-Saholiariliva N, Gunnell Y, Petit C, et al. Techniques for quantifying the accuracy of gridded elevation models and for mapping uncertainty in digital terrain analysis[J]. Progress in Physical Geography, 2011, 35(6):739-764.

[6] Chen C F, Yue T X. A method of DEM construction and related error analysis[J]. Computers & Geosciences, 2010, 36(6):717-725.

[7] Hannah M J. Error detection and correction in digital terrain models[J]. Photogrammetric Engineering & Remote Sensing, 1981, 47(1):63-69.

[8] Felicísimo A M. Parametric statistical method for error detection in digital elevation models[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 1994, 49(4):29-33.

[9] López C. Locating some types of random errors in digital terrain models[J]. International Journal of Geographical Information Science, 1997, 11(7):677-698.

[10] Liu H, Jezek K C, O'Kelly M E. Detecting outliers in irregularly distributed spatial data sets by locally adaptive and robust statistical analysis and GIS[J]. International Journal of Geographical Information Science, 2001, 15(8):721-741.

[11] Maronna R A, Martin R D, Yohai V J. Robust statistics: Theory and methods[M]. New York: Wiley Blackwell, 2006.

[12] Chen C F, Li Y Y. A robust multiquadric interpolation for DEM construction[J]. Mathematical Geosciences, 2013, 45(3): 297-319.

[13] Franke R, Hagen H. Least squares surface approximation using multiquadrics and parametric domain distortion[J]. Computer Aided Geometric Design, 1999, 16(3):177-196.

[14] Chen C F, Li Y Y. A robust method of thin plate spline and its application to DEM construction[J]. Computers & Geosciences, 2012(48):9-16.

[15] Yue T X, Chen C F, Li B L. An adaptive method of high accuracy surface modeling and its application to simulating elevation surfaces[J]. Transactions in GIS, 2010, 14(5):615-630.

[16] 吕言.数字地面模型中多面函数内插法的研究[J].武汉测绘学院学报, 1981(2):14-28.

[17] 龚健雅.关于DTM中多面函数内插法几个问题的研究[J].测绘通报, 1985(5):12-16.

[18] Hoerl A E, Kennard R W. Ridge regression: Biased estimation for nonorthogonal problems[J]. Technometrics, 1970, 12(1):55-67.

[19] Hoerl A E, Kannard R W, Baldwin K F. Ridge regression:some simulations[J]. Communications in Statistics, 1975, 4(2):105-123.

[20] Falivene O, Cabrera L, Tolosana-Delgado R, et al. Interpolation algorithm ranking using cross-validation and the role of smoothing effect: A coal zone example[J]. Computers & Geosciences, 2010, 36(4):512-519.

[21] Rousseeuw P J, Croux C. Alternatives to the median absolute deviation[J]. Journal of the American Statistical Association, 1993, 88(424):1273-1283.

[22] Rousseeuw P J, Debruyne M, Engelen S, et al. Robustness and outlier detection in chemometrics[J]. Critical Reviews in Analytical Chemistry, 2006, 36(3-4):221-242.

[23] Rousseeuw P J. Least median of squares regression[J]. Journal of the American Statistical Association, 1984(79):871-880.

[24] Rousseeuw P J, Van Driessen K. Computing LTS regression for large data sets[J]. Data Mining and Knowledge Discovery, 2006, 12(1):29-45.