Předchozí příspěvek Následující příspěvek
Aspekty přesných GPS aplikací pro speciální geodetické práce
Bureš J., Švábenský O., Weigel J.
Technologie GPS je charakteristická svou vysokou relativní přesností. V případech požadavků co nejvyšší přesnosti při speciálních aplikacích (např. při měření posunů a deformací) je nutné hledat možnosti modelování některých systematických vlivů a použití upravených měřících postupů. Mezi systematické chyby se řadí chyby kosmického systému (např. excentricity fázových center antén družic), dále vlivy prostředí, kterým se signál šíří (tj. ionosféry, troposféry), a systematické chyby přijímače (jako excentricity fázových center antén a jejich variace, vliv vícecestného šíření signálu). Nezanedbatelným je také vliv geometrické konfigurace družic. Při vyhodnocení GPS měření je obtížné jednotlivé vlivy separovat. Některé vlivy je možné ošetřit zaváděním dodatečných korekcí (soubory excentricit fázových center antén družic, globální modely ionosféry a troposféry, soubory excentricit fázových center přijímacích antén). Při kombinaci různých typů aparatur GPS se projeví zejména relativní diference poloh fázových center přijímačů, které jsou při přesných aplikacích nezanedbatelné.
Další možností zvýšení přesnosti je využití korelace dat blízkých GPS aparatur.
Vliv konstelace družic
Vliv konstelace družic je možné výrazně eliminovat opakováním kratších observací ve vhodném časovém odstupu. V tabulce č.1 jsou uvedeny některé charakteristiky pro výškovou složku hodinových řešení vektoru 0027-0021 sítě Sněžník. Vyrovnanou triádou se rozumí řešení tří jednohodinových vektorů v časovém odstupu dvou, čtyř a osmi hodin. Vnitřní přesnost je určena jako kvadratický průměr ze směrodatných odchylek vyrovnaných triád. Vnější přesnost je vyjádřena směrodatnou odchylkou výsledků triád. Je zřejmé, že rozdíl minimální a maximální hodnoty řešení hodinových řešení dosahuje až čtyř centimetrů v průběhu 24 hod. a vnitřní přesnost se výrazně liší od vnější přesnosti. Při kombinaci jednohodinových řešení vektorů po osmi hodinách se zde rozdíl maximálního a minimálního řešení snížil prakticky čtyřikrát, vnitřní přesnost se výrazně snížila, ale vnější přesnost se třikrát zvýšila a lépe koresponduje s vnitřní přesností.
Tabulka 1 : Vektor 0027 – 0021 (délka 7094 m , převýšení 430 m)
Výšková složka vektoru |
Hodinová řešení |
Vyrovnané triády z hodinových řešení |
||
po 2 hod |
po 4 hod |
po 8 hod |
||
Střední hodnota [m] |
7,4749 |
7,4724 |
7,4708 |
7,4746 |
Max – Min [m] |
0,0423 |
0.0262 |
0,0161 |
0,0108 |
Vnitřní přesnost [mm] |
0,8 |
4,2 |
5,2 |
5,4 |
Vnější přesnost [mm] |
10,1 |
7,4 |
4,2 |
3,6 |
Poměr přesností |
12,6 |
1,8 |
0,8 |
0,66 |
Graficky je průběh hodinových řešení spolu s řešením triád po osmi hodinách znázorněn v grafu na obr.1. Lomenou čarou je vyznačen průběh hodinových řešení. Průběh triád je znázorněn tečkami. Je zřejmé, že průběh triád má menší rozptyl (do jednoho centimetru) oproti samostatným hodinovým řešením.
Obr. 1
Vliv volby modelů prostředí
Vliv použitých modelů ionosféry a troposféry je rovněž nezanedbatelným. Diference výsledků při použití různých modelů pro výše uvedený vektor dosahovaly hodnot až několika centimetrů. Pro dlouhé observační doby jsou diference velmi malé, řádově v milimetrech. Bez použití modelu ionosféry mohou diference vůči řešení se standardním modelem dosahovat až centimetrových hodnot. Pro kratší observace (délky např. jedné hodiny) v nočních podmínkách jsou opět diference souřadnic při použití různých modelů malé, při zanedbání vlivu ionosféry vzrostou řádově do 5 mm. Naopak u kratších observací za denních podmínek diference souřadnic mohou dosahovat až dvoucentimetrových hodnot a nezavedení ionosférického modelu může znamenat až několikacentimetrové diference v souřadnicích. Podobná situace je i v důsledku volby různých modelů troposféry, kde diference zejména ve výškové složce dosahovaly pro daný vektor až pěticentimetrových hodnot.
Kombinace různých typů GPS aparatur
Při nasazení více GPS přijímačů je vhodné provést jejich srovnávací měření, jehož cílem je ověření vzájemné geometrie antén použitých aparatur. V našem případě jsme se zaměřili na testy excentricit fázových center antén při jejich různém nasměrování. Testy byly realizovány v rámci opakovaných kampaní v síti Sněžník. Ověřována byla možnost rychlé kontroly antén s využitím metody Stop&Go. Pro úplnost je třeba definovat pojem fázového centra jako bodu, ke kterému je měřena geometrická vzdálenost mezi družicí a přijímačem. Variace fázového centra se liší pro jednotlivé frekvence L1, L2 a závisí též na směru příjmu signálu GPS - je tedy funkcí jeho azimutu a zenitového úhlu.
Testovací měření byla realizována s GPS přijímači/anténami stejného typu (Leica) i různých typů (Leica, Ashtech, Trimble). Při měření bylo použito režimu metody Stop&Go, která byla zvolena záměrně neboť umožňuje v relativně krátkém čase určit polohu více bodů. Antény byly proměřovány v poloze otočení na sever, jih , východ a západ. Po společné inicializaci antén se poloha dalších bodů určuje pouze z několika málo epoch. Pro určení spolehlivé polohy stačí zaměření 2 až 3 epoch, přičemž zvětšování počtu epoch nepřináší výraznější zvýšení přesnosti. Relativní rozptyl těchto několika epoch se pohybuje okolo 1 mm a právě této skutečnosti bylo využito při proměřování antén.
Testována byla vždy nejméně trojice přijímačů/antén, které byly rozmístěny přibližně ve tvaru rovnostranného trojúhelníku, vždy s jednou stranou orientovanou na sever. Při testu více antén je vhodné, aby byly rozloženy rovnoměrně okolo jedné centrální antény. Vzdálenost antén byla minimální (cca do 1 m) z důvodu maximální eliminace společných systematických vlivů. Použití tří a více antén umožňuje určit jednoznačně jejich vzájemnou prostorovou polohu přímým klasickým měřením a tím pořídit nezávislý dostatečně přesný srovnávací etalon. Na začátku měření byla provedena společná inicializace všech antén při jejich stejné orientaci (na sever) pro vyřešení ambiguit. Po inicializaci byly proměřované antény otočeny na jih, přičemž se provedlo v této poloze zaznamenání několika určujících epoch a poté byly antény otočeny opět na sever. Při proměřování antén ve čtyřech polohách byly antény otáčeny ve směru světových stran podle schématu SJVZ. Celá procedura se opakovala tak dlouho až každá poloha antény byla určena 8 až 10 krát. Centrální (referenční) anténa zůstala po celou dobu testu orientována na sever. Testovací měření trvalo cca 1 hod. až 2 hod. Parametry observace anténních testů jsou uvedeny v tabulce č.2.
Tabulka č.2: Parametry observace anténních testů
Název |
Test. typy antén |
Délka inicial. |
Poč. urč. epoch |
Celk. délka testu |
Int. zázn. |
Elev. |
Polohy |
ATEST96 | 5xL | 5 min. | 60 |
55 min. | 5s |
15o |
S, J |
ATEST98 | 3xL, 2xT, 2xA | 72 min. | 12 |
1h 54 min. | 5s |
10o |
S, J |
ATEST99 | 2xL, 1xA | 30 min | 12 |
1h 43 min. | 5 s |
5o |
S,J,V,Z |
Naměřená data byla zpracována firemním software SKI 2.3 při použití jednotné elevační masky 15o, standardního modelu ionosféry, modelu troposféry Hopfieldové, řešení na frekvencích L1&L2. Podmínkou úspěšnosti výpočtu je vyřešení ambiguit při inicializaci antén. Z rozptylu řady měření příslušné polohy antény byla vypočtena výběrová směrodatná odchylka aritmetického průměru s, která činila v průměru do 1 mm pro horizontální složku a do 2 mm pro vertikální složku. Referenční souřadnicový etalon byl získán zpracováním klasického měření délek pásmem a převýšení nivelací. Přesnost etalonu byla lepší než 0,5 mm. Základní parametry zpracování dat jsou sestaveny v tabulce č.3.
Tabulka č.3: Parametry zpracování dat anténních testů
Název |
Elev. |
Model ionosf. |
Model troposf. |
řešení |
Počet určení |
sf |
sl |
s h |
ATEST96 | 15o |
Standard |
Hopfield |
L1&L2 |
11xS, 8xJ |
0,9 mm |
0,5 mm |
1,4 mm |
ATEST98 | 15o |
Standard |
Hopfield |
L1&L2 |
6xS,J |
0,8 mm |
0,9 mm |
1,7 mm |
ATEST99 | 15o |
Standard |
Hopfield |
L1&L2 |
9xS,J,V,Z |
0,8 mm |
0,5 mm |
1,1 mm |
Souřadnice WGS-84 byly transformovány s využitím meridiánového (M) a příčného (N) poloměru křivosti do roviny s bodem dotyku v referenční aparatuře. Rovinné souřadnice určené z klasických měření byly transformovány Weitbrechtovou transformací shodnostního typu do výše uvedeného systému roviných souřadnic. Výškové rozdíly byly vypočteny vzhledem k referenční anténě. Z rozdílu rovinných souřadnic GPS a z klasických měření byly vypočteny polohové a výškové diference. Hodnoty horizontálních diferencí vzhledem ke klasicky určené referenční poloze byly zakresleny do grafů na obr.2, 4 a 6, výškové diference vzhledem k nivelovanému převýšení jsou uvedeny v grafech na obr.3, 5 a 7. Číselně byly diference sestaveny do tabulek č.4 až 7. Rozdíly souřadnic při otočení antén na jih a na sever dosahují hodnot až 5 mm u antén stejného typu, což je z hlediska přesnosti určení této diference hodnota statisticky významná.
Vyhodnocení testu antén stejného typu - ATEST96
obr.2 obr.3
Tabulka č.4: Diference antén stejného typu vzhledem k LEICA SR399
Anténa |
Diference (N-Ref.) |
Diference (S-Ref.) |
Diference (S-N) |
||||||
d fe [m] |
d le [m] |
dhe [m] |
d fe [m] |
d le [m] |
dhe [m] |
d fe [m] |
d le [m] |
dhe [m] |
|
Leica SR399 (S3/E) | 0,0020 |
-0,0012 |
0,0000 |
-0,0020 |
0,0012 |
0,0000 |
0,0000 |
0,0000 |
0,0000 |
Leica SR299 (S2/A) | 0,0037 |
-0,0009 |
-0,0023 |
-0,0015 |
0,0014 |
-0,0043 |
-0,0051 |
0,0023 |
-0,0020 |
Leica SR299 (S1/B) | 0,0020 |
-0,0013 |
-0,0012 |
-0,0017 |
0,0016 |
-0,0063 |
-0,0037 |
0,0029 |
-0,0051 |
Leica SR299 (K1/C) | 0,0011 |
-0,0010 |
-0,0023 |
-0,0021 |
0,0014 |
-0,0017 |
-0,0031 |
0,0024 |
0,0005 |
Leica SR299 (K2/D) | 0,0023 |
-0,0014 |
-0,0043 |
-0,0013 |
0,0006 |
-0,0056 |
-0,0036 |
0,0021 |
-0,0014 |
Průměr | 0,0022 |
-0,0012 |
-0,0014 |
-0,0017 |
0,0012 |
-0,0033 |
-0,0039 |
0,0024 |
-0,0020 |
Stand. odch. průměru s | 0,0013 |
0,0006 |
0,0020 |
Vyhodnocení testu antén různého typu - ATEST98
Předmětem testu bylo 7 antén tří typů Leica SR299, SR399, Ashtech Z-XII (Marine, Geodetic) a Trimble 4000 SSi (Compact L1/L2 GND, Compact L1/L2).
obr.4 obr.5
Tabulka č.5: Diference antén různého typu vzhledem k LEICA SR399
Anténa |
Diference (N-Ref.) |
Diference (S-Ref.) |
Diference (S-N) |
||||||||
d fe [m] |
d le [m] |
dhe [m] |
d fe [m] |
d le [m] |
dhe [m] |
d fe [m] |
d le [m] |
dhe [m] |
|||
Leica SR399 Internal (S3/7) | 0,0020 |
-0,0012 |
0,0000 |
-0,0020 |
0,0012 |
0,0000 |
|||||
Leica SR299 Internal (S1/1) | 0,0002 |
-0,0003 |
0,0015 |
-0,0015 |
0,0031 |
-0,0023 |
-0,0018 |
0,0034 |
-0,0038 |
||
Trimble 4000SSi COMPACT L1/L2 GND (T1/2) | -0,0002 |
-0,0004 |
-0,0094 |
-0,0047 |
-0,0009 |
-0,0090 |
-0,0044 |
-0,0005 |
0,0004 |
||
Trimble 4000SSi COMPACT L1/L2 (T2/3) | 0,0005 |
-0,0016 |
0,0003 |
-0,0008 |
-0,0030 |
-0,0016 |
-0,0014 |
-0,0014 |
-0,0019 |
||
Leica SR299 Internal (S2/4) | 0,0022 |
-0,0029 |
0,0002 |
-0,0009 |
0,0014 |
-0,0011 |
-0,0031 |
0,0042 |
-0,0013 |
||
Ashtech Z-XII3 GEODETIC (A1/5) | -0,0040 |
-0,0034 |
-0,0157 |
-0,0049 |
-0,0005 |
-0,0175 |
-0,0009 |
0,0029 |
-0,0018 |
||
Ashtech Z-XII3 MARINE (A2/6) | -0,0025 |
-0,0016 |
-0,0007 |
-0,0011 |
-0,0008 |
-0,0020 |
0,0014 |
0,0008 |
-0,0013 |
||
Průměr | -0,0017 |
0,0016 |
-0,0016 |
||||||||
Stand. odch. průměru s | 0,0011 |
0,0012 |
0,0024 |
Vyhodnocení testu antén různého typu - ATEST99
Předmětem testu byly tři antény Leica SR299, SR399, Ashtech UZ-12 (Geodetic UE00).
obr.6 obr.7
Tabulka č.6: Diferen
ce poloh anténAnténa |
Diference (S-N) |
Diference (W-E) |
||||
d fe [m] |
d le [m] |
dhe [m] |
d fe [m] |
d le [m] |
dhe [m] |
|
Leica SR299 Internal (LS1) | ||||||
Ashtech UZ-12 GEODETIC UE00 | -0,0007 |
0,0022 |
-0,0003 |
-0,0007 |
-0,0003 |
0,0008 |
Leica SR399 Internal (LS3) | -0,0016 |
0,0017 |
-0,0003 |
-0,0021 |
-0,0015 |
0,0019 |
Průměr | -0,0011 |
0,0020 |
-0,0003 |
-0,0014 |
-0,0009 |
0,0014 |
Stand. odch. průměru s | 0,0010 |
0,0008 |
0,0016 |
0,0011 |
0,0005 |
0,0016 |
Tabulka č.7: Diference antén různého typu vzhledem k LEICA SR299
Diference (N-Ref.) |
Diference (S-Ref.) |
Diference (E-Ref.) |
Diference (W-Ref.) |
|||||||||
Anténa |
dfe [m] | dle [m] | dhe [m] | dfe [m] | dle [m] | dhe [m] | dfe [m] | dle [m] | dhe [m] | dfe [m] | dle [m] | dhe [m] |
Leica SR299 Internal (LS1) | 0,0014 |
-0,0016 |
0,0000 |
-0,0014 |
0,0016 |
0,0000 |
0,0000 |
0,0000 |
0,0000 |
0,0000 |
0,0000 |
0,0000 |
Ashtech UZ-12 GEODETIC UE00 | -0,0037 |
-0,0015 |
-0,0016 |
-0,0044 |
0,0007 |
-0,0019 |
-0,0035 |
-0,0009 |
-0,0031 |
-0,0042 |
-0,0012 |
-0,0023 |
Leica SR399 Internal (LS3) | -0,0003 |
-0,0017 |
0,0038 |
-0,0019 |
0,0000 |
0,0035 |
0,0005 |
-0,0001 |
0,0024 |
-0,0016 |
-0,0015 |
0,0044 |
Z výsledků vyhodnocení vyplynulo, že hodnoty diferencí při různé orientaci antén mohou dosahovat hodnot až 5 mm u antén stejného typu. Vzhledem ke krátké délce testu, kdy se konfigurace družic v průběhu testu nijak výrazně nemění, lze diference interpretovat jako vliv změny polohy fázového centra v důsledku změny azimutu natočení antény. Při proměření antén se diference jednotlivých poloh jeví jako systematické avšak různé pro každý typ antény. Vzájemný rozptyl poloh antén typů Ashtech a Trimble se jeví menší oproti typu Leica. Při porovnání absolutních rozdílů poloh těžišť při různých orientacích antén oproti referenční poloze (získané z terestrických měření) je poloha těžiště u aparatur Leica příznivější. Tento postřeh dokumentuje téměř symetrické rozložení poloh antén vzhledem k referenční poloze v testu antén stejného typu (ATEST96). Tuto skutečnost jistě nelze na základě pouze těchto měření zevšeobecňovat.
Vyhodnocení výškových diferencí antén různého typu (ATEST98) přineslo zajímavé zjištění, že antény s GP (Ground Plane) mají poněkud odlišné vlastnosti, což se projevilo zejména ve výškové složce. Anténa typu Trimble Compact L1/L2 GND a anténa Ashtech Geodetic vykázaly relativní výškovou diferenci 9 mm a 15 mm vůči ostatním typům bez GP.
Aplikace výsledků v lokální síti Sněžník
Vektory zaměřené různými typy aparatur v síti Sněžník vykazovaly při zpracování větší odchylky výškové složky než odchylky při použití aparatur stejného typu. Na vektorech zaměřených kombinací testovaných aparatur byly zavedeny korekce určené z testovacího měření.
Vypočítané hodnoty převýšení jednotlivých vektorů byly porovnány s hodnotami určenými ze spolehlivých referenčních souřadnic. Vyhodnocení bylo provedeno softwarem SKI 2.2, korekce byly zaváděny formou oprav hodnot výšek antén v datových souborech. Výsledky bez korekcí a se zavedením korekcí jsou uvedeny v tabulce č.8. Zavádění korekcí se jeví jako přínosné a ve většině případů došlo ke zlepšení výsledků.
Tabulka č.8: Výsledky zavedení diferencí z testovacích měření
VÝŠKOVÉ ODCHYLKY Z KAMPANĚ SNEZ98 |
|||
bez korekce |
s korekcí |
||
Bod |
Aparatura |
D h ref. - Dh SNEZ98 |
D h ref. - Dh SNEZ98 |
VYHL0027 | Leica SR 399 | ref. |
ref. |
VESE0031 | Ashtech Z-XII-Geodetic | 0,015 m |
0,001 m |
ROUD0030 | Trimble 4000SSi GND | 0,034 m |
0,018 m |
KLEP0026 | Ashtech Z-XII3 - Marine | 0,004 m |
0,002 m |
VLAS0120 | Leica SR299 | -0,001 m |
-0,001 m |
Korelace dat blízkých aparatur
Okamžité polohy současně GPS zaměřovaných bodů jsou v důsledku působení mnoha společných systematických vlivů více či méně silně korelovány. Těsnost korelací je závislá na rozličných faktorech, zejména na vzájemných vzdálenostech bodů a na jejich vzdálenosti od referenční stanice. Významné bývají na lokalitách menšího rozsahu, důležitým předpokladem jsou ovšem srovnatelné podmínky příjmu signálu. Na obr. 8 je vidět příklad korelací časového průběhu výšek dvou stanic SO, NS (znázorněné body grafu jsou řešení jednominutových dat vyhodnocených v čtvrthodinových intervalech), které byly počítány z referenční stanice TUBO vzdálené cca 20 km. Vzdálenost mezi oběma body je přibližně 18 m. Vypočtené koeficienty korelace pro daný případ byly 0,31 (L1 řešení), resp. 0,84 (L3 řešení).
obr.8
Korelace mezi výslednými parametry polohy jednotlivých bodu je možné využít ke zpřesnění výsledků GPS miření. K vyhodnocení je třeba použít data nejméně tří aparatur GPS. Jedna z nich je referenční (umístěná zpravidla na známém bodě), druhá je tzv. pseudoreferenční (bývá umisťována co nejblíže lokalitě měření) a třetí je vlastní měřící aparatura. Principem vylepšení výsledků je zavedení korekcí ze změn poloh pseudoreferenční aparatury pro aparaturu měřící. Zkušenosti ukazují, že hodnoty korekcí mohou dosahovat pro vertikální složku velikosti řádově až centimetrů, pro horizontální složku se pohybují většinou pod 1 cm. Tohoto způsobu zpřesnění výsledků bylo úspěšně použito při měření geometrického průběhu železniční koleje (Švábenský - Vitásek - Weigel , 1999).
Závěry
Využití korelace dat pro zavádění dodatečných korekcí prostorové polohy se ukazuje jako efektivní možnost zlepšování výsledků GPS kinematických a pseudokinematických metod. Pro správné určení korekcí je žádoucí zpracování dat totožné konstelace družic na fixní (pseudoreferenční) i pohyblivé stanici pro jednotlivé zaměřované epochy, což je možné dodatečně ošetřit při výpočtu.
Zpřesňovat výsledky lze dále prostřednictvím důsledné kalibrace antén GPS aparatur. Způsob kontroly antén s využitím metody Stop&Go se jeví jako vhodný, rychlý a realizovatelný i v polních podmínkách. Zavedení korekcí z anténních testovacích měření do jiných nezávislých měření při kombinaci různých typů aparatur přináší znatelné zlepšení výsledků.
Při použití různých atmosférických modelů mohou diference dosahovat zejména ve výškové složce několika centimetrů i pro relativně krátké délky několika kilometrů. Pro určitou lokalitu a zvolený způsob vyhodnocení. Zefektivnění postupu měření a zpřesnění výsledků lze rovněž dosáhnout kombinacemi několika kratších observací s vhodnými časovými rozestupy (6 až 8 hod.)
Tento referát byl zpracován v rámci řešených projektů CEZ: J22/98:261100007 a GAER 103/99/1633.
LITERATURA:
MADER, G. - Antenna Calibration Summary. Http://www.grdl.noaa.gov/GRD/GPS/Projects /ANTCAL/Files/summary.html
ROTHACHER M., MADER, G. - Combination of Antenna Phase Center Offsets and Variations, Antenna calibration set: IGS_01, FTP AIUB, 1996
ŠVÁBENSKÝ, O., VITÁSEK, J., WEIGEL, J. Měření geometrických parametrů koleje pomocí GPS. Sborník mezinárodní konference Geodézia a kartografia v doprave, Banská Štiavnica, 1999, (ISBN 80-968200-5-2), s. 39-44.
WEIGEL, J., ŠVÁBENSKÝ, O. Přesnost GPS výšek v lokálních sítích. Sborník (Seminář GPS a výšky) VUT v Brně, únor 1999, s. 55-62.
Předchozí příspěvek | Následující příspěvek |