COURS DE CALCULS

TOPOMETRIQUES

INGENIEUR TINA

1 /77
 SOMMAIRE
GENERALITES......................................................................................................................................................4
I. DEFINITION..............................................................................................................................................4
II. ORGANISATION DES CALCULS.................................................................................................................4
1. Généralités :........................................................................................................................................4
2. Organisation des calculs......................................................................................................................4
III. UNITES UTILISEES EN TOPOMETRIE.....................................................................................................5
1. Unités linéaires....................................................................................................................................5
2. Unités angulaires.................................................................................................................................6
3. Unités de surface (unités agraires)......................................................................................................6
4. Conversions angulaires........................................................................................................................6
CHAPITRE I : CALCUL DE DISTANCES ET PROBLEME D’ORIENTATION.................................................................8
I. CALCUL DE DISTANCES : RESOLUTION DE TRIANGLES.............................................................................8
1. Triangle rectangle................................................................................................................................8
2. Triangle quelconque............................................................................................................................8
II. LES DIFFERENTES ORIENTATIONS............................................................................................................9
1. Gisement d’une direction..................................................................................................................10
2. G0 de station.....................................................................................................................................11
CHAPITRE II : CALCUL D’UN CANEVAS PLANIMETRIQUE...................................................................................13
I. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR INTERSECTION.........................................13
1. Formules simples...............................................................................................................................13
2. Formule globale................................................................................................................................14
3. Formule de Hatt................................................................................................................................14
II. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR RELEVEMENT...........................................15
1. Principe.............................................................................................................................................15
2. Méthode géométrique ou relèvement italien...................................................................................16
3. Méthode barycentrique....................................................................................................................17
4. Méthode de Delambre......................................................................................................................18
III. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR DOUBLE RELEVEMENT.........................20
IV. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR RECOUPEMENT...................................22
V. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR TRILATERATION.......................................23
CHAPITRE III : CALCULER ET COMPENSER UN CHEMINEMENT PLANIMETRIQUE OU UNE POLYGONALE.........24
I. DEFINITIONS ET PRINCIPE.....................................................................................................................24
1. Un cheminement en antenne ou ouvert...........................................................................................24
2. Le cheminement encadré..................................................................................................................24

2 /77
 3. Le cheminement dit fermé ou bouclé :.............................................................................................25
4. Le point nodal :..................................................................................................................................25
II. CALCULS................................................................................................................................................25
1. Calcul d’un cheminement encadré....................................................................................................26
2. Calcul d’un cheminement bouclé ou fermé......................................................................................30
CHAPITRE IV : CALCUL DE SURFACES.................................................................................................................37
INTRODUCTION.............................................................................................................................................37
I. CALCUL DE SURFACES PAR DECOMPOSITION EN FIGURES GEOMETRIQUES SIMPLES..........................37
II. CALCUL DE SURFACES PAR COORDONNEES RECTANGULAIRES.............................................................39
III. CALCUL DE SURFACES PAR COORDONNEES POLAIRES (ANGLES, DISTANCES)..................................40
IV. CALCUL DE SURFACES PAR LA METHODE POLYGONALE (OU METHODE DE SARRON)......................42
V. REDRESSEMENT DE LIMITES..................................................................................................................44

3 /77
 CHAPITRE I : CALCUL DE DISTANCES ET PROBLEME
D’ORIENTATION
I. CALCUL DE DISTANCES : RESOLUTION DE TRIANGLES
Résoudre un triangle consiste à calculer ses données manquantes (angles ; côtés ; surface) tout
en faisant des contrôles.

1. Triangle rectangle
A
a=c . sinA
b=a . tanB=c . sinB
b c a b
c= =
sin A sin B
1 1 1
C S= ab . sinC = ac . sinB= ab . sinC
a B 2 2 2
2. Triangle quelconque
1er cas : Un côté et deux angles connus : Données : a ; C ; B
A A=200−(B+C)
a b c
2 R= = = (triangle inscrit dans un cercle)
b c sin ⁡( B+C ) sin ⁡B sin ⁡C
a²
2 S=
cotanB +cotanC
C a B 2 S=ab . sinC=ac . sinB =cb . sinA

2ème cas : Deux côtés et un angle connus : Données : c ; b ; A

A hC = [Link]
BC’ = c-
hc b . sinA
tanB= tanB=
[Link]
c BC ' c−b . cosA
b B’ hB
h B= [Link]
C’

BC’ = b-
hc hB
tanC= c . sinA
tanC=
[Link]
CB ' b−c . cosA
A+B+C = 200 gr
C a B
Contrôle
Remarque : si tan est inférieur à 0, alors l’angle est égal à + 200 gr

2S =
a b c
= =

[Link]
sin ⁡A sin ⁡B sin ⁡C
ème
cas : Trois côtés connus : Données : a ; b ; c
a² = b²+c²-[Link]
A 1ère méthode :
b² = a²+c²-[Link]
b c c² = a²+b²-
[Link]
2ème méthode : p =
C a B a+b+ c
2
4 /77
 A
2
=
√( p−b ) ( p−c)
sin Application
bc :
B
sin =
2 √
( p−a ) ( p−c)
ac
sin
C
2
=
√( p−a ) (p−b) S= √ p∗( p−a ) ( p−b )( p−c )
ab
N° 1 : Résoudre le triangle ABC sachant que B = 69.894 gr ; C = 51.312 gr ; BC=a = 315.712 m

N°2 : Résoudre le triangle ABC sachant que BC=a = 224.55 m ; AC=b = 251.86 m ; AB=c = 412.29
m

N°3 : Résoudre le triangle ABC sachant que BC=a = 151.46 m ; AC=b = 212.28 m ; C = 28.654 gr

N°4 : Résoudre le triangle ABC sachant que AC= b = 49.12 m ; BC= a = 32.81 m ; S = 281.52 m²

N°5 : Résoudre le triangle ABC sachant que AC= b = 51.02 m ; A = 127.100 gr ; S = 432.83 m²

N°6 : Résoudre le triangle ABC sachant que A= 51.200 gr ; B = 121.720 gr ; C= 27.080 gr ; S =

2989.12 m²

II. LES DIFFERENTES ORIENTATIONS

Il existe trois nord qui sont le nord géographique ; le nord magnétique et le nord de la
représentation plane.

- L’angle que fait le nord magnétique avec la direction AB est appelé azimut magnétique
( AZm )
NM

B
AZm

A
- L’angle que fait le nord géographique avec la direction AB est appelé azimut
géographique (AZg)
NG

AZg B

A
- L’angle que fait le nord de la représentation plane avec la direction AB est appelé
Gisement noté souvent G ou V.
Y
B
GAB =VAB

A X

5 /77
 1. Gisement d’une direction
Le gisement d’une direction est l’angle orienté compris entre le nord de la représentation plane
et la direction considérée. Il est compté positivement à partir des Y positifs dans le sens des
aiguilles d’une montre. Il varie de 0 à 400 gr.

