ステレオ投影

提供: testwiki
ナビゲーションに移動 検索に移動
地球を、北極を接点とする平面に、ステレオ投影した図
球面を球面の下の平面に、北極からステレオ投影する3次元の説明図

ステレオ投影(ステレオとうえい、テンプレート:Lang-en-short)は、球面平面投影する方法の一つである。ステレオ投影は複素解析学地図学結晶学写真術など様々な分野で重要である。

stereographic projection の訳語は分野によって異なる。ステレオ投影は主に物理学機械工学において用いられる。数学においては写像という意味で立体射影あるいはステレオグラフ射影、地図学では図法という意味で平射図法またはステレオ図法と呼ばれる。このように訳語が異なってはいるが、内容は全て同一視できる。

ステレオ投影は、数学的には写像として定義される。定義域は、球面から光源の一点を除いたところである。写像は滑らかかつ全単射である。また、等角写像、すなわち角度が保存される。一方、長さや面積は保存されない。これはとくに光源点付近では顕著である。

すなわち、ステレオ投影は、いくらかの避けられない妥協を含む、球面を平面に描く方法である。実際面では、コンピュータや、ウルフネットまたはステレオネットと呼ばれるテンプレート:仮リンクなどを使って、投影図が描かれる。

歴史

ステレオ投影は、ヒッパルコスクラウディオス・プトレマイオスに知られていたが、おそらくもっと早くから古代エジプトでも知られていた。これはもともと、平球投影(英語:planisphere projection[1])として知られていた[2]。プトレマイオスの著書 "テンプレート:仮リンク" は、ステレオ投影についてかかれた現存する最古の文書である。この投影の最も重要な使い方は星図を表すことであった。[2] 星座早見盤の英語 planisphere のように、今でもそのような図にこの言葉が使われる。

最初の世界地図は、1507年にグアルテリアス・ラドによって、ステレオ投影を基にそれぞれの半球を円盤に投影して描かれたと言われている[3]。ステレオ投影の赤道面への投影するという特徴は、17世紀18世紀東半球西半球の地図を描くのに利用された[4]

フランソワ・デギュイヨンが、彼の1613年の作品"Opticorum libri sex philosophis juxta ac mathematicis utiles"(哲学者と数学者に等しく役立つ6冊の光学の本)で、この投影にステレオ投影と名付けた[5]

数学的な定義

単位球面の北極から平面 テンプレート:Math への立体射影を表した断面図テンプレート:Mvar の像が テンプレート:Mvar である。

この節では、単位球面を北極から赤道を通る平面に投影する場合を扱う。その他の場合は後述

三次元空間 テンプレート:Math 内の単位球面は、テンプレート:Math と表すことができる。ここで、点 テンプレート:Math を"北極"とし、テンプレート:Mvar を球面の残りの部分とする。平面 テンプレート:Math は球の中心を通る。"赤道" はこの平面と球面の交線である。

テンプレート:Mvar 上の任意の点 テンプレート:Mvar に対して、テンプレート:Mvarテンプレート:Mvar を通る直線が一意的に存在し、この直線は平面 テンプレート:Math とちょうど一点 テンプレート:Mvar で交わる。テンプレート:Mvar立体射影による像を、その点 テンプレート:Mvar と定義する。

球面上の直交座標 テンプレート:Math と平面上の テンプレート:Math を用いると、立体射影とその逆写像は、次の式で与えられる。

(X,Y)=(x1z,y1z),
(x,y,z)=(2X1+X2+Y2,2Y1+X2+Y2,1+X2+Y21+X2+Y2).

球面上のテンプレート:仮リンク テンプレート:Math と平面上の極座標 テンプレート:Math を用いると、立体射影とその逆写像は

(R,Θ)=(sinφ1cosφ,θ),
(φ,θ)=(2arctan(1R),Θ)

となる。ただし、テンプレート:Math の場合は、テンプレート:Math と解釈する。

また、三角関数の等式を用いて、この式を書き直す方法がたくさんある。球面上の円柱座標 テンプレート:Math と平面上の極座標 テンプレート:Math を用いると、立体射影とその逆写像は、

