### Descripteur

Ajouter le résultat dans votre panier Affiner la recherche Interroger des sources externes

Etendre la recherche sur niveau(x) vers le bas

Progress towards a rigorous error propagation for total least-squares estimates / Burkhard Schaffrinin Journal of applied geodesy, vol 14 n° 2 (April 2020)

[article]inJournal of applied geodesy > vol 14 n° 2 (April 2020) . - pp 159 – 166

Titre : Progress towards a rigorous error propagation for total least-squares estimates Type de document : Article/Communication Auteurs : Burkhard Schaffrin, Auteur ; Kyle Snow, Auteur Année de publication : 2020 Article en page(s) : pp 159 – 166 Note générale : bibliographie Langues : Anglais ( eng)Descripteur : [Vedettes matières IGN] Géodésie spatiale

[Termes descripteurs IGN] équation non linéaire

[Termes descripteurs IGN] itération

[Termes descripteurs IGN] linéarisation

[Termes descripteurs IGN] matrice

[Termes descripteurs IGN] mesure de précision

[Termes descripteurs IGN] méthode des moindres carrés

[Termes descripteurs IGN] modèle d'erreur

[Termes descripteurs IGN] propagation d'erreurRésumé : (auteur) After several attempts at a formal derivation of the dispersion matrix for Total Least-Squares (TLS) estimates within an Errors-In-Variables (EIV) Model, here a refined approach is presented that makes rigorous use of the nonlinear normal equations, though assuming a Kronecker product structure for both observational dispersion matrices at this point. In this way, iterative linearization of a model (that can be established as being equivalent to the original EIV-Model) is avoided, which might be preferred since such techniques are based on the last iteration step only and, therefore, produce dispersion matrices for the estimated parameters that are generally too optimistic. Here, the error propagation is based on the (linearized total differential of the) exact nonlinear normal equations, which should lead to more trustworthy measures of precision. Numéro de notice : A2020-216 Affiliation des auteurs : non IGN Thématique : MATHEMATIQUE/POSITIONNEMENT Nature : Article nature-HAL : ArtAvecCL-RevueIntern DOI : 10.1515/jag-2019-0062 date de publication en ligne : 03/04/2020 En ligne : https://doi.org/10.1515/jag-2019-0062 Format de la ressource électronique : url article Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=94909 [article]Epipolar arrangement of satellite imagery by projection trajectory simplification / M. Wangin Photogrammetric record, vol 25 n° 132 (December 2010 - February 2011)

[article]inPhotogrammetric record > vol 25 n° 132 (December 2010 - February 2011) . - pp 422 - 436

Titre : Epipolar arrangement of satellite imagery by projection trajectory simplification Type de document : Article/Communication Auteurs : M. Wang, Auteur ; F. Hu, Auteur ; J. Li, Auteur Année de publication : 2010 Article en page(s) : pp 422 - 436 Note générale : Bibliographie Langues : Anglais ( eng)Descripteur : [Vedettes matières IGN] Acquisition d'image(s) et de donnée(s)

[Termes descripteurs IGN] barre de parallaxe

[Termes descripteurs IGN] courbe épipolaire

[Termes descripteurs IGN] image Ikonos

[Termes descripteurs IGN] image IRS

[Termes descripteurs IGN] image SPOT-HRG

[Termes descripteurs IGN] linéarisation

[Termes descripteurs IGN] lissage de courbe

[Termes descripteurs IGN] méthode des moindres carrés

[Termes descripteurs IGN] modèle stéréoscopique

[Termes descripteurs IGN] orbitographie

[Termes descripteurs IGN] trajectoireRésumé : (Auteur) This paper presents an efficient algorithm for approximate epipolar image generation from stereoscopic satellite images by projection trajectory simplification. In this algorithm, the original non-linear projection trajectory is simplified to a linear one by the least squares line fitting of sample points selected from the projection trajectory. The performance of the algorithm is evaluated using IRS-P5, SPOT5-HRG and IKONOS stereo-images acquired in different stereoscopic imaging modes over different terrain types. The results obtained indicate that the vertical parallax of each pair of conjugate image points on the epipolar images thus generated achieves accuracy at the sub-pixel level. Numéro de notice : A2010-534 Affiliation des auteurs : non IGN Thématique : IMAGERIE Nature : Article DOI : 10.1111/j.1477-9730.2010.00602.x date de publication en ligne : 22/12/2010 En ligne : https://doi.org/10.1111/j.1477-9730.2010.00602.x Format de la ressource électronique : URL article Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=30726 [article]## Réservation

