トモグラフィー

トモグラフィー（英: tomography）は、物理探査、医療診断等で用いられる逆解析技術の一つ。日本語訳は、断層映像法または断層影像法である。コンピュータを用いて処理することで、画像を構成する技術はコンピュータ断層撮影と呼ばれる。 その多くは、対象領域を取り囲む形で、走査線（線源と検出器）を配置し、内部の物性（音速、比抵抗、音響インピーダンス、密度など）の分布を調べる技術である。評価したい対象物によって、X線CT、地震波トモグラフィー、海洋音響トモグラフィーなどと呼ばれている。 本記事では、トモグラフ像の撮影と、復元について、原理と装置構成を説明する。トモグラフ像の撮影方法には、主に、平行ビーム光学系を用いる方法（図2参照）と、扇形ビーム（ファンビーム）光学系（図3参照）と円錐ビーム(コーンビーム)を用いる方法がある。 CT画像再構成法は解析的再構成法、代数的再構成法、統計的再構成法に大別され、逆投影法は解析的再構成法に分類され、逐次近似画像再構成法は代数的再構成法と統計的再構成法に分類される。 これまでCT画像再構成法の主流はフィルタ補正逆投影法（filtered back projection：FBP法）であったが、近年では画像ノイズ低減効果やアーチファクト低減効果が期待される逐次近似画像再構成法（iterative reconstruction：IR法）が増えつつある。 逆投影法（Back projection） 逆投影法では1回の計算で解(再構成像)が求まる。 逐次近似画像再構成法（Iterative reconstruction） 逐次近似法は最初に初期画像を仮定してこの画像から計算で作成した投影(順投影)と実測投影との整合性を反復計算によって高めていく手法で、反復計算により計算時間を多く必要とするものの、コンピュータの高速化に伴って統計雑音の性質、装置の分解能、被写体の滑らかさなどの事前情報などを式中に組み込める柔軟性や近年発展の著しい圧縮センシング理論を取り入れることにより徐々に増えつつある。 トモグラフィーの数学的な基礎はラドン変換と、ラドン逆変換である。ラドン変換は、トモグラフィーの基本原理であるばかりでなく、例えば、ハフ変換等にも応用される。応用範囲の広い数学的手法であるが、ここでは、トモグラフィーのモデル化という観点に重きをおいて説明する。 関数 μ ( x , y ) {\displaystyle \mu (x,y)} のラドン変換は、以下の式で与えられる。 p ( s , θ ) = − ∫ − ∞ ∞ μ ( s cos ⁡ θ − t sin ⁡ θ , s sin ⁡ θ + t cos ⁡ θ ) d t {\displaystyle p(s,\theta )=-{\int }_{-\infty }^{\infty }\mu (s\cos \theta -t\sin \theta ,s\sin \theta +t\cos \theta )\,dt} 即ち、 「 μ ( x , y ) {\displaystyle \mu (x,y)} のラドン変換の ( θ , s ) {\displaystyle (\theta ,s)} での値 p ( s , θ ) {\displaystyle p(s,\theta )} 」は、 「関数 μ ( x , y ) {\displaystyle \mu (x,y)} の直線 l [ θ , s ] ( t ) {\displaystyle {l}_{[\theta ,s]}(t)} に沿う線積分の値」である。但し、 l [ θ , s ] ( t ) {\displaystyle {l}_{[\theta ,s]}(t)} は、 l [ θ , s ] ( t ) = [ s cos ⁡ θ − t sin ⁡ θ s sin ⁡ θ + t cos ⁡ θ ] {\displaystyle {l}_{[\theta ,s]}(t)={\begin{bmatrix}s\cos \theta -t\sin \theta \\s\sin \theta +t\cos \theta \\\end{bmatrix}}} で定まる直線（tについての直線）である。ここで、上式をtについての直線とみなす際には、 θ {\displaystyle \theta } ,sは、固定されているものと考えるが、その際、 θ {\displaystyle \theta } は前記の直線の傾き角を表し、sは、前記の曲線と原点との間の距離を表すことに注意されたい。 本節では、座標(x,y)における被写体の吸収係数を μ ( x , y ) {\displaystyle \mu (x,y)} とおく。そのうえで、 吸収係数の位置依存性 μ ( x , y ) {\displaystyle \mu (x,y)} にラドン変換を施すことで、測定結果、即ち透過光によって得られた像 p ( s , θ ) {\displaystyle p(s,\theta )} が得られる(モデル化される)こと 測定結果にラドン逆変換を施すことで、 μ ( x , y ) {\displaystyle \mu (x,y)} が復元されること。 を説明する。 被写体を光線が透過した際に、透過光がどれだけ減衰するかを考えることで、上記のラドン変換が導出される。以下、その導出を行う。 ラドン変換を考える際、光線は幾何光学的な光を考える。即ち、光線は、極めて直進性がよく、吸収はされるが、回折や散乱をしないと考え、さらに反射もしないと考えてよいとする。例えばX線を、人体に透過させる場合には、このように考えて差しさわりない。幾何光学において、光線は直線で表される。光線の軌跡が、x-y断面上の直線 l {\displaystyle l} で表される場合について考える。 吸光が、ランベルト・ベールの法則に従うとすると、 前記光線の入射強度を I 0 {\displaystyle {I}_{0}} 、透過後の強度を I {\displaystyle I} 表記したとき、 I = I 0 exp ⁡ ( − ∫ μ ( x , y ) d l ) = I 0 exp ⁡ ( − ∫ − ∞ ∞ μ ( l ( t ) ) | l ˙ ( t ) | d t ) {\displaystyle I=I_{0}\exp \left({-\int \mu (x,y)\,dl}\right)=I_{0}\exp \left({-{\int }_{-\infty }^{\infty }\mu (l(t))\,|{\dot {l}}(t)|dt}\right)} が成り立つ。従って、光線lに沿った吸光度を p l {\displaystyle p_{l}} と表すと、 p l = ln ⁡ ( I / I 0 ) = − ∫ μ ( x , y ) d l = − ∫ − ∞ ∞ μ ( l ( t ) ) | l ˙ ( t ) | d t {\displaystyle p_{l}=\ln(I/I_{0})=-\int \mu (x,y)\,dl=-{\int }_{-\infty }^{\infty }\mu (l(t))\,|{\dot {l}}(t)|dt} 次に、x-y平面に対し、角度θをなす光束を考える。この光束の像について考察しよう。新たに s − t {\displaystyle s-t} 座標系を、 [ s t ] = [ cos ⁡ θ sin ⁡ θ − sin ⁡ θ cos ⁡ θ ] [ x y ] {\displaystyle {\begin{bmatrix}s\\t\end{bmatrix}}={\begin{bmatrix}\cos \theta &\sin \theta \\-\sin \theta &\cos \theta \\\end{bmatrix}}{\begin{bmatrix}x\\y\end{bmatrix}}} により定義する。即ち、s-t座標系は、x-y座標系を角度θだけ回転した座標系である。このとき、回転行列の性質から、 [ x y ] = [ cos ⁡ θ − sin ⁡ θ sin ⁡ θ cos ⁡ θ ] [ s t ] = [ s cos ⁡ θ − t sin ⁡ θ s sin ⁡ θ + t cos ⁡ θ ] {\displaystyle {\begin{bmatrix}x\\y\end{bmatrix}}={\begin{bmatrix}\cos \theta &-\sin \theta \\\sin \theta &\cos \theta \\\end{bmatrix}}{\begin{bmatrix}s\\t\end{bmatrix}}={\begin{bmatrix}s\cos \theta -t\sin \theta \\s\sin \theta +t\cos \theta \\\end{bmatrix}}} となる。今、上式において、sとθを固定すると、上式は、tを変数とする直線と見做せる。 [ x y ] = t [ − sin ⁡ θ cos ⁡ θ ] + [ s cos ⁡ θ s sin ⁡ θ ] {\displaystyle {\begin{bmatrix}x\\y\end{bmatrix}}=t{\begin{bmatrix}-\sin \theta \\\cos \theta \\\end{bmatrix}}+{\begin{bmatrix}s\cos \theta \\s\sin \theta \\\end{bmatrix}}} のように書くと、より直線らしく見えるであろう。 即ち、x-y平面に対し、角度θをなす光束は、以下の l [ θ , s ] ( t ) {\displaystyle {l}_{[\theta ,s]}(t)} で定まる直線を、すべてのsにわたって集めてきたものと考えられる。 l [ θ , s ] ( t ) = t [ − sin ⁡ θ cos ⁡ θ ] + [ s cos ⁡ θ s sin ⁡ θ ] {\displaystyle {l}_{[\theta ,s]}(t)=t{\begin{bmatrix}-\sin \theta \\\cos \theta \\\end{bmatrix}}+{\begin{bmatrix}s\cos \theta \\s\sin \theta \\\end{bmatrix}}} さらに、光線 l [ θ , s ] ( t ) {\displaystyle {l}_{[\theta ,s]}(t)} による吸光量を、 p ( s , θ ) {\displaystyle p(s,\theta )} と書くと、 p ( s , θ ) = − ∫ − ∞ ∞ μ ( s cos ⁡ θ − t sin ⁡ θ , s sin ⁡ θ + t cos ⁡ θ ) d t {\displaystyle p(s,\theta )=-{\int }_{-\infty }^{\infty }\mu (s\cos \theta -t\sin \theta ,s\sin \theta +t\cos \theta )\,dt} が成り立つ。 以上から、x軸に対し、角度θをなす平行光束による透過像（一次元透過像）のプロファイルが、ラドン変換によって与えられることが判った。この、 p ( s , θ ) {\displaystyle p(s,\theta )} が、CT撮影により測定される測定データである。 ラドン逆変換とは、実測された p ( s , θ ) {\displaystyle p(s,\theta )} から、 μ ( x , y ) {\displaystyle \mu (x,y)} を復元する作業のことを指すが、これを説明するためには2変数関数のフーリエ変換について知っておく必要があるので、簡単に復習する。 まず、 μ ( x , y ) {\displaystyle \mu (x,y)} のフーリエ変換とは、 μ