a) Le gisement peut être obtenu par le calcul si on connaît les coordonnées des extrémités de la
visée AB.
Y
tan GAB =
B B} over {B A
B’’ BB” = XB-XA =∆ XAB
B
YB
G B”A = YB-XA =∆ YAB
YA X B −X A ∆ X AB
tanG AB = =
A B’ Y B −Y A ∆ Y AB
XA XB
∆ X AB
GAB = arctg ( )
(0 gr) ∆ Y AB
Y

Quadrant IV : ∆ X<0 et ∆ Y>0 Quadrant I : ∆ X>0 et ∆ Y>0

∆X ∆X
g= arctg ( ¿<0 g= arctg ( ¿>0
∆Y ∆Y
GAB = g+400gr GAB =g B1
B4

(300gr) X (100gr)
A

B3 B2
Quadrant III : ∆ X<0 et ∆ Y<0 Quadrant II : ∆ X>0 et ∆ Y<0
∆X ∆X
g= arctg ( ¿>0 g= arctg ( ¿<0
∆Y ∆Y
GAB = g+200gr GAB = g+200gr
(200gr)

Tableau récapitulatif

N° quadrant ∆X ∆Y ∆X G
g= arctg ( ¿
∆Y
1 + + g>0 G=g
2 + - g<0 G=g +200gr
3 - - g>0 G= g +200gr
4 - + g<0 G=g + 400gr

6 /77
Application :
Pts X(m) Y(m)
A 789 050.55 313 030.90
B 786 786.62 315 309.88
C 786 450.78 313 065.25
D 791 753.90 311 665.23

Calculons les distances et les gisements AB ; BC ; CD ; AC ; AD et BD.

Résolution
∆X ∆Y D (m) ∆X G
g= arctg ( ¿
∆Y
AB
BC
CD
AD
AC
BD

Exercice :

Pts X(m) Y(m)

A 125 645.77 299 489.56
B 289 789.89 297 583.23
C 134 498.47 289 445.78
D 154 482.01 294 578.42

Calculer les gisements et les distances AB ; BC ; CD ; AC ; AD et BD.

b) Le gisement peut être obtenu à partir des observations :

GAM = G0 + Lm

2. G0 de station
Le G0 d’une station est le gisement de la graduation 0 du limbe.

Gisement du zéro du limbe pour une direction

Symbole : G0A
Le gisement du zéro du limbe est obtenu à partir du gisement
d'une direction connue et de la lecture sur le limbe horizontal
correspondant.
Point de station S (XS ; YS) ; Point visé connu A (XA ; YA)
GSA =arctg ¿

Lecture sur le limbe relative à A : LA

G0A = GSA - LA

7 /77
Gisement du zéro du limbe pour un tour d'horizon :
Symbole : G0 ou G 0moyen
Gisement moyen du zéro, origine de la graduation du limbe après réduction du tour d'horizon.
Il s'obtient par la moyenne arithmétique des gisements des zéros relatifs à plusieurs visées sur
des points connus.

G0A=GSA-LA, G0B=GSB-LB, G0C=GSC-LC

G 0 A + G0 B +G 0C
G 0moyen =
3
n

∑ G 0i
G 0moyen = i=1
n

Utilité du G 0moyen:

Le G 0moyen est le gisement qui, ajouté aux lectures

réduites d'un tour d'horizon, donne les gisements des
visées :

GSM =G 0 moyen+ LM (LM lecture réduite sur M)

Exercice

Station Points visés Lectures azimutales X (m) Y (m)

18 285 152.36 843 277.34
14 0.000 277 697.47 845 224.20
11 24.483 276 228.59 849 881.14
17 76.002 284 212.37 850 968.62
9 363.333 274 326.64 839 683.62

a) Calculer le G 0moyen de la station.

b) Calculer le gisement G18-8 sachant que la lecture azimutale réduite sur 8 est L8 =
398.578 gr.
c) Calculer les coordonnées rectangulaires du point 8 sachant que D18-8 = 315.45 m

8 /77
 CHAPITRE II : CALCUL D’UN CANEVAS PLANIMETRIQUE

I. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR

INTERSECTION

L’intersection est la méthode topographique qui permet de déterminer un point en le visant à

partir d’au moins 2 points.

M
M’

GAM

GBM

GAM
GAB
^
A I K
X
A H
^B GBM

Soient les deux points A et B de coordonnées connues avec la mesure des angles A et B. nous
devons déterminer les coordonnées du point M

1. Formules simples

Considérons le triangle AMB. En utilisant les relations des sinus, on a :

AM BM AB AM BM AB AB∗sinB
= = ou = = d ’ où AM =¿ et BM =¿
sinB sinA sinM sinB sinA sin ⁡(A + B) sin ⁡( A+ B)
X M = XA + X or X = AM*sin GAM
1
X M = XB + X or X = BM*sin GBM
2

et
Y M = YA + Y or Y = AM*cos GAM
1
Y M = YB + Y or Y = BM*cos GBM
2

X M = XA + AM*sin GAM
1
X M = XB + BM*sin GBM
2
et
Y M = YA + AM*cos GAM
1
Y M = YB + BM*cos GBM
2

Les deux valeurs permettent de faire un contrôle. Si ces valeurs concordent, le résultat final sera
la moyenne.
X M + XM et
Y M +Y M
XM = 1 2
YM = 1 2

2 2

9 /77
En pratique, les deux points connus en coordonnées constituent une condition nécessaire, mais
pas suffisante. Il faut obligatoirement stationner un troisième point.

2. Formule globale

En considérant le triangle AMB et en utilisant les relations des sinus, on a :

YM - YA = AM*cos GAM
Y M −Y A AM∗cos G AM Avec AM =¿
=
YB - YA = AB*cos GAB Y B−Y A AB∗cos G AB AB∗sinB
sin ⁡( A+ B)
Y M −Y A sinB∗cos G AM
=
Y B−Y A sin ⁡( A+ B)∗cos G AB

Y
M −¿Y A =(Y B −Y A )
sinB∗cos G AM
sin ( A +B )∗cos G AB
¿ et X M − X A=( Y M −Y A )∗tan G AM

3. Formule de Hatt

H est le projeté de M sur l’axe des X ; I est le projeté de B sur l’axe des X ; K est l’intersection de
(BM) et l’axe des X.

1er cas : AK = AH – KH (1)

2ème cas : AK = AI + IK (2)

Considérons le 1er cas : AK = AH – KH

 Expression de AH
AH
sin G AM =  AH = AM∗sin G AM =X M −X A (2)
AM
AH
tanG AM = AH = AM '∗tan G AM Or A M ' =Y M −Y A  AH =( Y M −Y A ) tan G AM (3)
AM '
(2)=(3)  AH =¿ X M − X A=( Y M −Y A ) tan G AM ( 4)

 Expression de KH
KH
tanG BM =  KH = A M '∗tanG BM =(Y ¿ ¿ M −Y A )∗tan GBM ¿ (5)
AM '
( 4 ) et ( 5 ) dans ( 1 ) donnent :

AK =( Y M −Y A ) tan G AM −¿ (Y ¿ ¿ M −Y A )∗tanG BM ¿

AK =( Y M −Y A )( tan G AM −tanG BM ) (6)

Considérons le 2ème cas : AK = AI + IK

10 /77
IK
tanG BM = IK =IB∗tan GBM avec IB=Y A−Y B
IB
 AI =X B −X A

AK =¿( X B− X A ¿−( Y B−Y A )∗tan G BM (7)

( 7 )=( 6 ) ( X B− X A ¿−( Y B−Y A )∗tan G BM =( Y M −Y A ) ( tan G AM −tan G BM )

( X B −X A )−( Y B−Y A )∗tan GBM

Y M −Y A =
tan G AM −tan G BM

X M − X A=( Y M −Y A ) tan G AM

Application

Soient deux points connus A et B tel que

XA = 782 333.32 XB = 785 489.74

m m

YA = 310 192.99 YB = 315 556.44

m m

L M =372.5104 gr L M =129.5077 gr

Calculer les coordonnées du point M par la méthode des formules, globale et de Hatt.

II. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR RELEVEMENT

1. Principe

Le relèvement est l’opération topométrique permettant de déterminer les coordonnées d’un

point en stationnant celui-ci pour viser au moins trois (3) points connus. Cependant en pratique,
trois (3) points constituent une condition nécessaire mais pas suffisante, alors il faut un 4ème
point connu en coordonnées rectangulaires.

2. Méthode géométrique ou relèvement italien

11 /77
Les points A, B et C sont des points connus. Le cercle passant par les points M, A et C coupe la
droite (BM) au point I. La détermination des coordonnées de I nous permet de calculer les
coordonnées de M en ramenant le problème à un simple calcul d’intersection.

Procédé de calcul

 Calcul du gisement AC

∆ X AC
G AC =arctg
∆ Y AC

 Calcul des gisements AI et CI :

G AI =G AC –

GCI =G AC ± 200+ α

 Détermination des coordonnées de I par la

formule de Hatt

( X A− X C )−( Y A−Y C )∗tanG AI

Y I −Y C =
tan G CI −tan G AI

X I − X C =( Y I −Y C ) tanG CI

 Calcul des gisements AM et CM.

Les points B, I et M étant alignés, on en déduit que : GBI =GBM

G AM =G BM −α
GCM =GBM +¿
 Calcul des coordonnées du point M.

( X C − X A ) −( Y C −Y A )∗tan GCM
Y M −Y A =
tan G AM −tan GCM

X M − X A=( Y M −Y A ) tan G AM

Actuellement, le relèvement se calcule par des méthodes numériques. La méthode de calcul

que nous venons de voir comporte une certaine ambigüité. En effet, si le point B est à l’intérieur
du cercle, il y aura confusion dans le calcul du gisement BI. Cependant, cette méthode est
toujours utilisée pour les calculs du relèvement double.

12 /77
 3. Méthode barycentrique

En visant les points connus, l’opérateur en station en M, fait un tour d’horizon.

Procédé de calcul

 Calculer les angles A, B, et C s’ils ne sont pas connus à partir des gisements.

A=G AC −G AB B=G BA−G BC

C=GCA −GCB

Contrôle : A + B + C = 200 gr

 Calculer les angles a, b et c à partir des lectures.

a=LC −LB b=L A− LC c=L B−L A

Contrôle : a + b + c = 400 gr

 Calculer les coefficients de chaque point.

1
mA=
cotan A−cotana

1
mB =
cotan B−cotanb

1
mC=
cotanC−cotanc

 Calculer les coordonnées du point M.

∑ mi . X i m A . X A + mB . X B +mC . X C
X M = i=1 n =
m A + mB +mC
∑ mi
i=1

∑ mi . Y imA . Y A + mB .Y B +mC .Y C
Y M = i=1n =
mA +mB +mC
∑
La formule est valable pour tous les m
cas de figure. Il suffit de tourner dans le sens
i
i=1
trigonométrique.

Les coordonnées approchées de M se calculent avec des combinaisons qui permettent une
meilleure détermination.

NB : On se contrôle en calculant le G0 de la station en considérant tous les points sur lesquels

l’on s’est relevé.
13 /77
 4. Méthode de Delambre
B (X ; Y)
A (X ; Y) Soient les gisements AM et BM.

D’après le schéma, GBM =G AM + α

α GCM =G AM + ¿

M
C (X ; Y)

Considérons le point M comme intersection à partir des points A et B avec des gisements pour
le moment inconnus G AM etGBM .

( X B −X A )−( Y B−Y A )∗tan GBM

Y M −Y A =
tan G AM −tan G BM

Et remarquant que :
1 cos G AM∗cos GBM cos G AM∗cos G BM
= =
tan G AM −tanG BM sin G AM∗cos G AM −sin GBM∗cos G AM sin ⁡(G AM −G BM )
Il vient que :

Y M −Y A =
[ ( X B−X A )−( Y B−Y A )∗tan G BM ]∗cos G AM∗cos GBM
−sin ⁡α

Avec sin ( G AM −G BM ) =sin (−α )=−sin α

Y M −Y A =
[ ( X B−X A )−
( Y B−Y A )∗sinG BM
cos G BM ]
∗cos G AM∗cos GBM

−sin ⁡α

( X B −X A ) cos GBM −( Y B−Y A )∗sin GBM

Y M −Y A = ∗cos G AM
−sin ⁡α

Opérant de la même manière en considérant M comme une intersection à partir des points A et
C avec des gisements pour le moment inconnus G AM etGCM =G AM + ¿

( X C − X A ) −( Y C −Y A )∗tan GCM
Y M −Y A =
tan G AM −tan G CM

Qui transformé comme précédemment devient :

14 /77
 Y M −Y A =
[ ( X C− X A ) −( Y C−Y A )∗tan GCM ]∗cos G AM∗cos GCM
−sin ⁡

Avec sin ( G AM −GCM )=sin ¿ ¿

Y M −Y A =
[ ( X C −X A ) −
( Y C −Y A )∗sin GCM
cos GCM ]
∗cos G AM∗cos GCM

−sin ⁡

( X C − X A ) cos GCM −( Y C −Y A )∗sin GCM

Y M −Y A = ∗cos G AM
−sin ⁡

On peut donc écrire :

( X B−X A ) cos GBM − ( Y B −Y A )∗sin GBM ( X C −X A ) cos GCM −( Y C −Y A )∗sinGCM

∗cos G AM = ∗cos G AM
−sin ⁡α −sin ⁡

( X B−X A ) cos ( G AM +α )−( Y B−Y A )∗sin ( G AM + α )

=( X C − X A ) cos ¿ ¿
sin ⁡α

( X B− X A ) cotg α∗cos G AM −( X B−X A ) sin G AM −( Y B −Y A ) cotg α∗sin G BM −( Y B −Y A ) cos G AM =¿

( X C −X A ) cotg∗cos G AM −( X C − X A ) sin G AM −( Y C −Y A ) cotg∗sin GCM −( Y C −Y A ) cos G AM

En regroupant tous les termes composés de sin G AM à part et tous les termes composés de
cos G AM à part, on obtient :

[ ( Y B−Y A ) cotg α −( Y C−Y A ) cotg−( X C−X B ) ] sin G AM =[ ( X B− X A ) cotg α −( X C−X A ) cotg+ ( Y C−Y B ) ] cos G AM

( X B− X A ) cotg α −( X C −X A ) cotg+ ( Y C −Y B )
tanG AM =
( Y B−Y A ) cotg α −( Y C −Y A ) cotg−( X C −X B )
Remarque : tanG AM =tan GMA or le gisement calculé est un gisement par relèvement d’où :

( X B −X A ) cotgα −( X C −X A ) cotg+ ( Y C −Y B )
tanG MA =
( Y B−Y A ) cotg α −( Y C −Y A ) cotg−( X C −X B )

III. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR DOUBLE

RELEVEMENT
B
A

α2 ❑2
❑1 α1
J
I M P ❑2 α1
❑1 15 /77
α1
Soit le point M dont on veut déterminer les coordonnées. De M, on ne peut voir que deux (2)
points connus A et D. Nous choisissons un point P situé à une distance donnée de M de manière
à pouvoir viser deux autres points B et C connus également.

La droite (MP) en la prolongeant va couper les deux (2) cercles l’un passant par les points A, D
et M au point I et l’autre passant par les points P, B et C au point J.

La solution du problème est donnée par la résolution du relèvement italien.

Procédé de calcul

a) Calculer les coordonnées de I

 Calculer gisement AD puis gisements AI et DI

∆ X AD
G AD=arctg
∆Y AD

G AI =G AD +❑1 et G DI =G DA−α 1

 Calculer les coordonnées de I par la formule de Hatt

( X A− X D )−( Y A −Y D )∗tan G AI
Y I −Y D=
tan G DI −tanG AI

X I − X D =( Y I −Y D ) tan G DI

b) Calculer les coordonnées de J

 Calculer gisement AD puis gisements AI et DI

∆ X BC
GBC =arctg
∆ Y BC

GBJ =GBC −❑2 et GCJ =GCB + α 2

 Calculer les coordonnées de I par la formule de Hatt

( X B− X C )−( Y B−Y C )∗tan GBJ

Y J −Y C =
tanGCJ −tan G BJ

16 /77
X J − X C =( Y J −Y C ) tan GCJ
 c) Calculer les coordonnées de
 Calculer gisement IJ et en déduire le gisement MP

Les points I, M, P et J sont alignés et connaissant les coordonnées de I et J, nous pouvons

calculer le gisement IJ

∆ X IJ
G IJ =arctg =GMP
∆ Y IJ

 Calculer les gisements AM et DM ; PB et PC

G MA=G MP +200+ α 1 G AM =G MA ±200

G MD=G MP +200−❑1 G DM =GMC ± 200

 Calculer les coordonnées de P et M par la formule de Hatt

G PB=GMP −α 2 et G PC =GMB +❑2

( X A − X D )−( Y A −Y D )∗tan G AM
Y M −Y D=
tan G DM −tan G AM

X M − X D =( Y M −Y D ) tan G DM

( X B− X C )−( Y B −Y C )∗tan GBP

Y P−Y C =
tan G CP −tanG BP

X P− X C =( Y P−Y C ) tan GCP

NB : Le contrôle se fera par la comparaison de la distance MP mesurée sur le terrain et le calcul

de la distance MP à partir des coordonnées trouvées.

17 /77
 IV. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR
RECOUPEMENT

M Le recoupement est le procédé topométrique

qui utilise simultanément l’intersection et le
relèvement par la détermination d’un point.

Soient A, B et C, 3 points connus. Seul le point B

C (inaccessible) est stationnable. On veut déterminer le point
inconnu M. on stationne B et on vise A, C, M et
B A
on fait les lectures LA, LC, LM. On stationne M, on
(stationnable) (inaccessible)
vise A et B et fait les lectures LA et LB.

Procédé de calcul

a) Calcul du G0moyen à partir des points connus

∆ X BA ∆ X BC
 GBA =arctg ; GBC =arctg
∆ Y BA ∆ Y BC

A C
A C G0 +G 0
D’où G 0 =GBA −L A ; G 0 =G BC −LC G 0moy =
2

b) Calcul des gisements BM et AM

GBM =G0 moy + L M ; G AM =G BM −^

AMB=GBM −(LB −L A )

c) Calcul des coordonnées de M par la formule de Hatt.

V. CALCUL DES COORDONNEES D’UN POINT DETERMINE PAR

TRILATERATION

M La trilatération est l’opération topométrique qui

consiste à mesurer uniquement la distance entre
D AM le point nouveau à déterminer et deux points

18 /77
DBM
A
connus (« dits anciens ». Ce procédé nécessitant la mesure de plusieurs distances s’appelle aussi
multilatération.

Procédé de calcul

1) Calcul du gisement AB et de la distance AB

2) Résolution du triangle ABM

3) en déduire les gisements AM et BM

4) Calcul des coordonnées du point M en utilisant l’une des méthodes de l’intersection.

Application.

19 /77
 CHAPITRE III : CALCULER ET COMPENSER UN CHEMINEMENT
PLANIMETRIQUE OU UNE POLYGONALE

I. DEFINITIONS ET PRINCIPE

Un cheminement planimétrique ou cheminement polygonal est une ligne brisée orientée dans
laquelle on connaît les longueurs des côtés successifs, et les angles que deux cotés consécutifs
font entre eux. Le but est de déterminer les coordonnées des sommets de la ligne brisée. Le
cheminement polygonal est aussi appelé par simplification polygonale. On distingue plusieurs
sortes de cheminements polygonaux :

