Elipsoid

Elipsoid Elipsoid je omezená kvadratická plocha. Někdy se chápe jako prostorové těleso tvořené množinou všech bodů, jejichž poloha vůči zadanému bodu (středu) splňuje podmínky dané následující nerovnicí. Pokud bychom znak ≤ nahradili znakem =, rovnici by splňovaly právě body na povrchu elipsoidu.

{x^2 \over a^2}+{y^2 \over b^2}+{z^2 \over c^2} \leq 1

kde a, b a c jsou konstantní kladná reálná čísla, určující délky poloos ve směru jednotlivých os. Uvedená definice předpokládá, že střed elipsoidu leží v počátku soustavy souřadnic a že osy elipsoidu jsou totožné s osami soustavy souřadnic. +more Pokud tomu tak není, je třeba nerovnici rozšířit o popis posunutí a otočení elipsoidu v prostoru.

Rovinnými řezy elipsoidu podél jednotlivých souřadnicových os jsou elipsy. Poloosy jednotlivých elips odpovídají poloosám elipsoidu.

Klasifikace elipsoidů

Elipsoid, jehož dvě poloosy jsou shodné, se nazývá rotační elipsoid nebo také sferoid. Rotační elipsoid lze také chápat jako těleso vzniklé rotací elipsy kolem jedné z jejich os. +more Při rotaci kolem hlavní osy se jedná o rotační elipsoid protáhlý, při rotaci elipsy kolem vedlejší osy jde o zploštělý rotační elipsoid.

Elipsoid, který má shodné všechny tři poloosy, je koule.

Budeme-li předpokládat, že a ≥ b ≥ c, potom: * jestliže a > b > c, jde o obecný (trojosý) elipsoid * jestliže a > b = c, jde o protáhlý (doutníkovitý) sferoid * jestliže a = b > c, jde o zploštělý (diskovitý) sferoid * jestliže a = b = c, jde o kouli

Soubor:elipsoid_trojosy321.png|Trojosý elipsoid s poměrem poloos 3:2:1. Soubor:elipsoid_protahly.png|Protáhlý rotační elipsoid. Soubor:elipsoid_zplostely.png|Zploštělý rotační elipsoid.

Vlastnosti

Elipsoid je jediná omezená kvadrika. Rovinným řezem je elipsa, jednobodová nebo prázdná množina.

V projektivní geometrii je ekvivalentní eliptickému paraboloidu a dvojdílnému hyperboloidu.

Objem libovolného elipsoidu je roven :\frac{4}{3}\pi abc\,\!, kde a, b, c jsou délky poloos.

Speciální případy

Objem protáhlého rotačního elipsoidu je :V = \frac{4}{3}\pi ab^2, kde a>b.

Objem zploštělého rotačního elipsoidu je :V = \frac{4}{3}\pi a^2 b, kde a>b.

Povrch elipsoidu se počítá podstatně složitěji. S použitím číselné výstřednosti elipsy \varepsilon =\frac{\sqrt{a^2-b^2}}{a} lze povrch protáhlého rotačního elipsoidu zapsat ve tvaru :P = 2\pi \left(b^2 + a b \frac{\operatorname{arcsin}\varepsilon}{\varepsilon}\right) Pro zploštělý elipsoid pak dostaneme výraz :P = 2\pi \left(a^2 + {b^2} \frac{\operatorname{artanh}\varepsilon}{\varepsilon}\right) = 2\pi \left(a^2 + \frac{b^2}{2\varepsilon} \ln\frac{1+\varepsilon}{1-\varepsilon}\right)

Povrch trojosého elipsoidu se počítá složitěji pomocí eliptických integrálů.

Zobecněný hyperelipsoid

Typicky pracujeme s elipsoidem v trojrozměrném prostoru. Uvedenou definici lze zobecnit pro n-rozměrný prostor následovně (i nadále platí, že osy elipsoidu jsou totožné s osami soustavy souřadnic):

\sum_{i=1}^n {x_i^2 \over r_i^2} \leq 1