(R,Θ)=(r1z,θ),
(r,θ,z)=(2R1+R2,Θ,R21R2+1)

となる。

性質

前の節で定義された立体射影は、“南極” テンプレート:Mathテンプレート:Math に、“赤道” を単位円に、南半球をその円の内側に、北半球をその円の外側に射影する。

この変換は テンプレート:Math では定義されない。この点の小さい近傍は、平面の テンプレート:Math から遠く離れた部分集合に射影される。テンプレート:Mvarテンプレート:Math に近づけるほど、像は テンプレート:Math より遠くなる。この性質により、一般に テンプレート:Math は射影によって平面の無限遠点に写ると言われ、また球面は “無限遠点” を付け加えて平面を完備化しているといわれる。この事実が、射影幾何学複素解析において有用である。単に位相的なレベルでは、それは球面が平面の一点コンパクト化同相である様子を示している。

平面上の10×10の正方格子は、球の上では歪んでいる。格子線は射影後も互いに垂直だが、格子内の正方形の面積は北極に近いほど縮んでいる。
平面上の半径 5 の極格子は、球の上で歪んでいる。格子曲線は射影後も互いに垂直だが、格子内の面積は北極に近いほど縮んでいる。

立体射影は等角写像、つまり曲線同士が互いに交わる角度を保存する(図参照)。しかし、面積は保存しない。一般に球面の領域の面積は平面に射影した面積に等しくない。テンプレート:Mvar-座標での面積要素は次のようになる: dA=4(1+X2+Y2)2dXdY. 単位円 テンプレート:Math では、面積の歪みは見られない。テンプレート:Math の近くでは4倍に歪められ、無限遠点の近くではいくらでも小さくなる。

角度と面積をともに保存する球面から平面への写像は存在しない。仮にあるとすれば、それは局所等長写像で、ガウス曲率を保存しているはずである。しかし球面と平面はガウス曲率が異なるので、これは不可能である。

立体射影の等角性は、いくつかの有用な幾何学的性質を意味する。光源点を「通らない」球面上の円は、平面上の円に射影される。光源点を「通る」球面上の円は、平面上の直線に射影される。このような直線は、無限遠点を通る円や無限大の半径を持つ円とみなされる事がある。

平面上のすべての直線は、立体射影の逆写像により球面上の円に写されると、光源点(無限遠点)で交わるようになる。平面上の平行線は、平面上では交わることはないが、球面上に投影されると光源点で接する。このように、平面上のすべての直線は、球面上のどこかで交わる—2点で横断的に交わるか、あるいは無限遠点で接するかである。(似たような性質は実射影平面でも成り立つが、交点の関係はそれとは違っている。)


球面上の等角航路を平面上に写した曲線は次の式で表される。

R=eθ/tan(B),

ここで テンプレート:Math は、等角航路の方位角である。よって等角航路は、等角螺旋になる。等角航路が常に子午線と同じ角度で交わるのと同じように、この螺旋は平面上の放射線(子午線の投影像)に常に同じ角度で交わる。

ウルフネット

ウルフネット(またはステレオネット)。手でステレオ投影図を書くときに使われる。

ステレオ投影は、前節で与えたような数式を用いて、コンピューターに算出させることができる。しかし、手でグラフを描くには、これらの数式は扱いにくい。代わりに、目的に特化して設計されたグラフ用紙を使うのが一般的である。このグラフ用紙を作るには、半球上に緯線と経線の格子を置き、それらを円盤の上にステレオ投影した曲線を描く。これをステレオネットまたはウルフネットという。ウルフネットの名は、ロシア人鉱物学者であるジョージ・ウルフ(George (Yuri Viktorovich) Wulff)[6] に由来する。

このネットの中心付近の区画と端の方の区画を比べることで、ステレオ投影の面積が歪む特性を見ることができる。この二つの区画は球面上では同じ面積を持っていた。円盤上では、端の方の区画は中心の区画の4倍近くの面積を面積を持っている。球の上の格子の目が細かければ細かいほど、その面積の比は4倍に近づく。

格子線を見ることで、この投影の正角性も見ることが出来る。球面上の緯線と経線は直角に交わるが、ウルフネット上のそれらの像も同じく直角で交わっている。