1. Un cheminement en antenne ou ouvert

C’est un cheminement qui partant d’un point connu en coordonnées avec une référence
connue aboutit à un point à déterminer. Ce cheminement est à éviter parce qu’il n’y a pas de
contrôles.

2. Le cheminement encadré
C’est un cheminement qui partant d’un point connu en coordonnées avec une référence
connue passe par des points à déterminer pour se refermer sur un point également connu en
coordonnées avec une référence également connue. Autrement dit, c’est une ligne polygonale
qui relie deux points connus en coordonnées, avec une orientation (Gisement) au départ et à
l'arrivée.

3. Le cheminement dit fermé ou bouclé :

C’est un cheminement qui passe par des points à déterminer pour se
refermer sur le même point de départ. Autrement dit, c’est un
cheminement qui a le point et la référence de départ confondus avec
le point et la référence d'arrivée.

20 /77
 4. Le point nodal :

Un point nodal N est un point d'arrivée inconnu

commun à 3 ou plusieurs cheminements polygonaux.
N
Les coordonnées brutes du point nodal peuvent être
calculées depuis chaque cheminement polygonal. On
effectue alors la moyenne pondérée (en fonction de la
précision de chaque cheminement) de ces coordonnées
brutes pour obtenir les coordonnées définitives du
point nodal. Les coordonnées des sommets de chaque
cheminement polygonal peuvent ensuite être calculées comme précédemment (cheminement
encadré).

Remarque : On peut se trouver face à des configurations beaucoup plus complexes dans
lesquelles les cheminements sont imbriqués les uns dans les autres. Le calcul s'effectue alors
soit par la méthode des moindres carrés, soit plus simplement en fixant un cheminement
polygonal principal calculé en premier puis des cheminements successifs qui s'appuient sur les
cheminements de niveaux supérieurs. Il est à noter que plus on descend dans le niveau du
cheminement, moins il sera précis.

II. CALCULS

Les calculs s’effectuent en deux phases : orientation et coordonnées des différents sommets de
la polygonale. On détermine de proche en proche les coordonnées de tous les sommets de la
polygonale. Pour cela, il faut connaître :

- les coordonnées du point de départ.

- le gisement origine au départ du cheminement appelé gisement de départ.

- les distances horizontales des côtés réduites à la projection, si la précision l'exige.

- les angles que font entre eux les différents côtés consécutifs depuis la direction de référence
jusqu'à celle de l'arrivée.

Si le cheminement est encadré, il faut en plus connaître :

- les coordonnées du point d'arrivée.

21 /77
- le gisement origine à l'arrivée du cheminement appelé gisement d’arrivée.

1. Calcul d’un cheminement encadré

a) Orientation du cheminement

 Transmission des gisements

En général, pour la transmission des gisements (détermination de l'orientation de chaque côté

du cheminement polygonal), on utilise les angles de gauche en chaque station.

Ces angles correspondent à ceux situés à gauche du sens de parcours du cheminement. Chaque
angle de gauche est donc égal à la différence entre les lectures angulaires sur les points avant et
arrière.

α i=l S S −l S S
i i+1 i i−1

Le gisement du côté " Si S i+1" se déduit de celui du côté " Si S i−1".

On a ainsi : GS S =GS S ± 200+ α i

i i+1 i i−1

Connaissant le gisement de référence au point de départ Gdépart , on peut donc calculer les
gisements de chaque côté du cheminement de proche en proche jusqu'au gisement de
référence au point d'arrivée Garrivée connu

 Ecart de fermeture angulaire ou fermeture angulaire

22 /77
On appelle fermeture angulaire notée fa la différence entre le gisement arrivée observé
(gisement calculé à partir des mesures entachées des imprécisions de mesure) et le gisement
arrivée connu ou calculé à partir à partir des coordonnées connues.

f a=Garrivée observé
−Garrivée connu

 Ecart type sur l’écart de fermeture angulaire

 Tolérance de fermeture angulaire

Les tolérances sont calculées à partir des précisions des appareils (écarts types sur les mesures).
On considère que sur chaque mesure angulaire on a un écart type σ α et sur les gisements de
départ et d'arrivée un écart type σ G.

L'écart type sur l'écart de fermeture angulaire vaut alors:

σ f =± √ 2 ×σ G2+(n+ 1)×σ α 2
a

On considérera que la précision sur les gisements de départ et d'arrivée est nulle. D’où :

L'écart type sur l'écart de fermeture angulaire :

σ f =± σ α √ n+1
a

et la tolérance sur fa :

8
T a=± σ f
3 a

8
T a=± σ √ n+1
3 α

 Compensation de l'écart de fermeture angulaire

En fonction des gammes de précision des appareils de mesure, des conditions dans lesquelles
ont été effectuées ces mesures, du cahier des charges,... on compare l'écart de fermeture
angulaire à la tolérance.

23 /77
Si fa n'est pas acceptable, soit on effectue une recherche de faute, soit on refait les mesures.

Si fa est acceptable, on répartit cet écart sur tous les angles mesurés et donc les gisements. La
compensation sur chaque angle vaut:

−f a
c a=
n+1

La compensation sur chaque gisement brut vaut :

−i× f a
ca =
i
n+1

b) Coordonnées rectangulaires des sommets

 Coordonnées brutes

Tous les côtés du cheminement polygonal sont dorénavant orientés puisque nous avons calculé
et compensé leurs gisements. Nous pouvons calculer les coordonnées planes des sommets
(stations) de proche en proche, à partir du point de départ, comme des points rayonnés. Nous
allons obtenir des coordonnées brutes des sommets et du point d'arrivée.

 Fermeture planimétrique ou écart de fermeture planimétrique

Les coordonnées du point d'arrivée B étant connues, nous pouvons donc les comparer avec
celles mesurées et calculées ci-dessus (coordonnées brutes).

Nous obtenons : - un écart de fermeture en X: f X= X B observé− X B calculé ou connu

- un écart de fermeture en Y: f Y =Y B observé −Y B calculé ou connu

Ces deux écarts forment un vecteur de fermeture planimétrique noté fp dont la norme vaut :

f p=√ f 2X + f 2Y

24 /77
  Ecart type sur l’écart de fermeture planimétrique ou linéaire.

L’étude des erreurs de fermeture planimétrique se fait sur un cheminement considéré

rectiligne. Nous distinguons deux erreurs de fermeture planimétrique qui sont :

 Erreur longitudinale ou écart de fermeture longitudinale σ D ou σ L .

Elle est provoquée par l’erreur sur les distances.

 Erreur transversale ou écart de fermeture transversale.

Elle est provoquée par l’erreur sur l’angle. Supposons que l’erreur commises sur les sommets
sont les mêmes.

 Tolérance de fermeture planimétrique ou tolérance planimétrique

8
Tp=± σ f
3 p

25 /77
 8 2 2
Tp=± √ σ +σ
3 T L

 Compensation de l'écart de fermeture planimétrique

Comme pour la fermeture angulaire, on compare l'écart de fermeture planimétrique à la

tolérance.

Si fp>Tp, fp n'est pas acceptable. Soit on effectue une recherche de faute, soit on refait les
mesures.

Si fp<Tp, fp est acceptable. On répartit les écarts de fermeture en x et en y sur les coordonnées
brutes.

Il existe plusieurs méthodes pour répartir ces écarts (compensation proportionnelle, par
rotation homothétie, à la cubique, par les moindres carrés, méthode Carrier, …).

Pour simplifier, on va utiliser la méthode de compensation proportionnelle à la longueur des

côtés.

Les compensations de chaque côté en X et en Y valent :

−Di
cX =
i n
×f X
∑ Di
i

−Di
cY =
i n
×fY
∑ Di
i

Remarque : Si les côtés sont de longueurs sensiblement égales, tous les côtés vont être
compensés de la même valeur.

2. Calcul d’un cheminement bouclé ou fermé

Deux cas se présentent : le cas où un côté du cheminement bouclé est orienté c’est-à-dire que
le gisement d’un côté est connu et le cas où la référence est un point qui ne fait pas parti du
cheminement bouclé. B

1er cas :

A 3

26 /77
Transmission de gisement

Elle se fait comme dans le cas du cheminement encadré. Le gisement de départ sera G AB .

NB : Les formules ci-après sont valables dans les deux cas du cheminement bouclé sauf au
niveau de la compensation angulaire.

 Fermeture angulaire :
n
fa = ∑ α i− ( n−2 ) ×200 (angles intérieurs)
i=1

n
= ∑Tolérance
fa  α i− ( n+2 ) × 200 (angles extérieurs)
angulaire :
8
T a=± σ √ n
i=1
3 α

 Compensation angulaire
−f a
c a= sur chaque angle
n
Si fa < Ta
−i× f a
ca = sur chaque gisement
i
n

 Coordonnées rectangulaires des sommets

 Fermeture planimétrique :

f X= ∑ ∆ X i
F p= √ f X 2+ f Y 2
f Y =∑ ∆ Y i

 Erreur longitudinale : σ L =± σ D √ n

 Erreur transversale : σ T =± L . σ α (rad )

√ n
3

 Tolérance planimétrique : T p=± √ σ L2 +σ T 2

27 /77
 EXERCICE D’APPLICATION
E
85.61 m 85.23 m

142.740 gr
D
A
115.701 gr
115.735 gr

79.39 m 107.157 gr 90.67 m

118.657 gr
+

B C
120.09 m

Calculer les coordonnées compensées des sommets de la polygonale ci-après, sachant que les
coordonnées du point A sont X=1 996.50 m et Y= 911.77 m, le gisement EA est 264.3633 gr ;
l’écart type angulaire sur une visée est ±3 mgr et l’écart type sur la distance est ±5 cm.

28 /77
 CORRECTION

Point
α i(gr) C α (gr ) α i (gr ) Gi (gr ) Di (m) ∆ X i (m) C X (m) ∆ X i (m) C Y (m) X (m) Y (m)
s i comp comp i i

E
264,3633
A 115,735 0,002 115,737 1996,50 911,77
180,1003 79,39 24,414 0,001 -75,543 -0,001
B 118,657 0,002 118,659 2020,92 836,23
98,7593 120,09 120,067 0,003 2,340 -0,003
C 107,157 0,002 107,159 2140,98 838,56
5,9183 90,67 8,417 0,001 90,278 -0,002
D 115,701 0,002 115,703 2149,40 928,84
321,6213 85,23 -80,362 0,001 28,393 -0,002
E 142,740 0,002 142,742 2069,04 957,23
264,3633 85,61 -72,544 0,001 -45,459 -0,002
A 1996,50 911,77
599,990 0,010 600,000 460,99 -0,007 0,007 0,010 -0,010

f a=599.990−( 5−2 ) ×200=−0.010 gr

8 f X=−0.007 m; f Y =−0.010 m
Ta= × 0.003 √2 × √5=±0.025 gr
3
F P= √ (−0.007 ) + (−0.010 ) =+ 0.112m
2 2
f a< Ta Compensation

−(−0.010) σ L =0.05 √ 5=± 0.112 m

Cα = =+0.002 gr
5
√
i

2 ×5
σ T =30 × √ 2 ×1.57 . 10−6 × 460.99 × =±0.056 m
3

8
T P=
3
√ 0.1122 +0.0382=±0.334 m
−f X . Di
CX = et
i
∑ Di 29 /77

−f Y . Di
C =
F P <¿ T P Compensation

30 /77
2ème cas :

B +
Conseil :

Sur le terrain, il est souhaitable de faire un tour d’horizon en A

B afin d’éviter des erreurs d’orientation de tout le
cheminement.
A C
L’angle extérieur A est la somme des deux angles extérieurs
RAB et ^
^ FAR .
F D

A=( l B −l R ) +(l R −l F ).
^

Pour le calcul du cheminement, ces deux angles apparaissent dans le tableau. Le premier angle
sera ^
RAB le second ^ ABC … et le dernier ^
FAR .

ATTENTION :

a. le gisement de départ sera et non G AR et G RA.

b. le 1er gisement calculé sera alors : G AB=G AR + ^
RAB=G AR + ( l A −l A B R
)
c. la transmission des gisements se fera selon la formule : Gi=Gi−1 +α g ± 200 gr

A partir du côté BC jusqu’à FA (dernier côté), tous les gisements se calculent selon cette
formule.
d. Dans la compensation angulaire, le nombre de nouveaux côtés correspond au nombre de
−f a
sommets : soit n. c a = Or dans le tableau, vous aurez (n+1) angles.
i
n
En effet, l’angle ^A est cindé en 2 angles : ^
RAB et ^
FAR . Ainsi, vous diviserez par 2 la
correction à apporter à l’angle ^
A . Soit

cα cα
^
RAB corrigé = ^
RAB brut + ^
FAR corrigé = ^
FAR brut +
2 2

EXERCICE D’APPLICATION

Calculer les coordonnées compensées des sommets de la polygonale ci-après, sachant que les
coordonnées du point A sont X=1 000.00 m et Y= 1000.00 m, le gisement AR est 350.0000 gr ;
l’écart type angulaire sur une visée est ±3 mgr et l’écart type sur la distance est ±5 cm.

