Teoretická geodézie 2 - úloha 1: Porovnání verzí

Aktuální verze z 30. 10. 2013, 07:52

Název úlohy

GPS - zpracování kódových měření

Zadání úlohy

Například v adresáři pub/igs/data/2013/142/ protokolu FTP na igs.ensg.ign.fr si stáhněte dle numerického zadání soubor ve formátu RINEX obsahující data měřená přijímačem GNSS stanice sítě EUREF EPN v rozmezí jednoho dne v GPS týdnu 1741 (povšimněte si významného výročí spadajícího na tento den). Na soubory je použita jak komprese pomocí SW compress (.Z), tak komprese Hatanaka (.RRd místo .RRo). WinZip formát .Z rozbalit neumí, ve Windows použijte program compress.exe (návod, decompr.bat), nebo gzip.exe. O kompresi Hatanaka se více dočtete zde, kde je také odkaz na SW na dekompresi pro různé OS (i pro Windows). Kompresi .tar umí rozbalit např. TotalCommander. Pokud se vám zdá, že váš rinex neodpovídá formátu souboru rinex v níže uvedeném manuálu rinex300.pdf, patrně jste neprovedli Hatanaka dekompresi.

Využijte dat z tohoto souboru a spočítejte polohu stanice pomocí pseudovzdáleností získaných z P2-kódu na druhé nosné vlně. Použijte data pro epochy v časech $t_{1}$ , $t_{2}$ , $t_{3}$ a využijte všechna měření na družice systému GPS. Kromě polohy stanice určete také opravu hodin přijímace pro každou epochu. Jako apriorní přibližné souřadnice stanice využijte souřadnice zadané v souboru RINEX. Polohu družic spočítejte ve svém programu za použití přesných efemerid z navigačního rinexu a opravu družicových hodin nalezněte tamtéž.

Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.

Podrobný popis formátu RINEX a SP3c naleznete v dokumentech:

Numerické zadání úlohy

číslo zadání	stanice EPN	epocha $t_{1}$			epocha $t_{2}$			epocha $t_{3}$
číslo zadání	stanice EPN	hod	min	sec	hod	min	sec	hod	min	sec
1	brux1420.13d	22	45	0.000000	23	0	0.000000	23	15	0.000000
2	bucu1420.13d	0	15	0.000000	0	30	0.000000	0	45	0.000000
3	bzrg1420.13d	1	0	0.000000	1	15	0.000000	1	30	0.000000
4	cags1420.13d	1	45	0.000000	2	0	0.000000	2	15	0.000000
5	cagz1420.13d	2	30	0.000000	2	45	0.000000	3	0	0.000000
6	ccj21420.13d	3	15	0.000000	3	30	0.000000	3	45	0.000000
7	cebr1420.13d	4	0	0.000000	4	15	0.000000	4	30	0.000000
8	cedu1420.13d	4	45	0.000000	5	0	0.000000	5	15	0.000000
9	chwk1420.13d	5	30	0.000000	5	45	0.000000	6	0	0.000000
10	chan1420.13d	6	15	0.000000	6	30	0.000000	6	45	0.000000
11	cit11420.13d	7	0	0.000000	7	15	0.000000	7	30	0.000000
12	chil1420.13d	7	45	0.000000	8	0	0.000000	8	15	0.000000
13	chpi1420.13d	8	30	0.000000	8	45	0.000000	9	0	0.000000
14	ffmj1420.13d	9	15	0.000000	9	30	0.000000	9	45	0.000000
15	chum1420.13d	10	0	0.000000	10	15	0.000000	10	30	0.000000
16	chur1420.13d	10	45	0.000000	11	0	0.000000	11	15	0.000000
17	chwk1420.13d	11	30	0.000000	11	45	0.000000	12	0	0.000000
18	clar1420.13d	12	15	0.000000	12	30	0.000000	12	45	0.000000
19	cnmr1420.13d	13	0	0.000000	13	15	0.000000	13	30	0.000000
20	coco1420.13d	13	45	0.000000	14	0	0.000000	14	15	0.000000
21	frdn1420.13d	14	30	0.000000	14	45	0.000000	15	0	0.000000
22	conz1420.13d	15	15	0.000000	15	30	0.000000	15	45	0.000000
23	coso1420.13d	16	0	0.000000	16	15	0.000000	16	30	0.000000
24	crao420.13d	1	15	0.000000	1	30	0.000000	1	45	0.000000
25	gode1420.13d	2	15	0.000000	2	30	0.000000	2	45	0.000000
26	cusv1420.13d	3	15	0.000000	3	30	0.000000	3	45	0.000000
27	lpal1420.13d	4	15	0.000000	4	30	0.000000	4	45	0.000000
28	dakr1420.13d	5	15	0.000000	5	30	0.000000	5	45	0.000000
29	darw1420.13d	6	15	0.000000	6	30	0.000000	6	45	0.000000
30	dgar1420.13d	7	15	0.000000	7	30	0.000000	7	45	0.000000
31	drao1420.13d	8	15	0.000000	8	30	0.000000	8	45	0.000000
32	dubo1420.13d	9	15	0.000000	9	30	0.000000	9	45	0.000000
33	geno1420.13d	11	15	0.000000	11	30	0.000000	11	45	0.000000
34	fair1420.13d	12	15	0.000000	12	30	0.000000	12	45	0.000000
35	koug1420.13d	13	15	0.000000	13	30	0.000000	13	45	0.000000