kde r_i je vektor poloos v jednotlivých rozměrech. Pro n>3 jde o abstraktní matematické těleso, které je obtížné vizualizovat. +more Pro n=2 jde o elipsu v rovině. Pro n=1 jde o úsečku na přímce.

Referenční elipsoid

Polární souřadnice

Klasické kartézské souřadnice neposkytují dostatečně intuitivní představu o poloze bodu na povrchu elipsoidu. Pokud víme, že se pohybujeme po povrchu, stačí nám totiž pouhé dvě souřadnice místo tří. +more Při využití elipsoidu jako referenčního tělesa pro povrch planety se proto obvykle používají polární (úhlové) souřadnice: zeměpisná šířka a zeměpisná délka. V mnoha případech nám totiž stačí pracovat s body na povrchu referenčního elipsoidu a odchylky způsobené nadmořskou výškou nebo lokální odchylkou gravitačního potenciálu zanedbáváme. Pokud bychom je nechtěli zanedbávat, museli bychom soustavu polárních souřadnic doplnit o třetí rozměr, výšku.

* Délka bodu A je úhel mezi rovinou xz (tj. rovinou obsahující osy x a z) a rovinou obsahující osu z a bod A. +more Označme tuto rovinu \delta. Vzhledem k tomu, že délka může být kladná i záporná, známe nejen úhel, ale i směr rotace od výchozí roviny (úhel je orientovaný) a máme navíc informaci o polorovině, na které se daný bod nachází (rovina \delta je rozdělena na poloroviny osou z). Na referenčním elipsoidu Země se délka určuje v rozsahu -180° až +180° (zápornou délku mají body, které mají zápornou kartézskou souřadnici y); u jiných těles se obvykle používá rozsah 0° až 360°. * Šířka je úhel mezi rovinou xy a přímkou r, která je podmnožinou roviny \delta a prochází středem elipsoidu i bodem A. Šířka se určuje v rozsahu -90° až +90° (zápornou šířku mají body, které mají zápornou kartézskou souřadnici z). * Výška je (případná) vzdálenost bodu A od povrchu elipsoidu a jde o reálné číslo, které je zdola omezeno hodnotou -a (opačné číslo k největší z poloos elipsoidu), shora omezené není. Jde-li o referenční elipsoid planety, udává se výška obvykle v metrech, popř. ve stopách. * Šířka geodetická je úhel mezi rovinou xy (tj. rovinou rovníku) a přímkou h. Přímka h je normálou k tečné rovině elipsoidu v daném bodě a obecně nemusí procházet středem elipsoidu. Pro počty s geodetickou šířkou musí platit poněkud jiné definice, než jaké jsou použity v následujících příkladech. Proto pozor, o kterou šířku se jedná. Viz např. standard WGS 84.

Převod polárních souřadnic na kartézské

Poznámka: Při nasazení níže uvedeného postupu v praxi je třeba si ověřit, zda software, s jehož pomocí počítáme funkci tangens, očekává úhel ve stupních, nebo v radiánech. Stupně lze přepočítat na radiány vzorcem radi\acute{a}ny = stupn\check{e} {\pi \over 180}.

Poznámka: Použití funkce tg nelze doporučit, v blízkosti singularity vede k chybám. Místo toho lze vystačit s funkcemi sin, cos a návod by bylo vhodné upravit.