R +
276,335 gr
85,453 gr 176,00 m
194,50 m B

199,792 gr A C 276,426 gr
31 /77
210,10 m
 Point
α i(gr) C α (gr ) α i (gr ) Gi (gr ) Di (m) ∆ X i (m) C X (m) ∆ X i (m) C Y (m) X (m) Y (m)
s i comp comp i i

R
350,0000
A 85,453 +0,002 85,455 1000,00 1000,00
35,4550 194,50 102,808 +0,007 165,108 +0,002
B 276,335 +0,003 276,338 1102,82 1165,11
111,7930 176,00 172,989 +0,006 -32,417 +0,001
C 276,426 +0,003 276,429 1275,81 1132,69
188,2220 208,50 38,355 +0,008 -204,942 +0,002
D 284,204 +0,003 284,207 1314,17 927,75
272,4290 183,00 -166,105 +0,007 -76,800 +0,001
E 277,773 +0,004 277,777 1148,08 850,96
350,2060 210,10 -148,082 +0,007 149,043 +0,002
A 199,792 +0,002 199,794 1000,00 1000,00
350,0000
R
1399,983 + 0.017 1400,000 972,100 -0,034 +0,034 -0,008 +0,008

f X=−0.034 m; f Y =−0.008 m
f a=1399.983−( 5+2 ) ×200=−0.017 gr
F P= √ (−0.034 ) + (−0.008 ) =+0.035 m
2 2

8
Ta= × 0.003 √2 × √5=±0.025 gr σ L =0.05 √ 5=± 0.112m
3

f a< Ta Compensation

−(−0.017)
σ T =0.003 √ 2 ×
π
200
5
√
× 972.10 × =± 0.084 m
3
Cα = =+0.003 gr 8
i
5 T P=
3
√ 2 2
0.112 +0.075 =±0.373 m

32 /77
 −f X . Di
F P <¿ T P Compensation CX = et
i
∑ Di
−f Y . D i
CY =
i
∑ Di

33 /77
 CHAPITRE IV : CALCUL DE SURFACES

INTRODUCTION

Le calcul de surface est l’une des opérations la plus importante en topographie.

Ces opérations vont servir à la vente des parcelles, à l’évaluation de l’impôt foncier et au
partage des terrains.

Les calculs de surfaces sont effectués soit d’après les relevés de terrain, soit d’après un plan
déjà établi. Les surfaces sont déterminées à partir de :

- Formules mathématiques usuelles ;

- Procédés topométriques (numériques, graphique…) ou
- par comparaison de différentes méthodes.

Ce présent chapitre fera l’objet des calculs de surface à partir de formules mathématiques
usuelles et procédés topométriques numériques.

I. CALCUL DE SURFACES PAR DECOMPOSITION EN FIGURES

GEOMETRIQUES SIMPLES

Soit un polygone dont on veut déterminer la superficie. Pour se faire, on va procéder à la

décomposition en figures géométriques simples qui sont : le triangle, le rectangle, le carré et le
trapèze.

Rappel
1 1 a² 1
Triangle : S= ab sin c= =√ p ( p−a )( p−b ) ( p−c )= b× h
2 2 ( cot B+ cotC ) 2

Rectangle : S= L×l

Carré : S=a ²

Trapèze : B+b
S= ×h
2

34 /77
 II. Y CALCUL DE SURFACES PAR COORDONNEES RECTANGULAIRES

Y1 1
Soient les trapèzes (1, X1, X3 ; 3) ; (1, 2, X2, X1) et (2, X2,
3 X3 ; 3)
Y3
Y2 2 Trapèze (1, X1, X3, 3)
2 S 1=( Y 1 +Y 3 ) ( X 1−X 3 )

X Trapèze (1, 2, X2, X1) 2 S 2=( Y 1 +Y 2 ) ( X 2−X 1 )

X3 X1 X2
Trapèze (2, X2, X3 ; 3) 2 S 3=( Y 3 +Y 2 ) ( X 2−X 3 )

2 S 1 ,2 ,3 =2 S1 +2 S 2−2 S3

2 S=( Y 1+Y 3 ) ( X 1−X 3 ) + ( Y 1+Y 2) ( X 2−X 1 ) −( Y 3 +Y 2 ) ( X 2− X 3 )

2 S=Y 1 X 1−Y 1 X 3 +Y 3 X 1−Y 3 X 3 +Y 1 X 2−Y 1 Y 1 +Y 2 X 2−Y 2 X 1−Y 3 X 2 +Y 3 X 3−Y 2 X 2+Y 2 X 3

2 S=Y 1 X 2−Y 1 X 3 −Y 2 X 1+Y 2 X 3 +Y 3 X 1−Y 3 X 2

2 S=Y 1 ( X ¿ ¿ 2−X 3)+Y 2 (X ¿ ¿3−X 1 )+Y 3 (X 1− X 2)¿ ¿

−2 S= X 1 (Y ¿ ¿ 2−Y 3 )+ X 2 (Y ¿ ¿ 3−Y 1)+ X 3 (Y 1−Y 2 )¿ ¿

Formule générale : n n
2 S=∑ Y i ( X ¿ ¿ i+ 1−X i−1 )=−¿ ∑ X i (Y ¿ ¿ i+1−Y i−1)¿ ¿ ¿
i=1 i=1

Application

Pts X (m) Y (m) X i +1−X i−1 Y i (X ¿ ¿ i+1−X i−1)¿ Y i+ 1−Y i−1 X i (Y ¿ ¿ i+1−Y i−1)¿
A 158.183 0.000
B 100.000 0.000
C 91.904 50.829
D 162.188 47.790
E 197.762 31.970
A 158.183 0.000
B 100.000 0.000
∑ ..................... ...................... .............................. ............................
S (m²) .............................. ............................

35 /77
 O
III. L1CALCUL DE SURFACES PAR COORDONNEES POLAIRES (ANGLES,
1
L2 DISTANCES)
d1
L3 2 S P 12=d1 . d 2 sin (¿ L2−L1)¿ 2 S P 23=d 2 . d 3 sin(¿ L3−L2 )¿
d2
2
P
d3 2 S P 13=d 1 . d 3 sin(¿ L3−L1)¿ 2 S 123 =2 S P 12+ 2 S P 23−2 S P 13
3
2 S=d1 . d 2 sin (¿ L2−L1)+d 2 . d 3 sin( ¿ L3 −L2)−d 1 . d 3 sin(¿ L3−L1)¿ ¿ ¿

2 S=d1 . d 2 sin (¿ L2−L1)+d 2 . d 3 sin(¿ L3 −L2)+d 1 . d 3 sin(¿ L1 −L3)¿ ¿ ¿

Formule générale :

n
2 S=∑ d i . d i+1 sin(¿ Li+ 1−Li) ¿
i=1

La formule est valable quel que soit la parcelle et quel que soit la position de la station, à
l’intérieur comme à l’extérieur du polygone.

Inconvénient : Il est très rare de voir d’une seule station tous les sommets de la parcelle.

Le rayonnement n’offre aucun contrôle. Si les distances sont longues, le chaînage perd sa
précision.

Avantages : rapidité des éditions, utilisation de formule simples.

NB : Il est déconseillé d’utiliser cette formule pour le foncier et le cadastre (parce qu’il n’y a pas
de contrôle) Sommet Gi (gr) Di (m)
s
Application 1 333.842 72.41
5 6
2 368.805 40.67
1 2 3 340.791 25.15
4 59.999 22.83
3 4 5 29.467 45.36
P 6 68.666 75.19
7 121.503 68.51
8 7 8 277.014 67.02

36 /77
Résolution

Station Sommets Di (m) Gi (gr) Di+ 1 . Di . sin(G¿ ¿i+1−Gi)¿

1 72.41 333.842
2 40.67 368.805

3 25.15 340.791

4 22.83 59.999

5 45.36 29.467
P
6 75.19 68.666

7 68.51 121.503

8 67.02 277.014

1 72.41 333.842
2 S (m²)
S (m²)

37 /77
 IV. CALCUL DE SURFACES
^ PAR LA METHODE POLYGONALE (OU METHODE
ab
DE SARRON)
3
b ^ a^c
4’’
bc
^
ab ^ La méthode de Sarron permet de calculer la
ab
surface de tout polygone de n côtés dont on
2 connaît les longueurs de (n-1) côtés et les
2’ h c valeurs des (n-2) angles.
b
4’
a ha Soit la parcelle 1, 2, 3,4 parcourue dans le sens
des gisements. On a les angles dirigés qui sont
3’
donnés par les vecteurs des côtés. Un angle
4
1 dirigé est un angle qui est le complément de
l’angle interne compris par les côtés.

' '
h a=4 3 +3 ' 4
' ' ' ^
4 3 =22 =b sin ab

3' 4=c sin ( 200− a^c ) =c sin a^c

^ c sin a^c
h a=b sin ab+

^ ) =c sin bc
h b=c sin ( 200−bc ^

^ c sin a^c ) +bc sin bc

