TR32-Database: Database of Transregio 32

[977] - Approximate continuation of harmonic functions in geodesy: A spline based least squares approach with regularization

All available metadata of the dataset are listed below. Some features are available, e.g. download of dataset or additional description file.

Jager, G., Kunoth, A., Schuh, W., 2013. Approximate continuation of harmonic functions in geodesy: A spline based least squares approach with regularization. Journal of Computational and Applied Mathematics, 237 (1), 62 - 82. DOI: 10.1016/
Title(s):Main Title: Approximate continuation of harmonic functions in geodesy: A spline based least squares approach with regularization
Description(s):Abstract: This paper is concerned with the approximate reconstruction of the earth’s potential field from geometric and gravimetric data. This is an ill-posed problem involving typically large amounts of data which are to be continued by a harmonic function. The standard approach in geodesy is based on spherical harmonics which are globally supported. Thus, a least squares approach for the data fitting yields a linear system of equations with a fully populated system matrix. This becomes computationally prohibitive for large amounts of data and, therefore, presents the biggest bottleneck for fast and efficient computations. Motivated by the early work [24], we propose in this paper an alternative and pose the harmonicity requirement on the continuation together with the data fitting as a minimization problem for a least squares functional with regularization involving the Laplacian. This approach enables the use of locally supported functions in the reconstruction for which we employ tensor products of cubic splines. The linear system resulting from the weighted least squares approach is therefore sparsely populated which allows for iterative solvers of complexity proportional to the total number of unknowns. We extensively study the choice of the regularization parameter balancing the data fit and the harmonicity requirement for both synthetic as well as earth potential data. We compare the results with discretizations using finite differences and finite elements for solving Laplace’s equation.
Identifier(s):DOI: 10.1016/
Responsible Party
Creator(s):Author: Gabriela Jager
Author: Angela Kunoth
Author: Wolf-Dieter Schuh
TR32 Topic:Other
Related Sub-project(s):C1, C7
Subject(s):CRC/TR32 Keywords: Regularization, Least Squares
GEMET: mathematical analysis
File Details
File Name:Jager_2013_JCAM.pdf
Data Type:Text
File Size:3185 kB (3.11 MB)
Date(s):Date Submitted: 2011-03-18
Available: 2012-07-20
Mime Type:application/pdf
Data Format:PDF
Download Permission:OnlyTR32
Download Information:Copyright © 2012 Elsevier B.V. All rights reserved.
General Access and Use Conditions:For internal use in TR32 only.
Access Limitations:For internal use in TR32 only.
Licence:TR32DB Data policy agreement
North:-no map data
Measurement Region:Other
Measurement Location:--Other--
Specific Informations - Publication
Article Type:Journal
Source:Journal of Computational and Applied Mathematics
Page Range:62 - 82
Metadata Details
Metadata Creator:Angela Kunoth
Metadata Created:2014-07-16
Metadata Last Updated:2014-07-16
Funding Phase:2
Metadata Language:English
Metadata Version:V40
Dataset Metrics
Page Visits:138
Metadata Downloads:0
Dataset Downloads:1
Dataset Activity