Disciplina: GEODEZIE SATELITARĂ Proiect

Pagina 1
DISCIPLINA: GEODEZIE SATELITARĂ

Student: Nume Prenume

Facultate.........................
An IV................................
Grupa..............................
Îndrumător proiect:........
Pagina 2
Tema proiectului:

Etapa 1 - Proiectarea măsurătorilor GPS

Date de intrare:
2 mașini, 5 GPS-uri, 5 operatori și 2 soferi, timp de staționare simultan într-o sesiune: 2 ore.
Se va considera de asemeni că distanța între 2 puncte este egală cu distanța din
coordonate + N în m, iar viteza medie de mers de la un punct la altul este de 60 km/oră.
Date de ieșire:
2 variante de proiect în care se va determina:
 lungimea totală parcursă de mijloacele auto,
 timpul total necesar măsurătorilor
 numărul de sesiuni necesare pentru finalizarea acestora.
Între cele 2 variante se va indica varianta optimă și se va justifica alegerea.

Etapa 2 - Transformarea coordonatelor geodezice (B,L,h) în coordonate geocentrice (X, Y, Z)

Etapa 3 - Unificarea determinărilor GPS într-o singură rețea
Etapa 4. - Selecția și rezolvarea triunghiurilor care compun rețeaua astfel încât să nu existe triunghiuri
care se intersectează.

Date de intrare: Vectorii 3D

Date de ieșire: Unghiurile triunghiurilor selectate și schita rețelei de sprijin

Etapa 5 - Problema geodezică directă

Pornind de la N latura, începând cu cea mai scurtă, singura considerată cunoscută, inclusiv
coordonatele geodezice ale unui punct de capăt, să se determine coordonatele geodezice ale
tuturor punctelor din rețea știind și unghiurile triunghiurilor determinate în etapa 3.

Date de intrare: Cel mai mic al N-lea vector 3D din rețea, coordonatele geodezice ale punctelor ce
formează vectorul și unghiurile triunghiurilor calculate la etapa 3.
Date de iesire: Coordonate geocentrice ale punctelor din rețea.

Etapa 6 - Problema geodezică inversă

Să se determine orientările și lungimile vectorilor care formează triunghiurile rețelei.

Date de intrare: Coordonatele geodezice ale punctelor din rețea

Date de ieșire: Lungimi și orientari ale laturilor triunghiurilor din rețea.

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

Disciplina: GEODEZIE SATELITARĂ Proiect

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

N - este numărul de ordine al fiecărui student, conform catalog

Pagina 3...n – Descriere metodologie

Etapa 1 - Proiectarea măsurătorilor GPS

Etape de realizat
 Întocmirea schiței rețelei și evidențierea sesiunilor, a traseelor si a echipamentelor utilizate.
Se vor întocmi 2 tabele cu ordinea punctelor de măsurat, a timpilor de mers și de măsurare, a
km-ilor parcurși și se vor compara și justifica parametri celor două soluții și alegerea solușiei
optime.

Etapa 2 - Transformarea coordonatelor geodezice (B,L,h) în coordonate geocentrice (X, Y, Z)

Descrierea procedurilor, formule folosite: Formula de transformare în grade cu zecimale, formula de

transformare în radiani, formule de transformare coordonate:
 Transformarea coordonatelor B L în B L cu zecimale
B transformat= Grade+(Minute +secunde/60)/60
L transformat= Grade+(Minute +secunde/60)/60
B- Latitudine; L- Longitudine
a cos B cos L c cos B cos L
X  ;
W V
a cos B sin L c cos B sin L
Y  ;
W V
a(1  e 2 ) sin B c(1  e 2 ) sin B
Z  .
W V
X t  ( N  H e ) cos B cos L;
Yt  ( N  H e ) cos B sin L;
 
Z t  N (1  e 2 )  H e sin B.

Procedură de lucru și calcul - exemplu

Date inițiale: Sesiunea 1
Nr. Cota
Nume grade minute secunde grade minute secunde
Crt. elipsoidală
1 CETATUIA_N-V 43 46 9.63288 25 43 39.61183 83.5321

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

Disciplina: GEODEZIE SATELITARĂ Proiect

2 HODIVOAIA 43 55 4.37604 25 46 46.6784 137.8064

3 MOVILA_BURNAS 43 49 52.8463 25 41 1.79537 119.3112
4 MOVILA_TIGANULUI 43 48 32.77372 25 48 7.39265 72.4333
5 PADUREA_BUCIUMENI 44 2 54.5071 25 53 18.23237 132.2562
6 SLOBOZIA_NORD-EST 43 52 15.90413 25 55 21.46693 73.5762
7 TOPORU_SUD 43 59 56.08529 25 38 59.86238 128.9227

Transformarea coordonatelor punctelor măsurate în sesiunea 1 în grade cu zecimale