Réserver ce document## Exemplaires (1)

Code-barres Cote Support Localisation Section Disponibilité 106-2010041 RAB Revue Centre de documentation En réserve 3L DisponibleThe curvilinear datum transformation / F.I. Okeke (1998)

Titre : The curvilinear datum transformation Type de document : Thèse/HDR Auteurs : F.I. Okeke, Auteur Editeur : Munich : Bayerische Akademie der Wissenschaften Année de publication : 1998 Collection : DGK - C Sous-collection : Dissertationen num. 481 Importance : 138 p. Format : 21 x 30 cm ISBN/ISSN/EAN : 978-3-7696-9521-2 Note générale : Bibliographie Langues : Anglais ( eng)Descripteur : [Vedettes matières IGN] Systèmes de référence et réseaux

[Termes descripteurs IGN] ellipsoïde de référence

[Termes descripteurs IGN] linéarisation

[Termes descripteurs IGN] méthode des moindres carrés

[Termes descripteurs IGN] modèle mathématique

[Termes descripteurs IGN] propagation d'erreur

[Termes descripteurs IGN] système de référence géodésique

[Termes descripteurs IGN] transformation de coordonnéesNuméro de notice : 67052 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Thèse étrangère Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=61694 ## Réservation

Réserver ce document## Exemplaires (2)

Code-barres Cote Support Localisation Section Disponibilité 67052-01 30.11 Livre Centre de documentation Géodésie Disponible67052-02 30.11 Livre Centre de documentation Géodésie DisponibleEllipsoidische und topographische Effekte im geodätischen Randwertproblem / Kurt Seitz (1997)

Titre : Ellipsoidische und topographische Effekte im geodätischen Randwertproblem Titre original : [Les effets ellipsoïdiques et topographiques dans le problème géodésique de valeurs limites] Type de document : Thèse/HDR Auteurs : Kurt Seitz, Auteur Editeur : Munich : Bayerische Akademie der Wissenschaften Année de publication : 1997 Collection : DGK - C Sous-collection : Dissertationen num. 483 Importance : 139 p. Format : 21 x 30 cm ISBN/ISSN/EAN : 978-3-7696-9523-6 Note générale : Bibliographie Langues : Allemand ( ger)Descripteur : [Vedettes matières IGN] Géodésie physique

[Termes descripteurs IGN] analyse harmonique

[Termes descripteurs IGN] approximation

[Termes descripteurs IGN] champ de pesanteur terrestre

[Termes descripteurs IGN] linéarisation

[Termes descripteurs IGN] potentiel de pesanteur terrestre

[Termes descripteurs IGN] problème des valeurs limitesRésumé : (Editeur) To determine both the external gravity potential W and the geometry of the earth's surface 5, various boundary value problems (bvp) can be formulated. They depend on the utilised observaWes L and whether the boundary is supposed to be known or unknown. If the geometry of 5 is already determined by the classical terrestrial techniques or by methods of satellite geodesy, then the fixed boundary value problem is under consideration. Otherwise the resulting bvp is of free type. The relation between the unknowns W, S and the observables L is given by boundary conditions. Generally, they are of non-linear structure.

This thesis focuses on the scalar free bvp and the fixed gravimetric bvp.

Analytical and numerical aspects require the linearisation of the boundary conditions. Therefore, suitable approximations have to be introduced for the unknowns. For the gravity potential, the normal gravity potential w is introduced. In the case of the scalar free boundary value problem, the boundary surface is approximated by the telluroid s, resulting from some telluroid mapping. In this context, the gravity potential W is substituted by the disturbing potential 6w := W w.

