(57)【特許請求の範囲】
【請求項１】
物体を透過した電離性電磁放射線のスペクトルを修正するための方法であって、
物体は、放射線源と検出器との間に配置され、放射線源は、前記物体に向けて、入射放射線と呼ばれる電離性電磁放射線を放出するために設けられ、
検出器は複数の画素を備え、各画素は、物体を透過して検出器に向かう放射線を検出し、そのスペクトルを獲得するために設けられ、透過放射線は、物体中での入射放射線の散乱より生じた散乱放射線と、１次放射線とを有し、
方法は、次のステップを有する。
ａ）線源と物体との間にマスクを挿入し、いくつかの画素により、物体を透過した第１放射線の第１スペクトルを獲得するステップ。マスクは、前記入射放射線の一部分を減衰することができる減衰要素を備え、マスクの検出器上への投影は、第１グループの画素の境界を定める。
ｂ）第１グループの画素の各画素について、前記マスクなしの場合に、物体を透過して検出器に到達する第２放射線の第２スペクトルを取得するステップ。
ｃ）前記第１グループの画素の各画素について、比較スペクトルを取得するために、第１スペクトルと第２スペクトルとを比較するステップ。
ｄ）第１グループの画素の各画素において、あらかじめ得た遷移行列を前記比較スペクトルに適用し、物体を透過して前記画素に到達した１次放射線を示す、いわゆる１次スペクトルの評価値を取得し、物体を透過した散乱スペクトルを示す散乱スペクトルを評価するステップ。
ｅ）検出器の全て又はいくつかの画素について、修正スペクトルを取得するために、ステップｄ）で評価した各散乱スペクトルから、第２スペクトル又は第１スペクトルを修正するステップ。
【請求項２】
ステップｅ）は、修正の前に、検出器の全ての画素について、散乱スペクトルを評価することを有する
請求項１記載の方法。
【請求項３】
ステップｅ）は、各画素について、第２スペクトル又は第１スペクトルから、評価した散乱スペクトルを減算することを有する
請求項１又は請求項２記載の方法。
【請求項４】
ステップｂ）において、第１グループの各画素について、
前記画素により第２スペクトルを獲得することによって、又は、
前記第１グループに属さない少なくとも１つの画素により獲得された、少なくとも１つの第１スペクトルに基づく評価によって、
第２スペクトルを取得する
請求項１乃至３のいずれか１項記載の方法。
【請求項５】
ステップａ）は複数の配置で実行され、各配置において、各画素について、第１スペクトルを取得するために、各配置は物体に対する検出器及び線源の位置と関連付けられ、各配置は更に第１グループの画素と関連付けられ、
ステップｂ）は、配置と関連付けられた第１グループの全て又はいくつかの画素について、他の配置により取得された第１スペクトルから第２スペクトルを決定することを有する
請求項１乃至４のいずれか１項記載の方法。
【請求項６】
各配置はパラメータと関連付けられ、第１パラメータと関連付けられた配置による少なくとも１つの第２スペクトルは、第１パラメータとは異なる第２パラメータと関連付けられた他の配置により取得された第１スペクトルから取得される
請求項５記載の方法。
【請求項７】
前記パラメータが、物体に対する線源及び／又は検出器の傾斜角である
請求項６記載の方法。
【請求項８】
ステップｄ）において、１次スペクトルの評価値は、各比較スペクトルによる前記遷移行列の行列積を有する
請求項１乃至７のいずれか１項記載の方法。
【請求項９】
各減衰要素は、それが曝される放射線の５〜８０％を減衰するように設けられる
請求項１乃至８のいずれか１項記載の方法。
【請求項１０】
マスクは表面に沿って延在し、各減衰要素は、１ｃｍよりも小さな距離により、他の減衰要素から離れている
請求項１乃至９のいずれか１項記載の方法。
【請求項１１】
ステップａ）〜ｅ）は複数の配置により実行され、各配置において、複数の画素について、修正スペクトルを取得するために、各配置は、物体に対する検出器及び線源の位置と関連付けられ、各配置の修正スペクトルは、物体の断層再構成を生成するために用いられる
請求項１乃至１０のいずれか１項記載の方法。
【請求項１２】
ステップａ）〜ｅ）は、次のステップにより補われる、
ｆ）少なくとも１つのエネルギー又はエネルギー範囲を選択するステップ。
ｇ）画像の各画素が、選択されたエネルギー内、又は、選択されたエネルギー範囲内の、検出器の画素と関連付けられた、修正スペクトルから取得したデータを有する、画像を生成するステップ。
請求項１乃至１０のいずれか１項記載の方法。
【請求項１３】
ステップｇ）において、画像の各画素は、選択されたエネルギー範囲内の前記修正スペクトルの積分値又は平均値と関連する情報項目を有する
請求項１１記載の方法。
【請求項１４】
マスクが、放射線源と物体との間に挿入される
請求項１乃至１３のいずれか１項記載の方法。
【請求項１５】
遷移行列は、複数のいわゆる校正計測を実行することにより取得され、各校正計測は、放射線源と検出器との間に既知の性質及び厚さを有する材料を挿入することにより実行される
請求項１乃至１４のいずれか１項記載の方法。
【請求項１６】
請求項１乃至１４のいずれか１項記載の方法のステップｃ）〜ｅ）を実行するための、マイクロプロセッサにより実行可能な指令を有する
情報記録媒体。
【請求項１７】
入射放射線と呼ばれる電離性電磁放射線を物体に向けて放出することができる放射線源と、
物体を透過して検出器に向かう放射線を検出し、そのスペクトルを獲得することができる画素を備える検出器と、
線源と物体との間に挿入されることができ、検出器上への投影が第１グループの画素の境界を定める、入射放射線の一部を減衰することができる減衰要素を備える、マスクと、
各画素により検出されたスペクトルを受け、請求項１乃至１４のいずれか１項記載の方法のステップｃ）〜ｅ）を実行することができるプロセッサと
を備える、物体の画像を生成するための装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の技術分野は、電離放射線、特に、Ｘ線の線源により照射された物体を透過した放射線のスペクトルの処理である。本発明の目的は、ピクセル型検出器（pixelated detector）を用いて複数のスペクトルを取得すること、及び、散乱放射線を示す成分の影響を限定するために、スペクトルの全て又はいくつかの修正を行うことである。その応用例は、医用画像又は非破壊検査での利用である。
【背景技術】
【０００２】
医療又は工業の分野でのＸ線による物体の検査が、非常に広範に利用されている。既存の方法は、放射線源と検出器との間に物体を配置して、線源を用いて物体を照射することにある。そして、検出器は、物体を透過した放射線の、一般に２次元の、画像を形成する。この画像は、線源により放出された放射線の、物体による減衰を示す。
【０００３】
物体を透過した放射線は、線源により放出された放射線の、物体による散乱の結果の成分を有する。そのことは、放射線のエネルギーが弱いとき、及び／又は、物体が大きな原子番号の材料からなるときに、より一層顕著になる。一般に散乱放射線と呼ばれるこの成分は、それだけが物体による減衰に間接的につながることから、画像の解釈を混乱させる。その上、１次放射線と呼ばれる非散乱放射線は、線源と検出器との間を直線軌道に沿って伝播するのに対して、散乱された放射線は、物体のどのような箇所からも発生し、発生のこの箇所からのその軌道は、異なる角度により分散される。それゆえ、検出器により計測された信号から散乱成分を抽出するために、画像を解釈することを目的として画像を処理する前に、この散乱成分を評価することが求められている。
【０００４】
線源と物体との間を直線軌道で伝播する、１次放射線と呼ばれる非散乱放射線を本質的に示す画像を得るために、これらの画像での散乱放射線の割合を評価し、削減することを試みる多くの方法が開発されてきた。
【０００５】
例えば、非特許文献１は、Ｘ線の線源と物体との間に、格子の形態を採る除去可能なマスクを挿入することを有する方法を記載する。非特許文献１は、そのようなマスクが、高い空間周波数において、１次放射線の顕著な空間変調を生じるのに対して、低い空間周波数において、散乱放射線へのその影響があまり顕著ではないことに基づいている。このように、周波数空間において、１次放射線と散乱放射線との間の区別を得ることは可能である。
【０００６】
その上、非特許文献２も、Ｘ線の線源と物体との間に除外可能なマスクを挿入することを有する方法を記載する。マスクは減衰を可能にし、マスクが線源と物体との間に配置されたときに、検出器により計測される放射線は、散乱放射線だけを示すと考えられる。言い換えると、このマスクは、全体として考えると、１次放射線の減衰を可能にする。
【０００７】
現在の検出器は、決して進歩を止めない品質を有する２次元の画像を得ることを可能にする。１つの最近の傾向は、スペクトル画像、すなわち、異なるエネルギー帯による画像の獲得を可能にする検出器の出現である。分光測定機能を有する半導体検出器にしばしば基づくこれらの検出器は、一般に２次元で得られる、獲得したデータにスペクトル次元を追加する。そのとき、物体を透過した、異なるエネルギー帯による放射線の画像を得ることが可能である。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】Zhu L, "Scatter Correction Method for X-ray CT Using Primary Modulation: Theory and Preliminary Results", IEEE Transactions on Medical Imaging, （米）, Dec 2006, Vol. 25, No. 12
【非特許文献２】Ning R, "X-ray Scatter Correction Algorithm for Cone Beam CT-Imaging", Med.Phys., （米）, May 2004, 31 (5)
【発明の概要】
【発明が解決しようとする課題】
【０００９】
現在、上記方法は、スペクトル画像解析（spectral imaging）とは関係しない。本発明者は、検出器により検出された放射線から１次放射線を抽出し、しかも、それをいくつかのエネルギー帯において同時に行うために、検出器により検出された放射線の効果的な修正を可能にする、スペクトル画像解析に基づく、代わりの解決方法を提案する。明細書の最後に記載した試験は、方法の効果を証明する。
【課題を解決するための手段】
【００１０】
本発明は、添付した請求項のいずれか１つに係る、物体を透過した電離性電磁放射線（ionizing electromagnetic radiation）のスペクトルを修正する方法である。より正確には、本発明の目的は、物体を透過した電離性電磁放射線のスペクトルを修正するための方法であり、
物体は、放射線源と検出器との間に配置され、放射線源は、上記物体に向けて、入射放射線と呼ばれる電離性電磁放射線を放出することが可能であり、
検出器は複数の画素を備え、各画素は物体を透過して検出器に向かう放射線を検出して、それからスペクトルを獲得することができ、透過放射線は、物体中での入射放射線の散乱により生じた散乱放射線と、１次放射線とを有し、
方法は、次のステップを有する。
ａ）線源と物体との間にマスクを挿入し、いくつかの画素により、物体を透過した第１放射線の第１スペクトルを獲得するステップ。マスクは、上記入射放射線の一部を減衰することができる減衰要素を備え、マスクの検出器上への投影は第１グループの画素の境界を定める。
ｂ）第１グループの画素の各画素について、上記マスクがない場合の、物体を透過して検出器に到達する第２放射線の第２スペクトルを取得するステップ。
ｃ）上記第１グループの画素の各画素について、比較スペクトルを取得するために、第１スペクトルと第２スペクトルとを比較するステップ。
ｄ）第１グループの画素の各画素において、あらかじめ得た遷移行列を上記比較スペクトルに適用し、物体を透過して上記画素に到達する１次放射線を示す、いわゆる１次スペクトルの評価値を得て、物体を透過した散乱放射線を示す、散乱スペクトルを評価するステップ。
ｅ）検出器の全て又はいくつかの画素について、修正スペクトルを得るために、ステップｄ）で評価した各散乱スペクトルから、第２スペクトル又は第１スペクトルを修正するステップ。
【００１１】
マスクは、好ましくは、放射線源と物体との間に挿入される。
【００１２】
実施の形態によると、
ステップａ）は複数の配置において実行され、各配置において及び各画素について第１スペクトルを得るために、各配置は、物体に対する検出器及び線源の位置と関連付けられ、各配置は、第１グループの画素とも関連付けられ、
配置と関連付けられた第１グループの画素の全て又はいくつかについて、ステップｂ）は、他の配置より得られる第１スペクトルから第２スペクトルを決定することを有する。
【００１３】
そのことは、除外可能なマスクを用いることを回避する。そして、画素が、１つの配置により第１グループの画素に属し、上記他の配置と関連付けられた第１グループの画素に属さず、放射線源に直接さらされるという事実が利用される。この他の配置によりそれらの画素が受ける放射線は、マスクがない場合に、上記配置において、それらが受ける放射線を示すものと考えられる。
【００１４】
この実施の形態によると、各配置はパラメータと関連付けられ、第１パラメータと関連付けられた配置による少なくとも１つの第２スペクトルは、第１パラメータとは異なる第２パラメータと関連付けられた他の配置により得られる第１スペクトルから得られる。このパラメータは、物体に対する線源及び／又は検出器の傾斜角であることができる。
【００１５】
断層撮影の（tomographic）実施の形態と呼ばれる１つの実施の形態によると、ステップａ）〜ｅ）は複数の配置により実行され、各配置において、複数の画素に関して修正スペクトルを得るために、各配置は、物体に対する検出器及び線源の位置と関連付けられ、各配置の修正スペクトルは、物体の断層再構成（tomographic reconstruction）を生成するために用いられる。そのことは、本発明を断層再構成に適用することを可能にする。
【００１６】
本発明の他の主題は、この出願で記載される方法を実行するための指令を有する情報記録媒体であり、これらの指令はマイクロプロセッサにより実行されることができる。
【００１７】
本発明の他の主題は、物体の画像を生成するための装置であり、当該装置は、
入射放射線と呼ばれる電離性電磁放射線を上記物体に向けて放出することができる放射線源と、
複数の画素を備え、各画素は、物体を透過して検出器に向かう放射線を検出し、それからスペクトルを獲得することができる、検出器と、
線源と物体との間に挿入されることができ、検出器上への投影が第１グループの画素の境界を定める、入射放射線の一部を減衰することができる減衰要素を備えるマスクと、
各画素により検出されたスペクトルを受け、この明細書に記載された方法のステップｃ）〜ｅ）を実行することができるプロセッサと
を備える。
【発明の効果】
【００１８】
本発明は、次の図面に基づいて、以下に記載される例示の実施の形態からより良く理解されるであろう。
【図面の簡単な説明】
【００１９】
【図１】図１（ａ）、（ｂ）は本発明の実行を可能にする装置を示す。図１（ｃ）は装置が設けられ、更に、図１（ａ）に描かれた、マスクを示す。
【図２】図２（ａ）は検出器応答行列（detector response matrix）を示す。図２（ｂ）、（ｃ）はそれぞれ、この応答行列の行、列を示す。
【図３】本発明の実施の形態の主要なステップを示す。
【図４】図４（ａ）は本発明で実行される遷移行列を示す。図４（ｂ）はこの遷移行列の列を示す。図４（ｃ）はこの遷移行列を得ることを可能にする方法のステップを示す。
【図５】図５（ａ）、（ｂ）は物体の断層再構成に適した本発明の実施の形態を示す。
【図６】図６（ａ）、（ｂ）、（ｃ）は物体の画像をシミュレートするモデル化の結果を示し、各画像は検出器の複数の画素により計測されたスペクトルに基づき、これらのスペクトルはそれぞれ、全放射線、１次放射線、本発明による１次放射線の評価値を示す。図６（ｄ）は図６（ａ）、（ｂ）、（ｃ）のそれぞれの水平プロファイルを示す。
【図７】図７（ａ）、（ｂ）、（ｃ）、（ｄ）は全放射線のスペクトル、１次放射線のスペクトル、本発明により評価された１次放射線のスペクトル、検出器の異なる画素により計測された１次放射線のスペクトルを示す。
【図８】図８（ａ）、（ｂ）、（ｃ）は異なるエネルギー範囲における断層再構成のモデル化を示し、断層再構成はそれぞれ、全放射線、１次放射線、本発明による１次放射線の評価値に基づいて生成される。
【図９】図９（ａ）、（ｂ）、（ｃ）は人体模型で行われた実験計測の画像を示し、各画像は検出器の複数の画素により計測されたスペクトルに基づき、これらのスペクトルはそれぞれ、全放射線、参考方法により評価された１次放射線、本発明による１次放射線の評価値を示す。図９（ｄ）は図９（ｂ）に示す水平線に沿ったプロファイルを示す。
【図１０】図１０（ａ）、（ｂ）、（ｃ）、（ｄ）は全放射線、参考方法により評価された１次放射線、及び、本発明により得られた１次放射線のエネルギースペクトルを示し、これらのスペクトルは、図９（ｂ）に示す４箇所に対応する。
【発明を実施するための形態】
【００２０】
図１（ａ）は、本発明に係る方法を実行する装置１の実施の形態を示す。
【００２１】
放射線源１１は、入射放射線と呼ばれる、電離性電磁放射線Ｉ０を物体１０に向けて放出する。物体１０は、放射線源１１と放射線検出器２０との間に配置される。放射線検出器は、検出平面Ｐと呼ばれる平面上に配置された画素２０ｉ，ｊを備える検出器である。係数ｉ，ｊは、検出平面の各画素の座標を示す。画素は行状に延在することができるが、一般には、それらは通常の２次元マトリックスに従って延在する。
【００２２】
物体１０は、活き活きとした（living）生体組織、例えば、動物又は人間の体の一部であることができる。装置は医用撮像装置である。物体は、工業用部品又は手荷物であることもでき、そのとき、装置は非破壊検査目的で用いられる。この例では、物体は人体模型の一部であって、母体（matrix）１００、２つの肺１０１、１０２、脊柱１０３、椎骨１０４、側面（dimension）１０５、腫瘍１０７、１０８、１０９を備える。
【００２３】
用語「電離性電磁放射線」は、１ｋｅＶよりも高い、そして、好ましくは、５ＭｅＶよりも低いエネルギーの光子からなる電磁放射線を示す。電離放射線のエネルギー範囲は、１ｋｅＶと２ＭｅＶとの間にあることができるが、それはしばしば、１ｋｅＶと、１５０ｋｅＶ又は３００ｋｅＶとの間に拡がる。電離放射線は、Ｘ線又はγ線であることができる。好ましくは、電離放射線源は多エネルギー性（poly-energetic）であって、入射放射線は、一般に数十ｋｅＶ、更には数百ｋｅＶに渡るエネルギー範囲により放出される。それは、特に、Ｘ線放出管である。
【００２４】
入射放射線Ｉ０を形成する光子の一部分は、物体と相互作用をすることなく、物体を通過して、検出器２０に到達する。これらの光子は、１次成分又は一次放射線ＩＰを形成する。入射放射線Ｉ０を形成する他の光子は、例えば、光電効果により、物体中で減衰される。最終的には、いくつかの光子はサンプル中で、コンプトン非弾性散乱タイプ又はレーリー弾性散乱タイプの散乱相互作用を経験する。非弾性、弾性にかかわらず、散乱は光子の向きに変化を生じさせる。
【００２５】
このように、線源１１により照射された物体１０は、透過放射線と呼ばれる放射線Ｉを検出器２０に向けて透過し、放射線Ｉは、
物体と相互作用することがなく、線源からの軌道が直線である、直接成分又は１次放射線ＩＰと、
物体中の入射放射線の散乱に基づく、散乱成分Ｉｄｉｆｆ又は散乱放射線と
を備える。
【００２６】
物体を透過した放射線Ｉは、検出器２０の画素に到達し、各画素はこの放射線の一部を検出する。物体を透過し、画素２０ｉ，ｊにより検出された放射線をＩｉ，ｊで示す。
【００２７】
先行技術と関連付けて説明したように、散乱放射線Ｉｄｉｆｆは計測の解釈を混乱させる。実際に、１次放射線ＩＰと比べて、散乱放射線は物体から検出器に向けて、いろいろな向きで伝播する。このように、検出器の各画素２０ｉ，ｊにより収集された放射線の一部は、放射線源１１から直接生じず、散乱現象に由来する。画像の解釈は、普通、検出器による入射放射線の減衰に基づき、減衰は、所定のエネルギー範囲における、入射放射線Ｉ０の強度に対する１次放射線ＩＰの強度の比により得られる。画像の良い解釈は、１次放射線ＩＰの強度の知識を前提とするが、物体を透過し、検出器により計測される放射線Ｉは、上記１次放射線ＩＰと散乱放射線Ｉｄｉｆｆとの合計値を有する。
【００２８】
各画素２０ｉ，ｊは放射線検出器を構成し、放射線検出器は、
物体を透過する放射線Ｉの光子と相互作用でき、シンチレータ（scintillator）タイプの材料、又は、好ましくは、ＣｄＴｅ、ＣｄＺｎＴｅタイプの、外気温での使用に対応する半導体材料である、検出器材料と、
検出器材料中で相互作用する各光子により預けられるエネルギーに依存し、又は、好ましくは、比例する強度を有する信号を生成することができる電子回路と、
取得期間と呼ばれる期間に検出した信号のエネルギーのＳｉ，ｊで示されるスペクトルを得ることができる、スペクトル測定回路（spectrometry circuit）と
を備える。
【００２９】
このように、画素が規則正しくマトリックス配列で配置されたときに、各画素は物体を透過した放射線のスペクトルＳｉ，ｊを生成することができる。検出器は、各画像が所定のエネルギー範囲ΔＥの各スペクトルの内容を示す、いくつかの画像を形成することができる。典型的には、各画像は、上記エネルギー帯における各スペクトルＳｉ，ｊの積分値又は平均値を有する。検出器が空間的にもスペクトル的にも解像される（resolved）ので、用語「スペクトル画像解析」が適合する。
【００３０】
また、入射放射線Ｉ０による照射の効果の下で、物体１０は透過放射線と呼ばれる放射線Iをピクセル型分光検出器（pixelated spectrometric detector）２０に向けて透過し、その各画素２０ｉ，ｊは、上記透過放射線Ｉを検出し、正しく検出した放射線Ｉｉ，ｊのエネルギースペクトルＳｉ，ｊを形成することができる。
【００３１】
用語「エネルギースペクトル」は、スペクトルの取得期間に検出した信号の強度ヒストグラムに対応する。強度ＡとエネルギーＥとの間の関係は、当業者に知られた原理による
Ｅ ＝ ｇ（Ａ）
のようなエネルギー校正関数（energy calibration function）により得ることができる。それゆえ、エネルギースペクトルＳｉ，ｊはベクトルであって、その各項Ｓｉ，ｊ（ｎ）は、∂Ｅが各チャネルのスペクトル幅であるときの、エネルギー範囲Ｅ±∂Ｅ／２において、画素２０ｉ，ｊにより検出される放射線の量を示す。
【００３２】
各エネルギースペクトルＳｉ，ｊは、雑音項内で（to within a noise term）、ＳＰｉ，ｊで示される１次放射線のスペクトルと、散乱放射線のスペクトルＳdiffｉ，ｊとの合計であると考えることができる。また、
Ｓｉ，ｊ ≒ ＳPｉ，ｊ ＋ Ｓdiffｉ，ｊ  （１）
である。記号「≒」は、雑音項内で等しいことを意味し、この雑音は特に、検出器に由来する、又は、２つの入射光子が同時に検出されたときに起こる、いわゆるスタッキング効果に由来することができる。
【００３３】
本発明の１つの目的は、散乱成分Ｓdiffｉ，ｊを削減し、Ｓ*ｉ，ｊ ≒ ＳPｉ，ｊとなるような修正スペクトルＳ*ｉ，ｊを得るために、各画素で計測されたスペクトルを修正することである。言い換えると、修正スペクトルＳ*ｉ，ｊは、画素２０ｉ，ｊに到達する１次放射線のスペクトルの評価値
【数１】
0006770866.tif 000002
に対応する。
【数２】
0006770866.tif 000003
【００３４】
装置は、更に、線源１１と検出器２０との間、この例では、線源１１と物体１０との間に配置され、好ましい配置を構成するマスク１５を備える。このマスクは、マスクが延在する表面１５Ｓ上に空間的に分散された減衰要素１５Ｘを備える。各減衰要素は、照射源により生成された入射放射線Ｉ０の一部分を部分的に減衰することができる。減衰要素は離散的に分散され、２つの隣り合う減衰要素は、上記減衰要素よりもより少なく減衰する。言い換えると、減衰要素は、減衰係数ａｔｔｘ１５、ａｔｔｘ’１５の離散的な空間分布を定義し、２つの隣り合う減衰要素１５Ｘ、１５Ｘ’の間では、減衰係数ａｔｔ０１５は、各減衰要素と関連する減衰係数ａｔｔｘ１５、ａｔｔｘ’１５よりも小さい。
【００３５】
用語「減衰係数」は当業者に知られている。減衰係数は、Ｉ０（Ｅ）がエネルギーＥで減衰要素１５Ｘ上に入射する入射放射線Ｉ０の強度を示し、ＩＸ（Ｅ）が上記エネルギーＥで減衰要素１５Ｘを透過する入射放射線ＩＸの強度を示す場合に、式
ａｔｔｘ１５（Ｅ） ＝ −ｌｎ［ＩＸ（Ｅ）／Ｉ０（Ｅ）］
で表される。
【００３６】
一般に、線源と検出器との間にマスクを挿入することは、マスクなしの配置と比べて、検出器から生じる散乱放射線を著しく変化させない方が良い。また、好ましくは、各減衰要素は、前に定義したように、減衰係数が、入射放射線Ｉ０が放出されるエネルギー範囲の１つのエネルギーにおいて、又は、このエネルギー範囲の平均エネルギーにおいて、０．０５と１．５との間にあることを示す。このように、散乱を無視することにより、各減衰要素は、好ましくは、線源により生成され、及び／又は、マスクの減衰要素の間に拡がる空間でマスクを通過する入射放射線Ｉ０の５〜８０％を減衰する。好ましくは、減衰係数は１よりも小さく、更には０．５よりも小さく、好ましくは０．３よりも小さい。このように、各減衰要素は、線源により生成された放射線の、又は、マスクの減衰要素の間を通過した放射線の、６０％よりも少ない、又は、４０％よりも少ない、好ましくは、３０％よりも少ない量をそれぞれ減衰する。線源により生成された放射線の５％の減衰に対応する、０．０５に等しい減衰係数以下では、本発明者は減衰は不十分であると考える。言い換えると、マスク１５は、それゆえ、減衰要素１５Ｘと、上記減衰要素の間に拡がる空間との間の減衰の差異（attenuation contrast）を得ることを可能にし、減衰要素は、上記空間を通過する放射線の５〜３０％、更には４０％、更にはそれ以上を減衰する。
【００３７】
加えて、又は、代わりに、マスクにより減衰された入射放射線の割合による充填率（filling factor）の結果の形態で、マスク１５の減衰全体（global attenuation）を定義することが可能であり、充填率は、照射線源１１により放出された入射放射線Ｉ０のエネルギー範囲のエネルギーにおいて、又は、この範囲の平均エネルギーにおいて決定される。充填率は、マスクの全表面領域に対する、減衰要素１５Ｘの組みにより占められるマスクの表面領域の割合に対応する。このように定義したマスクの減衰全体は、好ましくは、１％よりも大きく、１０％よりも小さい。このように、この条件を満足するマスクは０．０８に等しい充填率を有することができ、マスクの各要素１５Ｘは、前に定義したように、０．０８（８％）に等しいマスクの減衰全体を与える、入射放射線の１０％を減衰する、。
【００３８】
各減衰要素は、いかなる形状を有することもできるが、それが延在する表面１５Ｓの方向の少なくとも１つの大きさは、５ｍｍよりも小さく、好ましくは２ｍｍよりも小さく、更には、１ｍｍよりも小さい。前に説明した実施の形態の全てにおいて、マスクは、好ましくは、検出器の画素が延在する平面に平行なＸＹ平面上に延在する。
【００３９】
上記マスク上の２つの隣り合う減衰要素の間隔は、５ｍｍよりも小さく、好ましくは、１〜５ｍｍの間にあることができる。一般に、検出器２０上に投影した後の２つの隣り合う減衰要素の間隔は、有利には、１〜１０ｃｍの間にあって、好ましくは、５ｃｍよりも小さく、又は、３ｃｍよりも小さい。以下で説明するように、各減衰要素１５Ｘの検出器上への投影は基本陰影領域（elemental shadow zone）の境界を定める。各基本陰影領域は、中心点の周りに拡がる。有利には、２つの隣り合う基本陰影領域の複数の中心点の間の空間は、１〜１０ｃｍの間にあり、好ましくは、１〜５ｃｍの間にある。「投影」は、線源により放出された放射線の伝播の方向での投影を意味すると理解されるべきである。
【００４０】
マスクの例を図１（ｂ）に示す。各減衰要素１５Ｘは、平行六面体の（parallelepipedal）ブロックであり、マスクが延在する表面１５Ｓ上の各減衰要素の領域は、１ｍｍ×１ｍｍであり、各要素１５Ｘの中心の間隔はそれぞれ、第１方向Ｙにおいて３ｍｍであり、第１方向と直角の第２方向Ｘにおいて２．２ｍｍである。この例では、マスクを形成する材料はアルミニウムであり、その厚さは２ｍｍであり、入射放射線の１０％が減衰する、０．９〜１００ｋｅＶの減衰係数に対応する。その厚さは、マスクが延在する表面に直角の方向に伸びる。一般に、マスクは、１次放射線の十分な部分を減衰するような十分に減衰する材料で生成される。しかしながら、入射放射線が物体に到達する前に、入射放射線の顕著な散乱を生じる、例えば、鉛のような、重金属タイプの過剰に密度の高い材料は、好ましくは、避けるべきである。マスクに好ましい材料には、アルミニウム、銅、黒鉛、ホウ素が含まれる。
【００４１】
例えば、異なる減衰要素の間の不規則な間隔、又は、各減衰要素の不規則な幾何学的配置（geometry）を考えることにより、他の幾何学的配置を想像することができる。メッシュを定義する格子形状のマスクもまた可能であり、減衰要素は各メッシュの間に配置される。
【００４２】
各減衰要素の間の、マスクが延在する表面１５Ｓは、好ましくは、考えられるエネルギー範囲において光子を透過すると考えられる材料からなる。それは、薄い、プラスチック、紙、又は、軽金属、アルミニウム、鉄、銅、又は、自由な、空気が占める空間であることができる。このように、各減衰要素１５Ｘの間で、前に定義したように、減衰係数は、好ましくは、０．５よりも小さく、更に、０．２又は０．１よりも小さい。好ましくは、各減衰要素の間で、減衰係数は無視しても良い。
【００４３】
減衰要素の数は、検出器の観察領域に覆うように設定される。説明した例では、減衰要素は、１７×３要素、全５１要素のマトリックスにより分散される。
【００４４】
マスクは、除外又は固定できる。除外できるマスクを有することにより、図１（ｂ）に示すように、マスクなしのスペクトルを生成することが可能になる。
【００４５】
マスクが線源と検出器との間に挿入されたときに、マスクの、入射放射線Ｉ０の伝播の方向での、検出器上への投影は、上記伝播の方向における、各減衰要素１５Ｘに対して位置合わせされた検出器の画素２０Ｘｉ，ｊをグループ化する陰影領域の境界を定める。特に、前に説明したように、入射放射線の伝播の方向での各減衰要素１５Ｘの投影は、検出器上に基本陰影領域を形成する。基本陰影領域の組みは、上記陰影領域を構成する。陰影領域の画素は、Ｇ１で示される第１グループの画素を構成する。この第１グループＧ１は、線源と検出器との間に物体がなくても、計算により又は実験によりあらかじめ決定することができる。この第１グループに属さない画素
【数３】
0006770866.tif 000004
は、減衰要素１５Ｘにより減衰されていない放射線を受け、この第１グループに属する、上記画素と放射線源１１との間に伸びる線上に位置する各画素２０Ｘｉ，ｊ∈Ｇ１は、減衰要素１５Ｘにより減衰された放射線を受ける。
【００４６】
装置は更に演算部、又は、例えば、マイクロプロセッサのようなプロセッサ２１を備え、それは、検出器の画素により計測された各スペクトルＳｉ，ｊを処理することができる。特に、プロセッサは、この明細書で説明するスペクトル処理及び演算動作を実行するための指令のシーケンスを記憶したプログラム可能なメモリ２２に接続したマイクロプロセッサである。これらの指令は、プロセッサにより読み出し可能な、ハードディスク、ＣＤＲＯＭ又は他の種類のメモリの記録媒体上にバックアップをとることができる。プロセッサは、例えば、スクリーンのような、表示部２４に接続されることもできる。
【００４７】
検出器は、検出の不完全性を示す、検出器応答行列Ｄにより特徴付けられることができる。大きさＮ×Ｎ（Ｎは、検出器により形成される各スペクトルのチャネル数を示す）の、この行列を図２（ａ）に示す。この例では、Ｎ ＝ １５０であるが、一般にチャネル数は２よりも大きく、更に１０よりも大きく、数千に及ぶことができる。この行列の各項Ｄ（ｕ，ｖ）は、検出器に入射するエネルギーｖの光子が、検出器によって、エネルギーｕを有すると判断される可能性を示す。
【００４８】
言い換えると、図２（ｂ）に示すような行列の各行Ｄ（ｕ，・）は、エネルギーｕの光子を検出したときに、検出器に到達した光子のエネルギーｖの可能性の分散に対応する。図２（ｂ）は、行列Ｄの行７０を示す。同様に、図２（ｃ）に示すような行列の各列Ｄ（・，ｖ）は、検出器に到達した光子がエネルギーｖを有するときに、検出器により検出されたエネルギーｕの可能性の分散を示す。図２（ｃ）は、行列Ｄの列７０を示す。検出器のエネルギー解像度が高いほど、この行列が対角線行列になる傾向が大きくなる。完全な検出器の場合には、行列Ｄは恒等行列（identity matrix）である。
【００４９】
ここで、図３を参照して、本発明の第１の実施の形態に係る、方法の主要なステップを説明する。
【００５０】
（ステップ１００）  第１スペクトルの獲得
マスク１５が、線源１１と検出器２０との間に挿入される。各画素２０ｉ，ｊは、マスクの存在下で物体を透過した第１放射線と呼ばれる放射線Ｉ１に曝され、第１スペクトルＳ１ｉ，ｊと呼ばれる、そのスペクトルを獲得する。減衰要素１５Ｘの延長上に位置する画素２０ｉ，ｊ∈Ｇ１は第１グループＧ１に属し、他の画素
【数４】
0006770866.tif 000005
と比べて、減衰した放射線を検出する。第１グループの画素により獲得された第１スペクトルはＳ１ｉ，ｊ∈Ｇ１で示される。第１グループに属さない画素により獲得された第１スペクトルは
【数５】
0006770866.tif 000006
で示される。
【００５１】
（ステップ１１０）  第１グループＧ１の各画素２０ｉ，ｊ∈Ｇ１について、線源と物体との間に挿入されたマスクなしに、物体を透過した第２放射線Ｉ２を示す第２スペクトルＳ２ｉ，ｊ∈Ｇ１を取得すること
この第２スペクトルは、図１（ｂ）に示す、マスクが除外された配置により、実験的に取得することができる。この第２スペクトルは、第１グループＧ１に属さず、減衰要素により減衰されていない放射線を受ける第１スペクトル
【数６】
0006770866.tif 000007
から、評価により取得することができる。この評価は、例えば、線形補間又は多項式補間のような補間により行うことができる。それは、実験計測を行う必要なしに、マスクが存在しない場合に、第１グループの画素により受けられる第２スペクトルＳ２ｉ，ｊ∈Ｇ１を評価することを可能にする。
【００５２】
図５（ａ）、５（ｂ）に関連して、後で説明する他の実施の形態によると、各計測は対応する計測配置Ｃを有し、それは、物体に対する検出器及び／又は線源の位置を示す。この計測配置Ｃは２つの連続する計測の間で変化する。配置Ｃによる、検出器の第１グループの画素が受ける第２スペクトルＳ２，Ｃｉ，ｊ∈Ｇ１，Ｃは、上記画素が減衰要素の延長上になく、他の配置において第１グループＧ１，Ｃ’に属していない、他の配置Ｃ’において生成される第１スペクトル
【数７】
0006770866.tif 000008
から評価できる。この特別な場合については、図５（ａ）、５（ｂ）の記載において、後でより詳細に説明する。
【００５３】
（ステップ１２０）  第１グループＧ１の各画素２０ｉ，ｊ∈Ｇ１に関して、同じ画素について第１スペクトルＳ１ｉ，ｊ∈Ｇ１と第２スペクトルＳ２ｉ，ｊ∈Ｇ１とを比較することにより得られる、いわゆる比較スペクトルΔＳｉ，ｊ∈Ｇ１を決定すること
例えば、比較は減算であって、その場合に、
ΔＳｉ，ｊ∈Ｇ１ ＝ Ｓ２ｉ，ｊ∈Ｇ１ − Ｓ１ｉ，ｊ∈Ｇ１  （２）
である。スペクトルを減算することは、ベクトル減算（vector subtraction）タイプの、各チャネルの内容の減算を意味すると理解されるべきである。
【００５４】
（ステップ１３０）  上記画素と線源との間に配置されたマスクなしで、第１グループの各画素について、いわゆる１次スペクトル
【数８】
0006770866.tif 000009
を評価すること
そのようなスペクトルは、第１グループＧ１の画素に到達する第２スペクトルＳ２ｉ，ｊ∈Ｇ１の１次成分を示す。この評価は、このグループの各画素について、行列積
【数９】
0006770866.tif 000010
を介して、ステップ１２０で決定した比較スペクトルΔＳｉ，ｊ∈Ｇ１を検討することにより行われる。
【００５５】
Ｍは、あらかじめ決定され、ステップ１２０で決定した比較スペクトルと、マスクなしの場合に第１グループの画素２０ｉ，ｊ∈Ｇ１に到達した放射線Ｉ２の１次成分の評価
【数１０】
0006770866.tif 000011
との間の関係を得る、遷移行列である。
【００５６】
この評価は、本発明の重要な態様を構成する。その起源は、マスクあり及びマスクなしでそれぞれ第１グループの画素２０ｉ，ｊ∈Ｇ１が計測したスペクトルＳ１ｉ，ｊ∈Ｇ１、Ｓ２ｉ，ｊ∈Ｇ１を説明することにある。それは、マスクあり又はマスクなしでの散乱放射線のスペクトルＳdiffｉ，ｊ∈Ｇ１の変化を無視できるとする推定に基づく。それゆえ、その有効性は、前に説明したように、減衰し散乱しないマスクの利用を条件としている。
【００５７】
式（１）を用いて、Ｓ１ｉ，ｊ∈Ｇ１は、
Ｓ１ｉ，ｊ∈Ｇ１ ＝ Ｓ１Ｐｉ，ｊ∈Ｇ１ ＋ Ｓ１ｄｉｆｆｉ，ｊ∈Ｇ１  （４）
の形式で表現でき、Ｓ１Ｐｉ，ｊ∈Ｇ１及びＳ１ｄｉｆｆｉ，ｊ∈Ｇ１はそれぞれスペクトルＳ１ｉ，ｊ∈Ｇ１の直接又は散乱スペクトルを示す。
【００５８】
【数１１】
0006770866.tif 000012
であり、Ｄは前に説明した検出器応答行列を示し、Ｓ０は入射放射線Ｉ０のスペクトルを示し、ａｔｔ１０、ａｔｔ１５はそれぞれスペクトルの各エネルギーで画素２０ｉ，ｊに到達した放射線の物体１０及びマスク１５中での減衰を示す。
ａｔｔ１０ ＝ Σｑｌｑμｑ
ａｔｔ１５ ＝ｌ１５Ｘμ１５Ｘ
であり、ｌｑは物体１０の各要素１０ｑ中での移動距離であり、μｑはスペクトルの各エネルギーにおける、各要素１０ｑの線形減衰係数のベクトルであり、ｌ１５Ｘは画素２０ｉ，ｊと結合した減衰要素１５Ｘ中での移動距離であり、μ１５Ｘはスペクトルの各エネルギーにおける、線形減衰係数のベクトルである。×は行列積を示し、
【数１２】
0006770866.tif 000013
はアダマール積を示す。線形減衰係数μｑは次の定義
ｌｑμｑ（Ｅ） ＝ −ｌｎ［Ｉｑ（Ｅ）／Ｉ０，ｑ（Ｅ）］
を証明し、Ｉ０，ｑ、Ｉｑはそれぞれ厚さｌｑの物体の要素１０ｑに入射した放射線及び透過された放射線を示す。当業者に知られたそのような線形減衰係数は、エネルギー及び要素１０ｑを構成する材料に依存する。
【００５９】
スペクトルＳ２ｉ，ｊ∈Ｇ１は、減衰スクリーン１５なしの獲得に対応するので、
Ｓ２ｉ，ｊ∈Ｇ１ ＝ Ｓ２Ｐｉ，ｊ∈Ｇ１ ＋ Ｓ２ｄｉｆｆｉ，ｊ∈Ｇ１  （６）
【数１３】
0006770866.tif 000014
である。
【００６０】
このため、推定により、
Ｓ２ｄｉｆｆｉ，ｊ∈Ｇ１ ＝Ｓ１ｄｉｆｆｉ，ｊ∈Ｇ１
なので、
【数１４】
0006770866.tif 000015
である。このため、
【数１５】
0006770866.tif 000016
であり、１は１だけを備えるベクトルである。したがって、
【数１６】
0006770866.tif 000017
であり、このため、
ΔＳｉ，ｊ∈Ｇ１ ＝ Ｄ × ｗＳ２Ｐｉ，ｊ∈Ｇ１  （１２）
で、
ｗ ＝ （１−ｅ−ａｔｔ１５）  （１３）
である。
【００６１】
Ｗ−１で示す、大きさＮ×Ｎの正方行列を構成することが可能であり、この行列の項は、ゼロである対角線上に位置せず、対角線の各項ｎは、値
Ｗ−１（ｎ，ｎ） ＝ １／（１−ｅ−ａｔｔ１５（ｎ））’  （１４）
を有し、ａｔｔ１５（ｎ）はランク（rank）ｎのチャネルに対応するエネルギーでのスクリーン１５の減衰を示す。係数ｎはスペクトルのチャネルのランクを示し、それはエネルギー値Ｅに例えることができる。減衰は、減衰材料を構成する材料の線形減衰係数により乗算される減衰要素の厚さに対応する。
【００６２】
このように、画素２０Ｓｉ，ｊ∈Ｇ１について、比較スペクトルΔＳｉ，ｊ∈Ｇ１を決定する式（１３）（１４）に従い、あらかじめ定義した検出器応答行列Ｄ及び行列Ｗ−１を用いて、式
【数１７】
0006770866.tif 000018
により、１次スペクトル
【数１８】
0006770866.tif 000019
を評価することが可能である。
【００６３】
この式は、次のように書き換えることができ、
【数１９】
0006770866.tif 000020
それは、前に書いた式（３）に対応する。
【００６４】
ここで、Ｍは、
Ｍ ＝ Ｄ × Ｗ−１ × Ｄ−１  （１７）
のような遷移行列である。
【００６５】
検出器応答行列Ｄ及びマスク１５の減衰を考慮した行列から生成された遷移行列Ｍは、減衰あり及び減衰なしでそれぞれ決定されたスペクトルＳ１ｉ，ｊ∈Ｇ１、Ｓ２ｉ，ｊ∈Ｇ１の間の比較ΔＳｉ，ｊ∈Ｇ１から第１グループの画素２０ｉ，ｊ∈Ｇ１に減衰なしで到達する放射線の１次スペクトル
【数２０】
0006770866.tif 000021
を評価することを可能にする。この例では、比較はこれらの２つのスペクトルの間の差異である。比較スペクトルΔＳｉ，ｊから減衰がない１次スペクトルの評価に分析的に取り組むことを可能にすることから、この遷移行列を用いることは、本発明の重要な要素である。
【００６６】
遷移行列Ｍは、分析的に決定することができるが、本発明者は、以下で説明するように、それを実験的に評価することを可能にする方法を提案した。
【００６７】
（ステップ１４０）  第１グループの各画素２０ｉ，ｊ∈Ｇ１についての散乱スペクトル
【数２１】
0006770866.tif 000022
の評価
【００６８】
第１グループの各画素２０ｉ，ｊ∈Ｇ１について、Ｍ及びΔＳｉ，ｊ∈Ｇ１からスペクトル
【数２２】
0006770866.tif 000023
を評価することにより、（６）を用いて、このグループの各画素２０ｉ，ｊ∈Ｇ１に影響を及ぼす散乱スペクトルを評価することで可能である。
【数２３】
0006770866.tif 000024
【００６９】
（ステップ１５０）  検出器の各画素２０ｉ，ｊについての散乱スペクトルの
【数２４】
0006770866.tif 000025
の評価
第１グループの各画素２０ｉ，ｊ∈Ｇ１の散乱スペクトルを評価することにより、散乱放射線の空間的変化の頻度は小さいとする推定が基礎として用いられる。言い換えると、１つの画素から他の隣の画素への、散乱放射線のスペクトルの変動は、突然に変化しない。そして、第１グループＧ１の画素の散乱スペクトル
【数２５】
0006770866.tif 000026
の補間により、検出器の全ての画素の散乱スペクトル
【数２６】
0006770866.tif 000027
を評価することが可能である。第１グループに入らない画素については、これらの画素がマスクによる減衰を受けないので、散乱スペクトルはマスクあり及びマスクなしで当然同じである。
【００７０】
補間の前に、いくつかの任意の前処理動作を実行することが可能である。同じ減衰要素１５ｋの陰に配置された画素として２０ｉ，ｊ∈ｋを例に挙げる。これらの画素は核（kernel）ｋを形成し、核ｋの中心画素と呼ばれる１つの画素２０ｋが中心を構成する。この核を示す散乱スペクトル
【数２７】
0006770866.tif 000028
を評価し、それを中心画素に割り当てることが可能である。核ｋを示すこのスペクトルは
【数２８】
0006770866.tif 000029
であって、Ｃｋは核ｋの画素数を示す。スペクトル
【数２９】
0006770866.tif 000030
は、核ｋの画素の散乱スペクトルの評価の中心値を示す。補間は、各核ｋを示す各スペクトル
【数３０】
0006770866.tif 000031
を考慮して行われる。
【００７１】
他の任意の前処理動作は、検出器の全ての画素に到達する散乱スペクトルは同じ形状を有する、つまり、画素毎に同様であるとする推定に基づく、各核ｋを示す各スペクトル
【数３１】
0006770866.tif 000032
の拡大縮小（scaling）にある。この推定は、検出器の表面領域が小さい、典型的には数ｃｍ２又はそれよりも小さいときに、なお一層正当化される。この拡大縮小は、核ｋを示す散乱スペクトル
【数３２】
0006770866.tif 000033
の各チャネルｎにスカラーａｋ，ｎを適用することにあり、このスカラーは式
【数３３】
0006770866.tif 000034
により決定される。
【００７２】
Ｋは核ｋの数を示す。そして、各チャネルｎについて、
【数３４】
0006770866.tif 000035
で示される修正スペクトルを取得するような方法で、各核ｋでの散乱を示すスペクトル
【数３５】
0006770866.tif 000036
を修正することが可能である。
【数３６】
0006770866.tif 000037
【００７３】
２次元ＸＹに配置された検出器の画素を用いて、Ｘ方向について、次に他方のＹ方向について、又は、Ｘ方向及びＹ方向の両方について同時に、スペクトル
【数３７】
0006770866.tif 000038
【数３８】
0006770866.tif 000039
【数３９】
0006770866.tif 000040
に基づいて、用いる変数により、補間が行われる。補間はエネルギーチャネル毎に行われる。それは多項式であることができる。検討する例では、検出器は６４×６４０画素を備える。補間は、６４０画素の各行について４次多項式に従って、次に、６４画素の各列について２次多項式に従って行われる。
【００７４】
（ステップ１６０）  物体を透過した１次放射線のスペクトルを示す修正スペクトルＳ＊ｉ，ｊを取得するための、検出器の画素２０ｉ，ｊの全て又はいくつかについての第２スペクトルＳ２ｉ，ｊ（又は、第１グループに属さない画素についての第１スペクトルＳ１ｉ，ｊ）の修正
言い換えると、もしもＳＰｉ，ｊが物体を透過して画素２０ｉ，ｊ上に到達した１次放射線のスペクトルを示すならば、
【数４０】
0006770866.tif 000041
である。
【００７５】
このステップは単純に、各画素２０ｉ，ｊについて、第２スペクトルから（又は、第１グループに属さない画素についての第１スペクトルから）、散乱スペクトル
【数４１】
0006770866.tif 000042
の評価値を減算することにより行われる。
【００７６】
すなわち、第１グループの画素２０ｉ，ｊ∈Ｇ１について、
【数４２】
0006770866.tif 000043
であり、第１グループに属さない画素
【数４３】
0006770866.tif 000044
について、
【数４４】
0006770866.tif 000045
である。
【００７７】
そして、検出器の各画素について、物体を透過した１次放射線のスペクトルの評価値を示すスペクトルＳ＊ｉ，ｊが得られる。そして、各画素により検出された、１つ又は複数のエネルギー（Ｅ）における１次放射線の強度、すなわち、スペクトル画像の意味を示す画像Ｉｍ（Ｅ）を得ることが可能である。それは、計測が解釈されるこのスペクトル画像に基づく。
【００７８】
断層再構成の場合には、前に説明したように、方法は異なる配置Ｃで実行されるが、各配置は、断層再構成が生成されることに基づく、対応するスペクトル画像又は対応する複数のスペクトル画像ＩｍＣ（Ｅ）を有する。
【００７９】
どのような実施の形態であれ、異なるチャネルＥがグループ化されて、スペクトル帯（spectral band）ΔＥを構成することができる。そして、これらのスペクトル帯のそれぞれについて、画像Ｉｍ（ΔＥ）を取得することが可能である。
（遷移行列を得ること）
【００８０】
本発明の重要な要素の１つは、ステップ１３０で用いる遷移行列Ｍである。この行列は式（１７）による計算により取得することができる。しかしながら、本発明者は、遷移行列を実験的に得ることが好ましいと評価した。行列Ｍの各行Ｍ（ｃ，ｎ）は、
Ｍ（ｃ，ｃ） ＝ １／（１−ｅ−ａｔｔ１５（Ｃ））  （２５）
のような、対角項Ｍ（ｃ，ｃ）の一方の側での変形を促進する関数であると考えることができる。
【００８１】
このように、行Ｍ（ｃ，ｎ）の各項は、式
Ｍ（ｃ，ｎ） ＝ ［１／（１−ｅ−ａｔｔ１５（Ｃ））］ × ｆ（ｃ，ｎ，α１・・・αＰ）  （２６）
により説明することができ、ここで、ｆは変形関数であり、α１・・・αＰは変形関数ｆのパラメータである。
【００８２】
遷移行列を得るには、変形関数ｆの分析的表現（analytical expression）を定義する必要がある。そのために、シミュレーションが異なる校正配置（calibration configuration）により行われるが、各校正配置Ｃａｌｉｂは既知の性質及び厚さの材料に対応する。これらの配置は、校正の基礎（calibration base）を構成する。各校正配置Ｃａｌｉｂについて、第１グループＧ１に属する検出器の複数の画素を考慮して、図４（ｃ）に関連して説明するように、次のステップが行われる。
（ステップ９０）  １組みのパラメータ（α１・・・αＰ）によりパラメータ化された、関数ｆの分析形態の選択
（ステップ９２）  各校正配置Ｃａｌｉｂについての、１次放射線のスペクトルＳＣａｌｉｂ，Ｐｉ，ｊ∈Ｇ１の、シミュレーションによる決定
（ステップ９４）  各配置についての、マスクあり及びマスクなしで検出器に到達した放射線のスペクトルの、シミュレーション又は計測による決定、及び、これらの２つのスペクトルの間の差異を示す、比較スペクトルΔＳＣａｌｉｂｉ，ｊ∈Ｇ１の計算
（ステップ９６）  異なる校正配置に対応するスペクトルＳＣａｌｉｂ，Ｐｉ，ｊ∈Ｇ１を用いた、式
【数４５】
0006770866.tif 000046
によるパラメータ（α１・・・αＰ）の組みの決定。ここで、
【数４６】
0006770866.tif 000047
は、パラメータ（α１・・・αＰ）によりパラメータ化された遷移行列Ｍを式
【数４７】
0006770866.tif 000048
により実行することにより得られる１次スペクトルＳＣａｌｉｂ，Ｐｉ，ｊ∈Ｇ１の評価値であり、ψはスペクトルＳＣａｌｉｂ，Ｐｉ，ｊ∈Ｇ１の組みと、パラメータ（α１・・・αＰ）によりパラメータ化された遷移行列Ｍを用いて得られた、それらのそれぞれの評価値
【数４８】
0006770866.tif 000049
との間の誤差を示す関数である。それは、２次誤差（quadratic error）のタイプの関数であることができる。
（ステップ９８）  ステップ９６でパラメータが決定された関数ｆを用いて、及び、式（２６）を用いて、行列Ｍを取得すること
（他の実施の形態）
【００８３】
図５（ａ）、５（ｂ）に関連して説明する、１つの実施の形態によれば、物体に対する検出器及び／又は線源の位置は、２つの連続する計測の間で変化する。更に、各計測は対応する配置Ｃを有し、各配置は、前に定義したように、物体１０に関する検出器２０及び／又は線源１１の相対位置により特徴付けられる。各配置は、物体に関する検出器及び／又は線源の位置を示すパラメータＰＣにより特徴付けられる。このパラメータは、スカラー又はベクトルであることができる。図５（ａ）、５（ｂ）に示す例では、線源は検出器に対して固定され、線源及び検出器からなる組立物は、ＸＺ平面に直角なＹ軸上で、初期位置に対して角度Θで、物体の周りを回転する。この角度は、各計測配置のパラメータＰＣを構成する。
【００８４】
図５（ａ）に示す第１配置Ｃでは、Θ＝０である。画素２０ｌは、マスク１５の減衰要素１５Ｘにより遮蔽される。図５（ｂ）に示す第２配置Ｃ’では、Θ＝Θ’＞０である。画素２０ｌは、マスク１５により遮蔽されない。そのような実施の形態において、画素ｌについて、図５（ａ）に示す配置Ｃにより第１スペクトルＳ１，Ｃｌが得られ、第２スペクトルＳ２，Ｃｌが配置Ｃ’により得られる第１スペクトルＳ１，Ｃ’ｌに対応する。言い換えると、
Ｓ２，Ｃｌ ＝ Ｓ１，Ｃ’ｌ
である。この実施の形態は、線源と物体との間に固定マスクを保持することを可能とし、画素２０ｉ，ｊの第２スペクトルＳ２，Ｃｉ，ｊが、配置Ｃ’の第１スペクトルＳ１，Ｃ’ｉ，ｊの計測により得られる。第１配置Ｃでは、画素２０ｉ，ｊは第１グループＧ１，Ｃに属し、第２配置Ｃ’では、この画素は第１グループＧ１，Ｃ’に属さない。
【００８５】
好ましくは、これらの２つの配置は、第２配置で画素２０ｌにより計測されたスペクトルに十分に匹敵し、第１配置のマスクなしの放射線を示す。このように、好ましくは、配置Ｃ、Ｃ’にそれぞれ関連付けられたパラメータＰ、Ｐ’は類似する。「類似」は、それらの相対偏差（relative deviation）が所定のしきい値、例えば、１０又は２０％を超えないことを意味すると理解されるべきである。
【００８６】
１つの実施の形態によると、各配置において、複数の画素について、修正スペクトルＳ＊，Ｃｉ，ｊを取得するために、方法は、各配置が物体に対する検出器及び線源の関連位置をそれぞれ有するような、複数の配置Ｃにより実行される。各配置と関連付けられた修正スペクトルＳ＊，Ｃｉ，ｊは、物体１０の断層再構成を生成するために用いられる。各配置は、線源−検出器組立物及び物体により形成される角度を示す、パラメータＰＣ＝ΘＣによりパラメータ化されることができる。
【００８７】
再構成は、物体を形成する要素１０ｑの形状、及び、それらの線形減衰係数μｑ又は減衰を表す何らかの他の係数を再構成することを目的とする。それは、各エネルギーＥに対応するスペクトル画像ＩｍＣ（Ｅ）に基づいて、既知の逆転アルゴリズム（inversion algorithm）を実行することにより生成され、各画像は上記エネルギーでの放射線量を示し、この量は修正スペクトルＳ＊，Ｃｉ，ｊにより得られる。これらの修正スペクトルを用いることは、以下の例で説明するように、断層再構成の質を顕著に改善することを可能にする。
【００８８】
再構成の前に、スペクトル帯のそれぞれについてスペクトルＳ＊，Ｃｉ，ｊ（ΔＥ）を形成し、スペクトル帯に対応するスペクトル画像ＩｍＣ（ΔＥ）を取得するために、典型的には数十ｋｅＶのスペクトル帯ΔＥのスペクトルのグループ化を進めることができる。
（例）
【００８９】
次に、計算コード（computation code）により実行されるシミュレーションに基づく、Ｘ線撮影での応用例及び断層撮影での応用例に係る、本発明の例示的な実施の形態の説明を続ける。まず、遷移行列Ｍを得ることの例を説明する。
【００９０】
図４（ｃ）に関連して説明するように、遷移行列は１つ以上の材料を備える校正の基礎を用いて決定される。この例では、２つの材料（水及びアルミニウム）が検討され、検討された厚さは、
アルミニウムについては、ピッチ０．１ｃｍで０〜０．２ｃｍ
水については、ピッチ０．５ｃｍで０〜２０ｃｍ
である。関数ｆの解析形式は
【数４９】
0006770866.tif 000050
【数５０】
0006770866.tif 000051
のように選択される。
【００９１】
係数ｃ及びｎは、遷移行列の行及び列とそれぞれ関連付けられる。
【００９２】
その上、各校正配置について、比較スペクトルΔＳＣａｌｉｂｉ，ｊ∈Ｇ１を形成するために、実行された検出器の第１グループの画素に到達した１次スペクトルＳＣａｌｉｂ，Ｐｉ，ｊ∈Ｇ１はモデル化され、スペクトルの模擬計測（simulated measurement）はマスクあり及びマスクなしで実行された。そして、前に説明したように、値α１、α２が決定された。結果は、α１＝０．２であり、０．０４≦α２≦０．１８であった。図４（ａ）は取得した遷移行列Ｍを示し、図４（ｂ）はこの行列の座標ｎ＝４５の列を示す。
【００９３】
次に、図１（ａ）、１（ｂ）に示す装置を実行するシミュレーションの説明を続ける。物体は、人間の胴体をまねた人体模型である。それは、母体１００を備え、ＸＺ平面に平行な平面での、母体の断面は楕円形であり、Ｚ軸上及びＸ軸上での最大の大きさはそれぞれ２０ｃｍ及び４０ｃｍである。この母体１００は体を表し、水でできている。母体１００内部に配置された第１要素１０１及び第２要素１０２は肺を示す。ＸＺ平面に平行な平面での、それらの断面は楕円形であり、Ｚ軸上及びＸ軸上での最大の大きさはそれぞれ１１ｃｍ及び１５ｃｍである。これらの要素は空気からなる。第３要素１０３は脊柱を示す。ＸＺ平面に平行な平面での、その断面は直径３ｃｍの円盤の形状である。長方形の第４要素１０４は椎骨を示す。その構成材料は骨であり、国際放射線単位測定委員会（ＩＣＲＵ：International Commission on Radiation Units and Measurements）のデータベースによりモデル化される。最後に、楕円形の断面の環状の骨からなる母体の周辺１０５は側面を示す。その厚さは１ｃｍである。直径２ｃｍの実質的に球状の介在物（spherical inclusion）１０７、１０８、１０９は、第１要素及び第２要素の高さに配置され、肺を示す。これらの介在物はＰＭＭＡ（ポリメチルメタクリレート）からなり、癌性腫瘍を示す。
【００９４】
線源１１は、電圧１１０ｋＶを受ける、タングステンアノードを有するＸ線管である。検出器２０は（Ｘ軸上で）６４０画素×（Ｙ軸上で）６４画素であり、各画素は厚さ５ｍｍのＣｄＴｅを備える。検出器が延在するＸＹ平面での各画素の表面領域は、１ｍｍ×１ｍｍである。検出器はエネルギー分析され、各画素は、１ｋｅＶのエネルギーチャネルによりスペクトルを取得することを可能にする。用いられるマスク１５は、図２に示すものである。
【００９５】
図６（ａ）、６（ｂ）、６（ｃ）は、検出器の全ての画素上の、マスクなしでの全放射線、１次放射線、及び、図３に示す方法により本発明を実行することにより評価された１次放射線をそれぞれ示す。方法の実行中に、図２に示すように、線源から１０ｃｍの距離に位置するマスクをシミュレーションすることにより、第２スペクトルＳ２ｉ，ｊ∈Ｇ１が得られる。各画素２０ｉ，ｊに対応し、色コード（colour code）により示される量は、
【数５１】
0006770866.tif 000052
のような、減衰値Ａｉ，ｊである。ここで、
Ｓ０ｉ，ｊは、検出器が線源に直接さらされたときに、画素２０ｉ，ｊにより計測されるスペクトルである。それは、入射放射線Ｉ０のスペクトルを示す。
Ｓ’ｉ，ｊは、画素２０ｉ，ｊにより計測される、目的のスペクトルである。
図６（ａ）では、Ｓ’ｉ，ｊは、マスクなしで、物体を透過した放射線のスペクトルを示す。
Ｓ’ｉ，ｊ＝Ｓｉ，ｊ
図６（ｂ）では、Ｓ’ｉ，ｊは、マスクなしで、線源と検出器との間に配置された物体を透過した放射線の１次スペクトルを示す。
Ｓ’ｉ，ｊ＝ＳＰｉ，ｊ。
図６（ｃ）では、Ｓ’ｉ，ｊは、本発明により修正されたスペクトルを示し、それは、１次スペクトルＳｉ，ｊの評価値に対応する。
【数５２】
0006770866.tif 000053
【００９６】
各量Ａｉ，ｊは、物体中での全減衰（global attenuation）を示し、用語「全（global）」は、それがスペクトルの全てのエネルギーチャネルで決定されるという事実を説明する。そのことは、図示を簡単にする、スカラーにより目的の各スペクトルＳ’ｉ，ｊを示すことを可能にする。図６（ａ）、６（ｂ）、６（ｃ）は次のことを示す。
本発明を実行することによる、１次放射線のスペクトル
【数５３】
0006770866.tif 000054
の評価値（図６（ｃ））は、モデル化された１次放射線のスペクトルＳＰｉ，ｊ（図６（ｂ））と一致する。
１次放射線のスペクトル（図６（ｂ）、６（ｃ））は、全放射線のスペクトル（図６（ａ））よりも大きなコントラストを示す。１次放射線のスペクトルの分離は、異なる減衰を示す領域のより良い空間的な分離を可能にする。そのことは、密度の高い要素１０３、１０４（骨）に対応する、画像の中央部分で特に観ることができる。１次放射線のスペクトルに頼ることは、現実に近い減衰値（attenuation value）を取得することを可能にする。
【００９７】
これらの結果は、これらの画像のそれぞれの中心形状を同一の値ｙ＝３２で示す、図６（ｄ）で確かめられる。この図では、係数Ａ、Ｂ、Ｃはそれぞれ図６（ａ）、６（ｂ）、６（ｃ）を言及する。
【００９８】
図７（ａ）、７（ｂ）、７（ｃ）、７（ｄ）は、それぞれ、検出器の座標（ｉ，ｊ）＝（１００，１５）、（１９２，３８）、（３２，５５）、（３２０，３２）の画素２０ｉ，ｊ上の、全放射線のスペクトルＳｉ，ｊ、本発明によりモデル化された１次放射線のスペクトルＳＰｉ，ｊ、本発明により評価された１次放射線のスペクトル
【数５４】
0006770866.tif 000055
を示す。これらの画像のそれぞれにおいて、係数１、２、３はそれぞれ全放射線のスペクトルＳｉ，ｊ、本発明によりモデル化された１次放射線のスペクトルＳＰｉ，ｊ、本発明により評価された１次放射線のスペクトル
【数５５】
0006770866.tif 000056
を示す。ポイントのそれぞれにおける、本発明によりモデル化された１次放射線のスペクトルと、本発明により評価された１次放射線のスペクトルとの間の一致は、注目されるであろう。図７（ｄ）に示す、（ｉ，ｊ）＝（３２０，３２）の画素は、検出された放射線のほとんど全てが散乱放射線であり、１次放射線の寄与は無視できる点において注目に値する。この画素は、強く減衰する骨要素１０３、１０４に面して位置している。そのことは、この画素により計測される放射線が、基本的に、物体から発する散乱放射線からなることを意味している。本発明の目的及び効果は、ここにおいて判断される。
【００９９】
散乱スペクトルを修正するための方法の効果は、散乱線対１次線比（Scatter-to-Primary-Ratio）の略語ＳＰＲで示される指数により、しばしば計測される。
【０１００】
この指数の、積分表現（integral expression）と呼ばれる第１の表現は
ＳＰＲ ＝ ［ΣｎΣｉ，ｊｎＳｉ，ｊ（ｎ）／ΣｎΣｉ，ｊｎＳＰｉ，ｊ（ｎ）］ × １００  （３３）
である。
【０１０１】
用語「積分」は、指数が全ての画素及びエネルギーチャネルに基づくという事実を説明する。
【０１０２】
この指数はピクセル化される、すなわち、検出器の各画素２０ｉ，ｊと関連付けられることができ、その場合には、それは空間表現（spatial expression）と呼ばれる、第２の表現
ＳＰＲ（ｉ，ｊ） ＝ ［ΣｎｎＳｉ，ｊ（ｎ）／ΣｎｎＳＰｉ，ｊ（ｎ）］ × １００  （３４）
により説明される。
【０１０３】
それは、更に、検出器の全ての画素上のエネルギーの関数として表現されることができ、その場合には、スペクトル表現（spectral expression）と呼ばれる第３の表現
ＳＰＲ（ｎ） ＝ ［Σｉ，ｊｎＳｉ，ｊ（ｎ）／Σｉ，ｊｎＳＰｉ，ｊ（ｎ）］ × １００  （３５）
により説明される。
【０１０４】
その表現が積分、空間、スペクトルのいずれであっても、係数ＳＰＲが小さくなるほど、画素２０ｉ，ｊにおいて、又は、エネルギー範囲ｎにおいて、検出器の全ての上で散乱される割合（share）は小さくなる。
【０１０５】
更に、検出器の各画素２０ｉ，ｊについての評価の有効性を反映する、規格化微分指数（normalized differential indicator）Ｉｎｄ＿ｄｉｆｆにより１次放射線のスペクトルの評価の有効性を定量化することが可能である。前に説明した係数ＳＰＲと同様に、この指数Ｉｎｄ＿ｄｉｆｆは、それぞれ、表現
【数５６】
0006770866.tif 000057
【数５７】
0006770866.tif 000058
【数５８】
0006770866.tif 000059
により、積分、空間、スペクトルの形式で表現することができる。
【０１０６】
その表現が積分、空間、スペクトルのいずれであっても、係数Ｉｎｄ＿ｄｉｆｆが小さくなるほど、評価値は良くなる。
【０１０７】
表１は、検出器の全ての画素上での、上で説明した異なる指数を示す。
【表１】
0006770866.tif 000060
【０１０８】
検討した画像において、信号対１次積分比（signal to primary integral ratio）は２０％に近い。係数Ｉｎｄ＿ｄｉｆｆは平均してとても小さいことがわかり、そのことは、評価の良い品質を証明する。
【０１０９】
物体の周辺での、Ｚ軸に平行な回転軸での、計測装置（線源、マスク、検出器）の回転をシミュレーションすることにより、そして、１度の各角度間隔での獲得を実行することにより、図１に示す物体の断層再構成が生成された。再構成の前に、スペクトルのチャネルは、２３〜１１０ｋｅＶの間で、２１ｋｅＶ幅の４つのスペクトル帯にグループ化される。これらのスペクトル帯は、２３〜４４ｋｅＶ、４５〜６６ｋｅＶ、６７〜８８ｋｅＶ、８９〜１１０ｋｅＶである。図８（ａ）、８（ｂ）、８（ｃ）は、それぞれ、
全放射線のスペクトルＳｉ，ｊ
シミュレーションした１次放射線のスペクトルＳＰｉ，ｊ
各角度間隔で、マスクあり及びマスクなしの計測をシミュレーションした後に、本発明により評価された１次放射線のスペクトル
【数５９】
0006770866.tif 000061
に基づいて、再構成された物体の減衰の断面図を示す。
【０１１０】
各スペクトルは、１ｋｅＶのスペクトル解像度によりモデル化された。そして、前に説明した４つのスペクトル帯を取得するためにチャネルのグループ化が実行された。各スペクトル帯による再構成を図８（ａ）、８（ｂ）、８（ｃ）に示す。
【０１１１】
シミュレーションされた１次放射線に基づいて生成された再構成と、本発明により評価された１次放射線に基づいて生成された再構成との間の良い一致は、注目されるであろう。１次放射線を考慮することは、エネルギーが低いときに相対的により大きな効果を有することを観ることができる。
【０１１２】
人の胸部を示す擬人化された人体模型（擬人化された胸部人体模型（Anthropomorphic thorax phantom））を用いて、実験試験が行われた。実験条件は、以下の通りである。
線源：Ｘ線管YXLON Y.TU 160-D06  動作電圧１１０ｋＶ
検出器：  １２８×１画素MutliX ME 100エネルギー分解ＣｄＴｅ検出器  間隔０．８ｍｍ  厚さ３ｍｍ  線形検出器が１２８×４５１画素の２次元画像を獲得するために変形された。
直径２ｍｍ、高さ２ｍｍのアルミニウムシリンダの１１×２１のアレイからなるマスク、２つの隣り合うシリンダーの間の距離は５ｍｍである。
アルミニウムシリンダ及び直径２ｃｍのＰＭＭＡ（ポリメチルメタクリレート）を用いて得られた遷移行列Ｍ。アルミニウムシリンダの厚さは、０〜１６ｃｍの間で４ｃｍずつ変化する。ＰＭＭＡシリンダの厚さは、０〜３ｃｍの間で１ｃｍずつ変化する。
【０１１３】
図９（ａ）、９（ｂ）、９（ｃ）は、それぞれ、式（３２）に関連して説明したように、エネルギー範囲２５〜１１０ｋｅＶの範囲で、次のものに基づいて計算された減衰Ａｉ，ｊを示す。
各画素２０ｉ，ｊにより計測された、透過放射線のスペクトル
散乱放射線の評価値が減算された、透過放射線のスペクトル。この評価は、ビームストップ（Beam Stop）と呼ばれる、減衰帯を設けることによりなされる。この減衰帯に面して位置する画素が、散乱スペクトルを示すスペクトルを計測する。このスペクトルは、次に、各画素についての１次スペクトルの評価値を得るために、透過スペクトルから減算される。これは、参照修正方法（reference correction method）である。
各画素２０ｉ，ｊ上での、本発明により評価された１次放射線のスペクトル
【０１１４】
図９（ｂ）及び図９（ｃ）において良い一致を観ることができ、そのことは、本発明による修正の有効性を証明する。
【０１１５】
図９（ｂ）に、水平線を示す。図９（ｄ）は、この線に沿った、図９（ａ）〜９（ｃ）のそれぞれの上の強度プロファイルを示す。図９（ｂ）及び図９（ｃ）から得られるプロファイルの間に良い一致を観ることができる。この図において、係数Ａ、Ｂ、Ｃはそれぞれ図９（ａ）、９（ｂ）、９（ｃ）を言及する。
【０１１６】
図１０（ａ）、１０（ｂ）、１０（ｃ）、１０（ｄ）は、図９（ｂ）に示す各点１、２、３、４上の、透過放射線のスペクトル、参照方法により修正された透過放射線のスペクトル、本発明により修正された透過放射線のスペクトルを示す。これらの図面のそれぞれにおいて、係数１、２、３はそれぞれ、全放射線のスペクトルＳｉ，ｊ、参照方法を用いた１放射線のスペクトルＳＰｉ，ｊ、本発明により評価された１放射線のスペクトル
【数６０】
0006770866.tif 000062
を示す。
【０１１７】
本発明は、医療応用例のために、又は、より一般的に、上記物体の内部構造を調査することを目的とする、物体の非破壊検査において、電離放射線、特に、Ｘ線又はガンマ線を提供するスペクトル画像解析方法に適用することができる。物体は、例えば、荷物、工業製品、例えば、パイプラインのような設備の構造要素、核廃棄物などであることができる。
【０１１８】
本発明は、このように、散乱放射線の影響を限定して、放射線の１次成分の評価を可能にする。取得した画像の質、特に、空間解像度が、改善される。その結果は、より正確な結果であり、調査した物体に対してより共形（conformal）な結果である。
【０１１９】
方法を実行することは簡単であり、現存する装置に適合させることができる。その上、遷移行列をあらかじめ得ることができ、大きな計算時間を必要としないで、方法の速やかな実行を可能にする。その上、方法は、同業者費用（industrial rate）での実行に適している。
【図１】
0006770866.tif 000063
【図２】
0006770866.tif 000064
【図３】
0006770866.tif 000065
【図４】
0006770866.tif 000066
【図５】
0006770866.tif 000067
【図６】
0006770866.tif 000068
【図７】
0006770866.tif 000069
【図８】
0006770866.tif 000070
【図９】
0006770866.tif 000071
【図１０】
0006770866.tif 000072