WO2018179905A1

WO2018179905A1 - インテリアｃｔ画像生成方法

Info

Publication number: WO2018179905A1
Application number: PCT/JP2018/004603
Authority: WO
Inventors: 博幸工藤
Original assignee: 国立大学法人筑波大学
Priority date: 2017-03-27
Filing date: 2018-02-09
Publication date: 2018-10-04
Also published as: US20200029917A1; CN110520052A; EP3603520B1; EP3603520A4; US11457878B2; CN110520052B; JP7154611B2; JPWO2018179905A1; EP3603520A1

Abstract

先験情報が利用できない場合にも厳密解法により高精度な画像再構成が可能な汎用性の高いインテリアＣＴ画像生成方法を提供する。インテリアＣＴの画像生成方法であって、ＣＴ測定の幾何学系において撮像対象の内部の関心領域（ＲＯＩ）を通過する全ての量子ビームのみを測定してインテリアＣＴ投影データを取得する工程と、前記ＣＴ測定の幾何学系において前記撮像対象の外部にあるセグメントからの前記撮像対象の全体を通過する量子ビームを測定して部分的な全体投影データを取得する工程を含んでおり、更に、取得した前記インテリアＣＴ投影データと前記部分的な全体投影データに基づいて前記ＲＯＩを厳密に再構成する処理工程を含んでいる。

Description

インテリアＣＴ画像生成方法

　本発明は、物体内部における物理量分布の線積分値を測定してデータ処理により物理量分布を画像生成するＣＴの画像生成方法に関し、特に、インテリアＣＴの画像生成方法に関する。

　まず、インテリアＣＴ（コンピュータトモグラフィー）の原理について説明する。なお、以下では説明を分かりやすくするためＣＴ測定に使用する量子ビームがＸ線の場合を想定しているが、他の全ての量子ビームを用いたＣＴについても同様に適用可能なことは当業者であれば明らかである。ＣＴにおいて、その目的から、例えば心臓や乳房など小さな検査の関心領域（ＲＯＩ： Region of Interest）の画像のみで十分であり、ＲＯＩ外部の画像は不要である状況が想定される。通常のＣＴでは、このような状況であっても、図１（ａ）に示すように、ＲＯＩを通過するＸ線のみではなく、対象物全体をカバーするＸ線を照射して投影データを測定している。しかし、直感的にＲＯＩを通らない直線上の投影データはＲＯＩの情報を持っていないため、不要なことが予想される。そこで、このような無駄な測定をなくすため、図１（ｂ）に示すように、ＲＯＩを通過するＸ線のみを照射して投影データを測定してＲＯＩのみの画像を生成する方式のＣＴが、インテリアＣＴと呼ばれるものである。このようなインテリアＣＴには、（１）ＲＯＩ外部の被曝量の低減、（２）検出器やビーム幅の削減、（３）視野に収まらない大きな物体の撮影が可能になること、などの長所がある。また、マイクロ（ナノ）ＣＴなどの小視野を拡大して撮影するＣＴ装置では、試料が視野をはみ出してしまうインテリアＣＴの状況が頻繁に発生する。

　このように、インテリアＣＴにおいては、測定した投影データから画像を生成する画像再構成の問題は、以下のように定式化できる。図２に示すように、実線で示すＲＯＩ Ωを通る投影データは全て測定し、一方、点線で示すＲＯＩ Ωを通らない投影データは全て測定しないものとする。そして、このような不完全な投影データから、如何にしてＲＯＩ Ωの画像を数学的に厳密に再構成するかが、インテリアＣＴにおける画像再構成問題である。

　一方、インテリアＣＴに関する従来技術においては、まず、非特許文献[1]において、NattererはインテリアＣＴの画像再構成は解が一意に定まらず数学的に厳密な画像再構成は不可能であることを証明した。そこで、長年様々な近似解法が研究されてきたが、画像の低周波成分に誤差が発生して実用に至らなかった。例えば、代表的な画像再構成法として、各方向投影データの欠損部分を滑らかな関数で外挿してからフィルタ補正逆投影（ＦＢＰ： Filtered Backprojection）法により画像再構成を行う手法、不完全な投影データから代数的再構成（ＡＲＴ： Algebraic Reconstruction Technique）法や同時反復型逐次近似再構成（ＳＩＲＴ： Simultaneous Iterative Reconstruction Technique）法などの逐次近似法を適用する手法、などが研究された（非特許文献[2]～[4]）。

　図３に、インテリアＣＴにおいて発生する典型的なアーティファクトの例を示す。画像の低周波成分が歪むシェーディングアーティファクトやＲＯＩ周辺部の画像値が盛り上がるカッピングアーティファクトがよく発生する。これらの先行研究に対して、２００７年になりようやく数学的に厳密な画像再構成法が発見され、インテリアＣＴにおける画像再構成の新しい方向性が生まれた。これらの厳密解法の重要なキーポイントを一言で要約すると、『インテリアＣＴの投影データのみでは厳密な画像再構成は不可能だが、対象物に関するごく僅かな先験情報があれば数学的に厳密な画像再構成が可能になる』と述べることができる。

　更に、上述した厳密解法の発展について述べると、まず、２００７年のYeらの論文（非特許文献[5]）、２００８年のKudoらの論文（非特許文献[6]）、２００８年のWangらの特許（特許文献[1]）では、図４（ａ）に示すように、ＲＯＩ Ω内部の任意小領域B（先験情報領域）における画像値が既知という先験情報があれば、ＲＯＩ Ωの厳密な画像再構成が可能であることが証明された。次に、２００９年のYuらの論文（非特許文献[7]）と２０１４年のWangらの特許（特許文献[2]）では、図４（ｂ）に示すようにＲＯＩ Ω全体で画像値が区分的一様（完全な一定値を持つ有限個の領域で構成されていること）であれば、ＲＯＩ Ωの厳密な画像再構成が可能であることが証明された。更に、２０１０年のYangらの論文（非特許文献[8]）では、Yuらの手法を拡張して、図４（ｂ）に示すように、ＲＯＩ Ω全体で画像値が区分的多項式（有限次数の多項式で表される変化を持つ有限個の領域で構成されていること）であればＲＯＩ Ωの厳密な画像再構成が可能であることが証明された。最後に、２０１７年の工藤らの現在出願中の特許出願（ＰＣＴ出願）やKudoの論文（非特許文献[9]）では、Yuら、Wangら、Yangらの手法において必要な先験情報を大幅に削減することに成功して、ＲＯＩ Ω内部の任意小領域B（先験情報領域）において画像値が区分的一様または区分的多項式であれば、ＲＯＩ Ωの厳密な画像再構成が可能であることが証明された。

　加えて、非特許文献[10]では、KatsevichのＦＢＰ法と呼ばれるＦＢＰ型の構造を持つ画像再構成法が、非特許文献[11]では、KatsevichのＦＢＰ法におけるフィルタ補正投影データが対称性を持つことが、そして、非特許文献[12]では、位相ＣＴにも大きく異なる２つの原理が存在することが示されている。

米国特許第７，６９７，６５８号米国特許第８，８１１，７００号

Natterer F: The Mathematics of Computerized Tomography. Wiley, 1986 Ogawa K, Nakajima M, Yuta S: A reconstruction algorithm from truncated projections. IEEE Transactions on Medical Imaging 3: 34-40, 1984 Ohnesorge B, Flohr T, Schwarz K, Heiken JP, Bae KT: Efficient correction for CT image artifacts caused by objects extending outside the scan field of view. Medical Physics 27: 39-46, 2000 Ohyama N, Shiraishi A, Honda T, Tsujiuchi J: Analysis and improvement in region-of-interest tomography. Applied Optics 23: 4105-4110, 1984 Ye Y, Yu H, Wei Y, Wang G: A general local reconstruction approach based on a truncated Hilbert transform. International Journal of Biomedical Imaging 2007: Article ID 63634, 2007 Kudo H, Courdurier M, Noo F, Defrise M: Tiny a priori knowledge solves the interior problem in computed tomography. Physics in Medicine and Biology 53: 2207-2231, 2008 Yu H, Wang G: Compressed sensing based interior tomography. Physics in Medicine and Biology 54: 2791-2805, 2009 Yang J, Yu H, Jiang M, Wang G: High order total variation minimization for interior tomography. Inverse Problems, 26: Article ID 35013, 2010 Kudo H: Practical interior tomography. Proceedings of International Forum on Medical Imaging in Asia 2017, Paper No. O7-I-2, 2017 Katsevich A, Analysis of an exact inversion algorithm for spiral cone-beam CT. Physics in Medicine and Biology 47: 2583-2598, 2002 Noo F, Defrise M, Clackdoyle R, Kudo H: Image reconstruction from fan-beam projections on less than a short scan. Physics in Medicine and Biology 47: 2525-2546, 2002 百生敦: 位相コントラストＸ線イメージング. 放射光 10: 273-285, 1997