2 S 123 =2 S 124 +2 S 234 =a . h a+ b . hb=a ( b sin ab+ ^

^
2 S=ab sin ab+ac ^
sin a^c +bc sin bc

Remarque :

^ bc
a^c= ab+ ^; ^
ab=200− 2^ ; ^
bc=200− 3^

Pour un polygone de n côtés, le nombre de termes dans la formule de Sarron est : ( n−1 ) ( n−2 )
2
Le nombre de côtés utilisés est n-1 et celui des angles dirigés n-2.

Formule générale
^
2 S=ab sin ab+ac sin a^c +…+ a(n−1)sin ^
(a , n−1) n-2 termes

^ + …+¿ b(n−2)sin ^
^ bd sin bd
+bc sin bc+ (b , n−2)¿ n-2 termes

.
.
.
n(n-1)= 1 terme
+ ( n−2 ) (n−1)sin ^
( n−2 ) (n−1)
38 /77
Application :

^
bc D
^B=110.0506 AB=98.18 m
c ^
cd
C E ^
C=87.1931 gr BC=51.47 m

b ^
D=176.0324 gr gr CD=70.35 m
^
ab
DE=38.96 m
a A
B

^
2 S=ab sin ab+ac ^
sin a^c +ad sin ad
^ bd sin bd
+bc sin bc+ ^

^
+cd sin cd

ij^
Pts ij^ Di ^ 2Si
[Link]
^ a^c ^
ad
A ab
a
… … …^ … … … …...
ac. sin a^c ^
ad sin ad
ab sin ab
B ^
ab
b ^
bc ^
bd
^… … … …...
C … …bc
… ^
[Link] bc ^
bd sin bd
c
D ^
cd ^
cd
d
…E… … … … … … . ^
cd sin cd
2 S (m²) = ∑ 2 Si
S (m²)

39 /77
 V. C
REDRESSEMENT DE LIMITES
^
bc
b
A E
c x
Le ^ problème consiste à rendre
a ab c^x
B rectiligne une limite formée par
D
une ligne brisée (limite entre 2
parcelles, rectification de
chemin ....) de manière que les
parcelles considérées
conservent leur superficie respectives.

La solution la plus simple est de remplacer la limite considérée (ABCD par exemple) par une
droite issue de A et aboutissant en un point E.

Les surfaces hachurées de part et d’autre de AE doivent s’équilibrer c.-à-d. que la surface du
polygone ABCDE doit être nulle.

On choisit donc comme sens de parcours celui des gisements et on applique la méthode de
Sarron.

^
2 S=ab sin ab+ac sin a^c +ax sin a^x =0

^ bx sin bx
^ x=−¿ ¿
+bc sin bc+

+cx sin c^x

Application

^
ab
D ^ ^
a dx a = 14.05 m ab=147.39 gr
b d x ^
J b = 47.45 m bc=376.97 gr
I
^
bc c c = 40.70 m ^ =342.41 gr
cd
B
C d = 84.95 m ^
dx=142.45 gr
^
cd

^
2 S=ab sin ab+ac sin a^c +ax sin a^x+ ax sin a^x =0

^ bd sin bd
+bc sin bc+ ^ +¿ bx sin bx
^¿

^ + cx sin c^x
+cd sin cd
^
+dx sin dx

40 /77
 ^ ac sin a^c +bc sin bc+
−(ab sin ab+ ^ cd sin cd)
^
x=
^ +c sin c^x+ d sin dx
a sin a^x+ b sin bx ^

^
ij
Pts ^
ij Di 2Si
^
[Link]
^ a^c ^
ad a^x
I ab
a
…ab…
sin…+
^ … … ..ac
ab x sin a^c ^
ad sin ad (a sin a^x) x
A ^
ab
b ^
bc ^ ^
bx
bd
^
B… .. …+…
bc . … .bc ^
x sin bc ^
bd sin bd ^) x
( b sin bx
c
C ^ ^
cd c^x
cd
… … …+ … …dx ^
cd sin cd ( c sin c^x ) x
D ^
dx
^
dx
E x
… … …+ … … x ^)x
( d sin dx
∑
…… ……
ij sin ij^ … … … .

∑…( i …
… i^j )…
sin… x …x

− …………………………………… …
x= = =… … … … … … … … . .
… … … … … … … … … … … … … … ..

41 /77
 CHAPITRE VI : COMPENSATION PAR LA METHODE DES MOINDRES
CARRES
INTRODUCTION
Dans toute science d’observation, il est d’usage de multiplier les observations. Ainsi, en sciences
géographiques, nous effectuons toujours les mesures en nombre supérieur à celui des besoins de
la résolution mathématique du problème à résoudre. Cette manière de faire permet tout d’abord
de déceler et d’éliminer les fautes mais surtout l’expérience montre que cet usage conduit à une
amélioration des résultats. Cependant, l’observateur qui connaît ou détermine la précision des
mesures qu’il a faites est amené alors à répondre aux questions suivantes :

- à partir de l’ensemble des mesures effectuées, quelle valeur faut-il choisir pour la quantité
qu’il s’agit de déterminer ?
- quelle est la précision de la valeur adoptée ?

Malgré la redondance des observations dont il dispose, le calculateur souhaite en outre obtenir
une valeur unique pour la quantité qu’il se propose de déterminer. Cette exigence va le conduire
à modifier (compenser) les diverses observations.

Pour cela, il va substituer aux éléments observés des éléments conventionnels calculés formant
un système tel que :

- soit assurée l’homogénéité de l’ensemble tel que le système des valeurs adoptées soit
unique et satisfasse à toutes les conditions requises par la structure étudiée ;
- les valeurs substituées aux observations réelles soient voisines que possible dans leur
ensemble des valeurs observées et que les écarts entre valeurs observées et valeurs
compensées soient de l’ordre de grandeur des erreurs d’observations.

Un tel système forme un ensemble qui est dit ‘’compensé ‘’. Les compensations peuvent revêtir
différentes formes :

 empiriques (par expérience)

 mécaniques ou analogiques
 graphiques
 numériques.

Les formes numériques sont les plus rigoureuses. Elles peuvent entraîner de grandes quantités de
calculs qui sont aisément aujourd’hui abordés grâce aux traitements sur ordinateurs.

En sciences géographiques, la méthode numérique de compensation dérive des moindres carrés ;

ce qui n’est pas indispensable mais commode et admis. Il s’agit en définitive d’une méthode
simple, bien étudiée et connue dont on s’entend à dire qu’elle donne de bons résultats. On ne
voit pas actuellement quelle autre méthode pourrait valablement lui être opposée.

42 /77
 I. OBSERVATIONS INDIRECTES : METHODES DES MOINDRES CARRES

Hypothèse

Nous supposons que les observations sont supérieures en nombre d’inconnues. Toutes les
relations d’observations doivent être mises sous forme linéaire.

Pour la démonstration, nous allons utiliser des relations d’observations qui se limitent à trois
inconnues.

a i X + bi Y + c i Z+ k i=0

Les coefficientsa i, b i et c i sont à déterminer ; X, Y et Z sont des inconnues ; k i sont des valeurs
observées directement sur le terrain, elles sont donc entachées d’erreurs et la relation
d’observation est alors :

a i X + bi Y + c i Z+ k i=v i avec v i=résidu

a 1 X + b1 Y + c1 Z + k 1=v 1

a 2 X + b2 Y + c 2 Z+ k 2=v 2

a 3 X + b3 Y + c 3 Z+ k 3=v 3
. .
. .
. .

a n X + bn Y +c n Z +k n =v n

n
Nous allons rendre minimum la somme des carrés des résidus : ∑ v i =minimum
2

i=1

Pour rendre minimum la somme des carrés des résidus, il faut que la dérivée partielle par rapport
à chaque variable soit nulle.

¿ → n
2∑ ai (a ¿ ¿i X +b i Y + ci Z + k i)=0 ¿
i=1

n
2∑ bi (a ¿ ¿i X +b i Y + ci Z + k i)=0 ¿
i=1

Le développement donne :

∑ ai ai X +∑ a i b i Y +∑ ai ci Z +∑ ai k i=¿0
∑ bi ai X +∑ bi bi Y +∑ b i ci Z +∑ b i k i=¿0
∑ ci ai X +∑ c i b i Y +∑ ci c i Z+∑ c i k i=¿ 0

43 /77
 Avec

∑ ai ai=[ aa ] ∑ bi bi=[ bb ] ∑ ci c i=[ cc ]

∑ ai bi=∑ bi ai=[ ab ] ∑ bi ci =∑ c i b i=[ bc ] ∑ ci k i =[ ck ]
∑ ai ci =∑ c i ai=[ ac ] ∑ bi k i=[ bk ]
∑ ai k i=[ ak ]
Le système devient donc : [ aa ] X + [ ab ] Y + [ ac ] Z + [ ak ] =0
[ ba ] X + [ bb ] Y + [ bc ] Z+ [ bk ] =0
[ ca ] X + [ cb ] Y + [ cc ] Z + [ ck ] =0

On peut aussi écrire ce système sous forme matricielle :

[ ][ ][ ]
[ aa ] [ ab ] [ ac ] X [ ak ]
[ ba ] [ bb ] [ bc ] × Y + [ bk ] =0
[ ca ] [ cb ] [ cc ] Z [ ck ]

Ce système est appelé système des équations finales. Il a la particularité d’être un système
symétrique par rapport à la première diagonale et compte autant d’équations que d’inconnues.

La résolution de ce système donne une solution unique quel que soit le chemin suivi.

1. Formation pratique des équations finales

A partir des relations d’observations, nous créons le tableau des relations d’observations.

Tableau des relations d’observations

N° X Y Z ki ∑ ou S
1 a1 b1 c1 k1 s1 = a1 + b1+ c1+ k1
2 a2 b2 c2 k2 s2 = a2 + b2+ c2+ k2
. . . . . .
. . . . . .
. . . . . .
n an bn cn kn sn = an + bn+ cn+ kn

A l’aide du tableau des relations d’observations, nous formons les différentes équations finales :

44 /77
 a) Etape de la formation de la 1ère équation

Multiplicateur X Y Z ki ∑
a1 *.... a1 b1 c1 k1 s1
a2*.... a2 b2 c2 k2 s2
. . . . . .
. . . . . .
. . . . . .
an*.... an bn cn kn sn

On obtient la 1ère équation finale :

X Y Z ki ∑
a1a1 a1b1 a1c1 a1k1 a1s1
a2a2 a2b2 a2c2 a2k2 a2s2
. . . . .
. . . . .
. . . . .
anan anbn ancn ankn ansn
[ aa ] [ ab ] [ ac ] [ ak ] [ as ]
Remarque:
[ as ] = [ aa ]+[ ab ]+[ ac ] +[ ak ]. Cela sert de contrôle des résultats

b) Etape de la formation de la 2ème équation :

Multiplicateur Y Z ki ∑
b1*.... b1 c1 k1 s1
b2*.... b2 c2 k2 s2
. . . . .
. . . . .
. . . . .
bn*..... bn cn kn sn

On obtient la 2ème équation finale :

Y Z ki ∑
b1b1 b1c1 b1k1 b1s1
b2b2 b2c2 b2k2 b2s2
. . . .
. . . .
. . . .
bnbn bncn bnkn bnsn
[ bb ] [ bc ] [ bk ] [ bs ]
Remarque:
[ bs ] = [ ab ]+[ bb ]+[ bc ]+[ bk ] . Cela sert de contrôle des résultats

45 /77
 c) Etape de la formation de la 3ème équation

Multiplicateur Z ki ∑
c1*.... c1 k1 s1
c2*.... c2 k2 s2
. . . .
. . . .
. . . .
cn*..... cn kn sn

On obtient la 3ème équation finale :

Z ki ∑
c 1c 1 c 1k 1 c 1s 1
c 2c 2 c 2k 2 c 2s 2
. . .
. . .
. . .
cncn cnkn cnsn
[ cc ] [ ck ] [ cs ]
Remarque:
[ cs ] = [ ac ] +[ bc ]+[ cc ] +[ ck ]. Cela sert de contrôle des résultats

Tous ces différents résultats nous permettent d’établir la table des équations finales
suivante :
Table des équations finales :

X Y Z k ∑
[ aa ] [ ab ] [ ac ] [ ak ] [ as ]= [ aa ]+[ ab ]+[ ac ] +[ ak ]
[ ab ] [ bb ] [ bc ] [ bk ] [ bs ] = [ ab ]+[ bb ]+[ bc ]+[ bk ]
[ ac ] [ bc ] [ cc ] [ ck ] [ cs ]= [ ac ] +[ bc ]+[ cc ] +[ ck ]

Remarque :
A partir du tableau des relations d’observations, on peut directement calculer les
différents coefficients et passer à la table des équations finales.
Les coefficients ai, bi, ci et ki doivent être à 4 chiffres après la virgule.
Les valeurs à prendre en compte dans la colonne des sommes sont les résultats[ as ], [ bs ] et
[ cs ].

46 /77
 Application
Soit le système des relations d’observations ci-dessous :
1.0596 X + 0.7072Y + Z+ 0.83=v 1
0.6124 X−2.7157 Y + Z +0.80=v 2
−2.4554 X +0.7317 Y + Z+ 0.72=v 3
−0.2054 X−2.0941 Y + Z−1.54=v 4
0.6736 X + 1.6504 Y +Z — 0.80=v 5
Etablir le tableau des relations d’observations et la table des équations finale

Correction

Tableau des relations d’observations

N° X Y Z ki ∑
1 1.0596 0.7072 1 0.83
2 0.6124 −2.7157 1 0.80
3 −2.4554 0.7317 1 0.72
−0.2054 −2.0941 −1.54
4 1
0.6736 1.6504 — 0.80
5 1

a) 1ère équation finale

Multiplicateur X Y Z ki ∑
1.0596 1.0596 0.7072 1 0.83
0.6124 0.6124 −2.7157 1 0.80
−2.4554 −2.4554 0.7317 1 0.72
−0.2054 −0.2054 −2.0941 −1.54
1
0.6736 0.6736 1.6504 — 0.80
1

On obtient la 1ère équation finale :

X Y Z ki ∑

Contrôle

47 /77
 b) 2ème équation :

Multiplicateur Y Z ki ∑
0.7072 0.7072 1 0.83
−2.7157 −2.7157 1 0.80
0.7317 0.7317 1 0.72
−2.0941 −2.0941 −1.54
1
1.6504 1.6504 — 0.80
1

On obtient la 2ème équation finale :

Y Z ki ∑

Contrôle :

c) 3ème équation

Multiplicateur Z ki ∑
1 1 0.83
1 1 0.80
1 1 0.72
−1.54
1 1
— 0.80
1 1

On obtient la 3ème équation finale :

Z ki ∑

Contrôle

48 /77
 Table des équations finales :

X Y Z k ∑
8.0227 -1.1685 -0.3152 -0.6211 5.9179
15.5196 -1.7205 0.8458 13.4764
5.0000 0.0100 2.9743

2. Résolution du système des équations finales

L1 [ aa ] X + [ ab ] Y + [ ac ] Z + [ ak ] =0

L2 [ ab ] X + [ bb ] Y + [ bc ] Z+ [ bk ] =0

L3 [ ac ] X + [ cb ] Y + [ cc ] Z+ [ ck ] =0

Ce système peut être résolu par la méthode de substitution ou par la méthode des matrices.

Nous allons utiliser la méthode de substitution.

1ère étape : Nous allons exprimer l’inconnue X en fonction des inconnues Y et Z. pour cela, nous
allons considérer la ligne L1.

− [ ab ] [ ac ] [ ak ]
X= Y− Z−
[ aa ] [ aa ] [ aa ]

Résolvante de X.

Ensuite, nous allons reporter X dans les lignes L2 et L3.

