Géodésie

Transformation dans le Système R.G.F 93

de la base de données géographiques
de la Ville de Nice
Ludovic ANDRES

L’intérêt de travailler dans le nouveau système de Utilisation du système GPS

référence R.G.F.93 en vigueur en France est exposé
L’acquisition et la mise en œuvre de récepteurs GPS a mis en
à travers l’expérience de migration des données évidence une difficulté liée à son utilisation dans le cadre du
géographiques de la Ville de Nice : l’utilisation système N.T.F ou de sa densification locale : en effet, ce sys-
d’une méthode classique de changement de tème de positionnement permet d’obtenir aisément une pré-
système de coordonnées de l’ensemble de la base cision de l’ordre du centimètre quel que soit le lieu dès lors
de données géographiques rattachée au système que l’on travaille dans un référentiel précis et homogène.
L’établissement d’un canevas à l’aide de cette technologie
planimétrique N.T.F. est présentée. Cette technique
peut donc conduire à un réseau comportant de faibles défor-
permet de préserver au maximum la précision des mations. A contrario, la N.T.F. présente des déformations plus
données, et peut être mise en œuvre simplement à importantes et non homogènes. Convertir des coordonnées
l’aide de logiciels existants. Cette expérience peut obtenues par GPS dans la N.T.F. engendre donc une dégra-
être transposée à d’autres bases de données dation de la précision initiale des données. Les déformations
géographiques et doit inciter à migrer dans le de la N.T.F. étant importantes au niveau national, et non négli-
geables localement, il convient de les modéliser et d’en tenir
système de référence en vigueur.
compte dans le processus de transformation de coordonnées
afin d’insérer le point mesuré au mieux dans le réseau
déformé de la N.T.F. Une grille de conversion (GR3DF97A,

L
e décret d’application de l’article 53 de la loi
d’Aménagement et de Développement Durable (Décret 1997) modélisant ces déformations a été élaborée par l’IGN et
n° 2000-127, Loi n°95-115) a défini depuis le est disponible gratuitement sur son site internet. Cette grille
26 décembre 2000 le R.G.F.93 et sa projection associée le permet d’insérer les données dans la N.T.F. avec une erreur
Lambert 93 comme le nouveau système de référence en variant globalement entre 0 et 10 cm : elle convient donc a
vigueur. Ce système venait en remplacement de la Nouvelle priori pour transformer toutes les données de petite et
Triangulation de la France - N.T.F. (Levallois, 1988) dont la moyenne échelle, mais n’est pas adaptée pour les données à
faible précision de l’ordre de 1 cm/km, et des problèmes très grande échelle que sont les données topographiques et
d’échelle locaux rendaient inadapté aux exigences actuelles le canevas densifié de précision centimétrique. Cela nous
en terme de levés topographiques et générait également des imposait donc d’effectuer des mesures GPS complémentaires
difficultés liées à l’entretien du canevas (Kasser, 2001). Le nou- afin de modéliser la géométrie locale de notre canevas topo-
veau référentiel R.G.F.93 remédie à ces problèmes en offrant graphique dans le but de la prendre en compte dans le pro-
une précision accrue de l’ordre du centimètre sur l’ensemble cessus de transformation de coordonnées. Cette étape sup-
du pays. Il permet d’autre part de travailler de manière beau- plémentaire constituait d’une part une surcharge de travail et
coup plus simple et précise avec le système de positionne- d’autre part conduisait malgré tout, in fine, à une dégradation
ment par satellites GPS - Global Positioning system (Botton de la précision initiale des données. Cette solution n’était donc
et al., 1997). De très nombreux articles ont expliqué en détail pas satisfaisante et générait un sentiment de frustration chez
tous les avantages liés à ce nouveau référentiel, et l’on peut nos géomètres, attachés à effectuer de la belle ouvrage.
également obtenir simplement des informations à ce sujet sur
les sites internet de l’Institut Géographique National - IGN Maintenance et amélioration de la qualité
([Link] ou du Conseil National de l’Information du canevas
Géographique - C.N.I.G. ([Link]
Le canevas topographique de la Ville de Nice présente à
l’échelle de la commune des erreurs comprises entre 0 et 10
cm environ. Ces erreurs peuvent parfois générer des com-
Le choix du R.G.F.93 pensations importantes au regard des distances considé-
En ce qui concerne la Ville de Nice (VdN), différentes raisons rées. Il a donc été décidé d’exploiter au mieux les technolo-
ont conduit à adopter le système R.G.F.93. gies à notre disposition et de modifier nos méthodes de tra-

56 Revue XYZ • N°97 – 4e trimestre 2003

vail afin d’améliorer progressivement la qualité de notre La topographie
canevas. De plus le développement de notre canevas
La Ville de Nice dispose d’une base de données topogra-
conduit parfois à des travaux de polygonation importants et
phique dont l’élaboration a débuté dans les années 1920. Il
immobilisants. Tous cela nous a conduit à favoriser l’utili-
s’agissait à l’époque de plans topographiques relevés dans un
sation du système GPS et à opter pour un référentiel géo-
système local, nommé « Système Ville de Nice » et dont la tri-
désique plus précis.
angulation initiale de 122 sommets avait été conduite par
Monsieur Danger, ingénieur géomètre. Il est important de
Utilisation du R.G.F.93 et du R.G.P.
souligner que dès la création du service il avait été jugé néces-
L’Institut Géographique National ne maintient plus le réseau saire de rattacher tous les travaux à un même système de
N.T.F. qu’il a remplacé (Kasser, 2001) par le Réseau de Base référence planimétrique afin d’assurer un continuum géo-
Français - R.B.F. (1023 sites) formant l’accès opérationnel au graphique des données : le but recherché était d’en faciliter la
R.G.F.93 (23 sites). Les bornes détruites de l’ancien référentiel gestion (création, mise à jour, édition de plans), et de per-
ne sont donc plus remplacées. En sus du R.B.F., un Réseau de mettre de mieux coordonner les alignements de voirie sur
GPS Permanents (R.G.P.) a été développé par l’I.G.N et diffé- l’ensemble de la commune. Le dessin s’effectuait sur des
rents partenaires (Collectivités Territoriales, Universités, feuilles Invar de format A0.
CNRS, etc...). Il couvre l’ensemble du territoire français et per-
Avec la mise en place de la N.T.F. et la parution le 2 Juin 1948
met un accès aisé et précis au référentiel national R.G.F.93. Les
du décret interministériel rendant son application obligatoire,
travaux de canevas par GPS s’en trouvent donc simplifiés et
le Système Ville de Nice devenait obsolète. Aussi, au début
de qualité accrue.
des années 70 tous les plans ont été transformés manuelle-
ment dans le système N.T.F. L’objectif recherché était non seu-
Travaux topographiques de précision
lement de travailler dans un cadre légal, mais également de
D’importants projets d’infrastructure ont été lancés par la Ville faciliter les échanges de données avec les différents parte-
de Nice tels la construction d’une ligne de tramway traversant naires de la Ville de Nice (cabinets de géomètres, services de
la commune, ou le prolongement d’une autoroute urbaine. l’Etat, secteur privé, etc...). Tous les nouveaux levés étaient
Ces énormes chantiers sont soumis à des contraintes impor- donc réalisés directement dans la N.T.F. En réalité, il avait été
tantes en terme de précision des travaux topographiques affé- procédé à une densification locale de la N.T.F., densification
rents. L’exploitation du système G.P.S dans le R.G.F.93 per- qui s’est étoffée au fur et à mesure des nombreux travaux
mettait de répondre aux spécifications techniques liées à la topographiques.
topographie, de bâtir une ossature de canevas et d’intégrer
Les évolutions technologiques ont affecté naturellement les
toutes les données actuelles et futures dans la base de don-
nées topographique de la commune. Travailler dans l’ancien
système N.T.F. n’aurait aucunement permit d’atteindre simul-
méthodes de levé, mais également celles liées aux calculs
topographiques et aux dessins des plans : dès le début des ...
tanément ces objectifs. Les importants investissements réali-
sés en terme de topographie et liés à ces chantiers s’en trou-
vent ainsi capitalisés et les résultats obtenus pourront être
facilement exploités pour des études ou travaux ultérieurs.

Respect de la Loi
En tant qu’acteur local important, la Ville de Nice doit s’inté-
grer au cadre légal dans le véritable esprit de la loi, qui est de
s’appuyer sur le R.G.F.93 et d’abandonner la N.T.F.

De ce fait, en montrant l’exemple, et à travers les nombreux

échanges de données qu’elle réalise avec ses différents par-
tenaires, la Ville de Nice contribue à faire adopter ce nouveau
référentiel. Elle participe ainsi à l’effort national concernant la
gestion cohérente et homogène des informations géogra-
phiques sur l’ensemble du territoire.

Contexte et Historique
Avant de présenter la méthodologie utilisée pour réaliser la
conversion dans le système R.G.F.93, il est nécessaire de pré-
senter un bref historique ainsi que le contexte dans lequel Figure 1 : répartition géographique des relevés
s’opère cette transformation. topographiques sur la commune de Nice

Revue XYZ • N°97 – 4e trimestre 2003 57

 Géodésie

... années 80, des ressources informatiques ont été attribuées et

ont permis de mettre en place une base de données gra-
Le processus de transformation de la base de données peut
se décliner en deux étapes principales :
phique de plans topographiques. Environ 50 % des plans • il s’agit en premier lieu de déterminer la méthode de trans-
anciens jugés encore d’actualité ont été numérisés tandis que formation de coordonnées ainsi que ses paramètres associés
tous les nouveaux levés étaient directement saisis sur le logi- • en second lieu il conviendra d’appliquer la transformation
ciels de dessins qui fonctionnait alors sur un système central sur l’ensemble des couches constituant le SIGN à l’aide de
IBM. Aujourd’hui, la base de données topographique de la différents logiciels.
Ville de Nice s’étend sur plus de 80 % du corps de rue de la
commune. Elle comporte 1 300 000 points de semis et 13500
stations. La Figure 1 illustre la répartition géographique des Choix de la méthode de transformation et
relevés topographiques sur la commune de Nice. détermination des paramètres
Tous les points de la base de donnée sont renseignés d’une Comme on l’a vu précédemment, les données géographiques
altitude provenant d’un réseau de nivellement propre à la Ville de la Ville de Nice sont exprimées dans un système de coor-
de Nice, et qui constitue une densification de celui de l’IGN : à données correspondant à une densification du système de
ce propos, une conversion dans le système NGF - IGN69 référence géodésique N.T.F. Or nous souhaitons dorénavant
(Levallois et Maillard, 1970) à partir du système NGF - obtenir pour chacun des points et objets de la base, des coor-
Lallemand a eu lieu en 2002 (Andres, 2002) afin de travailler données dans le nouveau système R.G.F.93. C’est la transfor-
dans le système altimétrique en vigueur. La gestion entière- mation permettant d’exprimer les coordonnées de ces points
ment numérique des données s’effectue sous le logiciel dans un autre système géodésique (en l’occurrence dans le
Autodesk Map 5 en assurant un continuum géographique. R.G.F.93) qui constitue le changement de système géodésique.
Cette transformation peut s’appliquer sur des coordonnées
Les autres données du Système d’Information planes, géographiques ou bien cartésiennes. De nombreuses
Géographique méthodes mathématiques de transformation existent, et la
figure 2 (IGN, 1999) présente pour quelques modèles de trans-
Les données topographiques, ne forment qu’une partie du
formation un schéma synoptique permettant d’appréhender
Système d’Information Géographique de la Ville de Nice
le processus global de changement de système
(SIGN). Ce dernier comprend également de très nombreuses
autres informations dont l’acquisition et l’exploitation numé-
rique a débuté au début des années 1990. Aujourd’hui une
importante base de données regroupant plus 300 couches
correspondant à plus de 30 secteurs d’activité différents est
en service. La table 1 présente les principales couches d’in-
formation formant le SIGN.

Le fait que la base géographique représente un volume très

important (environ 60 Go), réparti en de très nombreuses
couches gérées dans différents formats correspondants à plu-
sieurs logiciels, génère des difficultés qu’il faudra appréhen-
der lors du changement de système de coordonnées. De plus,
Figure 2 : Processus de changement de système
son exploitation de manière continue et intensive aussi bien
géodésique : pour chaque type de coordonnées un
en consultation qu’en mise à jour ne saurait tolérer une inter-
modèle de transformation est présenté (source : Institut
ruption de service durable.
Géographique National)

Table 1 : Principales données utilisées et constituant le Système d’Information Géographique de la Ville de Nice
Echelle Vecteur Raster
Petite Scan 25 IGN, Scan 100 IGN
Moyenne Plan Urbain de la Ville de Nice
Périmètres administratifs : Canton, quartiers...
Bd Eclairage Publique Ville de Nice
Filaire de voie - Référentiel des adresses
Ilots et Iris INSEE
Grande Cadastre DGI précision Bd Ortho IGN
Plan d’Occupation des Sols Orthophoto VdN
Autres données Ville de Nice (habitat, nettoiement, sport, éducation, Orthophoto satellite 60cm
équipement, publics, réseaux, patrimoine, espaces verts, culture, etc...)
Très grande BD Topographique Ville de Nice et canevas géodésique
Alignements de voirie Ville de Nice

58 Revue XYZ • N°97 – 4e trimestre 2003

Les principales méthodes de transformation analysées à la Transformation standard N.T.F - W.G.S.84
Ville de Nice ont été les suivantes :
Il s’agit d’une transformation réalisant une translation en
• Transformation standard N.T.F.->W.G.S.84 avec les para-
3 dimensions entre les coordonnées géocentriques carté-
mètres uniques standards
siennes du système N.T.F. et celles du système W.G.S.84, ce
• Transformation en utilisant la grille de conversion planimé-
dernier système étant compatible au niveau métrique avec le
trique GR3DF97A
R.G.F.93. Les paramètres de translation (RTG n°14) sont illus-
• Transformation en utilisant une grille locale à la Ville de Nice
trés dans le système d’équations ci-dessous : pour inverser la
• Transformation par similitude 2D à 5 paramètres (transfor-
transformation, il suffit d’inverser le signe des valeurs de la
mation plane)
translation.
• Transformation par une similitude 3D à 7 paramètres
La précision moyenne de ce
modèle de transformation
Critères de choix
est d’environ 2 mètres, et
Il a été décidé, afin de conserver l’homogénéité de la base de peut aller jusqu’à 5 mètres
données, et de ne pas altérer sa topologie, de transformer selon les lieux (IGN, 1999).
l’ensemble des données avec le même modèle de transfor- Cela signifie que cette trans-
mation. En effet, de nombreuses couches de données ont été formation standard N.T.F.- W.G.S.84 est écartée d’emblée, car
créées en s’appuyant sur d’autres couches de référence : c’est l’erreur qui en résulterait ne serait pas compatible avec la pré-
le cas notamment du Plan d’Occupation des Sols (POS) saisi cision de nos données.
sur le plan cadastral ou bien des alignements établis à l’aide
des relevés topographiques. Transformation avec la grille de conversion
planimétrique GR3DF97A
De plus, la transformation retenue doit être en accord avec la
précision initiale des données : les couches de données les plus Cette grille mise au point par l’Institut géographique National
précises étant les données topographiques centimétriques, la est basée sur le même principe de transformation que celui
précision moyenne de la transformation doit donc s’en appro- présenté ci-dessus, c’est à dire une translation 3D entre sys-
cher. tème géodésique. Cependant, au lieu de disposer d’un seul
jeu de paramètres, de très nombreux jeux de paramètres (éta-
La mise en œuvre ultérieure de la transformation est également
blis par mesures GPS précises sur les sites R.B.F. / N.T.F.) ont
un critère très important : les différents logiciels SIG et DAO uti-
été interpolés dans une grille qui recouvre toute la France et
lisés devront permettre de réaliser la transformation retenue.
dont le pas est de 0,1° en latitude et longitude. La précision
Celle-ci devra être analysée sous deux angles distincts :
obtenue à l’aide de ce modèle de transformation est de
• Possibilité de réaliser la transformation de manière définitive
quelques centimètres (principalement de 0 à 10 centimètres
(la base de données est définitivement convertie)
selon les lieux).
• Possibilité de réaliser la transformation à la volée (les don-
Nous avons effectué des mesures GPS précises en Lambert
nées peuvent rester dans le système ancien et le logiciel les
93 des principales bornes N.T.F. de la commune de Nice. Nous
transforme à la volée au moment de l’affichage ou de l’ana-
avons ensuite comparé les coordonnées obtenues avec celles
lyse, dans le système choisi)
issues de la transformation des coordonnées N.T.F. de ces
La simplicité de la transformation est également à nos yeux
bornes à l’aide de la grille GR3DF97A et du logiciel CIRCE de
un facteur à prendre en compte : une transformation avec
l’IGN.
paramètres uniques conduit à une transformation en bloc
(chaque point est transformé de la même manière) : à préci- La table 2 présente les résultats obtenus pour ces sites ainsi
sion équivalente, elle conduira à une simplicité accrue et que les écarts constatés. On s’aperçoit, en ce qui concerne la
pourra permettre une réversibilité plus simple qui serait utile
si une période transitoire longue s’avérait nécessaire.
commune de Nice que les écarts en distance constatés entre
un point stationné “nativement” en R.G.F.93 et le même point ...
Borne
Ordre N.T.F N par GPS E par GPS N par CIRCE E par CIRCE DN DE
géodésique

Mt Chauve 1 6306012,125 1042471,044 6306012,145 1042471,055 0,020 0,011

Observatoire 1 6301373,633 1046291,064 6301373,574 1046291,139 0,059 0,075
Serena 3 6302898,705 1042601,069 6302898,671 1042601,095 0,034 0,026
Albert1er 3 6297711,936 1043971,146 6297711,807 1043971,163 0,129 0,017
Mt Alban 3 6298538,279 1046500,483 6298538,213 1046500,572 0,066 0,089
Cagnes 3 6292515,762 1035105,736 6292515,783 1035105,507 0,021 0,229
Moyenne 0,055 0,074

Table 2 : Evaluation de la précision de la grille GR3DF97A sur les bornes géodésiques N.T.F. de la Ville de Nice
(système de coordonnées : R.G.F.93/Lambert93)

Revue XYZ • N°97 – 4e trimestre 2003 59

 Géodésie

... transformé à partir de ses coordonnées N.T.F. s’échelonnent

entre 2 et 23 centimètres, avec une moyenne de 9 centimètres
environ. Cette précision constatée sur quelques points
confirme la précision moyenne de 5 à 10 cm annoncée par
l’IGN pour cette grille de transformation. Nos données topo-
graphiques étant de précision centimétrique, cette méthode
de transformation est donc insuffisante car elle génèrerait une
perte de précision de nos données.

Figure 3 : principe d’une similitude 3D entre deux

Transformation en utilisant une grille locale repères de l’espace
à la Ville de Nice
une rotation et à appliquer un facteur d’échelle, afin de passer
Le principe serait d’établir une grille de transformation sem-
d’un repère orthonormé de l’espace à un autre repère de
blable à la grille GR3DF97A, mais adaptée à notre géodésie
même type (figure 3)
locale. Cela aurait pour effet d’obtenir une précision améliorée,
puisque la grille serait calculée à partir de données mesurées Soit M un point de l’espace dont les coordonnées sont (X1, Y1,
par GPS et réparties sur l’ensemble de la commune. Z1) dans le repère (O1, i1, j1, k1) et (X2, Y2, Z2) dans le repère (O2,
Cependant, une analyse approfondie de notre canevas local, i2, j2). On peut écrire l’équation matricielle suivante :
permise grâce aux nombreuses mesures GPS effectuées dans
le système R.G.F.93, a permis de constater qu’il présente des
déformations en amplitude et en orientation qui varient de
façon très importante selon les lieux (figure 6). Pour que la grille
où les 7 paramètres de la transformation sont :
locale envisagée modélise de manière cohérente les déforma-
tion du canevas, il faudrait densifier de manière très importante • (TX, TY, TZ) qui représentent le vecteur translation entre le

• εX, εY, εZ, sont les rotations autour respectivement des axes
les observations à réaliser, pour qu’elles soient en adéquation repère (O1, i1, j1, k1) et le repère (O2, i2, jv, k2).
avec les polygonales réalisées depuis plus de vingt ans, poly-
gonales qu’il conviendrait d’identifier précisément. Outre un (O1, i1), (O1, j1), (O1, k1).
travail de préparation colossal, cela engendrerait une masse • E est le facteur d’échelle, défini de telle sorte que
d’observations très importante. Aussi cette solution n’a pas été ||i2||=||j2||=||k2||=E.||i1||

Les rotations εX, εY, εZ étant faibles, on peut faire l’approxi-

retenue : l’amélioration de précision que l’on obtiendrait n’étant
pas en rapport avec l’investissement en terme de temps passé
mation :
et de coût. De plus les versions de logiciels DAO et SIG utilisées
actuellement à la Ville de Nice, ne supportent pas encore une
transformation avec ce type de grille de conversion.

Transformation par similitude 2D à 5 paramètres

(transformation plane)
L’équation précédente peut s’exprimer de la manière suivante :
Cette méthode de transformation consiste à trouver le meilleur
jeu de paramètres pour une similitude à 5 paramètres (trans-
lation en X et Y, rotation en X et Y, facteur d’échelle), entre des
coordonnées planes de la N.T.F. (Lambert III). et les coordon-
nées planes du R.G.F.93 (Lambert 93). Le problème rencontré
avec cette méthode est que peu de logiciels l’implémentent.
Aussi, eu égard au contexte logiciel de la Ville de Nice, il était
difficile de mettre en application cette transformation à l’aide de
nos outils DAO et SIG. Il aurait fallu exporter toutes les données Pour un point considéré de l’espace, ce système comporte
pour les convertir à l’aide d’un processeur intégrant ce modèle 3 équations et 7 inconnues. Avec au minimum 3 points de
de transformation, et les réintégrer ensuite à leur format natif l’espace connus dans les deux systèmes on pourra donc cal-
dans leur logiciels respectifs. De plus, les calculs de paramètres culer les paramètres.
que nous avons réalisés selon ce modèle de transformation
En pratique, un nombre important de points vont être mesu-
présentaient des résultats moindres en terme de précision que
rés dans les deux systèmes et vont conduire à une sur-déter-
ceux obtenus avec la méthode évoquée ci-dessous.
mination de ce système d’équations : pour N points mesurés,
nous aurons 3N équations et 7 inconnues. Du fait des incerti-
Transformation par une similitude 3D
tudes liées aux mesures, les 3N équations apparaissent
à 7 paramètres uniques
comme étant toutes indépendantes. On peut résoudre de tels
Cette transformation connue également sous le nom de trans- systèmes linéaires par la méthode de l’approximation des
formation de Bursa Wolf consiste à effectuer une translation, moindres carrés :

60 Revue XYZ • N°97 – 4e trimestre 2003

en effet, si i est l’indice du point observé dans les deux sys- Le meilleur moyen d’accéder rapidement et avec précision à
tèmes (1<=i <=N), et en posant la référence R.G.F.93 est d’utiliser des techniques de position-
nement par GPS. La Ville de Nice est d’ailleurs équipée en
récepteurs GPS bi-fréquence, et possède une station perma-
nente GPS (Figure 4) intégrée au R.G.P.
Le territoire de la commune s’étend sur plus de 7200 hectares
avec un relief de type collinaire comportant des vallons orien-
tés Nord Sud. Les altitudes varient de 0 à 500 mètres. Nous
recherchons une précision d’environ un centimètre sur les
composantes planimétriques des coordonnées R.G.F.93. Pour
cela deux ou trois sessions d’observation GPS de 3 heures
sont prévues en mode statique (les résultats seront égale-
ment utilisés ultérieurement lors de la mise au point d’une
et le vecteur des 7 paramètres : grille de conversion altimétrique Hauteur Ellipsoïdale ->
Altitude IGN69).
Un maillage de 1,5 km a été appliqué sur la commune condui-
sant à une quarantaine de points observables en régie. Ces
points ont été choisis parmi les 13500 stations existantes sur
le terrain.
Le système que l’on cherche à résoudre peut s’écrire : M.p=b
Chaque point a été étudié avec précaution de manière à
Pour calculer le vecteur p (le meilleur jeux de 7 paramètres)
répondre aux critères suivants :
il faut minimiser la norme du vecteur erreur e=M.p-b, c’est à
• situation géographique proche de la position théorique don-
dire minimiser le carré de l’erreur quadratique moyenne mul-
nées par l’application du maillage
tipliée par 3N :
• le minimum de masque pour les observations GPS
3N.e2 = | M.p-b |2 = ( M.p-b )T . ( M.p-b )
• déport (et rattachement en Lambert III du point déporté) de
Le système M.p=b n’a pas de solution car il contient moins
la station si nécessaire pour éradiquer les masques
d’inconnues que d’équations. Par contre, le système
• contrôle de la qualité de chaque station (consultation du cal-
[Link]=MT.b, de 7 équations à 7 inconnues a en général une
cul de sa polygonale, et du levé topographique environnant)
solution si le déterminant de MT.M n’est pas nul.
• vérification sur le terrain de chaque station : pérennité de la
Il suffit donc de résoudre ce système pour déterminer le
station
meilleur jeu de 7 paramètres.
• nivellement de précision des stations à partir de repères de
nivellement IGN lorsque cela est possible (utile pour l’éta-
Détermination des paramètres de transformation
blissement d’une future grille de conversion altimétrique)
de la similitude 3D
Le nombre de points à observer étant conséquent, les obser-
Comme nous venons de le voir, il nous faut disposer de points
vations ont été rigoureusement planifiées et la stratégie d’ob-
de canevas possédant des coordonnées dans les deux sys-
servation affinée (choix des constellations favorables, choix
tèmes de référence N.T.F. et R.G.F.93. Les données topogra-
des antennes après vérification/calibration, stratégie d’occu-
phiques étant les données de référence les plus précises dis-
pation, disponibilité des moyens matériels et humains, etc...).
ponibles à la Ville de Nice, c’est cette couche d’informations
La figure 5 illustre une session d’observation réalisée en bord
qui a été utilisée pour établir les paramètres du modèle de
de mer.
transformation.
Toutes les calculs ont été effectués par rapport au pivot de
La base de données topographiques est originellement dans base que représente notre station permanente NICA, située
le système de référence N.T.F. Il ne reste qu’à déterminer pour géographiquement à moins de 10 km de tous les points à sta-
un certain nombre de points à choisir judicieusement leurs
coordonnées en R.G.F.93.
tionner. Plusieurs récepteurs GPS réalisant les mesures simul-
tanément, un calcul de fermetures internes et externes ...

Figure 4 : Photographies et architecture fonctionnelle de la station permanente GPS de la Ville de Nice (NICA).

Revue XYZ • N°97 – 4e trimestre 2003 61

 Géodésie

Point Résidu Résidu Résidu en

géodésique Est (cm) Nord (cm) distance (cm)
Albert 1er-D1 -3,79 -8,13 8,97
Brancolar 0,92 -2,32 2,50
Carras -2,25 1,77 2,86
Cimetiere de l'est 5,83 -0,67 5,87
Corniche des oliviers 2,34 1,12 2,59
Raiola -8,20 4,53 9,36
Durandy 0,75 2,29 2,41
Ecole Bellanda -1,69 1,37 2,18
Fabron 2,71 1,88 3,30
Ginestiere 3,98 -3,52 5,31
Jeanne D'arc -0,23 6,87 6,87
Figure 5 : Observation d’un point géodésique situé sur la Leclerc 4,09 -1,32 4,30
Les combes 2,95 -0,50 2,99
Promenade des Anglais - Matériel utilisé : Antenne
Lingostiere -0,72 -5,03 5,08
ASHTECH 701945E_M, Récepteur THALES 6502

... (Botton et al, 1997) était systématiquement effectué à l’issue

des traitements en mode statique. Certaines lignes de base
Madeleine-D1
Maeterlink
Magnan
4,42
-2,93
0,72
2,56
4,43
-0,96
5,11
5,31
1,20
présentant des incohérences entre sessions successives, il a Marmorini 2,17 -4,09 4,63
Mont Alban 2,94 -1,92 3,51
été nécessaire d’effectuer des observations complémentaires
Montreal -1,45 -3,60 3,88
afin de lever certains doutes. Les résidus planimétriques obte-
MtChauve-b -2,86 14,22 14,51
nus par cette méthode de transformation en bloc sur chacun NICA 7,78 1,56 7,93
des points sont présentés en table 3. Nikaï 5,97 -1,86 6,25
On peut s’apercevoir que les écarts en distance varient de Observatoire-D1 2,80 2,09 3,50
2 cm à 9 cm environ si l’on ne considère pas le point du Mont Pont du Var -3,23 -0,04 3,23
chauve qui semble présenter un résidu anormal : ce point (1er Pontremoli 6,31 -0,08 6,31
ordre N.T.F.), n’est d’ailleurs pas en bon accord avec l’autre Ponts jumeaux -3,07 -1,47 3,40
point N.T.F. 1er ordre de la commune (Observatoire). Port-D2 -2,56 -7,82 8,22
L’écart moyen de 5 cm constaté correspond à la précision Rimiez Lombardie -0,05 1,62 1,62
attendue de notre canevas, et l’écart type de 2,7 cm autour de Riquier -2,88 -0,71 2,97
Serena -2,14 6,00 6,37
cette moyenne peut être considéré comme faible.
St pancrace -2,89 6,76 7,35
La figure 6 propose quant à elle, une représentation gra-
St pierre de Feric -4,97 -4,71 6,85
phique de ces écarts : le caractère anisotrope des orientations St Roman -4,04 5,42 6,75
et des amplitudes des résidus apparaît de manière évidente. Verola 2,18 -3,24 3,90
Vinaigrier Gde Corniche -7,12 0,46 7,14
Moyenne des |Résidus| 3,14 2,86 5,13
Ecart type 3,85 4,38 2,68

Table 3 : Liste des résidus obtenus sur les points N.T.F.

observés par GPS dans le système R.G.F. 93. Ces résidus
sont obtenus en utilisant le modèle de transformation de
Bursa Wolf (similitude 3D à 7 paramètres) calculé.
Cela justifie pleinement le fait d’utiliser une transformation en
bloc comme celle fournie par le modèle de Bursa Wolf avec
un jeu unique de paramètres. La mise au point d’une grille de
transformation locale de type GR3DF97A, comme évoqué
précédemment, aurait nécessité une très importante densifi-
cation pour modéliser de manière réaliste les déformations
de notre canevas.
Les paramètres calculés pour passer du système N.T.F. den-
sifié Ville de Nice au R.G.F. 93 sont les suivants :
Tx = -166,1092 m Ty = -59,6894 m Tz = 321,8207 m
Rx = 0,27702 “ Ry = 0,014708 “ Rz = 0,24028 “
E = 1-0,41527. 10-6
L’établissement de ces paramètres aura finalement nécessité
Figure 6 : Représentation graphique des résidus issus approximativement :
du modèle de transformation 3D à 7 paramètres pour • 100 heures de préparation, repérage, matérialisation et rat-
les points géodésiques observés du canevas tachement éventuel des points

62 Revue XYZ • N°97 – 4e trimestre 2003

• 300 heures d’observation Outil DAO
• 50 heures de vidage et calcul des observations
Le logiciel Autodesk Map 5 présente des fonctionnalités SIG
• 50 heures d’études, mise en forme et calcul des paramètres.
permettant notamment d’interroger (visualisation ou conver-
Cela représente au total environ 500 heures de travail. sion définitive) un dessin dans un système différent de celui
Le nombre de sessions ainsi que la durée de chaque session dans lequel se trouvent les données originales. Ces transfor-
aurait pu être minimisé, mais effectuant ce travail de manière mations tiennent compte du changement de système géodé-
unique, nous avons souhaité prendre le maximum de pré- sique par la méthode de Molodensky ou bien par une trans-
cautions quant à la qualité des coordonnées obtenues par formation à 7 paramètres de Bursa Wolf. Toutes les projec-
GPS. De plus nous souhaitions obtenir une bonne précision tions et ellipsoïdes sont également paramétrables. Après
sur la hauteur ellipsoïdale afin de réutiliser ces résultats pour avoir effectué de nombreux tests et paramétrages pour les
analyser la grille de conversion altimétrique R.A.F.98 projections nous concernant, et les transformations géodé-
(Duquenne, 1998), ou en établir une localement si celle-ci siques souhaitées nous ne sommes pas parvenu à obtenir
s’avérait insuffisante. des résultats satisfaisants. Le problème a été signalé à
Autodesk qui reconnaît que ces fonctionnalités ne sont pas
pleinement opérationnelles en version Autodesk Map 5, et
Mise en œuvre de la transformation nous oriente sur les versions Autodesk Map 6 et ultérieures.
sur la base de données géographiques D’autre part aucun support de la grille de conversion
GR3DF97A n’existe dans cette version 5.
Les données de notre système d’information géographique
sont gérées principalement à l’aide des logiciels ArcView Les versions 6 et ultérieures présentent les mêmes fonction-
(ESRI) et Autodesk Map (AUTODESK). Nous allons étudier nalités que la version 5 sur ce sujet avec en plus le support de
successivement les possibilités offertes par ces produits pour la grille GR3DF97A. Après des paramétrages adéquats, il
mener à bien nos transformations. s’avère que les changements de systèmes géodésiques fonc-
tionnent parfaitement bien, notamment en ce qui concerne la
Outil SIG transformation par similitude à 7 paramètres de type Bursa
Wolf. Nous avons également testé avec succès la transfor-
Le logiciel ArcView 3.3 comporte de manière native des outils
mation utilisant la grille de conversion GR3DF97A (après para-
de projection mais ne comporte pas d’outils permettant de
métrage du symbole décimal comme caractère point).
changer de système géodésique. Nous n’avons pas trouvé
Il est nécessaire de préciser que ces transformations ne
sur le site d’ESRI France d’outils complémentaires permettant
concernent que les données de type vecteur manipulées par
de réaliser cela.
le produit Autodesk Map.
Pour pouvoir travailler dans le système R.G.F.93 (en
D’autre part il n’y a aucune possibilité d’exploiter ces fonc-
Lambert93 par exemple) il convient donc de transformer les
tionnalités pour des données de type raster : toutes les don-
données une fois pour toute dans ce système de référence, à
nées raster doivent donc obligatoirement être converties de
l’aide d’un autre outil.
manière définitive avec un outil spécifique dans le système de
Le logiciels Arcview 8.3 présente de nombreuses fonctionna- référence choisi.
lités concernant l’exploitation simultanée de données qui sont
Nous sommes donc forcés de convertir définitivement dans
dans des systèmes de référence différents.
le système de référence choisi les données :
La transformation par similitude à 7 paramètres de type Bursa
• Raster (à cause des outils DAO utilisés)
Wolf fonctionne à la volée ou permet une transformation défi-
• DAO (à cause des outils SIG utilisés)
nitive des données vecteur. Elle fonctionne également sur les
données raster mais uniquement à la volée et ne permet pas Nous avons donc décidé de transformer de manière définitive
de transformation définitive sur ce type de données sans le l’ensemble de nos données géographiques, qu’elles soient de
recours à une extension spécifique. Enfin, elle ne fonctionne type raster, ou vecteur (SIG et DAO)
pas sur les données CAO/DAO Autocad (DWG-DXF) ou Cette transformation a été planifiée sur une période d’une
Intergraph (DGN). semaine. Toutes les données dont la fréquence de mise à
jour est faible ont été préalablement transformées (données
Cette version 8.3 du logiciel ne supporte pas la grille de
raster et vecteur). Les autres ont été identifiées et transfor-
conversion GR3DF97A : ESRI France nous a donné la possi-
mées rapidement en concertation étroite avec les services
bilité d’évaluer la version 9 d’Arcview qui sera commerciali-
utilisateurs. Des actions d’information ont eu lieu réguliè-
sée courant 2004, pour laquelle nous avons pu constater le
rement depuis environ une années, et des séances de for-
support de la grille GR3DF97A ainsi que son bon fonctionne-
mations sont organisées pour les correspondants des dif-
ment.
férents services utilisateurs. Enfin, vis à vis des profession-
Puisqu’à la ville de Nice nous travaillons simultanément avec nels (géomètres, entreprises, acteurs locaux, ...), des for-
des données SIG et DAO, nous serons donc forcés de conver- mulaires d’information ont été envoyés par courrier, et une
tir nos données DAO une fois pour toute dans le système sou- réunion d’information sera organisée début 2004 sur l’inté-
haité pour pouvoir utiliser ces données DAO à l’aide du logi-
ciel SIG Arcview.
rêt de ce nouveau référentiel et sa mise en œuvre à la Ville
de Nice. ...
Revue XYZ • N°97 – 4e trimestre 2003 63
 Géodésie

... Conclusion Remerciements

L’adoption du nouveau référentiel R.G.F 93 apparaît comme un Je tiens à remercier tout spécialement mon collègue Bernard
choix inévitable et sage afin de rester en adéquation avec les Laugier, responsable des brigades topographiques, pour son
méthodes de travail modernes et le degré de précision qu’elles implication personnelle, le travail de grande qualité qu’il a fourni
permettent. La seule difficulté, de taille malgré tout, consiste et l’aide indispensable qu’il a apporté à ce projet. Je félicite
en la transformation des données existantes. Ce travail impor- d’autre part tous les topographes de la Direction Information
tant et coûteux doit cependant être réalisé : la difficulté de cette Géographique (brigades de topographie et des alignements)
tache est liée à la précision et à la complexité des données à pour leur contribution active à toutes les observations GPS qui
transformer. D’ores et déjà notre étude montre qu’il peut être ont été réalisées ainsi que les collègues qui ont participé à la
aisé de transformer des données de précision décimétrique à transformation de la base de données géographiques. Je remer-
l’aide de modèles établis tel la grille de conversion GR3DF97A cie enfin Messieurs Michel Kasser (IGN) et Denis Delerba (Ville
qu’il est possible de mettre en œuvre simplement à l’aide de de Nice) pour leurs commentaires, Monsieur Alain Harmel (IGN)
logiciels existants. Pour des données de haute précision, notre pour ses conseils, ainsi que les sociétés AUTODESK France,
expérience prouve que cela est réalisable dès lors que la ESRI France, THALES Navigation, AEC Informatique (St Laurent
volonté de travailler dans le système légal existe. du Var, 06), Digitech International (Bagneux, 92), et Latitude
Peu de témoignages sur ce type d’expérience existent en Géosystèmes (Deuil la Barre, 95) pour leur collaboration.
France. Aussi le changement de système de coordonnées réa-
lisée ici doit avoir une action de démystification : elle doit inci- Références
ter les gestionnaires de bases de données géographiques tels ANDRES L., 2002, Conversion dans le système altimétrique IGN69
des collectivités territoriales, des services de l’Etat, des de la base de données topographiques de la Ville de Nice, XYZ
concessionnaires de réseaux, ou des cabinets de géomètres n°91, Juin 2002
demeurés dans la N.T.F. ou dans des systèmes locaux à ARRÊTÉ INTERMINISTÉRIEL du 2 juin 1948
migrer dans le système légal R.G.F.93. ●
fixant les conditions d’exécution et de publication des levés de
plans entrepris par les services publics, Journal Officiel.
Contacts BOTTON S., DUQUENE, F., EGELS Y., EVEN M, WILLIS P., 1997,
Ludovic ANDRES - Ingénieur en Chef GPS localisation et navigation, collection. CNIG-GPSD, éditions
Direction Information Géographique Hermes.
membre du GT SIG TOPO de l’AITF DECRET N° 2000-1276, 26 Décembre 2000, Décret portant
Mairie de Nice application de l’article 89 de la loi n° 95-115 du 4 février 1995
3 rue de la Terrasse -06364 NICE Cedex 4 modifiée d’orientation pour l’aménagement et le développement
Tel : 04 97 13 25 53 du territoire relatif aux conditions d’exécution et de publication
Fax : 04 97 13 29 24 des levés de plans entrepris par les services publics
Email : [Link]@[Link] DUQUENNE, H., 1998, QGF98 - a new solution for the quasi-
geoid in France, Proceedings of the Second Continental Workshop
Bernard Laugier
on the Geoid in Europe. - Reports of the Finnish Geodetic
Direction Information Géographique Institute, [Link] and J. Adam, editors, Masala, Finland.
[Link]@[Link]
GR3DF97A, 1997, Grille de paramètres de transformation de
Denis Delerba coordonnées - notice d’utilisation, Service de Géodésie et de
Responsable de la Direction Information Géographique Nivellement, Institut Géographique National.
Animateur National GT SIG TOPO de l’AITF IGN, 1999, Document d’information générale sur la transforma-
[Link]@[Link] tion de coordonnées, Institut Géographique National, Service de
Géodésie et de Nivellement.
ABSTRACT KASSER, M., 2001, Systèmes R.G.F. 93 et Lambert 93 : de petits
problèmes à court terme, de nombreuses simplifications à long
The City of Nice operates a large geographic database in terme, Géomètre, avril 2001.
a former french geodetic reference system named N.T.F. LEVALLOIS, J.J., 1988, Mesurer la Terre - 300 ans de géodésie
that has been replaced by a new one, the R.G.F.93. The française, Association Française de Topographie, Paris, 1988.
advanges of working in the R.G.F.93 system are pointed
LEVALLOIS J.J. AND MAILLARD J.,1970, Le nouveau réseau de
out through the data migration project that has been
nivellement de 1er ordre du territoire français, conséquences pra-
carried out. The coordinate conversion process that has
tiques et scientifiques. IGN
been selected aims to preserve the data precision, and
its implementation involving the use of common LOI N° 95-115, 4 Février 1995, Loi d’orientation pour l’aména-
algorithms and softwares is presented. Finaly, such gement et le développement du territoire, article 89
geographic database coordinate conversions are RT/G n°14, Nouveaux systèmes géodésiques utilisables en France
encouraged in order to support the R.G.F.93 system. (WGS84, ED87), Rapport Technique, Institut Géographique
National.

64 Revue XYZ • N°97 – 4e trimestre 2003