@@ Řádek 2: / Řádek 2: @@
 GPS - zpracování kódových měření
 ==Zadání úlohy==
-Například v [ftp://igs.ensg.ign.fr/pub/igs/data/2009/099 adresáři /pub/igs/data/2009/099/ protokolu FTP na igs.ensg.ign.fr] si stáhněte dle numerického zadání soubor ve formátu RINEX obsahující data měřená přijímačem GNSS stanice sítě EUREF EPN v rozmezí jednoho dne v GPS týdnu 1526. Na soubory je použita jak komprese pomocí SW ''compress'' (''.Z''), tak komprese ''Hatanaka'' (''.RRd'' místo ''.RRo''). WinZip formát ''.Z'' rozbalit neumí, ve Windows použijte program [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/compress.exe compress.exe] ([ftp://athena.fsv.cvut.cz/VG/VYG1/SW/compress.txt návod], [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/decompr.bat decompr.bat]), nebo [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/gzip.exe gzip.exe]. O kompresi ''Hatanaka'' se více dočtete [ftp://terras.gsi.go.jp/software/RNXCMP zde], kde je také odkaz na SW na dekompresi pro různé OS (i pro Windows). Kompresi ''.tar'' umí rozbalit např. ''TotalCommander''.
+Například v [ftp://igs.ensg.ign.fr/pub/igs/data/2013/142/ adresáři pub/igs/data/2013/142/] protokolu FTP na igs.ensg.ign.fr si stáhněte dle numerického zadání soubor ve formátu RINEX obsahující data měřená přijímačem GNSS stanice sítě EUREF EPN v rozmezí jednoho dne v GPS týdnu 1741 (povšimněte si významného výročí spadajícího na tento den). Na soubory je použita jak komprese pomocí SW ''compress'' (''.Z''), tak komprese ''Hatanaka'' (''.RRd'' místo ''.RRo''). WinZip formát ''.Z'' rozbalit neumí, ve Windows použijte program [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/compress.exe compress.exe] ([ftp://athena.fsv.cvut.cz/VG/VYG1/SW/compress.txt návod], [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/decompr.bat decompr.bat]), nebo [ftp://athena.fsv.cvut.cz/VG/VYG1/SW/gzip.exe gzip.exe]. O kompresi ''Hatanaka'' se více dočtete [http://facility.unavco.org/software/preprocessing/preprocessing.html#hatanaka zde], kde je také odkaz na SW na dekompresi pro různé OS (i pro Windows). Kompresi ''.tar'' umí rozbalit např. ''TotalCommander''. Pokud se vám zdá, že váš rinex neodpovídá formátu souboru rinex v níže uvedeném manuálu rinex300.pdf, patrně jste neprovedli Hatanaka dekompresi.
 Využijte dat z tohoto souboru a spočítejte polohu stanice pomocí pseudovzdáleností získaných z P2-kódu na druhé nosné vlně.
 Použijte data pro epochy v časech <math>t_1</math>, <math>t_2</math>, <math>t_3</math> a využijte všechna měření na družice systému GPS. Kromě polohy stanice určete také opravu hodin přijímace pro každou epochu.
 Jako apriorní přibližné souřadnice stanice využijte souřadnice zadané v souboru RINEX.
-Polohu družic a opravu družicových hodin spočítejte na základě souboru formátu SP3 obsahujícího přesné efemeridy družic.
+Polohu družic spočítejte ve svém programu za použití přesných efemerid z navigačního rinexu a opravu družicových hodin nalezněte tamtéž.
+<!--
 Výpočet polohy družice v okamžiku vyslání signálu vypočítejte pomocí kvadratické interpolace ze třech etap souboru SP3, časově nejbližších zadaným epochám <math>t_1</math>, <math>t_2</math>, <math>t_3</math>.
 Soubory SP3 dohledejte na internetových stránkách sítě [http://www.epncb.oma.be EUREF EPN] nebo na stránkách služby [http://igscb.jpl.nasa.gov IGS] (pomůcka - oproti měřeným datům - observačním souborům, se soubory SP3 nacházejí v adresáři ''products'', ne ''data'', a jsou řazeny ne podle roků a dnů v roce, ale podle čísla GPS týdne).
+-->
 Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.
@@ Řádek 18: / Řádek 21: @@
 *[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/sp3c.txt sp3c.txt]
 *[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/gps_zerocode.pdf gps_zerocode.pdf]
+*[http://www.colorado.edu/engineering/ASEN/asen6090/rinex_files.html rinex_files]
 == Numerické zadání úlohy ==
+{| class = "border"
+!rowspan="2" valign="top"| číslo zadání
+!rowspan="2" valign="top"|  stanice EPN
+!colspan="3"| epocha <math>t_1</math>
+!colspan="3"| epocha <math>t_2</math>
+!colspan="3"| epocha <math>t_3</math>
+|-
+!halign="right"|hod!!min!!sec!!hod!!min!!sec!!hod!!min!!sec
+|-
+| 1 ||   brux1420.13d  || 22 || 45 ||    0.000000 || 23 || 0 ||    0.000000 || 23 || 15 ||    0.000000
+|-
+| 2 ||   bucu1420.13d  || 0 || 15 ||    0.000000 || 0 || 30 ||    0.000000 || 0 || 45 ||    0.000000
+|-
+| 3 ||   bzrg1420.13d  || 1 || 0 ||    0.000000 || 1 || 15 ||    0.000000 || 1 || 30 ||    0.000000
+|-
+| 4 ||   cags1420.13d  || 1 || 45 ||    0.000000 || 2 || 0 ||    0.000000 || 2 || 15 ||    0.000000
+|-
+| 5 ||   cagz1420.13d  || 2 || 30 ||    0.000000 || 2 || 45 ||    0.000000 || 3 || 0 ||    0.000000
+|-
+| 6 ||   ccj21420.13d  || 3 || 15 ||    0.000000 || 3 || 30 ||    0.000000 || 3 || 45 ||    0.000000
+|-
+| 7 ||   cebr1420.13d  || 4 || 0 ||    0.000000 || 4 || 15 ||    0.000000 || 4 || 30 ||    0.000000
+|-
+| 8 ||   cedu1420.13d  || 4 || 45 ||    0.000000 || 5 || 0 ||    0.000000 || 5 || 15 ||    0.000000
+|-
+| 9 ||   chwk1420.13d  || 5 || 30 ||    0.000000 || 5 || 45 ||    0.000000 || 6 || 0 ||    0.000000
+|-
+| 10 ||   chan1420.13d  || 6 || 15 ||    0.000000 || 6 || 30 ||    0.000000 || 6 || 45 ||    0.000000
+|-
+| 11 ||  cit11420.13d  || 7 || 0 ||    0.000000 || 7 || 15 ||    0.000000 || 7 || 30 ||    0.000000
+|-
+| 12 ||  chil1420.13d  || 7 || 45 ||    0.000000 || 8 || 0 ||    0.000000 || 8 || 15 ||    0.000000
+|-
+| 13 ||  chpi1420.13d  || 8 || 30 ||    0.000000 || 8 || 45 ||    0.000000 || 9 || 0 ||    0.000000
+|-
+| 14 ||  ffmj1420.13d  || 9 || 15 ||    0.000000 || 9 || 30 ||    0.000000 || 9 || 45 ||    0.000000
+|-
+| 15 ||  chum1420.13d  || 10 || 0 ||    0.000000 || 10 || 15 ||    0.000000 || 10 || 30 ||    0.000000
+|-
+| 16 ||  chur1420.13d  || 10 || 45 ||    0.000000 || 11 || 0 ||    0.000000 || 11 || 15 ||    0.000000
+|-
+| 17 ||  chwk1420.13d  || 11 || 30 ||    0.000000 || 11 || 45 ||    0.000000 || 12 || 0 ||    0.000000
+|-
+| 18 ||  clar1420.13d  || 12 || 15 ||    0.000000 || 12 || 30 ||    0.000000 || 12 || 45 ||    0.000000
+|-
+| 19 ||  cnmr1420.13d  || 13 || 0 ||    0.000000 || 13 || 15 ||    0.000000 || 13 || 30 ||    0.000000
+|-
+| 20 ||  coco1420.13d  || 13 || 45 ||    0.000000 || 14 || 0 ||    0.000000 || 14 || 15 ||    0.000000
+|-
+| 21 ||  frdn1420.13d  || 14 || 30 ||    0.000000 || 14 || 45 ||    0.000000 || 15 || 0 ||    0.000000
+|-
+| 22 ||  conz1420.13d  || 15 || 15 ||    0.000000 || 15 || 30 ||    0.000000 || 15 || 45 ||    0.000000
+|-
+| 23 ||  coso1420.13d  || 16 || 0 ||    0.000000 || 16 || 15 ||    0.000000 || 16 || 30 ||    0.000000
+|-
+| 24 ||  crao420.13d  || 1 || 15 ||    0.000000 || 1 || 30 ||    0.000000 || 1 || 45 ||    0.000000
+|-
+| 25 ||  gode1420.13d  || 2 || 15 ||    0.000000 || 2 || 30 ||    0.000000 || 2 || 45 ||    0.000000
+|-
+| 26 ||  cusv1420.13d  || 3 || 15 ||    0.000000 || 3 || 30 ||    0.000000 || 3 || 45 ||    0.000000
+|-
+| 27 ||  lpal1420.13d  || 4 || 15 ||    0.000000 || 4 || 30 ||    0.000000 || 4 || 45 ||    0.000000
+|-
+| 28 ||  dakr1420.13d  || 5 || 15 ||    0.000000 || 5 || 30 ||    0.000000 || 5 || 45 ||    0.000000
+|-
+| 29 ||  darw1420.13d  || 6 || 15 ||    0.000000 || 6 || 30 ||    0.000000 || 6 || 45 ||    0.000000
+|-
+| 30 ||  dgar1420.13d  || 7 || 15 ||    0.000000 || 7 || 30 ||    0.000000 || 7 || 45 ||    0.000000
+|-
+| 31 ||  drao1420.13d  || 8 || 15 ||    0.000000 || 8 || 30 ||    0.000000 || 8 || 45 ||    0.000000
+|-
+| 32 ||  dubo1420.13d  || 9 || 15 ||    0.000000 || 9 || 30 ||    0.000000 || 9 || 45 ||    0.000000
+|-
+| 33 ||  geno1420.13d  || 11 || 15 ||    0.000000 || 11 || 30 ||    0.000000 || 11 || 45 ||    0.000000
+|-
+| 34 ||  fair1420.13d  || 12 || 15 ||    0.000000 || 12 || 30 ||    0.000000 || 12 || 45 ||    0.000000
+|-
+| 35 ||  koug1420.13d  || 13 || 15 ||    0.000000 || 13 || 30 ||    0.000000 || 13 || 45 ||    0.000000
+|}
+{{Teoretická geodézie}}