[ ab ] [−[ ab ]
[ aa ]
Y−
[ ac ]
[ aa ]
Z−
[ ak ]
[ aa ] ]
+ [ bb ] Y + [ bc ] Z + [ bk ] =0

][
]]
− ab ]
[ [ ac ] [ ak ]
[ ac Y− Z− + bc ] Y + [ cc ] Z+ [ ck ] =0
[
[
aa ] [ aa ] [ aa

[[ bb ] −[ ab ]
] [
[ ab ]
[ aa ]
Y + [ bc ]− [ ab ]
] [
[ ac ]
[ aa ]
Z + [ bk ]− [ ab ]
[ ak ]
[ aa ]
=0
]
[ bb 1 ] [ bc 1 ] [ bk 1 ]

[[ bc ]− [ ac ]
[ ab ]
[ aa ] ] [
Y + [ cc ]− [ ac ]
[ ac ]
[ aa ] ] [
Z+ [ ck ] −[ ac ]
[ ak ]
[ aa ]
=0
]
[ bc 1 ] [ cc 1 ] [ ck 1 ]

49 /77
On obtient un nouveau système symétrique :

L2’ [ bb 1 ] Y + [ bc 1 ] Z+ [ bk 1 ]=0

L3’ [ bc 1 ] Y + [ cc 1 ] Z + [ ck 1 ] =0
2ème étape : Nous allons exprimer l’inconnue Y en fonction de l’inconnues Z. Pour cela, nous allons
considérer la ligne L2’.

−[ bc 1 ] [ bk 1 ]
Y= Z−
[ bb 1 ] [ bb 1 ]

Nous reportons Y dans la ligne L3’.

[ bc 1 ] [ [[
− bc 1 ]
bb 1 ]
Z−
[ bk 1 ]
[ bb 1 ] ]
+ [ cc 1 ] Z + [ ck 1 ]=0

[ [ cc 1 ] − [ bc 1 ] ] [
[ bc 1 ]
[ bb 1 ]
Z + [ ck 1 ] −[ bc 1 ]
]
[ bk 1 ]
[ bb 1 ]
=0

[ ck 2 ]
[ cc 2 ]
On obtient enfin :

[ cc 2 ] Z+ [ ck 2 ] =0

D’où :
− [ ck 2 ]
Z=
[ cc 2 ]

Le système proposé devient donc :

[ aa ] X + [ ab ] Y + [ ac ] Z + [ ak ] =0
[ bb 1 ] Y + [ bc 1 ] Z+ [ bk 1 ]=0
[ cc 2 ] Z+ [ ck 2 ] =0

La dernière équation donne la valeur de l’inconnue Z et l’on calcule de proche en proche les
autres inconnues.

50 /77
 Rappel des expressions suivantes des différents coefficients :
[ ab ] [ ak ]
[ bb 1 ] =[ bb ]− [ ab ] [ ck 1 ] =[ ck ] −[ ac ]
[ aa ] [ aa ]

[ ac ] [ bc 1 ]
[ bc 1 ]= [ bc ]− [ ab ] [ cc 2 ] =[ cc 1 ]− [ bc 1 ]
[ aa ] [ bb 1 ]
[ ak ] [ bk 1 ]
[ bk 1 ] =[ bk ] −[ ab ] [ ck 2 ] =[ ck 1 ] − [ bc 1 ]
[ aa ] [ bb 1 ]

[ ac ]
[ cc 1 ] =[ cc ] −[ ac ]
[ aa ]

Le tableau général de substitution appelé tableau de Doolittle est le suivant :

Tableau de Doolittle :
L0
X Y Z k ∑
L1 [ aa ] [ ab ] [ ac ] [ ak ] [ as ]

−[ ab ] −[ ac ] −[ ak ] −[ as ]
L2 −1 [ aa ] [ aa ] [ aa ] [ aa ]

L3 Report [ bb ] [ bc ] [ bk ] [ bs ]

[ ab ] [ ac ] [ ak ] [ as ]
L4 −[ ab ] −[ ab ] −[ ab ] −[ ab ]
S2 (X) [ aa ] [ aa ] [ aa ] [ aa ]

L5 Résolvante [ bb 1 ] [ bc 1 ] [ bk 1 ] [ bs 1 ]

−[ bc 1 ] −[ bk 1 ] −[ bs 1 ]
L6 −1 [ bb 1 ] [ bb1 ] [ bb 1 ]

L7 Report [ cc ] [ ck ] [ cs ]

[ ac ] [ ak ] [ as ]
L8 −[ ac ] −[ ac ] −[ ac ]
S3 (X) [ aa ] [ aa ] [ aa ]

[ bc 1 ] [ bk 1 ] [ bs 1 ]
L9 −[ bc 1 ] −[ bc 1 ] −[ bc 1 ]
S3 (Y) [ bb 1 ] [ bb 1 ] [ bb 1 ]

L10 Résolvante [ cc 2 ] [ ck 2 ] [ cs 2 ]

−[ ck 2 ] −[ cs 2 ]
L11 −1 [ cc 2 ] [ cc 2 ]

51 /77
Calcul des racines : Vérification des résultats :

− [ ck 2 ] −[ cs 2 ]
Z= Z' = =Z−1
[ cc 2 ] [ cc 2 ]

−[ bc 1 ] [ bk 1 ] −[ bc 1 ] [ bs 1 ]
Y= Z− Y '= Z '− =Y −1
[ bb 1 ] [ bb 1 ] [ bb1 ] [ bb 1 ]

− [ ab ] [ ac ] [ ak ] −[ ab ] [ ac ] [ as ]
X= Y− Z− X '= Y '− Z '− =X−1
[ aa ] [ aa ] [ aa ] [ aa ] [ aa ] [ aa ]

Principe de fonctionnement du tableau de DOOLITTLE pour 3 inconnues

Détermination des différents coefficients de chaque ligne du tableau :

L0 = entête du tableau : inconnues à déterminer ; les constantes ki et les différentes sommes

L1 = coefficients de la 1ère équation finale se trouvant dans la table des équations finales

L1
L2 =
−[ aa ]

L3=coefficients de la 2 ème équation finale se trouvant dans la table des équations finales

L4 =[ ab ] × L2

L5=L3 + L4

L5
L6 =
−[ bb 1 ]

L7=coefficients de la 3 ème équation finale se trouvant dans la table des équations finales

L8= [ ac ] × L2

L9= [ bc 1 ] × L6

L10=L7 + L8 + L9

L1
L11=
−[ cc 2 ]

52 /77
 Tableau de DOOLITTLE pour la détermination de 2 inconnues
L0
X Y k ∑
L1 [ aa ] [ ab ] [ ak ] [ as ]

−[ ab ] −[ ak ] −[ as ]
L2
−1 [ aa ] [ aa ] [ aa ]
L3 Report [ bb ] [ bk ] [ bs ]

[ ab ] [ ak ] [ as ]
L4 −[ ab ] −[ ab ] −[ ab ]
S2 (X) [ aa ] [ aa ] [ aa ]

L5 Résolvante [ bb 1 ] [ bk 1 ] [ bs 1 ]

−[ bk 1 ] −[ bs 1 ]
L6 −1 [ bb1 ] [ bb 1 ]

Calcul des racines : Vérification des résultats :

−[ bk 1 ] −[ bs 1 ]
Y= Y '= =Y −1
[ bb 1 ] [ bb 1 ]

− [ ab ] [ ak ] −[ ab ] [ as ]
X= Y− X'= Y '− =X −1
[ aa ] [ aa ] [ aa ] [ aa ]

Principe de fonctionnement du tableau de DOOLITTLE pour 2 inconnues

Détermination des différents coefficients de chaque ligne du tableau :

L0 = entête du tableau : inconnues à déterminer, les constantes ki et les différentes sommes

L1 = coefficients de la 1ère équation finale se trouvant dans la table des équations finales

L1
L2 =
−[ aa ]

L3=coefficients de la 2 ème équation finale se trouvant dans la table des équations finales

L4 =[ ab ] × L2

L5=L3 + L4

L5
L6 =
−[ bb 1 ]

53 /77
 APPLICATIONS

APPLICATION 1

Soit le système des relations d’observations ci-dessous :

1.0596 X + 0.7072Y + Z+ 0.83=0
0.6124 X−2.7157 Y + Z +0.80=0
−2.4554 X +0.7317 Y + Z+ 0.72=0
−0.2054 X−2.0941 Y + Z−1.54=0
0.6736 X + 1.6504 Y +Z — 0.80=0
Déterminer X, Y et Z par la méthode des moindres carrés.

RESOLUTION

 Tableau des relations d’observations

N° X Y Z ki ∑
1 1,0596 0,7072 1 0,83 3,5968
2 0,6124 -2,7157 1 0,8 -0,3033
3 -2,4554 0,7317 1 0,72 -0,0037
4 -0,2054 -2,0941 1 -1,54 -2,8395
5 0,6736 1,6504 1 -0,8 2,524

 Table des équations finales

X Y Z ki ∑
8,0227 -1,1685 -0,3152 -0,6211 5,9179
15,5196 -1,7205 0,8458 13,4764
5,0000 0,0100 2,9743

54 /77
  Tableau de DOOLITTLE

X Y Z ki ∑
8,0227 -1,1685 -0,3152 -0,6211 5,9179
-1,0000 0,1457 0,0393 0,0774 -0,7376
15,5196 -1,7205 0,8458 13,4764
-0,1702 -0,0459 -0,0905 0,8620
15,3494 -1,7664 0,7554 14,3384
-1,0000 0,1151 -0,0492 -0,9341
5,0000 0,0100 2,9743
-0,0124 -0,0244 0,2325
-0,2033 0,0869 1,6501
4,7843 0,0725 4,8569
-1,0000 -0,0152 -1,0152

Z=−0.0152

Y =0.1151∗(−0.0152 )−0.0492=−0.0510

X =0,1457∗(−0.0510 ) +0,0393∗(−0.0152)=0.0694

Vérification

Z '=−1.0152

Z−1=−0.0152−1−¿−1.0152

Y '=0.1151∗(−1.0152 )−0,9341=−1.0510

Y −1=−0.0510−1−¿−1.0510

X ' =0,1457∗(−1.0510 ) +0,0393∗(−1.0152)=0.9306

X −1=0.0694−1−¿ 0.9306

APPLICATION 2

Soit le système des relations d’observations ci-dessous :

−1.9353 X +0.6522 Y =0
−0.8909 X −2.2826 Y =0
3.4980 X −2.1462Y −4.16=0
1.8474 X+ 1.7924 Y +0.42=0

Déterminer X, Y et Z par la méthode des moindres carrés.

55 /77
RESOLUTION

 Tableau des relations d’observations

ki
Point X Y ∑
A -1,9353 0,6522 0,00 -1,2831
B -0,8909 -2,2826 0,00 -3,1736
C 3,4980 -2,1462 -4,16 -2,8116
D 1,8474 1,7924 0,42 4,0611

 Table des équations finales

X Y ki ∑
20,1880 -3,4246 -13,7854 2,9780
13,4545 9,6906 19,7205

 Tableau de DOOLITTLE

X Y ki ∑
20,1880 -3,4246 -13,7854 2,9780
-1,0000 0,1696 0,6829 -0,1475
13,4545 9,6906 19,7205
-0,5809 -2,3385 0,5052
12,8736 7,3521 20,2257
-1,0000 -0,5711 -1,5711

Y =−0,5711

X =0,1696∗(−0,5711 )−0,6829=0,5860

Vérification

Y '=−1,5711

Y −1=−0,5711−1−¿−1.5711

X ' =0,1696∗(−1.5711) −0,1475=−0,4140

X −1=0,5860−1−¿−0,4140

56 /77
 II. METHODE DE VARIATION DES COORDONNEES

Compenser une figure géodésique, c’est substituer au réseau approché obtenu un réseau
géométrique bien défini s’appuyant sur l’ensemble des points.

Notation

 Points anciens : ce sont des points dont les coordonnées sont connues ; elles ne peuvent
être modifiées. On les note en général A, B, C ....
 Points nouveaux : ce sont les points à déterminer ; on les note I ; II ; III ….
 Coordonnées des points approchés : on les note X0 ; Y0.
 Corrections à apporter aux coordonnées approchées : on les note dX ; dY.
 δX et δY : inconnues, auxiliaires se rapportant à dX et dY.
 Gobs : gisement observé : Gobs =G 0+lecture L
∆X
 Gcal : gisement calculé : Gcal=arctan
∆Y
 G 0: gisement du zéro du limbe
 G 0 prov : G0 provisoire ou G 0moyen
 dG 0: correction à apporter au G 0 prov pour avoir le G 0déf (G0 définitif )
 δG0 : inconnue, auxiliaire se rapportant à dG 0.

1. INTERSECTION

Y
Visée sur M
Gobs

G
M 0 point approché
Gcal dX dY
v
G0 Dcal M point définitif
L
0

La visée de A vers M donne :

X M −X A X M −X A + dX ∆ X 0+ dX
G=arctan =arctan 0
=arctan
Y M −Y A Y M −Y A +dY
0
∆ Y 0 +dY

Le développement de Taylor limité au 1er degré nous permet de linéariser la fonction arctan qui
nous donne :

57 /77
 ( ) ( )
' '
∆ X0 ∆ X0 ∆ X0
G=arctan + arctan . dX+ ¿ arctan . dY ¿
∆Y 0 ∆Y 0 X ∆Y0 Y

∆ X0
Or G=arctan =Gcal
∆Y 0

( )
'
∆ X0 1 1 ∆Y 0 ∆Y0
arctan = × = = 2

( )
2 2 2
∆Y0 X ∆ X0 ∆ Y0 ∆ X 0 +∆ Y0 Dcal
1+
∆Y0

( ) ( )
'
∆ X0 1 −∆ X 0 ∆ X0 −∆ X 0
arctan = × = =

( )
2 2 2 2 2
∆Y0 Y ∆ X0 ∆Y0 ∆ X 0 +∆ Y0 D cal
1+
∆Y0

D’où

∆Y 0 ∆ X0
G=G cal+ 2
dX − 2
dY =Gobs + v
Dcal D cal

∆Y0 ∆ X0
2
dX− 2
dY −(G obs−G ¿ ¿ cal)=v ¿
D cal Dcal

Harmonisation des unités.

Le résidu v étant la quantité à ajouter au gisement observé Gobs pour avoir le gisement définitif G,
nous allons l’exprimer en seconde centésimale. Pour cela, Gobs −Gcal doit être en seconde
1
centésimale. Cela entraîne que les unités linéaires doivent être multipliées par ̏
sin 1

D’où :
1
̏
∆Y0
2
sin 1 D cal( ∆ X0
dX− 2 dY −(Gobs −G ¿ ¿ cal)̏ =v ̏ ¿
Dcal )
( )
3 3
∆Y 0 ( km ) . 10 ∆ X 0 ( km ) .10
636620 dX − dY −(Gobs −G ¿ ¿ cal)̏ =v ̏ ¿
D2cal ( km ) . 106 D2cal ( km ) . 106

6,36620. 10
2
( ∆ Y 0 ( km )
2
D ( km )
cal
dX −
∆ X 0 ( km )
D 2cal ( km ) )
dY −(G obs−G ¿ ¿ cal)̏ =v ̏ ¿

6,36620
( 10 ∆ Y 0 ( km )
2
D ( km )
cal
dX−
10 ∆ X 0 ( km )
2
D ( km )
cal
dY −
) (Gobs −Gcal )̏ v ̏
10
=
10

6,36620
( 10 cos G cal
Dcal ( km )
dX−
10 sin Gcal
D cal ( km )
dY −
(Gobs −Gcal )̏ v ̏
10
=
10 )
Posons:
(Gobs −Gcal )̏ 3
δY =6,36620 dY ḱ = =10 ∗(Gobs −Gcal )
10
58 /77
δX =6,36620 dX

10cos Gcal −10 sinG cal

a= b=
Dcal ( km ) Dcal ( km )
La relation d’observation d’une vise d’intersection s’écrit donc :
vi ̏
aδX +bδY −k i ̏ =
10

 Calcul des coefficients a, b et –ki

Points Gobs ( gr) Gcal ( gr ) Dcal (km) a b −k i ̏

approch anciens
é
A
B
M
C
D
NB:

 Les distances Dcal (km) seront données à 6 chiffres après la virgule ;

 les gisements à 5 chiffres,
 les coefficients a et b à 4 chiffres et
 –ki à 2 chiffres après la virgule.
 Calcul des vi

Avant compensation : v i ̏ =10 k i ̏ après compensation : v i ̏ =10(aδX +bδY −k i ̏ )

vi ̏ vi ² vi ̏ vi ²
v1 v1
v2 v2
v3 v3
v4 v4
∑vi ² ∑vi ²
Précision angulaire :

√ √
n
n

∑ vi ² ∑ vi ²
i
σα̏ = i σα̏ =
n n
Calcul des rayons avant et après compensation : v i ̏ × Di (km)
r i ( cm )= cal

6,36620

Avant compensation : après compensation

r i (cm) ri ² r i (cm) ri ²
r1 r1
59 /77
r2 r2
r3 r3
 A B

Précision sur le rayon moyen quadratique :

D
C

√ √
n
n

∑ ri ² ∑ ri ²
i
σ r (cm)= i σ r (cm)=
n n

METHODOLOGIE :

- Calculer les coordonnées approchées du point M à partir de deux points anciens ;

- Déterminer les gisements observés à partir du calcul de gisements de points anciens et
des lectures angulaires puis déterminer les gisements calculés et les distances calculées et
enfin calculer les coefficients a, b et –k ;
- Etablir le tableau des relations d’observations ;
- Etablir la table des équations finales ;
- Etablir le tableau de Doolittle et calculer les dX et dY
- Calculer les coordonnées définitives
- Calculer les vi avant et après compensation et en déduire la précision angulaire
- Calculer les différents rayons et en déduire la précision
sur le rayon moyen quadratique

APPLICATION

Stations Points Lectures X (m) Y (m)

visés
A B 0.0000 785 489.74 315 558.44
M 45.4548
B M 0.0000
A 57.5437 782 333.32 310 192.99
a) Calculer les coordonnées C D 0.0000 789 761.91 309 051.37
approchées du point M à partir M 76.2305
des points A et B. D M 0.0000
C 39.7009 789 050.55 313 030.98
b) Calculer les coordonnées
définitives du point M
c) Calculer les précisions angulaires et les précisions sur le rayon moyen quadratique avant et
après compensation

