1. clear all
2. p_elipsoidy_GRS80;
3.  
4. %dane
5. %punkt C
6. a_fi_d = 53;
7. a_fi_m = 38;
8. a_fi_s = 0;
9.  
10. a_la_d = 21;
11. a_la_m = 22;
12. a_la_s = 0;
13.  
14. %punkt D
15. b_fi_d = 53;
16. b_fi_m = 21;
17. b_fi_s = 51.9190;
18.  
19. b_la_d = 21;
20. b_la_m = 9;
21. b_la_s = 25.8088;
22.  
23. %obliczenia dla 1992
24. %zamiana na radiany
25. a_fi = stopnie_d_to_rad(stopnie_to_stopnie_d( a_fi_d, a_fi_m, a_fi_s ));
26. a_la = stopnie_d_to_rad(stopnie_to_stopnie_d( a_la_d, a_la_m, a_la_s ));
27.  
28. b_fi = stopnie_d_to_rad(stopnie_to_stopnie_d( b_fi_d, b_fi_m, b_fi_s ));
29. b_la = stopnie_d_to_rad(stopnie_to_stopnie_d( b_la_d, b_la_m, b_la_s ));
30.  
31. %wyznaczenie poludnika strefowego dla ukladu 1992
32. L_0 = 19;
33. L_0 = stopnie_d_to_rad(L_0);
34.  
35. %wyznaczenie wsp na płaszczyznie gaussa
36. [a_x_gk,a_y_gk] = bl_gk(a_fi,a_la,L_0,a,e2);
37. [b_x_gk,b_y_gk] = bl_gk(b_fi,b_la,L_0,a,e2);
38.  
39. %zamiana na wspołrzędne w układzie 1992
40. [a_x_1992, a_y_1992] = cechowanie_1992(a_x_gk, a_y_gk);
41. [b_x_1992, b_y_1992] = cechowanie_1992(b_x_gk, b_y_gk);
42.  
43. %obliczenie odwrotne
44. [a_x_gk_test, a_y_gk_test] = odcechowanie_1992(a_x_1992, a_y_1992);
45. [b_x_gk_test, b_y_gk_test] = odcechowanie_1992(b_x_1992, b_y_1992);
46.  
47. [a_fi_test, a_la_test] = gk_bl(a_x_gk_test, a_y_gk_test, L_0, a, e2);
48. [b_fi_test, b_la_test] = gk_bl(b_x_gk_test, b_y_gk_test, L_0, a, e2);
49.  
50. %zamiana na stopnie
51. [a_fi_d_test, a_fi_m_test, a_fi_s_test] = rad_to_degres(a_fi_test);
52. [a_la_d_test, a_la_m_test, a_la_s_test] = rad_to_degres(a_la_test);
53.  
54. [b_fi_d_test, b_fi_m_test, b_fi_s_test] = rad_to_degres(b_fi_test);
55. [b_la_d_test, b_la_m_test, b_la_s_test] = rad_to_degres(b_la_test);
56.  
57. %raport po 1992
58. disp([' ']);
59. disp(['________________________________________________________________________________']);
60. disp(['Raport obliczeniowy programu GeoKał produkcji Kałowego Softwaru. /Michał Kałuski']);
61. disp(['Układy GK - przeliczenie współrzędnych geodezyjnych na prostokątne w układzie PUWG1992 oraz PUWG2000']);
62.  
63. disp([' ']);
64. disp(['Zadanie 1']);
65.  
66. disp([' ']);
67. disp(['Dane wejściowe:']);
68.  
69. disp(['c_fi = ',num2str(a_fi_d), ' d ' , num2str(a_fi_m), ' m ', num2str(a_fi_s, 8), ' s']);
70. disp(['c_la = ',num2str(a_la_d), ' d ' , num2str(a_la_m), ' m ', num2str(a_la_s, 8), ' s']);
71. disp([' ']);
72. disp(['d_fi = ',num2str(b_fi_d), ' d ' , num2str(b_fi_m), ' m ', num2str(b_fi_s, 8), ' s']);
73. disp(['d_la = ',num2str(b_la_d), ' d ' , num2str(b_la_m), ' m ', num2str(b_la_s, 8), ' s']);
74.  
75. disp([' ']);
76. disp(['Wyniki 1992:']);
77.  
78. disp(['c_x_gk = ',num2str(a_x_gk),]);
79. disp(['c_y_gk = ',num2str(a_y_gk),]);
80. disp([' ']);
81. disp(['d_x_gk = ',num2str(b_x_gk),]);
82. disp(['d_y_gk = ',num2str(b_y_gk),]);
83. disp([' ']);
84.  
85. disp(['c_x_1992 = ',num2str(a_x_1992),]);
86. disp(['c_y_1992 = ',num2str(a_y_1992),]);
87. disp([' ']);
88. disp(['d_x_1992 = ',num2str(b_x_1992),]);
89. disp(['d_y_1992 = ',num2str(b_y_1992),]);
90.  
91. disp([' ']);
92. disp(['Kontrola:']);
93. disp(['c_fi_odwrotne = ',num2str(a_fi_d_test), ' d ' , num2str(a_fi_m_test), ' m ', num2str(a_fi_s_test, 8), ' s']);
94. disp(['c_la_odwrotne = ',num2str(a_la_d_test), ' d ' , num2str(a_la_m_test), ' m ', num2str(a_la_s_test, 8), ' s']);
95. disp([' ']);
96. disp(['d_fi_odwrotne = ',num2str(b_fi_d_test), ' d ' , num2str(b_fi_m_test), ' m ', num2str(b_fi_s_test, 8), ' s']);
97. disp(['d_la_odwrotne = ',num2str(b_la_d_test), ' d ' , num2str(b_la_m_test), ' m ', num2str(b_la_s_test, 8), ' s']);
98.  
99.  
100. disp([' ']);
101.  
102. %obliczenia dla 2000
103. %poludnik zerowy musi byc w stopniach dziesietnych zeby zadzialal IF
104. a_la_dz = rad_to_stopnie_d(a_la);
105. b_la_dz = rad_to_stopnie_d(b_la);
106.  
107. %wyznaczenie południka osiowego i numeru strefy
108. [ a_nr_strefy, a_L_0 ] = poszukiwanie_strefy_2000( a_la_dz );
109. [ b_nr_strefy, b_L_0 ] = poszukiwanie_strefy_2000( b_la_dz );
110.  
111. % do dalszych funkcji potrzebujemy juz poludnika zerowego ale w radianach
112. a_L_0_d = a_L_0;
113. b_L_0_d = b_L_0;
114. a_L_0 = stopnie_d_to_rad(a_L_0);
115. b_L_0 = stopnie_d_to_rad(b_L_0);
116.  
117. %wyznaczenie wsp na plaszczyznie gaussa
118. [a_x_gk,a_y_gk] = bl_gk(a_fi,a_la,a_L_0,a,e2);
119. [b_x_gk,b_y_gk] = bl_gk(b_fi,b_la,b_L_0,a,e2);
120.  
121. %zamiana na wspołrzędne w układzie 2000
122. [a_x_2000, a_y_2000] = cechowanie_2000(a_x_gk, a_y_gk, a_nr_strefy);
123. [b_x_2000, b_y_2000] = cechowanie_2000(b_x_gk, b_y_gk, b_nr_strefy);
124.  
125. %obliczenie odwrotne
126. [ a_x_gk_test, a_y_gk_test ] = odcechowanie_2000( a_x_2000, a_y_2000, a_nr_strefy );
127. [ b_x_gk_test, b_y_gk_test ] = odcechowanie_2000( b_x_2000, b_y_2000, b_nr_strefy );
128.  
129. [a_fi_test, a_la_test] = gk_bl(a_x_gk_test, a_y_gk_test, a_L_0, a, e2);
130. [b_fi_test, b_la_test] = gk_bl(b_x_gk_test, b_y_gk_test, b_L_0, a, e2);
131.  
132. %zamiana na stopnie
133. [a_fi_d_test, a_fi_m_test, a_fi_s_test] = rad_to_degres(a_fi_test);
134. [a_la_d_test, a_la_m_test, a_la_s_test] = rad_to_degres(a_la_test);
135.  
136. [b_fi_d_test, b_fi_m_test, b_fi_s_test] = rad_to_degres(b_fi_test);
137. [b_la_d_test, b_la_m_test, b_la_s_test] = rad_to_degres(b_la_test);
138.  
139. %raport po 2000
140. disp(['Wyniki 2000:']);
141.  
142. disp(['poludnik osiowy punktu C = ',num2str(a_L_0_d), ' d']);
143. disp(['numer strefy punktu C = ',num2str(b_nr_strefy)]);
144. disp([' ']);
145. disp(['poludnik osiowy punktu D = ',num2str(b_L_0_d), ' d']);
146. disp(['numer strefy punktu D = ',num2str(a_nr_strefy)]);
147. disp([' ']);
148.  
149. disp(['c_x_gk = ',num2str(a_x_gk),]);
150. disp(['c_y_gk = ',num2str(a_y_gk),]);
151. disp([' ']);
152. disp(['d_x_gk = ',num2str(b_x_gk),]);
153. disp(['d_y_gk = ',num2str(b_y_gk),]);
154. disp([' ']);
155.  
156. disp(['c_x_2000 = ',num2str(a_x_2000),]);
157. disp(['c_y_2000 = ',num2str(a_y_2000),]);
158. disp([' ']);
159. disp(['d_x_2000 = ',num2str(b_x_2000),]);
160. disp(['d_y_2000 = ',num2str(b_y_2000),]);
161.  
162. disp([' ']);
163. disp(['Kontrola:']);
164. disp(['c_fi_odwrotne = ',num2str(a_fi_d_test), ' d ' , num2str(a_fi_m_test), ' m ', num2str(a_fi_s_test, 8), ' s']);
165. disp(['c_la_odwrotne = ',num2str(a_la_d_test), ' d ' , num2str(a_la_m_test), ' m ', num2str(a_la_s_test, 8), ' s']);
166. disp([' ']);
167. disp(['d_fi_odwrotne = ',num2str(b_fi_d_test), ' d ' , num2str(b_fi_m_test), ' m ', num2str(b_fi_s_test, 8), ' s']);
168. disp(['d_la_odwrotne = ',num2str(b_la_d_test), ' d ' , num2str(b_la_m_test), ' m ', num2str(b_la_s_test, 8), ' s']);
169.  
170. %zadanie 2
171. %wyznaczenie sasiedniego pasa odwzorowawczego
172. if b_nr_strefy == 5
173.     b_nr_strefy_sasiedni = 6;
174. elseif a_nr_strefy == 8
175.     b_nr_strefy_sasiedni = 7;
176. else
177.     b_nr_strefy_sasiedni = b_nr_strefy + 1;
178. end;
179.  
180. if b_nr_strefy_sasiedni == 5
181.     b_L_0 = 15;
182. elseif b_nr_strefy_sasiedni == 6
183.     b_L_0 = 18;
184. elseif b_nr_strefy_sasiedni == 7
185.     b_L_0 = 21;
186. elseif b_nr_strefy_sasiedni == 8
187.     b_L_0 = 24;
188. end;
189.  
190. b_L_0_d = b_L_0;
191. b_L_0 = stopnie_d_to_rad(b_L_0);
192.  
193. %wyznaczenie wsp na plaszczyznie gaussa
194. [b_x_gk,b_y_gk] = bl_gk(b_fi,b_la,b_L_0,a,e2);
195.  
196. %zamiana na wspołrzędne w układzie 2000
197. [b_x_2000, b_y_2000] = cechowanie_2000(b_x_gk, b_y_gk, b_nr_strefy_sasiedni);
198.  
199. %raport końcowy
200. disp([' ']);
201. disp(['Zadanie 2']);
202. disp([' ']);
203.  
204. disp(['Wyniki:']);
205. disp(['nowy południk osiowy = ', num2str(b_L_0_d)]);
206. disp(['nowa strefa = ', num2str(b_nr_strefy_sasiedni)]);
207. disp(['d_x_2000 w pasie sasiednim = ',num2str(b_x_2000),]);
208. disp(['d_y_2000 w pasie sasiednim = ',num2str(b_y_2000),]);
209. disp([' ']);
210. disp(['Współrzędne w układzie 1992 z oczywistych względów są takie same jak w Zadaniu 1']);
211.  
212. disp([' ']);
213. disp(['________________________________________________________________________________']);
214.