152VYG1 Vyšší geodézie 1 - úloha 1

Termín odevzdání

20.10.

Název úlohy

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

Zadání úlohy

Například v adresáři /pub/igs/data/2011/075/ 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 1627. 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 "okénkový" program 7-zip, 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.

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 a opravu družicových hodin spočítejte na základě souboru formátu SP3 obsahujícího přesné efemeridy družic. 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 $t_{1}$ , $t_{2}$ , $t_{3}$ .

Soubory SP3 dohledejte na internetových stránkách sítě EUREF EPN nebo na stránkách služby 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.

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

Numerické zadání úlohy

číslo zadání	student(ka)	stanice EPN	epocha $t_{1}$			epocha $t_{2}$			epocha $t_{3}$
číslo zadání	student(ka)	stanice EPN	hod	min	sec	hod	min	sec	hod	min	sec
1		abpo0990.09o	0	15	0.000000	0	30	0.000000	0	45	0.000000
2		adis0990.09o	1	00	0.000000	1	15	0.000000	1	30	0.000000
3		aira0990.09o	1	30	0.000000	1	45	0.000000	2	0	0.000000
4		ajac0990.09o	2	15	0.000000	2	30	0.000000	2	45	0.000000
5		albh0990.09o	3	15	0.000000	3	30	0.000000	3	45	0.000000
6		algo0990.09o	4	00	0.000000	4	15	0.000000	4	30	0.000000
7		dlft0990.09o	4	45	0.000000	5	00	0.000000	5	15	0.000000
8		alrt0990.09o	5	30	0.000000	5	45	0.000000	6	00	0.000000
9		dgar0990.09o	6	15	0.000000	6	30	0.000000	6	45	0.000000
10		antc0990.09o	7	0	0.000000	7	15	0.000000	7	30	0.000000
11		areq0990.09o	7	45	0.000000	8	00	0.000000	8	15	0.000000
12		artu0990.09o	8	30	0.000000	8	45	0.000000	9	00	0.000000
13		aspa0990.09o	9	15	0.000000	9	30	0.000000	9	45	0.000000
14		crfp0990.09o	10	00	0.000000	10	15	0.000000	10	30	0.000000
15		baie0990.09o	10	45	0.000000	11	00	0.000000	11	15	0.000000
16		bako0990.09o	11	30	0.000000	11	45	0.000000	12	00	0.000000
17		bamf0990.09o	12	15	0.000000	12	30	0.000000	12	45	0.000000
18		barh0990.09o	13	0	0.000000	13	15	0.000000	13	30	0.000000
19		bdos0990.09o	13	45	0.000000	14	00	0.000000	14	15	0.000000
20		bill0990.09o	14	30	0.000000	14	45	0.000000	15	00	0.000000
21		bjfs0990.09o	15	15	0.000000	15	30	0.000000	15	45	0.000000
22		blyt0990.09o	13	00	0.000000	13	15	0.000000	13	30	0.000000
23		bogi0990.09o	16	45	0.000000	17	00	0.000000	17	15	0.000000
24		bogt0990.09o	17	30	0.000000	17	45	0.000000	18	00	0.000000
25		bran0990.09o	18	15	0.000000	18	30	0.000000	18	45	0.000000
26		crao0990.09o	19	00	0.000000	19	15	0.000000	19	30	0.000000
27		brew0990.09o	19	45	0.000000	20	00	0.000000	20	15	0.000000
28		brft0990.09o	20	30	0.000000	20	45	0.000000	21	00	0.000000
29		coyq0990.09o	21	15	0.000000	21	30	0.000000	21	45	0.000000
30		brst0990.09o	22	00	0.000000	22	15	0.000000	22	30	0.000000
31		brus0990.09o	22	45	0.000000	23	00	0.000000	23	15	0.000000
32		bucu0990.09o	0	15	0.000000	0	30	0.000000	0	45	0.000000

Vyšší geodézie 1