RESOLUTION

1. Calcul des coordonnées approchées du point M à partir des points A et B

60 /77
  Calcul des angles A et B
^
A=l M −l B=¿ ……………….. ^B=l A −l M =¿……………………..

 Calcul des gisements G AB , G AM et GBM

X B− X A
G AB=arctan ⁡( )=… … … … … … … … …. .
Y B−Y A
G AM =G AB + ^
A=… … … … … … … … …. ^
G BM =GBA− B=… …………………..

 Calcul des coordonnées approchées du point M

X B− X A−( Y B −Y A )∗tan G BM
Y M =Y A + =… … … … … … … … .
0
tan G AM −tan GBM

X M = X A + ( Y M −Y A )∗tan G AM =… … … … … … … …. .
0

2. Calcul des coordonnées définitives du point M

 Calcul des gisements observés

^ M −l D =¿………………..
C=l ^
D=l C −l M =¿……………………..

^
GCM =GCD + C=… ……………………. G DM =G DC − ^
D=… … … … … … … … ..

 Calcul des coefficients a, b et -k.

vi ̏
La relation d’observation d’une visée d’intersection est : aδX +bδY −k i ̏ =
10

Points Gobs (gr) Gcal (gr ) Dcal (km) a b −k i ̏

appro. anc.
M A
B
C
D

−10 sinG cal 10cos Gcal 3

−ḱ i=10 ∗(Gcal−G obs )
b= a=
Dcal ( km ) Dcal ( km )

Gcal=arctan ⁡ ( )∆X
∆Y
Dcal =√ ∆ X 2 +∆ Y 2 ⁡
61 /77
  Tableau des relations d’observation

N° δX δY -ki ∑
1 ….…………. ….…………. ….…………. ….………….

2 ….…………. ….…………. ….…………. ….………….

… … … … …. . … … … … …. . … … … … …. . … … … … …. .
3
… … … … …. . … … … … …. . … … … … …. . … … … … …. .
4

 Table des équations finales

δX δY -ki ∑

 Tableau du tableau de DOOLITTLE et résolution

δX δY -ki ∑
…………. ………… ……. ………… …… ………… …….

… … … … …. ………… ……… … … … … …. ………… …….

Report ………… ……. ………… ……. ………… …….

S2 (X) ………… ……… ………… ……… ………… ………

Résolvante ………… ……. ………… ……. ………… …….

………… ……… ………… ……… ………… ………

Racines

δ Y = ………….

δ X=… … … … … … … … … … .=… … …. .

Vérification
'
δ Y =… … … … … δ Y −1=… … … … …

62 /77
 '
δ X =… … … … … … … … … … … ..=… … … … … .
δ X−1=… … … … … … … … … …=… … … … …

 Calcul des variations dX et dY

δX
dX (m)= =… … … … .
6,36620

δY
dY (m)= =… … … … .
6,36620

 Coordonnées définitives du point M

X M = X M +dX =… … … … … … … … …
0

Y M =Y M + dY =… … … … … … … … …
0

3. Calcul des précisions angulaires et les précisions sur le rayon moyen quadratique avant et
après compensation
 Précisions angulaires avant et après compensation

Avant compensation : v i ̏ =10 k i ̏ après compensation : v i ̏ =10(aδX +bδY −k i ̏ )

vi ̏ vi ² vi ̏ vi ²

∑vi ² ∑vi ²

Précision angulaire avant compensation : Précision angulaire après compensation :

√ √
n n

∑ vi ² ∑ vi ²
i i
σα̏ = =… … … … … … σα̏ = =… … … … … … .
n n

 Calcul des rayons avant et après compensation :

v i ̏ × Di (km)
r i ( cm )= cal

6,36620

Avant compensation : après compensation

r i (cm) ri ² r i (cm) ri ²

∑r i ² ∑r i ²
63 /77
 Précision sur le rayon moyen quadratique :

Avant compensation Après compensation

√ √
n
n
∑ ri ²
∑ ri ² σ r (cm)= i
=… … … … …
σ r (cm)= i
=… … … … … n
n

2. RELEVEMENT

Les gisements G AM de l’intersection et de G MA de relèvement, visée de M vers A diffèrent de 200

gr et par conséquent :

( )
'
∆ X0 −∆ Y 0
( )
'
arctan = 2 ∆ X0 ∆ X0
∆Y0 X et arctan = 2
D cal ∆ Y 0 Y Dcal

En outre, le théodolite en station au point M génère une variable supplémentaire, le G0 tel que :

G 0=G 0 prov +dG 0=G 0moy + dG0 .

On aura donc Gobs −Gcal =(G 0+ L)−Gcal=G0 moy +dG 0+ L−Gcal

Gobs −Gcal =dG 0+G 0moy −(G¿¿ cal−L)=dG 0+G 0moy −G 0i ¿

( Gobs −Gcal ) ̏ dG 0 ̏ ( G 0moy −G 0i ) ̏

= +
10 10 10

En posant que :
( G 0 moy−G 0i ) ̏ dG 0 ̏
3 δG 0̏ =
ḱ i= =10 ∗(G moy−G 0i ) 10
10
D’où avec les notations précédentes, la relation d’observation d’une visée de relèvement
devient :
vi ̏
−aδX −bδY +δG 0̏ −k i ̏ =
10

 Calcul des coefficients - a, - b et - ki

Points Gcal ( gr ) Dcal (km) G 0i (gr ) G 0moy ( gr) −a −b −k i ̏

approché anciens
A
B
M ...........
C
D

64 /77
Remarque : Le coefficient de δG 0̏ est 1

 Calcul des vi

Avant compensation : v i ̏ =10.k i ̏ après compensation : v́ i=10(−aδX−bδY +δG 0 ̏ −k i ̏ )

vi ̏ vi ² vi ̏ vi ²
v1 v1
v2 v2
v3 v3
v4 v4
∑vi ² ∑
vi ² ̏
Précision angulaire :

√ √
n
n

∑ vi ² ∑ vi ²
i
σα̏ = i σα̏ =
n n
v i ̏ × Di (km)
Calcul des rayons avant et après compensation : r i ( cm )= cal

6,36620

Avant compensation : après compensation

r i (cm) ri ² r i (cm) ri ²
r1 r1
r2 r2
r3 r3
r4 r4
∑r i ² ∑r i ²

Précision sur le rayon Précision sur le rayon

moyen quadratique : moyen quadratique :

√ √
n n
∑ ri ² ∑ ri ²
i
σ r (cm)= σ r (cm)= i
n n

METHODOLOGIE :

- Calculer les coordonnées approchées du point M à partir de deux points anciens en

utilisant l’une des trois méthodes de calcul de relèvement
- Calculer les gisements et les distances calculées ; les G0i et le G0moyen et les coefficients –a,
–b et k ;
- Etablir le tableau des relations d’observations ;
- Etablir la table des équations finales ;
- Etablir le tableau de Doolittle et calculer les dX, dY et dG0
- Calculer les coordonnées définitives et le G0 définitif ;

65 /77
 - Calculer les vi avant et après compensation et en déduire la précision angulaire
- Calculer les différents rayons et en déduire la précision sur le rayon moyen quadratique

APPLICATION

Sachant que l’opérateur a stationné en un point M et a visé successivement les points connus A,
B et C, on se propose de calculer les coordonnées définitives de ce point M.

A B Station Points Lectures X (m) Y (m)

visés azimutales (gr)
M
A 150,6978 12271,85 4683,26

B 237,5126 12819,31 4737,39

D
M C 331,2076 12994,39 4138,31
C D 64,7218 12103,72 4092,39

A 150,6997
1. Calculer les coordonnées
approchées du point M à partir des points A, B et C

2. Calculer les coordonnées définitives du point M et le G0 définitif

3. Calculer les précisions angulaires et le rayon moyen quadratique avant et après

compensation

RESOLUTION

1. Calcul des coordonnées approchées du point M à partir des points A et B.

Angles au Angles au Coefficients

Pts Gi (gr) Li (m)
sommet centre mi
A
B

A
B
Σ

1
m i=
cotan Angle au sommet I −cotan angle au centre i

66 /77
 n

∑ mi . X i
X M 0= i=1 n =… … … … … … … … … …
∑ mi
n

∑ mi . Y i
i=1

Y M 0= i=1n =… … … … … … … … … …. .
∑ mi
i=1
2. Calcul des coordonnées définitives du point M
 Calcul des coefficients – a, – b et – ki
vi ̏
La relation d’observation d’une visée de relèvement s’écrit : −aδX −bδY +δG 0̏ −k i ̏ =
10
Points Gcal ( gr ) Dcal (km) G 0i (gr ) G 0moy ( gr) −a −b −k i ̏
appro Anc.
A
B
M ..................
C
D

10sin Gcal −10 cos Gcal 3

−ḱ =10 ∗(G 0i−Gmoy )
−b= −a=
Dcal ( km ) D cal ( km )

 Tableau des relations d’observation

N° δX δY δG 0 -ki ∑
1 ….…………. ….…………. … … … … …. . ….…………. ….………….
G cal=arctan ⁡
2
∆X
… . … …∆Y…….
( ) Dcal =√ ∆ X 2 +∆ Y 2 ⁡
….…………. ….………….
G 0i=Gcal −Li
….…………. ….………….

