トモグラフィー

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動先: 案内検索

トモグラフィー: tomography[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [注釈 1] [注釈 2] [注釈 3] [注釈 4] は、物理探査医療診断等で用いられる逆解析技術の一つ。日本語訳は、断層映像法または断層影像法である。

その多くは、対象領域を取り囲む形で、走査線(線源と検出器)を配置し、内部の物性音速比抵抗、音響インピーダンス密度など)の分布を調べる技術である。評価したい対象物によって、X線CT地震波トモグラフィー海洋音響トモグラフィーなどと呼ばれている。

Figure1: CT撮影(人体)の様子[注釈 1]

概要[編集]

本記事では、トモグラフ像の撮影と、復元について、原理と装置構成を説明する。 トモグラフ像の撮影方法には、主に、平行ビーム光学系を用いる方法(図2参照)と、扇形ビーム(ファンビーム)光学系(図3参照)を用いる方法がある [注釈 5] [注釈 6]

平行ビーム光学系を用いたトモグラフ像の撮影と復元[編集]

被写体のx-y断面を考える。トモグラフィーの数学的な基礎はラドン変換と、ラドン逆変換である。ラドン変換は、トモグラフィーの基本原理であるばかりでなく、 例えば、ハフ変換等にも応用される [16][17] [18] [19] 応用範囲の広い数学的手法であるが、ここでは、トモグラフィーのモデル化という観点に重きをおいて説明する。