As shown in previous studies (Heck and Seitz, 1993, 1995; Seitz et al., 1994), the non-linear terms in the boundary condition can amount to maximum values of ±200 10~8ma~2 in case of the scalar free bvp, and ±40 10~8ms~2 for the fixed gravimetric problem. As a consequence of omitting the non-linear terms, the equipotential surfaces of the solution for the disturbing potential are shifted in vertical direction. Utilising a normal gravity field of Somigiiana-Pizzetti-type (e.g. GRS80), this shift takes values up to 4mm. If the gravitational portion of the normal gravity field is represented by a truncated series of spherical harmonics (e.g. Nv = 20), this effect can be reduced to 2mm. In the vicinity of the earth's physical surface, the accuracy of present day global solutions of the earth's gravity potential varies in the scope of dm-m. Therefore, the non-linear effects in the boundary condition can be neglected for the purpose of this thesis.

f

The linear boundary condition is now based on a linear operator D acting on the disturbing potential 6w. Applying the evaluation operator E,, the resulting linear functional is restricted to the telluroid. In case of the fixed bvp, the evaluation operator es is applied to D{Sw} and restricts the boundary condition to the physical earth's surface. The coefficients of the differential operator D are functional of the normal gravity potential w. Introducing an approximation wa of the normal gravity potential tu, the coefficients of D can be represented in terms of a Taylor series expansion.

For the scalar free boundary value problem as well as for the fixed gravimetric bvp the second order approximation of the differential operator D is derived. Thereby, the approximation wa of the normal gravity potential is represented by the potential of a Helmert's normal spheroid. This potential is symmetrical with respect to the mean earth's equatorial plane and independent of the geocentric longitude. In addition to the centrifugal part, this normal potential contains the isotropic term fj?/r and, furthermore, two gravitational terms, which are proportional to the zonal harmonic coefficients «/2 and «/4.

The analytical representation of the differential operator is carried out on the basis of wa. The extension of wa by the coefficient J± leads to the second order approximation of D. Further simplifications of wa result in the first order approximation and the isotropic approximation of D{6w} (fundamental equation of physical geodesy). This evaluation is done for a general normal potential and a Somigliana-Pizzetti reference field respectively.

The numerical studies clearly show that the error in the boundary condition, caused by isotropic approximation, exceeds the accuracy of relative gravimetric measurements by the multiple. It is also larger, by a factor of 5-10, than the non-linear terms in the corresponding boundary condition. According to this result, the ellipsoidal terms must not be neglected, as it is done in the isotropic approximation. If a Somigiiana-Pizzetti-field is used as normal field w, the first order approximation of the linear boundary condition produces a maximum error of ±0.5 10~8ms~2. The analytical

approximation of the differential operator finds its qualitative end, if the second order ellipsoidal terms, which are caused by the zonal coefficient J^, are taken into account. The errors due to the resulting second order approximation of the ellipsoidal terms in the linear boundary condition can be neglected in practice. The residuals are less than ±0.01 10~8ms~2.

The properties are different, when the normal field w contains tessera! and sectorial terms. In the present thesis, the normal gravity potential is given in terms of spherical harmonics, with a maximum degree Nv = 20. The apposition of the zonal coefficient J^ in the approximated normal field wa gives no better results for the approximated ellipsoidal terms. The maximum residuals are ±7 10~8m.s~2 for both levels of analytical approximation, the first order and the second order approximation. Even if second order tesseral and sectorial terms are included in w0, the accuracy of the approximation cannot be improved. According to the magnitude of the low degree spherical harmonic coefficients, all terms (e.g. up to degree and order 4) have to be included into the analytical approximation of the ellipsoidal terms in the boundary condition. The necessary analytical and numerical expenditure cannot be justified. Because of the slowly decreasing spectrum of this normal gravity potential, it must be concluded that the maximum degree Nv of the spherical harmonic representation has to be in accordance with the order of analytical approximation.