　上述した全ての厳密解法においては、数学的に厳密な画像再構成を可能にするためには、対象物に関する値が既知、区分的一様なことが既知、区分的多項式なことが既知、などの先験情報が必要になる。

　一方、現実のＣＴイメージングにおいて先験情報を獲得する方法としては、例えば、（１）同一患者の以前のＣＴ画像や他モダリティで撮影した画像からの推定、（２）他の患者の同じ部位を撮影したＣＴ画像からの推定、（３）上記の特許出願（ＰＣＴ出願）で提案された近似解法で再構成したアーティファクトを含む画像からユーザが手動または画像解析で推定するなど、多様な方法が考えられる。そして、これらの先験情報が得られるイメージング状況であれば、厳密解法により従来の近似解法よりはるかに高精度な画像再構成が可能である。しかし、撮影前に対象物に関する先験情報が既知であるという状況は必ずしも多くなく、一般性に欠けるものと言っても過言ではない。

　そこで、本発明は、上述した従来技術における課題、即ち、先験情報が利用できない場合においても、厳密解法により高精度な画像再構成が可能となるより実用的で汎用性にも優れたインテリアＣＴ画像生成方法を提供することをその目的とする。

　上記の目的を達成するため、本発明によれば、まず、インテリアＣＴの画像生成方法であって、ＣＴ測定の幾何学系において撮像対象の内部の関心領域（ＲＯＩ）を通過する全ての量子ビームのみを測定してインテリアＣＴ投影データを取得する工程と、前記ＣＴ測定の幾何学系において前記撮像対象の外部にあるセグメントからの前記撮像対象の全体を通過する量子ビームを測定して部分的な全体投影データを取得する工程を含んでおり、更に、取得した前記インテリアＣＴ投影データと前記部分的な全体投影データに基づいて前記ＲＯＩを厳密に再構成する処理工程を含んでいることを特徴とするインテリアＣＴの画像生成方法が提供される。

　本発明では、上記に記載のインテリアＣＴの画像生成方法において、前記セグメントは、前記撮像対象を取り囲む曲線に対応する複数の点を含んだ奇数個（１、３、５、…）のセグメントであることが好ましく、更には、前記インテリアＣＴ投影データ取得工程を、360度円軌道ファンビームで行い、前記部分的な全体投影データの取得を、前記インテリアＣＴ投影データを取得する前記360度円軌道ファンビームの円軌道内に含まれている前記セグメントからファンビームで行うことが、又は、前記インテリアＣＴ投影データ取得工程を、ファンビームショートスキャンで行い、前記部分的な全体投影データの取得を、前記インテリアＣＴ投影データを取得する前記ファンビームショートスキャンの円弧軌道内に含まれている前記セグメントからファンビームで行うことが、又は、前記インテリアＣＴ投影データ取得工程を、180度平行ビームスキャンで行い、前記部分的な全体投影データの取得を、前記インテリアＣＴ投影データを取得する前記180度平行ビームスキャンの軌道内に含まれている前記セグメントから平行ビームで行うことが好ましい。

　また、本発明では、上記に記載のインテリアＣＴの画像生成方法において、前記セグメントはその角度が、最低でも一つ以上の投影データが含まるように設定されており、更に、上記に記載のインテリアＣＴの画像生成方法において、前記ＲＯＩの再構成処理工程は、解析的画像再構成法、逐次近似画像再構成法、統計的画像再構成法の何れか、又は、それらの組み合わせにより実行されることが好ましく、更に、アクティブコリメータを利用して、前記インテリアＣＴ投影データと前記部分的な全体投影データを取得する量子ビームの開口角度を制御すること、又は、前記部分的な全体投影データとして、位置合わせのためのScout-Viewスキャン投影データを利用すること、又は、前記部分的な全体投影データを、前記インテリアＣＴ投影データよりも低い分解能の他の量子ビームを照射する測定を利用して取得することが好ましく、更に、量子ビームを対象物に照射して測定される前記インテリアＣＴ投影データ及び前記全体投影データが、量子ビームと物体の相互作用により発生する、吸収、位相シフト、散乱、回折、屈折の少なくとも何れか一つの空間的な物理量分布の、前記量子ビームが通過する直線上の線積分になっていることが好ましい。

　加えて、本発明によれば、上記の目的を達成するため、インテリアＣＴの画像生成方法であって、ＣＴ測定の幾何学系において撮像対象の内部の関心領域（ＲＯＩ）を通過する全ての量子ビームのみを測定してトランケーションを含んだインテリアＣＴ投影データを取得する工程と、前記ＣＴ測定の幾何学系において前記撮像対象の外部にあるセグメントからの前記撮像対象の全体を通過する量子ビームを測定してトランケーションを含まない部分的な全体投影データを取得する工程を含んでおり、更に、取得したトランケーションを含んだ前記インテリアＣＴ投影データとトランケーションを含まない前記部分的な全体投影データに基づいて前記ＲＯＩを厳密に再構成する処理工程を含んでいるインテリアＣＴの画像生成方法が提供される。

　また、本発明では、上記に記載のインテリアＣＴの画像生成方法において、前記ＲＯＩの再構成処理工程は、解析的画像再構成法、逐次近似画像再構成法、統計的画像再構成法の何れか、又は、それらの組み合わせにより実行されることが好ましく、又は、前記トランケーションを含んだ前記インテリアＣＴ投影データにおけるヒルベルト変換の計算不可能性を、前記トランケーションを含まない前記部分的な全体投影データを用いて計算可能にすることにより、前記ＲＯＩを再構成することが好ましい。

　上述した本発明によれば、先験情報が利用できない場合においても厳密解法により高精度な画像再構成が可能となる、実用性や汎用性においてより優れたインテリアＣＴ画像生成方法を提供するという優れた効果を発揮する。

通常のＣＴとの比較により本発明が関わるインテリアＣＴについて説明する図である。インテリアＣＴにおける画像再構成問題の定義を説明する図である。インテリアＣＴにおいて発生する典型的なアーティファクトの例を示す図である。画像値の先験情報が既知である、又は、画像値が区分的一様又は区分的多項式であればＲＯＩの厳密な画像再構成が可能であることを説明する図である。本発明の画像再構成方法である先験情報を利用することなくインテリアＣＴの数学的に厳密な画像再構成を実現する手法を説明する図である。本発明の画像再構成方法を一般化する場合の考え方を説明する図である。上記画像再構成方法を、180度平行ビーム、ファンビームショートスキャン、多角形軌道ファンビームスキャンに適用した場合の説明図である。通常のＣＴとインテリアＣＴの場合における360度円軌道ファンビームＣＴのデータ収集の様子を示す図である。非インテリアＣＴとインテリアＣＴの投影データに画像再構成法を適用した際の、違いや計算の破綻について検討した結果を示す図である。２つの対向位置の投影データ座標に対する、KatsevichのＦＢＰ法におけるフィルタ補正投影データの対称性を説明する図である。全体投影データの角度範囲Eが大きい場合の状況を説明する図である。全体投影データの角度範囲Eが任意の小さな円弧セグメントの場合の状況を説明する図である。Ｘ線源軌道、物体サポート、ＲＯＩ Ωが原点を中心とする円でない一般的な場合を考えた系１と系２の設定を説明する図である。本発明の有効性を示すために行った数値シミュレーション実験の結果を示す図である。本発明の有効性を示すために行った数値シミュレーション実験の結果を示す図である。本発明の有効性を示すために行った追加の数値シミュレーション実験の結果を示す図である。全体投影データを測定するセグメントを複数の小円弧セグメントに分割して設定する原理を説明する図である。回折格子を用いたＸ線位相顕微鏡やTalbot型干渉計による方法における幾何学的関係を示す図である。本発明の画像再構成方法を利用した装置の一例である、一般的なＸ線ＣＴ装置の外観構成を示す図である。上記Ｘ線ＣＴ装置の内部構成の一例を示す図である。