# Mějme bod A daný délkou \lambda, šířkou \phi a výškou v. # y = x \mbox{tg } \lambda\,\. +more Zvláštním případem je \| \lambda \| = 90, kdy tangens není definován. Tehdy platí x=0\,\. # z = (\sqrt {x^2+y^2}) \mbox{tg } \phi, kde \sqrt {x^2 + y^2} je poloměr \phi-té rovnoběžky. Zvláštním případem je \| \phi \| = 90, kdy tangens není definován. Tehdy platí x=y=0\,\. # Nyní již můžeme popsat polopřímku r, která vychází ze středu elipsoidu a prochází bodem A. Výše uvedené vztahy mezi x, y a z vycházejí z úhlů definujících polopřímku a platí proto pro každý její bod. Na polopřímce tedy leží i bod A_1=[1;\mbox{tg }\lambda;(\sqrt{1+(\mbox{tg }\lambda)^2})\mbox{tg }\phi]. # Potřebujeme zjistit, v jaké vzdálenosti od středu elipsoidu leží bod, ve kterém polopřímka r protíná jeho povrch. Nejdříve tedy musíme zjistit souřadnice tohoto průsečíku. Hledáme takový skalár k_Z, že bod A_Z = k_Z \cdot A_1 leží na povrchu elipsoidu, tj. splňuje rovnici {{X_Z^2+y_Z^2}\over{a^2}}+{{Z_Z^2}\over{b^2}}=1. # Dosazením do této rovnice získáme k_Z = \sqrt{{A^2b^2}\over{b^2x_1^2+b^2y_1^2+a^2z_1^2}}. # Vzdálenost bodu A_Z od středu elipsoidu \|A_Z\|=\sqrt{x_Z^2+y_Z^2+z_Z^2}. # Nyní hledáme takové k, že pro bod A, jehož souřadnice hledáme, platí A=k\cdot A_1. k a k_Z jsou ve stejném poměru, v jakém jsou vzdálenosti A a A_Z od středu, tedy k = k_Z {{\|A_Z\|+v} \over {\|A_Z\|}}. Souřadnice bodu A pak už triviálně získáme pomocí k a souřadnic A_1.

Příklad

Mějme bod určený polárními souřadnicemi

* \phi = 50° s. š. = +50° = 0.872664626 radiánů * \lambda = 14° 30' v. d. = +14,5° = 0.253072741 radiánů * v = 250 m

Předpokládejme, že polární souřadnice jsou vztaženy k povrchu referenčního elipsoidu WGS-84, tedy

* a = 6 378 137 m * b = 6 356 752,3 m

Potom

* A_1 = [1; 0.258618; 1.230963] * k_Z = 3 961 374 * A_Z = [3961374; 1024481; 4876303] * \|A_Z\| = 6365562 * {{\|A_Z\|+v} \over {\|A_Z\|}} = 1.000039 * A = [3961529; 1024521; 4876495]

Tedy

* x = 3 961 529 m * y = 1 024 521 m * z = 4 876 495 m

Převod kartézských souřadnic na polární

Poznámka: Použití funkce arctg je chybné, mj. nerespektuje znaménko výsledku. +more Návod by bylo vhodné upravit, ale bude složitější. Inspirací může být implementace funkce atan2 systému Matlab.

Předpokládejme, že pracujeme se zploštělým sferoidem, např. s referenčním elipsoidem Země. Pro obecný elipsoid lze postup zobecnit.

# Mějme bod A a jeho kartézské souřadnice x, y a z. # Polopřímku r, která vychází ze středu elipsoidu a prochází bodem A, lze parametricky vyjádřit jako [0;0;0]+k \cdot A, kde k>0. +more # Protože bod A může mít nenulovou výšku a nemusí tedy ležet přímo na povrchu elipsoidu, potřebujeme najít bod A_Z, který leží na průniku polopřímky r s povrchem elipsoidu. # Tento bod musí současně splňovat rovnici pro polopřímku i rovnici pro povrch elipsoidu {{X_Z^2+y_Z^2} \over {a^2}}+{{Z_Z^2} \over {b^2}}=1. # Vyjádříme-li souřadnice A_Z pomocí k_Z a souřadnic A, dostaneme {{K_Z^2x^2+k_Z^2y^2} \over {a^2}}+{{K_Z^2z^2} \over {b^2}}=1. # {k_Z^2}={{A^2b^2}\over{b^2x^2+b^2y^2+a^2z^2}} # k_Z=\pm\sqrt{{A^2b^2}\over{b^2x^2+b^2y^2+a^2z^2}} # Dostaneme 2 koeficienty, jeden kladný a jeden záporný, protože přímka protne povrch elipsoidu na dvou místech. Jak už ale bylo zmíněno, zajímá nás pouze polopřímka vedoucí ze středu elipsoidu směrem k bodu A, využijeme tedy pouze kladný z obou koeficientů. # S pomocí koeficientu k dostaneme souřadnice bodu A_Z=k_Z\cdot[x;y;z]. # Délka vektoru z bodu A_Z do bodu A odpovídá hledané nadmořské výšce: v=\pm\sqrt{(x-x_Z)^2+(y-y_Z)^2+(z-z_Z)^2}. Směr vektoru nám prozradí, zda má být výška kladná, nebo záporná. Totéž lze vyčíst i z koeficientu k - pokud je větší než 1, leží bod A_Z na polopřímce r až za bodem A, a výška bodu A je tedy záporná. V opačném případě je výška kladná. # Délka \lambda=\mbox{tg}^{-1}{y\over{x}} # Šířka \phi=\mbox{tg }^{-1}{z\over\sqrt{x^2+y^2}}.