… … … … …. . … … … … …. . ….…………. … … … … …. . … … … … …. .
3
… … … … …. . … … … … …. . … … … … …. . … … … … …. . … … … … …. .
4

 Table des équations finales

67 /77
 δX δY δG 0 -ki ∑

 Tableau du tableau de DOOLITTLE et résolution

δX δY δG 0 -ki ∑
…………. ………… ……. ………… ……. ………… …… ………… …….

… … … … …. ………… ……… ………… ……. … … … … …. ………… …….

Report ………… ……. ………… ……. ………… ……. ………… …….

S2 (X) ………… ……… ………… ……. ………… ……… ………… ………

Résolvante ………… ……. ………… ……. ………… ……. ………… …….

………… ……… ………… ……. ………… ……… ………… ………

Report ………… ……. ………… ……. ………… …….

S3 (X) ………… ……. ………… ……. ………… …….

S3 (Y) ………… ……. ………… ……. ………… …….

Résolvante ………… ……. ………… ……. ………… …….

………… ……. ………… ……. ………… …….

Racines

δG 0 = ………….

δ Y =… … … … … … … … … … … … … … … … … ..=… … …. .

δ X=… … … … … … … … … … … … … … … … … … … … … … … … … … …=… … … . .

Vérification

68 /77
δG 0' =… … … … … δG 0−1=… … … … … … … … .=… … … … … . .
'
δ Y =… … … … … … … … … … … … … … … …=… … … … … .
δ Y −1=… … … … … … … … …=… … … … …
'
δ X =… … … … … … … … … … … … … … … … … … … … … … ..=… … … … … .
δ X−1=… … … … … … … … …=… … … … …

 Calcul des variations dX, dY et dG0

δX δY
dX (m)= =… … … … . dY (m)= =… … … … .
6,36620 6,36620

dG 0(dmgr)=δ G 0× 10=… … … … .

 Coordonnées définitives du point M et le G0 définitif

X M = X M +dX =… … … … … … … … …
0

Y M =Y M + dY =… … … … … … … … …
0

G 0definitif =G 0moy + dG 0=… … … … … … … … …

3. Calcul des précisions angulaires et le rayon moyen quadratique avant et après compensation
 Calcul des précisions angulaires avant et après compensation

Avant compensation : v i ̏ =10.k i ̏ après compensation : v́ i=10(−aδX−bδY +δG 0 ̏ −k i ̏ )

vi ̏ vi ² vi ̏ vi ²

∑vi ² ∑vi ² ̏
Précision angulaire avant compensation : Précision angulaire après compensation

√ √
n
n

∑ vi ² ∑ vi ²
i
σα̏ = i
=… … … … . . σα̏ = =… … … … … .
n n

v i ̏ × Di (km)
r i ( cm )= cal

6,36620 69 /77
Calcul des rayons avant et après compensation :

Avant compensation : après compensation

r i (cm) ri ² r i (cm) ri ²

∑r i ² ∑r i ²

Précision sur le rayon moyen Précision sur le rayon moyen

quadratique : quadratique :

√ √
n
n

∑ ri ² ∑ ri ²
i
σ r (cm)= i
=… … … … σ r (cm)= =… … … … .
n n

70 /77
 3. RECOUPEMENT

Il faut au total au moins 6 visées dont 3 visées d’intersection et 3 visées de relèvement

Relations d’observations :
vi ̏
 Visées d’intersection : aδX +bδY −k i ̏ =
10
vi ̏
 Visées de relèvement : −aδX −bδY +δG 0̏ −k i ̏ =
10

 Calcul des coefficients a, b, -a, -b et –ki

Points Gobs (gr) Gcal (gr ) Dcal (km) G 0i (gr ) G 0moy (gr) a ou b ou δG 0̏ −ki ̏
appro anciens -a -b
Relèvement
........... 1
1
1
M
Intersection
0
0
0

3
−ki ̏ =10 ∗(G0 i−Gmoy ) pour ≤relèvement
NB:
3
−ki ̏ =10 ∗(Gcal−Gobs ) pour l '∩¿

 Calcul des vi

NB : le calcul des v i tient compte de l’opération effectuée soit le relèvement, soit l’intersection.

Avant compensation : v i ̏ =10 k i ̏ après compensation : v i ̏ =10(aδX +bδY −k i ̏ ) ou

v́ i=−10(aδX +bδY −δG 0 ̏ +k i ̏ )

vi ̏ vi ² vi ̏ vi ²
v1 v1
v2 v2
v3 v3
v4 v4
∑vi ² ∑vi ²
Précision angulaire :

√ √
n
n

∑ vi ² ∑ vi ²
i
σα̏ = i σα̏ =
n n
71 /77
 v i ̏ × Di (km)
Calcul des rayons avant et après compensation : r i ( cm )= cal

6,36620

Avant compensation : après compensation

r i (cm) ri ² r i (cm) ri ²
r1 r1
r2 r2
r3 r3
r4 r4
∑r i ² ∑r i ²
Précision sur le rayon
moyen quadratique :

√
n

∑ ri ²

√
n

∑ ri ² σ r (cm)= i
n
i
σ r (cm)=
n

METHODOLOGIE :

- Calculer les gisements observés à partir du calcul de gisements de points anciens et des
lectures angulaires ;
- Calculer les coordonnées approchées du point M soit par intersection ou soit par
relèvement
- Déterminer les gisements et distances calculés ;
- Calculer les G0i et le G0moyen
- Calculer les coefficients a, b et –k ;
- Etablir le tableau des relations d’observations ;
- Etablir la table des équations finales ;
- Etablir le tableau de Doolittle et calculer les dX, dY et dG0
- Calculer les coordonnées définitives et G0 définitif
- Calculer les vi avant et après compensation et en déduire la précision angulaire
- Calculer les différents rayons et en déduire la précision sur le rayon moyen quadratique

72 /77
 4. MULTILATERATION
La Multilatération est la détermination de points par la mesure des distances. Les notions
d’intersection et de relèvement n’existent pas ici. Tout point nouveau doit recevoir au moins 3
visées (distances).
Y
M 0 point approché

Gcal dX dY
Dcal
M point définitif
D

D= AM étant la distance au point définitif, il vient :

D= √ ( X M −X A ) ² + ( Y M −Y A ) ²= √ [( X
2 2
M0 − X A ) + dX ] + [ ( Y M −Y A ) + dY ]
0

Soit en développant :

√( X √(X √( X
' '
D= M0 −X A ) ²+ ( Y M −Y A ) ²+ ( M0 −X A ) ²+ ( Y M −Y A ) ² ) .dX + ( M0 − X A ) ² + ( Y M −Y A ) ² ) . dY
0 0 dX 0 dY

Or √ ( X M − X A ) ²+ ( Y M −Y A ) ² ¿ D cal
0 0

2 ( XM − X A) Dcal .sin Gcal

( √ ( X M −X A ) ²+( Y M −Y A ) ² )dX =
'
= =sin Gcal
0

0 0
2 √( X M0 −X A ) ²+ ( Y M −Y A ) ²
0
Dcal

2 ( Y M −Y A ) Dcal .cos Gcal

(√ ( X
'
−X A ) ²+ ( Y M −Y A ) ² )dY = = =cos Gcal
0

√( X
M0 0
2 M0 −X A ) ²+ ( Y M −Y A ) ² Dcal
0

D’où D=D obs + v D =Dcal +sin Gcal . dX +cos G cal . dY

Soit : sin Gcal . dX +cos Gcal . dY −( Dobs −Dcal )=v D

Posons :

a=sin Gcal b=cos Gcal k i=Dobs −Dcal

Alors la relation d’observation d’une mesure de distance s’écrit : adX + bdY −k i=v i

 Calcul des coefficients a, b et –ki

73 /77
 Points Gcal ( gr ) Dcal (m) Dobs (m) a b −k i (m)
approché anciens

NB :

 Les distances Dcal (m) seront données à 2 chiffres après la virgule ;

 Les gisements à 4 chiffres,
 Les coefficients a et b à 4 chiffres et
 –ki à 2 chiffres après la virgule.

 Calcul des vi

Avant compensation : v i=k i après compensation : v i=adX +bdY −k i

vi vi ² vi vi ²
v1 v1
v2 v2
v3 v3
v4 v4
∑vi ² ∑vi ²

Précision sur la distance :

√ √
n
n

∑ vi ² ∑ vi ²
i
σ D (m)= i σ D (m)=
n n

METHODOLOGIE :

- Calculer un gisement observé à partir de deux points anciens et la distance entre ces
points ;
- Résoudre le triangle formé par les deux points anciens utilisés ci-dessus et le point à
déterminer ; il s’agira de calculer les 3 angles par la formule des cosinus
- Déterminer les deux gisements des 2 directions formées par le point à déterminer et les
2 points anciens ;
- Calculer les coordonnées approchées du point M à partir des deux points anciens ;
- Déterminer tous les gisements et toutes les distances calculées
- Calculer les coefficients a, b et –k ;
- Etablir le tableau des relations d’observations ;
- Etablir la table des équations finales ;
- Etablir le tableau de Doolittle et calculer les dX et dY
74 /77
 - Calculer les coordonnées définitives
- Calculer les vi avant et après compensation et en déduire la précision sur la distance ou
précision linéaire.

APPLICATION

En vue de déterminer les coordonnées du point M, un technicien supérieur géomètre, muni

d’une chaîne a mesuré les distances suivantes :

AM =232,989 m; BM = 224,851 m; CM = 127,442 m; DM = 218,701 m

Coordonnées des points d’appui A B

Points X (m) Y (m)

A 8 911,95 5 122,79 M
B 9 212,36 5 103,07
C 9 086,44 4 816,49
D 8 845,62 4 884,46 D

1. Calculer les coordonnées approchées du point M à partir des points A et D C

2. Calculer les coordonnées définitives du point M
3. Calculer les précisions linéaires avant et après compensation.
RESOLUTION

1. Calcul des coordonnées approchées du point M à partir des points A et D

 Calcul de G AD et AD.

X D −X A D AD =√ ∆ X 2 +∆ Y 2 ⁡=… … … …
G AD=arctan ⁡( )=… … … … … … … … … . .
Y D −Y A

 Calcul des angles du triangle ADM

( )
2 2 2
^ A M + A D −D M
A=arccos =… … … … … … … …
2 AM . AD

( )
2 2 2
^ DA + DM − AM
D=arccos =… … … … … … … . .
2 DA . DM

( )
2 2 2
^ MA + MD −AD
M =arccos =… … … … … … … …
2 MA . MD

Contrôle : ^ ^+M
A+ D ^ =… … … … … … … … … .

 Calcul de G AM et G DM

G AM =G AD − ^
A=… … … … … … … … … . G DM =G DA + ^
D=… … … … … … … … ..

75 /77
  Calcul des coordonnées approchées du point M

X D −X A −( Y D−Y A )∗tan G DM
Y M =Y A + =… … … … … … … … .
0
tan G AM −tanG DM

X M = X A + ( Y M −Y A )∗tan G AM =… … … … … … … …. .
0

2. Calcul des coordonnées définitives du point M

 Calcul des coefficients a, b et -ki.

La relation d’observation d’une visée d’intersection est : adX + bdY −k i=v i

Points Gcal ( gr ) Dcal (m) Dobs (m) a b −k i( m)

appro anc
A
B
M
C
D

Gcal=arctan ⁡ ( )
∆X
∆Y
Dcal =√ ∆ X 2 +∆ Y 2 ⁡

a=sin Gcal b=cos Gcal −k =Dcal−D obs

 Tableau des relations d’observation

N° dX dY -ki ∑
1 ….…………. ….…………. ….…………. ….………….

2 ….…………. ….…………. ….…………. ….………….

… … … … …. . … … … … …. . … … … … …. . … … … … …. .
3
… … … … …. . … … … … …. . … … … … …. . … … … … …. .
4

 Table des équations finales

dX dY -ki ∑

76 /77
  Tableau du tableau de DOOLITTLE et résolution

dX dY -ki ∑
…………. ………… ……. ………… …… ………… …….

… … … … …. ………… ……… … … … … …. ………… …….

Report ………… ……. ………… ……. ………… …….

S2 (X) ………… ……… ………… ……… ………… ………

Résolvante ………… ……. ………… ……. ………… …….

………… ……… ………… ……… ………… ………

Racines

d Y = ………….

d X =… … … … … … … … … … .=… … … . .

Vérification
'
d Y =… … … … … d Y −1=… … … … …
'
d X =… … … … … … … … … … … ..=… … … … … .
d X −1=… … … … … … … … … …=… … … … …

 Coordonnées définitives du point M

X M = X M +dX =… … … … … … … … …
0

Y M =Y M + dY =… … … … … … … … …
0

3. Calcul des précisions linéaires avant et après compensation

 Précisions linéaires avant et après compensation

Avant compensation : v i=k i après compensation : v i=adX +bdY −k i

77 /77
 vi vi ²
vi vi ²

∑vi ² ∑vi ²

Précisions linéaires :

Avant compensation : Après compensation :

√ √
n n

∑ vi ² ∑ vi ²
i i
σ D= =… … … … … … σ D= =… … … … … … .
n n

78 /77