　以下、本発明の実施の形態について詳細に説明するが、それに先立ち、本発明者等が提案するインテリアＣＴ厳密解法について述べる。

＜本発明のインテリアＣＴ厳密解法＞
　本発明では、上述したように従来技術における対象物に関する先験情報を利用することなく、インテリアＣＴの数学的に厳密な画像再構成を実現する手法が提案されており、これについて以下に詳述する。

　インテリアＣＴにおいて解の一意性（画像再構成問題の解が一つに定まること）が成り立ち、厳密な画像再構成を可能にするためには、インテリアＣＴ投影データに加えて、何らかの付加的な情報が必要である。従来の厳密解法では、既述のように、対象物に関する先験情報を利用して解が一意に定まるようにしているが、本発明ではこれに代わるものとして、インテリアＣＴでは通常は測定しない投影データ、即ち、ＲＯＩ Ωを通らない余分な投影データを測定して、この補足データを利用する。即ち、提案する厳密解法の重要な特徴は、ＲＯＩ Ωを通らない必要最小限の余分な投影データを測定することにある。

　なお、上記の手法はＣＴの多様な幾何学系（Ｘ線源の軌道や光源からのＸ線の出し方等）にも適用できる一般性を有するが、以下では、最も実用的であるＸ線源が円軌道上を360度動くファンビームＣＴについてその原理を説明する。

　図５（ａ）に示すように、360度の円軌道のファンビームＣＴにおいて、ＲＯＩ Ωを通過するＸ線のみを照射して投影データを測定するインテリアＣＴの状況を考える。もちろん、インテリアＣＴ投影データだけでは数学的に厳密なＲＯＩ Ωの画像再構成は不可能であり、図５（ｂ）に示すように、付加的な情報として、円軌道の一部の円弧セグメントEから対象物全体をカバーするＸ線を照射して、全体の投影データを測定する。この部分的な全体投影データ（以下、単に「全体投影データ」とも言う）をインテリアＣＴ投影データに加味して、厳密なＲＯＩ Ωの再構成ができるようにする。この場合に問題になるのは、大きく以下の２点である。

（１）全体投影データを測定する円弧セグメントEは、どの程度短くできるか（もしも広範囲のEが必要であれば、非インテリアＣＴに近づいてしまい、提案手法の長所は小さいと思われる）。
（２）円弧セグメントEの全体投影データがあれば、厳密なＲＯＩ Ωの画像再構成が可能であることが数学的に証明できるか。

　本発明では、上記（１）の問題に対しては、驚くことに、Eは任意の小さな円弧セグメントでよいことが証明できた。

　また、上記（２）の問題に対しては、過去の対象物に関する先験情報を用いたインテリアＣＴ厳密解法の研究で解の一意性を示すツールとして使用された微分逆投影（ＤＢＰ： Differentiated Backprojection）法とトランケーションヒルベルト変換を組み合わせた画像再構成法（非特許文献[5]-[9]）では証明できないが、ＦＢＰ法におけるフィルタリング処理にトランケーションヒルベルト変換を導入した新しい画像再構成法を用いることで証明できた（詳細は後述する）。ここで、本発明で明らかになった360度円軌道ファンビームＣＴにおけるインテリアＣＴ画像再構成問題の解の一意性をまとめると、次のようになる。

［解の一意性］
　インテリアＣＴ投影データに加えて任意の円弧セグメントE（いくら小さくともよい）の全体スキャン（左右のトランケーションなし）投影データがあれば、インテリアＣＴの画像再構成の解は一意である。

　この新しいインテリアＣＴの厳密解法について、以下において、まず証明することに成功した解の一意性の結果を述べ、次に、その証明を示す。

＜得られたインテリアＣＴにおける解の一意性＞
　ここでは、後に説明するＦＢＰ法のフィルタリング処理にトランケーションヒルベルト変換を導入した新しい画像再構成法に基づき証明することに成功した、インテリアＣＴ画像再構成における解の一意性の結果をまとめて述べる。なお、記号の定義として、対象画像（物体）をf(x,y)、ＲＯＩをΩで表す。

（１）360度円軌道ファンビームＣＴの場合のインテリアＣＴの解の一意性について
　図５（ａ）に示す360度円軌道ファンビームＣＴにおいて、ＲＯＩ Ωを通過するＸ線のみを照射して投影データ測定を行うインテリアＣＴの状況を考える。なお、インテリアＣＴ投影データのみでは画像再構成問題の解は一意に定まらないので、図５（ｂ）に示すように、円軌道上の任意の小円弧セグメントEから対象物の全体をカバーするＸ線を照射して全体投影データを測定する。このとき、次の解の一意性が成り立つ。

［解の一意性（360度円軌道ファンビームスキャン）］
　ＲＯＩ Ωを通る全ての投影データに加えて、任意の小円弧セグメントE（いくら小さくともよい）から全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

　なお、この一意性の問題点は、ここでは、360度円軌道ファンビームＣＴに特化しているため、他の幾何学系（例えば、180度平行ビームスキャン、ファンビームショートスキャン、非円軌道を用いたファンビームスキャン）を用いたＣＴに適用できないことである。そこで、任意の幾何学系の場合におけるインテリアＣＴの解の一意性について、以下に述べる。

（２）任意の幾何学系の場合のインテリアＣＴの解の一意性
　そこで、360度円軌道ファンビームＣＴの場合の解の一意性を、任意の幾何学系で測定した投影データに適用可能なように一般化する。

　まず、図６に基づいて、一般化する際の考え方を述べる。任意の幾何学系を用いた場合においても、360度円軌道ファンビームＣＴの場合においても、物体の直線上の線積分値の集合を測定している点は共通であり、両者の投影データの間には座標変換（投影データ間の座標変換は、ＣＴ分野では、リビニングと呼ばれる）の関係が成立する。そこで、任意の幾何学系で測定した投影データを、一旦、図６（ａ）に示すように、仮想的な360度円軌道ファンビームＣＴの投影データに座標変換して、そのデータが360度円軌道ファンビームＣＴの場合における一意性の条件を満足しているか調べ、満足していれば解の一意性が成立すると判断するのが、基本となる考え方である。この考え方で任意の幾何学系で測定した投影データに適用可能な解の一意性を導出すると、最終的に、次の結論が得られる。

［解の一意性（任意の幾何学系）］
　以下の２つの両方の条件が満足されるように投影データが測定されていれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

　（条件１）ＲＯＩ Ωを通る全ての投影データが測定されていること（インテリアＣＴの測定条件）。

　（条件２）物体外部にある任意の（物体を内部に含むある円に対応する）小円弧セグメントEを通る全ての直線上の投影データ（部分的な全体投影データ）が測定されていること（解を一意に定めるために余分に測定する投影データの条件）。

　ただし、小円弧セグメントEの長さはいくら短くともよく、どの場所にあってもよい（但し、物体を内部に含むある円に対応するセグメントであることは必要）。

　上述した２つの幾何学的条件の意味を、図６（ｂ）に示す。

　更に、上述の一般化した解の一意性が有効な３つの事例を以下に述べる。いずれも、360度円軌道ファンビームＣＴの場合の解の一意性では、厳密なＲＯＩ Ωの再構成可能性が証明できない例になっている。

（ａ）180度平行ビームスキャンの場合
　ＣＴのデータ収集法の中で、最も基本的な180度平行ビームスキャンでインテリアＣＴを実施する場合を考える。動径をr、角度をθとして投影データをp(r,θ)（投影角度範囲-π/2≦θ＜π/2）で表す。いま、-ε≦θ≦ε（εは小角度）の角度範囲では（トランケーションなしの）全体投影データが測定され、それ以外ではインテリアＣＴ投影データしか測定されないとする。このとき、小円弧セグメントEを図７（ａ）に示すようにとれば上述の解の一意性の条件を満足していることが分かり、ＲＯＩ再構成の解は一意である。

［解の一意性（180度平行ビーム）］
　-π/2≦θ＜π/2の角度範囲のインテリアＣＴ投影データに加え任意の小角度範囲E（いくら小さくともよい）で全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

　なお、この一意性は、360度円軌道ファンビームCTの場合の解の一意性において円軌道の半径を無限にした極限と解釈することができ、180度平行ビームスキャンにおいても360度円軌道ファンビームＣＴと同様な解の一意性が成立する。なお、この結果は任意の幾何学系に一般化した一意性を用いて初めて証明できる。

