(57)【特許請求の範囲】
【請求項１】
被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素を含む検出器と、
前記複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いた演算を行うことにより、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理を行う情報処理部と、
を有する放射線撮像システム。
【請求項２】
前記情報処理部は、推定された前記物質の厚みと密度に関する情報に基づいて、前記物質の画素値を生成する画素値生成部を更に有することを特徴とする請求項１に記載の放射線撮像システム。
【請求項３】
前記情報処理部は、
前記複数の画素のうちの任意の画素の画素値から前記任意の画素の画素値の平均を算出する第１算出部と、
前記任意の画素の画素値から算出された前記任意の画素の画素値の分散と、前記第１算出部によって算出された前記任意の画素の画素値の平均と、に基づいて、前記任意の画素の放射線量子のエネルギーの平均値を算出する第２算出部と、
前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数を用いた演算を行うことにより、前記任意の画素の画素値の平均を推定する第１推定部と、
前記２つの物質の減弱係数を用いた演算を行うことにより、前記任意の画素の放射線量子のエネルギーの平均値を推定する第２推定部と、
前記第１推定部で推定された前記任意の画素の画素値の平均と、前記第２推定部で推定された前記任意の画素の放射線量子のエネルギーの平均値と、を用いた演算を行うことにより、前記被写体を構成する互いに異なる２つの物質のうち少なくとも一方の厚みと密度に関する情報を推定する第３推定部と、
を有することを特徴とする請求項２に記載の放射線撮像システム。
【請求項４】
前記第３推定部は、前記算出された前記任意の画素の画素値の平均と前記第１推定部で推定された前記任意の画素の画素値の平均とを近似計算し、且つ、前記算出された前記任意の画素の放射線量子のエネルギーの平均値と前記第２推定部で推定された前記任意の画素の放射線量子のエネルギーの平均値とを近似計算することにより、前記被写体を構成する物質の厚みと密度に関する情報を推定することを特徴とする請求項３に記載の放射線撮像システム。
【請求項５】
前記第３推定部は、反復法による数値的解法を用いて前記被写体を構成する物質の厚みと密度に関する情報を推定することを特徴とする請求項４に記載の放射線撮像システム。
【請求項６】
前記第１推定部は、前記２つの物質のうちの第１構成物質の厚さをｄ１、前記第１構成物質の密度をρ１、前記放射線のエネルギーがＥ［ｋｅｖ］の際の前記第１構成物質の質量減弱係数をμ（Ｅ）１、前記２つの物質のうちの第２構成物質の厚さをｄ２、前記第２構成物質の密度をρ２、前記放射線のエネルギーがＥ［ｋｅｖ］の際の前記第２構成物質の質量減弱係数をμ（Ｅ）２、エネルギーを画素値に変換するための任意の係数をβ、とすると、前記任意の画素の画素値の平均ｉＡｖｅを、
【数１】
0006676338.tif 000017
を用いた演算を行うことにより推定することを特徴とする請求項５に記載の放射線撮像システム。
【請求項７】
前記第２推定部は、前記任意の画素の放射線量子のエネルギーの平均値ｅＡｖｅを、
【数２】
0006676338.tif 000018
を用いた演算を行うことにより推定することを特徴とする請求項６に記載の放射線撮像システム。
【請求項８】
前記第３推定部は、前記第１構成物質の厚みと密度に関する情報を第１情報、前記第２構成物質の厚みと密度に関する情報を第２情報、ニュートン法にて反復ｎ回目の第１情報をρ１ｄ１ｎ、反復ｎ回目の第２情報をρ２ｄ２ｎ、とすると、
【数３】
0006676338.tif 000019
【数４】
0006676338.tif 000020
【数５】
0006676338.tif 000021
を用いた演算を行うことにより、反復して前記第１情報及び前記第２情報を推定することを特徴とする請求項７に記載の放射線撮像システム。
【請求項９】
前記第３推定部は、反復ｍ回目の第１情報ρ１ｄ１ｍ、及び、反復ｍ回目の第２情報をρ２ｄ２ｍが収束しているか否かを判定する収束判定を行うことを特徴とする請求項８に記載の放射線撮像システム。
【請求項１０】
前記画素値生成部は、前記第１構成物質の厚みと密度に関する情報である第１情報に基づいて前記第１構成物質の画素値である第１画素値を生成する第１画素値生成部と、前記第２構成物質の厚みと密度に関する情報である第２情報に基づいて前記第２構成物質の画素値である第２画素値を生成する第２画素値生成部と、有することを特徴とする請求項９に記載の放射線撮像システム。
【請求項１１】
前記画素値生成部は、前記第１画素値をＩ１、前記第２画素値をＩ２とすると、
【数６】
0006676338.tif 000022
及び
【数７】
0006676338.tif 000023
の少なくとも一方を用いた演算を行うことにより、前記物質の画素値を生成することを特徴とする請求項１０に記載の放射線撮像システム。
【請求項１２】
前記画素値生成部は、前記第１画素値をＩ１、前記第２画素値をＩ２、前記放射線の実効エネルギーをＥｅｆｆ、とすると、
【数８】
0006676338.tif 000024
及び
【数９】
0006676338.tif 000025
の少なくとも一方を用いた演算を行うことにより、前記物質の画素値を生成することを特徴とする請求項１０に記載の放射線撮像システム。
【請求項１３】
前記画素値生成部は、任意の単色放射線エネルギーをＥｍｏｎｏ、任意のスペクトルを持ったエネルギーの画素値をＩ（Ｅｍｏｎｏ）と設定し、
【数１０】
0006676338.tif 000026
を用いた演算を行うことにより、前記任意のスペクトルを持ったエネルギーの画素値を生成することを特徴とする請求項１０に記載の放射線撮像システム。
【請求項１４】
前記画素値生成部は、前記第１画素値をＩ１、前記第２画素値をＩ２、前記放射線の実効エネルギーをＥｅｆｆ、前記算出された前記任意の画素の画素値の平均をＩＡｖｅ、とすると、
【数１１】
0006676338.tif 000027
及び
【数１２】
0006676338.tif 000028
の少なくとも一方を用いた演算を行うことにより、前記物質の画素値を生成することを特徴とする請求項１０に記載の放射線撮像システム。
【請求項１５】
前記算出された放射線量子のエネルギーの平均値は、前記任意の画素の画素値の時間的及び／或いは空間的な変化量を用いて算出されることを特徴とする請求項１から１４のいずれか１項に記載の放射線撮像システム。
【請求項１６】
被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いて、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理を行う、ことを特徴とする放射線画像の情報処理装置。
【請求項１７】
被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いて、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理を行う、ことを特徴とする放射線画像の情報処理方法。
【請求項１８】
被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いて、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理をコンピュータに行わせる、ことを特徴とする放射線画像の情報処理を実行するプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線撮像システム、放射線画像の情報処理装置、放射線画像の情報処理方法、及び、そのプログラムに関する。
【背景技術】
【０００２】
放射線（Ｘ線）による医療画像診断や非破壊検査に用いる撮像装置として、半導体材料によって形成された平面検出器（Ｆｌａｔ  Ｐａｎｅｌ  Ｄｅｔｅｃｔｏｒ、以下ＦＰＤ）を用いた放射線撮像装置が知られている。このような放射線撮像装置は、例えば医療画像診断において、静止画や動画などのデジタル撮像装置として用いられうる。
【０００３】
ＦＰＤとしては、例えば、積分型のセンサおよびフォトンカウンティング型のセンサがある。積分型のセンサは、放射線の入射により発生した電荷の総量を計測する。それに対し、フォトンカウンティング型のセンサは、入射した放射線のエネルギー（波長）を識別し、複数のエネルギーレベルの各々について放射線の検出回数をカウントする。即ち、フォトンカウンティング型のセンサは、エネルギー分解能を有するため、物質の弁別や仮想的に単一のエネルギーの放射線で撮影した場合の画像の生成、骨密度の測定などへの応用が期待できる。しかしながら、入射する放射線量子の個数は膨大であるため、これを個別に数えるには高い動作速度が必要となる。このため、大面積のＦＰＤでフォトンカウンティング型のセンサを実現することは困難であった。
【０００４】
そこで、特許文献１では、所定領域ごとに平均画像濃度情報と画像濃度の分散情報とを用いて、放射線量子の個数やエネルギーの平均値を推定することにより、エネルギー分解能を有する放射線撮像装置が提案されている。即ち、特許文献１には、所定領域ごとに平均画像濃度情報と画像濃度の分散情報とを用いて、放射線量子の個数やエネルギーの平均値を推定し、放射線量子の個数と放射線量子のエネルギーの平均値という２種類の画像情報を取得する情報処理方法が開示されている。特許文献１の方式では、フォトンカウンティング型のセンサと比較して低い動作速度であっても、エネルギー分解能を有するセンサが実現できる。
【０００５】
一方、特許文献２には、エネルギーサブトラクション法という技術が開示されている。このエネルギーサブトラクション法は、２種類のエネルギーの放射線で撮影し、それぞれのエネルギーで撮影された２つの画像を取得し、所望の演算を施した後に２つの画像を差分処理することで、減弱係数の異なる２つの物質を分離した画像を生成できる。これは、放射線が物質中を透過する際に減衰する程度を示す減弱係数が物質ごとに異なり、この減弱係数は放射線のエネルギーに依存するという現象を利用したものである。また、特許文献２では、同様の現象を利用した、骨密度を測定するための技術である、二重エネルギーＸ線吸収測定法（ＤＥＸＡ法）という技術が開示されている。しかしながら、特許文献２に記載されているエネルギーサブトラクション法やＤＥＸＡ法では、２種類のエネルギーでの２回の放射線撮影を必要とする。そのため、エネルギーを切り替えている際の被写体の動きによるアーチファクトの発生や、放射線エネルギーの高速切替を要する、といった課題がある。このような課題に対しては、１種類のエネルギーでの１回の放射線撮影で物質の弁別が可能な特許文献１の処理方法の方が優位である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００９−２８５３５６号公報
【特許文献２】特開２０１３−２３６９６２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１では、取得された放射線量子の個数と放射線量子のエネルギーの平均値という２種類の画像情報から、放射線画像を構成する２つの物質を分離しうる情報の取得方法は開示されていない。そこで本発明は、放射線のエネルギーを変化させた撮影を行わずとも、被写体の放射線画像を構成する２つの物質を分離しうる情報を取得することが可能な技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の放射線撮像システムは、被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素を含む検出器と、前記複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いた演算を行うことにより、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理を行う情報処理部と、を有する。本発明の放射線画像の情報処理装置は、被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いて、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理を行う。本発明の放射線画像の情報処理方法は、被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いて、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理を行う。本発明の放射線画像の情報処理を実行するプログラムは、被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素のうちの任意の画素の画素値から算出された前記任意の画素の画素値の平均と前記任意の画素の画素値から算出された前記任意の画素の画素値の分散とに基づいて算出された前記任意の画素の放射線量子のエネルギーの平均値と、前記被写体を構成する物質のうち互いに異なる２つの物質の減弱係数と、を用いて、前記被写体を構成する物質の厚みと密度に関する情報を推定する処理をコンピュータに行わせる。
【発明の効果】
【０００９】
本発明により、放射線のエネルギーを変化させた撮影を行わずとも、被写体の放射線画像を構成する２つの物質を分離しうる情報を取得することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】放射線撮像システムの機能構成例を示す模式図
【図２】処理フローを示すフローチャート
【図３】一例としての骨と脂肪の質量減弱係数
【図４】放射線撮像システムの機能構成例を示す模式図
【図５】処理フローを示すフローチャート
【図６】放射線撮像システムの模式的ブロック図
【発明を実施するための形態】
【００１１】
以下、本発明の実施形態について、添付の図面を参照して具体的に説明する。なお、本明細書では、放射性崩壊によって放出される粒子（光子を含む）の作るビームであるα線、β線、γ線などの他に、同程度以上のエネルギーを有するビーム、例えばＸ線や粒子線、宇宙線なども、放射線に含まれるものとする。
【００１２】
（第１の実施形態）
まず、図１、図２を用いて、第１の実施形態の放射線撮像システムに係る構成と処理フローを説明する。図１は、第１の実施形態の放射線撮像システムの機能構成例を示す模式図である。図２は、第１の実施形態の処理フローを示すフローチャートである。
【００１３】
放射線撮像システムは、放射線撮像装置１０、コンピュータ１３、放射線制御装置１２、及び、放射線発生装置１１を含み得る。放射線発生装置１０１は被写体に放射線を照射する。放射線撮像装置１０は、被写体を透過して入射された放射線に応じた画素値を取得するための複数の画素を含む検出器を含み、検出器によって被写体を通過して入射した放射線に応じた画素値が得られる。本発明の情報処理部及び／又は情報処理装置である、コンピュータ１３は、画素値から検出器に到達した放射線量子のエネルギーの平均値を推定し、被写体を構成する物質の密度と厚さに関する情報の推定を行う。この推定については、後で詳細に説明する。また、コンピュータ１３は、コンピュータ１３の制御卓（不図示）を介して撮影者（不図示）から入力された撮影情報に基づいて、放射線撮像装置１０及び放射線制御装置１２に制御信号を与える。放射線制御装置１２は、コンピュータ１３からの制御信号を受けて、放射線発生装置１１の放射線源（不図示）から放射線を出射する動作や照射野絞り機構（不図示）の動作の制御を行う。放射線制御装置１２によって制御された放射線発生装置１１から出射された放射線に応じて、放射線撮像装置１０の検出器は当該放射線に応じた画像信号を出力する。出力された画像信号は、信号処理部１０５によってオフセット補正等の画像処理がなされた後、コンピュータ１３に伝送される。ここで、伝送には、公知の無線通信や有線通信が適用され得る。伝送された画像信号は、制御コンピュータ１３によって必要な画像処理がなされた後、コンピュータ１３の表示部（不図示）に表示され得る。ここで、画素値は画像信号を構成するものである。
【００１４】
コンピュータ１３は、その機能構成として、第１算出部１３１、第２算出部１３２、第１推定部１３３、第２推定部１３４、第３推定部１３５を具備する。
【００１５】
Ｓ２０１において、被写体を介して放射線撮像装置に放射線を所定期間照射して複数枚のデジタル画像信号をコンピュータ１３が取得する。
【００１６】
次にＳ２０２及びＳ２０３における演算の原理を説明する。ここで、所定期間に照射される放射線は一定であり、被写体は動きがないものと仮定する。得られたデジタル画像信号のうち、任意の一画素を選択し、選択した画素から得られたデジタル信号（以下、画素値を称する）は、理想的には一定となるはずだが、実際には画素値の時系列にばらつきが生じる。このばらつきには、量子ノイズが含まれる。量子ノイズは、単位時間当たりの放射線量子の量子数（例えばＸ線フォトンの数）がばらつくことによって生じる。この放射線量指数のばらつきは、離散的な事象に対する単位時間当たりの生起確率として鑑みると、所与の時間間隔で発生する離散的な事象を数える特定の確率変数を持つ離散確率分布であるポアソン分布に従う。ポアソン分布では、定数λ＞０に対し、自然数を値にとる確率変数が所望の条件を満たすとき、この確率変数はパラメータλのポアソン分布に従うという。すなわち、画素値の平均値が同じであっても、エネルギーが大きい放射線量子で形成された画像のほうが、画素値の分散が大きくなる。これを利用して、Ｘ線フォトン等の放射線量子のエネルギーを推定することができる。
【００１７】
以下に、式を用いて放射線量子のエネルギーを推定する方法を説明する。まず、放射線撮像装置にＴ（Ｔは２以上の自然数）回の放射線照射を行い、放射線撮像装置からＴ枚のデジタル画像信号が取得されたものとする。ここで、ｔ（ｔは２以上Ｔ以下の自然数）枚目のデジタル画像信号のある画素の画素値をＩ（ｔ）とし、当該画素に到達して吸収された放射線量子の量子数の合計をＮ個、放射線量子のエネルギーの平均値をＥＡｖｅとすると、以下の式（１）が成り立つ。
【００１８】
ＥＡｖｅ×Ｎ＝ΣＩ（ｔ）  ・・・（１）
式（１）より、１枚のデジタル画像信号の当該画素に到達して吸収された放射線量子の量子数の相加平均をｎＡｖｅとすると、以下の式（２）で表される。
【００１９】
ｎＡｖｅ＝Ｎ／Ｔ＝ΣＩ（ｔ）／ＥＡｖｅ／Ｔ  ・・・（２）
また、式（１）より、１枚のデジタル画像信号の当該画素に到達して吸収された放射線量子の量子数の標本分散をｎＶａｒとすると、以下の式（３）で表される。
【００２０】
ｎＶａｒ＝Σ［｛Ｉ（ｔ）／ＥＡｖｅ−ｎＡｖｅ｝２］／Ｔ  ・・・（３）
ここで、ポアソン分布では、期待値及び分散がパラメータλに等しい。また、サンプル数が大きくなるほど、相加平均は期待値に近づき、標本分散は分散に近づく。そこで、サンプル数を十分大きく（好ましくは無限大）として、放射線量子の量子数の相加平均ｎＡｖｅと放射線量子の量子数の標本分散をｎＶａｒとが等しいものと近似すると、式（２）と式（３）が等しいという仮定により以下の式（４）が導かれる。
【００２１】
ＥＡｖｅ＝Σ｛Ｉ（ｔ）２｝／Σ｛Ｉ（ｔ）｝−Σ｛Ｉ（ｔ）｝／Ｔ  ・・・（４）
このようにして、任意のｔ枚目のデジタル画像信号のある画素の画素値Ｉ（ｔ）より、当該画素に到達して吸収された放射線量子のエネルギーの平均値ＥＡｖｅを推定して算出することができる。
【００２２】
また、画素値Ｉ（ｔ）の相加平均をＩＡｖｅとすると、画素値Ｉ（ｔ）の相加平均ＩＡｖｅは、放射線量子の量子数の相加平均ｎＡｖｅを用いて以下の式（５）で表される。
【００２３】
ＩＡｖｅ＝ｎＡｖｅ×ＥＡｖｅ  ・・・（５）
また、画素値の標本分散をＩＶａｒとすると、画素値の標本分散ＩＶａｒは放射線量子の量子数の標本分散をｎＶａｒより以下の式（６）で表される。
【００２４】
ＩＶａｒ＝ｎＶａｒ×ＥＡｖｅ２  ・・・（６）
従って、当該画素に到達して吸収された放射線量子のエネルギーの平均値Ｅは、以下の式（７）でも表される。
【００２５】
ＥＡｖｅ＝ＩＶａｒ  ／ＩＡｖｅ・・・（７）
Ｓ２０２において、第１算出部１３１は、下記式（８）より、任意の画素の画素値Ｉ（ｔ）を用いて、任意の画素の画素値の標本分散ＩＶａｒを算出する。なお、ここでは、分散として標本分散を用いたが、不偏分散を用いてもよい。また、画素値の平均として、画素値I（ｔ）の相加平均ＩＡｖｅを用いたが、本発明はそれに限定されるものではない。
【００２６】
【数１】
0006676338.tif 000002
次に、Ｓ２０３において、第２算出部１３２は、式（７）より算出された下記式（９）より、任意の画素の画素値の標本分散ＩＶａｒを用いて、任意の画素の放射線量子のエネルギーの平均値ＥＡｖｅを算出する。ここで、αは、画素値とエネルギーの単位を変換する任意の定数である。なお、ここでは、分散として標本分散を用いたが、不偏分散を用いてもよい。
【００２７】
【数２】
0006676338.tif 000003
以降のステップでは、Ｓ２０２及びＳ２０３で算出された任意の画素の放射線量子のエネルギーの平均値ＥＡｖｅと、任意の画素の画素値と、を用いた演算を行う。この演算を行うことによって、被写体を構成する物質の厚みと密度に関する情報を推定する。なお、ここでは説明を簡略化するために、被写体の例としての人体を構成する物質の中で骨（第１構成物質）以外の物質を脂肪（第２構成物質）としている。なお、第１構成物質と第２構成物質は互いに異なるものである。これは、骨以外の物質とは、脂肪、筋肉、臓器、水分などであるが、これらは骨と比較して近い質量減弱係数を持つと考えられるためである。図３に、一例として、骨と脂肪の質量減弱係数を示す。なお、以下では骨と脂肪の質量減弱係数を用いて説明を行うが、診断用途や被写体の構成物質に応じて、任意の２つの質量減弱係数を用い得る。ここで、放射線のエネルギーがＥ［ｋｅｖ］の際の骨の質量減弱係数をμ（Ｅ）１、骨の厚さをｄ１、骨の密度をρ１とする。また、放射線のエネルギーがＥ［ｋｅｖ］の際の脂肪の質量減弱係数をμ（Ｅ）２、脂肪の厚さをｄ２、脂肪の密度をρ２、被写体を通過する前の放射線の各エネルギーの放射線量子の量子数の割合（エネルギースペクトル）をｎ（Ｅ）とする。Ｓ２０４では、第１推定部１３３は、以下の式（１０）を用いた演算を行うにより任意の画素の画素値の相加平均ｉＡｖｅを推定する。なお、ここで、βはエネルギーを画素値に変換するための任意の係数である。なお、画素値の平均として、画素値I（ｔ）の相加平均ｉＡｖｅを推定したが、本発明はそれに限定されるものではない。
【００２８】
【数３】
0006676338.tif 000004
なお、ｎ（Ｅ）は市販のスペクトロメーター等を用いて予め測定しておけば良い。また、簡略的にエネルギースペクトルを導出できる簡易式が公表されているため、それを用いれば放射線を照射する際の管電圧や付加フィルター等の条件に基づいてエネルギースペクトルを算出することも可能である。また、βは、例えば放射線撮像装置に用いられるシンチレータの特性やシステムゲインなどの特性に基づいて決定すればよい。
【００２９】
次に、Ｓ２０５では、第２推定部１３４は、Ｓ２０４で推定された任意の画素の画素値の相加平均ｉＡｖｅを用いて、以下の式（１１）を用いた演算を行うにより任意の画素の任意の画素の放射線量子のエネルギーの平均値ｅＡｖｅを推定する。
【００３０】
【数４】
0006676338.tif 000005
以上により、取得された画素値Ｉ（ｔ）に基づく放射線量子のエネルギーの平均値ＥＡｖｅと、取得された画素値Ｉ（ｔ）に基づく画素値の相加平均ＩＡｖｅと、が得られる。また、以上により、推定された画素値の相加平均ｉＡｖｅと、推定された放射線量子のエネルギーの平均値ｅＡｖｅと、が得られる。そして、Ｓ２０６で、第３推定部１３５は、Ｓ２０５までに得られた各値に基づいて、第１構成物質（骨）の厚さと密度に関する第１情報ρ１ｄ１、及び、第２構成物質（脂肪）の厚さと密度に関する第２情報ρ２ｄ２を推定する。ここでは、第３推定部１３５は、取得された各値と推定された各値とを近似計算することにより、第１情報ρ１ｄ１、及び、第２情報ρ２ｄ２を推定する。すなわち、第３推定部１３５は、下記式（１２）の非線形連立方程式を解く演算を行うことで、第１情報ρ１ｄ１、及び、第２情報ρ２ｄ２を推定する。ここで、第１情報は、第１構成物質（骨）の厚さと密度を乗算した値であり、第２情報は、第２構成物質（脂肪）の厚さと密度を乗算した値である。なお、第１情報ρ１ｄ１の値の単位はｇ／ｃｍ２となる。これは、１ｃｍ２あたりの骨量、つまり、骨密度として利用することが可能である。
【００３１】
【数５】
0006676338.tif 000006
なお、式（１２）に示す方程式系は非線形のため、代数的解法で解くことは困難である。そこで、本実施形態では、第３推定部１３５がニュートン法等の反復法による数値的解法を用いて解く方法を示す。
【００３２】
まず、第３推定部１３５は、第１情報ρ１ｄ１の初期値、及び、第２情報ρ２ｄ２の初期値を定める。ここで、それぞれの物質の厚さと密度がおおよそわかっている場合には、その値を初期値にすると反復は早く収束するが、それぞれの物質の厚さと密度がわからない場合には、０等の適当な値を入れても演算は可能である。ニュートン法にて反復ｍ回目の第１情報をρ１ｄ１ｍ、反復ｍ回目の第２情報をρ２ｄ２ｍ、とすると、下記式（１３）を用いた演算を行うにより反復して求めていくことができる。
【００３３】
【数６】
0006676338.tif 000007
なお、ここで、ｉＡｖｅｍは、反復ｍ回目の推定された画素値の相加平均ｉＡｖｅで、ｅＡｖｅｍは、反復ｍ回目の推定された放射線量子のエネルギーの平均値であり、式（１０）及び式（１１）によりそれぞれ、以下の式（１４）及び式（１５）で示される。
【００３４】
【数７】
0006676338.tif 000008
【００３５】
【数８】
0006676338.tif 000009
そして、Ｓ２０７で、第３推定部１３５は、反復ｍ回目の第１情報ρ１ｄ１ｍ、及び、反復ｍ回目の第２情報をρ２ｄ２ｍが収束しているか否かを判定する収束判定を行う。収束していないと判定された場合には、Ｓ２０４に戻り、Ｓ２０４〜Ｓ２０６の各ステップを反復する。収束していると判定された場合には、第３推定部１３５は、反復ｍ回目の第１情報ρ１ｄ１ｍ、及び、反復ｍ回目の第２情報ρ２ｄ２ｍを出力し、処理が終了される。なお、収束判定に際して、反復ｍ回目の第１情報ρ１ｄ１ｍと反復ｍ−１回目の第１情報ρ１ｄ１ｍ−１の差分を予め定めた閾値と比較することによって判定してもよい。また、反復ｍ回目の第２情報ρ２ｄ２ｍと反復ｍ−１回目の第２情報ρ２ｄ２ｍ−１の差分を予め定めた閾値と比較することによって判定してもよい。また、予め定めた回数の反復がなされた時点で処理終了としてもよく、閾値を用いた判定と反復の回数とを用いた判定との組み合わせを用いてもよい。
【００３６】
以上のステップＳ２０１からＳ２０７までの処理を行うことにより、被写体を構成する２つの物質の厚さと密度に関する情報である、第１情報及び第２情報を推定することが可能となる。
【００３７】
本実施形態によれば、放射線のエネルギーを変化させた撮影を用いなくても被写体を構成する互いに異なる２つの物質のうちの少なくとも１つの物質の厚さと密度に関する情報を推定することが可能となる。そのため、一般撮影や透視撮影において、複雑な放射線撮影システムにも適用可能な、放射線画像を構成する２つの物質を分離しうる情報が取得できる。
【００３８】
なお、本実施形態では、被写体として人体を用い、人体における骨と脂肪についての情報を推定する方法を例としたが、本発明ではこれに限定されない。例えば、造影剤と人体についての情報を推定し、造影剤の分布の可視性を向上させ、造影剤使用量を低減する事が可能となる。また、ガイドワイヤーと人体についての情報を推定する事により、ガイドワイヤーの可視性を向上させ、手術中の患者の安全や医師の負担を低減する事が可能となる。このように、本発明は物質の種類に限定されることなく、任意の互いに異なる二種類の物質のうちの少なくとも一種類の物質についての厚さと密度に関する情報を得る事が可能となる。
【００３９】
また、本実施形態では、骨と脂肪のそれぞれの密度と厚さに関する情報を取得するとしたが、骨および脂肪の密度が別の方法で入手もしくは推定できる場合は、それぞれの厚さであるｄ１、ｄ２を求める方法としてもよい。この場合、質量減弱係数は密度を乗算する事によって、線減弱係数となる。つまり、本実施形態で用いた式の質量減弱係数μを線減弱係数ρとして、演算を行えばよく、本発明は、減弱係数を用いることで行い得る。
【００４０】
また、本実施形態では、時系列の複数の画素値を用いて標本分散と相加平均を算出し放射線量子のエネルギーの平均値を推定したが、本発明ではこれに限定されない。例えば、Ｘ軸を列数、Ｙ軸を行数とした２次元空間の配列位置を用いて２次元に配列された複数の画素のうちの任意の画素の画素値に対して、周辺の複数の画素の画素値を用いて任意の画素の画素値の標本分散と相加平均を算出し放射線量子のエネルギーの平均値を推定してもよい。このようにする事で、一枚の静止画像からエネルギー画像を算出する事が可能である。
【００４１】
また、本実施形態の画素値の放射線量子のエネルギーの平均値ＥＡｖｅは、任意の画素の画素値から算出された相加平均と標本分散を用いたが、本発明はそれに限定されるものではない。以下に説明するように、任意の画素の画素値に基づいて放射線量子のエネルギーの平均値が算出されればよく、例えば、任意の画素の画素値の時間的及び／或いは空間的な変化量を用いて算出してもよい。
【００４２】
実際の放射線撮影においては、複数枚の画像を撮影（放射線撮像装置によってデジタル画像信号を複数回取得）している間に、被写体が動いてしまうことがある。例えば、心臓などの動く臓器を撮影する場合や、手術中に透視撮影を行う場合などである。被写体に動きがあると、放射線撮像装置が複数枚の画像データを出力する間に、ある画素に到達する放射線量子の量子数の例であるＸ線フォトンの数が変化することとなる。すなわち、ポアソン分布のパラメータλが変化してしまうこととなる。そのため、エネルギーの平均値を用いて生成した画像にアーチファクトが発生することとなり、診断性能を低下させてしまう。
【００４３】
そこで、任意の画素における放射線量子のエネルギーの平均値を推定するに際して、当該任意の画素の画素値の時間的及び／或いは空間的な変化量を用いて、任意の画素における放射線量子のエネルギーを推定することが望ましい。ここでいう時間的な画素値の変化量とは、前記任意の画素を特定したフレームと異なるフレームにおける前記任意の画素に対応した画素（同一位置或いはその近接位置の画素）の画素値と前記任意の画素の画素値の差分を意味する。また、空間的な画素値の変化量とは、前記任意の画素を特定したフレーム内において前記任意の画素に隣接或いは近接に位置した画素の画素値と前記任意の画素の画素値との差分を意味する。時間的及び空間的な変化量とは上記の混在を意味する。なお、前記任意の画素の画素値或いはそれと差分を取る画素の画素値は単一の場合、複数の場合の両方を含む。複数の画素の場合は、その画素値はそれらの画素の代表値（例えば、任意の画素値に対してリカーシブフィルタ処理された値や平均化処理等）とする。なお、時間的な変化量を求める際に用いる前記異なるフレームは、時間軸上に隣接するフレームであることが好ましいが、効果を損なわない程度に離れたものであってもよい。また、空間的な変化量における、任意の画素と異なる画素とは、隣接していることが好ましいが、効果を損なわない程度に互いに離れていても良い。このように、離れていても効果を損なわない範囲を、信号処理に用いられる全画素値のうちの一部である任意の範囲とする。この変化量を用いることにより、前記任意の画素のエネルギーの平均値の推定における誤差が抑制され得る。
【００４４】
より具体的に説明すれば、本発明においては、画素値の標本分散を任意の画素値の時間的及び／又は空間的な変化量の二乗の１／２倍とみなして任意の画素における放射線量子のエネルギーを近似する、という概念に基づいている。代表的には、任意の画素の画素値と、前記任意の画素を特定したフレームとは異なるフレーム内における前記任意の画素値と同一位置の画素の画素値の差分の二乗の１／２倍、とみなして、近似する、という概念に基づいている。そして、近似された任意の画素における放射線量子のエネルギーを用いて平均化し、任意の画素における放射線量子のエネルギーの平均値を算出する。任意の画素における放射線量子のエネルギーは、ポアソン分布のパラメータλと等しい期待値が変化した瞬間にしか、大きな誤差を生じない。従って、上記アーチファクトを抑制することが可能となる。
【００４５】
また、例えば、フォトンカウンティング型のセンサを検出器として用い、複数のエネルギーレベルの各々について任意の画素毎に検出回数をカウントすることによって得られた放射線量子のエネルギーから、放射線量子のエネルギーの平均値を推定してもよい。その場合、本発明の任意の画素の画素値は、フォトンカウンティング型のセンサの任意の画素の画素値も含み得る。
【００４６】
（第２の実施形態）
本実施形態では、第１の実施形態で取得した被写体を構成する２つの物質の厚さと密度に関する情報を用いて、２つの物質ごとの画素値の生成方法を説明する。以下に図４及び図５を用いて、一例として人体の骨（第１構成物質）とそれ以外の第２構成物質に分離した画像の生成方法を説明する。図４は、第２の実施形態の放射線撮像システムの機能構成例を示す模式図である。図５は、第２の実施形態の処理フローを示すフローチャートである。なお、第１の実施形態と同じ機能構成及び処理ステップについては同じ番号を付与し、詳細な説明は省略する。
【００４７】
図４の画素値生成部１４は、被写体を構成する物質の厚みと密度に関する情報（第１情報及び第２情報の少なくとも一方）に基づいて、それぞれの物質毎の物質画素値（第１画素値及び第２画素値を生成する。画素値生成部１４は、第１画素値生成部１４１と、第２画素値生成部１４２を有する。第１画素値生成部１４１では、骨の密度と厚さに関する情報（第１情報）に基づいて骨の画素値（第１画素値）を生成する。また、第２画素値生成部１４２は、脂肪の密度と厚さに関する情報（第２情報）に基づいて脂肪の画素値（第２画素値）を生成する。なお、ここでは説明を簡略化するために、第１の実施形態と同様に、被写体を構成する物質の中で骨以外の物質を可視化した画像を脂肪画像としている。これは、骨以外の物質とは、脂肪、筋肉、臓器、水分などであるが、これらは骨と比較して近い質量減弱係数を持つと考えられるためである。
【００４８】
図５のステップＳ２０８では、ステップＳ２０１からＳ２０７までの処理で得られた第１情報に基づいて、以下の式（１６）を用いた演算を行うことにより、骨の画素値（第１画素値）Ｉ１を算出する。
【００４９】
【数９】
0006676338.tif 000010
ステップＳ２０９では、ステップＳ２０１からＳ２０７までの処理で得られた第２情報に基づいて、以下の式（１７）を用いた演算を行うことにより、脂肪の画素値（第２画素値）Ｉ２を算出する。
【００５０】
【数１０】
0006676338.tif 000011
以上の処理を行うことにより、被写体を構成する２つの物質ごとの物質画素値を生成することが可能となる。
【００５１】
なお、本実施形態では、被写体として人体を用い、人体における骨と脂肪についての画素値を生成する方法を例としたが、本発明ではこれに限定されない。例えば、造影剤と人体についての情報を推定し、造影剤の分布の可視性を向上させ、造影剤使用量を低減する事が可能となる。また、造影剤の画素値を作成すれば、デジタルサブトラクションアンギオグラフィーの画像と同等の画像が造影剤注入前の画像の撮影を行う事なく作成する事ができるため、造影中の被写体と放射線撮像装置の相対的な位置変化などにも対応可能である。また、ガイドワイヤーと人体についての情報を推定する事により、ガイドワイヤーの可視性を向上させ、手術中の患者の安全や医師の負担を低減する事が可能となる。このように、本発明は物質の種類に限定されることなく、任意の二種類の物質についての画素値を生成する事が可能となる。また、二種類の画素値が必要でなく、骨の画素値のみ、また、造影剤の画素値のみを要求する場合では、ステップＳ２０８及びステップＳ２０９の一方を省略する事が可能である。なお、式（１６）及び式（１７）で得られた画素値は、放射線撮像装置で仮想的に骨のみ、又は、脂肪のみを撮影した画素値に相当する。
【００５２】
また、各画素値の生成は、照射放射線のエネルギースペクトルの平均（実効エネルギー）Ｅｅｆｆを用いて、以下の式（１８）及び式（１９）を用いた演算を行うことによって算出してもよい。
【００５３】
【数１１】
0006676338.tif 000012
【００５４】
【数１２】
0006676338.tif 000013
このように、照射放射線の実効エネルギーを用いることにより、計算を簡略化することが可能となる。
【００５５】
また、任意の単色放射線エネルギーＥｍｏｎｏを設定することで、以下の式（２０）を用いた演算を行うことにより、仮想的に任意のスペクトルを持ったエネルギーの画素値を生成することができる。
【００５６】
【数１３】
0006676338.tif 000014
例えば、μ１をヨード造影剤の質量減弱係数、μ２を人体の質量減弱係数、ρ１ｄ１をヨード造影剤の密度と厚さの積、ρ２ｄ２を人体の密度と厚さの積、とすると、ヨード造影剤とそれ以外の物質のコントラストを向上させたい場合を鑑みる。この場合、ヨード造影剤は３３．２ｋｅＶにＫ吸収端を持っているため、Ｅｍｏｎｏを３３．２ｋｅＶに設定して計算を行うとよい。
【００５７】
また、照射放射線の実効エネルギーを用い、取得された画素値Ｉ（ｔ）に基づく画素値の相加平均ＩＡｖｅを用いて、以下の式（２１）及び式（２２）を用いた演算を行うことにより、各画素値の生成を行ってもよい。
【００５８】
【数１４】
0006676338.tif 000015
【００５９】
【数１５】
0006676338.tif 000016
また、得られた情報μ１ρ１ｄ１もしくはρ１ｄ１を骨の画素値（第１画素値）として、μ２ρ２ｄ２もしくはρ２ｄ２を脂肪の画素値（第２画素値）としてもよい。μρｄを表示する事で、それぞれの物質による放射線の減弱の度合いを可視化できる。また、ρｄを表示する事でそれぞれの物質の面密度を可視化できる。例えばρ１ｄ１は骨密度の分布を可視化する事ができる。さらに、骨および脂肪の密度が別の方法で入手もしくは推定できる場合は、それぞれの厚さであるｄ１、ｄ２をそれぞれ骨の画素値（第１画素値）、脂肪の画素値（第２画素値）としてもよい。このようにすることで、各物質の厚さの二次元分布を可視化する事が可能となる。
【００６０】
以下に、本発明に使用する画素値を取得するのに好適な放射線撮像装置及び放射線撮像システムについて説明する。
【００６１】
まず、図６を用いて放射線撮像システムを説明する。図６は、放射線撮像システムの模式的ブロック図である。なお、図１及び図４で説明した構成と同じ構成については同じ番号を付与し、詳細な説明は省略する。
【００６２】
検出器１０１は、放射線又は光を電気信号に変換する画素を２次元行列状に複数備えた画素アレイ１０２と、画素アレイ１０２を駆動する駆動回路１０３と、駆動された画素アレイ１０２からの電気信号を画像信号として出力する出力回路１０４と、を含み得る。画素アレイ１０２は、入射された放射線に応じた画素値を取得するために放射線に応じた電気信号を出力する複数の画素が配列されており、好ましくは２次元行列状に複数の画素が配列されている。複数の画素はそれぞれ、光電変換素子と、画素回路部と、を含み得る。光電変換素子は、放射線を光に変換するシンチレータで変換された光を電荷に変換する素子であり、本実施形態では、光電変換素子として、シリコン基板等の半導体基板に設けられたフォトダイオードが用いられている。ただし、本発明はそれに限定するものではない。例えば、ガラス基板などの絶縁基板上に配置されたアモルファスシリコンの光電変換素子や、シンチレータを用いずに放射線を直接電荷に変換する変換素子を用いてよい。放射線撮像装置１０の制御部１０７は、制御コンピュータ１３からの制御信号を受けて、放射線撮像装置１０の各部の制御を行う。放射線制御装置１２によって制御された放射線発生装置１１から出射された放射線に応じて、放射線撮像装置１０の検出器１０１は当該放射線に応じた画像信号を出力する。出力された画像信号は、信号処理部１０５によってオフセット補正等の画像処理がなされた後、コンピュータ１３に伝送される。ここで、伝送には、公知の無線通信や有線通信が適用され得る。伝送された画像信号は、コンピュータ１３によって必要な画像処理がなされた後、コンピュータ１３の表示部（不図示）に表示され得る。
【００６３】
なお、本発明は、例えば、システム、装置、方法、プログラム若しくは記憶媒体等としての実施態様を採ることもできる。具体的には、複数の機器から構成されるシステムに適用してもよいし、また、一つの機器からなる装置に適用してもよい。また、上述した処理は、好適にはプログラムを用いて処理されることが望ましいが、処理の全て、あるいは一部を、回路を用いて実施してもよい。また、コンピュータ１３ではなく、信号処理部１０５によって行われてもよく、信号処理部１０５とコンピュータ１３の両方を活用して行ってもよい。すなわち、本発明の情報処理部及び／又は情報処理装置は、信号処理部１０５、コンピュータ１３、及び、信号処理部１０５及びコンピュータ１３、の少なくとも１つが相当する。
【００６４】
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給する。そして、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ、ＧＰＵ等）がプログラムを読み出して実行する処理である。
【符号の説明】
【００６５】
１０１  検出器
１０５  信号処理部
１３  コンピュータ
【図１】
0006676338.tif 000029
【図２】
0006676338.tif 000030
【図３】
0006676338.tif 000031
【図４】
0006676338.tif 000032
【図５】
0006676338.tif 000033
【図６】
0006676338.tif 000034