Příklad

Mějme bod A určený kartézskými souřadnicemi

* x = 1 113 547 m * y = -4 823 300 m * z = 4 008 571 m

Hledáme polární souřadnice tohoto bodu vztažené k povrchu referenčního elipsoidu WGS-84, tedy

* a = 6 378 137 m * b = 6 356 752,3 m

Potom

* \phi = 39 * \lambda = -77 * k_Z = 0,999994 * A_Z = [1113540; -4823270; 4008546] * \|A_Z\| = 6365562 * A-A_Z = [6,992792; -30,2891; 25,17282] * \|A-A_Z\| = 40

Tedy

* φ = 39° s. š. * λ = 77° z. d. * v = 40 m

(Pohledem na mapu pak zjistíme, že tento bod leží na předměstí Washingtonu.)

Vzdálenost dvou bodů na povrchu elipsoidu

Vzdálenost dvou bodů na povrchu (nebo v blízkosti povrchu) elipsoidu nás zajímá zejména u referenčního elipsoidu, chceme-li určit vzdálenost dvou bodů na povrchu Země či jiné planety.

Po přímce

Nejkratší vzdáleností dvou bodů v prostoru je délka úsečky, kterou tyto body vymezují, tedy vzdálenost měřená na přímce. Pro praktické užití na referenčním elipsoidu má ovšem takto zjištěná vzdálenost zásadní nevýhodu, totiž že spojnice obou bodů typicky prochází pod povrchem elipsoidu a neodpovídá tedy vzdálenosti, kterou musí při přesunu mezi oběma body urazit např. +more letadlo. Pouze u malých vzdáleností je tento rozdíl zanedbatelný.

Vzdálenost bodů A a B lze přímo vypočítat z jejich kartézských souřadnic, stačí tedy převést polární souřadnice na kartézské podle výše uvedeného postupu.

d=\sqrt{(x_B-x_A)^2+(y_B-y_A)^2+(z_B-z_A)^2}

Příklad

Mějme body A (Praha) a B (Washington) jako ve výše uvedených příkladech, tedy

* A = [3961529; 1024521; 4876495] * B = [1113547; -4823300; 4008571]

Potom

* (B-A) = [-2847982,876; -5847821,155; -867924,0585] * d = 6 562 112 m

Po elipse

Vzdálenost vedená po povrchu elipsoidu by byla délka úseku elipsy, která prochází oběma body a má střed shodný se středem elipsoidu. Tato elipsa je průnikem povrchu elipsoidu a roviny definované středem elipsoidu a oběma body, jejichž vzájemnou vzdálenost měříme. +more Problém je, že neznáme poloosy této elipsy; víme pouze, že leží někde mezi největší a nejmenší poloosou elipsoidu. Můžeme však aspoň provést horní a dolní odhad vzdálenosti, pokud budeme místo elipsy postupně uvažovat kružnice s poloměry odpovídajícími největší a nejmenší poloose elipsoidu. Tento postup navíc nezohledňuje případnou výšku bodů A a B nad povrchem elipsoidu.

Při cestování po dráze dané elipsou se středem uprostřed elipsoidu se soustavně mění azimut, tedy úhel mezi směrem pohybu a severem. Přesun z bodu A do bodu B lze také provést po loxodromě (anglicky rhumb line), tedy čáře protínající všechny poledníky pod stejným úhlem; taková trasa však bude delší než trasa po elipse. +more I tento rozdíl je pro blízké body zanedbatelný.