（ｂ）ファンビームショートスキャンの場合
　図７（ｂ）に示すファンビームショートスキャンの場合を考える。円軌道上のＸ線源の位置をβ∈[-π/2-α_max,π/2+α_max)（α_maxはショートスキャンの条件から決まるオーバースキャン角度、非特許文献[11]）、直線検出器上の座標をuとしてファンビーム投影データをg(u,β)で表す。いま、-ε≦β≦ε（εは小角度）の角度範囲で（トランケーションなしの）全体投影データが測定され、それ以外ではインテリアＣＴ投影データしか測定されないとする。このとき、小円弧セグメントEを図７（ｂ）に示すようにとれば、上述の解の一意性の条件を満足していることが分かり、ＲＯＩ再構成の解は一意である。

［解の一意性（ファンビームショートスキャン）］
　-π/2-α_max≦β＜π/2+α_maxの角度範囲のインテリアＣＴ投影データに加えて、任意の小角度範囲E（いくら小さくともよい）で全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

（ｃ）多角形軌道を用いたファンビームスキャン：
　図７（ｃ）に示す正五角形軌道を用いたファンビームスキャンを考える。正五角形軌道上のＸ線源の位置をβ∈[-π,π)（βは正五角形の中心から軌道上の点を見た方位角）、直線検出器上の座標をuとしてファンビーム投影データをg(u,β)で表す。今、-ε≦β≦ε（εは小角度で正五角形軌道の一辺がEになるように決定）の角度範囲では（トランケーションなしの）全体投影データが測定され、それ以外ではインテリアＣＴ投影データしか測定されないとする。

　このとき、小円弧セグメントEを図７（ｃ）に示すようにとれば、上述の解の一意性の条件を満足していることが分かり、ＲＯＩ再構成の解は一意である。

［解の一意性（多角形軌道を用いたファンビームスキャン）］
　-π≦β＜πの角度範囲のインテリアＣＴ投影データに加えて、任意の小角度範囲E（いくら小さくともよい）で全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

　もちろん、上述した３つの事例は、任意の幾何学系を用いた解の一意性が有効であることを示す一端であり、本発明の手法により解の一意性が証明できるＣＴのデータ収集法の範囲はかなり広いことは、当業者であれば明らかであろう。

（３）測定の実施例
　上述したように、本発明の特徴は、通常のインテリアＣＴ投影データに加えて部分軌道Eから全体投影データを測定して、厳密な画像再構成を可能にする点にあり、このような測定を実現する実施例としては、各ＣＴ装置のハードウェア構成やイメージング状況に応じて、多様なものが考えられる。以下では、典型的な５つの実施例（ａ）～（ｅ）について述べる。

（ａ）純粋に２回の測定を行う
　１回目の投影データ測定は小角度範囲Eから物体全てをカバーするＸ線を用いて行い、２回目の測定はＲＯＩを通過するＸ線のみを照射してインテリアＣＴ投影データ全てが測定されるように行う。

（ｂ）位置合わせのScout-Viewスキャン投影データを全体投影データに利用する
　マイクロ（ナノ）CT装置などでは、撮影を行う前に見たいＲＯＩを上手く視野内に収める位置決めの目的で、物体全てが視野内に入るような低分解能のScout-Viewスキャンが行われる。このScout-Viewスキャンの投影データを、部分軌道Eの全体投影データとして利用することができる。

（ｃ）複数Ｘ線源を用いた装置
　近年、複数のＸ線源を用いた医療用Ｘ線ＣＴ装置やマイクロ（ナノ）ＣＴ装置が実用されている。これらの装置では、複数のＸ線源のうち１つを低分解能で全体投影データを測定するように、他のＸ線源でＲＯＩ Ωを高分解能で撮影するように設定すれば、低分解能投影データを部分軌道Eの全体投影データとして利用することができる。

（ｄ）アクティブコリメータを利用する方法
　開口角度を動的に変化可能なアクティブコリメータをＸ線源の前に設置して、部分軌道Eではトランケーションがないように開口角度を大きく、それ以外では開口角度を小さくしてＲＯＩを通過するＸ線のみを照射するように制御し、投影データ測定を行う。

（ｅ）以前に撮影した画像から計算で全体投影データを推定する方法
　同一の患者の以前に撮影したＣＴ画像や他の患者の同じ部位を撮影したＣＴ画像が利用できる状況であれば、これらの画像から部分軌道Eの投影データのうちトランケーションにより、未観測のデータを計算で推定することができる。

＜具体的な画像再構成法＞
　上述したように、ある幾何学系で測定した投影データが解の一意性の条件を満足していれば、具体的な画像再構成法としては（数学的に正しい方法であれば）どのような手法を用いても正しく画像再構成が可能である。したがって、本発明に対応する画像再構成法の選択肢は無数にある。具体的な例としては、以下の代表的な画像再構成法の全てが利用できる。

（１）解析的画像再構成法
　まず、360度円軌道ファンビームＣＴの場合には、後述する解の一意性を証明するのに用いたＦＢＰ法におけるフィルタリング処理にトランケーションヒルベルト変換を導入した新しい画像再構成法を用いることができる。360度円軌道ファンビームＣＴ以外の幾何学系の場合には、上述の画像再構成法を測定に用いた各々の幾何学系に拡張した直接再構成法、または、一旦360度ファンビーム投影データに座標変換（リビニング）した後に、上述の画像再構成法を適用するリビニング法を用いることができる。もちろん、本発明で証明に利用した画像再構成法と異なる解析的再構成法が将来発見されれば、それで置き換えることも可能であろう。

（２）逐次近似画像再構成法
　画像の画素値を一列に並べた画像ベクトルと測定したインテリアＣＴ投影データを一列に並べた投影データベクトルを、各々、x, bとすると、画像再構成問題は線形方程式Ax=bを解く問題と定式化できる。この線形方程式をＡＲＴ法、ＳＩＲＴ法、同時反復型ＡＲＴ（ＳＡＲＴ： Simultaneous Algebraic Reconstruction Technique）法などのＣＴ画像再構成分野における代表的な反復法を用いて解く。

（３）統計的画像再構成法
　画像xから計算される予測投影データAxと実測投影データbの距離distance(Ax,b)を雑音の統計的性質に基づいて定義しておき、この距離を凸最適化の手法を用いて最小化する。代表的な距離としては、重み付き最小２乗（ＷＬＳ： Weighted Least Squares）法、放射型ＣＴのPoisson分布対数尤度、透過型ＣＴのPoisson分布対数尤度などがある。距離の最小化に用いる凸最適化の反復法も多数のものが存在する。

　即ち、解の一意性が保証されているインテリアＣＴ測定データであれば、具体的な画像再構成法の選択には制限はなく（数学的に正しい方法であれば）、何れでもよいということである。

＜360度ファンビームＣＴの場合の解の一意性の証明＞
　以降では、本発明の基礎となっており、その実用可能性が最も高いと思われる360度円軌道ファンビームＣＴの場合の解の一意性を証明する。なお、任意の幾何学系に一般化する方法は既に説明している。なお、証明は、解の一意性の条件を満足する際に具体的に厳密な画像再構成が可能な手法を、２００２年に非特許文献[10]で発見されたKatsevichのＦＢＰ法と呼ばれるＦＢＰ型の構造を持つ画像再構成法に基づき示すことで行う。

（１）基礎事項
　まず、証明で考える問題設定と使用する用語・記号の定義を行い、証明に使用する画像再構成法であるKatsevichのＦＢＰ法について説明する。図８に、通常のＣＴとインテリアＣＴの場合における360度円軌道ファンビームＣＴのデータ収集の様子を示す。画像化の対象となる物体をf(x,y)で表す。そして、物体のサポート（対象物が存在する領域）は原点を中心とする半径aの円、ＲＯＩ Ωは原点を中心とする半径dの円の状況を考える（ただし、a>d）。また、Ｘ線源が動く軌道は原点を中心とする半径Rの円a (β)=(Rcos β,Rsin β)（β∈[-π,π)）とする（ただし、βはＸ線源の位置を表す変数でR>a>d）。そして、各Ｘ線源の位置βから測定されたファンビーム投影データを直線検出器上の座標uを用いて表す。検出器は対象物の外部に設置するが、画像再構成の研究でよく行われるように、実検出器と平行な原点を通る仮想検出器を考え、その座標uを用いて投影データを表す。即ち、ファンビーム投影データを(u,β)の２変数関数としてg(u,β)で表す。ただし、ＲＯＩ Ω、物体サポート、Ｘ線源軌道の全てが原点を中心とする円と仮定するのは、数式の論述を複雑にしないためであり、Ｘ線源軌道の内部に物体サポートが含まれ、物体サポートの内部にＲＯＩ Ωが含まれる幾何学的配置であれば、以下の証明は積分範囲を読み替えるなどの少ない変更で容易に拡張できることを述べておく（この点は証明の後に再度述べる）。

