(57)【要約】      （修正有）【課題】開口直径の大きなレンズなどの集向装置を使用し、正確なCT再構成が可能となるようにする。【解決手段】集向装置を用いて取得したプロジェクションデータ群３８から、MLEM法を用いて観察領域の輝度あるいは対数吸収量の分布を再構成する方法であって、各プロジェクション素４３に配置されたプロジェクションデータを、当該プロジェクション素４３からレンズ４０の開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビーム４８に沿って再構成領域４１に配置するように、各プロジェクション素４３に配置されたプロジェクションデータと、各プロジェクション素４３の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビーム４８が通過する再構成領域４１の全ての正方ピクセルに、それぞれの正方ピクセルのCT値とビーム通過面積との積に応じて分配する。【選択図】図７
2011102758.tif 000002
【特許請求の範囲】
【請求項１】
観察領域の輝度分布をレンズなどの集向装置を用いて取得したプロジェクションデータ群から、観察領域の輝度分布を二次元再構成する方法であって、
前記プロジェクションデータ群（積算輝度値）を、プロジェクションデータエリアの各プロジェクション素に配置し、前記プロジェクション素に配置されたプロジェクションデータを、そのバックプロジェクションビームが通過する再構成領域の全ての正方ピクセルに、それぞれの正方ピクセルのCT値（再構成輝度値）とビーム通過面積との積に応じて分配するMLEM法を用いて前記観察領域の輝度分布を再構成する方法であり、
（ａ）各プロジェクション素に配置されたプロジェクションデータを、当該プロジェクション素からレンズなどの集向装置の開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビームに沿って前記再構成領域に配置するように、各プロジェクション素に配置されたプロジェクションデータと、各プロジェクション素の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビームが通過する前記再構成領域の全ての正方ピクセルに、それぞれの正方ピクセルのCT値（再構成輝度値）とビーム通過面積との積に応じて分配するステップを有して、前記観察領域の輝度分布を再構成する方法。
【請求項２】
（ｂ）各プロジェクション素の開口分布を、レンズ素ごとの通過プロジェクションビームに関して、通過ピクセルの再構成輝度と通過面積との積の合計値を個別合計値として算出し、また全てのレンズ素に関して総合計を行って総合計値を算出し、総合計値に対する個別合計値の割合として、算出するステップを有し、前記（ａ）、（ｂ）のステップを繰り返すことにより、前記観察領域の輝度分布を再構成することを特徴とする請求項１に記載の方法。
【請求項３】
観察領域の輝度分布をレンズなどの集向装置を用いて取得したプロジェクションデータ群から、観察領域の輝度分布を三次元再構成する方法であって、
前記プロジェクションデータ群（積算輝度値）を、プロジェクションデータエリアの各プロジェクション素に配置し、前記プロジェクション素に配置されたプロジェクションデータを、そのバックプロジェクションビームが通過する再構成領域の全ての立方ボクセルに、それぞれの立方ボクセルのCT値（再構成輝度値）とビーム通過体積との積に応じて分配するMLEM法を用いて前記観察領域の輝度分布を再構成する方法であり、
（ａ）各プロジェクション素に配置されたプロジェクションデータを、当該プロジェクション素からレンズなどの集向装置の開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビームに沿って前記再構成領域に配置するように、各プロジェクション素に配置されたプロジェクションデータと、各プロジェクション素の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビームが通過する前記再構成領域の全ての立方ボクセルに、それぞれの立方ボクセルのCT値（再構成輝度値）とビーム通過体積との積に応じて分配するステップを有して、前記観察領域の輝度分布を再構成する方法。
【請求項４】
（ｂ）各プロジェクション素の開口分布を、レンズ素ごとの通過プロジェクションビームに関して、通過ピクセルの再構成輝度と通過体積との積の合計値を個別合計値として算出し、また全てのレンズ素に関して総合計を行って総合計値を算出し、総合計値に対する個別合計値の割合として、算出するステップを有し、前記（ａ）、（ｂ）のステップを繰り返すことにより、前記観察領域の輝度分布を再構成することを特徴とする請求項３に記載の方法。
【請求項５】
観察領域の吸収量の分布をレンズなどの集向装置と背景拡散放射装置を用いて取得したプロジェクションデータ群（積算吸収量の対数）から、観察領域の吸収量の分布を二次元再構成する方法であって、
前記プロジェクションデータ群（積算吸収量の対数）を、プロジェクションデータエリアの各プロジェクション素に配置し、前記プロジェクション素に配置されたプロジェクションデータを、そのバックプロジェクションビームが通過する再構成領域の全ての正方ピクセルに、それぞれの正方ピクセルのCT値（再構成対数吸収量）とビーム通過面積との積に応じて分配するMLEM法を用いて前記観察領域の対数吸収量の分布を再構成する方法であり、
（ａ）各プロジェクション素に配置されたプロジェクションデータを、当該プロジェクション素からレンズなどの集向装置の開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビームに沿って前記再構成領域に配置するように、各プロジェクション素に配置されたプロジェクションデータと、各プロジェクション素の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビームが通過する前記再構成領域の全ての正方ピクセルに、それぞれの正方ピクセルのCT値（再構成対数吸収量）とビーム通過面積との積に応じて分配するステップを有して、前記観察領域の対数吸収量の分布を再構成する方法。
【請求項６】
（ｂ）各プロジェクション素の開口分布を、レンズ素ごとの通過プロジェクションビームに関して、通過ピクセルの再構成対数吸収量と通過面積との積の合計値を個別合計値として算出し、また全てのレンズ素に関して総合計を行って総合計値を算出し、総合計値に対する個別合計値の割合として、算出するステップを有し、前記（ａ）、（ｂ）のステップを繰り返すことにより、前記観察領域の対数吸収量の分布を再構成することを特徴とする請求項５に記載の方法。
【請求項７】
観察領域の対数吸収量の分布をレンズなどの集向装置と背景拡散放射装置を用いて取得したプロジェクションデータ群（積算吸収量の対数）から、観察領域の吸収量の分布を三次元再構成する方法であって、
前記プロジェクションデータ群（積算吸収量の対数）を、プロジェクションデータエリアの各プロジェクション素に配置し、前記プロジェクション素に配置されたプロジェクションデータを、そのバックプロジェクションビームが通過する再構成領域の全ての立方ボクセルに、それぞれの立方ボクセルのCT値とビーム通過体積との積に応じて分配するMLEM法を用いて前記観察領域の対数吸収量の分布を再構成する方法であり、
（ａ）各プロジェクション素に配置されたプロジェクションデータを、当該プロジェクション素からレンズなどの集向装置の開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビームに沿って前記再構成領域に配置するように、各プロジェクション素に配置されたプロジェクションデータと、各プロジェクション素の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビームが通過する前記再構成領域の全ての立方ボクセルに、それぞれの立方ボクセルのCT値（再構成対数吸収量）とビーム通過体積との積に応じて分配するステップを有して、前記観察領域の対数吸収量の分布を再構成する方法。
【請求項８】
（ｂ）各プロジェクション素の開口分布を、レンズ素ごとの通過プロジェクションビームに関して、通過ピクセルの再構成対数吸収量と通過体積との積の合計値を個別合計値として算出し、また全てのレンズ素に関して総合計を行って総合計値を算出し、総合計値に対する個別合計値の割合として、算出するステップを有し、前記（ａ）、（ｂ）のステップを繰り返すことにより、前記観察領域の対数吸収量の分布を再構成することを特徴とする請求項７に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、二次元あるいは三次元CT（コンピュータ断層）再構成のように観察領域の物理量分布を再構成する方法に関するものである。
【背景技術】
【０００２】
特許文献１にあるように、Ｘ線を用いる透過型のCT（コンピュータ断層）装置においては、小さなＸ線源から被検体（人体）に向けて線ビーム、シートビームあるいはコーンビームを放射し、被検体（人体）を減衰透過したビームの強度を直接的にＸ線センサーで捕らえることで、プロジェクション（投射）データを得る．このとき、レンズ等の集光装置は用いない。
【０００３】
また、非特許文献１では、主として、観察領域からの放射を再構成する放射型CTについて記述されている。すなわち、非特許文献１では、観察領域から放射される可視光の三次元CT再構成を行っており、計測手法としては、40個のレンズを有する特殊カメラにより、乱流火炎を多方向から同時撮影を行い、得られた40枚の二次元輝度分布から、観察領域の三次元輝度分布の再構成をMLEM（DBP）アルゴリズムにより実行している。
【０００４】
非特許文献１にあるようにCT再構成アルゴリズムには、フィルタードバックプロジェクション法、代数的手法、分配バックプロジェクション（DBP）法（MLEM法）などがある。
【０００５】
以下に、二次元の場合の、プロジェクションの取得手順と、MLEM法によるプロジェクションデータからの再構成手順を記述する。
【０００６】
図１に、プロジェクション取得手順の説明図を示す。図において、１は観察領域、２はレンズ、３は一次元光センサー素子列、４はカメラ、５は観察領域中の任意小領域、６はデータ取得・ＣＴ演算コンピュータである。ここでは、レンズはピンホールレンズとして扱われる。また、プロジェクション取得は、最低でも180゜に渡り、なるべく多方向である方が望ましいが、ここでは、図１に示すように４方向であるとして説明を行う。
【０００７】
放射型CTにおけるプロジェクションデータの取得は、通常の「写真撮影」と等しい。すなわち、観察領域１にのみ輝度が分布しているとし、観察領域中のある小領域５から全方向へ放射される光の一部は、ピンホールレンズ２により光センサー素子列３に結像する。観察領域１に存在するすべての輝度からの結像により、観察領域のプロジェクションデータが光センサー素子列３によって獲得され、カメラ台数分のプロジェクションデータ群がデータ取得・ＣＴ演算コンピュータ６により記録される。
【０００８】
つぎに、MLEM法によるプロジェクションデータからの輝度分布の再構成手順の説明を行う。図２に再構成手順の説明図を示す。この再構成は、データ取得・ＣＴ演算コンピュータ６の内部で仮想的に演算処理される。
【０００９】
図２において、６はデータ取得・ＣＴ演算コンピュータ、７はプロジェクションデータ群、８はプロジェクションデータエリア、９はレンズ中心、１０は再構成領域、１１は再構成領域を分割する細かな正方ピクセルの例、１２はプロジェクション素の例、１３はプロジェクションデータの分布のイメージの例、１４は任意のm番目のプロジェクション素、１５はプロジェクション素１４に配置されたプロジェクションデータ（値P(m)）、１６はプロジェクション素１４からレンズ中心９を通過して伸びるバックプロジェクションビームである。
【００１０】
データ取得・ＣＴ演算コンピュータ６には、上記プロジェクションデータ群７が記録されており、はじめに、仮想空間に、プロジェクション取得時（図１参照）と同じ配置となるように、プロジェクションデータエリア８を一次元光センサー列３の位置に配置し、レンズ中心９をピンホールレンズ２の位置に配置し、再構成領域１０を観察領域１の領域に配置する。再構成領域１０は通常、細かな正方ピクセル１１に分割しておく。正方ピクセルの数をＳとし、各正方ピクセルには、１からＳまでの番地を付けておく。また、プロジェクションデータエリア８も、光センサー素子列３の各光センサー領域に対応するように、プロジェクション素１２に分割しておく。プロジェクション素の数はＭとし、各プロジェクション素には、１からＭまでの番地をつけておく。
【００１１】
データ取得・ＣＴ演算コンピュータ６に記録されているプロジェクションデータ群７を、プロジェクションデータエリア８の各プロジェクション素１２に配置する。ここで、プロジェクションデータの分布のイメージ１３を、プロジェクションデータエリアに沿うように、プロジェクションデータエリアに垂直の距離で図示した。
【００１２】
さて、プロジェクションデータ群７を再構成領域１０の各正方ピクセル（例えば１１）に再配置し、各正方ピクセルの持つＣＴ値（ｓ番目の正方ピクセルのＣＴ値をＣ(s)とする）を再構成する工程に移る。
【００１３】
MLEM法（DBP法）は、繰り返し演算法である。最初の繰り返し演算の前に、初期値として、再構成領域１０中の全ての正方ピクセルに、等しいCT値を配置しておく。
【００１４】
MLEM法（DBP法）では、あるm番目のプロジェクション素１４に配置されたプロジェクションデータ１５（値はP(m)）を、プロジェクション素１４からレンズ中心９を通過して伸びるバックプロジェクションビーム１６に沿って、下記手順により再構成領域１０の各正方ピクセルに配置する。
【００１５】
はじめに、当該のバックプロジェクションビーム１６が通過する全ての通過正方ピクセルを選び出す。その個数がＮ(m)個であるとする。１番目からＮ(m)番目の通過正方ピクセルの番地を、それぞれ、s(m,1)からs(m,N(m))とする。つぎに、バックプロジェクションビーム１６とその通過正方ピクセルが重なる面積A(m,1)〜A(m,N(m))を演算し、その通過正方ピクセルに配置されているCT値C(s(m,1))〜C(s(m,N(m))）を用いて、次の演算を行い、m番目のプロジェクション素１４から伸びるバックプロジェクションビーム１６がｉ番目に通過する正方ピクセルに分配するCT値の分配量 C(m,i)を算出する。
【００１６】
C(m,i)＝P(m)×（A(m,i)×C(s(m,i))）／（ A(m,1)×C(s(m,1))＋・・・＋A(m,N(m))×C(s(m,N(m))）
すなわちMLEM法（DBP法）は、プロジェクションデータP(m)を、そのバックプロジェクションビームが通過する全ての正方ピクセルに、それぞれの正方ピクセルのCT値とビーム通過面積との積に応じて分配する方法である。
【００１７】
全てのプロジェクション素１４について上記の分配量を演算し、全ての正方ピクセルについて分配量を合計することで、全ての正方ピクセルについての新しいCT値が計算される。
【００１８】
以上で、１回の繰り返し演算が終了する。繰り返し演算を適切な回数だけ繰り返すことで、再構成領域１０の各正方ピクセルのCT値の分布が完成し、プロジェクションデータ群からの再構成領域の輝度分布が再構成されたことになる。
【００１９】
なお、三次元の場合も、二次元の場合とほぼ同様に行うことができることが、非特許文献１に記載されている。
【００２０】
また、均一の拡散背景放射が放射されている場合の観察領域の吸収量分布を再構成する場合も、単位長さあたりの吸収量の対数を用いることで、上記の方法が使用できる。すなわち、輝度の代わりに、単位長さ当たりの吸収率の対数を使用して、背景放射経路に沿った単位長さ当たりの吸収量の対数の和として、再構成の演算を行う。
【先行技術文献】
【特許文献】
【００２１】
【特許文献１】特開２００４−１１３８０２号公報
【非特許文献】
【００２２】
【非特許文献１】石野洋二郎，大岩紀生，「乱流火炎形状のスキャンレス三次元CT計測  −40眼カメラによる多方向同時撮影法と三次元CT法による乱流火炎の瞬間輝度分布の三次元計測−」日本燃焼学会誌，第46巻137号，pp.142-152，2004年
【発明の概要】
【発明が解決しようとする課題】
【００２３】
特許文献１にあるように、Ｘ線を用いる透過型のCT（コンピュータ断層）装置においては、小さなＸ線源から被検体（人体）に向けて線ビーム、シートビームあるいはコーンビームを放射し、被検体（人体）を減衰透過したビームの強度を直接的にＸ線センサーで捕らえることで、プロジェクション（投射）データを得る。このとき、レンズ等の集光装置は用いないため、プロジェクションデータ収集時には、被写界深度は考慮しない。
【００２４】
ところが一方で、非特許文献１では、主として、観察領域から放射される可視光の三次元CT再構成のプロジェクションデータの獲得の際に、開口を持つレンズを使用している。この場合には本来、被写界深度を考慮する必要がある。すなわち、レンズ開口直径が大きく、観察領域が大きく、観察領域とレンズまでの距離が近い場合には、焦点面の前後の観察領域のプロジェクション像にはボケが生じ、正確なCT再構成が不可能となる。
【００２５】
本発明は上記点に鑑みて、開口直径の大きなレンズを使用し獲得した、焦点面の前後の観察領域の像にボケが生じているプロジェクションデータを使用して、正確なCT再構成が可能となるようにすることを目的とする。
【課題を解決するための手段】
【００２６】
上記目的を達成するため、請求項１に記載の発明では、観察領域の輝度分布をレンズを用いて取得したプロジェクションデータ群から、観察領域の輝度分布を二次元再構成する方法であって、
前記プロジェクションデータ群を、プロジェクションデータエリアの各プロジェクション素に配置し、前記プロジェクション素に配置されたプロジェクションデータを、そのバックプロジェクションビームが通過する再構成領域の全ての正方ピクセルに、それぞれの正方ピクセルのCT値とビーム通過面積との積に応じて分配するMLEM法を用いて前記観察領域の輝度分布を再構成する方法であり、
（ａ）各プロジェクション素に配置されたプロジェクションデータを、当該プロジェクション素からレンズの開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビームに沿って前記再構成領域に配置するように、各プロジェクション素に配置されたプロジェクションデータと、各プロジェクション素の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビームが通過する前記再構成領域の全ての正方ピクセルに、それぞれの正方ピクセルのCT値とビーム通過面積との積に応じて分配するステップを有して、前記観察領域の輝度分布を再構成する方法を特徴としている。
【００２７】
請求項２に記載の発明では、請求項１に記載の方法において、
（ｂ）各プロジェクション素の開口分布を、レンズ素ごとの通過プロジェクションビームに関して、通過ピクセルの再構成輝度と通過面積との積の合計値を個別合計値として算出し、また全てのレンズ素に関して総合計を行って総合計値を算出し、総合計値に対する個別合計値の割合として、算出するステップを有し、前記（ａ）、（ｂ）のステップを繰り返すことにより、前記観察領域の輝度分布を再構成することを特徴としている。
【００２８】
請求項３に記載の発明では、観察領域の輝度分布をレンズを用いて取得したプロジェクションデータ群から、観察領域の輝度分布を三次元再構成する方法であって、
前記プロジェクションデータ群を、プロジェクションデータエリアの各プロジェクション素に配置し、前記プロジェクション素に配置されたプロジェクションデータを、そのバックプロジェクションビームが通過する再構成領域の全ての立方ボクセルに、それぞれの立方ボクセルのCT値とビーム通過体積との積に応じて分配するMLEM法を用いて前記観察領域の輝度分布を再構成する方法であり、
（ａ）各プロジェクション素に配置されたプロジェクションデータを、当該プロジェクション素からレンズの開口領域を分割したレンズ素を通過し屈折して伸びる多数のバックプロジェクションビームに沿って前記再構成領域に配置するように、各プロジェクション素に配置されたプロジェクションデータと、各プロジェクション素の開口分布とから分乗プロジェクションデータを得、この分乗プロジェクションデータを、そのバックプロジェクションビームが通過する前記再構成領域の全ての立方ボクセルに、それぞれの立方ボクセルのCT値とビーム通過体積との積に応じて分配するステップを有して、前記観察領域の輝度分布を再構成する方法を特徴としている。
【００２９】
請求項４に記載の発明では、請求項３に記載の方法において、
（ｂ）各プロジェクション素の開口分布を、レンズ素ごとの通過プロジェクションビームに関して、通過ピクセルの再構成輝度と通過体積との積の合計値を個別合計値として算出し、また全てのレンズ素に関して総合計を行って総合計値を算出し、総合計値に対する個別合計値の割合として、算出するステップを有し、前記（ａ）、（ｂ）のステップを繰り返すことにより、前記観察領域の輝度分布を再構成することを特徴としている。
【００３０】
また、均一の拡散背景放射が放射されている場合の観察領域の吸収量分布を再構成する場合も、単位長さあたりの吸収量の対数を用いることで、上記の方法が使用できる。すなわち、輝度の代わりに、単位長さ当たりの吸収率の対数を使用して、背景放射経路に沿った単位長さ当たりの吸収量の対数の和として、再構成の演算を行う。
【図面の簡単な説明】
【００３１】
【図１】従来方法が適用できる場合のプロジェクション取得手順の説明図である。
【図２】従来方法（MLEM法，DBP法）の再構成手順の説明図である。
【図３】開口が大きく観察領域に近接したレンズを用いたプロジェクション取得手順の説明図である。
【図４】非合焦による「ボケ」の発生を説明する説明図である。
【図５】非合焦による「ボケ」の発生を詳細に説明する詳細説明図である。
【図６】非合焦時の光センサーが感知する光線の経路を説明する説明図である。
【図７】第１実施形態の再構成方法を説明する説明図である。
【図８】第１実施形態の再構成アルゴリズムのフローチャートである。
【図９】第２実施形態の再構成方法を説明する説明図である。
【図１０】第２実施形態の再構成アルゴリズムのフローチャートである。
【発明を実施するための形態】
【００３２】
（第１実施形態）
以下に、開口が大きく、観察領域に近接したレンズを用いたプロジェクションの取得手順と、本実施形態のＣＴアルゴリズムによるプロジェクションデータからの再構成手順を、いずれも、二次元の場合について記述する。
【００３３】
図３に、開口が大きく、観察領域に近接したレンズを用いたプロジェクション取得手順の説明図を示す。図３において、１は観察領域、２はレンズ、３は一次元光センサー素子列、４はカメラ、５は観察領域中の任意小領域、６はデータ取得・ＣＴ演算コンピュータである。レンズ２の位置および焦点距離は、一次元光センサー素子列３に対して観察領域１の中心が合焦するように設定されている。すなわち、観察領域１の中心からの光はレンズによって光センサー素子列３上の、観察領域中心とレンズ中心とを結ぶ直線上の点に焦点を結ぶ。ここで図３は、レンズ２の中心位置と観察領域１の中心との距離に比べて、レンズ２の直径および観察領域１の大きさが無視できない場合である。また、プロジェクション取得は、最低でも180゜に渡る方が望ましく、なるべく多方向である方が望ましいが、ここでは、図３に示すように４方向であるとして説明を行う。
【００３４】
放射型CTにおけるプロジェクションデータの取得は、通常の「写真撮影」と等しい。すなわち、観察領域１にのみ輝度が分布しているとし、観察領域中のある小領域５から全方向へ放射される光の一部は、レンズ２により光センサー素子列３に導かれる。ただし、ここで図３は、レンズ２の中心位置と観察領域１の中心との距離に比べて、レンズ２の直径および観察領域１の大きさが無視できない場合であるため、小領域５の位置によって、小領域５からの光は、レンズ２により光センサー素子列３に合焦せず、光センサー素子列３上の広範な領域に導かれ、非合焦状態による「ボケ」が生じる。
【００３５】
非合焦による「ボケ」の発生を、詳しく説明すると以下のようになる。説明図を図４に示す。
【００３６】
図４は、観察領域１にある一つの小領域５からの光を、レンズ２の開口部全体で捕らえた様子を示している。図４（ａ）は、小領域５とレンズ２が生み出す合焦面１７が一次元光センサー素子列３と一致し、小領域５とレンズ２と一次元光センサー素子列３とが合焦状態に配置されている場合を示しており、光は一次元光センサー素子列３の上で結像する。その結果、図４（ａ）の下部に示すように、プロジェクションデータ１８は、鋭いピークを持つ分布となる。
【００３７】
一方、図４（ｂ）は、小領域５が図４（ａ）の位置よりレンズ２に近い場合であり、小領域５とレンズ２が生み出す合焦面１７が一次元光センサー素子列３より遠くに位置する。光は本来ならば合焦面１７上に結像するため、一次元光センサー素子列３上では広範に広がった輝度分布となり、図４（ｂ）の下部に示すように、プロジェクションデータ１８も広がった分布となる。
【００３８】
図４（ｃ）は、小領域５が図４（ａ）の位置よりレンズ２から遠い場合であり、小領域５とレンズ２が生み出す合焦面１７が一次元光センサー素子列３よりレンズ近くに位置する。光は合焦面１７上に結像するため、一次元光センサー素子列３上では広範に広がった輝度分布となり、図４（ｃ）の下部に示すように、プロジェクションデータ１８も広がった分布となる。
【００３９】
以上が、観察領域１にある一つの小領域５からの光を、レンズ２の開口部全体で捕らえたときの、合焦状態および非合焦状態の場合の、プロジェクション獲得の様子である。
【００４０】
つぎに図４と同条件の場合について、レンズ各部を通過する光線について詳細に説明する。説明図を図５に示す。
【００４１】
図５（ａ）は、小領域１９とレンズ２が生み出す合焦面１７が一次元光センサー素子列３と一致し、小領域１９とレンズ２と一次元光センサー素子列３とが合焦状態に配置されている場合を示しており、小領域１９からの光は一次元光センサー素子列３の上で結像する。
【００４２】
その結果、図５（ａ）の下部に示すように、プロジェクションデータ２０は、鋭いピークを持つ分布となる。ここで、小領域１９からの光のうち、レンズ２の中心部を通過した光線２１は一次元光センサー素子列３上の結像位置に到達する。また、レンズ２の左部分を通過した光線２２も一次元光センサー素子列３上の結像位置に到達し、さらに、レンズ２の右部分を通過した光線２３も一次元光センサー素子列３上の結像位置に到達する。すなわち、レンズ２を通過する全ての光線は、一次元光センサー素子列３上の結像位置に到達する。
【００４３】
一方、図５（ｂ）は、小領域２４が図５（ａ）の位置よりレンズ２に近い場合であり、小領域２４とレンズ２が生み出す合焦面１７が一次元光センサー素子列３より遠くに位置する。光は本来ならば合焦面１７上に結像するため、一次元光センサー素子列３上では広範に広がった輝度分布となり、図５（ｂ）の下部に示すように、プロジェクションデータ２５も広がった分布となる。ここで、小領域２４からの光のうち、レンズ２の中心部を通過した光線２６は一次元光センサー素子列３上の結像の中央位置に到達するが、レンズ２の左部分を通過した光線２７は、一次元光センサー素子列３上の結像の中央位置より左側に到達する。また、レンズ２の右部分を通過した光線２８は、一次元光センサー素子列３上の結像の中央位置より右側に到達する。この他、他のレンズ部分を通過する光線も、一次元光センサー素子列３上の結像の中央位置以外に到達する。図４（ｂ）で示した非合焦状態の場合のプロジェクション分布の広がりは、このようにして発生する。
【００４４】
図５（ｃ）は、小領域２９が図５（ａ）の位置よりレンズ２から遠い場合であり、小領域２９とレンズ２が生み出す合焦面１７が一次元光センサー素子列３よりレンズ近くに位置する。光は合焦面１７上に結像するため、一次元光センサー素子列３上では広範に広がった輝度分布となり、図５（ｃ）の下部に示すように、プロジェクションデータ３０も広がった分布となる。ここで、当該の光束のうち、レンズ２の中心部を通過した光線３１は一次元光センサー素子列３上の結像の中央位置に到達するが、レンズ２の左部分を通過した光線３２は、一次元光センサー素子列３上の結像の中央位置より右側に到達する。また、レンズ２の右部分を通過した光線３３は、一次元光センサー素子列３上の結像の中央位置より左側に到達する。この他、他のレンズ部分を通過する光線も、一次元光センサー素子列３上の結像の中央位置以外に到達する。図４（ｃ）で示した非合焦状態の場合のプロジェクション分布の広がりは、このようにして発生する。
【００４５】
さらに、図４（ａ）〜（ｃ）および図５（ａ）〜（ｃ）で使用した、位置の異なる３つの小領域１９、２４、２９が同時に存在し、等しい強さの光を放射している場合について、レンズ各部を通過する光線について詳細に記述する。説明図を図６に示す。
【００４６】
図６（ａ）は、位置の異なる３つの小領域１９、２４、２９に輝度が同時に存在する場合であるので、図５（ａ）〜（ｃ）のレンズ各部を通過する光線を合成して、得られる図である。図６（ａ）のように、一次元光センサー素子列３上には、３つの像（２０、２５、３０）を合成したプロジェクション分布３４が得られる。
【００４７】
ここで、任意の光センサー素子に到達する光の経路について考える。一次元光センサー素子列３上の任意の光センサー素子３５に到達する光線のみを、図６（ｂ）に図示する。説明を分かりやすくするため、図６（ｂ）では、結像の中央位置付近の光センサー素子３５を例にとっている。
【００４８】
光センサー素子３５には、小領域１９からの光は、光線２１、２２、２３に代表されるように、レンズ２の開口の全面を通って到達する。一方、小領域２４からの光は、レンズ２の開口のうち極中央部を通る光のみが到達する。同様に、小領域２９からの光も、レンズ２の開口のうち極中央部を通る光のみが到達する。以上のように、光センサー素子３５に到達し、プロジェクションデータとして記録される光は、レンズ２を均一な分布で通過してくるのではなく、レンズ開口領域３６に対して開口分布３７のような分布で通過してくることが分かる。
【００４９】
したがって、任意の光センサー素子（たとえば３５）で獲得したプロジェクションデータ値を、バックプロジェクションビームを用いて再構成領域に分配させる場合、はじめに、当該の光センサー素子に対して開口分布３７を算出しておき、当該光センサー素子３５からレンズ開口全面に向けて多数の分割バックプロジェクションビームを逆投影させればよい。ただし、このとき、当該光センサー素子３５が獲得したプロジェクションデータ値は、各分割バックプロジェクションビームに、開口分布に比例させて分乗させて、この分割プロジェクションデータの値を、前述したMLEM法の手順に従い、分割バックプロジェクションビームが通過する小領域に分配すればよい。
【００５０】
開口が大きく、観察領域に近接したレンズを用いて獲得したプロジェクションデータ群から、観察領域の輝度分布を再構成する本実施形態の再構成手順の説明を行う。説明図として図７を示す。
【００５１】
この再構成は、データ取得・ＣＴ演算コンピュータ３７の内部で仮想的に演算処理される。データ取得・ＣＴ演算コンピュータ３７には、開口が大きく、観察領域に近接したレンズを用いて獲得したプロジェクションデータ群３８が記録されている。
【００５２】
はじめに、データ取得・ＣＴ演算コンピュータ３７の内部の仮想空間に、プロジェクションデータ群３８の取得時（図３参照）と同じ配置となるように、プロジェクションデータエリア３９を一次元光センサー列３（図３参照）の位置に配置し、レンズ４０をレンズ２（図３参照）の位置に配置し、再構成領域４１を観察領域１（図３参照）の領域に配置する。再構成領域４１は通常、細かな正方ピクセル（例えば４２）に分割しておく。正方ピクセルの数をＳとし、各正方ピクセルには、１からＳまでの番地を付けておく。また、プロジェクションデータエリア３９も、光センサー素子列３（図３参照）の各光センサー領域に対応するように、プロジェクション素４３に分割しておく。プロジェクション素の数はＭとし、各プロジェクション素には、１からＭまでの番地をつけておく。
【００５３】
データ取得・ＣＴ演算コンピュータ３７に記録されているプロジェクションデータ群３８を、プロジェクションデータエリア３９の各プロジェクション素４３に配置する。ここで、プロジェクションデータの分布のイメージ４４を、プロジェクションデータエリアに沿うように、プロジェクションデータエリアに垂直の距離で図示した。
【００５４】
さて、プロジェクションデータ群３８を再構成領域４１の各正方ピクセル（例えば４２）に再配置し、各正方ピクセルの持つＣＴ値（ｓ番目の正方ピクセルのＣＴ値をＣ(s)とする）を再構成する工程に移る。
【００５５】
本実施形態の再構成アルゴリズムは、繰り返し演算法であるMLEM法（DBP法）（前述参照）を基礎とする。本実施形態の再構成アルゴリズムは、繰り返し演算法である。最初の繰り返し演算の前に、初期値として、再構成領域４１の全ての正方ピクセル（例えば４２）に、等しいCT値を配置しておく。また、各プロジェクション素ごとに存在する開口分布Ｂも、各プロジェクション素ごとに、対応するレンズに関する合計が１となるように均一の値を代入しておく。
【００５６】
本実施形態のアルゴリズムでは、あるm番目のプロジェクション素（例えば４５）に配置されたプロジェクションデータ４６（値はP(m)）を、プロジェクション素４５からレンズ４０の開口領域を細かく分割したレンズ素（例えば４７）を通過し屈折して伸びる多数のバックプロジェクションビーム（例えば４８）に沿って、再構成領域４１に下記手順により配置する。ただし、各バックプロジェクションビーム（例えば４８）には、当該プロジェクション素（ここでは４５）に配置されたプロジェクションデータ４６（値はP(m)）を、当該のプロジェクション素（ここでは４５）に関する開口分布４９に比例させて分乗させる。この分乗プロジェクションデータの値を、前述したMLEM法の手順に従い、ビームが通過する正方ピクセル（例えば４２）に分配する。
【００５７】
この工程を具体的に以下に記述する。
【００５８】
合計でＭ個存在するプロジェクション素のうち、代表として、m番目のプロジェクション素４５に配置されたプロジェクションデータ４６（値はP(m)）を、再構成領域４１の正方ピクセル（例えば４２）に配置することを考える。各レンズを１番目からＥ番目までの合計Ｅ個のレンズ素（例えば４７）に分割しておく。m番目のプロジェクション素４５に関する、e番目のレンズ素の開口分布をB(m,e)とする。開口分布B(m,e)の決定方法は後述するが、プロジェクション素（例えば４５）毎の、対応するレンズにおけるB(m,e)の合計は１である。
【００５９】
m番目のプロジェクション素４５と対応するレンズのe番目のレンズ素（例えば４７）とに関するバックプロジェクションビーム（例えば４８）に分乗させる分乗プロジェクションデータQ(m,e)は、Q(m,e)＝P(m)×B(m,e)となる。
【００６０】
つぎに、当該のバックプロジェクションビーム（たとえば４８）が通過する全ての通過正方ピクセルを選び出す。その個数がＮ(m,e)個であるとする。１番目からＮ(m,e)番目の通過正方ピクセルの番地を、それぞれ、s(m,e,1)からs(m, e,N(m,e))とする。つぎに、バックプロジェクションビーム（たとえば４８）とその通過正方ピクセルが重なる面積A(m,e,1)〜A(m,e,N(m))を算出し、その通過正方ピクセルに配置されているCT値C(s(m,e,1))〜C(s(m,e,N(m,e))）を用いて、次の演算を行い、m番目のプロジェクション素４５と対応するレンズのe番目のレンズ素（例えば４７）とに関するバックプロジェクションビーム（例えば４８）がｉ番目に通過する通過正方ピクセル（番地はs(m,e,i)）に分配するCT値の分配量 C(m,e,i)を得る。
【００６１】
C(m,e,i)＝Q(m,e)×（A(m,e,i)×C(s(m,e,i))）／（ A(m,e,1)×C(s(m,e,1))＋・・・＋A(m,e,N(m,e))×C(s(m,e,N(m,e))）
すなわち本実施形態のアルゴリズムでは、分乗プロジェクションデータQ(m,e)を、そのバックプロジェクションビームが通過する全ての正方ピクセルに、それぞれの正方ピクセルのCT値とビーム通過面積との積に応じて分配する。
【００６２】
１番目からＭ番目までの全てのプロジェクション素（例えば４５）と該当するレンズの１番目からＳ番目までの全てのレンズ素（例えば４７）について上記の分配量を算出し、全ての正方ピクセルについて分配量 C(m,e,i)を合計することで、全ての正方ピクセル（例えば４２）についての新しいCT値が計算される。
【００６３】
また、次回の繰り返し演算で使用する新しい開口分布B(m,e)は、次式で算出する。ただし、各正方ピクセルのCT値は、上記で計算された最新のものを用いる。
【００６４】
B(m,e)＝［A(m,e,1)×C(s(m,e,1))＋・・・＋A(m,e,N(m,e))×C(s(m,e,N(m,e))］／｛［A(m,e=1,1)×C(s(m,e=1,1))＋・・・＋A(m, e=1,N(m, e=1))×C(s(m, e=1,N(m, e=1))］＋・・・＋［A(m,e=E,1)×C(s(m, e=E,1))＋・・・＋A(m, e=E,N(m, e=E))×C(s(m, e=E,N(m, e=E))］｝
すなわち、あるプロジェクション素（例えば４５）の開口分布B(m,e)は、はじめに、レンズ素ごとの通過プロジェクションビーム（例えば４８）に関して、通過ピクセルの再構成輝度と通過面積との積の合計値を個別合計値として算出しておき、つぎに、全てのレンズ素に関してする総合計を行い総合計値を算出しておくと、開口分布B(m,e)は、総合計値に対する個別合計値の割合としてそのレンズ素の開口分布値が算出される。
【００６５】
以上で、１回の繰り返し演算が終了する。繰り返し演算を適切な回数だけ繰り返すことで、再構成領域１０（図３参照）の各正方ピクセルのCT値の分布が完成し、プロジェクションデータ群からの再構成領域の輝度分布が再構成される。
【００６６】
上記した本実施形態の再構成アルゴリズムのフローチャートを図８に示す。
【００６７】
以上述べたように、本実施形態は、二次元CT再構成の場合であり、再構成領域は平面、多次元光センサーは一次元光センサー、再構成領域要素は正方平面ピクセルとなる。
【００６８】
なお、プロジェクションデータ群の獲得の際には、多数のカメラを用意してもよいし、観察対象を回転させてもよい。
【００６９】
なお、プロジェクションデータ群の獲得の際には、多数のカメラもしくは多レンズカメラを用いても良く、１台のカメラを円周方向に位置を変えて用いても良い。また、カメラを固定し、観察対象を回転させても良い。この場合、再構成時には、観察対象を固定し、多数のカメラが存在するかのように再構成場を構成する。
（第２実施形態）
第２実施形態は、三次元CT再構成の場合であり、再構成領域は立体、多次元光センサー３は二次元光センサー、再構成領域要素５は立体ボクセルとなる。
【００７０】
以下に、開口が大きく、観察領域に近接したレンズを用いたプロジェクションの取得手順と、本実施形態のＣＴアルゴリズムによるプロジェクションデータからの再構成手順を、いずれも、三次元の場合について記述する。
【００７１】
開口が大きく、観察領域に近接したレンズを用いたプロジェクション取得手順は、図３と同様に行う。放射型CTにおけるプロジェクションデータの取得は、通常の「写真撮影」と等しい。
【００７２】
第２実施形態の場合、図３において、１は観察領域、２はレンズ、３は二次元光センサー素子列、４はカメラ、５は観察領域中の任意小領域（ボクセル）、６はデータ取得・ＣＴ演算コンピュータである。
【００７３】
また、プロジェクション取得は、最低でも180゜に渡る方が望ましく、なるべく多方向である方が望ましいが、ここでは、図３に示すように４方向であるとして説明を行う。方向の配置は、同一面内には限らず、立体的な方法を設定してもよい。
【００７４】
プロジェクションデータ群の獲得における「ボケ」の発生原理は、第１実施形態に二次元の場合を詳述したが、三次元の場合である第２実施形態においても「ボケ」の発生原理は変わらない。
【００７５】
開口が大きく、観察領域に近接したレンズを用いて獲得したプロジェクションデータ群から、観察領域の輝度分布を再構成する本実施形態の再構成手順の説明を行う。説明図として図９を示す。
【００７６】
この再構成は、データ取得・ＣＴ演算コンピュータ５０の内部で仮想的に演算処理される。データ取得・ＣＴ演算コンピュータ５０には、開口が大きく、観察領域に近接したレンズを用いて獲得したプロジェクションデータ群５１が記録されている。
【００７７】
はじめに、データ取得・ＣＴ演算コンピュータ５０の内部の仮想空間に、プロジェクションデータ群５１の取得時（図３参照）と同じ配置となるように、プロジェクションデータエリア５２を二次元光センサー列（例えば図３の３番、参照）の位置に配置し、レンズ５３をレンズ（例えば図３の２番、参照）の位置に配置し、再構成領域５４を観察領域（例えば図３の１番、参照）の領域に配置する。再構成領域５４は通常、細かな立方ボクセル（例えば５５）に分割しておく。立方ボクセルの数をＳとし、各立方ボクセルには、１からＳまでの番地を付けておく。また、プロジェクションデータエリア５２も、光センサー素子列（例えば図３の３番、参照）の各光センサー領域に対応するように、プロジェクション素５６に分割しておく。プロジェクション素の数はＭとし、各プロジェクション素には、１からＭまでの番地をつけておく。
【００７８】
データ取得・ＣＴ演算コンピュータ５０に記録されているプロジェクションデータ群５１を、プロジェクションデータエリア５２の各プロジェクション素５６に配置する。ここで、プロジェクションデータの分布のイメージ５７を、プロジェクションデータエリアに沿うように、プロジェクションデータエリアに垂直の距離で図示した。
【００７９】
さて、プロジェクションデータ群５１を再構成領域５４の各立方ボクセル（例えば５５）に再配置し、各立方ボクセルの持つＣＴ値（ｓ番目の立方ボクセルのＣＴ値をＣ(s)とする）を再構成する工程に移る。
【００８０】
本実施形態の再構成アルゴリズムは、繰り返し演算法であるMLEM法（DBP法）（前述参照）を基礎とする。本実施形態の再構成アルゴリズムは、繰り返し演算法である。最初の繰り返し演算の前に、初期値として、再構成領域５４の全ての立方ボクセル（例えば５５）に、等しいCT値を配置しておく。また、各プロジェクション素ごとに存在する開口分布Ｂも、各プロジェクション素ごとに、対応するレンズに関する合計が１となるように均一の値を代入しておく。
【００８１】
本実施形態のアルゴリズムでは、あるm番目のプロジェクション素（例えば５８）に配置されたプロジェクションデータ５９（値はP(m)）を、プロジェクション素５８からレンズ５３の開口領域を細かく分割したレンズ素（例えば６０）を通過し屈折して伸びる多数のバックプロジェクションビーム（例えば６１）に沿って、再構成領域５４に下記手順により配置する。ただし、各バックプロジェクションビーム（例えば６１）には、当該プロジェクション素（ここでは５８）に配置されたプロジェクションデータ５９（値はP(m)）を、当該のプロジェクション素（ここでは５８）に関する開口分布６２に比例させて分乗させる。この分乗プロジェクションデータの値を、前述したMLEM法の手順に従い、ビームが通過する立方ボクセル（例えば５５）に分配する。
【００８２】
この工程を具体的に以下に記述する
合計でＭ個存在するプロジェクション素のうち、代表として、m番目のプロジェクション素５８に配置されたプロジェクションデータ５９（値はP(m)）を、再構成領域５４の立方ボクセル（例えば５５）に配置することを考える。各レンズを１番目からＥ番目までの合計Ｅ個のレンズ素（例えば６０）に分割しておく。m番目のプロジェクション素５８に関する、e番目のレンズ素の開口分布をB(m,e)とする。開口分布B(m,e)の決定方法は後述するが、プロジェクション素（例えば５８）毎の、対応するレンズにおけるB(m,e)の合計は１である。
【００８３】
m番目のプロジェクション素５８と対応するレンズのe番目のレンズ素（例えば６０）とに関するバックプロジェクションビーム（例えば６１）に分乗させる分乗プロジェクションデータQ(m,e)は，Q(m,e)＝P(m)×B(m,e)となる。
【００８４】
つぎに、当該のバックプロジェクションビーム（たとえば６１）が通過する全ての通過立方ボクセルを選び出す。その個数がＮ(m,e)個であるとする。１番目からＮ(m,e)番目の通過立方ボクセルの番地を、それぞれ、s(m,e,1)からs(m, e,N(m,e))とする。つぎに、バックプロジェクションビーム（たとえば６１）とその通過立方ボクセルが重なる体積A(m,e,1)〜A(m,e,N(m))を算出し、その通過立方ボクセルに配置されているCT値C(s(m,e,1))〜C(s(m,e,N(m,e))）を用いて、次の演算を行い、m番目のプロジェクション素５８と対応するレンズのe番目のレンズ素（例えば６０）とに関するバックプロジェクションビーム（例えば６１）がｉ番目に通過する通過立方ボクセル（番地はs(m,e,i)）に分配するCT値の分配量 C(m,e,i)を得る。
【００８５】
C(m,e,i)＝Q(m,e)×（A(m,e,i)×C(s(m,e,i))）／（ A(m,e,1)×C(s(m,e,1))＋・・・＋A(m,e,N(m,e))×C(s(m,e,N(m,e))  ）
すなわち本実施形態のアルゴリズムでは、分乗プロジェクションデータQ(m,e)を、そのバックプロジェクションビームが通過する全ての立方ボクセルに、それぞれの立方ボクセルのCT値とビーム通過体積との積に応じて分配する。
【００８６】
１番目からＭ番目までの全てのプロジェクション素（例えば５８）と該当するレンズの１番目からＳ番目までの全てのレンズ素（例えば６０）について上記の分配量を算出し、全ての立方ボクセルについて分配量 C(m,e,i)を合計することで、全ての立方ボクセル（例えば５５）についての新しいCT値が計算される。
【００８７】
また、次回の繰り返し演算で使用する新しい開口分布B(m,e)は、次式で算出する。ただし、各立方ボクセルのCT値は、上記で計算された最新のものを用いる。
【００８８】
B(m,e)＝［A(m,e,1)×C(s(m,e,1))＋・・・＋A(m,e,N(m,e))×C(s(m,e,N(m,e))］／｛［A(m,e=1,1)×C(s(m,e=1,1))＋・・・＋A(m, e=1,N(m, e=1))×C(s(m, e=1,N(m, e=1))］＋・・・＋［A(m,e=E,1)×C(s(m, e=E,1))＋・・・＋A(m, e=E,N(m, e=E))×C(s(m, e=E,N(m, e=E))］｝
すなわち、あるプロジェクション素（例えば５８）の開口分布B(m,e)は、はじめに、レンズ素ごとの通過プロジェクションビーム（例えば６１）に関して、通過ピクセルの再構成輝度と通過体積との積の合計値を個別合計値として算出しておき、つぎに、全てのレンズ素に関してする総合計を行い総合計値を算出しておくと、開口分布B(m,e)は、総合計値に対する個別合計値の割合としてそのレンズ素の開口分布値が算出される。
【００８９】
以上で、１回の繰り返し演算が終了する。繰り返し演算を適切な回数だけ繰り返すことで、再構成領域（例えば図３の１０番参照）の各立方ボクセルのCT値の分布が完成し、プロジェクションデータ群からの再構成領域の輝度分布が再構成される。
【００９０】
上記した本実施形態の再構成アルゴリズムのフローチャートを図１０に示す。
【００９１】
なお、プロジェクションデータ群の獲得の際には、多数のカメラもしくは多レンズカメラを用いても良く、１台のカメラを円周方向あるいは球周方向に位置を変えて用いても良い。また、カメラを固定し、観察対象を回転させても良い。この場合、再構成時には、観察対象を固定し、多数のカメラが存在するかのように再構成場を構成する。
（第３実施形態）
第１実施形態と同様に、二次元再構成を行うが、観察対象の吸収量の分布を求めるため、観察領域に対して、レンズの反対側に拡散背景放射を行う装置（たとえば，放射種類が可視光ならば、拡散光を放つ衝立で、可能）を配置し、そこからの放射に対して、観察対象が吸収した吸収後の透過量分布を近接した大開口のレンズによりプロジェクションデータ群を得る。このプロジェクションデータ群から、第１実施形態で使用した再構成アルゴリズムにより、観察領域の吸収量分布を得る。ただし、プロジェクションデータがもたらす吸収量は背景放射経路に沿った単位長さ当たりの吸収量の積であるので、第１実施形態で扱った輝度の代わりに、単位長さ当たりの吸収率の対数を使用して、背景放射経路に沿った単位長さ当たりの吸収量の対数の和として、再構成の演算を行う。
（第４実施形態）
第２実施形態と同様に、三次元再構成を行うが、観察対象の吸収量の分布を求めるため、観察領域に対して、レンズの反対側に拡散背景放射を行う装置（たとえば，放射種類が可視光ならば、拡散光を放つ衝立で、可能）を配置し、そこからの放射に対して、観察対象が吸収した吸収後の透過量分布を近接した大開口のレンズによりプロジェクションデータ群を得る。このプロジェクションデータ群から、第２実施形態で使用した再構成アルゴリズムにより、観察領域の吸収量分布を得る。ただし、プロジェクションデータがもたらす吸収量は背景放射経路に沿った単位長さ当たりの吸収量の積であるので、第２実施形態で扱った輝度の代わりに、単位長さ当たりの吸収率の対数を使用して、背景放射経路に沿った単位長さ当たりの吸収量の対数の和として、再構成の演算を行う。
【符号の説明】
【００９２】
１                観察領域
２                レンズ
３                センサー素子列
４                カメラ
５                観察領域中の任意小領域
６、３７、５０    データ取得・ＣＴ演算コンピュータ
７、３８、５１    プロジェクションデータ群
８、３９、５２    プロジェクションデータエリア
１０、４１、５４  再構成領域
１１、４２        正方ピクセル
１２、４３、５６  プロジェクション素
１３、４４、５７  プロジェクションデータの分布のイメージ
１６、４８、６１  バックプロジェクションビーム
５５              立方ボクセル
【図１】
2011102758.tif 000003
【図２】
2011102758.tif 000004
【図３】
2011102758.tif 000005
【図４】
2011102758.tif 000006
【図５】
2011102758.tif 000007
【図６】
2011102758.tif 000008
【図７】
2011102758.tif 000009
【図８】
2011102758.tif 000010
【図９】
2011102758.tif 000011
【図１０】
2011102758.tif 000012