TY - GEN

T1 - A total least-squares approach in two stages for semivariogram modeling of aeromagnetic data

AU - Felus, Yaron A.

AU - Schaffrin, Burkhard

PY - 2005

Y1 - 2005

N2 - Semivariogram analysis and estimation procedures are important for geostatisticians, Geographical Information Systems experts and earth scientists. A mathematical model, selected from a list of admissible functions, is fitted to the vector of empirical semivariogram values. Least-squares techniques (ordinary, weighted, and generalized) are often used in this process to compute the parameters within this mathematical model. However, least-squares fitting techniques consider the lag-distances (on the abscissa) of the empirical semivariogram graph as fixed or error-free; this is an imprecise assumption since these lag-distances are also measured quantities and include random errors. Here, a novel approach is presented that treats the empirical semivariogram lag-distances on the abscissa and the semivariogram values symmetrically under the assumption of random errors in both variables. This approach employs the relatively new Total Least-Squares (TLS) estimation technique, introduced in 1980 by Golub and van Loan, and provides a more accurate fitting with special emphasis on the crucial behavior near the origin. The main formulas and concepts of the TLS approach are described, along with its adaptation to the spatial context. Experiments of fitting semivariogram models to aeromagnetic data collected in Antarctica are presented, and various quality indicators are studied to compare the Weighted Least-Squares (WLS) with the Total Least-Squares methods. These studies demonstrate the general superiority of the TLS approach over the common WLS approach.

AB - Semivariogram analysis and estimation procedures are important for geostatisticians, Geographical Information Systems experts and earth scientists. A mathematical model, selected from a list of admissible functions, is fitted to the vector of empirical semivariogram values. Least-squares techniques (ordinary, weighted, and generalized) are often used in this process to compute the parameters within this mathematical model. However, least-squares fitting techniques consider the lag-distances (on the abscissa) of the empirical semivariogram graph as fixed or error-free; this is an imprecise assumption since these lag-distances are also measured quantities and include random errors. Here, a novel approach is presented that treats the empirical semivariogram lag-distances on the abscissa and the semivariogram values symmetrically under the assumption of random errors in both variables. This approach employs the relatively new Total Least-Squares (TLS) estimation technique, introduced in 1980 by Golub and van Loan, and provides a more accurate fitting with special emphasis on the crucial behavior near the origin. The main formulas and concepts of the TLS approach are described, along with its adaptation to the spatial context. Experiments of fitting semivariogram models to aeromagnetic data collected in Antarctica are presented, and various quality indicators are studied to compare the Weighted Least-Squares (WLS) with the Total Least-Squares methods. These studies demonstrate the general superiority of the TLS approach over the common WLS approach.

UR - http://www.scopus.com/inward/record.url?scp=43149088577&partnerID=8YFLogxK

M3 - ???researchoutput.researchoutputtypes.contributiontobookanthology.conference???

AN - SCOPUS:43149088577

SN - 0973422025

SN - 9780973422023

T3 - GIS and Spatial Analysis - 2005 Annual Conference of the International Association for Mathematical Geology, IAMG 2005

SP - 215

EP - 220

BT - GIS and Spatial Analysis - 2005 Annual Conference of the International Association for Mathematical Geology, IAMG 2005

T2 - 2005 Annual Conference of the International Association for Mathematical Geology: GIS and Spatial Analysis, IAMG 2005

Y2 - 21 August 2005 through 26 August 2005

ER -