To solve the fixed or scalar free boundary value problem by the aid of a harmonic analysis, the respec--tive boundary condition must be related to a geometrically defined surface of revolution, symmetrical to the mean earth's rotation axis. In this context, a sphere K 9 k with radius r = a and the surface E 9 e of an ellipsoid of revolution is considered. The explained algorithm used for the harmonic ana-lysis results in the spherical harmonic spectrum, whether the (reduced) boundary data are given on K 3 k or E 9 e.

The analytical continuation of the boundary condition is realised by a formal Taylor series expansion of the evaluation operator E, or es respectively. Additionally, the optimal choice of the Taylor point has to be discussed. The Taylor point can be chosen at the auxiliary surface on which the data are analyticaly continued (sphere or surface of an ellipsoid). The Taylor point can also be situated at the original boundary surface. In case of the fixed bvp it can be situated at the physical earth's surface S 9 P, and on the telluroid s 3 p in case of the scalar free boundary value problem. The order of the Taylor series expansion of the evaluation operator is set up in accordance with the error level of the second order analytical approximation of the boundary condition. By this consideration the absolute error level of 10~8ms~2, which is aspired for the analytical approximation of the differential operator, is assigned to the analytically continued boundary condition. For each boundary condition we end up with a representation, where the analytical side contains an isotropic term and so-called ellipsoidal and topographical constituents. The isotropic term refers to the surface K or E. The ellipsoidal and topographical terms in the boundary condition are caused by the anisotropy of the normal potential and by the difference between the boundary surface and the auxiliary surface (K or E), on which the boundary condition is continued. They have to be evaluated either on the original boundary surface (S 3 P - fixed bvp; s 9 p - scalar free bvp) or on the auxiliary surface (K 9 k or E B e) according to the choice of the Taylor point. The ellipsoidal and topographical terms can be applied as a reduction to the original boundary data. The result is a spherical boundary value problem on the sphere K (Neumann-problem in the case of the fixed bvp), or rather a quasi-isotropic bvp with an ellipsoidal boundary surface. If the complete ellipsoidal and topographical terms are neglected, the constant radius approximation will be obtained. The reduction terms are functionals of the disturbing potential Sw we solve for. Therefore, an iterative procedure is required.

To guarantee an absolute accuracy of 1 lQ~8ms~2 for the continuation terms (ellipsoidal and topogra--phical terms), the partial derivatives of the disturbing potential are considered up to the 8th order. The partial derivatives of

j%- are considered up to the 4th order, with respect to the geocentric distance r. The numerical studies confirmed that this high degree of evaluation is necessary in the framework of the analytical continuation to the sphere. Otherwise, the level of accuracy obtained for the ellipsoidal terms cannot be reached for the continuation terms. If the analytical continuation is performed to the surface of an ellipsoid, the Taylor series can be truncated after the 5th order terms without anyloss of accuracy. The accuracy of approximation is not influenced by the choice of the Taylor point. The extremal values of the whole ellipsoidal and topographical part of the linear boundary condition, that has to be modelled by the analytical representation, decreases from ±100 10~5ms~2 in the case of a sphere to ±20 10~5ms~2, if the boundary operator is continued to the surface of an ellipsoid. These terms are identical with the neglected terms in case of the spherical approximation on K or the isotropic approximation on E, respectively. They have the same order of magnitude as the boundary data itself. Undoubtedly, this kind of approximation, where all ellipsoidal and topographical terms are neglected, leads to no high-precision solution of the fixed or scalar free boundary value problem.

First of all, the error of approximation depends on the order of the Taylor series, which are setup for the analytical continuation of the boundary condition. Therefore, the surface of an ellipsoid should be favoured above the sphere. The modelling of the normal field has a minor effect. With regard to the analytical handling, a normal field of SomigJiana-Pizzetti-type should be preferred to a truncated sphe--rical harmonic axpansion. The error in the developed analytical representation of the ellipsoidal and topographical terms reaches maximum values of ±1 10~8ms~2, if the GRS80 (Somigliana-Pizzetti-field) is used as a normal field. If the normal gravity potential is modelled by a spherical harmonic expansion, (OSU91alf) truncated at degree Nv = 20, the approximation errors slightly increase up to ±7 10-8ms-2.