Nr. Crt. Nume punct Latitudine Longitudine Cota elipsoidala
1 CETATUIA_N-V 43.769342 25.72766995 83.5321
2 HODIVOAIA 43.917882 25.77963289 137.8064
3 MOVILA_BURNAS 43.831346 25.68383205 119.3112
4 MOVILA_TIGANULUI 43.809104 25.80205351 72.4333
5 PADUREA_BUCIUMENI 44.048474 25.88839788 132.2562
6 SLOBOZIA_NORD-EST 43.871084 25.92263036 73.5762
7 TOPORU_SUD 43.998913 25.64996177 128.9227

Transformarea coordonatelor geodezice în coordonate geocentrice – Sesiunea 1

Nr. Crt. Nume punct W X Y Z
1 CETATUIA_N-V 0.998397 4155974.213 2002607.813 4389678.114
2 HODIVOAIA 0.998388 4143896.892 2001421.054 4401619.242
3 MOVILA_BURNAS 0.998393 4153231.242 1997371.983 4394675.230
4 MOVILA_TIGANULUI 0.998395 4150611.235 2006667.396 4392859.647
5 PADUREA_BUCIUMENI 0.998381 4131020.910 2004882.333 4412056.418
6 SLOBOZIA_NORD-EST 0.998391 4142089.713 2013312.850 4397827.645
7 TOPORU_SUD 0.998384 4142776.266 1989329.426 4408094.252

Notă: Identic și pentru sesiunea 2 și 3

Etapa 3 - Unificarea determinărilor GPS într-o singură rețea

Descrierea procedurii, justificarea metodei folosite

Coordonate inițiale
Sesiunea 1 X Y Z
CETATUIA_N-V 4155974.213 2002607.813 4389678.115
HODIVOAIA 4143896.892 2001421.054 4401619.242
MOVILA_BURNAS 4153231.243 1997371.983 4394675.230
MOVILA_TIGANULUI 4150611.235 2006667.396 4392859.647
PADUREA_BUCIUMENI 4131020.911 2004882.333 4412056.418
SLOBOZIA_NORD-EST 4142089.736 2013312.802 4397827.645
TOPORU_SUD 4142776.266 1989329.426 4408094.253
Sesiunea 2 X Y Z
BUDENI_VEST 4113692.266 2015082.139 4423470.483
CUCURUZU_2 4134139.292 1994324.610 4413906.356
LA_CORNET 4103452.270 2031135.969 4425607.563
MOVILA_BRANISTE 4131212.128 2021177.785 4404404.521
PADUREA_BUCIUMENI 4131020.913 2004882.338 4412056.422
SLOBOZIA_NORD-VEST 4142089.749 2013312.804 4397827.646

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

Disciplina: GEODEZIE SATELITARĂ Proiect

VLAD_TEPES 4117463.570 2017461.311 4418938.390

Sesiunea 3 X Y Z
BUDENI_VEST(GR56) 4113690.226 2015081.427 4423468.757
CUCURUZU(GR57) 4108423.588 1976164.694 4445790.215
GHIMPATI(GR55) 4114418.112 1990712.298 4433784.141
MIRSA(GR52) 4120718.178 1971711.834 4436438.423
PODU_DOAMNEI(GR53) 4118515.386 1984164.984 4432927.147
STOENESTI(GR51) 4134137.641 1994324.030 4413905.145
TUFANILOR(GR54) 4122967.378 1992565.050 4425060.467

Diferente sesiunea 2 si 3
2.0402 0.7125 1.7266

Diferente sesiunea 1 si 2 Transformat

2.0378 0.7068 1.7230
2.0270 0.7107 1.7262
Medie 2.0324 0.70875 1.7246

Nume punct X Y Z
PADUREA_BUCIUMENI -0.0024 -0.0057 -0.0036
S1 - S2 SLOBOZIA_NORD-EST -0.0132 -0.0018 -0.0004
Media -0.0078 -0.0038 -0.0020

Nume punct x y z
S2-S3
BUDENI_VEST 2.0324 0.7087 1.7246

Coordonatele X,Y,Z unificate într-o singură rețea

Nr. Crt. Nume punct X Y Z

1 CETATUIA_N-V 4155974,2126 2002607,8127 4389678,1148
2 HODIVOAIA 4143896,8919 2001421,0543 4401619,2417
3 MOVILA_BURNAS 4153231,2425 1997371,9827 4394675,2304
S1
4 MOVILA_TIGANULUI 4150611,2350 2006667,3958 4392859,6468
5 PADUREA_BUCIUMENI 4131020,9105 2004882,3327 4412056,4182
6 SLOBOZIA_NORD-EST 4142089,7358 2013312,8023 4397827,6451
7 TOPORU_SUD 4142776,2660 1989329,4260 4408094,2525

Nr. Crt. Nume punct X Y Z

1 BUDENI_VEST 4113692,2581 2015082,1355 4423470,4811
2 CUCURUZU_2 4134139,2837 1994324,6058 4413906,3536
3 LA_CORNET 4103452,2626 2031135,9652 4425607,5613
S2
4 MOVILA_BRANISTE 4131212,1199 2021177,7815 4404404,5188
5 PADUREA_BUCIUMENI 4131020,9051 2004882,3347 4412056,4198
6 SLOBOZIA_NORD-VEST 4142089,7412 2013312,8004 4397827,6435
7 VLAD_TEPES 4117463,5624 2017461,3071 4418938,3879

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

