(57)【要約】【課題】Ｘ線診断装置及び画像処理装置において、目的部位の診断に適したＤＥ画像を表示すること。【解決手段】Ｘ線診断装置１０は、被検体に対して、第１線質を用いたＸ線撮影と、第１線質より小さい第２線質を用いたＸ線撮影とを行なうように制御するＸ線撮影制御部３１と、第１線質を用いたＸ線撮影によって生成される第１画像と、第２線質を用いたＸ線撮影によって生成される第２画像とを記憶する画像メモリ１３と、第１画像内の代表値群としての画素値群を取得すると共に、第２画像内の代表値群としての画素値群を取得する代表値演算部３２２と、第１画像内の代表値群間の画素値の幅を第１幅として算出すると共に、第２画像内の代表値群間の画素値の幅を第２幅として算出する幅演算部３２３と、第１幅及び第２幅の比を算出する比演算部３２４と、比を重み係数とし、第１画像及び第２画像を基にデュアルエナジー画像を生成する係数処理部３３と、を有する。【選択図】  図４
2010240054.tif 000002
【特許請求の範囲】
【請求項１】
被検体に対して、第１線質を用いたＸ線撮影と、前記第１線質より小さい第２線質を用いたＸ線撮影とを行なうように制御する制御部と、
前記第１線質を用いたＸ線撮影によって生成される第１画像と、前記第２線質を用いたＸ線撮影によって生成される第２画像とを記憶する記憶手段と、
前記第１画像における観察対象組織の画素値の分布範囲を表した第１範囲と、前記第２画像における観察対象組織の画素値の分布範囲を表した第２範囲とを求める算出手段と、
前記第１範囲及び前記第２範囲に基づいて重み係数を求め、前記重み係数を用いて前記第１画像及び前記第２画像の差分画像を生成する画像生成手段と、
を有することを特徴とするＸ線診断装置。
【請求項２】
前記画像生成手段は、
前記第１画像内の代表値群としての画素値群を取得すると共に、前記第２画像内の代表値群としての画素値群を取得する取得手段と、
前記第１画像内の前記代表値群間の画素値の幅を第１幅として算出すると共に、前記第２画像内の前記代表値群間の画素値の幅を第２幅として算出する幅算出手段と、
前記第１幅及び前記第２幅の比を前記重み係数として算出する比算出手段と、
を有することを特徴とする請求項１に記載のＸ線診断装置。
【請求項３】
前記取得手段は、前記第１画像を構成する画素値群に基づくヒストグラムと、前記第２画像を構成する画素値群に基づくヒストグラムとを基に、前記代表値群をそれぞれ取得する構成とすることを特徴とする請求項２に記載のＸ線診断装置。
【請求項４】
前記取得手段は、前記第１画像のヒストグラムを基に前記代表値群としての直接線のピークにおける画素値及び前記被検体内の最高画素値をそれぞれ取得し、前記第２画像のヒストグラムを基に前記代表値群としての直接線のピークにおける画素値及び前記被検体内の最高画素値をそれぞれ取得し、
前記幅算出手段は、前記第１画像に基づく前記直接線の画素値及び前記最高画素値の幅と、前記第２画像に基づく前記直接線の画素値及び前記最高画素値の幅とをそれぞれ求める構成とすることを特徴とする請求項３に記載のＸ線診断装置。
【請求項５】
前記取得手段は、前記第１画像のヒストグラムを基に前記代表値群としての前記被検体内の最高画素値及び最低画素値をそれぞれ取得し、前記第２画像のヒストグラムを基に前記代表値群としての前記被検体内の最高画素値及び最低画素値をそれぞれ取得し、
前記幅算出手段は、前記第１画像に基づく前記最高画素値及び前記最低画素値の幅と、前記第２画像に基づく前記最高画素値及び前記最低画素値の幅とをそれぞれ求める構成とすることを特徴とする請求項３に記載のＸ線診断装置。
【請求項６】
前記第１画像と、前記第２画像とで同一であるとみなされる部分にそれぞれ関心領域を設定する関心領域設定手段をさらに有し、
前記画像生成手段は、前記関心領域設定手段によって設定された前記第１画像の前記関心領域内の画素値の平均値を前記代表値群の一つとして設定し、前記第２画像の前記関心領域内の画素値の平均値を前記代表値群の一つとして設定する構成とすることを特徴とする請求項２乃至５のうちいずれか一項に記載のＸ線診断装置。
【請求項７】
関心領域を入力する関心領域入力手段をさらに有し、
前記画像生成手段は、前記関心領域入力手段を介して入力された、前記第１画像の前記関心領域内の画素値の平均値を前記代表値群の一つとして設定し、前記第２画像の前記関心領域内の画素値の平均値を前記代表値群の一つとして設定する構成とすることを特徴とする請求項２乃至５のうちいずれか一項に記載のＸ線診断装置。
【請求項８】
被検体に対する第１線質を用いたＸ線撮影によって生成される第１画像と、前記第１線質より小さい第２線質を用いたＸ線撮影によって生成される第２画像とを記憶する記憶手段と、
前記第１画像における観察対象組織の画素値の分布範囲を表した第１範囲と、前記第２画像における観察対象組織の画素値の分布範囲を表した第２範囲とを求める算出手段と、
前記第１範囲及び前記第２範囲に基づいて重み係数を求め、前記重み係数を用いて前記第１画像及び前記第２画像の差分画像を生成する画像生成手段と、
を有することを特徴とする画像処理装置。
【請求項９】
前記画像生成手段は、
前記第１画像内の代表値群としての画素値群を取得すると共に、前記第２画像内の代表値群としての画素値群を取得する取得手段と、
前記第１画像内の前記代表値群間の画素値の幅を第１幅として算出すると共に、前記第２画像内の前記代表値群間の画素値の幅を第２幅として算出する幅算出手段と、
前記第１幅及び前記第２幅の比を前記重み係数として算出する比算出手段と、
を有することを特徴とする請求項８に記載の画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デュアルエナジーサブトラクション方式を用いたＸ線撮影に使用されるＸ線診断装置及び画像処理装置に関する。
【背景技術】
【０００２】
Ｘ線診断装置によるＸ線撮影法の一つであるデジタルサブトラクション法の中に、同一撮影位置で管電圧を変えて線質の異なるＸ線による２枚の画像を撮影し、各画像の濃度値に重み係数（ｗ値）を掛けて（重み係数を加えたり、下駄を履かせたりする等も可能）、対応する画素毎に引算を行なうことにより、撮影画像の中から例えばＸ線吸収の大きい骨等を消去した画像や、逆にＸ線吸収の小さい軟組織等を消去した画像を得るデュアルエナジー（ＤＥ）サブトラクション法と呼ばれるものがある。
【０００３】
Ｘ線診断装置は、重み係数を変化させることで、軟部組織を除去した骨だけの画像や、骨を除去した軟部組織だけのＤＥ画像を取得することができる。Ｘ線診断装置では、照射Ｘ線のエネルギーに対応する重み係数がテーブルとして予め決定されている。
【０００４】
デュアルエナジーサブトラクション方式を用いたＸ線診断装置として、３極Ｘ線管を用いて、フィラメント電流を低管電圧撮影時に許容最大電流を流す電流値に設定し、高管電圧撮影時にはグリッドに負の電圧を印加して、高管電圧に対する許容最大電流以下になるように制御し、低管電圧撮影時にはグリッドに印加する電圧を０にして低管電圧に対する許容最大電流を流すようにする装置がある（例えば、特許文献１参照。）。
【０００５】
また、デュアルエナジーサブトラクション方式を用いたＸ線診断装置として、特許文献２のような装置がある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００８−７３１１５号公報
【特許文献２】特開２００２−３５９７８１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
予め決定された重み係数のテーブルを用いてＤＥ画像を生成すると、被検体の体厚等の個人差によって、輝度のばらついたＤＥ画像が生成されるので、診断者の技量によって診断がばらつく可能性があった。
【０００８】
本発明は、上述のような事情を考慮してなされたもので、目的部位の診断に適したＤＥ画像を表示できるＸ線診断装置及び画像処理方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係るＸ線診断装置は、上述した課題を解決するために、被検体に対して、第１線質を用いたＸ線撮影と、前記第１線質より小さい第２線質を用いたＸ線撮影とを行なうように制御する制御部と、前記第１線質を用いたＸ線撮影によって生成される第１画像と、前記第２線質を用いたＸ線撮影によって生成される第２画像とを記憶する記憶手段と、前記第１画像における観察対象組織の画素値の分布範囲を表した第１範囲と、前記第２画像における観察対象組織の画素値の分布範囲を表した第２範囲とを求める算出手段と、前記第１範囲及び前記第２範囲に基づいて重み係数を求め、前記重み係数を用いて前記第１画像及び前記第２画像の差分画像を生成する画像生成手段と、を有する。
【００１０】
本発明に係る画像処理装置は、上述した課題を解決するために、被検体に対する第１線質を用いたＸ線撮影によって生成される第１画像と、前記第１線質より小さい第２線質を用いたＸ線撮影によって生成される第２画像とを記憶する記憶手段と、前記第１画像における観察対象組織の画素値の分布範囲を表した第１範囲と、前記第２画像における観察対象組織の画素値の分布範囲を表した第２範囲とを求める算出手段と、前記第１範囲及び前記第２範囲に基づいて重み係数を求め、前記重み係数を用いて前記第１画像及び前記第２画像の差分画像を生成する画像生成手段と、を有する。
【発明の効果】
【００１１】
本発明に係るＸ線診断装置及び画像処理装置によると、目的部位の診断に適したＤＥ画像を表示できる。
【図面の簡単な説明】
【００１２】
【図１】ＤＥ画像の生成方法を示す概念図。
【図２】本実施形態のＸ線診断装置の全体構成を示す概略図。
【図３】画像処理装置の構成を示す図。
【図４】第１実施形態のＸ線診断装置の機能を示すブロック図。
【図５】高エネルギー画像、低エネルギー画像のそれぞれについての画素値のヒストグラムを示す図。
【図６】低エネルギー画像のヒストグラム上の代表値としての直線性のピークにおける画素値の求め方を説明するための図。
【図７】低エネルギー画像のヒストグラム上の代表値としての胸部の最高・最低画素値の求め方を説明するための図。
【図８】高エネルギー画像に基づく代表値群間の幅の求め方を説明するための図。
【図９】低エネルギー画像に基づく代表値群間の幅の求め方を説明するための図。
【図１０】第２実施形態のＸ線診断装置の機能を示すブロック図。
【図１１】第３実施形態のＸ線診断装置の機能を示すブロック図。
【発明を実施するための形態】
【００１３】
本発明に係るＸ線診断装置及び画像処理装置の実施形態について、添付図面を参照して説明する。
【００１４】
まず、ＤＥ（デュアルエナジー）画像の生成方法について説明する。
【００１５】
物質（透過係数μ、厚さｘ）を透過した透過Ｘ線のエネルギー（線質、透過力）Ｉは、Ｘ線管から照射された単色Ｘ線のエネルギーＩ０を用いて次の式のように表せる。
【数１】
2010240054.tif 000003
【００１６】
また、透過Ｘ線のエネルギーＩに比例する値をＬＯＧ変換した信号ｓは、次の式（１）のように表せる。
【数２】
2010240054.tif 000004
【００１７】
【数３】
2010240054.tif 000005
【００１８】
さらに、高エネルギーＩ０Ｈの単色Ｘ線が軟部組織（透過係数μａＨ、厚さｘａ）と骨（透過係数μｂＨ、厚さｘｂ）とが混在する被検体を透過してＬＯＧ変換された信号ＳＨと、低エネルギーＩ０Ｌの単色Ｘ線を軟部組織（透過係数μａＬ、厚さｘａ）と骨（透過係数μｂ、厚さｘｂ）とが混在する被検体を透過してＬＯＧ変換された信号ＳＬとは、次の式（２）及び（３）のように表せる。なお、透過係数は入射されるＸ線のエネルギーに依存する。
【数４】
2010240054.tif 000006
【００１９】
重み係数ｗをつけて、上記式（２）及び（３）の差分をとると次の式（４）のように表せる。
【数５】
2010240054.tif 000007
【００２０】
ここで、次の式（５）に示すような重み係数ｗが求められれば、上記式（４）を基に骨（ｘｂ）がキャンセルされて軟部組織のＤＥ画像が得られる一方、次の式（６）に示すような重み係数ｗが求められれば、上記式（４）を基に軟部組織（ｘａ）がキャンセルされて骨のＤＥ画像が得られる。
【数６】
2010240054.tif 000008
【００２１】
図１は、ＤＥ画像の生成方法を示す概念図である。
【００２２】
図１に示すように、高エネルギーのＸ線（高い管電圧によるＸ線）を用いた撮影によって生成された高エネルギー画像ＩＨと、低エネルギーのＸ線（低い管電圧によるＸ線）を用いた撮影によって生成された低エネルギー画像ＩＬとのうち一方、例えば低エネルギー画像ＩＬのデータに重み係数ｗを乗じて差分処理することで、差分画像としてのＤＥ画像ＩＤＥが生成される。
【００２３】
本実施形態のＸ線診断装置は、上記式（５）のμｂＨとμｂＬとの比を精度よく求めて重み係数ｗを求め、骨がキャンセルされた軟部組織のＤＥ画像ＩＤＥを精度よく生成するものである。また、本実施形態のＸ線診断装置は、上記式（６）のμａＨとμａＬとの比を精度よく求めて重み係数ｗを求め、軟部組織がキャンセルされた骨のＤＥ画像ＩＤＥを精度よく生成するものである。
【００２４】
図２は、第１実施形態のＸ線診断装置の全体構成を示す概略図である。
【００２５】
図２は、第１実施形態のＸ線診断装置１０を示す。そのＸ線診断装置１０は、Ｘ線管１１、Ｘ線検出装置１２、画像メモリ１３、画像処理装置１４、コントローラ１５及び高電圧供給器１６を備える。
【００２６】
Ｘ線管１１は、高電圧供給器１６から高電圧電力の供給を受けて、この高電圧電力の条件に応じて被検体（患者Ｐの所定部位）に向かってＸ線を曝射する。Ｘ線管２１は、被検体を撮影するための撮影Ｘ線を曝射することができる。
【００２７】
Ｘ線検出装置１２は、Ｘ線管１１の出射側に設けられ、被検体を透過したＸ線を検出する。Ｘ線検出装置１２は、平面検出器（ＦＰＤ：ｆｌａｔ  ｐａｎｅｌ  ｄｅｔｅｃｔｏｒ）１２ａ、ゲートドライバ１２ｂ及びデータ生成部１２ｃを備える。なお、Ｘ線検出装置１２は、ＣＲ（ｃｏｍｐｕｔｅｄ  ｒａｄｉｏｇｒａｐｈｙ）系や、Ｉ．Ｉ．（ｉｍａｇｅ  ｉｎｔｅｎｓｉｆｉｅｒ）−ＴＶ系であってもよい。
【００２８】
平面検出器１２ａは、被検体を透過したＸ線を電荷に変換して蓄積する。ゲートドライバ１２ｂは、平面検出器１２ａに蓄積された電荷をＸ線画像信号として読み出す。データ生成部１２ｃは、読み出された電荷を画像データに変換する。データ生成部１２ｃは、平面検出器１２ａから読み出された電荷を電圧に変換する電荷・電圧変換器（図示しない）と、電荷・電圧変換器の出力をデジタル信号に変換するＡ／Ｄ変換器（図示しない）と、平面検出器１２ａからライン単位でパラレルに読み出される画像信号をシリアルな信号に変換するパラレル・シリアル変換器（図示しない）とを備える。
【００２９】
画像メモリ１３は、データ生成部１２ｃからのＸ線画像を記憶する。
【００３０】
図３は、画像処理装置１４の構成を示す図である。
【００３１】
画像処理装置１４は、ＣＰＵ２１、包括メモリ２２、ＨＤＤ（ｈａｒｄ  ｄｉｓｋ  ｄｒｉｖｅ）２３、通信制御装置２４、ディスプレイ２５及び入力装置２６を有する。
【００３２】
ＣＰＵ２１は、半導体で構成された電子回路が複数の端子を持つパッケージに封入されている集積回路（ＬＳＩ）の構成をもつ制御装置である。ＣＰＵ２１は、医師及び技師等のオペレータによって入力装置２６が操作等されることにより指令が入力されると、包括メモリ２２に記憶しているプログラムを実行する。又は、ＣＰＵ２１は、ＨＤＤ２３のＨＤ（ｈａｒｄ  ｄｉｓｋ）に記憶しているプログラム、ネットワークＮから転送され通信制御装置２４で受信されてＨＤにインストールされたプログラム、又は記録媒体用のドライブ（図示しない）に装着された記録媒体から読み出されてＨＤにインストールされたプログラムを、包括メモリ２２にロードして実行する。
【００３３】
包括メモリ２２は、ＲＯＭ（ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）及びＲＡＭ（ｒａｎｄｏｍ  ａｃｃｅｓｓ  ｍｅｍｏｒｙ）等の要素を兼ね備える構成をもつ記憶装置である。包括メモリ２２は、ＩＰＬ（ｉｎｉｔｉａｌ  ｐｒｏｇｒａｍ  ｌｏａｄｉｎｇ）、ＢＩＯＳ（ｂａｓｉｃ  ｉｎｐｕｔ／ｏｕｔｐｕｔ  ｓｙｓｔｅｍ）及びデータを記憶したり、ＣＰＵ２１のワークメモリやデータの一時的な記憶に用いたりする記憶装置である。
【００３４】
ＨＤＤ２３は、磁性体を塗布又は蒸着した金属のディスクであるＨＤを着脱不能で内蔵する構成をもつ読み取り装置である。ＨＤは、画像処理装置１４にインストールされたプログラム（アプリケーションプログラムの他、ＯＳ（ｏｐｅｒａｔｉｎｇ  ｓｙｓｔｅｍ）等も含まれる）や、データを記憶する。また、ＯＳに、ユーザに対する情報の表示にグラフィックを多用し、基礎的な操作を入力装置２６によって行なうことができるＧＵＩ（ｇｒａｐｈｉｃａｌ  ｕｓｅｒ  ｉｎｔｅｒｆａｃｅ）を提供させることもできる。
【００３５】
通信制御装置２４は、各規格に応じた通信制御を行なう。通信制御装置２４は、ネットワークＮに接続することができる機能を有しており、これにより、Ｘ線診断装置１０は、通信制御装置２４からネットワークＮ網に接続することができる。
【００３６】
ディスプレイ２５は、液晶ディスプレイやＣＲＴ（ｃａｔｈｏｄｅ  ｒａｙ  ｔｕｂｅ）等によって構成される。ディスプレイ２５は、ビデオ信号に基づいて、ＤＥ画像を種々のパラメータの文字情報や目盛等と共に表示する機能を有する。
【００３７】
入力装置２６としては、オペレータによって操作が可能なキーボード及びマウス等によって構成される。入力装置２６による操作に従った入力信号は、ＣＰＵ２１に送られる。
【００３８】
図２に示すコントローラ１５は、図示しないＣＰＵ及びメモリを含んでいる。コントローラ１５は、画像処理装置１４のＣＰＵ２１からの指示に従って、高電圧供給器１６の動作を制御する。
【００３９】
高電圧供給器１６は、コントローラ１５による制御の下、Ｘ線管１１に高電圧電力を供給する。
【００４０】
図４は、第１実施形態のＸ線診断装置の機能を示すブロック図である。
【００４１】
図３に示すＣＰＵ２１がプログラムを実行することによって、Ｘ線診断装置１０は、Ｘ線撮影制御部３１、係数設定部３２及び係数処理部３３として機能する。なお、各部３１乃至３３は、Ｘ線診断装置１０の機能として備えられるものとして説明するが、Ｘ線診断装置１０にハードウェアとして備えられるものであってもよい。
【００４２】
Ｘ線撮影制御部３１は、同一の被検体に対して、高エネルギーＸ線及び低エネルギーＸ線を用いて、２枚のＸ線画像（高エネルギー画像ＩＨ、低エネルギー画像ＩＬ）を撮影するようにコントローラ１５を制御する機能を有する。
【００４３】
係数設定部３２は、Ｘ線撮影制御部３１による制御によって画像メモリ１３に記憶された高エネルギー画像ＩＨ、低エネルギー画像ＩＬを基に、上記式（４）に付加された重み係数ｗを設定する機能を有する。係数設定部３２は、ヒストグラム生成部３２１、代表値演算部３２２、幅演算部３２３及び比演算部３２４を有する。
【００４４】
ヒストグラム生成部３２１は、Ｘ線撮影制御部３１の制御によって生成される高エネルギー画像ＩＨ、低エネルギー画像ＩＬのそれぞれについての画素値のヒストグラム（図５に図示）を生成する機能を有する。なお、図５に示すヒストグラムについて、一般的に、出力される画像データの画素値は、線量に対してＬＯＧ変換されている。よって、これらのヒストグラムの横軸の画素値は、線量のＬＯＧに比例する値となっている。
【００４５】
図４に示す代表値演算部３２２は、高エネルギー画像ＩＨのヒストグラム上の代表的な画素値（代表値）群を求める一方、低エネルギー画像ＩＬのヒストグラム上の代表値群を求める機能を有する。代表値としては、ヒストグラム上の直接線のピークにおける画素値、被検体内の最高画素値、及び被検体内の最低画素値等が挙げられる。
【００４６】
図６，図７は、低エネルギー画像ＩＬのヒストグラム上の代表値群の求め方の一例としての胸部のヒストグラムを示す図である。図６は、低エネルギー画像ＩＬのヒストグラム上の代表値としての直線性のピークにおける画素値の求め方を説明するための図である。図７は、低エネルギー画像ＩＬのヒストグラム上の代表値としての胸部の最高・最低画素値の求め方を説明するための図である。
【００４７】
図６に示すように、ヒストグラムの高画素値側（右側）からみて閾値Ａ１を超える最初のピークにおける画素値や、最も頻度の高いピークにおける画素値を、直接線のピークにおける画素値とする。
【００４８】
図７に示すように、ヒストグラムの直接線を除いた部分が、胸部のヒストグラムと言える。そして、胸部のヒストグラムの高画素値側からみて最初に頻度が閾値Ａ２を超える頻度に相当する画素値を最高画素値とする。また、図７に示すように、胸部のヒストグラムの低画素値側（左側）からみて最初に閾値Ａ３を超える頻度に相当する画素値を最低画素値とする。
【００４９】
また、図４に示す幅演算部３２３は、高エネルギー画像ＩＨのヒストグラムに基づく代表値群間の幅と、低エネルギー画像ＩＬのヒストグラムに基づく代表値群間の幅とを求める機能を有する。
【００５０】
図８は、高エネルギー画像ＩＨに基づく代表値群間の幅の求め方を説明するための図である。図９は、低エネルギー画像ＩＬに基づく代表値群間の幅の求め方を説明するための図である。
【００５１】
図８，図９に示すように、胸部内の最高画素値は胸部内の比較的薄い軟部組織を表しているとみなせ、最低画素値は胸部内の比較的厚い軟部組織を表しているとみなせる。なお、ヒストグラムの最低画素値には、厳密には骨等の画素値も入ってしまうが、ヒストグラムの最低画素値を胸部の軟部組織の厚い部分と判断してもほとんど影響はない。
【００５２】
ここで、幅演算部３２３は、高エネルギー画像ＩＨのヒストグラムに基づく代表値群間の幅と、低エネルギー画像ＩＬのヒストグラムに基づく代表値群間の幅とを求める第１例として、高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値と比較的薄い軟部組織の画素値との幅と、低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値と比較的薄い軟部組織の画素値との幅を求める。
【００５３】
胸部の高エネルギー画像ＩＨ及び低エネルギー画像ＩＬの中で、比較的薄い軟部組織の画素値を胸部内の最高画素値とすると、幅演算部３２３は、各ヒストグラムにより、直接線の画素値と比較的薄い軟部組織の画素値との幅をそれぞれ求めることができる。一方、比較的薄い軟部組織（厚さｘａ１、透過係数μａ）とすると、上記式（１）により比較的薄い軟部組織を透過することによって減衰した信号ｓ´は次の式（７）のように表せる。
【数７】
2010240054.tif 000009
【００５４】
ここで、比較的薄い軟部組織により減衰した信号ｓ´、及び、直接線の画素値と比較的薄い軟部組織の画素値との幅は同等とみなすことができる。よって、幅演算部３２３は、上記式（７）の透過係数μａをμａＨとして、高エネルギー画像ＩＨのヒストグラムを基に直接線の画素値と比較的薄い軟部組織の画素値との幅（μａＨ×ｘａ１）を算出し、また、上記式（７）の透過係数μａをμａＬとして、低エネルギー画像ＩＬのヒストグラムを基に直接線の画素値と比較的薄い軟部組織の画素値との幅（μａＬ×ｘａ１）を算出する。
【００５５】
比演算部３２４は、幅演算部３２３によって算出された高エネルギー画像ＩＨのヒストグラムに基づく代表値群間の幅と、低エネルギー画像ＩＬのヒストグラムに基づく代表値群間の幅とを基に、比を求める機能を有する。例えば、比演算部３２４は、次の式（８）に示すように、高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値及び比較的薄い軟部組織の画素値の幅（μａＨ×ｘａ１）と、低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値及び比較的薄い軟部組織の画素値の幅（μａＬ×ｘａ１）との比を求める。
【数８】
2010240054.tif 000010
【００５６】
上記式（６），（８）によると、比演算部３２４は、高エネルギー画像ＩＨに基づく直接線の画素値及び比較的薄い軟部組織の画素値の幅（μａＨ×ｘａ１）と、低エネルギー画像ＩＬに基づく直接線の画素値及び比較的薄い軟部組織の画素値の幅（μａＬ×ｘａ１）との比を求めることで、軟部組織がキャンセルされた骨のＤＥ画像ＩＤＥを生成するための重み係数ｗを求めることができる。
【００５７】
また、幅演算部３２３は、高エネルギー画像ＩＨのヒストグラムに基づく代表値群間の幅と、低エネルギー画像ＩＬのヒストグラムに基づく代表値群間の幅とを求める第２例として、高エネルギー画像ＩＨのヒストグラムに基づく比較的薄い軟部組織の画素値と比較的厚い軟部組織の画素値との幅と、低エネルギー画像ＩＬのヒストグラムに基づく比較的薄い軟部組織の画素値と比較的厚い軟部組織の画素値との幅を求める。
【００５８】
胸部の高エネルギー画像ＩＨ及び低エネルギー画像ＩＬの中で、比較的厚い軟部組織の画素値を胸部での最低画素値とすると、幅演算部３２３は、各ヒストグラムにより、比較的薄い軟部組織と厚い軟部組織の画素値の幅を求めることができる。一方、比較的薄い軟部組織の厚さをｘａ１、比較的厚い軟部組織の厚さをｘａ２、軟部組織の透過係数をμａとすると、比較的薄い軟部組織を透過することによって減衰した信号と比較的厚い軟部組織を透過することによって減衰した信号との差Δｓ´は、次の式（９）のように表される。
【数９】
2010240054.tif 000011
【００５９】
なお、比較的厚い軟部組織と比較的薄い軟部組織とでは、厳密には同一の透過係数とならないこともある。しかし、重要なのは骨と軟部組織の透過係数の差であるので、軟部組織同士で同一の透過係数としても差し支えない。
【００６０】
ここで、比較的薄い軟部組織を透過することによって減衰した信号と比較的厚い軟部組織を透過することによって減衰した信号との差Δｓ´、及び、比較的薄い軟部組織の画素値と比較的厚い軟部組織の画素値との幅は同等とみなすことができる。よって、幅演算部３２３は、上記式（９）の透過係数μａをμａＨとして、高エネルギー画像ＩＨのヒストグラムを基に比較的薄い軟部組織の画素値と比較的厚い軟部組織の画素値との幅（μａＨ×（ｘａ２−ｘａ１））を算出し、また、上記式（９）の透過係数μａをμａＬとして、低エネルギー画像ＩＬのヒストグラムを基に比較的薄い軟部組織の画素値と比較的厚い軟部組織の画素値との幅（μａＬ×（ｘａ２−ｘａ１））を算出する。
【００６１】
比演算部３２４は、次の式（１０）に示すように、高エネルギー画像ＩＨのヒストグラムに基づく比較的薄い軟部組織の画素値及び比較的厚い軟部組織（観察対象組織）の分布範囲の画素値の幅（μａＨ×（ｘａ２−ｘａ１））と、低エネルギー画像ＩＬのヒストグラムに基づく比較的薄い軟部組織（観察対象組織）の分布範囲の画素値及び比較的厚い軟部組織の画素値の幅（μａＬ×（ｘａ２−ｘａ１））との比を求める。
【数１０】
2010240054.tif 000012
【００６２】
上記式（６），（１０）によると、比演算部３２４は、高エネルギー画像ＩＨに基づく比較的薄い軟部組織の画素値及び比較的厚い軟部組織の画素値の幅（μａＨ×（ｘａ２−ｘａ１））と、低エネルギー画像ＩＬに基づく比較的薄い軟部組織の画素値及び比較的厚い軟部組織の画素値の幅（μａＬ×（ｘａ２−ｘａ１））との比を求めることで、軟部組織がキャンセルされた骨のＤＥ画像ＩＤＥを生成するための重み係数ｗを求めることができる。
【００６３】
係数処理部３３は、低エネルギー画像ＩＬのデータに係数設定部３２によって設定された重み係数ｗを乗じて高エネルギー画像ＩＨから差分することで、ＤＥ画像ＩＤＥを生成する機能を有する。係数処理部３３によって生成される軟部組織がキャンセルされた骨のＤＥ画像ＩＤＥは、ディスプレイ２５に表示される。
【００６４】
なお、係数設定部３２は、軟部組織を透過することによって減衰した信号ｓ´を基に軟部組織がキャンセルされた骨のＤＥ画像ＩＤＥを生成するための重み係数ｗ（上記式（６））を生成するものとして説明したが、その場合に限定されるものではない。係数設定部３２は、骨を透過することによって減衰した信号ｓ´を基に骨がキャンセルされた軟部組織のＤＥ画像ＩＤＥを生成するための重み係数ｗ（上記式（５））を生成してもよい。
【００６５】
第１実施形態のＸ線診断装置１０によると、予め決定された重み係数ｗのテーブルを用いずに個々の被検体のＸ線画像を基に適宜重み係数ｗを求めるので、目的部位の診断に適したＤＥ画像ＩＤＥを表示できる。
【００６６】
図１０は、第２実施形態のＸ線診断装置の機能を示すブロック図である。なお、第２実施形態のＸ線診断装置の構成は、図２，図３に示す第１実施形態のＸ線診断装置の構成と同様であるので、説明を省略する。
【００６７】
図３に示すＣＰＵ２１がプログラムを実行することによって、第２実施形態のＸ線診断装置１０Ａは、Ｘ線撮影制御部３１、係数設定部３２Ａ、係数処理部３３及びＲＯＩ（ｒｅｇｉｏｎ  ｏｆ  ｉｎｔｅｒｅｓｔ）設定部３４として機能する。なお、各部３１乃至３４は、Ｘ線診断装置１０Ａの機能として備えられるものとして説明するが、Ｘ線診断装置１０Ａにハードウェアとして備えられるものであってもよい。係数設定部３２Ａは、ヒストグラム生成部３２１、代表値演算部３２２Ａ、幅演算部３２３Ａ及び比演算部３２４Ａを有する。
【００６８】
ＲＯＩ設定部３４は、高エネルギーのＸ線を用いた撮影によって生成された高エネルギー画像ＩＨと、低エネルギーのＸ線を用いた撮影によって生成された低エネルギー画像ＩＬとで同一であるとみなされる部分にそれぞれ関心領域としてのＲＯＩを設定する機能を有する。
【００６９】
係数設定部３２Ａは、ＲＯＩ設定部３４によって設定された高エネルギー画像ＩＨのＲＯＩ内の画素値の平均値等を代表値群の一つとして設定し、低エネルギー画像ＩＬのＲＯＩ内の画素値の平均値等を代表値群の一つとして設定する。
【００７０】
係数設定部３２Ａの代表値演算部３２２Ａは、高エネルギー画像ＩＨのＲＯＩ内の画素値の平均値と高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値とを演算し、また、低エネルギー画像ＩＬのＲＯＩ内の画素値の平均値と低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値とを演算する。
【００７１】
幅演算部３２３Ａは、高エネルギー画像ＩＨのＲＯＩ内の画素値の平均値と高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値との幅と、低エネルギー画像ＩＬのＲＯＩ内の画素値の平均値と低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値との幅を求める。
【００７２】
比演算部３２４Ａは、高エネルギー画像ＩＨのＲＯＩ内の画素値の平均値と高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値との幅と、低エネルギー画像ＩＬのＲＯＩ内の画素値の平均値と低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値との幅との比を求めることで、ＲＯＩで設定された部分がキャンセルされるような重み係数ｗを求めることができる。
【００７３】
なお、図１０に示すＸ線診断装置１０Ａにおいて、図４に示すＸ線診断装置１０と同一機能には同一符号を付して説明を省略する。
【００７４】
第２実施形態のＸ線診断装置１０Ａによると、予め決定された重み係数ｗのテーブルを用いずに個々の被検体のＸ線画像を基に適宜重み係数ｗを求めるので、目的部位の診断に適したＤＥ画像ＩＤＥを表示できる。
【００７５】
図１１は、第３実施形態のＸ線診断装置の機能を示すブロック図である。なお、第３実施形態のＸ線診断装置の構成は、図２，図３に示す第１実施形態のＸ線診断装置の構成と同様であるので、説明を省略する。
【００７６】
図３に示すＣＰＵ２１がプログラムを実行することによって、第３実施形態のＸ線診断装置１０Ｂは、Ｘ線撮影制御部３１、係数設定部３２Ｂ、係数処理部３３及びＲＯＩ設定部３４として機能する。なお、各部３１乃至３４は、Ｘ線診断装置１０Ｂの機能として備えられるものとして説明するが、Ｘ線診断装置１０Ｂにハードウェアとして備えられるものであってもよい。係数設定部３２Ｂは、ヒストグラム生成部３２１、代表値演算部３２２Ｂ、幅演算部３２３Ｂ及び比演算部３２４Ｂを有する。
【００７７】
係数設定部３２Ｂは、高エネルギー画像と低エネルギー画像とで単一であると思われる部分内に入力装置２６を介して手動ＲＯＩがそれぞれ入力されることで、高エネルギー画像ＩＨの手動ＲＯＩ内の画素値の平均値等を代表値群の一つとして設定し、低エネルギー画像ＩＬの手動ＲＯＩ内の画素値の平均値等を代表値群の一つとして設定する。
【００７８】
係数設定部３２Ｂの代表値演算部３２２Ｂは、高エネルギー画像ＩＨの手動ＲＯＩ内の画素値の平均値と高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値とを演算し、また、低エネルギー画像ＩＬの手動ＲＯＩ内の画素値の平均値と低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値とを演算する。
【００７９】
幅演算部３２３Ｂは、高エネルギー画像ＩＨの手動ＲＯＩ内の画素値の平均値と高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値との幅と、低エネルギー画像ＩＬの手動ＲＯＩ内の画素値の平均値と低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値との幅を求める。
【００８０】
比演算部３２４Ｂは、高エネルギー画像ＩＨの手動ＲＯＩ内の画素値の平均値と高エネルギー画像ＩＨのヒストグラムに基づく直接線の画素値との幅と、低エネルギー画像ＩＬの手動ＲＯＩ内の画素値の平均値と低エネルギー画像ＩＬのヒストグラムに基づく直接線の画素値との幅との比を求めることで、手動ＲＯＩ内の同一であるとみなされる部分がキャンセルされるような重み係数ｗを求めることができる。
【００８１】
なお、図１１に示すＸ線診断装置１０Ｂにおいて、図４に示すＸ線診断装置１０と同一機能には同一符号を付して説明を省略する。
【００８２】
第３実施形態のＸ線診断装置１０Ｂによると、予め決定された重み係数ｗのテーブルを用いずに個々の被検体のＸ線画像を基に適宜重み係数ｗを求めるので、目的部位の診断に適したＤＥ画像ＩＤＥを表示できる。
【符号の説明】
【００８３】
１０，１０Ａ，１０Ｂ  Ｘ線診断装置
１３  画像メモリ
１４  画像処理装置
２１  ＣＰＵ
３１  Ｘ線撮影制御部
３２，３２Ａ，３２Ｂ  係数設定部
３３  係数処理部
３４  ＲＯＩ設定部
３２１  ヒストグラム生成部
３２２，３２２Ａ，３２２Ｂ  代表値演算部
３２３，３２３Ａ，３２３Ｂ  幅演算部
３２４，３２４Ａ，３２４Ｂ  比演算部
【図２】
2010240054.tif 000013
【図３】
2010240054.tif 000014
【図４】
2010240054.tif 000015
【図５】
2010240054.tif 000016
【図６】
2010240054.tif 000017
【図７】
2010240054.tif 000018
【図８】
2010240054.tif 000019
【図９】
2010240054.tif 000020
【図１０】
2010240054.tif 000021
【図１１】
2010240054.tif 000022
【図１】
2010240054.tif 000023