UNIVERSITATEA DE STIINTE AGRONOMICE SI MEDICINA

VETERINARA BUCURESTI, FACULTATEA DE IMBUNATATIRI

FUNCIARE SI INGINERIA MEDIULUI
Specializarea masuratori terestre si cadastru

PROIECT GEODEZIE
ANUL III

Student: TUDORACHE ALINA

Grupa 9302
Nr. ordine 22
 Tema proiect

Etapa 1. Transformarea coordonatelor geodezice (B,L,h) in coordonate geocentrice (X,

Y, Z)

Etapa 2. Unificarea determinarilor GPS intr-o singura retea

Etapa 3. Selectia si rezolvarea triunghiurilor care compun reteaua astfel incat sa nu

existe triunghiuri care se intersecteaza.
Se deseneaza triunghiuri prin puncte doar folosind vectorii existenti dar fara a se
suprapune. Vectorii neutilizati se vor sterge din lista.
Date de intrare: Vectorii 3D
Date de iesire: Unghiurile triunghiurilor selectate si schita retelei

Etapa 4. Problema geodezica directa

Primul pas este determinarea traseului de calcul (cu rosu in imaginea de mai sus). Pentru
asta din toti vectorii se selecteaza vectorii care compun triunghiurile retelei si se
ordoneaza crescator, vectorul cu numar de ordine egal cu NRO (NumaR de Ordine)
devine vectorul de baza. De la acest vector se determina un traseu de parcurgere a retelei
astfel incat sa trec o singura data prin fiecare punct si sa ajung inapoi in acelasi punct de
plecare. Pornind de la latura de bază, singura considerata cunoscuta, inclusiv
coordonatele geodezice ale unui punct de capat, sa se determine coordonatele geodezice
ale tuturor punctelor din retea folosind formulele dela problema geodezica directa (se
folosesc orientarile si lungimile reduse la orizont
Date de intrare: Vectorul de bază, coordonatele geodezice ale punctului de capat, traseul
de calcul, lungimile si orientarile vectorilor din traseul de calcul
Date de iesire: Coordonate geodezice ale punctelor din retea

Etapa 5 Problema geodezica inversa

Sa se determine orientarile si lungimile vectorilor care formeaza triunghiurile retelei
Date de intrare: Coordonatele geodezice ale punctelor din retea
Date de iesire: Lungimi si orientari ale laturilor triunghiurilor din retea.
Parametri Elipsoidul Krasovski Elipsoidul WGS-84
a 6378245.000 m 6378137.000 m
f 0.003352329 0.003352810
b 6356863.024m 6356752.314 m
e2 0.006693422 0.006694380
e 2 0.006738525 0.006739497
E 521825.488 m 521854.012 m
c 6399698.902 m 6399593.626 m

Din fisierele initiale de date in format text (*.rep) se extrag datele in 3 fisiere, un fisier
[Link], un alt fisier [Link] si un alt fisier [Link], apoi acestea se citesc in Excel.

FISIER BLH-DATE EXTRASE

grade latit B in
nr crt. nume grade minute secunde zecimale
1 Cetatuia N-v 43 46 9.63288 43.76934247
2 Hodivoaia 43 55 4.37604 43.91788223
3 Movila Burnas 43 49 52.8463 43.83134619
4 Movila Tiganului 43 48 32.77372 43.80910381
Padurea
5 Buciumeni 44 2 54.5071 44.04847419
6 Slobozia N-v 43 52 15.90413 43.87108448
7 Toporu Sud 43 59 56.08529 43.99891258

Etapa 1. Transformarea coordonatelor geodezice in coordonate geocentrice

Avem date de intrare cunoscute: coordonate geodezice din fiecare sesiune S1, S2,
S3 date in anexa 1,2,3
Trebuie sa determinam coordonatele geocentrice pentru fiecare sesiune
Pentru a transforma coordonatele geodezice BLH in ccordonate geocentrice vom
parcurge urmatoarele etape
a) Cu ajutorul datelor din fiecare sesiune transformam latitudinea si
longitudinea in grade cu zecimale cu ajutorul urmatoarelor formule
B transformat= ((grade*60+minute)*60+secunde)/3600
L transformat= ((grade*60+minute)*60+secunde)/3600
B- Latitudine; L-Longitudine
Astfel ca am extras datele din sesiunea I fiind expimate in grade minute secunde,
urmand ca in tabelele 1 si 2, latitudinea si longitudinea punctelor sa fie
transformate in grade cu zecimale cu ajutorul formulelor de mai sus