When no ellipsoidal and topographical reductions are applied to the scalar gravity anomalies and when this data is analysed on a sphere, the resulting effect in the solution of the disturbing potential increases to ±40m2s~2. This corresponds to a vertical shift up to ±4m of the equipotential surfaces in the vicinity of the earth's surface. If the unreduced scalar gravity anomalies are analysed on the surface of an ellipsoid, the corresponding effect will be ±2m.

The proposed iterative procedure for solving the linear boundary value problem is investigated with respect to its convergence behaviour. When the analytical continuation to the surface of an ellipsoid is applied, a strong convergence can be observed. Already after 6 iteration steps, the spectrum of the solved disturbing potential does not change any more. The spectrum of the residuals represents the error of the solution in the frequency domain. It is the difference between the solved and the given spectrum, which was used to calculate the boundary data. The residuals illustrate the error in the space domain. The resulting vertical shift of the equipotential surfaces in the vicinity of the earth's physical surface reaches a maximum of ±4mm.

In case of the analytical continuation to a sphere, the iterative approach diverges. At first, the residual spectrum is decreasing within the scope of low and medium frequencies. Simultaneously, the error in the high frequency coefficients is increasing. The spectrum deteriorates in the whole frequency domain after about 5 iteration steps.

Only if the reduction terms (ellipsoidal and topographical terms) are generated in the first iteration step by the use of the given (exact) disturbing potential, convergence will be obtained. This was done in order to check the algorithms. For that reason it is supposed that the inaccurate initial solution is responsible for the divergent behaviour of the iterative process. It can be presumed that the analytical continuation to a sphere reacts more sensibly upon errors in the initial solution, than the analytical continuation to the surface of an ellipsoid. In this context, the problem of downward continuation is crucial, requiring further investigations.

(Editeur) To determine both the external gravity potential W and the geometry of the earth's surface 5, various boundary value problems (bvp) can be formulated. They depend on the utilised observaWes L and whether the boundary is supposed to be known or unknown. If the geometry of 5 is already determined by the classical terrestrial techniques or by methods of satellite geodesy, then the fixed boundary value problem is under consideration. Otherwise the resulting bvp is of free type. The relation between the unknowns W, S and the observables L is given by boundary conditions. Generally, they are of non-linear structure.

This thesis focuses on the scalar free bvp and the fixed gravimetric bvp.

Analytical and numerical aspects require the linearisation of the boundary conditions. Therefore, suitable approximations have to be introduced for the unknowns. For the gravity potential, the normal gravity potential w is introduced. In the case of the scalar free boundary value problem, the boundary surface is approximated by the telluroid s, resulting from some telluroid mapping. In this context, the gravity potential W is substituted by the disturbing potential 6w := W ? w.

As shown in previous studies (Heck and Seitz, 1993, 1995; Seitz et al., 1994), the non-linear terms in the boundary condition can amount to maximum values of ±200 10~8ma~2 in case of the scalar free bvp, and ±40 10~8ms~2 for the fixed gravimetric problem. As a consequence of omitting the non-linear terms, the equipotential surfaces of the solution for the disturbing potential are shifted in vertical direction. Utilising a normal gravity field of Somigiiana-Pizzetti-type (e.g. GRS80), this shift takes values up to 4mm. If the gravitational portion of the normal gravity field is represented by a truncated series of spherical harmonics (e.g. Nv = 20), this effect can be reduced to 2mm. In the vicinity of the earth's physical surface, the accuracy of present day global solutions of the earth's gravity potential varies in the scope of dm-m. Therefore, the non-linear effects in the boundary condition can be neglected for the purpose of this thesis.

f

The linear boundary condition is now based on a linear operator D acting on the disturbing potential 6w. Applying the evaluation operator E,, the resulting linear functional is restricted to the telluroid. In case of the fixed bvp, the evaluation operator es is applied to D{Sw} and restricts the boundary condition to the physical earth's surface. The coefficients of the differential operator D are functional of the normal gravity potential w. Introducing an approximation wa of the normal gravity potential tu, the coefficients of D can be represented in terms of a Taylor series expansion.

