| Tento článek má několik problémů. Prosím pomozte vylepši to nebo diskutovat o těchto problémech na internetu diskusní stránka. (Zjistěte, jak a kdy tyto zprávy ze šablony odebrat) | Téma tohoto článku nemusí splňovat požadavky Wikipedie obecný pokyn k notabilitě. Pomozte prosím určit notabilitu citováním spolehlivé sekundární zdroje to jsou nezávislý tématu a poskytnout jeho významné pokrytí nad rámec pouhé triviální zmínky. Pokud nelze určit významnost, je pravděpodobné, že článek bude sloučeny, přesměrovánnebo smazáno. Najít zdroje: "Cesty zemského řezu" – zprávy · noviny · knihy · učenec · JSTOR (Únor 2020) (Zjistěte, jak a kdy odstranit tuto zprávu šablony) |
| tento článek případně obsahuje původní výzkum. Prosím vylepši to podle ověřování vznesené nároky a přidání vložené citace. Výroky sestávající pouze z původního výzkumu by měly být odstraněny. (Leden 2020) (Zjistěte, jak a kdy odstranit tuto zprávu šablony) |
(Zjistěte, jak a kdy odstranit tuto zprávu šablony) |
Cesty zemského řezu jsou cesty na Zemi definované křižovatkou a referenční elipsoid a letadlo. Mezi běžné příklady pozemských řezů patří velká elipsa a normální řezy. Tato stránka poskytuje sjednocující přístup ke všem zemským sekcím a jejich přidruženým geodetické problémy.
Nepřímý problém
Nepřímým problémem pro zemské úseky je: daný dvěma body
a
na povrchu referenčního elipsoidu najděte délku,
, krátkého oblouku sféroidního úseku z
na
a také najít odjezd a příjezd (odkaz na skutečný sever) azimuty té křivky,
a
. Nechat
mít geodetickou šířku
a zeměpisná délka
(k = 1,2). Tento problém je nejlépe vyřešit pomocí analytická geometrie v ECEF souřadnice
a
být souřadnice ECEF dvou bodů, vypočítané pomocí diskutovaných geodetických transformací ECEF tady.
Rovina řezu
Chcete-li definovat rovinu řezu, vyberte libovolný třetí bod
není na řádku z
na
. Výběr
být na povrchu normální v
definuje normální sekci na
. Li
je původ, pak zemská část je velká elipsa. (Počátek by byl kolineární se 2 antipodálními body, takže v takovém případě musí být použit jiný bod). Protože existuje nekonečně mnoho možností
, výše uvedený problém je ve skutečnosti třídou problémů (jeden pro každé letadlo). Nechat
být dán. Chcete-li dát rovnici roviny do standardního tvaru,
, kde
, vyžaduje komponenty a jednotkový vektor,
, kolmé k rovině řezu. Tyto komponenty lze vypočítat následovně: Vektor z
na
má komponenty
a vektor z
na
má komponenty
. Proto,
=
×
), kde
je jednotkový vektor ve směru
. Zde je použita orientační konvence
ukazuje nalevo od cesty. Pokud tomu tak není, předefinujte
= -
. Nakonec lze parametr d pro rovinu vypočítat pomocí Tečkovaný produkt z
s vektorem od počátku do libovolného bodu v rovině, například
, tj. d =
. Rovnice roviny (ve vektorové podobě) je tedy
⋅
= d, kde
je vektor polohy z (x, y, z).
Azimut
Zkoumání transformace ENU na ECEF ukazuje, že souřadnice ECEF jednotkového vektoru směřujícího na východ v kterémkoli bodě elipsoidu jsou:
=
, jednotkový vektor směřující na sever je
=
a jednotkový vektor směřující nahoru je
=
. Vektor tangenta k cestě je:
takže východní složka
je
a severní složka je
. Azimut lze proto získat z a dvou argumentová arkustangensová funkce,
=
. Tuto metodu použijte u obou
a
dostat
a
.
Sekce elipsa
(Netriviální) průsečík roviny a elipsoidu je elipsa. Proto délka oblouku,
, na úseku cesty z
na
je eliptický integrál které lze vypočítat s jakoukoli požadovanou přesností pomocí zkrácené řady. Než to lze provést, musí být definována elipsa a vypočítány limity integrace.
a nechte
.Je-li p = 0, pak je řez vodorovným poloměrem
, který nemá řešení, pokud
.
Pokud p> 0, pak Gilbertson[1] ukázaly, že souřadnice ECEF středu elipsy jsou
, kde
,
poloviční hlavní osa je
, ve směru
, a poloviční vedlejší osa je
, ve směru
, který nemá řešení, pokud
.
Délka oblouku
Polární forma vzhledem ke středu pro rovnici elipsy je
, kde
, se týká excentricity elipsy, nikoli excentricity sféroidu (viz elipsa ). Nechť P je bod na elipsě a
, pak vektor z
na
má komponenty
. Pomocí argumentu podobného tomu pro azimut výše, pojďme
, pak
, a
, a
. Tímto způsobem získáme středové úhly
a
souhlasí s
a
resp. Je třeba dbát na to, aby to bylo zajištěno
≤
≤
. Pak délka oblouku podél elipsy je dán vztahem
=
Střídání
výše do tohoto vzorce, provedení označených operací, za použití jednoho dalšího výrazu než Gilbertsonova vyjádření a přeskupení, má za následek
, kde
![{ displaystyle { begin {aligned} Delta theta & = theta _ {2} - theta _ {1}, [6pt] Delta s2 & = sin (2 theta _ {2}) - sin (2 theta _ {1}), [6pt] Delta s4 & = sin (4 theta _ {2}) - sin (4 theta _ {1}), [6pt] Delta s6 & = sin (6 theta _ {2}) - sin (6 theta _ {1}), [6pt] {c_ {0}} & = 1 + e ^ {2} (4096 + 3328e ^ {2} + 2880e ^ {4}) / 16384, [6pt] {c_ {1}} & = e ^ {2} (512 + 384e ^ {2} + 380e ^ {4}) / 4096, [6pt] {c_ {2}} & = - e ^ {4} (64 + 80e ^ {2}) / 16384, [6pt] {c_ {3}} & = - 60e ^ { 6} / 12288. [6pt] end {zarovnáno}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c52c63c881c5545de1cf625b9e3249fae6c74a11)
Případně rozšíření pro Poledníkový oblouk lze zde použít nahrazením excentricity sféroidu excentricitou elipsy řezu.
Přímý problém
Přímý problém je uveden
, vzdálenost,
a odchodový azimut,
, najít
a příjezdový azimut,
.
Rovina řezu
Vytvořte tečný vektor v
,
, kde
a
jsou jednotkové vektory směřující na sever a na východ
. Vyberte vektor,
, definovat rovinu řezu, věnovat pozornost orientaci. Dodržujte to
nesmí být v rozsahu {
} (jinak by letadlo bylo tečné k Zemi v
, takže by nevznikla žádná cesta). Normální vektor
=
×
), dohromady s
definuje rovinu.
Lokalizovat 
Toto je 2-d problém v rozpětí {
}, který bude vyřešen pomocí výše uvedeného vzorce délky oblouku. Základním přístupem je použití Newton-Raphsonovy iterace k dosažení
. Základem odhadu je, že vektor polohy libovolného bodu na elipsě řezu lze vyjádřit pomocí vektoru polohy středu a středového úhlu jako
Chcete-li získat počáteční odhad
, nechť
,
= Central_Angle
(viz část délky oblouku výše),
,
.
Nyní inicializujte
=
a opakujte následující kroky:
![{ displaystyle { begin {aligned} s & = ArcLength ({ theta _ {1}}, { theta _ {2}}), [6pt] Err & = delta -s, [6pt] s '({ theta}) & = { frac {b ^ {*}} {(1-e ^ {2} cos ^ {2} theta)}} { sqrt { frac {(1- ( 2-e ^ {2}) e ^ {2} cos ^ {2} theta} {1-e ^ {2} cos ^ {2} theta}}}, [6pt] Delta theta & = { frac {Err} {s '({ theta _ {2}})}}, [6pt] theta _ {2} & = theta _ {2} + Delta theta, [6pt] end {zarovnáno}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/50e91c74a79053a8bb6803a76a34abf06158c7ee)
odejít, když 
Obvykle nejsou nutné více než tři iterace, i když téměř antipodální případy mohou být problematické. Nakonec
, a
= ECEF_to_Geo
pomocí Bowringova algoritmu z roku 1985,[2] nebo algoritmus tady.
Alternativně lze použít inverzi řady délky oblouku, aby se zabránilo iteracím.
Azimut
Azimut lze získat stejnou metodou jako nepřímý problém:
=
, kde dolní index 2 označuje vyhodnocení přidruženého množství při
.
Příklady
Velká elipsa
Nechat
být původ, takže
= vektor polohy
. Výše uvedený přístup poskytuje alternativu k přístupu ostatních, jako je Bowring.[3]
Normální sekce
Normální část v
je určeno nájmem
=
(povrch normální v
). Výše uvedený přístup poskytuje alternativu k přístupu ostatních, jako je Bowring.[4]
Střední normální část
Střední normální část z
na
je určeno pronajmutím
=
. To je dobrá aproximace geodetické z
na
pro letectví nebo plachtění.
Třída sekcí
Třídu sekcí lze představit otáčením
o spojování akordů
a
To vše lze vyřešit výše uvedeným jediným přístupem.
Křižovatky
Nechť jsou uvedeny dvě roviny řezu:
⋅
=
, a
⋅
=
. Za předpokladu, že dvě roviny nejsou rovnoběžné, je průsečík na obou rovinách. Proto ortogonální k oběma normálům, tj. Ve směru
.
Od té doby
a
nejsou kolineární
,
,
je základem pro
. Proto existují konstanty
a
tak, že přímka dvou rovin je dána vztahem
=
+
+ t
, kde t je nezávislý parametr.
Jelikož je tato čára v obou rovinách řezu, vyhovuje oběma:
+
(
·
) =
, a
(
·
) +
=
.
Řešení těchto rovnic pro
a
dává
[1 - (
] =
-
(
·
), a
[1 - (
] =
-
(
·
).
Definujte "úhel vzepětí",
tím, že
=
·
.Pak
=
, a
=
.
Na křižovatce, kterou máme
=
+ t
, kde
=
+
.Proto:
=
+ t
,
=
+ t
, a
=
+ t
, kde
=
+
,
=
+
, a
=
+
.a
=(
,
,
), pro i = 1,2,3.
Chcete-li najít průsečík této přímky se zemí, zapojte rovnice přímky
, dostat
, kde
=
,
=
,
=
.
Proto linka protíná Zemi v
. Li
, pak neexistuje žádná křižovatka. Li
, pak je čára tečná k Zemi v
(tj. řezy se protínají v tomto jediném bodě).
Dodržujte to
od té doby
a
nejsou kolineární. Zapojení do
=
+ t
, dává průsečíky pozemských řezů.
Příklady
Maximální nebo minimální zeměpisná šířka
na trase pozemského úseku lze nalézt zrušení indexů v daném úseku;
,
a nastavení
, aby
. Pak vyřešte pro
takhle
.
Od té doby
, a
, musíme mít
. Zapojení do
=
, dává průsečíky pozemských řezů. Případně stačí nastavit
.
Maximální nebo minimální délka
na trase pozemského úseku lze nalézt zrušení indexů v daném úseku;
,
a nastavení
, kde
je zeměpisná délka, kterou je třeba vyřešit
.
Případně stačí nastavit
.
Reference