SESIUNEA 1

TABEL1 TRANSFORMAREA IN GRADE ZECIMALE A LATITUDINII B

B transformat= ((grade*60+minute)*60+secunde)/3600
Grade, minute, secunde –extrase din sesiunea I

nr crt. nume grade minute secunde grade latit B in zecimale

1 Cetatuia N-v 43 46 9.63288 43.76934247
2 Hodivoaia 43 55 4.37604 43.91788223
3 Movila Burnas 43 49 52.8463 43.83134619
4 Movila Tiganului 43 48 32.77372 43.80910381
5 Padurea Buciumeni 44 2 54.5071 44.04847419
6 Slobozia N-v 43 52 15.90413 43.87108448
7 Toporu Sud 43 59 56.08529 43.99891258

TABEL 2- TRANSFORMAREA IN GRADE ZECIMALE A

LONGITUDINII L
L transformat= ((grade*60+minute)*60+secunde)/3600
Grade, minute, secunde –extrase din sesiunea I

grade long L in
nr crt. nume grade minute secunde zecimale
1 Cetatuia N-v 25 43 39.61183 25.72766995
2 Hodivoaia 25 46 46.6784 25.77963289
3 Movila Burnas 25 41 1.79537 25.68383205
4 Movila Tiganului 25 48 7.39265 25.80205351
5 Padurea Buciumeni 25 53 18.23237 25.88839788
6 Slobozia N-v 25 55 21.46693 25.9226297
7 Toporu Sud 25 38 59.86238 25.64996177
 b) Transformarea in radiani dupa formula
Y=(pi*X)/180
Unde Pi=3.14159
X-grade calculate
Y- coresponendul radiani

TABEL 3 TRANSFORMAREA IN RADIANI A LATITUDINII B dupa

formula Y=(pi*X)/180

transformare in radiani
nr. crt nume grade cu zecimale latitudine
1 Cetatuia N-v 43.76934247 0.763918492
2 Hodivoaia 43.91788223 0.766510998
3 Movila Burnas 43.83134619 0.76500066
4 Movila Tiganului 43.80910381 0.764612458
5 Padurea Buciumeni 44.04847419 0.768790256
6 Slobozia N-v 43.87108448 0.765694224
7 Toporu Sud 43.99891258 0.767925243

TABEL 4 TRANSFORMAREA IN RADIANI A LONGITUDII L dupa formula

Y=(pi*X)/180

transformare in radiani
nr. crt nume grade cu zecimale longitudine
1 Cetatuia N-v 25.72766995 0.44903217
2 Hodivoaia 25.77963289 0.449939094
3 Movila Burnas 25.68383205 0.448267055
4 Movila Tiganului 25.80205351 0.450330407
5 Padurea Buciumeni 25.88839788 0.451837399
6 Slobozia N-v 25.9226297 0.452434857
7 Toporu Sud 25.64996177 0.447675908

URMEAZA TABELE LANDSCAPE sesiunea 1

Dupa aceleasi etape vom efectua calculele pentru sesiunea 2 si 3

SESIUNEA 2

TABEL1 TRANSFORMAREA IN GRADE ZECIMALE A LATITUDINII B

B transformat= ((grade*60+minute)*60+secunde)/3600
Grade, minute, secunde –extrase din sesiunea I

grade latit B in
nr crt. nume grade minute secunde zecimale
1 Budeni vest 44 11 30,4609 44,19179469
2 Cucuruzu 2 44 4 17,93891 44,0716497
3 La Cornet 44 13 7,41055 44,21872515
4 Movila Braniste 43 57 11,64981 43,95323606
5 Padurea Buciumeni 44 2 54,50708 44,04847419
6 Slobozia N-v 43 52 15,90385 43,8710844
7 Vlad Tepes 44 8 5,05903 44,13473862

TABEL 2- TRANSFORMAREA IN GRADE ZECIMALE A