For the scalar free boundary value problem as well as for the fixed gravimetric bvp the second order approximation of the differential operator D is derived. Thereby, the approximation wa of the normal gravity potential is represented by the potential of a Helmert's normal spheroid. This potential is symmetrical with respect to the mean earth's equatorial plane and independent of the geocentric longitude. In addition to the centrifugal part, this normal potential contains the isotropic term fj?/r and, furthermore, two gravitational terms, which are proportional to the zonal harmonic coefficients «/2 and «/4.

The analytical representation of the differential operator is carried out on the basis of wa. The extension of wa by the coefficient J± leads to the second order approximation of D. Further simplifications of wa result in the first order approximation and the isotropic approximation of D{6w} (fundamental equation of physical geodesy). This evaluation is done for a general normal potential and a Somigliana-Pizzetti reference field respectively.

The numerical studies clearly show that the error in the boundary condition, caused by isotropic approximation, exceeds the accuracy of relative gravimetric measurements by the multiple. It is also larger, by a factor of 5-10, than the non-linear terms in the corresponding boundary condition. According to this result, the ellipsoidal terms must not be neglected, as it is done in the isotropic approximation. If a Somigiiana-Pizzetti-field is used as normal field w, the first order approximation of the linear boundary condition produces a maximum error of ±0.5 10~8ms~2. The analytical

approximation of the differential operator finds its qualitative end, if the second order ellipsoidal terms, which are caused by the zonal coefficient J^, are taken into account. The errors due to the resulting second order approximation of the ellipsoidal terms in the linear boundary condition can be neglected in practice. The residuals are less than ±0.01 10~8ms~2.

The properties are different, when the normal field w contains tessera! and sectorial terms. In the present thesis, the normal gravity potential is given in terms of spherical harmonics, with a maximum degree Nv = 20. The apposition of the zonal coefficient J^ in the approximated normal field wa gives no better results for the approximated ellipsoidal terms. The maximum residuals are ±7 10~8m.s~2 for both levels of analytical approximation, the first order and the second order approximation. Even if second order tesseral and sectorial terms are included in w0, the accuracy of the approximation cannot be improved. According to the magnitude of the low degree spherical harmonic coefficients, all terms (e.g. up to degree and order 4) have to be included into the analytical approximation of the ellipsoidal terms in the boundary condition. The necessary analytical and numerical expenditure cannot be justified. Because of the slowly decreasing spectrum of this normal gravity potential, it must be concluded that the maximum degree Nv of the spherical harmonic representation has to be in accordance with the order of analytical approximation.

To solve the fixed or scalar free boundary value problem by the aid of a harmonic analysis, the respec--tive boundary condition must be related to a geometrically defined surface of revolution, symmetrical to the mean earth's rotation axis. In this context, a sphere K 9 k with radius r = a and the surface E 9 e of an ellipsoid of revolution is considered. The explained algorithm used for the harmonic ana-lysis results in the spherical harmonic spectrum, whether the (reduced) boundary data are given on K 3 k or E 9 e.

The analytical continuation of the boundary condition is realised by a formal Taylor series expansion of the evaluation operator E, or es respectively. Additionally, the optimal choice of the Taylor point has to be discussed. The Taylor point can be chosen at the auxiliary surface on which the data are analyticaly continued (sphere or surface of an ellipsoid). The Taylor point can also be situated at the original boundary surface. In case of the fixed bvp it can be situated at the physical earth's surface S 9 P, and on the telluroid s 3 p in case of the scalar free boundary value problem. The order of the Taylor series expansion of the evaluation operator is set up in accordance with the error level of the second order analytical approximation of the boundary condition. By this consideration the absolute error level of 10~8ms~2, which is aspired for the analytical approximation of the differential operator, is assigned to the analytically continued boundary condition. For each boundary condition we end up with a representation, where the analytical side contains an isotropic term and so-called ellipsoidal and topographical constituents. The isotropic term refers to the surface K or E. The ellipsoidal and topographical terms in the boundary condition are caused by the anisotropy of the normal potential and by the difference between the boundary surface and the auxiliary surface (K or E), on which the boundary condition is continued. They have to be evaluated either on the original boundary surface (S 3 P - fixed bvp; s 9 p - scalar free bvp) or on the auxiliary surface (K 9 k or E B e) according to the choice of the Taylor point. The ellipsoidal and topographical terms can be applied as a reduction to the original boundary data. The result is a spherical boundary value problem on the sphere K (Neumann-problem in the case of the fixed bvp), or rather a quasi-isotropic bvp with an ellipsoidal boundary surface. If the complete ellipsoidal and topographical terms are neglected, the constant radius approximation will be obtained. The reduction terms are functionals of the disturbing potential Sw we solve for. Therefore, an iterative procedure is required.