　通常のＣＴでは、物体全てをカバーするようにＸ線を照射するため、投影データはg(u,β)（-A≦u≦A, β∈[-π,π)）であるが（ただし、A=aR/(R²-a²)^1/2）、インテリアＣＴではＲＯＩ Ωを通過するＸ線のみが照射されるため、投影データは検出器座標uに関してより狭いg(u,β)（-D≦u≦D, β∈[-π,π)）の範囲（ただし、D= dR/(R²-d²)^1/2）でしか測定されない（d<aのとき明らかにD<Aである）。この各方向投影データg(u,β)の左右が欠損することをトランケーションと呼び、インテリアＣＴの画像再構成を困難にしている要因である。

　本発明では、通常の非インテリアＣＴを対象として２００２年に非特許文献[10]でKatsevichにより発見されたＦＢＰ法を証明に用いる。この画像再構成法の具体的な処理手順を、以下にまとめて示す。

［Step 1］微分フィルタ

［Step 2］ヒルベルト変換フィルタ

［Step 3］逆投影

　ただし、p.vは積分のCauchyの主値を表し、式(2)で表される積分変換はヒルベルト変換（Hilbert Transform）と呼ばれる。通常のＦＢＰ法では、再構成画像のぼけを補正するフィルタリングはランプフィルタにより一回の畳み込み積分で行われるが、KatsevichのＦＢＰ法では、ランプフィルタを微分とヒルベルト変換に分解してStep 1で微分をStep 2でヒルベルト変換を計算している点が異なる。このような特殊なＦＢＰ法を証明に用いる理由は、追って明らかになるが、ヒルベルト変換の巧妙な性質を利用するのが目的である。

（２）解の一意性の証明
　KatsevichのＦＢＰ法は、通常の非インテリアＣＴを想定した画像再構成法である。そこで、非インテリアＣＴとインテリアＣＴの投影データにこの画像再構成法を適用した際に、どのような違いがあり、インテリアＣＴではどこで計算が破綻するかを考えてみる。その結果をまとめて図９に表の形で示す。

　図９における（ａ）欄及び（ｂ）欄においては、微分は近傍データのみで計算できる局所演算であり、積分は積分範囲の全データが存在して初めて計算できることに注意して、上記の式(1)-(3)により計算を行ったときに各関数が計算できる範囲を示している。式(1)の微分は局所演算であるため、非インテリアＣＴでも、インテリアＣＴでも、何れの場合でも問題なく計算できる。しかし、式(2)のヒルベルト変換は、-A≦u≦Aの範囲でg(u,β)が測定されている非インテリアＣＴでは問題なく計算できるが、投影データの左右にトランケーションがあり-D≦u≦Dの範囲でしかg(u,β)が測定されていないインテリアＣＴでは計算できない。即ち、インテリアＣＴでは、ヒルベルト変換の箇所でＦＢＰ法の計算が破綻して続行できなくなる。

　そこで、本発明では、以下に述べるように、小さなＸ線源の角度範囲E={β|-ε≦β≦ε}（εは小角度）において（トランケーションなしの）全体投影データが測定されていれば、上述の計算不可能なヒルベルト変換の箇所を正しく計算してフィルタ補正された投影データg^F(u,β)（-D≦u≦D, β∈[-π,π)）を求め、逆投影してＲＯＩ Ωの数学的に厳密な画像再構成が可能なことを示す。

　図９における（ｃ）欄に、その新手法で計算できる各関数の範囲を示す。なお、以降では簡単のため、全体投影データが測定されるＸ線源の角度範囲Eは-ε≦β≦εとしているが、Eが他の角度範囲でも証明は全く同じである。さて、それを証明する。まず、証明に使用するヒルベルト変換の重要な２つの性質を述べる。一つ目は、『フィルタ補正投影データの対称性』と呼ばれる性質である。式(1),(2)で微分フィルタとヒルベルト変換フィルタを作用させたg^F(u,β)（いわばランプフィルタを作用させた投影データ）は、次の補題１に示す対称性を持つ。

［補題１］
　図１０（ａ）に示すように、(u₀,β₀)と(u₁,β₁)を２つの対向位置の投影データ座標とするとき、KatsevichのＦＢＰ法におけるフィルタ補正投影データg^F(u,β)は次の対称性を持つ。

　補題１の証明はNooらの非特許文献[11]を参照されたい。式(4)のような対称性が測定された元投影データg(u,β)について成立するのは自明だが、KatsevichのＦＢＰ法におけるフィルタリング後の投影データg^F(u,β)についても成立するのが、この驚きの性質である。次に、二つ目の性質は『トランケーションヒルベルト変換の計算可能性』と呼ばれるものである。いま、図１０（ｂ）に示すように、サポート（関数のゼロでない区間）[-A,A]に含まれる一部区間[-D,D]でしか観測されていない関数g(u)があり、その次式で定義されるヒルベルト変換を計算したい状況を考える。