LONGITUDINII L
L transformat= ((grade*60+minute)*60+secunde)/3600
Grade, minute, secunde –extrase din sesiunea I

grade long L in
nr crt. nume grade minute secunde zecimale
1 Budeni vest 26 5 52,1185 26,09781069
2 Cucuruzu 2 25 45 10,2213 25,75283925
3 La Cornet 26 20 4,53256 26,33459238
4 Movila Braniste 26 4 11,9756 26,06999322
5 Padurea Buciumeni 25 53 18,23255 25,88839793
6 Slobozia N-v 25 55 21,46675 25,92262965
7 Vlad Tepes 26 6 13,6029 26,10377858
TABEL 3 TRANSFORMAREA IN RADIANI A LATITUDINII B dupa
formula Y=(pi*X)/180

transformare in radiani
nr. crt nume grade cu zecimale latitudine
1 Budeni vest 44,19179469 0,771291668
2 Cucuruzu 2 44,0716497 0,769194744
3 La Cornet 44,21872515 0,771761693
4 Movila Braniste 43,95323606 0,767128038
5 Padurea Buciumeni 44,04847419 0,768790256
6 Slobozia N-v 43,8710844 0,765694222
7 Vlad Tepes 44,13473862 0,770295853

TABEL 4 TRANSFORMAREA IN RADIANI A LONGITUDII L dupa formula

Y=(pi*X)/180
transformare in radiani
nr. crt nume grade cu zecimale longitudine
1 Budeni vest 26,09781069 0,455492339
2 Cucuruzu 2 25,75283925 0,449471457
3 La Cornet 26,33459238 0,459624956
4 Movila Braniste 26,06999322 0,455006833
5 Padurea Buciumeni 25,88839793 0,4518374
6 Slobozia N-v 25,92262965 0,452434856
7 Vlad Tepes 26,10377858 0,455596499
 Sesiune 3

TABEL1 TRANSFORMAREA IN GRADE ZECIMALE A LATITUDINII B

B transformat= ((grade*60+minute)*60+secunde)/3600
Grade, minute, secunde –extrase din sesiunea III

nr crt. nume grade minute secunde grade latit B in zecimale

1 Budeni Vest (Gr56) 44 11 30.46925 44.19179701
2 Cucuruzu (Gr 57) 44 28 19.75796 44.47215499
3 Ghimpati (Gr 55) 44 19 16.20276 44.32116743
4 Mirsa (Gr 52) 44 21 15.69538 44.35435983
5 Podu Doamnei (Gr 53) 44 18 37.25556 44.31034877
6 Stoenesti (Gr 51) 44 4 17.94992 44.07165276
7 Tufanilor (Gr 54) 44 12 41.52875 44.21153576

TABEL 2- TRANSFORMAREA IN GRADE ZECIMALE A

LONGITUDINII L
L transformat= ((grade*60+minute)*60+secunde)/3600
Grade, minute, secunde –extrase din sesiunea III

nr crt. nume grade minute secunde grade long L in zecimale

1 Budeni Vest (Gr56) 26 5 52.13011 26.09781392
2 Cucuruzu (Gr 57) 25 41 15.72395 25.6877011
3 Ghimpati (Gr 55) 25 49 10.08004 25.81946668
4 Mirsa (Gr 52) 25 34 13.90281 25.57052856
5 Podu Doamnei (Gr 53) 25 43 23.64477 25.72323466
6 Stoenesti (Gr 51) 25 45 10.23006 25.75284168
7 Tufanilor (Gr 54) 25 47 37.48232 25.79374509

Transformarea in radiani dupa formula

Y=(pi*X)/180
Unde Pi=3.14159
X-grade calculate
Y- coresponendul radiani
TABEL 3 TRANSFORMAREA IN RADIANI A LATITUDINII B dupa
formula Y=(pi*X)/180

nr. crt nume grade cu zecimale transformare in radiani latitudine

1 Budeni Vest (Gr56) 44.19179701 0.771291709
2 Cucuruzu (Gr 57) 44.47215499 0.776184874
3 Ghimpati (Gr 55) 44.32116743 0.773549647
4 Mirsa (Gr 52) 44.35435983 0.774128963
5 Podu Doamnei (Gr 53) 44.31034877 0.773360826
6 Stoenesti (Gr 51) 44.07165276 0.769194798
7 Tufanilor (Gr 54) 44.21153576 0.771636215