To guarantee an absolute accuracy of 1 lQ~8ms~2 for the continuation terms (ellipsoidal and topogra--phical terms), the partial derivatives of the disturbing potential are considered up to the 8th order. The partial derivatives of

j%- are considered up to the 4th order, with respect to the geocentric distance r. The numerical studies confirmed that this high degree of evaluation is necessary in the framework of the analytical continuation to the sphere. Otherwise, the level of accuracy obtained for the ellipsoidal terms cannot be reached for the continuation terms. If the analytical continuation is performed to the surface of an ellipsoid, the Taylor series can be truncated after the 5th order terms without anyloss of accuracy. The accuracy of approximation is not influenced by the choice of the Taylor point. The extremal values of the whole ellipsoidal and topographical part of the linear boundary condition, that has to be modelled by the analytical representation, decreases from ±100 10~5ms~2 in the case of a sphere to ±20 10~5ms~2, if the boundary operator is continued to the surface of an ellipsoid. These terms are identical with the neglected terms in case of the spherical approximation on K or the isotropic approximation on E, respectively. They have the same order of magnitude as the boundary data itself. Undoubtedly, this kind of approximation, where all ellipsoidal and topographical terms are neglected, leads to no high-precision solution of the fixed or scalar free boundary value problem.

First of all, the error of approximation depends on the order of the Taylor series, which are setup for the analytical continuation of the boundary condition. Therefore, the surface of an ellipsoid should be favoured above the sphere. The modelling of the normal field has a minor effect. With regard to the analytical handling, a normal field of SomigJiana-Pizzetti-type should be preferred to a truncated sphe--rical harmonic axpansion. The error in the developed analytical representation of the ellipsoidal and topographical terms reaches maximum values of ±1 10~8ms~2, if the GRS80 (Somigliana-Pizzetti-field) is used as a normal field. If the normal gravity potential is modelled by a spherical harmonic expansion, (OSU91alf) truncated at degree Nv = 20, the approximation errors slightly increase up to ±7 10-8ms-2.

When no ellipsoidal and topographical reductions are applied to the scalar gravity anomalies and when this data is analysed on a sphere, the resulting effect in the solution of the disturbing potential increases to ±40m2s~2. This corresponds to a vertical shift up to ±4m of the equipotential surfaces in the vicinity of the earth's surface. If the unreduced scalar gravity anomalies are analysed on the surface of an ellipsoid, the corresponding effect will be ±2m.

The proposed iterative procedure for solving the linear boundary value problem is investigated with respect to its convergence behaviour. When the analytical continuation to the surface of an ellipsoid is applied, a strong convergence can be observed. Already after 6 iteration steps, the spectrum of the solved disturbing potential does not change any more. The spectrum of the residuals represents the error of the solution in the frequency domain. It is the difference between the solved and the given spectrum, which was used to calculate the boundary data. The residuals illustrate the error in the space domain. The resulting vertical shift of the equipotential surfaces in the vicinity of the earth's physical surface reaches a maximum of ±4mm.

In case of the analytical continuation to a sphere, the iterative approach diverges. At first, the residual spectrum is decreasing within the scope of low and medium frequencies. Simultaneously, the error in the high frequency coefficients is increasing. The spectrum deteriorates in the whole frequency domain after about 5 iteration steps.

