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

Z GeoWikiCZ
 
(Není zobrazeno 8 mezilehlých verzí od 2 dalších uživatelů.)
Řádek 1: Řádek 1:
==Název úlohy==
==Název úlohy==
GPS - výpočet polohy stanice pomocí dvojitých diferencí
GPS - výpočet polohy stanice pomocí dvojitých diferencí
==Zadání úlohy==
==Zadání úlohy==
K dispozici máte observační a navigační soubory ve formátu RINEX pořízené na referenční stanici a na určované stanici současně. Vypočítejte polohu připojovaného (určovaného) bodu relativně k referenční stanici (pro polohu referenční stanice vstupující do vyrovnání použijte přibližné souřadnice přijímače referenční stanice uvedené v observačním souboru formátu RINEX). Pracujte s pseudovzdálenostmi získanými na zadaném typu kódového měření v časových epochách <math>t_1</math>, <math>t_2</math>, <math>t_3</math> observačních souborů RINEX z referenční i určované stanice a výpočet proveďte zpracováním dvojitě diferencovaných měření. Při výpočtu berte v úvahu všechny GPS družice společné pro referenční a připojovaný bod v dané epoše.
==Numerické zadání==
Numerický výpočet v LS2024 není vzhledem k okolnostem povinný. V případě zájmu bude numerické zadání včetně dat potřebných pro výpočet poskytnuto na vyžádání (v takovém případě se výpočet stane povinně závazným).
<!--
Numerické zadání (obsahující název observačního RINEXu referenční stanice, název observačního RINEXu určované stanice, 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/doubledif/zadani v souboru '''tg2_2022_u4_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.


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/doubledif/data.
-->
Pro výpočet uvažujte geocentrickou gravitační konstantu <math>GM = 398600.5 km^3s^{-2}</math> a úhlovou rychlost rotace Země <math>\omega_E = 7292115.1467*10^{-11} s^{-1}</math>. Jako rychlost světla použijte hodnotu <math>c = 299792458.0 m/s</math>.
==Dokumenty ke stažení==
Stručný popis zpracování dvojitě diferencovaných kódových měření je uveden v pdf-dokumentu (gps_zerocode_stare.pdf), který máte k dispozici ze 3.úlohy.
<!--
Poznámky ke zpracování kódových měření naleznete v souboru [http://athena.fsv.cvut.cz:8000/TG2/code/gps_zerocode_stare.pdf gps_zerocode.pdf]. (jde zatím o původní neaktualizovaný soubor s chybami, které však pozorný čtenář může odhalit; aktualizovaný a opravený soubor bude doplněn)
Podrobný popis formátu RINEX verze 3 (jak observačních, tak navigačních souborů) naleznete [http://athena.fsv.cvut.cz:8000/TG2/code/rinex304.pdf zde] (např. na stranách 16-17 je uveden podrobný seznam identifikátorů jednotlivých typů GPS observací, na dalších stranách následují identifikátory typů observací ostatních GNSS systémů).
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.
V případě zájmu si též můžete stáhnout a použít funkci [http://athena.fsv.cvut.cz:8000/TG2/code/readrinexo.m readrinexo.m] pro načtení celého observačního RINEXu (verze 3.04) pro účely 3.úlohy TG2 (tj. budou načteny pouze observace všech GPS družic ve všech epochách v RINEXu obsažených, zatímco observace ostatních GNSS systémů jsou ignorovány). Její použití však předpokládá znalost a využití nejen numerických a řetězcových matic, ale i datového typu buňka (cell). Vzhledem k tomu, že ve vašem zadání máte zpracovávat data pouze ze tří zadaných epoch, je asi jednodušší vyhledat si data trojice epoch ručně.
Máte k dispozici observační a navigační soubory formátu RINEX, pořízené pro body 4003 a 4004.
-->
<!--
==Zadání úlohy==
Máte k dispozici observační a navigační soubory formátu RINEX, pořízené pro body 4003 a 4004.
Máte k dispozici observační a navigační soubory formátu RINEX, pořízené pro body 4003 a 4004.


Vypočítejte polohu připojovaného bodu 4001/4003 relativně k bodu 4004 (danou polohu bodu 4004 vstupující do vyrovnání použijte přibližné souřadnice přijímače udané v observačním souboru formátu RINEX). Pracujte s pseudovzdálenostmi získanými v P2-kódu ve časových epochách <math>t_1</math>,<math>t_2</math>,<math>t_3</math> observačního souboru RINEX a výpočet proveďte zpracováním dvojitě diferencovaných měření.  
Vypočítejte polohu připojovaného bodu 4001/4003 relativně k bodu 4004 (danou polohu bodu 4004 vstupující do vyrovnání použijte přibližné souřadnice přijímače udané v observačním souboru formátu RINEX). Pracujte s pseudovzdálenostmi získanými v P2-kódu ve časových epochách <math>t_1</math>,<math>t_2</math>,<math>t_3</math> observačního souboru RINEX a výpočet proveďte zpracováním dvojitě diferencovaných měření.  
Při výpočtu berte v úvahu družice společné pro referenční a připojovaný bod.
Při výpočtu berte v úvahu družice společné pro referenční a připojovaný bod.
<!--
Polohu družic 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 <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). Souřadnice družic lze též z přesných efemerid, uvedených v navigačním rinexu.
-->
Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.
Jako rychlost světla použijte hodnotu c = 299792458.0 m/s.


Řádek 18: Řádek 43:
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/gps_zerocode.pdf gps_zerocode.pdf]
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/gps_zerocode.pdf gps_zerocode.pdf]
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/DD_dodatek.pdf DD_dodatek.pdf]
*[ftp://athena.fsv.cvut.cz/VG/VYG1/texty/DD_dodatek.pdf DD_dodatek.pdf]


==Numerické zadání==
==Numerické zadání==
Řádek 30: Řádek 54:
||  epocha <math>t_3</math> [s]
||  epocha <math>t_3</math> [s]
|-               
|-               
| Bárta  Aleš || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 49885 || 50485 || 51085
| Adamcová Alice || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 49815 || 50385 || 51080
|-               
|-               
|Belecká  Sabina||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 33310 || 33910 || 34510
|Beck Jakub||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 33315 || 33915 || 34110
|-               
|-               
|Běloch  Lukáš ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 48150 || 48750 || 49350
|Bláhová Taťána ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 48145 || 48755 || 49250
|-               
|-               
|Brachtl  Ondřej ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 33365 || 33965 || 34565
|Boušek Martin ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 33165 || 33865 || 34515
|-               
|-               
|Brňovják  Vojtěch ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvutz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42000 || 42600 || 43200
|Horník Martin ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvutz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42010 || 42605 || 43210
|-               
|-               
|Hoffmann  Marek ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 36265 || 36865 || 37465
|Korec Pavel ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 36165 || 36825 || 37265
|-               
|-               
|Hošek  Lukáš ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 45440 || 46040 || 46640
|Krauz Tomáš ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 45420 || 46045 || 46130
|-               
|-               
|Hrubý  Jan ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 37290 || 37890 || 38490
|Kučera Jakub ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 37190 || 37895 || 38195
|-               
|-               
|Kačerovský  Jonáš ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 49265 || 49865 || 50465
|Kučerová Adéla ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 49165 || 49875 || 50435
|-               
|-               
|Kličková  Sabina ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 35300 || 35900 || 36500
|Martoš Samuel ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 35100 || 35910 || 36505
|-               
|-               
|Králová  Tereza ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42040 || 42640 || 43240
|Pelikán Václav ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42140 || 42645 || 43340
|-               
|-               
|Křížek  Vojtěch ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 48590 || 49190 || 49790
|Sochor Miroslav ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 48190 || 49155 || 49780
|-               
|-               
|Křížová  Monika ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 49415 || 50015 || 50615
|Straka Tomáš ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 49215 || 50115 || 50645
|-               
|-               
|Kučera  Petr ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 38770 || 39370 || 39970
|Štropová Barbara ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 38170 || 39375 || 39175
|-               
|-               
|Kučera  Jan ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 40560 || 41160 || 41760
|Tichá Anna ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 40160 || 41135 || 41265
|-               
|-               
|Kuzmanov  Pane ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 41155 || 41755 || 42355
|Tomášek Filip ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 41055 || 41760 || 42255
|-               
|-               
|Masák  Václav ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 35275 || 35875 || 36475
|Trojan Jakub ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 35175 || 35810 || 36420
|-               
|-               
|Mrkvička  Jiří ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 35510 || 36110 || 36710
|Váchová Hana ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 35500 || 36185 || 36210
|-               
|-               
|Mužík  František ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42230 || 42830 || 43430
|Váňa Ondřej ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42030 || 42835 || 43220
|-               
|-               
|Novák  Michal ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 45620 || 46220 || 46820
|Novák  Michal ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 45620 || 46220 || 46820
Řádek 77: Řádek 101:
|-               
|-               
|Zimmerhaklová  Kateřina ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42590 || 43190 || 43790
|Zimmerhaklová  Kateřina ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 42590 || 43190 || 43790
<!--
|-               
|-               
| ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 47590 || 48190 || 48790
| ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 47590 || 48190 || 48790
Řádek 118: Řádek 141:
|-
|-
| M.Zikmundová ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 46450 || 46780 || 48980
| M.Zikmundová ||  [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4004.zip 4004] || [ftp://athena.fsv.cvut.cz/VG/VYG1/zadani/uloha3/4003.zip 4003] || 46450 || 46780 || 48980
|}
-->
-->
|}


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

Aktuální verze z 22. 4. 2024, 17:07

Název úlohy

GPS - výpočet polohy stanice pomocí dvojitých diferencí

Zadání úlohy

K dispozici máte observační a navigační soubory ve formátu RINEX pořízené na referenční stanici a na určované stanici současně. Vypočítejte polohu připojovaného (určovaného) bodu relativně k referenční stanici (pro polohu referenční stanice vstupující do vyrovnání použijte přibližné souřadnice přijímače referenční stanice uvedené v observačním souboru formátu RINEX). Pracujte s pseudovzdálenostmi získanými na zadaném typu kódového měření v časových epochách Nelze pochopit (SVG, alternativně PNG (MathML lze povolit skrze prohlížečový plugin): Neplatná odpověď („Math extension cannot connect to Restbase.“) od serveru „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t_1} , Nelze pochopit (SVG, alternativně PNG (MathML lze povolit skrze prohlížečový plugin): Neplatná odpověď („Math extension cannot connect to Restbase.“) od serveru „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t_2} , Nelze pochopit (SVG, alternativně PNG (MathML lze povolit skrze prohlížečový plugin): Neplatná odpověď („Math extension cannot connect to Restbase.“) od serveru „https://wikimedia.org/api/rest_v1/“:): {\displaystyle t_3} observačních souborů RINEX z referenční i určované stanice a výpočet proveďte zpracováním dvojitě diferencovaných měření. Při výpočtu berte v úvahu všechny GPS družice společné pro referenční a připojovaný bod v dané epoše.

Numerické zadání

Numerický výpočet v LS2024 není vzhledem k okolnostem povinný. V případě zájmu bude numerické zadání včetně dat potřebných pro výpočet poskytnuto na vyžádání (v takovém případě se výpočet stane povinně závazným).

Pro výpočet uvažujte geocentrickou gravitační konstantu Nelze pochopit (SVG, alternativně PNG (MathML lze povolit skrze prohlížečový plugin): Neplatná odpověď („Math extension cannot connect to Restbase.“) od serveru „https://wikimedia.org/api/rest_v1/“:): {\displaystyle GM = 398600.5 km^3s^{-2}} a úhlovou rychlost rotace Země Nelze pochopit (SVG, alternativně PNG (MathML lze povolit skrze prohlížečový plugin): Neplatná odpověď („Math extension cannot connect to Restbase.“) od serveru „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \omega_E = 7292115.1467*10^{-11} s^{-1}} . Jako rychlost světla použijte hodnotu Nelze pochopit (SVG, alternativně PNG (MathML lze povolit skrze prohlížečový plugin): Neplatná odpověď („Math extension cannot connect to Restbase.“) od serveru „https://wikimedia.org/api/rest_v1/“:): {\displaystyle c = 299792458.0 m/s} .

Dokumenty ke stažení

Stručný popis zpracování dvojitě diferencovaných kódových měření je uveden v pdf-dokumentu (gps_zerocode_stare.pdf), který máte k dispozici ze 3.úlohy.