ウルフネットに点を描くときの書き方

ウルフネットの使い方の例を説明する。 まずウルフネットが描かれた薄い紙を2枚用意し、片方をもう片方の上に重ね、互いの中心を揃えて鋲で止める。 仮に、下側の単位半球面上の点(0.321, 0.557, -0.766)を描くとする。この点はx軸正の方向から60° 反時計回り(またはy軸正の方向から30° 時計回り)の方向にあり、z = 0の水平面より50° 下である。これらの角度を知れば、次の4 ステップで描ける。

  1. ここでは格子の間隔は10° である。格子を使い、点(1, 0)から60° 反時計回り(または点(0, 1)から30° 時計回り)のネットの端に印を付ける。
  2. 上側のネットを回し、付けた印を下側のネットの点(1, 0)に合わせる。
  3. 下側の格子を使い、印を付けた点から中心に向かって50° のところに、点を打つ。
  4. 上側のネットをさっきとは逆方向に回し、下側のネットと揃えた元の位置に戻す。ステップ3 で打った点が、目標の点のステレオ投影である。

60° や 50° のような切りのいい数字ではない角度の点を描くには、近い格子の間を補間しなければならない。10°より目の細かいネットの方が使いやすく、格子間隔が2° のものが一般的である。

このステレオ投影点を元にして2点間の球面上の中心角を見つけるには、ウルフネットをその上に被せて、2点が同じ経線に載るまたは近くなるまで互いの中心を合わせて回す。その経線に沿って格子線を数えることで、中心角を測ることができる。

別な定式化と一般化

単位球面の北極から平面 テンプレート:Math の上への立体射影、ここではその断面を示す。

北極 テンプレート:Math から南極 テンプレート:Math における接平面 テンプレート:Math への立体射影を定義する文献もある[7]。この射影によって得られる座標 テンプレート:Mvar は前記の節で記述した赤道射影 (equatorial projection) のちょうど二倍の値を与える。例えば、この射影により赤道は原点を中心とする半径 テンプレート:Math の円に写される。赤道射影の場合に赤道面に沿った無限小面積で歪みが無くなるのに対して、この極-接平面射影の場合では南極における無限小面積が歪んでいない。

別な文献 テンプレート:Harvtxt では半径 テンプレート:Math の球面と平面 テンプレート:Math で考える。この場合の定義式は

(x,y,z)(ξ,η)=(x12z,y12z),
(ξ,η)(x,y,z)=(ξ1+ξ2+η2,η1+ξ2+η2,1+ξ2+η22+2ξ2+2η2)

となる。

テンプレート:Math から平面 テンプレート:Mvar への立体射影、ここではその断面を示す。

一般に、球面上の任意の点 テンプレート:Math から

という条件を満たす任意の平面 テンプレート:Mvar の上への立体射影を定義することができる。テンプレート:Mvar がこれら条件を満たす限りにおいて、テンプレート:Math でない任意の点 テンプレート:Math に対して テンプレート:Math を通る直線と テンプレート:Mvar との交わりは、テンプレート:Mathテンプレート:Mvar の上への立体射影として定義されるべき点 テンプレート:Math ただ一つである[8]

ここまでで述べた立体射影の定式化は何れも本質的に同じ性質を有する。これらは何れも射影点を除き至る所定義された滑らかな全単射(微分同相)であり、共形(等角)であって面積を保存しない。

より一般に、テンプレート:Math-次元ユークリッド空間 テンプレート:Math 内の テンプレート:Mvar-次元超球面 テンプレート:Mvar に対して立体射影を考えることができる。テンプレート:Mvar の点 テンプレート:Mathテンプレート:Math 内の超平面 テンプレート:Mvar に対して、点 テンプレート:Math の立体射影 テンプレート:Math は直線 テンプレート:Mathテンプレート:Mvar との交点とする。球面上の直交座標系 テンプレート:Math および平面上の直交座標系 テンプレート:Math に関して、テンプレート:Math からの射影は

Xi=xi1x0(i=1,,n)

で与えられる。テンプレート:Math と置けば、逆写像は

x0=s21s2+1,xi=2Xis2+1(i=1,,n)

