(57)【特許請求の範囲】
【請求項１】
放射線源から放出され被検体を透過した放射線を検出器アレイで検出して放射線強度分布を求める放射線計測方法であって、
応答関数評価装置が、前記検出器アレイ内のクロストークに関する応答関数を求め、
クロストーク補正装置が、規格化された応答関数および補正行列を用いて、クロストークが含まれた検出器アレイの信号出力から、ベイズの定理に基づく推定式を用いて検出器アレイ入射前のクロストークが含まれない放射線強度分布を求めることを特徴とする放射線計測方法。
【請求項２】
前記応答関数評価装置は、シミュレーション計算または実測によって前記応答関数を求めることを特徴とする請求項１に記載の放射線計測方法。
【請求項３】
放射線を発生する放射線源装置と、
被検体を透過した前記放射線を検出して検出信号を出力する検出器アレイと、
前記検出器アレイ内のクロストークに関する応答関数を求める応答関数評価装置と、
規格化された応答関数および補正行列を用いて、前記検出信号から、ベイズの定理に基づく推定式を用いてクロストークの影響を補正した放射線強度分布を求めるクロストーク補正装置を備えることを特徴とする放射線計測装置。
【請求項４】
前記応答関数評価装置は、シミュレーション計算または実測によって前記応答関数を求めることを特徴とする請求項３に記載の放射線計測装置。
【請求項５】
Ｘ線を発生するＸ線源と、
被検体を透過した前記Ｘ線を検出して検出信号を出力する検出器アレイと、
前記検出器アレイ内のクロストークに関する応答関数を求める応答関数評価装置と、
規格化された応答関数および補正行列を用いて、前記検出信号から、ベイズの定理に基づく推定式を用いてクロストークの影響を補正したＸ線強度分布を求めるクロストーク補正装置と、
前記Ｘ線強度分布を用いて前記被検体の透過像を作成する透視像作成装置を備えることを特徴とするＸ線透過像撮影装置。
【請求項６】
前記応答関数評価装置は、シミュレーション計算または実測によって前記応答関数を求めることを特徴とする請求項５に記載のＸ線透過像撮影装置。
【請求項７】
Ｘ線を発生するＸ線源と、
被検体を透過した前記Ｘ線を検出して検出信号を出力する検出器アレイと、
前記検出器アレイ内のクロストークに関する応答関数を求める応答関数評価装置と、
規格化された応答関数および補正行列を用いて、前記検出信号から、ベイズの定理に基づく推定式を用いてクロストークの影響を補正したＸ線強度分布を求めるクロストーク補正装置と、
前記Ｘ線強度分布を用いて前記被検体のＣＴ画像を作成するＣＴ画像再構成装置を備えることを特徴とするＸ線ＣＴ装置。
【請求項８】
前記応答関数評価装置は、シミュレーション計算または実測によって前記応答関数を求めることを特徴とする請求項７に記載のＸ線ＣＴ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の検出器チャネルにより構成される検出器アレイを用いたＸ線ＣＴ装置に関するものであり、特に、検出器チャネル間のクロストークを補正するために好適な放射線計測方法、放射線計測装置、Ｘ線透過像撮影装置及びＸ線ＣＴ装置に関するものである。
【背景技術】
【０００２】
Ｘ線ＣＴ装置などの放射線計測システムにおいて、撮影画像の空間分解能を向上しようとした場合には、検出器チャネルを稠密化して放射線の角度分布を詳細なピッチで取得することが考えられる。検出器チャネル間を稠密化するためには検出器チャネル間に従来配置していた遮蔽材の厚みを減じたり、無くしたりすることが不可欠となるが、このような場合において放射線がクロストークして真の信号以外の検出信号を周囲の検出器チャネルが計測してしまい、かえって分解能が悪化してしまうというトレードオフの課題があった。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第４１９４１２８号
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許第４１９４１２８号で提案された方法は、クロストークを表現した補正行列を作成し、近似行列を作成して補正を実行しようと試みている。しかしながら、一般的に行列演算のみで逆問題を解決することは困難であって近似的な解放を用いる必要があり、得られた解に恣意性が含まれてしまう。この恣意性は特に産業用のＸ線ＣＴ装置などで、ＣＴ画像を元にした形状計測などの解析を行う場合にトレーサビリティが欠如してしまうという課題も生じる、
本発明の目的は、稠密化した検出器チャネル間で顕在化するクロストークに対して、科学的根拠に基づいてチャネル間クロストークを補正できる放射線計測方法、放射線計測装置、Ｘ線透過像撮影装置及びＸ線ＣＴ装置を提供することにある。
【課題を解決するための手段】
【０００５】
上記目的を達成する本発明の特徴は、放射線源から放出され被検体を透過した放射線を検出器アレイで検出して放射線強度分布を求める放射線計測方法であって、応答関数評価装置が検出器アレイ内のクロストークに関する応答関数を求め、クロストーク補正装置が規格化された応答関数および補正行列を用いて、クロストークが含まれた検出器アレイの信号出力から、ベイズの定理に基づく推定式を用いて検出器アレイ入射前のクロストークが含まれない放射線強度分布を求めることにある。
【０００６】
また、上記目的を達成する他の特徴は、放射線を発生する放射線源装置と、被検体を透過した放射線を検出して検出信号を出力する検出器アレイと、検出器アレイ内のクロストークに関する応答関数を求める応答関数評価装置と、規格化された応答関数および補正行列を用いて、検出信号から、ベイズの定理に基づく推定式を用いてクロストークの影響を補正した放射線強度分布を求めるクロストーク補正装置を備える放射線計測装置にある。
【０００７】
また、上記目的を達成するさらに他の特徴は、Ｘ線を発生するＸ線源と、被検体を透過したＸ線を検出して検出信号を出力する検出器アレイと、検出器アレイ内のクロストークに関する応答関数を求める応答関数評価装置と、規格化された応答関数および補正行列を用いて、検出信号から、ベイズの定理に基づく推定式を用いてクロストークの影響を補正したＸ線強度分布を求めるクロストーク補正装置と、Ｘ線強度分布を用いて被検体の透過像を作成する透視像作成装置を備えるＸ線透過像撮影装置にある。
【０００８】
また、上記目的を達成するさらに他の特徴は、Ｘ線を発生するＸ線源と、被検体を透過したＸ線を検出して検出信号を出力する検出器アレイと、検出器アレイ内のクロストークに関する応答関数を求める応答関数評価装置と、規格化された応答関数および補正行列を用いて、検出信号から、ベイズの定理に基づく推定式を用いてクロストークの影響を補正したＸ線強度分布を求めるクロストーク補正装置と、Ｘ線強度分布を用いて被検体のＣＴ画像を作成するＣＴ画像再構成装置を備えるＸ線ＣＴ装置にある。
【発明の効果】
【０００９】
本発明によれば、検出器チャネル間のクロストークを補正することができため、空間分解能がより向上した放射線計測方法及び放射線計測装置を実現することができるようになる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施の一形態に係る放射線計測方法における計測の流れを示す図である。（実施例１）
【図２】本発明の実施の一形態に係る放射線計測方法の基本アルゴリズムを示す図である。
【図３】本発明の実施の一形態に係る放射線計測方法において応答関数をシミュレーション計算により評価する場合の構成を示す図である。
【図４】本発明の第２の実施例に係る放射線計測方法におい応答関数を実際の測定により評価する場合の構成を示す図である。
【図５】本発明の第３の実施例に係る放射線計測装置の基本構成を示す図である。
【図６】本発明の第４の実施例に係るＸ線ＣＴ装置の基本構成を示す図である。
【発明を実施するための形態】
【００１１】
本発明の実施例を以下に説明する。
【００１２】
＜実施例１＞
本発明の好適な実施例である実施例１の放射線計測方法を、図１を用いて以下に説明する。
【００１３】
図１は、本発明の実施の一形態に関わる放射線計測装置の主要構成を示す図である。本実施例の放射線計測装置は、Ｘ線（例えばＸ線ファンビーム）を放出するＸ線源装置１０と、被検体４０を透過したＸ線を検出する複数のＸ線検出器（検出器チャンネル）を有する検出器アレイ２１と、各々の検出器チャンネルからの出力信号を測定する多チャンネル測定回路３１と、応答関数評価装置９０と、検出器アレイ２１からの出力信号に基づいてクロストークを考慮した放射線強度分布を求めるクロストーク補正装置９２とを備える。図１には、検出器アレイ２１に入射する前のＸ線強度分布５１と、検出器アレイ２１からの信号出力から求めたＸ線強度分布５２を示す。
【００１４】
Ｘ線源装置１０から放出されたＸ線ファンビーム１１が被検体４０を透過して透過経路に応じたＸ線強度分布５１を形成する。検出器アレイ２１に入射したＸ線の強度分布に応じた検出器アレイの信号出力５２が得られる。しかしながら、検出器アレイ２１に入射したＸ線の一部は散乱等の相互作用により、検出器チャネル間でクロストークＸ線１２となり、隣接する検出器チャネルに回り込んで、検出器チャンネルはクロストークＸ線１２による信号出力５３を生成する。
【００１５】
さて、検出器アレイの検出器チャネルｄｎに入射するＸ線強度φ（Ｘｍ）によって、信号Ｓ（ｄｎ）だけでなく周囲の検出器チャネルｄｎ−１，  ｄｎ−２，  ・・・，  ｄ１ 、およびｄｎ＋１，  ｄｎ＋２，・・・，  ｄＮにクロストーク信号出力を付与する。つまり、検出器アレイ入射前のＸ線強度分布φ（ｘｊ）（ｊ＝１〜Ｍ）がＳ（ｄｎ）（ｉ＝１〜Ｎ）の信号出力で観測される。この過程は、検出器アレイ中の唯一の検出器チャネルｄｉにＸ線φ（ｘｊ）が入射したときにもたらされる信号強度ｒｉｊを、Ｎ行Ｍ列の構成要素とした応答関数行列Ｒで表現することができる。なお、応答関数評価装置９０が、検出器アレイの幾何体系と入射Ｘ線源のエネルギースペクトルを元にして計算機シミュレーションによって応答関数Ｒを算出する。
【００１６】
図３に、応答関数評価装置９０で行う計算機シミュレーションの幾何体系を示している。Ｘ線強度分布６１を持つ仮想的なペンシル状のＸ線ビーム１３を作成し、特定の位置Ｘｊに入力する。このとき得られる検出器アレイの信号出力６２が応答関数行列ＲのＭ列の成分となる。応答関数評価装置９０は、このシミュレーション計算を位置Ｘｊ（ｊ＝１〜Ｍ）まで行い、各列成分を組み合わせることで応答関数行列Ｒが得られる。
【００１７】
前記までの各関数を用いて、クロストーク補正装置９２において、図１の構成で得られる検出器アレイの信号出力Ｓ（ｄｉ）は次式（１）の関係で表すことができる。
【００１８】
Ｓ（ｄｉ）＝Ｒ・φ（ｘｊ）      （１）
しかしながら、条件付き確率に関するベイズの定理に基づいた補正手法を適用するため、応答関数行列Ｒの列方向の和が１となるように規格化された規格化応答関数行列Ｒ’および補正行列Ｂを用いて次式（２）に示すように演算を行う。
【００１９】
Ｒ＝Ｒ’・Ｂ                    （２）
このとき補正行列ＢはＭ行Ｍ列の正方行列で、その行列成分ｂｉｊはディラックのデルタ関数を用いて次式（３）となる。
【００２０】
ｂｉｊ＝δｉｊ×εｊ            （３）
成分εｊは次式（４）となる。
【００２１】
εｊ＝Σｒｉｊ  （ｉ＝１〜Ｎ）  （４）
図２は詳細なデータフローと計算式を示す図である。
【００２２】
図２において、応答関数評価装置９０は応答関数Ｒを出力可能であり、応答関数評価装置９０から出力された応答関数Ｒは演算部９４に供給される。演算部９４には応答関数評価装置９０からの応答関数Ｒが供給され、規格化応答関数Ｒ’と補正行列Ｂを演算し、ベイズ推定式演算部９５および除算部９６に供給する。また多チャンネル測定回路３１からの信号出力５２もベイズ推定式演算部９５に供給される。
【００２３】
ベイズ推定式演算部９５は、検出器アレイの信号出力Ｓ（ｄｊ）と、規格化応答関数Ｒ’とを用いて、次に示すベイズの定理に基づいた推定式（５）を演算し、クロストークが補正されたＸ線強度分布ｅｓｔ（ｘｊ）（ｌ＋１）を得る。
【００２４】
【数１】
0006560990.tif 000002
【００２５】
ベイズ推定式演算部９５が演算したｅｓｔ（ｘｊ）（ｌ＋１）は除算部９６に供給される。
【００２６】
除算部９６には、補正行列Ｂも供給されており、除算部９６は補正行列の成分ε’ｊを用いて次式（６）に示すように、ｅｓｔ（ｘｊ）（ｌ＋１）をε’ｊで除算して、検出器アレイ入射前のＸ線強度分布φ（ｘｊ）（ｌ＋１）を算出する。
【００２７】
【数２】
0006560990.tif 000003
【００２８】
なお、演算部９４とベイズ推定式演算部９５と除算部９６はクロストーク補正装置９２に備えられている。
【００２９】
つまり、応答関数評価装置９０は、ベイズの定理に基づいた推定式（５）を用いて、複数回の繰り返し推定計算を用い、その後ε’ｊで除算することによって、検出器アレイ内でクロストークを生じる直前の検出器アレイ入射前のＸ線強度分布φ（ｘｊ）（ｊ＝１〜Ｍ）を算出することができる。
【００３０】
以上のように、本実施例では、クロストーク補正装置９２は、検出器間のクロストークを表現した応答関数を規格化し、その規格化した応答関数および補正行列を算出する。算出した規格化応答関数と多チャンネル測定回路３１から供給される検出器アレイ２１からの検出信号５２とを用いてベイズの定理に基づいた式で繰り返し推定した演算結果を得る。クロストーク補正装置９２は補正行列の対角成分で上記演算結果を除算してクロストークが補正された放射線強度分布を求めている。
【００３１】
本実施例によれば、予めシミュレーション計算によって求めたクロストークを表現する応答関数行列を作成しておき、実際のクロストークした検出器チャネルの信号出力に対して、条件付き確率に関するベイズの定理に基づく推定計算を行うことで、クロストークする前の放射線の角度分布を詳細なピッチで得ることができるようになる。
【００３２】
また、本実施例によれば、規格化された応答関数および補正行列を用いて、クロストークが含まれた検出器アレイからの検出信号から、ベイズの定理に基づく推定式を用いて検出器アレイ入射前のクロストークが含まれない放射線強度分布を求めることができる。つまり、検出器チャンネル間で生じるクロストークによる影響を補正することができるため、空間分解能が向上した放射線計測方法及び放射線計測装置を実現できるようになる。このようなクロストークは検出器チャネル間の遮蔽材を設置しない検出器アレイや、薄い遮蔽材のみを設置した検出器アレイを用いた場合に特に生じやすいため、このような検出器アレイを用いる場合に特に有効である。
【００３３】
＜実施例２＞
第２の実施例を、図４を用いて説明する。第１の実施例では、応答関数評価装置９０がシミュレーション計算で応答関数行列Ｒを求める例を説明したが、本実施例では実測値を用いて応用関数行列Ｒを求める。本実施例について、第１の実施例と異なる構成を中心に、以下に説明する。
【００３４】
図４に示すように、Ｘ線源装置１０と検出器アレイ２１の間に、１チャンネルコリメータ４１を設置する。１チャンネルコリメータ４１は、Ｘ線源装置１０から放出されたＸ線がいずれか一つの検出器チャンネルに入射されるように孔部（貫通穴）をもつ。１チャネルコリメータ４１を検出器アレイ２１の前面に配置することで、Ｘ線源装置１０から放出されたＸ線ファンビーム１１が検出器チャネルｄｎ以外に入射しないように構成され、検出器アレイ２１に入射する前のＸ線強度分布が図４に示すような、ペンシル状のＸ線強度分布６１となる。本実施例では、検出器アレイの信号出力６２を検出器アレイ３２で実測することにより応答関数行列Ｒが得られる。得られた応答関数行列は第１の実施例と同様にベイズの定理に基づいた推定計算に用いる。また、実施例１と同様の方法により、検出器アレイ内でクロストークを生じる直前の検出器アレイ入射前のＸ線強度分布φ（ｘｊ）（ｊ＝１〜Ｍ）を算出する。
【００３５】
本実施例によれば、検出器アレイ２１を使った実測値を用いてクロストークを表現する応答関数行列を作成し、実際のクロストークした検出器チャネルの信号出力に対して、条件付き確率に関するベイズの定理に基づく推定計算を行うことで、クロストークする前の放射線の角度分布を詳細なピッチで得ることができるようになる。
【００３６】
本実施例によれば、規格化された応答関数および補正行列を用いて、クロストークが含まれた検出器アレイからの検出信号から、ベイズの定理に基づく推定式を用いて検出器アレイ入射前のクロストークが含まれない放射線強度分布を求めることができる。つまり、検出器チャンネル間で生じるクロストークによる影響を補正することができるため、空間分解能が向上した放射線計測方法及び放射線計測装置を実現できるようになる。このようなクロストークは検出器チャネル間の遮蔽材を設置しない検出器アレイや、薄い遮蔽材のみを設置した検出器アレイを用いる場合に特に生じやすいため、このような検出器アレイを用いる場合に特に有効である。
【００３７】
＜実施例３＞
第３の実施例を、図５を用いて説明する。本実施例では、Ｘ線源装置１０と検出器アレイ２１の間に被検体４０を配置して、被検体４０を透過してくるＸ線強度信号を検出器アレイ２１で測定し、Ｘ線透過像を撮影する透過像撮影装置について説明する。
【００３８】
本実施例の透過像撮影装置は、Ｘ線（例えばＸ線ファンビーム）を放出するＸ線源装置１０と、被検体４０を透過したエックス線を検出する複数のＸ線検出器（検出器チャンネル）を有する検出器アレイ２１と、各々の検出器チャンネルからの出力信号を測定する多チャンネル測定回路３１と、応答関数評価装置９０と、検出器アレイ２１からの出力信号に基づいてクロストークを考慮した放射線強度分布を求めるクロストーク補正装置９２と、クロストーク補正装置９２からの放射線強度分布の情報に基づいて透視画像を生成する透過像作成装置８０を備える。
【００３９】
本実施例の透視画像撮影装置は、検出器アレイ２１内における検出器チャネル間のクロストークＸ線１２による信号が多チャンネル測定回路３１から出力された段階で、応答関数評価装置９０およびクロストーク補正装置９２を用いて補正する。クロストーク補正装置９２からクロストークが補正されたＸ線強度分布を透過像作成装置８０に供給することにより空間分解能が向上した透過像が得られるようになる。本実施例の応答関数評価装置９０は、実施例１又は実施例２のいずれかの方法で応答関数行列Ｒを求めればよい。
【００４０】
本実施例によれば、クロストークを表現する応答関数行列を作成し、実際のクロストークした検出器チャネルの信号出力に対して、条件付き確率に関するベイズの定理に基づく推定計算を行うことで、クロストークする前の放射線の角度分布を詳細なピッチで得ることができ、より空間分解能が向上した透過像を得ることができるようになる。
【００４１】
本実施例によれば、規格化された応答関数および補正行列を用いて、クロストークが含まれた検出器アレイからの検出信号から、ベイズの定理に基づく推定式を用いて検出器アレイ入射前のクロストークが含まれない放射線強度分布を求めることができる。つまり、検出器チャンネル間で生じるクロストークによる影響を補正することができるため、空間分解能が向上した透過像撮影方法及び透過像撮影装置を実現できるようになる。このようなクロストークは検出器チャネル間の遮蔽材を設置しない検出器アレイや、薄い遮蔽材のみを設置した検出器アレイを用いる場合に特に生じやすいため、このような検出器アレイを用いる場合に特に有効である。
【００４２】
＜実施例４＞
第４の実施例を、図６を用いて説明する。本実施例では、Ｘ線ＣＴ装置のスキャナテーブル６０上に被検体４０を搭載し、被検体４０を透過してくるＸ線の強度信号を検出器アレイ２１で測定し、Ｘ線ＣＴ画像を撮影するＸ線ＣＴ装置について説明する。
【００４３】
本実施例のＸ線ＣＴ装置は、Ｘ線（例えばＸ線ファンビーム）を放出するＸ線源装置１０と、被検体４０をのせるためのスキャナテーブル６０と、被検体４０を透過したエックス線を検出する複数のＸ線検出器（検出器チャンネル）を有する検出器アレイ２１と、各々の検出器チャンネルからの出力信号を測定する多チャンネル測定回路３１と、応答関数評価装置９０と、検出器アレイ２１からの出力信号に基づいてクロストークを考慮した放射線強度分布を求めるクロストーク補正装置９２と、クロストーク補正装置９２からの放射線強度分布の情報に基づいてＣＴ画像を生成するＣＴ画像再構成演算装置８１を備える。
【００４４】
本実施例のＸ線ＣＴ装置は、検出器アレイ２１内における検出器チャネル間のクロストークＸ線１２による信号が多チャンネル測定回路３１から出力された段階で、応答関数評価装置９０およびクロストーク補正装置９２を用いて補正する。クロストーク補正装置９２からクロストークが補正されたＸ線強度分布をＣＴ画像再構成演算装置８１に供給することにより空間分解能が向上したＣＴ画像が得られるようになる。本実施例の応答関数評価装置９０は、実施例１又は実施例２のいずれかの方法で応答関数行列Ｒを求めればよい。
【００４５】
本実施例によれば、クロストークを表現する応答関数行列を作成し、実際のクロストークした検出器チャネルの信号出力に対して、条件付き確率に関するベイズの定理に基づく推定計算を行うことで、クロストークする前の放射線の角度分布を詳細なピッチで得ることができ、より空間分解能が向上したＸ線ＣＴ画像を得ることができるようになる。
【００４６】
本実施例によれば、規格化された応答関数および補正行列を用いて、クロストークが含まれた検出器アレイからの検出信号から、ベイズの定理に基づく推定式を用いて検出器アレイ入射前のクロストークが含まれない放射線強度分布を求めることができる。つまり、検出器チャンネル間で生じるクロストークによる影響を補正することができるため、空間分解能が向上したＸ線ＣＴ画像撮影方法及びＸ線ＣＴ画像装置を実現できるようになる。このようなクロストークは検出器チャネル間の遮蔽材を設置しない検出器アレイや、薄い遮蔽材のみを設置した検出器アレイを用いる場合に特に生じやすいため、このような検出器アレイを用いる場合に特に有効である。
【符号の説明】
【００４７】
１０  Ｘ線源装置
１１  Ｘ線ファンビーム
１２  検出器チャネル間のクロストークＸ線
１３  ペンシル状のＸ線ビーム
２１  検出器アレイ
３１  多チャンネル測定回路
４０  被検体
４１  １チャネルコリメータ
５１  検出器アレイ入射前のＸ線強度分布
５２  検出器アレイの信号出力
５３  クロストークＸ線による信号出力
５４  クロローク補正された検出器アレイの信号出力
６０  スキャナテーブル
６１  ペンシル状のＸ線の強度分布
６２  ペンシル状のＸ線の強度分布による検出アレイの信号出力
８０  透過像作成装置
８１  ＣＴ画像再構成演算装置
９０  応答関数評価装置
９２  クロストーク補正装置
９４  演算部
９５  ベイズ推定式演算部
９６  除算部
【図１】
0006560990.tif 000004
【図２】
0006560990.tif 000005
【図３】
0006560990.tif 000006
【図４】
0006560990.tif 000007
【図５】
0006560990.tif 000008
【図６】
0006560990.tif 000009