155TG2 Teoretická geodézie 2 / úloha 3: Porovnání verzí

Z GeoWikiCZ
VizuálníWikitext
mBez shrnutí editace
Řádek 1: Řádek 1:
==Název úlohy==
==Název úlohy==
GPS - zpracování kódových měření
GPS - zpracování kódových měření
==Zadání úlohy==
Využijte dat ze zadaného souboru GNSS observací (*.yyo) a spočítejte polohu stanice pomocí pseudovzdáleností získaných ze zadaného typu GNSS observací (zadán typ kódového měření). 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í zadaného typu na všechny družice systému GPS v uvedených třech epochách. 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 uvedené v observačním RINEXu. Polohu družic v okamžiku vyslání signálu spočítejte ve svém programu z úlohy 1 za použití vysílaných efemerid z navigačního RINEXu; opravu družicových hodin naleznete tamtéž.
==Numerické zadání==
Numerické zadání (obsahující název observačního RINEXu, zadaný typ kódového měření a trojici výpočetních epoch) naleznete v adresáři http://athena.fsv.cvut.cz:8000/TG2/code/zadani v souboru '''tg2_2022_u3_xx.m''', kde '''xx''' je číslo zadání. Číslo zadání studenta odpovídá číslování uvedenému na stránkách cvičení TGD2. (zatím není připraveno)
Observační RINEXy i jim příslušné navigační RINEXy ze shodného dne (se shodným DOY), které jsou pro vás již staženy a dekomprimovány z permanentní GNSS sítě CZEPOS, naleznete v adresáři http://athena.fsv.cvut.cz:8000/TG2/code/data.
Pro výpočet uvažujte geocentrickou gravitační konstantu <math>GM = 398600.5 km^3s^{-2}</math>.
Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.
==Dokumenty ke stažení==
Poznámky ke zpracování kódových měření naleznete v souboru [http://athena.fsv.cvut.cz:8000/TG2/code/gps_zerocode.pdf gps_zerocode.pdf]. (dokument zatím není k dispozici, bude doplněno)
Podrobný popis formátu RINEX verze 3 (jak observačních, tak navigačních souborů) naleznete v [http://athena.fsv.cvut.cz:8000/TG2/code/rinex304.pdf zde]
Výpočet poloh družic v daném čase z navigační zprávy již znáte z 1.úlohy a je popsán v souboru [http://athena.fsv.cvut.cz:8000/TG2/code/poznamky_druzice.pdf poznamky_druzice.pdf], který již též máte z 1.úlohy.
<!--
==Zadání úlohy==
==Zadání úlohy==
Například v adresáři [ftp://igs.ensg.ign.fr/pub/igs/data/2016/103/ pub/igs/data/2016/103/] 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ýdne 1892<!--(povšimněte si významného výročí spadajícího na tento den)-->. Případně při nefunkčnosti FTP na webu: https://webigs.ign.fr/gdc/en/data/search vyberte svou stanici a datum (103. den roku 2016: 12.4.2016) a data si stáhněte z košíku (Basket).  
Například v adresáři [ftp://igs.ensg.ign.fr/pub/igs/data/2016/103/ pub/igs/data/2016/103/] 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ýdne 1892. Případně při nefunkčnosti FTP na webu: https://webigs.ign.fr/gdc/en/data/search vyberte svou stanici a datum (103. den roku 2016: 12.4.2016) a data si stáhněte z košíku (Basket).  
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 na stránkách [http://unavco.org/software/preprocessing/preprocessing.html#hatanaka unavco.org], 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.  
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 na stránkách [http://unavco.org/software/preprocessing/preprocessing.html#hatanaka unavco.org], 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.  


Řádek 9: Řádek 32:
Jako apriorní přibližné souřadnice stanice využijte souřadnice zadané v souboru RINEX.  
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éž.
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>.
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).
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.
Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.
<!--
'''Poznámka, zejména pro kruh 17 (1. cvičení): Dejte si pozor na správný počet sloupců v matici A, díky 3 epochám měření je třeba vypočítat 3 chyby přijímače, takže matice A bude mít pro tento případ 6 sloupců - viz vzorec 12 na straně 5 gps-zerocode.pdf ''' -->


Podklady k úloze:  
Podklady k úloze:  
* popis výpočtu [ftp://athena.fsv.cvut.cz/VG/VYG1/texty/gps_zerocode.pdf gps_zerocode.pdf]
* popis výpočtu [ftp://athena.fsv.cvut.cz/VG/VYG1/texty/gps_zerocode.pdf gps_zerocode.pdf]
* popis formátu rinex [ftp://athena.fsv.cvut.cz/VG/VYG1/texty/rinex300.pdf rinex300.pdf]
* popis formátu rinex [ftp://athena.fsv.cvut.cz/VG/VYG1/texty/rinex300.pdf rinex300.pdf]
<!-- *[http://www.colorado.edu/engineering/ASEN/asen6090/rinex_files.html rinex_files] -->


== Numerické zadání úlohy ==
== Numerické zadání úlohy ==
Řádek 41: Řádek 62:
| Boušek Martin ||  cags1030.16d  || 10 || 45 ||    0.000000 || 11 || 0 ||    0.000000 || 15 || 15 ||    0.000000  
| Boušek Martin ||  cags1030.16d  || 10 || 45 ||    0.000000 || 11 || 0 ||    0.000000 || 15 || 15 ||    0.000000  
|-
|-
|Horník Martin ||  cagz1030.16d  || 11 || 45 ||    0.000000 || 13 || 30 ||    0.000000 || 17 || 0 ||    0.000000  
| Horník Martin ||  cagz1030.16d  || 11 || 45 ||    0.000000 || 13 || 30 ||    0.000000 || 17 || 0 ||    0.000000  
|-
|-
| Korec Pavel ||  ccj21030.16d  || 13 || 15 ||    0.000000 || 14 || 30 ||    0.000000 || 16 || 45 ||    0.000000  
| Korec Pavel ||  ccj21030.16d  || 13 || 15 ||    0.000000 || 14 || 30 ||    0.000000 || 16 || 45 ||    0.000000  
Řádek 67: Řádek 88:
| Trojan Jakub ||  chwk1030.16d  || 11 || 45 ||    0.000000 || 13 || 45 ||    0.000000 || 17 || 45 ||    0.000000  
| Trojan Jakub ||  chwk1030.16d  || 11 || 45 ||    0.000000 || 13 || 45 ||    0.000000 || 17 || 45 ||    0.000000  
|-
|-
|Váchová Hana ||  clar1030.16d  || 8 || 45 ||    0.000000 || 13 || 45 ||    0.000000 || 18 || 45 ||    0.000000  
| Váchová Hana ||  clar1030.16d  || 8 || 45 ||    0.000000 || 13 || 45 ||    0.000000 || 18 || 45 ||    0.000000  
|-
|-
|Váňa Ondřej ||  cnmr1030.16d  || 3 || 45 ||    0.000000 || 12 || 15 ||    0.000000 || 19 || 30 ||    0.000000  
| Váňa Ondřej ||  cnmr1030.16d  || 3 || 45 ||    0.000000 || 12 || 15 ||    0.000000 || 19 || 30 ||    0.000000  
<!--
|-
|-
| ||  coco1030.16d  || 11 || 45 ||    0.000000 || 14 || 15 ||    0.000000 || 19 || 45 ||    0.000000  
| ||  coco1030.16d  || 11 || 45 ||    0.000000 || 14 || 15 ||    0.000000 || 19 || 45 ||    0.000000  
Řádek 117: Řádek 137:
|-
|-
|  ||  madr1030.16d  || 13 || 15 ||    0.000000 || 13 || 30 ||    0.000000 || 13 || 45 ||    0.000000  
|  ||  madr1030.16d  || 13 || 15 ||    0.000000 || 13 || 30 ||    0.000000 || 13 || 45 ||    0.000000  
|}
-->
-->
|}




{{Teoretická geodézie}}
{{Teoretická geodézie}}

Verze z 5. 4. 2022, 11:29

Název úlohy

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

Zadání úlohy

Využijte dat ze zadaného souboru GNSS observací (*.yyo) a spočítejte polohu stanice pomocí pseudovzdáleností získaných ze zadaného typu GNSS observací (zadán typ kódového měření). Použijte data pro epochy v časech , , a využijte všechna měření zadaného typu na všechny družice systému GPS v uvedených třech epochách. 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 uvedené v observačním RINEXu. Polohu družic v okamžiku vyslání signálu spočítejte ve svém programu z úlohy 1 za použití vysílaných efemerid z navigačního RINEXu; opravu družicových hodin naleznete tamtéž.

Numerické zadání

Numerické zadání (obsahující název observačního RINEXu, zadaný typ kódového měření a trojici výpočetních epoch) naleznete v adresáři http://athena.fsv.cvut.cz:8000/TG2/code/zadani v souboru tg2_2022_u3_xx.m, kde xx je číslo zadání. Číslo zadání studenta odpovídá číslování uvedenému na stránkách cvičení TGD2. (zatím není připraveno)

Observační RINEXy i jim příslušné navigační RINEXy ze shodného dne (se shodným DOY), které jsou pro vás již staženy a dekomprimovány z permanentní GNSS sítě CZEPOS, naleznete v adresáři http://athena.fsv.cvut.cz:8000/TG2/code/data.

Pro výpočet uvažujte geocentrickou gravitační konstantu . Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.

Dokumenty ke stažení

Poznámky ke zpracování kódových měření naleznete v souboru gps_zerocode.pdf. (dokument zatím není k dispozici, bude doplněno)

Podrobný popis formátu RINEX verze 3 (jak observačních, tak navigačních souborů) naleznete v zde

Výpočet poloh družic v daném čase z navigační zprávy již znáte z 1.úlohy a je popsán v souboru poznamky_druzice.pdf, který již též máte z 1.úlohy.

Citováno z „https://geo.fsv.cvut.cz/index.php?title=155TG2_Teoretická_geodézie_2_/_úloha_3&oldid=39447“