となる。さらにより一般に、テンプレート:Mvar射影空間 テンプレート:Math 内の(非特異)二次超曲面とする。つまり、テンプレート:Mvarテンプレート:仮リンク テンプレート:Mvar に関する非特異二次形式 テンプレート:Math の零点の軌跡である。テンプレート:Mvar 上の任意の点 テンプレート:Mathテンプレート:Math を含まない テンプレート:Math 内の超平面 テンプレート:Mvar を固定すれば、点 テンプレート:Math の立体射影は、直線 テンプレート:Mathテンプレート:Mvar との唯一の交点と定められる。前と同じく、この立体射影は共形かつ「小さい」集合の外側で可逆である。この立体射影は二次超曲面を有理超曲面として与えるものである[9]。この構成は代数幾何学およびテンプレート:仮リンクにおいて役割を果たす。

数学での応用

複素解析

テンプレート:Main 任意の立体射影において球面上の一点(射影点)は欠けてしまうけれども、相異なる二点に関する二つの立体射影を用いれば全球を写すことができる。言葉を変えれば、球面は平面からの二つの立体図式テンプレート:仮リンク(立体射影の逆写像)によって被覆することができる。二つの媒介変数表示は球面上で同じ向きを持つように選ぶことができる。併せて、これらは球面を向き付けられた曲面(二次元多様体)として記述する。

この構成は複素解析において特に著しい意義を持つ。実平面上の点 テンプレート:Math複素数 テンプレート:Math と同一視すれば、北極から赤道面への立体射影は

ζ=x+iy1z,
(x,y,z)=(2e(ζ)1+ζ¯ζ,2m(ζ)1+ζ¯ζ,1+ζ¯ζ1+ζ¯ζ)

となる。同様に、別の複素座標 テンプレート:Math を取れば

ξ=xiy1+z,
(x,y,z)=(2e(ξ)1+ξ¯ξ,2m(ξ)1+ξ¯ξ,1ξ¯ξ1+ξ¯ξ)

が南極から赤道面への立体射影を定義する。テンプレート:Mvar-座標から テンプレート:Mvar-座標への遷移写像は テンプレート:Math および テンプレート:Math で与えられ、テンプレート:Math のとき テンプレート:Math であり、逆もまた然りである。これにより複素数に対するエレガントで有用な無限大の概念がえられ、これは実際に有理型函数の理論の全容をリーマン球面へ引き写すものになっている。単位球面上の標準リーマン計量はリーマン球面上のフビニ–スタディ計量に一致する。

平面と直線の可視化

面心立方格子の8つあるうちの4つの<111>結晶面の間で傾きを変えていくアニメーション。平面(この図の中では束になった曲線)は常に同じ角度で交わる。.

3次元空間の原点を通るすべての直線の集合は実射影平面と呼ばれる空間をなす。この空間は、3次元空間に埋め込む事が出来ないので、可視化するのが難しい。

しかし以下のようにして、それを円板として「ほぼ」可視化することが出来る。原点を通るすべての直線は、南半球 z ≤ 0 と交わる。その点は立体射影により円板上の点に射影することができる。水平な直線は、赤道上にある2点で南半球に交わる。そのどちらの点も、この円板に投影することができる。これに関しては、円板の境界上の対蹠点は、同一の直線を表すとして理解する(商位相空間を参照)。よって、原点を通る直線の任意の集合が、円板内の点の集合として、ほぼ完全に投影される。

また、原点を通る各々の平面は単位球面と大円(これを平面のtraceと呼ぶ)で交わる。この円は立体射影により、円に投影される。よって、この投影は、円板中の円弧として、平面を可視化することが出来る。コンピューターが普及する前には、大円のステレオ投影をするには、ビームコンパスを使って書く必要があるほど大きな半径の弧を書くことがあった。現在はコンピューターにより簡単にそれらを描くことができる。

さらに、任意の平面は、原点を通りその平面に垂直な直線を唯一つ持つ。その直線を平面の極線と呼ぶ。その直線は、その他の原点を通る直線と同様に円板の上に点を打つことができる。よって、立体射影は任意の平面も円板内の点として可視化することができる。たくさんの平面をプロットするには、平面の極線をプロットしていく方が、平面の軌跡をプロットしていくより、整然とした図を作ることができる。

