《積分幾何学》

【せきぶんきかがく (integral geometry)】

“ビュッフォンの針”のように図形に関係した確率を幾何確率といい, これらの理論的な部分は, 積分幾何学(integral geometry)を基礎にしている. 積分幾何学は幾何確率に応用可能であるばかりでなく, もっと広く図形に関する様々な局面で基礎となり得ると考えられる. そこで積分幾何学の基礎概念と直線に関する主公式にかぎって述べることにする. 詳しくは[1], [2]を参照されたい. 積分幾何学の基礎概念とは合同変換によって不変な測度を求めることにあるが, 例を直観では分かりづらい直線の集合に採り, 以下で議論を展開する.

図1：直線に下ろした垂線 スタイル検討

　図1のように, 直線 $g\,$ があって, 原点からこれに下ろした垂線の長さを $p\,$ , 垂線と 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x\, } 軸との角度を $\theta \,$ とすると, この直線は

$x\cos \theta +y\sin \theta =p\,$　　　　　$(1)\,$

と表される. そこで, この p と \theta で, 直線の集合(たとえばある領域と交わる直線の集合など)の測度を

$m(X)=\int _{X}f(p,\theta )\;{\mbox{d}}p\;{\mbox{d}}\theta \,$　　　　　$(2)\,$

とおき, 以下, 合同変換で不変の条件を考えよう. 直線 (1)で表されるような, 直線の集合 $X\,$ が, 合同変換

${\begin{array}{lll}x'&=&x\cos \alpha -\sin \alpha +a\\y'&=&x\sin \alpha +y\cos \alpha +b\end{array}}\,$　　　　　$(3)\,$

によって,

$x'\cos \theta '+y'\sin \theta '=p'\,$

なる直線の集合 X' に変換されたものとする. すると, 上式より

${\begin{array}{lll}\theta &=&\theta '-\alpha ,\\p&=&p'-a\cos \theta '-b\sin \theta '\end{array}}\,$　　　　　$(4)\,$

が得られる. ここで, 一様な直線の分布を考えているので, $X\,$ と $X'\,$ の測度は等しく, すなわ ち

$\int _{X}f(p,\theta )\;{\mbox{d}}p\;{\mbox{d}}\theta =\int _{X}f(p',\theta ')\;{\mbox{d}}p'\;{\mbox{d}}\theta '\,$

でなければならない. ところが, 式 (4)から $(p',\theta ')\,$ より $(p,\theta )\,$ への変数変換のヤコビアンは1なので, 上式の右辺を変形して

$\int _{X}f(p,\theta )\;{\mbox{d}}p\;{\mbox{d}}\theta =\int _{X}f(p',\theta ')\;{\mbox{d}}p\;{\mbox{d}}\theta$

と表すことができ, これは $X\,$ をどのようにとっても成立しなければならないので,

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f(p,\theta)=c\, } 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ( \, } 定数構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ) \, }

が得られる. そこで $c=1\,$ とし, ${\mbox{d}}p\;{\mbox{d}}\theta \,$ を ${\mbox{d}}G\,$ とおいて, これを直線の集合 $G\,$ の密度を表わすとすれば, 直線の集合の測度が式 (2)より

$m(X)=\int _{X}{\mbox{d}}G\,$ 　　　　　$(5)\,$

と導かれる. ここで, 図2のような滑らかな曲線ABがあり, その長さを 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L \, } とする. AB上の点Pでこの曲線と交わる直線を 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g\, } , 曲線ABのPにおける接線と $g\,$ とのなす角度を 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \varphi \, } , AからPまでの曲線上の長さを 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s \, } とおく. すると, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (p,\theta) \, } から 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (s,\varphi) \, } への変数変換におけるヤコビアンは $|\sin \varphi |\,$ となり, 曲線ABと交わる直線のすべては, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 0\leq s\leq L, 0\leq\varphi \leq \pi \, } の範囲を考えればよい. ゆえに

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{0}^{L}\! \int_{0}^{\pi}\sin\varphi \; \mbox{d}s \;\mbox{d}\varphi=2L\, }

図２：曲線と直線 スタイル検討

となるので, この 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 2L \, } を曲線 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle AB \, } と交わる直線の集合の測度 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \textstyle \int\mbox{d}p \; \mbox{d}\theta \, } としてよいかにみえる. ところが上記の積分は曲線 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle AB\, } に沿って行なわれたので, この曲線と二つ以上の点で交わる直線については, それぞれの交点で計算されている. したがって, 上記の 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 2L \, } という値は求めたい測度よりも大きい数値であり, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n(p,\theta) \, } でこの重複を表現すれば

$\int {n}(p,\theta )\;{\mbox{d}}p\;{\mbox{d}}\theta =2L\,$

が得られる. ところで, 曲線ABを滑らかなものとして議論してきたが, 滑らかな曲線を有限個つなげたものや, 不連続なものについても, ほとんどいたるところで接線は存在するので, 上式は成立する. そして対象となる曲線と交わらない直線については, $n(p,\theta )=0\,$ と解釈し,

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int{n}\mbox{d}G=2L\, } 　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (6)\, }

が導かれる. さて, これをもとに, 簡単な図形と交わる直線の集合の測度を求めよう. まず上記の曲線ABが直線分であれば, ほとんどいたるところ 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n=1 \, } なので, 長さ $L\,$ の線分と交わる直線の集合 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle X \, } の測度は

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m(X)=\int\mbox{d}G=2L\, } 　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (7)\, }

図３：構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C\, } と交わる確率 スタイル検討

となる. 次に曲線が凸閉曲線の場合, ほとんどいたるところ 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n=2 \, } が成立するので, 長さ 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L \, } の凸閉曲線と交わる直線の集合 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle X \, } の測度は

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle m(X)=\int\mbox{d}G=L\, } 　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (8)\, }

となる. これを用いると, 以下のような確率を簡単に求めることができる. 図3のように凸閉曲線 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C_0,C \, } があって, 長さをそれぞれ, $L_{0},L\,$ とする. 曲線 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C_0,C \, } に交わる直線の集合を 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle X_0,X \, } とすれば, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C_0 \, } に交わる直線が, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C\, } にも交わる確率 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P \, } は, 式(8)より

$P=m(X)/m(X_{0})=L/L_{0}\,$　　　　　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (9)\, }

となる.

---

参考文献

[1] L.A.Santaló, Integral Geometry and Geometric Probability, Addison-Wesley, Reading, Massachusetts, 1976.

[2] 腰塚武志, 「積分幾何学について(1)～(4)」, 『オペレーションズ・リサーチ』, 21 (1976), 524-529, 591-596, 654-659, 711-717.