　もちろん、通常式(5)を計算するには全区間[-A,A]のg(u)が必要であるが、一部区間[-D,D]でしか観測されていなくとも、ヒルベルト変換した結果Hg(u')が[-D,D]に含まれる任意の小区間u'∈P⊂[-D,D]（先験情報区間）で既知であれば、Hg(u')は[-D,D]全体で一意に定まり正しく計算できるというのがこの性質である。それを以下の補題２にまとめておく。

［補題２］
　サポートが[-A,A]の関数g(u)が[-D,D]⊂[-A,A]のみで観測され、かつ、g(u)のヒルベルト変換Hg(u')が[-D,D]に含まれる任意の小区間u'∈P⊂[-D,D]で既知であれば、Hg(u')は観測区間[-D,D]全体で一意に定まり観測データg(u)（-D≦u≦D）から正しく計算可能である。

　補題２の証明とその具体的な計算方法は、Kudoらの非特許文献[6]を参照されたい。非特許文献[6]では、凸射影法と呼ばれる先験情報を拘束条件として用いる反復法により、上述のトランケーションヒルベルト変換を計算している。別の非反復な計算方法として、行列の特異値分解を用いた擬似逆行列計算でも計算可能である。以上の２つの性質を用いて、目的の『X線源の小角度範囲-ε≦β≦ε（εは小角度）において（トランケーションなしの）全体投影データが測定されていれば、上述の計算不可能なヒルベルト変換の箇所を正しく計算してフィルタ補正された投影データg^F(u,β)（-D≦u≦D, β∈[-π,π)）を一意に求めることができる』事実を示す。証明は以下に述べる２段階の手順で行う。まず、第１ステップでは、全体投影データの角度範囲-ε≦β≦εが十分に大きい（εが大きい）場合について証明を行い、次に、第２ステップで、全体投影データの角度範囲-ε≦β≦εはいくらでも小さくできることを証明する。

（ａ）第１ステップ（全体投影データの角度範囲Eが大きい場合）
　まず、図１１に示すような状況を考える。図１１において、実線は全体投影データg(u,β)（-A≦u≦A）が測定されているＸ線源軌道の角度範囲E={β| -ε≦β≦ε}を表し、点線はインテリアＣＴ投影データg(u,β)（-D≦u≦D）のみが測定されている軌道の角度範囲N=[-π,π)＼Eを表す（ただし、小角度εはＲＯＩ Ωの境界に接する垂直線が軌道から切り取る円弧セグメントから決まる）。この場合、角度範囲Eのg(u,β)は全体投影データであるため、式(1),(2)で微分とヒルベルト変換したg^F(u,β)（-A≦u≦A）が問題なく計算できる。しかし、インテリアＣＴ投影データしか測定されていない角度範囲Nでは、微分は計算できるがヒルベルト変換は計算できない状況になる。これに対して、補題１と補題２を上手く組み合わせて用いると、次の手順によりg^F(u,β)（-D≦u≦D）が正しく計算できる。

　各β₀∈Nのトランケーション投影データに対して、以下のStep 1とStep 2を行う。

［Step 1］（ヒルベルト変換先験情報構築）
　β₀に対向する図１１中のB-Cの範囲の微分とヒルベルト変換が計算されているg^F(u,β)をβ₀側に折り返して集め補題1の対称性を用いると、g^F(u,β₀)の先験情報が構築できる（即ち、小区間u∈P⊂[-D,D]においてg^F(u,β₀)が既知という状況が構築できる）。

［Step 2］（トランケーションヒルベルト変換）
　Step 1で構築したg^F(u,β₀)（u∈P）を先験情報に用いて、補題２のトランケーションヒルベルト変換の手法により微分フィルタをかけたデータg^D(u,β₀)（-D≦u≦D）から、フィルタ補正されたデータg^F(u,β₀)（-D≦u≦D）を計算する。

　全体投影データの角度範囲Eが大きい図１１に示す設定の場合には、（空集合でない）先験情報区間Pを全てのβ∈Nについて構築でき、この方法で全てのβ∈Nについてg^F(u,β) （-D≦u≦D）を正しく計算できる。よって、フィルタ補正された投影データg^F(u,β)（-D≦u≦D, β∈[-π,π)）を式(3)により逆投影して、ＲＯＩ Ωの数学的に厳密な画像再構成が可能である。

（ｂ）第２ステップ（全体投影データの角度範囲Eが任意の小さな円弧セグメントの場合）
　第１ステップの証明だけでは、解の一意性の結論で述べた全体投影データの角度範囲が任意の小円弧セグメントEでよいことは示せない。このはるかに強い結果を示すには、もう一工夫が必要である。どのようにしてそれを証明するかを、図１２に示す。

　第１ステップで述べた対向データを利用してヒルベルト変換可能なＸ線源軌道の角度範囲を拡大する操作を、複数回繰り返すのが基本アイデアである。

　図１２の例では、元々ヒルベルト変換可能な角度範囲Eからスタートして、一回、第１ステップの操作を行うとヒルベルト変換可能な角度範囲E₁は真ん中の図まで増大し、もう一回、第１ステップの操作を行うと全てのβ∈[-π,π)についてヒルベルト変換可能なようになる。この例では２回の操作で十分だが、一般にEが空集合でない限り、無限回の操作で必ず360度角度範囲β∈[-π,π)についてヒルベルト変換できる状態に到達することが幾何学的に証明できる。よって、式(3)の逆投影によりＲＯＩ Ωの画像を生成するのに十分なフィルタ補正投影データg^F(u,β)（-D≦u≦D, β∈[-π,π)）が計算できるので、ＲＯＩ Ωを厳密に正しく再構成できる。

　最終的に得られた解の一意性に関する結論を定理の形でまとめると、以下のようになる。

［定理］
　図８（ｂ）の設定（Ｘ線源軌道は原点中心の半径Rの円、物体サポートは原点中心の半径aの円、ＲＯＩ Ωは原点中心の半径dの円）の360度円軌道ファンビームＣＴにおいて、ＲＯＩ Ωを通る全ての投影データに加えて任意の小円弧セグメントE（いくら小さくともよい）から全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

　上述の定理では、証明に現れる数式が複雑になり過ぎないように、Ｘ線源軌道、物体サポート、ＲＯＩ Ωが原点を中心とする円の場合を考えたが、この証明を拡張して以下の２つの系を証明することは容易である。具体的には、証明に用いる画像再構成法の積分範囲などが多少変わってくるが、証明の本質と手順は同じである。

［系１］
　図１３（ａ）の設定（Ｘ線源軌道は半径Rの円、物体サポートは円軌道の内部に含まれる任意の円、ＲＯＩ Ωは物体サポートの内部に含まれる任意の円）の360度円軌道ファンビームＣＴにおいて、ＲＯＩ Ωを通る全ての投影データに加えて任意の小円弧セグメントE（いくら小さくともよい）から全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

［系２］
　図１３（ｂ）の設定（Ｘ線源軌道は半径Rの円、物体サポートは円軌道の内部に含まれる任意の凸形状領域、ＲＯＩ Ωは物体サポートの内部に含まれる任意の凸形状領域）の360度円軌道ファンビームＣＴにおいて、ＲＯＩ Ωを通る全ての投影データに加えて任意の小円弧セグメントE（いくら小さくともよい）から全体投影データを測定すれば、ＲＯＩ Ωで画像f(x,y)は一意に定まり、Ωの厳密な再構成が可能である。

＜シミュレーション実験結果＞
　図１４及び図１５に、本発明の有効性を示すために行った数値シミュレーション実験の結果を示す。実験では、人体頭部のＣＴ画像を模擬した数値ファントムを用いて画像再構成を行った。360度あたりのファンビーム投影データ数（View数）は1,000であり、全体投影データを測定する小円弧セグメントE内の投影データ数MがM=300,M=100,M=33の３パターンについて、画像再構成を行った。なお、ここで、M=300はEの角度範囲が108度、M=100はEの角度範囲が36度、M=33はEの角度範囲が12度の場合に相当する。また、上述のセグメントEが１個の場合に加えて、セグメントEが３個の場合についても実験を行った。Eが３個のセグメントの場合は、M=94×3個（34度×3個）とM=25×3個（9度×3個）の設定とした。図１５（ａ）は、小円弧セグメントEが１個の場合の再構成画像であり、図１５（ｂ）は小円弧セグメントEが３個の場合の再構成画像である。　　　

　画像再構成に用いた手法は、解の一意性を証明するのに使用したKatsevichのＦＢＰ法におけるフィルタリング処理にトランケーションヒルベルト変換を導入して厳密に画像再構成できるようにした厳密解法である。また、比較のため、全ての投影データが全体投影データの場合の通常のＦＢＰ法、全ての投影データがインテリアＣＴ投影データの場合に最もよく使用されるローカルＦＢＰ法（投影データの欠損部分を滑らかな関数で外挿してからＦＢＰ法で画像再構成）についても実験を行った。

　これらの図からも明らかなように、ローカルＦＢＰ法では、ＲＯＩ周辺部の画像値が盛り上がるカッピングアーティファクトが強く発生しているのに対して、本発明のセグメントEの個数が１個の手法ではM=300とM=100についてはほぼ完全にＲＯＩを画像再構成することができている。なお、Eの角度範囲が小さなM=33の場合については、M=300やM=100と比較すると低周波数成分のアーティファクトは多少増大しているが、ローカルＦＢＰ法と比較すると、はるかに誤差が少ないことが分かる。また、セグメントEの個数が３個の場合の方が同等な角度範囲のEを用いる１個の場合と比較して、よりアーティファクトが少ない画像再構成が実現できる傾向が見られた。

＜追加シミュレーション実験結果＞
　図１６に、本発明の有効性を示すために行った追加の数値シミュレーション実験の結果を示す。実験では、医療用Ｘ線ＣＴ装置で測定した人体頭部のＣＴ実画像を数値ファントムとして画像再構成を行った。180度平行ビームスキャンを想定して、180度あたりの投影データ数（View数）は800と設定して、上述のセグメントEが１個、全体投影データを測定する小円弧セグメントE内の投影データ数MがM=1の場合（即ち全体投影データの数が僅か1個のみ）という最も難しい場合について、画像再構成を行った。なお、ここで、M=1はEの角度範囲が0.225度の場合に相当する。図１６（ａ）は原画像、図１６（ｂ）は全体投影データを用いない場合のローカルＦＢＰ法による再構成画像、図１６（ｃ）は全体投影データを用いない場合の逐次近似法による再構成画像、図１６（ｄ）は１個の全体投影データを用いる本発明の手法による再構成画像である。

　画像再構成に用いた手法は、段落［００５９］で述べた逐次近似画像再構成法である。具体的には、画像の画素値を一列に並べた画像ベクトルと測定したインテリアＣＴ投影データと１個の全体投影データを一列に並べた投影データベクトルを、各々、x, bとして、画像再構成問題を線形方程式Ax=bを解く問題に定式化して、この線形方程式をＣＴ画像再構成分野における代表的な反復法であるＡＲＴ法を用いて解いた。ただし、実装にあたっては、反復が進むにつれて計算誤差が増大する影響を抑制するため、トータルバリエーション（ＴＶ： Total Variation）による正則化（非特許文献[7]）を弱く施した。

　この図から明らかなように、全体投影データを用いない通常のインテリアＣＴにおいては、ローカルＦＢＰでも逐次近似法でもＲＯＩ周辺部の画素値が盛り上がるカッピング効果と画像の直流成分シフトによるアーティファクトが強く発生して脳内の構造が見えないほど画質が劣化しているのに対して、１個のセグメントEから僅か１個の全体投影データを測定して補足的に用いる本発明の手法ではＲＯＩをほぼ完全に画像再構成することができている。本実験では、全体投影データに誤差や雑音が全くない場合を想定したが、全体投影データの測定精度が高い場合には、本実験結果に見られるように１個のセグメントEから測定した僅か１個の全体投影データのみを補足して、インテリアＣＴにおいて発生するアーティファクトをほぼ完全に除去することができると結論づけられる。

＜複数の全体投影データを測定する円弧セグメントを用いた方法＞
　これまでの議論では、画像再構成問題の解を一意に定めるために全体投影データを測定する円弧セグメントEが単一の場合のみを考えてきた。この場合でも、上述の証明の第２ステップで行ったヒルベルト変換可能な角度範囲を拡大する操作を複数回逐次的に繰り返すことで、全体投影データを測定するセグメントEをいくらでも小さくできることは数学的に証明できるが、この第２ステップの操作回数が多いと画像再構成が複雑になり安定性（解が一意に定まりやすく雑音や計算誤差に敏感に影響されないこと）が悪くなることが危惧され、上述のシミュレーション実験でもEが小さいとアーティファクトが増大することが確認された。

　そこで、この問題点を緩和するため、ヒルベルト変換可能な角度範囲を拡大する操作が一回だけ（証明の第１ステップだけ）でも全体投影データを測定するセグメントEを十分に小さくできる方法として、Eを複数の小円弧セグメントに分割して設定する方法を考案した。

　図１７に、本手法の原理を示す。全体投影データを測定する円弧セグメントEが単一の場合には、ヒルベルト変換可能な角度範囲を拡大する操作を一回で済ませるには、（証明で述べたように）図１７（ａ）に示すようにEはかなり大きくとることが必要である。これに対して、図１７（ｂ），（ｃ）に示すように、全体投影データを測定する角度範囲Eを奇数n個（n=3,5,…）の等間隔に配置した小円弧セグメントに分割して設定すれば、どの方向のインテリアＣＴ投影データβ∈Nも対向データから一部が分かっている先験情報区間u∈Pが存在する状況を構築することがはるかに容易になり、トータルのEの角度範囲を大幅に小さくすることが可能である（図１７（ａ）と図１７（ｂ），（ｃ）の比較から明らかである）。即ち、全体投影データを測定するセグメントEを分割して配置することが可能であれば、より小さなEでより安定に画像再構成を行うことができる。

＜位相Ｘ線ＣＴへの拡張＞
　これまでは、主に物体内部のＸ線吸収係数分布を画像化する吸収ＣＴを想定して説明を行ってきたが、近年Ｘ線が物体を通過する際の位相シフト量分布を画像化する位相ＣＴの実用化が進んでいる。本発明で述べた解の一意性は、かかる位相ＣＴの場合にも全く同じように成立することが示される。以下にその理由を簡潔に述べる。

　位相ＣＴにも以下の大きく異なる２つの原理が存在する（非特許文献[12]）。

（１）単色Ｘ線を用いたBonse-Hart型干渉計による方法
　この方法では、測定した干渉縞パターンから位相復元処理を行って得られるデータは、位相シフト量空間分布の線積分値になる。よって、任意の幾何学系(γ,η)で測定した投影データp(γ,η)を仮想360度円軌道ファンビームＣＴの投影データg(u,β)に座標変換することができ、これまでの議論と同様に変換されたデータg(u,β)が解の一意性を満足することが厳密にＲＯＩ再構成可能な条件となる。したがって、解の一意性は全く同じであり、本発明の全ての結論が成立する。

（２）回折格子を用いたＸ線位相顕微鏡やTalbot型干渉計による方法
　この方法では、測定した干渉縞パターンやモアレ縞パターンから位相復元処理を行って得られる投影データは、位相シフト量空間分布の線積分値を各光線に垂直な方向に一次微分したデータ（微分位相投影データと呼ばれる）になる。いま、任意の幾何学系(γ,η)で測定した微分位相投影データをp^D(γ,η)で表し、各データが微分される方向の単位法線ベクトルをn(γ,η)とすると、各データp^D(γ,η)と同じ直線に対応するKatsevichのＦＢＰ法で現れる式(1)の微分フィルタ投影データg^D(u,β)の間には、次の関係式が成り立つ。

　ただし、a'(β)≡(d/dβ)a(β)=(-Rsin β,Rcos β)である。式(6)の証明は、Nooらの非特許文献[11]を参照されたい。なお、図１８には、式(6)の幾何学的関係を示す。

　したがって、式(6)を用いて任意の幾何学系で測定した微分位相投影データp^D(γ,η)を仮想360度円軌道ファンビームＣＴにおけるKatsevichのＦＢＰ法の微分投影データg^D(u,β)に変換することができ、これまでの議論と同様に、変換されたデータg^D(u,β)が解の一意性を満足することが厳密にＲＯＩ再構成可能な条件となる。したがって、解の一意性は全く同じであり本発明の全ての結論が成立する。

　以上に詳細に述べたように、本発明は、先験情報が既知であるという状況でない場合でも厳密解法によって従来の近似解法よりはるかに高精度な画像再構成が可能な、より一般性のあるインテリアＣＴ画像生成方法であり、当該方法は、物体内部における物理量分布の線積分値を測定してデータ処理により物理量分布を画像生成する原理に基づくあらゆるＣＴイメージング装置においても適用可能である。

　なお、ＣＴとは、通常、Ｘ線吸収係数分布を画像生成する吸収Ｘ線ＣＴのことを指す場合が多く、そこで、以下には、本発明の画像生成方法を適用する装置の最も一般的な例として、Ｘ線を用いて被検体の内部の断面像を得るＸ線ＣＴ装置について、図面を参照しながら概略を説明する。

　図１９は、上述した画像生成方法を利用して物体内部における物理量分布の線積分値を測定してデータ処理により物理量分布を画像生成する一般的なＸ線ＣＴ装置の全体外観構成を示す。即ち、Ｘ線ＣＴ装置は、以下にも述べるＸ線照射部等の構成要素を収納すると共に、その中央部には被検体が位置される略円筒状の空洞部を備えたガントリー部１と、その上面に被検体を載置する天板（クレードル）４を備えた基台部２と、そして、データ処理装置であるコンピュータ（ここでは図示せず）、得られた画像などを表示するためのディスプレイ装置５や必要な入力を行うためのキーボード６等を含むコンソール部３を備えている。

　ガントリー部１のハウジング内部やコンソール部３内には、図２０にも示すように、Ｘ線ＣＴ装置を構成する構成要素が設けられている。その一例として、図示のように、ハウジングの内部には、サンプルに対してＸ線を扇状に照射するＸ線発生装置１０と、当該装置から照射されて被検体を透過したＸ線を検出する円弧状のＸ線検出装置２０が、ここでは図示しないが、例えばリング状のフレーム上に取り付けられている。

　一方、Ｘ線発生装置１０とＸ線検出装置２０との間の空間には、その上面に被検体を載置（セッティング）する天板３０（図１９の符号４に対応）が設けられている。なお、その一部にＸ線発生装置１０とＸ線検出装置２０が取り付けられる部材は、回転駆動部５０を介して、例えば、上記ガントリー部１の内部に設けられたモータ等の回転駆動機構により、所定の方向に所定の回転速度で回転する（図の矢印を参照）。一方、被検体を載置するための天板３０は、上記Ｘ線発生装置１０とＸ線検出装置２０の回転面の略中央部の円筒状の空間に対向して配置され、サンプル載置台移動部６０により移動される。更に、ガントリー部１には、上記Ｘ線発生装置１０に対して高電圧を発生して供給するためのＸ線高電圧部４０、モータの回転制御によりＸ線発生装置１０とＸ線検出装置２０が取り付けられる部材を回転駆動する回転駆動部５０などが設けられている。

　上述したＸ線検出装置２０からの検出信号は、データ収集部７０に入力されて画像データとして収集され、更に、画像再生部７５において、サンプルの内部の断面画像または３次元画像等として再生される。なお、図中の符号７６は、画像再生部７５においてサンプルの内部の断面画像または３次元画像を再生する際に使用される記憶装置（画像メモリ）である。また、当該画像再生部７５により再生されたサンプルの内部の断面画像または３次元画像は、例えば、液晶表示装置等により構成された画像表示部８０（図１９の符号５に対応）において表示される。なお、この画像表示部８０に、所謂、タッチパネル（図示せず）を組み込むことによれば、当該画像表示部８０により、装置の操作に必要な入力を行うことが可能となる。しかしながら、本発明はこれに限定されることなく、当該装置は、当該タッチパネルに代え、キーボード（図１９の符号６に対応）やテンキーやマウス等を備えてもよい。

　図中の符号９０は、上述したＸ線ＣＴ装置を構成する各部の動作を制御するための制御部（図１９の符号３に対応）を示している。より具体的には、例えば、中央演算処理装置（ＣＰＵ）、そして、ＲＡＭやＲＯＭ等の記憶装置（メモリ）、更には、ＨＤＤ等の外部記憶装置等により構成されており、そして、記憶装置内に格納した各部の動作を制御するためのソフトウェアやファームウェアに基づいて必要な制御を実行する。

　そして、上述した本発明であるインテリアＣＴの画像生成方法は、上記Ｘ線ＣＴ装置を構成する画像再生部７５において、例えば、ソフトウェアとしてＲＡＭやＲＯＭ等の記憶装置（メモリ）内に格納され、中央演算処理装置（ＣＰＵ）によって実行される。

　なお、本発明が応用できる技術はこれに止まらず、例えば、Ｘ線を照射した際の位相シフト分布の線積分データから位相シフト分布を画像生成する位相Ｘ線ＣＴ、体内に投与した放射性薬剤分布の画像を生成する核医学イメージング装置であるＰＥＴ（ポジトロンエミッションＣＴ）やＳＰＥＣＴ（単光子放射型ＣＴ）、超音波・マイクロ波・音波・地震波などの波動を用いたＣＴ、電子線ＣＴ、投影データからの画像再構成を利用したＭＲＩ（磁気共鳴イメージング）などにも応用可能である。即ち、本発明において「物体」または「画像」とは、画像化する物理量の空間分布のことを指し、「投影データ」とは、その直線上の線積分値を表す測定データのことを指す。

　また、位相シフト、量子ビームの位相シフト、回折、又は回折を含む投影データの数値は、光学素子の追加あるいはその位置変更により取得した単数又は複数の投影データのセットから抽出され、当該抽出された前記量子ビームの位相シフト、回折、又は回折を含む前記投影データの数値を用いて画像を再構成することも可能である。

　以上には、本発明の種々の実施例になるインテリアＣＴの画像生成方法について詳細に述べた。しかしながら、本発明は、上述した実施例のみに限定されるものではなく、様々な変形例が含まれることは明らかである。例えば、上記した実施例は本発明を分かりやすく説明するためにシステム全体を詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能であろう。

産業上の利用性

　本発明は、物体内部における物理量分布の線積分値を測定してデータ処理により物理量分布を画像生成する画像再構成方法、特に、インテリアＣＴの画像再構成方法を提供する。

　１…ガントリー部、３…コンソール部、４、３０…天板、１０…Ｘ線発生装置、２０…Ｘ線検出装置、９０…制御部

Claims

　インテリアＣＴの画像生成方法であって、
　ＣＴ測定の幾何学系において撮像対象の内部の関心領域（ＲＯＩ）を通過する全ての量子ビームのみを測定してインテリアＣＴ投影データを取得する工程と、
　前記ＣＴ測定の幾何学系において前記撮像対象の外部にあるセグメントからの前記撮像対象の全体を通過する量子ビームを測定して部分的な全体投影データを取得する工程を含んでおり、更に、
　取得した前記インテリアＣＴ投影データと前記部分的な全体投影データに基づいて前記ＲＯＩを厳密に再構成する処理工程を含んでいることを特徴とするインテリアＣＴの画像生成方法。
　請求項１のインテリアＣＴの画像生成方法において、前記セグメントは、前記撮像対象を取り囲む曲線に対応する複数の点を含んだ奇数個（１、３、５、…）のセグメントであることを特徴とするインテリアＣＴの画像生成方法。
　請求項２のインテリアＣＴの画像生成方法において、前記インテリアＣＴ投影データ取得工程を、360度円軌道ファンビームで行い、前記部分的な全体投影データの取得を、前記インテリアＣＴ投影データを取得する前記360度円軌道ファンビームの円軌道内に含まれている前記セグメントからファンビームで行うことを特徴とするインテリアＣＴの画像生成方法。
　請求項２のインテリアＣＴの画像生成方法において、前記インテリアＣＴ投影データ取得工程を、ファンビームショートスキャンで行い、前記部分的な全体投影データの取得を、前記インテリアＣＴ投影データを取得する前記ファンビームショートスキャンの円弧軌道内に含まれている前記セグメントからファンビームで行うことを特徴とするインテリアＣＴの画像生成方法。
　請求項２のインテリアＣＴの画像生成方法において、前記インテリアＣＴ投影データ取得工程を、180度平行ビームスキャンで行い、前記部分的な全体投影データの取得を、前記インテリアＣＴ投影データを取得する前記180度平行ビームスキャンの軌道内に含まれている前記セグメントから平行ビームで行うことを特徴とするインテリアＣＴの画像生成方法。
　請求項３、４、５の何れか一項のインテリアＣＴの画像生成方法において、前記セグメントはその角度が、最低でも一つ以上の投影データが含まれるように設定されていることを特徴とするインテリアＣＴの画像生成方法。
　請求項１のインテリアＣＴの画像生成方法において、前記ＲＯＩの再構成処理工程は、解析的画像再構成法、逐次近似画像再構成法、統計的画像再構成法の何れか、又は、それらの組み合わせにより実行されることを特徴とするインテリアＣＴの画像生成方法。
　請求項１のインテリアＣＴの画像生成方法において、アクティブコリメータを利用して、前記インテリアＣＴ投影データと前記部分的な全体投影データを取得する量子ビームの開口角度を制御することを特徴とするインテリアＣＴの画像生成方法。
　請求項１のインテリアＣＴの画像生成方法において、前記部分的な全体投影データとして、位置合わせのためのScout-Viewスキャン投影データを利用することを特徴とするインテリアＣＴの画像生成方法。
　請求項１のインテリアＣＴの画像生成方法において、前記部分的な全体投影データを、前記インテリアＣＴ投影データよりも低い分解能の他の量子ビームを照射する測定を利用して取得することを特徴とするインテリアＣＴの画像生成方法。
　請求項１のインテリアＣＴの画像生成方法において、量子ビームを対象物に照射して測定される前記インテリアＣＴ投影データ及び前記全体投影データが、量子ビームと物体の相互作用により発生する、吸収、位相シフト、散乱、回折、屈折の少なくとも何れか一つの空間的な物理量分布の、前記量子ビームが通過する直線上の線積分になっていることを特徴とするインテリアＣＴの画像生成方法。
　インテリアＣＴの画像生成方法であって、
　ＣＴ測定の幾何学系において撮像対象の内部の関心領域（ＲＯＩ）を通過する全ての量子ビームのみを測定してトランケーションを含んだインテリアＣＴ投影データを取得する工程と、
　前記ＣＴ測定の幾何学系において前記撮像対象の外部にあるセグメントからの前記撮像対象の全体を通過する量子ビームを測定してトランケーションを含まない部分的な全体投影データを取得する工程を含んでおり、更に、
　取得したトランケーションを含んだ前記インテリアＣＴ投影データとトランケーションを含まない前記部分的な全体投影データに基づいて前記ＲＯＩを厳密に再構成する処理工程を含んでいることを特徴とするインテリアＣＴの画像生成方法。
　請求項１２のインテリアＣＴの画像生成方法において、前記ＲＯＩの再構成処理工程は、解析的画像再構成法、逐次近似画像再構成法、統計的画像再構成法の何れか、又は、それらの組み合わせにより実行されることを特徴とするインテリアＣＴの画像生成方法。
　請求項１２のインテリアＣＴの画像生成方法において、前記トランケーションを含んだ前記インテリアＣＴ投影データにおけるヒルベルト変換の計算不可能性を、前記トランケーションを含まない前記部分的な全体投影データを用いて計算可能にすることにより、前記ＲＯＩを再構成することを特徴とするインテリアＣＴの画像生成方法。