* Pro změření úseku kružnice mezi body A a B potřebujeme znát úhel \gamma mezi přímkami a a b, přičemž přímka a prochází bodem A, přímka b prochází bodem B a obě procházejí středem kružnice. * Velikost úhlu \gamma dokážeme určit pomocí goniometrických funkcí, pokud budeme znát polohu bodu K, ve kterém přímka c, spojující body A a B, protíná kolmici na ni vedenou ze středu elipsoidu S. +more ** Přímku c lze popsat jako C=A+k\cdot (B-A), kde k je reálné číslo a (B-A) je vektor vedoucí z bodu A do bodu B. Pro každý bod C přímky c existuje skalár k takový, že uvedená rovnice platí. Náš bod K má ležet na přímce c, budeme tedy hledat jemu odpovídající k. ** Přímku v, která prochází body S a K, lze popsat jako V=S+l\cdot (K-S). Protože bod S je střed elipsoidu a počátek soustavy souřadnic, platí S=[0;0;0], a tedy V=l\cdot K. ** Protože dále v má být kolmá na c, musí platit, že skalární součin jejich směrových vektorů je roven nule. Naším směrovým vektorem pro přímku c je (B-A), pro přímku v je to (zatím neznámé) K, tedy x_{(B-A)}x_K+y_{(B-A)}y_K+z_{(B-A)}z_K=0\,\. ** Dosadíme-li do právě uvedené rovnice vyjádření bodu K jako prvku přímky c, dostaneme rovnici o jedné neznámé k: x_{(B-A)}(x_A+kx_{(B-A)})+y_{(B-A)}(y_A+ky_{(B-A)})+z_{(B-A)}(z_A+kz_{(B-A)})=0\,\. ** Úpravou rovnice dostaneme k=-{{X_A(x_B-x_A)+y_A(y_B-y_A)+z_A(z_B-z_A)}\over{(x_B-x_A)^2+(y_B-y_A)^2+(z_B-z_A)^2}}. ** Získané k dosadíme do rovnice K=A+k\cdot (B-A) a dostaneme souřadnice bodu K. * Změříme úsečky AK, KB a KS. Např. \|AK\|=\sqrt{(x_K-x_A)^2+(y_K-y_A)^2+(z_K-z_A)^2}. * Vypočítáme úhel \gamma=\mbox{tg}^{-1}{{\|AK\|}\over{\|KS\|}}+\mbox{tg}^{-1}{{\|KB\|}\over{\|KS\|}}. * Jestliže jsme úhel \gamma dostali ve stupních, převedeme ho na radiány: radi\acute{a}ny = stupn\check{e} {\pi \over 180}. * Délku úseku kružnice vypočítáme jako d=\gamma r\,\. , kde r je poloměr kružnice.

Příklad

Mějme body A (Praha) a B (Washington) jako ve výše uvedených příkladech, tedy

* A = [3961529; 1024521; 4876495] * B = [1113547; -4823300; 4008571]

a poloosy referenčního elipsoidu

* a = 6 378 137 m * b = 6 356 752,3 m

Potom

* (B-A) = [-2847982,876; -5847821,155; -867924,0585] * k = 0,499427705 * K = [2539167,869; -1896042,727; 4443029,462] * |AK| = 3 277 300,419 m * |KB| = 3 284 811,352 m * |SK| = 5 457 367,707 m * γ = 1,082625728 rad = 62,03°

Pro kružnice s maximálním (a) a minimálním (b) poloměrem pak získáme horní a dolní odhad vzdálenosti po elipse:

* da = 6 905 135 m * db = 6 881 984 m

Jak vidíme, chyba způsobená tím, že se pohybujeme po kružnici místo po elipse, je relativně malá (23 km), zatímco rozdíl mezi vzdáleností po kružnici a vzdáleností po přímce je poměrně značný (320, resp. 343 km).

Odkazy

Související články

Sféroid * Geoid * Geometrický útvar * Kvadrika (kvadratická plocha) * Koule * Elipsa

Externí odkazy

Kategorie:Oblá tělesa Kategorie:Navigace Kategorie:Kartografie Kategorie:Plochy