この様な作業は、後で述べるように、結晶学や地学で方向のデータを可視化するのに使われる。

その他可視化

正十二面体のシュレーゲル図。

立体射影は多胞体の可視化にも使われる。テンプレート:仮リンクでは、Rn + 1n 次元の多胞体は、n 次元球面に射影され、それが Rn に立体射影される。Rn + 1 から Rn への低次元化は、多胞体の可視化と理解を簡単にしてくれる。

他分野での応用

地図学

ステレオ投影は地球の地図を描く際にも使われる。極付近を描くことが多いが、その他の点をとることもある。

地図学では、球面である地球を平面の地図に投影する図法の一つとして、ステレオ投影のことを平射図法(平射方位図法)またはステレオ図法と言う。角度(形)と面積を保ったまま球を平面に投影する地図は無い、という事実が地図学の基本的問題である。一般的に、統計学は積分をする傾向があるために、正積投影は統計向きである。一方、正角投影航海目的に好ましい。平射図法は正角投影に分類される。また、地球のどちらかの極を中心に投影したとき、経線は原点から放射状にのび、緯線は原点を中心とする円になるという、方位図法としての特性が得られる。すなわち平射図法は正角方位図法である。

地球を赤道半径 a離心率 e回転楕円体とするとき、緯度 φ経度 λ の点をXY平面に投影する式は以下のとおりである[10]。すなわち、座標原点を極点にとり、極点から赤道へ向かう方向を正方向とした中央経線をX軸に設定し、当該中央経線の経度を λ0 とするとき、

x=r(φ)cos(λλ0),y=r(φ)sin(λλ0)
r(φ)=2aexp(q(φ))(1+e)1+e(1e)1e

に投影する。ただし、q(φ) は、緯度 φ に対する等長緯度である。

以前は正角性に注目して使用されることも多かったが、近年ではメルカトル図法横メルカトル図法ガウス・クリューゲル図法)、斜軸メルカトル図法ランベルト正角円錐図法などにほとんど置き換わり、ごく狭い国(オランダなど)や北極南極などに限られるようになった。ただし天気図については、気圧気温の傾斜を正角で観察する必要がある上に、北半球南半球全体で考察することが増えたため、平射図法の使用が増えている。また「全ての小円が必ず円になる」という性質に着目して、クレーターが多い面の全図に使われた例[11] がある。 テンプレート:Seealso

結晶学

テンプレート:Main

[111]方向ダイヤモンド格子結晶の極点図形

結晶学では、3次元空間内の結晶軸や結晶面の方向は、例えばX線電子線回折パターンの解釈で中心となる幾何学的な概念である。 これらの方向は、上記の線と面の可視化の節で述べたように、可視化することができる。 それは、結晶軸と結晶面に対する極は、北半球に交わり、それゆえにステレオ投影を用いてプロットできる。 極点のプロットは、極点図形と呼ばれる。

地学

テンプレート:節スタブ

写真

平射投影を用いた球面パノラマ投影

脚注

  1. 訳注:本文中にもあるように、planisphereは「星図」の古い呼び名なので、星図投影とも読める。
  2. 2.0 2.1 Snyder (1993).
  3. According to (Snyder 1993), although he acknowledges he did not personally see it
  4. Snyder (1989).
  5. According to (Elkins, 1988) who references Eckert, "Die Kartenwissenschaft", Berlin 1921, pp 121--123
  6. Wulff, George, Untersuchungen im Gebiete der optischen Eigenschaften isomorpher Kristalle: Zeits. Krist.,36, l-28 (1902)
  7. Cf. Apostol (1974) p. 17.
  8. Cf. Pedoe (1988).
  9. Cf. Shafarevich (1995).
  10. 河瀬和重 (2014): Lambert正角円錐図法及びその極限としての平射図法の座標換算式に係る包括的導出に関する研究, 平成25年度調査研究年報, 国土地理院技術資料A4-No.12, 80–83
  11. 月の地形図 国土地理院サイト

参考文献

関連項目

外部リンク

テンプレート:Normdaten テンプレート:Sci-stub