Geodezie
Probleme geodezice de bazaPROBLEME GEODEZICE DE BAZA1.Scopul si importanta problemelor geodezice de bazaPe baza observatiilor executate pe suprafata terestrǎ, proiectate pe elipsoidul de referintǎ, urmeaza sa stabilim prin calcul pozitiile punctelor geodezice de prim ordin intr-un anumit sistem de coordonate utilizat pe elipsoid. Pentru ca sistemul coordonatelor geodezice prezinta avantaje deosebite comparativ cu celelalte sisteme, privind calculele necesare si modul unitar de determinare a punctelor geodezice, indiferent de pozitia lor pe suprafata elipsoidului, in practica geodezica este utilizat in mod curent. Problemele geodezice de baza se refera la sistemul coordonatelor geografice elipsoidale si constau in urmatoarele : a. Problema geodezica directa, determina coordonatele geodezice B2, L2 ale unui punct , S2 (fig.1) si azimutul A2 in acest punct in functie de coordonatele geodezice ale punctului S1(B1, L1) si elementele geodezice de legatura (distanta s12 intre punctele S1 si S2 si azimutul A1 in punctul P1). Aceasta problema este similarǎ cu problema cunoscuta la topografie, privind determinarea coordonatelor (x2, y2) ale unui punct S2 in functie de coordonatele (x1, y1) ale punctului S1,distanta s12 si orientarea θ12. Problema geodezica directa, dat fiind scopul care-l urmareste, mai poarta denumirea si de problema transportului de coordonate geodezice. b. Problema geodezica inversa, utilizeaza ca marimi coordonatele geodezice, B1, L1 si B2, L2 ale punctelor S1 si S2 cu ajutorul cǎrora se determinǎ azimutele A1,A2 si lungimea liniei geodezice s12. Cu problema geodezica directa se calculeazǎ coordonatele unor puncte geodezice ce formeaza o retea geodezica de ordinul I, iar cu problema geodezica inversa se calculeazǎ elementele geodezice initiale (distante si azimute) si se verificǎ calculele efectuate la problema geodezica directa. 2.Precizia de calcul a coordonatelor geodezicePrecizia de calcul, ca in orice problema geodezica, urmareste ca erorile de calcul sa fie de cca. 10 ori mai mici fata de erorile medii ale marimilor masurate. Evident, aceasta este in functie si de distantele intre punctele geodezice, respectiv cu cat distantele cresc precizia calculelor are o semnificatie deosebita. Se poate arata ca, in triangulatia de ordinul I in care intervin distante geodezice, (in medie de 30 4o km) este necesar ca aproximatia de calcul pentru coordonatele geodezice B si L sa fie de ± 0 , pentru azimute de ± 0,001, iar pentru distante de ± 0,001m. De asemenea din cauza distantelor geodezice relativ mici, la rezolvarea problemelor geodezice se aplica metode in care se accepta unele aproximatii cum ar fi: dezvoltari in serie, inlocuirea suprafetei elipsoidului cu sfera de raza medie, etc. Rezolvarea problemei geodezice directeIn functie de lungimea liniei geodezice se pot distinge urmǎtoarele metode de rezolvare a problemelor geodezice directe: metoda dezvoltǎrii in serie (metoda Legendre), pentru distante mici s<60 km; metoda argumentelor medii(metoda Gauss), pentru distante medii 60<s<600km; metoda Bessel, pentru distante geodezice mari <600<20000 km. De mentionat, cǎ in prezenta lucrare se vor trata numai metodele de rezolvare a problemelor geodezice pentru distante scurte (s<60 km), urmand ca persoanele interesate sǎ regǎseascǎ celelalte metode in [10] 12], [13], [24] etc. 1.Metoda dezvoltarilor in serieFiind cunoscuta pozitia punctului initial S1(B1,L1), lungimea liniei geodezice s si azimutul A1 al directiei S1S2 (fig.1), se pune problema determinǎrii coordonatelor geodezice B2 si L2, ale punctului S2 si azimutul invers A2. Deoarece diferentele (B2 B1), (L2 L1) si (A2 A1 vor fi functie numai de lungimea liniei geodezice s se acceptǎ urmǎtoarele dezvoltǎri in serie Maclaurain:
(1)
Fig.1.Problema geodezicǎ directǎ Derivatele din egalitatile (1.) se calculeaza in punctul S1, deci cu marimile cunoscute in acest punct. Pentru calculul derivatelor de ordinul I, ne vom referi la fig.2, in care ds este un element liniar al liniei geodezice S1S2. Elementele cunoscute in punctul S1 (latitudine, longitudine si azimut) sunt notate cu B, L, A, urmand ca expresiile derivatelor de ordinal I sǎ fie stabilite in mod geometric. Astfel din triunghiurile S1SS si OSS1 rezultǎ derivata latitudinii B in raport cu lungimea infinitezimalǎ a liniei geodezice : (2)
Fig.2. Metoda Legendre Din triunghiurile S1SS si
(3) Din triunghiurile S1SS si PSS rezultǎ derivata azimutului (A) in raport cu ds: (4)
Recapituland, expresiile derivatelor de ordinul intai se prezintǎ sub forma:
(5)
Inlocuind in (5) valorile razelor de curburǎ, conform relatiilor (l.34) si (1.36) se obtine:
(6)
Derivatele de ordin superior ordinului intai se obtin derivand succesiv in raport cu s expresiile (6) .Astfel, derivatele de ordinul doi se prezintǎ sub forma : (7)
In urma prelucrǎrii relatiilor (7.) rezulta :
(8)
in care s-a notat : (9) Pentru celelalte derivate se procedeaza analog. Revenind la relatia (l.), in care se inlocuiesc expresiile derivatelor se obtin relatiile definitive : (10) (11) (12) in care :
(13) Relatiile (10. (11.) si (12.) pot fi scrise mai simplu notand coeficientii marimilor u , v , u2, v2, cu a, b, c . ,de indici diferiti. Pentru elipsoidul de referinta Krasovski acesti coeficienti se extrag din tabele in functie de latitudinea B1 a punctului P1. 2.Metoda inlocuirii suprafetei elipsoidului cu sfera GaussConsideram punctele P1 si P2 (fig.a) pe o sfera de raza medie a triunghiului sferic P1PP2 .Se cunosc coordonatele geodezice ale punctului P1( B1, L1=0), azimutul directiei P1 P2 (A1) si lungimea liniei geodezice s. Se vor determina coordonatele punctului P2 (B2, L2) si azimutul inversA2. Pentru a usura scrierea relatiilor se alege un triunghi sferic a cǎror elemente sunt reprezentate in fig. b. Pentru rezolvarea problemei scriem formulele a trei elemente adiacente intr-un triunghi sferic oarecare (fig.3 b) si obtinem :
(13)
Fig. a.Triunghiul sferic pe sferǎ; b. Notatiile in triunghiul sferic. Triunghiul sferic ABC identificat cu triunghiul sferic P1P2P conduce la egalitatile : A = A1 a = 900 B2 B = 1800-A1 b = 900 B1 C = L2 c = s si cu ajutorul acestora :
(15) Introducand egalitatile (l5.) in (1) se obtine I II III IV Prin impartirea relatiilor I si II; III si IV se obtine:
(17) Rezolvand sistemul de ecuatii (17) se obtin necunoscutele si L2 in functie de coordonatele geodezice ale punctului P1 si lungimea geodezicǎ s. Azimutul invers A2=+200g. Pentru a afla latitudinea geodezicǎ a punctului P2 (B2) se impart relatiie I si III ; II si IV: (18)
Cele doua relatii (18.) determina aceeasi marime si anume latitudinea B2 a punctului P2. Este de precizat ca metoda descrisa conduce la aceleasi valori cu cele obtinute prin utilizarea relatiilor (10.), (11.), (12.), dar transformate, in sensul ca semiaxele suprafetei elipsoidului vor fi egale, deci a = b si in consecinta : excentricitatea e = e = 0
deci N =M= R 4. Rezolvarea problemei geodezice inverse.Metoda inlocuirii suprafetei elipsoidului cu sfera GaussIn limita aproximatiilor admisibile, pentru distante mici, se poate asimila suprafata elipsoidului cu suprafata unei sfere de razǎ medie Gauss. Ne vom referi la figurile 4 a si 4 b , pe baza carora putem, in prima etapa, scrie formulele semisumei si semidiferentei sinusului si cosinusului in triunghiurile sferice:
(20)
Prin identificarea triunghiurilor P1PP2 si CAB rezulta : A = L2 a = s B = 180- b = 90 B1 C = A1 c = 90 B2
Fig.4.a.Triunghiul pe sferǎ; b. Triunghiul sferic corespondent si cu ajutorul acestora se calculeazǎ:
(22)
Inlocuind egalitatile (22.) in (20.) se obtine : I II (23) III IV Impartind egalitatile I la II si III la IV obtinem :
(24) In relatiile (24.), coordonatele punctelor P1 si P2 fiind date, rezulta ca B1, B2 siL2 sunt cunoscute; deci prin rezolvarea sistemului de ecuatii pot fi determinate azimutele A1 si . Impartind egalitatile I la III , si II la IV se obtine: (25)
Una din relatiile (25) se utilizeazǎ pentru calculul lungimii liniei geodezice dintre punctele P1 si P2, iar cealaltǎ pentru verificare.
|