Elipsoid

uzavřená kvadratická plocha

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.

Elipsoid

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. 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ů editovat

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. 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

Vlastnosti editovat

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

 ,

kde a, b, c jsou délky poloos.

Speciální případy editovat

Objem protáhlého rotačního elipsoidu je

 ,

kde  .

Objem zploštělého rotačního elipsoidu je

 ,

kde  .

Povrch elipsoidu se počítá podstatně složitěji. S použitím číselné výstřednosti elipsy   lze povrch protáhlého rotačního elipsoidu zapsat ve tvaru

 

Pro zploštělý elipsoid pak dostaneme výraz

 

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

Zobecněný hyperelipsoid editovat

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):

 

kde   je vektor poloos v jednotlivých rozměrech. Pro   jde o abstraktní matematické těleso, které je obtížné vizualizovat. Pro   jde o elipsurovině. Pro   jde o úsečku na přímce.

Referenční elipsoid editovat

Související informace naleznete také v článku Referenční elipsoid.

Polární souřadnice editovat

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ří. 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   je úhel mezi rovinou   (tj. rovinou obsahující osy   a  ) a rovinou obsahující osu   a bod  . Označme tuto rovinu  . 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   je rozdělena na poloroviny osou  ). 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  ); u jiných těles se obvykle používá rozsah 0° až 360°.
  • Šířka je úhel mezi rovinou   a přímkou  , která je podmnožinou roviny   a prochází středem elipsoidu i bodem  . Šířka se určuje v rozsahu -90° až +90° (zápornou šířku mají body, které mají zápornou kartézskou souřadnici  ).
  • Výška je (případná) vzdálenost bodu   od povrchu elipsoidu a jde o reálné číslo, které je zdola omezeno hodnotou   (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   (tj. rovinou rovníku) a přímkou  . Přímka   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é editovat

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  .

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.

  1. Mějme bod   daný délkou  , šířkou   a výškou  .
  2.  . Zvláštním případem je  , kdy tangens není definován. Tehdy platí  .
  3.  , kde   je poloměr  -té rovnoběžky. Zvláštním případem je  , kdy tangens není definován. Tehdy platí  .
  4. Nyní již můžeme popsat polopřímku  , která vychází ze středu elipsoidu a prochází bodem  . Výše uvedené vztahy mezi  ,   a   vycházejí z úhlů definujících polopřímku a platí proto pro každý její bod. Na polopřímce tedy leží i bod  .
  5. Potřebujeme zjistit, v jaké vzdálenosti od středu elipsoidu leží bod, ve kterém polopřímka   protíná jeho povrch. Nejdříve tedy musíme zjistit souřadnice tohoto průsečíku. Hledáme takový skalár  , že bod   leží na povrchu elipsoidu, tj. splňuje rovnici  .
  6. Dosazením do této rovnice získáme  .
  7. Vzdálenost bodu   od středu elipsoidu  .
  8. Nyní hledáme takové  , že pro bod  , jehož souřadnice hledáme, platí  .   a   jsou ve stejném poměru, v jakém jsou vzdálenosti   a   od středu, tedy  . Souřadnice bodu   pak už triviálně získáme pomocí   a souřadnic  .

Příklad editovat

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

  •   = 50° s. š. = +50° = 0.872664626 radiánů
  •   = 14° 30' v. d. = +14,5° = 0.253072741 radiánů
  •   = 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

  •  
  •  
  •  
  •  
  •  
  •  

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í editovat

Poznámka: Použití funkce arctg je chybné, mj. nerespektuje znaménko výsledku. 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.

  1. Mějme bod   a jeho kartézské souřadnice  ,   a  .
  2. Polopřímku  , která vychází ze středu elipsoidu a prochází bodem  , lze parametricky vyjádřit jako  , kde  .
  3. Protože bod   může mít nenulovou výšku a nemusí tedy ležet přímo na povrchu elipsoidu, potřebujeme najít bod  , který leží na průniku polopřímky   s povrchem elipsoidu.
  4. Tento bod musí současně splňovat rovnici pro polopřímku i rovnici pro povrch elipsoidu  .
  5. Vyjádříme-li souřadnice   pomocí   a souřadnic  , dostaneme  .
  6.  
  7.  
  8. 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  , využijeme tedy pouze kladný z obou koeficientů.
  9. S pomocí koeficientu   dostaneme souřadnice bodu  .
  10. Délka vektoru z bodu   do bodu   odpovídá hledané nadmořské výšce:  . Směr vektoru nám prozradí, zda má být výška kladná, nebo záporná. Totéž lze vyčíst i z koeficientu   - pokud je větší než 1, leží bod   na polopřímce   až za bodem  , a výška bodu   je tedy záporná. V opačném případě je výška kladná.
  11. Délka  
  12. Šířka  

Příklad editovat

Mějme bod   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

  •  
  •  
  •  
  •  
  •  
  •  
  •  

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 editovat

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 editovat

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ř. letadlo. Pouze u malých vzdáleností je tento rozdíl zanedbatelný.

Vzdálenost bodů   a   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.

 

Příklad editovat

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 editovat

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. 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. I tento rozdíl je pro blízké body zanedbatelný.

  • Pro změření úseku kružnice mezi body   a   potřebujeme znát úhel   mezi přímkami   a  , přičemž přímka   prochází bodem  , přímka   prochází bodem   a obě procházejí středem kružnice.
  • Velikost úhlu   dokážeme určit pomocí goniometrických funkcí, pokud budeme znát polohu bodu  , ve kterém přímka  , spojující body   a  , protíná kolmici na ni vedenou ze středu elipsoidu  .
    • Přímku   lze popsat jako  , kde   je reálné číslo a   je vektor vedoucí z bodu   do bodu  . Pro každý bod   přímky   existuje skalár   takový, že uvedená rovnice platí. Náš bod   má ležet na přímce  , budeme tedy hledat jemu odpovídající  .
    • Přímku  , která prochází body   a  , lze popsat jako  . Protože bod   je střed elipsoidu a počátek soustavy souřadnic, platí  , a tedy  .
    • Protože dále   má být kolmá na  , musí platit, že skalární součin jejich směrových vektorů je roven nule. Naším směrovým vektorem pro přímku   je  , pro přímku   je to (zatím neznámé)  , tedy  .
    • Dosadíme-li do právě uvedené rovnice vyjádření bodu   jako prvku přímky  , dostaneme rovnici o jedné neznámé  :  .
    • Úpravou rovnice dostaneme  .
    • Získané   dosadíme do rovnice   a dostaneme souřadnice bodu  .
  • Změříme úsečky AK, KB a KS. Např.  .
  • Vypočítáme úhel  .
  • Jestliže jsme úhel   dostali ve stupních, převedeme ho na radiány:  .
  • Délku úseku kružnice vypočítáme jako  , kde   je poloměr kružnice.

Příklad editovat

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 editovat

Související články editovat

Externí odkazy editovat