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

Z GeoWikiCZ
 
(Není zobrazeno 7 mezilehlých verzí od stejného uživatele.)
Řá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 67: Řádek 91:
|-               
|-               
|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
|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 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 , , 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 a úhlovou rychlost rotace Země . Jako rychlost světla použijte hodnotu .

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.