Descripteur
Documents disponibles dans cette catégorie (86)
Ajouter le résultat dans votre panier
Visionner les documents numériques
Affiner la recherche Interroger des sources externes
Etendre la recherche sur niveau(x) vers le bas
Une exploration des signaux en ondelettes / Stéphane Mallat (2000)
Titre : Une exploration des signaux en ondelettes Type de document : Monographie Auteurs : Stéphane Mallat, Auteur Editeur : Palaiseau : Ecole Polytechnique EP Année de publication : 2000 Importance : 637 p. Format : 17 x 24 cm ISBN/ISSN/EAN : 978-2-7302-0733-1 Note générale : Bibliographie Langues : Français (fre) Descripteur : [Vedettes matières IGN] Traitement du signal
[Termes IGN] approximation
[Termes IGN] ondelette
[Termes IGN] principe d'incertitude de Heisenberg
[Termes IGN] traitement du signal
[Termes IGN] transformation de FourierIndex. décimale : 24.20 Traitement du signal Résumé : (Auteur) Traiter un signal, de la musique à l'image, est avant tout affaire de représentation. Des sinusoïdes aux ondelettes, il s'agit de trouver les structures élémentaires qui permettent de révéler le contenu utile d'un signal. En partant de la transformée de Fourier, l'ouvrage montre que la construction de représentations localisées en temps et fréquence est un jeu de pavage organisé autour du principe d'incertitude. Les transformées et bases orthogonales d'ondelettes, de paquets d'ondelettes et de cosinus, trouvent leurs applications au travers de la théorie de l'approximation et des statistiques, en passant par le chemin des algorithmes rapides. Compression audio ou vidéo, débruitage, analyse de singularités, de multifractales ou de transitoires, sont autant d'applications qui nous conduisent à la frontière des mathématiques comprises à ce jour. Des logiciels disponibles sur Internet permettent au lecteur d'appliquer les algorithmes et les théorèmes du livre, et de développer son intuition sur des exemples. Numéro de notice : 18974 Affiliation des auteurs : non IGN Thématique : IMAGERIE/POSITIONNEMENT Nature : Monographie Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=55523 Réservation
Réserver ce documentExemplaires(1)
Code-barres Cote Support Localisation Section Disponibilité 18974-01 24.20 Livre Centre de documentation Physique Disponible Approximation hilbertienne / M. Atteia (1999)
Titre : Approximation hilbertienne : splines, ondelettes, fractales Type de document : Monographie Auteurs : M. Atteia, Auteur ; J. Gaches, Auteur Editeur : Les Ulis : EDP Sciences Année de publication : 1999 Collection : Grenoble Sciences Importance : 158 p. Format : 21 x 30 cm ISBN/ISSN/EAN : 978-2-7061-0860-0 Note générale : Bibliographie Langues : Français (fre) Descripteur : [Vedettes matières IGN] Analyse numérique
[Termes IGN] approximation
[Termes IGN] B-Spline
[Termes IGN] espace de Hilbert
[Termes IGN] fonction spline
[Termes IGN] objet fractal
[Termes IGN] ondelette
[Termes IGN] Ploynôme
[Termes IGN] transformation de FourierRésumé : (Editeur) Splines, ondelettes et fractales sont des outils mathématiques très utilisés dans la pratique. L'objectif de l'ouvrage est d'expliquer l'origine de ces outils. Les premiers chapitres rappellent les connaissances de base telles les espaces de Hilbert et les séries de Fourier. Puis, utilisant les espaces classiques de Beppo Lévi et de Sobolev, on introduit la notion de noyau hilbertien. Ces noyaux hilbertiens permettent enfin d'élaborer les théories des splines, ondelettes et fractales. L'ouvrage n'est pas un simple traité : chaque chapitre est suivi d'exercices qui constituent un complément pour la compréhension des concepts et l'utilisation des outils. Numéro de notice : 18860 Affiliation des auteurs : non IGN Thématique : MATHEMATIQUE Nature : Monographie Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=55495 Réservation
Réserver ce documentExemplaires(1)
Code-barres Cote Support Localisation Section Disponibilité 18860-01 23.40 Livre Centre de documentation Mathématiques Disponible Segmentation d'images aériennes texturées par champs de Markov / Guillaume Verza (1999)
Titre : Segmentation d'images aériennes texturées par champs de Markov Type de document : Mémoire Auteurs : Guillaume Verza, Auteur Editeur : Champs/Marne : Université de Marne-la-Vallée Année de publication : 1999 Importance : 73 p. Format : 21 x 30 cm Note générale : Bibliographie
Mémoire de fin de stage, DEA Sciences de l'information géographiqueLangues : Français (fre) Descripteur : [Vedettes matières IGN] Traitement d'image
[Termes IGN] approximation
[Termes IGN] champ aléatoire de Markov
[Termes IGN] classification dirigée
[Termes IGN] classification non dirigée
[Termes IGN] image aérienne
[Termes IGN] ondelette
[Termes IGN] processus stochastique
[Termes IGN] résolution multiple
[Termes IGN] segmentation d'image
[Termes IGN] texture d'imageIndex. décimale : DSIG Mémoires du master 2 IG, du master 2 SIG, de l'ex DEA SIG Résumé : (Auteur) Ce rapport traite de la segmentation d'images aériennes, en noir et blanc, texturées. Segmenter une image consiste à créer une partition de l'image en régions. Pour cela, il est nécessaire de modéliser les textures afin de pouvoir calculer des paramètres appropriés. Une classification est ensuite réalisée sur ces paramètres pour déterminer les différentes régions de l'image. Les champs de Markov sont en effet bien adaptés au traitement d'images car ils permettent de prendre en compte les interactions entre pixels et sont à l'origine de nombreuses techniques de relaxation performantes. Deux modèles de segmentation sont proposés. Le premier est non supervisé. Il modèlise les textures comme des champs markoviens gaussiens anisotropes caractérisés par des variances conditionnelles directionnelles facilement estimables. La classification se fait ensuite, suivant chaque direction, par relaxation aux différents niveaux d'une approximation multirésolution de l'image. L'initialisation de la classification est réalisée au niveau le plus haut par un algorithme de "Fuzzy-C-means" avec un terme d'entropie. Les classifications obtenues pour chaque direction sont fusionnées en fin de traitement. Le deuxième modèle de segmentation est supervisé. Les textures sont modélisées par les interactions entre pixels pour un voisinage choisi par l'opérateur. La classification se fait par relaxation à partir des classes fixées par l'opérateur, en tenant compte de tous les paramètres à la fois. Ces deux modèles permettent d'aboutir à des résultats à peu près comparables. Ils segmentent correctement les zones homogènes et les micro-textures. Par contre ils sont tous les deux inefficaces pour segmenter les macro-textures. Des modèles plus performants pourraient être obtenus en rajoutant aux niveaux de gris des informations de ... Numéro de notice : 11051 Affiliation des auteurs : non IGN Thématique : IMAGERIE Nature : Mémoire Master 2 IG Organisme de stage : INRIA Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=49476 Réservation
Réserver ce documentExemplaires(2)
Code-barres Cote Support Localisation Section Disponibilité 11051-01 DSIG Livre Centre de documentation Travaux d'élèves Disponible 11051-02 K314 Livre LASTIG Dépôt en unité Exclu du prêt Approximation des profondeurs d'une scène urbaine à partir d'images à grande échelle / A. Benzenati (1997)
Titre : Approximation des profondeurs d'une scène urbaine à partir d'images à grande échelle Type de document : Mémoire Auteurs : A. Benzenati, Auteur Editeur : Paris : Université de Paris 6 Pierre et Marie Curie Année de publication : 1997 Importance : 54 p. Format : 21 x 30 cm Note générale : Bibliographie
DESS de mathématiques appliquéesLangues : Français (fre) Descripteur : [Vedettes matières IGN] Traitement d'image
[Termes IGN] appariement d'images
[Termes IGN] approximation
[Termes IGN] détection du bâti
[Termes IGN] milieu urbain
[Termes IGN] photographie à grande échelle
[Termes IGN] scène
[Termes IGN] stéréoscopieIndex. décimale : DESS Divers DESS et masters actuels professionnels ou spécialisés Résumé : (Auteur) Cette étude propose une méthode pour obtenir une approximation de profondeurs d'une scène prise en milieu urbain, à partir d'images aériennes à grande échelle. Un algorithme de correspondance, appelé corrélateur, permet de restituer les profondeurs à partir d'images à petite ou moyenne échelle avec des résultats décevants pour les images à grande échelle. L'approximation par la méthode des moindres carrés et la mise en place d'un algorithme de sélection de régions obtenues à partir d'une image des contours, permettent d'obtenir une description satisfaisante de la scène étudiée. Numéro de notice : 67716 Affiliation des auteurs : non IGN Thématique : IMAGERIE Nature : Mémoire DESS divers Organisme de stage : MATIS (IGN) Permalink : https://documentation.ensg.eu/index.php?lvl=notice_display&id=53624 Réservation
Réserver ce documentExemplaires(1)
Code-barres Cote Support Localisation Section Disponibilité 67716-01 DESS Livre Centre de documentation En réserve Mezzanine Disponible Ellipsoidische 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 IGN] analyse harmonique
[Termes IGN] approximation
[Termes IGN] champ de pesanteur terrestre
[Termes IGN] linéarisation
[Termes IGN] potentiel de pesanteur terrestre
[Termes IGN] problème des valeurs limitesIndex. décimale : 30.40 Géodésie physique Ré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 documentExemplaires(2)
Code-barres Cote Support Localisation Section Disponibilité 28009-01 30.40 Livre Centre de documentation Géodésie Disponible 28009-02 30.40 Livre Centre de documentation Géodésie Disponible Geodetic theory today / Fernando Sanso (1995)PermalinkInfluence de la rugosité de surface en radiométrie micro-onde des sols nus / L. Laguerre (1995)PermalinkAnalyse der Form und Lage von Objekten im Raum / E. Drixler (1993)PermalinkDifferentielle Verschiebungen und Drehstreckungen in dreidimensionalen Koordinatensystemen / D. Ehlert (1991)PermalinkRaumbezogene Informationssysteme und digitale Geländemodelle / D. Frisch (1991)PermalinkDigital filters / R.W. Hamming (1989)PermalinkAutomation in photogrammetric block adjustment systems / Tapani Sarjakoski (1988)PermalinkA scattering model for perfectly conducting random surfaces - 1 : model development / A.K. Fung in International Journal of Remote Sensing IJRS, vol 8 n° 11 (November 1987)PermalinkAnalyse und Synthese des räumlich und zeitlich variablen Schwerefeldes mittels vektorieller Kugelflächenfunktio- nen / M. Lindlohr (1987)PermalinkDie Gaußsche Ellipsoid-Abbildung mit komplexer Arithmetik und numerischen Näherungsverfahren / H. Glasmacher (1987)Permalink