Only if the reduction terms (ellipsoidal and topographical terms) are generated in the first iteration step by the use of the given (exact) disturbing potential, convergence will be obtained. This was done in order to check the algorithms. For that reason it is supposed that the inaccurate initial solution is responsible for the divergent behaviour of the iterative process. It can be presumed that the analytical continuation to a sphere reacts more sensibly upon errors in the initial solution, than the analytical continuation to the surface of an ellipsoid. In this context, the problem of downward continuation is crucial, requiring further investigations.Numéro de notice : 28009 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Thèse étrangère Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=63356 ## Réservation

Réserver ce document## Exemplaires (2)

Code-barres Cote Support Localisation Section Disponibilité 28009-01 30.40 Livre Centre de documentation Géodésie Disponible28009-02 30.40 Livre Centre de documentation Géodésie DisponibleBestimmung von Modellparametern der Erde durch Analyse ihrer Drehbewegung / H. Fröhlich (1994)

Titre : Bestimmung von Modellparametern der Erde durch Analyse ihrer Drehbewegung Titre original : [Estimation des paramètres de modèle de la Terre au travers de l'analyse de sa rotation] Type de document : Thèse/HDR Auteurs : H. Fröhlich, Auteur Editeur : Munich : Bayerische Akademie der Wissenschaften Année de publication : 1994 Collection : DGK - C Sous-collection : Dissertationen num. 420 Importance : 102 p. Format : 21 x 30 cm Note générale : Bibliographie Langues : Allemand ( ger)Descripteur : [Vedettes matières IGN] Géodésie physique

[Termes descripteurs IGN] correction différentielle

[Termes descripteurs IGN] linéarisation

[Termes descripteurs IGN] modèle mathématique

[Termes descripteurs IGN] problème des valeurs limites

[Termes descripteurs IGN] rotation de la Terre

[Termes descripteurs IGN] théorème de LiouvilleRésumé : (Auteur) The rotation of the Earth is modeled by means of the linearized EULER-LIOUVlLLE equation. The different components described physically or by indices are combined in a joint model of the excitation function. Pa-rameters are estimated by differential correction. Parts of the excitation function not explicitly considered in the model can be derived from the residuals. The linearized EULER-LIOUVlLLE equation is solved as a bound-ary value problem; the formulas according to the method of infinitely many variables are derived. Results of the estimation of parameters are discussed by using their accuracies and amount and temporal course of the residuals. Numéro de notice : 28067 Affiliation des auteurs : non IGN Thématique : POSITIONNEMENT Nature : Thèse étrangère Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=63414 ## Réservation

Réserver ce document## Exemplaires (1)

Code-barres Cote Support Localisation Section Disponibilité 2806701 30.40 Livre Centre de documentation Géodésie DisponibleLösung des fixen geodätischen Randwertproblems mit Hilfe der Randelementmethode / R. Klees (1992)PermalinkRelativistic motion of gyroscopes and space gradiometry / E. Gill (1992)PermalinkTowards a local split window method over land surface / François Beckerin International Journal of Remote Sensing IJRS, vol 11 n° 3 (March 1990)PermalinkSpot et la photogrammétrie analytique : solution de quelques problèmes - partie 1 : modélisation de l'imagerie Spot et spatiotriangulation / Philippe Hottierin Bulletin [Société Française de Photogrammétrie et Télédétection], n° 117 (Janvier 1990)PermalinkFestschrift Friedrich Ackermann zum 60. Geburtstag / Institut für Photogrammetrie der Universität Stuttgart (Stuttgart, Allemagne) (1989)PermalinkMéthode des moindres carrés / Henri Duquenne (1986)PermalinkEin dreidimensionales geometrisch-physikalisches Modell für konventionelle Beobachtungen / J. Zaiser (1984)PermalinkEine approximative Lösung der allgemeinen linearen Geodätischen Randwertaufgabe durch Reihenentwicklungen nach Kugelfunktionen / E.H. Knickmeyer (1984)PermalinkAdvanced physical geodesy / Helmut Moritz (1980)PermalinkThéorie des erreurs à deux composantes : Bruit et silence, 2. Tome 2, Systèmes de mesure pluridimensionnels à configuration margée associée à une "figure" ou la compensation des mesures infiniment précises / Philippe Hottier (1980)Permalink