Figure2: 平行ビーム照射光学系によるトモグラフ像撮影原理;
被写体と、透過光との角度を、とするような、平行ビーム照射光学系を考える。この光学系による投影像は、ある種の線積分の結果と見做すことが出来る。前記の線積分は、被写体をビームが貫通した際に生じた減衰量(attenuation)を表している。図中の符号はそれぞれ、以下の通り (1)被写体, (2)平行ビーム光源, (3)スクリーン, (4)透過光, (5)平行ビーム光源、スクリーンの軌道, (6)平行ビーム光源、スクリーンの軌道の中心, (7)透影像(一次元画像;をそれぞれあらわす。

関数のラドン変換は、以下の式で与えられる。

即ち、 「のラドン変換の での値 」は、 「関数の直線に沿う線積分の値」である。 但し、は、

で定まる直線(tについての直線)である。ここで、上式をtについての直線とみなす際には、 ,sは、固定されているものと考えるが、 その際、は前記の直線の傾き角を表し、sは、前記の曲線と原点との間の距離を表すことに注意されたい。


本節では、被写体の、座標(x,y)における吸収係数をとモデル化する。 そのうえで、

  1. 測定結果、即ち透過光によって得られた像のシリーズが、にラドン変換を施すことで得られた関数として表現される(モデル化される)こと
  2. 測定結果にラドン逆変換を施すことで、が復元されること。

を説明する。


平行ビーム光学系によるトモグラフ像撮影のモデル化[編集]

被写体を光線が透過した際に、透過光がどれだけ減衰するかを考えることで、上記のラドン変換が導出される。以下、その導出を行う。

ラドン変換を考える際、光線は幾何光学的な光を考える。即ち、光線は、極めて直進性がよく、吸収はされるが、回折散乱をしないと考え、さらに反射もしないと考えてよいとする。例えばX線を、人体に透過させる場合には、このように考えて差しさわりない。 幾何光学において、光線は直線で表される。光線の軌跡が、x-y断面上の直線で表される場合について考える。

吸光が、ランベルト・ベールの法則に従うとすると、 前記光線の入射強度を、透過後の強度を表記したとき、

が成り立つ。従って、光線lに沿った吸光度と表すと、

次に、x-y平面に対し、角度θをなす光束を考える。この光束の像について考察しよう。 新たに座標系を、

により定義する。即ち、s-t座標系は、x-y座標系を角度θだけ回転した座標系である。このとき、回転行列の性質から、

となる。今、上式において、sとθを固定すると、上式は、tを変数とする直線と見做せる。

のように書くと、より直線らしく見えるであろう。

即ち、x-y平面に対し、角度θをなす光束は、以下ので定まる直線を、すべてのsにわたって集めてきたものと考えられる。

さらに、光線による吸光量を、と書くと、

が成り立つ。

以上から、x軸に対し、角度θをなす平行光束による透過像(一次元透過像)のプロファイルが、ラドン変換によって与えられることが判った。この、が、CT撮影により測定される測定データである。


平行ビーム光学系によりトモグラフ像の、復元[編集]

二次元フーリエ変換についての補足[編集]

ラドン逆変換とは、実測されたから、を復元する作業のことを指すが、これを説明するためには2変数関数のフーリエ変換について知っておく必要があるので、簡単に復習する。

まず、のフーリエ変換とは、

である。のことを、 と書く場合もある。 ここで、""は、関数と関数の積(たんなる掛け算)を意味する。

先のと、に対し、以下の等式が成立する。これを、フーリエ逆変換と呼ぶ。

即ち、関数をフーリエ変換した後、フーリエ逆変換すれば、元の関数に戻る。 ここで、""は、関数と関数の積(たんなる掛け算)を意味する。

ラドン逆変換[編集]

実測されたから、今度はを復元することを考える。この復元操作は、数学的に、 以下の2ステップで行われる。

  • ラドン逆変換のステップ(1)

を計算する。

  • ラドン逆変換のステップ(2)

を計算する。

上記2ステップの計算を実施することを、ラドン逆変換という。以下、上記ステップにて、μ(x,y)が示せることを示す。

証明第一段階[編集]

(1) の証明

を、変数sについてフーリエ変換(一変数関数としてフーリエ変換)したものを、とする。即ち、

とする。これは、上記ステップ(1)の変換に他ならない。 当たり前のことだが、このは、とは別物である。そもそも定義が異なる。

今、の定義式、即ち、

を、上式に代入すると、




である。

今、2変数ベクトル値関数と、を、それぞれ、

と定めると、明らかに、

である。さらに、

である。従って、


である。

さらに、上式を、

に注意して積分の変数変換を施すと、




一方で、

の、(u,v) に を代入すると、


従って、

が判る。

証明第二段階[編集]

(2)の証明

第一段階の結論、すなわち、

より、以下の等式が、任意の(x,y)に対して成り立つ。

上式の右辺に、以下の変数変換:

を施すと、 積分の変数変換の公式から、


一方、二次元のフーリエ逆変換を考えると、

であるため、

を得る。

扇形ビーム光学系によるトモグラフ像の測定と復元[編集]

特に医療機器の場合には、図3のような扇形ビームを用いる場合が多い。
Figure3:扇型ビーム照射光学系によるトモグラフ像撮影原理。 ;

コンピュータ断層撮影[編集]

X線放射線以外にも電磁波などを利用して調査対象を走査コンピュータを用いて処理することで、調査対象の内部構造の画像を構成する技術。物質固有の透過率、誘電率、伝搬性などの差異を基に数値化して画像を取得する。

地震波トモグラフィー[編集]

地震波の伝播時間を用いて、地球内部の3次元速度構造を求める手法。

海洋音響トモグラフィー[編集]

海洋中の音波の伝播時間を測定して、音波の伝わった海洋の内部の状態を調べる手法。

ミュー粒子トモグラフィー[編集]

1kmの岩盤でも透過する宇宙から降り注ぐミュー粒子を使用して内部構造を調べる手法。火山マグマの上昇やピラミッドの内部構造の調査[20]原子炉の炉心の調査に使用される[21][22][23]

ジオトモグラフィ[編集]

ジオトモグラフィ(地盤断層撮影)は地質構造を調べる目的で使用される。弾性波トモグラフィと比抵抗トモグラフィが含まれる[24]

弾性波トモグラフィ[編集]

弾性波トモグラフィは弾性波を使用して地質構造を可視化する[24]

音響トモグラフィ[編集]

音響トモグラフィは高周波数の音響波を用いて、地盤や構造物の内部構造、岩盤やコンクリート内部の亀裂、樹木内部の空洞を非破壊で可視化する[25]

比抵抗トモグラフィ[編集]

比抵抗トモグラフィは地中にの電気伝導度を元に地下構造を可視化する[24]

光コヒーレンストモグラフィ[編集]

細胞レベルの高速・高分解能で深さ5mm未満の断層画像化に使用される[26]

電磁誘導トモグラフィ[編集]

電磁誘導トモグラフィは電磁波の伝播経路の伝播特性から内部構造を可視化する[27]

電気インピーダンストモグラフィ[編集]

電気インピーダンストモグラフィは電気インピーダンスを元に内部構造を可視化する[28][29]

電気容量トモグラフィ[編集]

電気容量トモグラフィは電気容量を元に内部構造を可視化する[30][31]

磁気誘導トモグラフィ[編集]

分析対象内に送信コイルによって誘導電流を生じさせてそれによって発生した磁場を検出することで対象の内部構造を可視化する[32]

熱音響トモグラフィ[編集]

熱音響トモグラフィは加熱によって生じた音響を検出して内部構造を可視化する[33]

中性子トモグラフィ[編集]

中性子トモグラフィとは中性子を利用して内部構造を可視化する[34]

テラヘルツ波トモグラフィ[編集]

テラヘルツ波(周波数100GHz~ 10THz、波長30um~3mm)を利用して内部構造を可視化する[35]

量子トモグラフィ[編集]

量子トモグラフィとは未知の量子的対象を完全に同定する手法の総称であり、量子情報処理実験において重要な役割を担っている[36]

電離層トモグラフィ[編集]

電離層トモグラフィとは電離層内での電波の伝播を利用して構造を可視化する[37]

水蒸気トモグラフィ[編集]

多数のGPS観測点でグローバル・ポジショニング・システム(GPS)からの電波を受信する事で大気圏を伝播する電波の水蒸気による衛星視線方向の大気遅延量を利用して大気中の水蒸気の空間分布を算出する[38][39]

脚注・参考文献[編集]

脚注[編集]

  1. ^ a b 図1のような医療用CT撮影機器は、医療機器に該当するため薬事法の規制をうける。従って、それぞれの製品毎に必ず添付文書が必ず存在する(薬事法上の一般名称は、全身用X線CT診断装置)。添付文書は 医薬品医療機器総合機構のデータベース[1]から検索できる。本記事では、東芝メディカルシステムズ Asteion TSX-021Bの添付文書[2]と、日立メディコ製の全身用X線CT診断装置 SCENARIA[3] の添付文書[4]を参考にした。
  2. ^ いくつかの病院のホームページ上で、CTスキャン装置の写真を見ることができる。 一例→[5](麦島内科の例)
  3. ^ 図面の豊富なその他の特許として、例えば、次のようなものがある。 アメリカ合衆国特許第6,879,657号アメリカ合衆国特許第6,574,296号アメリカ合衆国特許第6,775,346号アメリカ合衆国特許第7,215,734号
  4. ^ A. M. Cormackの2論文は、AIPの重要論文とされている。[6]
  5. ^ プローブとするビームは、(主にX線)、磁場に加え、 電子線(平均自由行程が短いため、真空中に限る)変わり種としては、ミュオン(山の断層写真の撮影例がある)等がある。
  6. ^ 少なくとも医療機器の場合には、扇形ビームをベースとした方法がほとんどである。さらにメーカー独自の改良をが加わっている。最近の技術動向は、特許庁特許事務所等が作成したパテントマップ等からある程度解読可能である。 特許庁作成版としては平成15年版[7] あるいは平成23年版[8]の パテントマップが公開されている。本文では15年版を引用している(23年版は 過去の動向や歴史、基盤技術に関する分析が簡潔するため)本シリーズの前シリーズには、 「技術分野別特許マップ」[9] があり、基本特許の特定などの有益な情報が充実していたが、「技術分野別特許マップ」シリーズに比べ、15年版ですら そういった基礎的な情報は大幅にプアになっていて、産業統計に近い状況になっている。 尚、「技術分野別特許マップ」シリーズでは、CT関連を直接扱ったものはない。

参考文献[編集]

  1. ^ Avinash Kak; Malcolm Slaney (1988). Principles of Computerized Tomographic Imaging. IEEE Press. ISBN 9780879421984. http://www.slaney.org/pct/pct-toc.html. 
  2. ^ 篠原 広行、坂口 和也、橋本 雄幸「Excelによる画像再構成入門 (画像再構成シリーズ)」、医療科学社、1993年2月ISBN 9784860033736
  3. ^ 斎藤 恒雄「画像処理アルゴリズム (アルゴリズム・シリーズ)」、近代科学社、1993年2月ISBN 9784764902053
  4. ^ フーリエ解析(13): フーリエ変換の医療分野への応用例 [10]
  5. ^ 梅垣 寿春「情報数理の基礎―関数解析的展開 (Information & Computing)」、サイエンス社1993年7月ISBN 9784781907079
  6. ^ 河田 聡、南 茂夫「「科学計測のための画像データ処理―パソコン EWS活用による画像計測&処理技術 」」、CQ出版1994年4月ISBN 9784789830300
  7. ^ Johann Radon, Uber die Bestimmung von Funktionen durch ihre Integralwerte l?ngs gewisser Mannigfaltigkeiten, Computed tomography (Cincinnati, Ohio, 1982) Proc. Sympos. Appl. Math., vol. 27, Amer. Math. Soc., Providence, R.I., 1982, pp. 71?86 (German). MR 692055 (84f:01040)
  8. ^ A. M. Cormack;"Representation of a function by its line integrals, with some radiological applications" J. Appl. Phys. 34, 2722-2727 (1963) [11]
  9. ^ Hounsfield GN; Computerised transverse axial scanning (tomography) I. Description of system. Br J Radiol 46: 1016-1022, 1973.[12]
  10. ^ Hounsfield GNのノーベル賞講演 [13]
  11. ^ Godfrey Newbold Hounsfield US4,115,698 [14](Godfrey Newbold Hounsfieldによる特許)。 同特許のパテントファミリー等は、Espacenet等より参照可[15]
  12. ^ 藤田保健衛生大学辻岡勝美研究室のページに、有益なリソースが多数ある。[16]
  13. ^ Steven W. Smith;"The Scientist and Engineer's Guide to Digital Signal Processing" [17]
  14. ^ 日立メディコ CTシステム関連技術 [18]
  15. ^ 特許庁総務部技術調査課  「医用画像診断装置に関する特許出願動向報告」 平成15年5月8日 [19]
  16. ^ 国際電気通信基礎技術研究所による特開平05-012438[20]
  17. ^ (新エネルギー・産業技術総合開発機構)即効型地域新生コンソーシアム研究開発 柔軟変形物ハンドリング用ビジョンチップの研究開発報告書 [21]および、立命館大学講義ノート [22]
  18. ^ MATLAB解説記事より[23]
  19. ^ [24]
  20. ^ 東大地震研:浅間山の 内部構造再現 素粒子使い立体的に 毎日新聞(2010年3月9日)
  21. ^ Muon scans confirm complete reactor meltdown at Fukushima
  22. ^ Muon scans confirm complete reactor meltdown at Fukushima Reactor #1
  23. ^ Our Next Two Steps for Fukushima Daiichi Muon Tomography
  24. ^ a b c 弾性波トモグラフィ
  25. ^ 地中を”見える化”する音響トモグラフィ地盤調査
  26. ^ 医療におけるフォトニクス技術の展開
  27. ^ 電磁誘導トモグラフィの3次元解析ソフトウェアの開発
  28. ^ Electrical Impedance Tomography:EITの紹介
  29. ^ 磁気を併用した生体電気インピーダンス CT の開発研究
  30. ^ Electrical Impedance Tomography:EITの紹介
  31. ^ 磁気を併用した生体電気インピーダンス CT の開発研究
  32. ^ 磁気誘導トモグラフィ
  33. ^ 熱音響トモグラフィ法および熱音響トモグラフ
  34. ^ 中性子断層撮影法の基礎
  35. ^ 実時間2次元テラヘルツ断層イメージング
  36. ^ 量子情報関東Student Chapter
  37. ^ 電離層電子密度トモグラフィ
  38. ^ GPS meteorology
  39. ^ GPS MET COM -GPSによる水蒸気観測-

関連項目[編集]