Disciplina: GEODEZIE SATELITARĂ Proiect

Nr. Crt. Nume punct X Y Z

1 BUDENI_VEST 4113692,2581 2015082,1355 4423470,4811
2 CUCURUZU 4108425,6203 1976165,4031 4445791,9391
3 GHIMPATI 4114420,1440 1990713,0072 4433785,8656
S3
4 MIRSA 4120720,2105 1971712,5432 4436440,1471
5 PODU_DOAMNEI 4118517,4186 1984165,6926 4432928,8715
6 STOENESTI 4134139,6733 1994324,7384 4413906,8699
7 TUFANILOR 4122969,4104 1992565,7585 4425062,1919

Etapa 4. - Selecția și rezolvarea triunghiurilor care compun rețeaua astfel încât să nu existe triunghiuri
care se intersectează

Descrierea procedurii, formule de calcul, justificare tehnică, schița rețelei

1 1-2
2 1-17
3 1-4
4 2-4

 Schița rețelei se va realiza în Autocad

 Plecând de la valorile date ale vectorilor rețelei, vom calcula unghiurile triunghiurilor formate
de acestia cu ajutorul Teoremei Cosinusului.
 Se vor ordona vectorii după lungime de la cel mai scurt la cel mai lung, se va alege un vector de
bază

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

 Disciplina: GEODEZIE SATELITARĂ Proiect

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()

Se va înmulți cu 180/PI() deoarece valorile din excel sunt în radiani

Se vor reduce la orizont distanțele date cu formula :

d = D*cos V;

Etapa 5. Problema geodezică directă

Pornind de la latura cea mai scurtă, singura considerată cunoscută, inclusiv coordonatele
geodezice ale punctelor de capăt, să se determine coordonatele geodezice ale tuturor
punctelor din rețea știind și unghiurile triunghiurilor determinate în etapa 3.

Date de intrare: Cel mai mic vector 3D din rețea, coordonatele geodezice ale punctelor ce
formează vectorul și unghiurile triunghiurilor calculate la etapa 3.
Date de ieșire: Coordonate geocentrice ale punctelor din rețea

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,88839793
HODIVOAIA

Formul e Folosite
a) Transformarea coordonatelor B L în B L cu zecimale
B transformat=(((Grade*60+Minute)*60+Secunde)/3600)
L transformat=(((Grade*60+Minute)*60+Secunde)/3600)
b) Transformarea distanțelor
s = (S*180) / (PI()*Ra)
unde
s-distanța în grade
S-distanța în metri
Ra-Raza de curbură a unei laturi de azimut A

Ra = (M*N) / (N*cos2A+M*sin2A)
M= (a*(1-e2)) / W3
N= a / W

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

Disciplina: GEODEZIE SATELITARĂ Proiect

c) Calcule Azimut, Latitudine și Diferența de Longitudine

 90 o  1  s
 A  A1
sin
tg 2 2
 tg  2
 2 2 90  1  s
o
sin
 2

 90  1  s
o
cos
 A2  2 A1 2
tg 2  tg 
2 90  1  s
o
 cos
 2

Unde
A2=azimutul invers;  2 =Diferenta de longitudine
A1=azimutul cunoscut; 1 =Latitudinea cunoscută
s=distanța cunoscută transformată în 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
A  2
cos 2
90  2
o
2 90o  1  s
tg   tg
2 A  2 2
cos 2
2
Cu una din formulele de mai sus determinăm latitudinea punctului nou
 2 -latitudinea care vrem să o determinam
Celelalte elemente fiind déjà cunoscute

Etapa 6 - Problema geodezică inversă

Să se determine orientările și lungimile vectorilor care formează triunghiurile rețelei

Date de intrare: Coordonatele geodezice ale punctelor din rețea
Date de iesire: Lungimi și orientări ale laturilor triunghiurilor din rețea.

Plecând de la coordonatele geodezice ale punctelor din rețea calculate anterior s-au determinat
Lungimile și Orientările laturilor triunghiurilor din rețea

VECTORUL Latitudine Longitudine

5 -7 44.04847419 25.88839793

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ

Disciplina: GEODEZIE SATELITARĂ Proiect

Formule Folosite:
Se va calcula orientarea și orientarea inversă 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 Inversă
 2 = Latitudinea punctului 2 în formula de mai sus
1 = Latitudinea punctului 1 în formula de mai sus

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

Ambele formule de mai sus rezolvă aceiași problemă și anume:

Calculul distanței dintre cele 2 puncte de coordonate cunoscute (B,L) rezultatul fiind o
distanță ca unghi.

Cu ajutorul formulei de mai jos se va transforma ulterior distanța ca unghi în distanță în

metri.

S = (2*PI()*Ra*s) / (360)
unde
s-distanța în grade
S-distanța în metri
Ra-Raza de curbură a

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

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

N= a / W

Coordonator proiect: Conf. Dr. Ing. Nicolae Ion BĂBUCĂ