TABEL 4 TRANSFORMAREA IN RADIANI A LONGITUDII L dupa formula

Y=(pi*X)/180

grade cu transformare in radiani

nr. crt nume zecimale longitudine
1 Budeni Vest (Gr56) 26.09781392 0.455492396
2 Cucuruzu (Gr 57) 25.6877011 0.448334583
3 Ghimpati (Gr 55) 25.81946668 0.450634324
4 Mirsa (Gr 52) 25.57052856 0.446289538
5 Podu Doamnei (Gr 53) 25.72323466 0.44895476
6 Stoenesti (Gr 51) 25.75284168 0.449471499
7 Tufanilor (Gr 54) 25.79374509 0.450185398
Etapa 2. Unificarea determinarilor GPS intr-o singura retea
Descrierea procedurii:in urma suprapunerii determinarilor GPS, rezulta puncte de
legattura. Se identifica, se face diferenta citirilor, media si repartitia mediei.
SESIUNEA 1
X-corectat Y-corectat Z-corectat
Crt nume

1 Cetatuia N-v 4155974.224 2002607.825 4389677.992

2 Hodivoaia 4143896.911 2001421.068 4401619.124
3 Movila Burnas 4153231.262 1997372.018 4394675.11
4 Movila Tiganului 4150611.234 2006667.435 4392859.529
5 Padurea Buciumeni 4131020.916 2004882.344 4412056.301
6 Slobozia N-E 4142089.752 2013312.846 4397827.528
7 Toporu Sud 4142776.291 1989329.451 4408094.312

SESIUNEA 2
Crt nume X-corectat Y-corectat Z-corectat

1 Budeni vest 4113692.333 2015082.17 4423470.572

2 Cucuruzu 2 4134139.346 1994324.635 4413906.432
3 La Cornet 4103452.348 2031136.006 4425607.593
4 Movila Braniste 4131212.227 2021177.833 4404404.589
Padurea
5 Buciumeni 4131021.025 2004882.392 4412056.483
6 Slobozia N-v 4142089.846 2013312.85 4397827.658
7 Vlad Tepes 4117463.621 2017461.334 4418938.473

SESIUNEA 3
Crt nume X-corectat Y-corectat Z-corectat

1 Budeni Vest (Gr56) 4113690.267 2015081.449 4423468.804

2 Cucuruzu (Gr 57) 4108423.63 1976164.717 4445790.264

3 Ghimpati (Gr 55) 4114418.155 1990712.321 4433784.184

4 Mirsa (Gr 52) 4120718.217 1971711.856 4436438.464

5 Podu Doamnei (Gr 53) 4118515.428 1984165.006 4432927.194

6 Stoenesti (Gr 51) 4134137.681 1994324.051 4413905.194

7 Tufanilor (Gr 54) 4122967.421 1992565.073 4425060.514

si-s2 nume x y z
Padurea Buciumeni 0.109 0.048 0.182
Slobozia n-v 0.094 0.004 0.13
Media
 s2-s3 nume x y z
Budeni Vest 2.023 0.703 1.755

Etapa 2. Calculul Unghiurilor triunghiurilor selectate si schita retelei

- Am realizat in Autocad 2007 schita retelei

- Plecand de la valorile date ale vectorilor retelei am calculat unghiurile

triunghiurilor formate de acestia cu ajutorul Teoremei Cosinusului.

- Mi-am ordonat vectorii dupa lungime de la cel mai scurt la cel mai lung
alegand ca vector de baza al 19-lea.

Cos a = (B2+C2-A2) / 2*B*C rezultand a=arccos((B2+C2-A2) / 2*B*C)*180/PI()

Cos b = (A2+C2-B2) / 2*A*C rezultand b=arccos(A2+C2-B2) / 2*A*C)*180/PI()

Cos c = (A2+B2-C2) / 2*A*B rezultand c=arccos((A2+B2-C2) / 2*A*B)*180/PI()

Am inmultit cu 180/PI() deoarece valorile din excel sunt in radiani

Am redus la orizont distantele date cu formula :

d = D*cos V;

D
V

