Tema 3
Tema 3
Tema 3
1
etc.), geodezie matematică (Jordan 1958 etc.), geodezie geometrică (Heiskanen 1969), obiectul
de studiu rămânând, în principiu, acelaşi.
Deoarece au fost utilizaţi de-a lungul timpului mai mulţi elipsoizi de rotaţie ca referinţă,
o problemă importantă pentru geodezie o reprezintă transcalculul de coordonate de pe un
elipsoid pe altul.
Pentru prezentarea formulelor şi dezvoltărilor ulterioare este necesară prezentarea
parametrilor geometrici prin care se poate defini, un elipsoid de rotaţie (figura 3.1). Parametrii
geometrici ai elipsoidului de rotaţie sunt:
– semiaxa mare (raza ecuatorială); (3.1)
– semiaxa mică; (3.2)
– turtirea geometrică; (3.3)
Z
St
P
He
( 3) normala la
S S
elipsoid
G
E΄
B S ( 2) E
E
(1)
S quL
S0ati
on
.3
X x
P΄
2
Elipsoidul de rotaţie poate fi definit prin doi parametri geometrici, dintre care unul
trebuie să fie liniar (lungime). Primii trei parametri se mai numesc şi parametri principali. În
lucrări mai dezvoltate se utilizează şi alţi parametri geometrici.
Între parametrii geometrici ai elipsoidului de rotaţie se pot stabili cu uşurinţă următoarele
relaţii principale de legătură:
; (3.8)
; (3.9)
; (3.10)
; (3.11)
; (3.12)
; (3.13)
(3.14)
; (3.15)
. (3.16)
Elipsoidul de referinţă, este elipsoidul folosit la un moment dat, într-o ţară sau în mai
multe ţări, pentru rezolvarea problemelor geodezice. Acesta este un elipsoid de rotaţie cu turtire
mică la poli. Reamintind că parametrii geometrici ai elipsoidului de rotaţie au fost definiţi
anterior, în Tabelul 3.1. se prezintă valorile numerice ale parametrilor a şi f pentru elipsoizii de
referinţă care au fost utilizaţi în decursul anilor în ţara noastră, precum şi pentru elipsoidul
caracteristic Sistemului Mondial de Referinţă WGS-84 folosit actualmente în geodezia cu
sateliţi.
Tabelul 3.1. Elipsoizi de referinţă
Denumirea elipsoidului de Anul Semiaxa mare, a
Turtirea numerică f
referinţă determinarii [m]
Bessel 1841 6337397.115 1:299.1528
Clarke 1880 6378243.000 1:293.5
Hayford 1909 6378388.000 1:297.0
Krasovski 1942 6378245.000 1:298.3
WGS-84 1984 6378137.000 1:298.257223563
3
3.2 Sisteme de coordonate utilizate frecvent în geodezia elipsoidală
Zz N
cercul paralelului
P X t , Yt , Z t
S
t
e
punctului S H
S3 B , L , S t
X ,Y , Z
S e
B , L, H 0
G
) 2
O B S ( 2S t Y
E΄
S
St1
1
L SS E
t
S0
x
X
meridianul origine P
meridianul punctelor St şi S
(al punctului Greenwich)
4
În geodezia cu sateliţi se determină coordonatele carteziene Xt Yt Zt ale punctului St de
pe suprafaţa fizică:
; ; . (3.17)
În geodezia elipsoidală se operează cu puncte raportate pe elipsoidul de referinţă, de
genul punctului S ale cărui coordonate carteziene sunt:
; ; .
(3.18)
Sistemul de coordonate geodezice B. L.
Sistemul de coordonate utilizat este acelaşi (X Y Z). Punctele St şi S au însă coordonate
diferite:
a) Latitudinea geodezică B – este unghiul format de normala la elipsoid cu planul
ecuatorului (sensul pozitiv de măsurare a latitudinii B este de la ecuator spre normală). Rezultă
că există latitudini nordice şi respectiv latitudini sudice.
b) Longitudinea geodezică L – este unghiul diedru format de meridianul punctului S t
considerat şi meridianul origine.
Meridianul unui punct reprezintă secţiunea determinată prin elipsoid de către un plan care
trece prin punctul considerat şi prin axa polilor .
Ca meridian origine este acceptat unanim meridianul geodezic al observatorului
astronomic Greenwich. Acest observator este utilizat şi în astronomie (îndeplineşte rolul de
origine pentru longitudinile astronomice şi respectiv pentru timp – aşa numitul timp universal
notat TU).
Atenţie: Coordonatele geodezice B, L nu definesc poziţia în spaţiu a punctului St ci doar
a normalei la elipsoid. Pentru definirea în spaţiu a punctului St mai este necesară o mărime:
c) Altitudinea geodezică :
. (3.19)
Se pot menţiona următoarele localităţi situate la limite extreme pentru teritoriul ţării
noastre:
limita nordică: Otaci ;
limita sudică: Giurgiulești ;
limita de est: Palanca ;
limita de vest: Criva .
Forma și întinderea teritoriului Republicii Moldova
Teritoriul Republicii Moldova se înscrie într-o bandă, care merge aproximativ pe direcția
sud-nord:
5
Pe latitudine:
Pe direcția est-vest:
Punctul S, ca proiecţie pe elipsoid a punctului St în lungul normalei are aceleaşi
coordonate geodezice B, L cu cele ale punctului St. însă He = 0.
Ecuaţiile parametrice ale elipsoidului de rotaţie:
(3.20)
Prin normala într-un punct oarecare de pe suprafaţa elipsoidului pot trece o infinitate
de plane. denumite plane normale (deoarece conţin ). Aceste plane intersectează elipsoidul
după o infinitate de secţiuni. denumite secţiuni normale.
Secţiunile care nu trec prin normală (care nu conţin normala) se numesc secţiuni înclinate
(ex.: secţiunea perpendiculară pe axa polilor, care intersectează elipsoidul după cercul paralel al
punctului respectiv).
Dintre secţiunile normale posibile una are raza de curbură maximă şi una are raza de
curbură minimă; acestea se numesc secţiuni normale principale şi au proprietatea de a fi
perpendiculare între ele.
Secţiunea meridiană. Secţiunea meridiană (care, aşa cum se va arăta în continuare are
raza de curbură minimă) este reprezentată de meridianul punctului S.
Secţiunea meridiană are forma unei elipse. fiind obţinută din intersecţia planului meridian
cu elipsoidul de rotaţie. Raza de curbură a acestei elipse se notează cu M:
. (3.21)
. (3.21)
Expresia razei de curbură a elipsei meridiane se poate determina prin următoarea relație:
6
. (3.22)
. (3.24)
(3.25)
(3.26)
. (3.28)
Din cele de mai sus rezultă semnificaţia geometrică a razei de curbură a primului vertical
N = SS0 în Fig. 3.2.
Se observă că şi raza de curbură a primului vertical are o variaţie de la ecuator la pol.
crescând în mărime odată cu creşterea argumentului B:
La ecuator. unde B = 0. rezultă:
(3.29)
(3.30)
Numai la pol cele două raze de curbură principale M şi N sunt egale. În celelalte situaţii
raportul dintre razele de curbură ale secţiunilor normale principale va fi:
7
. (3.30)
Rezultă N M. motiv pentru care raza de curbură N a primului vertical se mai numeşte şi
marea normală.
Raza de curbură a secţiunii meridiane (M) precum şi raza de curbură a secţiunii primului
vertical (N) se pot calcula în funcţie de parametrii elipsoidului a, e2 şi de latitudinea B a
punctului S (sunt funcţii de poziţia punctului pe un anumit elipsoid dat).
Raza de curbură a unei secţiuni normale oarecare în funcţie de razele de curbură
principale şi azimutul geodezic A
P
paralelul punctului S N
primul
A (c)
vertical
r
S
O
E
B E
S0
meridianul punctului S
P
Figura 3.3. Razele de curbură a secţiunilor normale
Se consideră un punct S pe elipsoid (situat pe un anumit meridian şi pe un anumit
paralel). Prin acest punct se poate duce normala la suprafaţă (o unică posibilitate). Prin
normală trec o infinitate de plane normale. care intersectează elipsoidul după secţiuni normale.
fiecare având o altă rază de curbură.
Azimutul A este unghiul format de curba c (directia ) cu direcţia nord a
meridianului. Se poate demonstra formula Euler de calcul a razei de curbură a unei secţiuni
normale. de azimut A:
. (3.31)
Din relaţia de mai sus se observă că raza de curbură a unei secţiuni normale oarecare este
exprimată în funcţie de azimutul său A şi. în cazul elipsoidului de rotaţie. de curburile secţiunii
meridianului (M) şi. respectiv. a primului vertical (N).
Raza de curbură medie Gauss
8
Media aritmetică a razelor de curbură ale secţiunilor normale care trec printr-un punct
situat pe elipsoid, atunci când numărul acestor secţiuni tinde către infinit, se numeşte rază medie
de curbură sau rază medie Gauss, notată R:
. (3.32)
Pentru a realiza media aritmetică din relaţia (3.32) utilizăm următoarea schiţă ajutătoare:
A = 0g
Ai
A
A = 100g
S
. (3.33)
. (3.34)
lor pe elipsoid
Din figura 3.5, cu ajutorul mărimilor determinate mai înainte, se pot determina:
coordonatele punctului St (Xt Yt Zt) situat pe suprafaţa fizică. Acestea intervin în
operaţiile specifice geodeziei cu sateliţi;
coordonatele punctului S (X. Y. Z) proiectat pe elipsoidul de referinţă în lungul
normalei la elipsoid.
Reamintind că semnificaţia geometrică a razei de curbură a primului vertical este N =
SS0, rezultă cu uşurinţă următoarele relaţii de calcul:
a. Pentru puncte proiectate pe elipsoid (S) rezultă prin utilizarea noilor notaţii
introduse:
9
(3.35)
b. Pentru puncte situate pe suprafaţa fizică (St) rezultă analog:
(3.36)
Zz
S(B)
S(3)
rx
90+B
z
O B x
S
S0
10
Pentru o suprafaţă oarecare coordonatele carteziene ale unui punct sunt funcţie de
coordonatele geografice elipsoidice:
, , (3.37)
În cazul elipsoidului de rotaţie coordonata z este funcţie numai de latitudinea punctului
.
Elementul liniar al unei curbe oarecare situate pe o suprafaţă poate fi determinat printr-o
relaţie de forma:
(3.38)
Pentru exprimarea elementului liniar al curbei în funcţie de coordonatele geografice
elipsoidice şi de cele carteziene se va diferenţia dx, dy şi dz, obţinând:
(3.38)
(3.40)
; ; (3.47)
ds = elementul de arc
Ţinând cont de expresiile lui dx, dy şi dz se poate scrie:
; ; (3.48)
; ; (3.49)
12
; ; (3.50)
Notând cu unghiul dintre liniile de coordonate generale şi , se poate
scrie: (3.51)
Înlocuind expresiile cosinuşilor directori se obţine:
; (3.52)
(3.53)
(3.55)
(3.56)
În cazul unui sistem ortogonal de coordonate este îndeplinită condiţia F=0, adică
, iar pentru elipsoidul de rotaţie terestru: şi
(3.56)
13
Dacă suprafaţa terestră se aproximează cu o sferă de rază medie Gauss, elementul de
arie devine:
(3.57)
Pentru suprafeţe mici se aplică relaţia:
(3.58)
; ; ;
Înlocuind în relaţia cosA şi ţinând cont de expresiile coeficienţilor E, F şi G se va
obţine:
(3.60)
14
(3.61)
şi (3.62)
Atunci când se consideră un domeniu infinit mic, lungimea unui arc de meridian şi de
paralel se determină cu relaţiile:
şi rezultă că:
şi (3.63)
(3.64)
Din relaţia tangentei rezultă o relaţie foarte importantă pentru teorema lui Clairaut.
15