Etapa 3. Problema geodezica directa

Pornind de la NRO latura incepand cu cea mai scurta, singura considerata

cunoscuta, inclusiv coordonatele geodezice ale punctelor de capat, sa se determine
 coordonatele geodezice ale tuturor punctelor din retea stiind si unghiurile
triunghiurilor determinate in etapa 3
Date de intrare: Cel mai mic al NRO vector 3D din retea, coordonatele geodezice
ale punctelor ce formeaza vectorul si unghiurile triunghiurilor calculate la etapa 3
Date de iesire: Coordonate geocentrice ale punctelor din retea

Azimutul Padurea Buciumeni

Date de intrare Distanta Grade Minute Secunde VECTORUL B L
PADUREA_BUCIUMENI-
16932,3955 211 3 33,38 5-7 44,04847419 25,888
HODIVOAIA

Formul e Folosite
a) Transformarea coordonatelor B L in B L cu zecimale
B transformat=(((Grade*60+Minute)*60+Secunde)/3600)
L transformat=(((Grade*60+Minute)*60+Secunde)/3600)

b) Transformarea distantelor

s = (S*180) / (PI()*Ra) unde

s-distanta in grade
S-distanta in metri
Ra-Raza de curbura a unei laturi de azimut A

Ra = (M*N) / (N*cos2A+M*sin2A)

M= (a*(1-e2)) / W3

N= a / W

c) Calcule Azimut, Latitudine si Diferenta de Longitudine

 90o  1  s
 A  A
sin
tg 2 2
 tg 1  2
 2 2 90 o
 1  s
sin
 2

 90  1  s
o
cos
 A2  2 A 2
tg 2  tg 1 
2 90o  1  s
 cos
 2
 Unde
A2=azimutul invers; 2 =Diferenta de longitudine
A1=azimutul cunoscut; 1 =Latitudinea cunoscuta
s=distanta cunoscuta transformata in valori unghiulare

Rezolvând sistemul de ecuaţii constituit mai sus, rezultă necunoscutele  2 şi

A2 , celelalte elemente fiind cunoscute.

A  2
sin 2
90o  2 2  tg 90  1  s
o
tg 
2 A  2 2
sin 2
2
A2  2
cos
90o  2 2 90o  1  s
tg   tg
2 A2  2 2
cos
2
Cu una din formulele de mai sus determinam latitudinea punctului nou
 2 -latitudinea care vrem sa o determinam
Celelalte elemente fiind déjà cunoscute

Etapa 4 Problema geodezica inversa

Sa se determine orientarile si lungimile vectorilor care formeaza triunghiurile retelei

Date de intrare: Coordonatele geodezice ale punctelor din retea
Date de iesire: lungimi si orientari ale laturilor triunghiurilor din retea.

Plecand de la coordonatele geodezice ale punctelor din retea calculate anterior s-au
determinat Lungimile si Orientarile laturilor triunghiurilor din retea
VECTORUL Latitudine Longitudine

5-7 44,04847419 25,88839793

Formule Folosite:
-Am calculat orientarea si orientarea inversa cu formulele:
    1
 A1  A2 cos 2
 2
tg  tg 
 2 2   
sin 2 1
 2

   1
sin 2
tg A2  A1  tg   2
 2 
2 cos 2  1

 2

Unde A1 = Orientarea
A2 = Orientarea Inversa
 2 = Latitudinea punctului 2 in formula de mai sus
1 = Latitudinea punctului 1 in formula de mai sus

A2  A1 
sin 
s 2  ctg 2   
tg  1

2 sin A2  A1 2 
2 

A2  A1 
cos
s 2 2  1 
tg   ctg
2 cos A2  A1 2 

2 

Ambele formule de mai sus rezolva aceiasi problema si anume:

Calculul distantei dintre cele 2 puncte de coordonate cunoscute (B,L)
rezultatul fiind o distanta ca unghi

Cu ajutorul formulei de mai jos am transformat ulterior distanta ca unghi in distanta in

metri.

S = (2*PI()*Ra*s) / (360) unde

s-distanta in grade
S-distanta in metri
Ra-Raza de curbura a

Ra = (M*N) / (N*cos2A+M*sin2A)

M= (a*(1-e2)) / W3

N= a / W