JP2021176357A - 医用情報処理装置、医用情報処理方法及び医用情報処理プログラム - Google Patents

Abstract

【課題】物質分離能を向上させること。【解決手段】実施形態の医用情報処理装置は、取得部と、生成部と、出力部とを備える。取得部は、異なるＸ線エネルギーに対応する複数のデータセットに基づいて生成された、複数の指標に応じた複数の第１軸を有する座標系における物質の分布を示す第１分布情報を取得する。生成部は、前記第１分布情報において設定した複数の領域に基づいて前記複数の第１軸とは異なる複数の第２軸を前記第１分布情報において設定し、前記複数の第２軸に応じて前記第１分布情報の変換を行なうことで、前記複数の第２軸を有する座標系における物質の分布を示す第２分布情報を生成する。出力部は、前記第２分布情報に基づく強調画像を出力する。【選択図】図４Ａ

Description

本明細書等に開示の実施形態は、医用情報処理装置、医用情報処理方法及び医用情報処理プログラムに関する。

従来、異なるエネルギーのＸ線を用いた撮影を行ない、物質弁別を行なう技術が知られている。例えば、２種類の異なるエネルギーを用いるデュアルエナジー（Ｄｕａｌ−Ｅｎｅｒｇｙ：ＤＥ）収集では、各Ｘ線エネルギーに対応する２つのデータセットを用いて、予め設定した２つの基準物質それぞれの存在量に基づく基準物質画像を生成することが可能である。

ここで、撮影対象となる被検体の生体組織の組成には個人差があり、年齢や病気の有無等によって変動する場合もある。即ち、生体組織の実際の組成と設定した基準物質との間には乖離が生じる場合があり、この乖離に起因して物質分離能が低下してしまう場合がある。

特開２０１３−８１７６５号公報 国際公開第２０１６／１５８２３４号 特開昭６０−１６３６４１号公報 特開２００７−４４５２０号公報

本明細書等に開示の実施形態が解決しようとする課題の一つは、物質分離能を向上させることである。ただし、本明細書等に開示の実施形態により解決される課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を、本明細書等に開示の実施形態が解決する他の課題として位置付けることもできる。

実施形態の医用情報処理装置は、取得部と、生成部と、出力部とを備える。取得部は、異なるＸ線エネルギーに対応する複数のデータセットに基づいて生成された、複数の指標に応じた複数の第１軸を有する座標系における物質の分布を示す第１分布情報を取得する。生成部は、前記第１分布情報において設定した複数の領域に基づいて前記複数の第１軸とは異なる複数の第２軸を前記第１分布情報において設定し、前記複数の第２軸に応じて前記第１分布情報の変換を行なうことで、前記複数の第２軸を有する座標系における物質の分布を示す第２分布情報を生成する。出力部は、前記第２分布情報に基づく強調画像を出力する。

図１は、第１の実施形態に係る医用情報処理システムの構成の一例を示すブロック図である。 図２は、第１の実施形態に係るＸ線ＣＴ装置の構成の一例を示すブロック図である。 図３は、第１の実施形態に係る基準物質画像及び第１分布情報を示す図である。 図４Ａは、第１の実施形態に係る第２分布情報の生成処理について説明するための図である。 図４Ｂは、第１の実施形態に係る第２軸の設定方法の一例を示す図である。 図４Ｃは、第１の実施形態に係る変換処理の一例を示す図である。 図５は、第１の実施形態に係る第２軸の設定方法の一例を示す図である。 図６は、第１の実施形態に係る第２軸の設定方法の一例を示す図である。 図７は、第１の実施形態に係る第２軸の設定方法の一例を示す図である。 図８は、第１の実施形態に係る医用情報処理装置の処理の一連の流れを説明するためのフローチャートである。

以下、添付図面を参照して、医用情報処理装置、医用情報処理方法及び医用情報処理プログラムの実施形態について詳細に説明する。

（第１の実施形態）
本実施形態では、図１に示すように、医用情報処理装置２０を含んだ医用情報処理システム１を例として説明する。図１は、第１の実施形態に係る医用情報処理システム１の構成の一例を示すブロック図である。また、本実施形態では、図１に示すＸ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）装置１０により収集されたデータセットを用いて物質弁別処理を行なう場合を例として説明する。Ｘ線ＣＴ装置１０及び医用情報処理装置２０は、ネットワークＮＷを介して相互に接続される。

医用情報処理装置２０は、Ｘ線ＣＴ装置１０によって収集されたデータセットに基づく各種の処理を実行する。例えば、医用情報処理装置２０は、図１に示すように、メモリ２１、ディスプレイ２２、入力インタフェース２３及び処理回路２４を備える。

メモリ２１は、例えば、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。例えば、メモリ２１は、医用情報処理装置２０に含まれる回路がその機能を実現するためのプログラムを記憶する。また、メモリ２１は、Ｘ線ＣＴ装置１０から取得したデータセットを記憶する。なお、メモリ２１は、医用情報処理装置２０とネットワークＮＷを介して接続されたサーバ群（クラウド）により実現されることとしてもよい。

ディスプレイ２２は、各種の情報を表示する。例えば、ディスプレイ２２は、処理回路２４によって生成される強調画像を表示する。また、例えば、ディスプレイ２２は、入力インタフェース２３を介してユーザから各種の指示や設定等を受け付けるためのＧＵＩ（Ｇｒａｐｈｉｃａｌ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）を表示する。例えば、ディスプレイ２２は、液晶ディスプレイやＣＲＴ（Ｃａｔｈｏｄｅ Ｒａｙ Ｔｕｂｅ）ディスプレイである。ディスプレイ２２は、デスクトップ型でもよいし、医用情報処理装置２０本体と無線通信可能なタブレット端末等で構成されることにしても構わない。

入力インタフェース２３は、ユーザからの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路２４に出力する。例えば、入力インタフェース２３は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、音声入力回路等により実現される。なお、入力インタフェース２３は、医用情報処理装置２０本体と無線通信可能なタブレット端末等で構成されることにしても構わない。また、入力インタフェース２３は、マウスやキーボード等の物理的な操作部品を備えるものだけに限られない。例えば、医用情報処理装置２０とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を処理回路２４へ出力する電気信号の処理回路も入力インタフェース２３の例に含まれる。

処理回路２４は、制御機能２４ａ、指定機能２４ｂ、取得機能２４ｃ、生成機能２４ｄ及び出力機能２４ｅを実行することで、医用情報処理装置２０全体の動作を制御する。ここで、指定機能２４ｂは、指定部の一例である。また、取得機能２４ｃは、取得部の一例である。また、生成機能２４ｄは、生成部の一例である。また、出力機能２４ｅは、出力部の一例である。

例えば、処理回路２４は、制御機能２４ａに対応するプログラムをメモリ２１から読み出して実行することにより、入力インタフェース２３を介してユーザから受け付けた各種の入力操作に基づいて、指定機能２４ｂ、取得機能２４ｃ、生成機能２４ｄ及び出力機能２４ｅといった各種の機能を制御する。

また、例えば、処理回路２４は、指定機能２４ｂに対応するプログラムをメモリ２１から読み出して実行することにより、被検体における関心領域を指定する。また、例えば、処理回路２４は、取得機能２４ｃに対応するプログラムをメモリ２１から読み出して実行することにより、第１分布情報を取得する。また、例えば、処理回路２４は、生成機能２４ｄに対応するプログラムをメモリ２１から読み出して実行することにより、第１分布情報に基づいて第２分布情報を生成する。また、例えば、処理回路２４は、出力機能２４ｅに対応するプログラムをメモリ２１から読み出して実行することにより、第２分布情報に基づく強調画像を出力する。なお、指定機能２４ｂ、取得機能２４ｃ、生成機能２４ｄ及び出力機能２４ｅによる処理については後述する。

図１に示す医用情報処理装置２０においては、各処理機能がコンピュータによって実行可能なプログラムの形態でメモリ２１へ記憶されている。処理回路２４は、メモリ２１からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、プログラムを読み出した状態の処理回路２４は、読み出したプログラムに対応する機能を有することとなる。

なお、図１においては単一の処理回路２４にて、制御機能２４ａ、指定機能２４ｂ、取得機能２４ｃ、生成機能２４ｄ及び出力機能２４ｅが実現するものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路２４を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。また、処理回路２４が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。

また、処理回路２４は、ネットワークＮＷを介して接続された外部装置のプロセッサを利用して、機能を実現することとしてもよい。例えば、処理回路２４は、メモリ２１から各機能に対応するプログラムを読み出して実行するとともに、医用情報処理装置２０とネットワークＮＷを介して接続されたサーバ群（クラウド）を計算資源として利用することにより、図１に示す各機能を実現する。

次に、Ｘ線ＣＴ装置１０について、図２を用いて説明する。図２は、第１の実施形態に係るＸ線ＣＴ装置１０の構成の一例を示すブロック図である。例えば、Ｘ線ＣＴ装置１０は、架台装置１１０と、寝台装置１３０と、コンソール装置１４０とを有する。

図２においては、非チルト状態での回転フレーム１１３の回転軸又は寝台装置１３０の天板１３３の長手方向をＺ軸方向とする。また、Ｚ軸方向に直交し、床面に対し水平である軸方向をＸ軸方向とする。また、Ｚ軸方向に直交し、床面に対し垂直である軸方向をＹ軸方向とする。なお、図２は、説明のために架台装置１１０を複数方向から描画したものであり、Ｘ線ＣＴ装置１０が架台装置１１０を１つ有する場合を示す。

架台装置１１０は、Ｘ線管１１１と、Ｘ線検出器１１２と、回転フレーム１１３と、Ｘ線高電圧装置１１４と、制御装置１１５と、ウェッジ１１６と、コリメータ１１７と、ＤＡＳ１１８とを有する。

Ｘ線管１１１は、熱電子を発生する陰極（フィラメント）と、熱電子の衝突を受けてＸ線を発生する陽極（ターゲット）とを有する真空管である。Ｘ線管１１１は、Ｘ線高電圧装置１１４からの高電圧の印加により、陰極から陽極に向けて熱電子を照射することで、被検体Ｐに対し照射するＸ線を発生する。

Ｘ線検出器１１２は、Ｘ線管１１１から照射されて被検体Ｐを通過したＸ線を検出し、検出したＸ線量に対応した信号をＤＡＳ１１８へと出力する。Ｘ線検出器１１２は、例えば、Ｘ線管１１１の焦点を中心とした１つの円弧に沿ってチャンネル方向（チャネル方向）に複数の検出素子が配列された複数の検出素子列を有する。Ｘ線検出器１１２は、例えば、チャネル方向に複数の検出素子が配列された検出素子列が列方向（スライス方向、roｗ方向）に複数配列された構造を有する。

例えば、Ｘ線検出器１１２は、グリッドと、シンチレータアレイと、光センサアレイとを有する間接変換型の検出器である。シンチレータアレイは、複数のシンチレータを有する。シンチレータは入射Ｘ線量に応じた光子量の光を出力するシンチレータ結晶を有する。グリッドは、シンチレータアレイのＸ線入射側の面に配置され、散乱Ｘ線を吸収するＸ線遮蔽板を有する。なお、グリッドはコリメータ（１次元コリメータ又は２次元コリメータ）と呼ばれる場合もある。光センサアレイは、シンチレータからの光量に応じた電気信号に変換する機能を有し、例えば、フォトダイオード等の光センサを有する。なお、Ｘ線検出器１１２は、入射したＸ線を電気信号に変換する半導体素子を有する直接変換型の検出器であっても構わない。

回転フレーム１１３は、Ｘ線管１１１とＸ線検出器１１２とを対向支持し、制御装置１１５によってＸ線管１１１とＸ線検出器１１２とを回転させる円環状のフレームである。例えば、回転フレーム１１３は、アルミニウムを材料とした鋳物である。なお、回転フレーム１１３は、Ｘ線管１１１及びＸ線検出器１１２に加えて、Ｘ線高電圧装置１１４やウェッジ１１６、コリメータ１１７、ＤＡＳ１１８等を更に支持することもできる。更に、回転フレーム１１３は、図２において図示しない種々の構成を更に支持することもできる。以下では、架台装置１１０において、回転フレーム１１３、及び、回転フレーム１１３と共に回転移動する部分を、回転部とも記載する。

Ｘ線高電圧装置１１４は、変圧器（トランス）及び整流器等の電気回路を有し、Ｘ線管１１１に印加する高電圧を発生する高電圧発生装置と、Ｘ線管１１１が発生するＸ線に応じた出力電圧の制御を行なうＸ線制御装置とを有する。高電圧発生装置は、変圧器方式であってもよいし、インバータ方式であってもよい。なお、Ｘ線高電圧装置１１４は、回転フレーム１１３に設けられてもよいし、図示しない固定フレームに設けられても構わない。

制御装置１１５は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等を有する処理回路と、モータ及びアクチュエータ等の駆動機構とを有する。制御装置１１５は、入力インタフェース１４３からの入力信号を受けて、架台装置１１０及び寝台装置１３０の動作制御を行なう。例えば、制御装置１１５は、回転フレーム１１３の回転や架台装置１１０のチルト、寝台装置１３０の動作等について制御を行なう。一例を挙げると、制御装置１１５は、架台装置１１０をチルトさせる制御として、入力された傾斜角度（チルト角度）情報により、Ｘ軸方向に平行な軸を中心に回転フレーム１１３を回転させる。なお、制御装置１１５は架台装置１１０に設けられてもよいし、コンソール装置１４０に設けられてもよい。

ウェッジ１１６は、Ｘ線管１１１から照射されたＸ線量を調節するためのＸ線フィルタである。具体的には、ウェッジ１１６は、Ｘ線管１１１から被検体Ｐへ照射されるＸ線が予め定められた分布になるように、Ｘ線管１１１から照射されたＸ線を減衰させるＸ線フィルタである。例えば、ウェッジ１１６は、ウェッジフィルタ（ｗｅｄｇｅ ｆｉｌｔｅｒ）やボウタイフィルタ（ｂｏｗ−ｔｉｅ ｆｉｌｔｅｒ）であり、所定のターゲット角度や所定の厚みとなるようにアルミニウム等を加工して作製される。

コリメータ１１７は、ウェッジ１１６を透過したＸ線の照射範囲を絞り込むための鉛板等であり、複数の鉛板等の組み合わせによってスリットを形成する。なお、コリメータ１１７は、Ｘ線絞りと呼ばれる場合もある。また、図２においては、Ｘ線管１１１とコリメータ１１７との間にウェッジ１１６が配置される場合を示すが、Ｘ線管１１１とウェッジ１１６との間にコリメータ１１７が配置される場合であってもよい。この場合、ウェッジ１１６は、Ｘ線管１１１から照射され、コリメータ１１７により照射範囲が制限されたＸ線を透過して減衰させる。

ＤＡＳ１１８は、Ｘ線検出器１１２が有する各検出素子によって検出されるＸ線の信号を収集する。例えば、ＤＡＳ１１８は、各検出素子から出力される電気信号に対して増幅処理を行なう増幅器と、電気信号をデジタル信号に変換するＡ／Ｄ変換器とを有し、検出データを生成する。ＤＡＳ１１８は、例えば、プロセッサにより実現される。

ＤＡＳ１１８が生成したデータは、回転フレーム１１３に設けられた発光ダイオード（Ｌｉｇｈｔ Ｅｍｉｔｔｉｎｇ Ｄｉｏｄｅ：ＬＥＤ）を有する送信機から、光通信によって、架台装置１１０の非回転部分（例えば、固定フレーム等。図２での図示は省略している）に設けられた、フォトダイオードを有する受信機に送信され、コンソール装置１４０へと転送される。ここで、非回転部分とは、例えば、回転フレーム１１３を回転可能に支持する固定フレーム等である。なお、回転フレーム１１３から架台装置１１０の非回転部分へのデータの送信方法は、光通信に限らず、非接触型の如何なるデータ伝送方式を採用してもよいし、接触型のデータ伝送方式を採用しても構わない。

寝台装置１３０は、スキャン対象の被検体Ｐを載置、移動させる装置であり、基台１３１と、寝台駆動装置１３２と、天板１３３と、支持フレーム１３４とを有する。基台１３１は、支持フレーム１３４を鉛直方向に移動可能に支持する筐体である。寝台駆動装置１３２は、被検体Ｐが載置された天板１３３を、天板１３３の長軸方向に移動する駆動機構であり、モータ及びアクチュエータ等を含む。支持フレーム１３４の上面に設けられた天板１３３は、被検体Ｐが載置される板である。なお、寝台駆動装置１３２は、天板１３３に加え、支持フレーム１３４を天板１３３の長軸方向に移動してもよい。

コンソール装置１４０は、メモリ１４１と、ディスプレイ１４２と、入力インタフェース１４３と、処理回路１４４とを有する。なお、コンソール装置１４０は架台装置１１０とは別体として説明するが、架台装置１１０にコンソール装置１４０又はコンソール装置１４０の各構成要素の一部が含まれてもよい。

メモリ１４１は、例えば、ＲＡＭ、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。例えば、メモリ１４１は、被検体Ｐに対するスキャンを実行することで収集される各種のデータを記憶する。また、例えば、メモリ１４１は、Ｘ線ＣＴ装置１０に含まれる回路がその機能を実現するためのプログラムを記憶する。なお、メモリ１４１は、Ｘ線ＣＴ装置１０とネットワークＮＷを介して接続されたサーバ群（クラウド）により実現されることとしてもよい。

ディスプレイ１４２は、各種の情報を表示する。例えば、ディスプレイ１４２は、ユーザからの各種操作を受け付けるためのＧＵＩや、被検体Ｐから収集された各種の画像を表示する。例えば、ディスプレイ１４２は、液晶ディスプレイやＣＲＴディスプレイである。ディスプレイ１４２は、デスクトップ型でもよいし、コンソール装置１４０本体と無線通信可能なタブレット端末等で構成されることにしても構わない。

入力インタフェース１４３は、ユーザからの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路１４４に出力する。例えば、入力インタフェース１４３は、スキャン条件の設定やスキャン開始の指示等を入力する操作をユーザから受け付ける。例えば、入力インタフェース１４３は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、操作面へ触れることで入力操作を行なうタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、音声入力回路等により実現される。なお、入力インタフェース１４３は、架台装置１１０に設けられてもよい。また、入力インタフェース１４３は、コンソール装置１４０本体と無線通信可能なタブレット端末等で構成されることにしても構わない。また、入力インタフェース１４３は、マウスやキーボード等の物理的な操作部品を備えるものだけに限られない。例えば、コンソール装置１４０とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を処理回路１４４へ出力する電気信号の処理回路も入力インタフェース１４３の例に含まれる。

処理回路１４４は、制御機能１４４ａ、収集機能１４４ｂ及び出力機能１４４ｃを実行することで、Ｘ線ＣＴ装置１０全体の動作を制御する。

例えば、処理回路１４４は、制御機能１４４ａに対応するプログラムをメモリ１４１から読み出して実行することにより、入力インタフェース１４３を介してユーザから受け付けた各種の入力操作に基づいて、収集機能１４４ｂ、出力機能１４４ｃといった各種の機能を制御する。

また、例えば、処理回路１４４は、収集機能１４４ｂに対応するプログラムをメモリ１４１から読み出して実行することにより、被検体Ｐに対するスキャンを実行する。例えば、収集機能１４４ｂは、Ｘ線高電圧装置１１４を制御することにより、Ｘ線管１１１に高電圧を供給する。これにより、Ｘ線管１１１は、被検体Ｐに対し照射するＸ線を発生する。また、収集機能１４４ｂは、寝台駆動装置１３２を制御することにより、被検体Ｐを架台装置１１０の撮影口内へ移動させる。また、収集機能１４４ｂは、ウェッジ１１６の位置、及び、コリメータ１１７の開口度及び位置を調整することで、被検体Ｐに照射されるＸ線の分布を制御する。また、収集機能１４４ｂは、制御装置１１５を制御することにより回転部を回転させる。また、収集機能１４４ｂによってスキャンが実行される間、ＤＡＳ１１８は、Ｘ線検出器１１２における各検出素子からＸ線の信号を収集し、検出データを生成する。また、収集機能１４４ｂは、ＤＡＳ１１８から出力された検出データに対して、前処理を施す。例えば、収集機能１４４ｂは、ＤＡＳ１１８から出力された検出データに対して、対数変換処理やオフセット補正処理、チャンネル間の感度補正処理、ビームハードニング補正等の前処理を施す。なお、前処理を施した後のデータについては生データとも記載する。また、前処理を施す前の検出データ及び前処理を施した後の生データを総称して、投影データとも記載する。

また、例えば、処理回路１４４は、出力機能１４４ｃに対応するプログラムをメモリ１４１から読み出して実行することにより、被検体Ｐから収集された各種のデータを出力する。例えば、出力機能１４４ｃは、被検体Ｐに対するスキャンを実行することで収集されたデータを、ネットワークＮＷを介して、医用情報処理装置２０に送信する。また、例えば、出力機能１４４ｃは、ディスプレイ１４２における表示の制御を行なう。

図２に示すＸ線ＣＴ装置１０においては、各処理機能がコンピュータによって実行可能なプログラムの形態でメモリ１４１へ記憶されている。処理回路１４４は、メモリ１４１からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、プログラムを読み出した状態の処理回路１４４は、読み出したプログラムに対応する機能を有することとなる。

なお、図２においては単一の処理回路１４４にて、制御機能１４４ａ、収集機能１４４ｂ及び出力機能１４４ｃが実現するものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路１４４を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。また、処理回路１４４が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。

また、処理回路１４４は、ネットワークＮＷを介して接続された外部装置のプロセッサを利用して、機能を実現することとしてもよい。例えば、処理回路１４４は、メモリ１４１から各機能に対応するプログラムを読み出して実行するとともに、Ｘ線ＣＴ装置１０とネットワークＮＷを介して接続されたサーバ群（クラウド）を計算資源として利用することにより、図２に示す各機能を実現する。

以上、Ｘ線ＣＴ装置１０及び医用情報処理装置２０を含んだ医用情報処理システム１の構成例について説明した。以下、医用情報処理システム１において行なわれる処理について詳細に説明する。

物質弁別処理を行なう場合、Ｘ線ＣＴ装置１０における収集機能１４４ｂは、被検体Ｐに対するスキャンを行なって、異なるＸ線エネルギーに対応する複数のデータセットを収集する。例えば、収集機能１４４ｂは、被検体Ｐに対するデュアルエナジーのスキャンを行なって、第１のＸ線エネルギーに対応する第１のデータセットと、第１のＸ線エネルギーよりも高い第２のＸ線エネルギーに対応する第２のデータセットとを収集する。

本実施形態では、デュアルエナジー収集が行われ、２つの異なるＸ線エネルギーに対応する第１のデータセット及び第２のデータセットが収集された場合を例として説明する。なお、デュアルエナジー収集の方式については特に限定されるものではない。例えば、収集機能１４４ｂは、Ｓｌｏｗ−ｋＶスイッチングや、Ｆａｓｔ−ｋＶスイッチング、積層型検出器（デュアルレイヤー）、デュアルソース方式、スプリット方式といった方式により、デュアルエナジー収集を実行することができる。

Ｓｌｏｗ−ｋＶスイッチングの方式でデュアルエナジー収集を実行する場合、収集機能１４４ｂは、第１のＸ線エネルギーのＸ線を用いて１回目のスキャンを行なった後、第２のＸ線エネルギーのＸ線を用いて２回目のスキャンを行なう。これにより、収集機能１４４ｂは、１回目のスキャンにおいて第１のＸ線エネルギーに対応する第１のデータセットを収集し、２回目のスキャンにおいて第２のＸ線エネルギーに対応する第２のデータセットを収集することができる。なお、収集機能１４４ｂは、Ｘ線高電圧装置１１４を制御し、Ｘ線管１１１に供給する管電圧を制御することで、第１のＸ線エネルギーと第２のＸ線エネルギーとの切り替えを行なうことができる。

Ｆａｓｔ−ｋＶスイッチングの方式でデュアルエナジー収集を実行する場合、収集機能１４４ｂは、１回のスキャンの中でビューごとにＸ線のエネルギーを変化させて、第１のデータセット及び第２のデータセットを収集する。即ち、スキャン中、Ｘ線管１１１は図２に示したＺ軸を回転軸として被検体Ｐの周囲を回転し、照射角度を変えながら被検体ＰにＸ線を照射する。ここで、収集機能１４４ｂは、照射角度（ビュー）ごとに、Ｘ線のエネルギーを第１のＸ線エネルギーと第２のＸ線エネルギーとの間で高速に切り替えることで、第１のデータセット及び第２のデータセットを収集する。

デュアルレイヤーの方式でデュアルエナジー収集を実行する場合、Ｘ線ＣＴ装置１０は、Ｘ線検出器１１２として、積層型検出器を備える。例えば、Ｘ線検出器１１２は、第１の層１１２ａと、第２の層１１２ｂとから構成され、Ｘ線管１１１から照射されたＸ線を分光して検出する。例えば、第１の層１１２ａは第１のＸ線エネルギーのＸ線を検出し、第２の層１１２ｂは、第１の層１１２ａを透過した第２のＸ線エネルギーのＸ線を検出する。この場合、収集機能１４４ｂは、第１の層１１２ａからの出力に基づいて第１のデータセットを収集し、第２の層１１２ｂからの出力に基づいて第２のデータセットを収集することができる。

デュアルソースの方式でデュアルエナジー収集を実行する場合、Ｘ線ＣＴ装置１０は、Ｘ線管１１１として、第１のＸ線管１１１ａ及び第２のＸ線管１１１ｂを備える。また、Ｘ線ＣＴ装置１０は、Ｘ線検出器１１２として、第１のＸ線管１１１ａから照射されたＸ線を検出する第１のＸ線検出器１１２ｃと、第２のＸ線管１１１ｂから照射されたＸ線を検出する第２のＸ線検出器１１２ｄとを備える。この場合、収集機能１４４ｂは、第１のＸ線管１１１ａから第１のＸ線エネルギーのＸ線を照射させ、第２のＸ線管１１１ｂから第２のＸ線エネルギーのＸ線を照射させることで、第１のＸ線検出器１１２ｃからの出力に基づいて第１のデータセットを収集し、第２のＸ線検出器１１２ｄからの出力に基づいて第２のデータセットを収集することができる。

スプリット方式でデュアルエナジー収集を実行する場合、Ｘ線ＣＴ装置１０は、ウェッジ１１６として、Ｘ線管１１１から照射されたＸ線をエネルギーの異なる複数のＸ線に分割するフィルタを備える。例えば、Ｘ線ＣＴ装置１０は、ウェッジ１１６として、Ｘ線フィルタ１１６ａ及びＸ線フィルタ１１６ｂを備える。この場合、収集機能１４４ｂは、Ｘ線管１１１から照射されてＸ線フィルタ１１６ａを透過したＸ線の検出結果に基づいて第１のデータセットを収集し、Ｘ線管１１１から照射されてＸ線フィルタ１１６ｂを透過したＸ線の検出結果に基づいて第２のデータセットを収集することができる。

次に、出力機能１４４ｃは、第１のデータセット及び第２のデータセットを医用情報処理装置２０に送信する。また、取得機能２４ｃは、第１のデータセット及び第２のデータセットを受信して、メモリ２１に記憶させる。なお、出力機能１４４ｃが送信する第１のデータセット及び第２のデータセットは、前処理前のデータでもよいし前処理後のデータでもよい。前処理前の状態で第１のデータセット及び第２のデータセットを送信する場合には、取得機能２４ｃが第１のデータセット及び第２のデータセットに対する前処理を実行することができる。

次に、取得機能２４ｃは、第１のデータセット及び第２のデータセットに基づいて、物質の分布を示す第１分布情報を生成する。例えば、取得機能２４ｃは、基準物質として第１物質及び第２物質を定義し、第１物質の存在量に関する指標に応じた軸、及び、第２物質の存在量に関する指標に応じた軸を有する座標系における物質の分布を示す散布図を生成する。以下、第１物質として水（Ｈ２Ｏ）が定義され、第２物質としてカルシウム（Ｃａ）が定義されている場合について説明する。水の存在量に関する指標は、第１指標の一例である。また、カルシウムの存在量に関する指標は、第２指標の一例である。

例えば、取得機能２４ｃは、まず、第１のデータセット及び第２のデータセットに基づく物質弁別処理を行ない、撮影対象部位に存在する水及びカルシウムを分離する。例えば、取得機能２４ｃは、第１のデータセット及び第２のデータセットにおける画素ごとに、水及びカルシウムそれぞれの存在量に関する指標を推定する。一例を挙げると、取得機能２４ｃは、画素ごとに、水に対するカルシウムの存在比率と、カルシウムに対する水の存在比率とをそれぞれ推定する。また、取得機能２４ｃは、水成分に対応する単色Ｘ線の投影データセット及びカルシウム成分に対応する単色Ｘ線の投影データセットをそれぞれ生成する。

次に、取得機能２４ｃは、水成分に対応する単色Ｘ線の投影データセットに基づいて水成分に対応するＣＴ画像（ボリュームデータ）を再構成し、カルシウム成分に対応する単色Ｘ線の投影データセットに基づいてカルシウム成分に対応するＣＴ画像を再構成する。次に、取得機能２４ｃは、再構成したＣＴ画像に対する画像処理を行なって、ＭＰＲ画像等の表示用画像を生成する。例えば、取得機能２４ｃは、水成分に対応するＣＴ画像に基づいて基準物質画像Ｉｍ１１を生成し、カルシウム成分に対応するＣＴ画像に基づいて基準物質画像Ｉｍ１２を生成する。

なお、投影データの段階で物質弁別処理を行なうものとして説明したが、実施形態はこれに限定されるものではない。例えば、取得機能２４ｃは、まず、第１のデータセットに基づくＣＴ画像と、第２のデータセットに基づくＣＴ画像とをそれぞれ再構成する。次に、取得機能２４ｃは、第１のデータセットに基づくＣＴ画像及び第２のデータセットに基づくＣＴ画像に基づいて水及びカルシウムを分離し、水成分に対応するＣＴ画像及びカルシウム成分に対応するＣＴ画像をそれぞれ生成する。そして、取得機能２４ｃは、水成分に対応するＣＴ画像に基づいて基準物質画像Ｉｍ１１を生成し、カルシウム成分に対応するＣＴ画像に基づいて基準物質画像Ｉｍ１２を生成する。

次に、出力機能２４ｅは、基準物質画像をディスプレイ２２に表示させる。例えば、出力機能２４ｅは、図３に示すように、水成分に対応する基準物質画像Ｉｍ１１をディスプレイ２２に表示させる。なお、図３は、第１の実施形態に係る基準物質画像Ｉｍ１１及び第１分布情報を示す図である。

ここで、指定機能２４ｂは、被検体Ｐにおける関心領域を指定する。例えば、図３においては被検体Ｐにおける３つの椎骨が基準物質画像Ｉｍ１１に含まれており、指定機能２４ｂは、これら３つの椎骨のうち中心の椎骨に対応する領域を、関心領域として指定する。

なお、指定機能２４ｂは、基準物質画像Ｉｍ１１を参照したユーザからの入力操作を受け付けることにより関心領域を指定してもよいし、基準物質画像Ｉｍ１１から骨領域を自動抽出することで関心領域を指定してもよい。

或いは、指定機能２４ｂは、基準物質画像Ｉｍ１１以外の画像に基づいて関心領域を指定してもよい。例えば、取得機能２４ｃは、第１のデータセットに基づいてＣＴ画像を再構成し、再構成したＣＴ画像に対する画像処理を行なって、図３の基準物質画像Ｉｍ１１と同じ位置のＭＰＲ画像を生成する。なお、第１のデータセットに基づくＭＰＲ画像については、低エネルギー画像とも記載する。同様に、取得機能２４ｃは、第２のデータセットに基づいてＣＴ画像を再構成し、再構成したＣＴ画像に対する画像処理を行なって、高エネルギー画像を生成する。そして、指定機能２４ｂは、低エネルギー画像や高エネルギー画像等に基づいて、関心領域を指定する。

次に、取得機能２４ｃは、指定された関心領域内の画素について、物質の分布を示す第１分布情報を生成する。例えば、取得機能２４ｃは、水の存在量に関する指標に応じて横軸を定義し、カルシウムの存在量に関する指標に応じて横軸を定義した座標系において、関心領域内の各画素をプロットすることにより、図３の右図に示す散布図Ｓｃ１１を生成する。なお、図３の散布図Ｓｃ１１における横軸及び縦軸は、第１軸の一例である。

散布図Ｓｃ１１においては、各画素のプロットを白点で示す。即ち、より水らしい画素は散布図Ｓｃ１１においてより右側にプロットされ、よりカルシウムらしい画素は散布図Ｓｃ１１においてより上側にプロットされる。なお、散布図Ｓｃ１１において横軸と縦軸のスケールは異なっている。具体的には、散布図Ｓｃ１１における縦軸は、横軸と比較して圧縮されている。

図３に示すように、関心領域（３つの椎骨のうち中心の椎骨に対応する領域）には、浮腫が生じている。浮腫に相当する領域には他の正常領域と比較して水成分がより多く存在するため、水成分に対応する基準物質画像Ｉｍ１１を用いれば、浮腫の観察を行なうことが可能である。

具体的には、基準物質画像Ｉｍ１１の正常領域における各画素は、散布図Ｓｃ１１における領域Ｒ１１にプロットされる。また、基準物質画像Ｉｍ１１の浮腫領域における各画素は、散布図Ｓｃ１１における領域Ｒ２１にプロットされる。なお、図３においては、領域Ｒ１１の中心及び領域Ｒ２１の中心を、それぞれ「×」のマークで示す。図３に示すように、領域Ｒ１１と領域Ｒ２１とは横軸方向にずれているため、領域Ｒ１１に対応する正常領域と領域Ｒ２１に対応する浮腫領域とは、基準物質画像Ｉｍ１１において色が異なることとなる。より具体的には、基準物質画像Ｉｍ１１における正常領域及び浮腫領域は、領域Ｒ１１及び領域Ｒ２１の中心間の距離Ｄに応じて平均画素値に差が生じるため、基準物質画像Ｉｍ１１を参照したユーザは、浮腫領域を正常領域と区別して観察することができる。

ここで、理論値としては、散布図Ｓｃ１１における領域Ｒ１１及び領域Ｒ２１は円形となるものである。例えば、図３に示すように散布図Ｓｃ１１の縦軸を圧縮している場合、理論値としては、領域Ｒ１１及び領域Ｒ２１は、横軸と平行な長軸を有する楕円形となるものである。

しかしながら、被検体Ｐの生体組織の組成によって、散布図Ｓｃ１１は理論値通りとならない場合が多い。例えば、浮腫領域の主要構成物質が純粋な水ではなく血液であること、骨密度には個人差があること等から、設定した基準物質と実際の生体組織の組成との間には乖離が生じる。この結果、領域Ｒ１１及び領域Ｒ２１は、散布図Ｓｃ１１の横軸と平行な長軸を有する楕円形とはならず、例えば長軸が散布図Ｓｃ１１の横軸に対して傾いた楕円形になることがある。

このように、領域Ｒ１１及び領域Ｒ２１が散布図Ｓｃ１１において斜め方向に分布することによって、領域Ｒ１１及び領域Ｒ２１の中心間の距離Ｄが小さくなり、正常領域と浮腫領域との間のコントラストが小さくなる場合がある。即ち、領域Ｒ１１及び領域Ｒ２１が散布図Ｓｃ１１において斜め方向に分布することによって、物質分離能が低下し、浮腫領域を観察しにくくなる場合がある。

そこで、医用情報処理装置２０は、散布図Ｓｃ１１等の第１分布情報に基づいて第２分布情報を更に生成し、第２分布情報に基づく強調画像を出力することで、物質分離能を向上させる。

例えば、生成機能２４ｄは、まず、散布図Ｓｃ１１の横軸及び縦軸とは異なる複数の軸を、散布図Ｓｃ１１において設定する。即ち、生成機能２４ｄは、複数の第１軸とは異なる複数の第２軸を、第１分布情報において設定する。ここで、生成機能２４ｄは、領域Ｒ１１と領域Ｒ２１との間の分離能が高くなるように、第２軸の設定を行なう。

一例を挙げると、生成機能２４ｄは、図４Ａに示すように、散布図Ｓｃ１１において、散布図Ｓｃ１１の横軸及び縦軸とは異なる軸Ａ１１及びＡ２１を設定する。なお、図４Ａは、第１の実施形態に係る第２分布情報の生成処理について説明するための図である。

具体的には、生成機能２４ｄは、散布図Ｓｃ１１において、領域Ｒ１１及び領域Ｒ２１を設定する。例えば、生成機能２４ｄは、散布図Ｓｃ１１において、基準物質画像Ｉｍ１１の正常領域の画素に対応するプロットのうち所定の割合以上のプロットが含まれる楕円形であって、面積が最小となる楕円形となるように、領域Ｒ１１を設定する。同様に、生成機能２４ｄは、散布図Ｓｃ１１において、基準物質画像Ｉｍ１１の浮腫領域の画素に対応するプロットのうち所定の割合以上のプロットが含まれる楕円形であって、面積が最小となる楕円形となるように、領域Ｒ２１を設定する。

次に、生成機能２４ｄは、領域Ｒ１１及び領域Ｒ２１のいずれかを選択する。以下では、領域Ｒ１１を選択したものとして説明する。次に、生成機能２４ｄは、選択した領域Ｒ１１の長軸に沿って軸Ａ１１を設定する。例えば、生成機能２４ｄは、楕円形である領域Ｒ１１の円周と、領域Ｒ１１の長軸との交点となる２点を結ぶ直線に沿って、軸Ａ１１を設定する。なお、軸Ａ１１は、第２軸及び第３軸の一例である。

また、生成機能２４ｄは、物質の存在量の変化に応じた軸を、軸Ａ２１として設定する。例えば、生成機能２４ｄは、カルシウムの存在量（骨密度）を変化させた複数の対象物について生成された、物質の分布を示す分布情報を取得する。ここで、生成機能２４ｄは、散布図Ｓｃ１１の横軸及び縦軸と同じ軸を有する座標系における物質の分布を示す分布情報を取得する。また、生成機能２４ｄは、取得した分布情報において、複数の対象物それぞれに対応する複数の対応領域を設定する。また、生成機能２４ｄは、設定した複数の対応領域に基づいて、軸Ａ２１を設定する。なお、軸Ａ２１は、第２軸及び第４軸の一例である。

ここで、軸Ａ２１の設定について、図４Ｂを用いてより詳細に説明する。図４Ｂは、第１の実施形態に係る第２軸の設定方法の一例を示す図である。

例えば、まず、骨密度を変化させた３つの対象物について、Ｘ線ＣＴ装置１０によるスキャンが実行される。例えば、収集機能１４４ｂは、骨密度の低いサンプルＳａ１、サンプルＳａ１より骨密度の高いサンプルＳａ２、サンプルＳａ２より骨密度の高いサンプルＳａ３を埋め込んだファントムについてデュアルエナジーのスキャンを実行し、異なるＸ線エネルギーに対応する２つのデータセットを収集する。

次に、生成機能２４ｄは、サンプルＳａ１、サンプルＳａ２及びサンプルＳａ３についての散布図を生成する。例えば、取得機能２４ｃは、水の存在量に関する指標に応じて横軸を定義し、カルシウムの存在量に関する指標に応じて横軸を定義した座標系において、サンプルＳａ１に対応する各画素、サンプルＳａ２に対応する各画素、及び、サンプルＳａ３に対応する各画素をそれぞれプロットすることにより、散布図を生成する。

次に、生成機能２４ｄは、生成した散布図においてサンプルＳａ１に対応する領域Ｒ３１を設定する。例えば、生成機能２４ｄは、生成した散布図において、サンプルＳａ１に対応する各画素のプロットのうち所定の割合以上のプロットが含まれる楕円形であって、面積が最小となる楕円形となるように、図４Ｂに示す領域Ｒ３１を設定する。同様に、生成機能２４ｄは、サンプルＳａ２に対応する領域Ｒ３２及びサンプルＳａ３に対応する領域Ｒ３３をそれぞれ設定する。領域Ｒ３１、領域Ｒ３２及び領域Ｒ３３は、対応領域の一例である。

次に、生成機能２４ｄは、領域Ｒ３１、領域Ｒ３２及び領域Ｒ３３に基づいて、図４Ｂに示す軸Ａ３を設定する。一例を挙げると、生成機能２４ｄは、まず、領域Ｒ３１の中心を特定する。なお、領域Ｒ３１の中心とは、楕円形の中心（長軸と短軸の交点）であってもよいし、サンプルＳａ１に対応する各画素のプロットの平均値であってもよい。同様に、生成機能２４ｄは、領域Ｒ３２及び領域Ｒ３３それぞれの中心を特定する。

次に、生成機能２４ｄは、領域Ｒ３１、領域Ｒ３２及び領域Ｒ３３の中心に沿って、軸Ａ３を設定する。例えば、生成機能２４ｄは、領域Ｒ３１、領域Ｒ３２及び領域Ｒ３３それぞれの中心座標について、最小二乗法によって一次関数を近似させることにより、軸Ａ３を設定する。更に、生成機能２４ｄは、散布図Ｓｃ１１上で軸Ａ３を平行移動させ、領域Ｒ１１の中心を通るように配置することで、軸Ａ２１を設定する。

なお、出力機能２４ｅは、散布図Ｓｃ１１や、散布図Ｓｃ１１において設定された軸Ａ１１及び軸Ａ２１をディスプレイ２２に表示させることとしてもよい。また、生成機能２４ｄは、設定した軸Ａ１１及び軸Ａ２１について、ユーザによる調整を受け付けることとしてもよい。

軸Ａ１１及び軸Ａ２１を設定した後、生成機能２４ｄは、軸Ａ１１及び軸Ａ２１に応じて散布図Ｓｃ１１の変換を行なうことで、図４Ａに示す散布図Ｓｃ２１を生成する。散布図Ｓｃ２１は、軸Ａ１１及び軸Ａ２１を有する座標系における物質の分布を示す散布図であり、第２分布情報の一例である。

ここで、生成機能２４ｄによる散布図Ｓｃ１１の変換処理について、図４Ｃを用いて説明する。図４Ｃは、第１の実施形態に係る変換処理の一例を示す図である。例えば、生成機能２４ｄは、まず、散布図Ｓｃ２１の原点となる座標を、散布図Ｓｃ１１において決定する。例えば、生成機能２４ｄは、図４Ｃに示すように、散布図Ｓｃ１１において軸Ａ１１と軸Ａ２１との交点となる座標（ｘ_０，ｙ_０）を、散布図Ｓｃ２１の原点となる座標として決定する。換言すると、生成機能２４ｄは、領域Ｒ１１の中心を、散布図Ｓｃ２１の原点となる座標として決定する。

次に、生成機能２４ｄは、軸Ａ１１の成分について投影を行なう。例えば、生成機能２４ｄは、散布図Ｓｃ１１上の任意の点である点Ｐ（ｘ，ｙ）について、軸Ａ２１に沿って延長した場合に軸Ａ１１と交差する点ａ（ｘ_ａ，ｙ_ａ）を求める。また、生成機能２４ｄは、点Ｐ（ｘ，ｙ）と点ａ（ｘ_ａ，ｙ_ａ）との間の距離（ｘ_ａ’）を求め、点Ｐ（ｘ，ｙ）を変換した後の横軸座標値（ｘ_ａ’，０）とする。

同様に、生成機能２４ｄは、軸Ａ２１の成分について投影を行なう。例えば、生成機能２４ｄは、点Ｐ（ｘ，ｙ）について、軸Ａ１１に沿って延長した場合に軸Ａ２１と交差する点ｂ（ｘ_ｂ，ｙ_ｂ）を求める。また、生成機能２４ｄは、点Ｐ（ｘ，ｙ）と点ｂ（ｘ_ｂ，ｙ_ｂ）との間の距離（ｙ_ｂ’）を求め、点Ｐ（ｘ，ｙ）を変換した後の縦軸座標値（０，ｙ_ｂ’）とする。これにより、点Ｐ（ｘ，ｙ）は、点Ｐ（ｘ’，ｙ’）に座標変換される。また、生成機能２４ｄは、点Ｐ（ｘ，ｙ）と同様にして散布図Ｓｃ１１上の各ピクセルを変換することにより、散布図Ｓｃ２１を生成することができる。

なお、図４Ｃの変換処理はあくまで一例であり、種々の変形が可能である。例えば、生成機能２４ｄは、軸Ａ１１と軸Ａ２１との交点となる座標（ｘ_０，ｙ_０）以外の任意の座標を原点として、変換処理を行なうことが可能である。

更に、出力機能２４ｅは、散布図Ｓｃ２１に基づく基準物質画像Ｉｍ２１を出力する。例えば、出力機能２４ｅは、散布図Ｓｃ２１上にプロットされている各画素について、散布図Ｓｃ２１の横軸座標に応じて画素値を割り当てることにより、基準物質画像Ｉｍ２１を生成する。また、出力機能２４ｅは、生成した基準物質画像Ｉｍ２１をディスプレイ２２に表示させる。基準物質画像Ｉｍ２１は、水成分に対応する基準物質画像であり、第２分布情報に基づく強調画像の一例である。

ここで、散布図Ｓｃ２１においては、散布図Ｓｃ１１の場合と比較して、領域Ｒ１１と領域Ｒ２１との中心間の距離Ｄが拡大している。従って、基準物質画像Ｉｍ２１においては、正常領域と浮腫領域との間のコントラストが向上しており、基準物質画像Ｉｍ２１を参照したユーザは、浮腫領域をより容易に観察することができる。

なお、出力機能２４ｅは、基準物質画像Ｉｍ２１とともに、散布図Ｓｃ１１や散布図Ｓｃ２１の表示を行なうこととしてもよい。また、生成機能２４ｄは、軸Ａ１１や軸Ａ２１といった第２軸について、ユーザによる調整を受け付けてもよい。この場合、出力機能２４ｅは、再設定された第２軸に基づいて基準物質画像Ｉｍ２１を順次更新しながらディスプレイ２２に表示させることができる。

また、図４Ａ〜図４Ｃにおいては、散布図Ｓｃ１１において設定された複数の領域のうち領域Ｒ１１が選択され、領域Ｒ１１との関係において第２軸が設定されるものとして説明した。しかしながら実施形態はこれに限定されるものではなく、生成機能２４ｄは、複数の領域から領域Ｒ２１を選択し、領域Ｒ２１との関係において第２軸を設定することとしてもよい。

また、第２軸の設定方法は、図４Ａ〜図４Ｃに示した手法に限定されるものではない。以下、第２軸の設定方法の変形例について、図５、図６及び図７を用いて説明する。図５、図６及び図７は、第１の実施形態に係る第２軸の設定方法の一例を示す図である。

図５においては、散布図Ｓｃ１１において領域Ｒ１２及び領域Ｒ２２が設定されているものとして説明する。領域Ｒ１２及び領域Ｒ２２は、図４Ａに示した領域Ｒ１１及び領域Ｒ２１と同様にして設定することができる。図５に示す場合、生成機能２４ｄは、まず、領域Ｒ１２及び領域Ｒ２２のいずれかを選択する。以下では、領域Ｒ１２を選択したものとして説明する。

次に、生成機能２４ｄは、選択した領域Ｒ１２の長軸に沿って軸Ａ１２を設定する。また、生成機能２４ｄは、軸Ａ１２と直交するように、軸Ａ２２を更に設定する。例えば、生成機能２４ｄは、軸Ａ１２と直交し、且つ、領域Ｒ１２の中心を通るように、軸Ａ２２を設定する。なお、軸Ａ１２は、第２軸及び第３軸の一例である。また、軸Ａ２２は、第２軸及び第４軸の一例である。

そして、生成機能２４ｄは、軸Ａ１２及び軸Ａ２２に応じて散布図Ｓｃ１１の変換を行ない、図５に示す散布図Ｓｃ２２を生成する。また、出力機能２４ｅは、散布図Ｓｃ２２に基づいて基準物質画像Ｉｍ２２を生成し、ディスプレイ２２に表示させる。図５に示すように、散布図Ｓｃ２２においては、散布図Ｓｃ１１の場合と比較して、領域Ｒ１２と領域Ｒ２２との中心間の距離Ｄが拡大している。従って、基準物質画像Ｉｍ２２においては、正常領域と浮腫領域との間のコントラストが向上しており、基準物質画像Ｉｍ２２を参照したユーザは、浮腫領域をより容易に観察することができる。

次に、図６について説明する。図６においては、散布図Ｓｃ１１において領域Ｒ１３及び領域Ｒ２３が設定されているものとして説明する。領域Ｒ１３及び領域Ｒ２３は、図４Ａに示した領域Ｒ１１及び領域Ｒ２１と同様にして設定することができる。

次に、生成機能２４ｄは、領域Ｒ１３及び領域Ｒ２３を繋ぐ曲線を生成する。ここで、生成機能２４ｄは、入力インタフェース２３を介してユーザからの入力操作を受け付けることにより曲線を生成してもよいし、散布図Ｓｃ１１を画像処理して曲線を自動生成してもよい。例えば、生成機能２４ｄは、領域Ｒ２３において長軸と円周が交差する点Ｘ１と、領域Ｒ２３の中心と、領域Ｒ１３の中心との３点を最短且つ連続的に結ぶように、曲線を生成する。なお、生成機能２４ｄは、端部（点Ｘ１、領域Ｒ１３の中心）の曲率を維持したまま曲線を延長することで、外側領域まで曲線を外挿することができる。また、生成機能２４ｄは、自動生成した曲線について、ユーザによる調整を受け付けることとしてもよい。そして、生成機能２４ｄは、生成した曲線に沿って、軸Ａ１３を設定する。

更に、生成機能２４ｄは、軸Ａ１３と交差するように軸Ａ２３を設定する。例えば、生成機能２４ｄは、軸Ａ２３として、図６に示す軸Ａ２３１及び軸Ａ２３２を設定する。一例を挙げると、生成機能２４ｄは、軸Ａ１３と直交し且つ領域Ｒ２３の中心を通るように、軸Ａ２３１を設定する。また、生成機能２４ｄは、軸Ａ１３と直交し且つ領域Ｒ１３の中心を通るように、軸Ａ２３２を設定する。なお、軸Ａ１３は、第２軸及び第３軸の一例である。また、軸Ａ２３１及び軸Ａ２３２は、第２軸及び第４軸の一例である。

そして、生成機能２４ｄは、軸Ａ１３と、軸Ａ２３１及び軸Ａ２３２とに応じて散布図Ｓｃ１１の変換を行ない、図６に示す散布図Ｓｃ２３を生成する。具体的には、生成機能２４ｄは、軸Ａ１３を直線となるように変形させて散布図Ｓｃ２３の横軸とする。また、生成機能２４ｄは、散布図Ｓｃ２３の横軸と直交するように、軸Ａ２３１及び軸Ａ２３２の傾きを変化させる。更に、生成機能２４ｄは、軸Ａ１３との間の距離を維持し、且つ、軸Ａ２３１との間の距離と軸Ａ２３２との間の距離との比率を維持するように、散布図Ｓｃ１１上の各ピクセルを座標変換することで、散布図Ｓｃ２３を生成する。なお、かかる変換処理を行なう場合には、散布図Ｓｃ１１において、軸Ａ２３１と軸Ａ２３２との交点の近傍領域が大幅に拡大されることとなる。そこで、生成機能２４ｄは、拡大率についての閾値を設け、閾値以上に拡大される領域を散布図Ｓｃ２３から省略することとしてもよい。

また、出力機能２４ｅは、散布図Ｓｃ２３に基づいて基準物質画像Ｉｍ２３を生成し、ディスプレイ２２に表示させる。図６に示すように、散布図Ｓｃ２３においては、散布図Ｓｃ１１の場合と比較して、領域Ｒ１３と領域Ｒ２３との中心間の距離Ｄが拡大している。従って、基準物質画像Ｉｍ２３においては、正常領域と浮腫領域との間のコントラストが向上しており、基準物質画像Ｉｍ２３を参照したユーザは、浮腫領域をより容易に観察することができる。

なお、図６においては、散布図Ｓｃ１１が２つの領域（領域Ｒ１３及び領域Ｒ２３）を含むものとして説明した。即ち、図６においては、２物質（例えば正常領域及び浮腫領域など）のケースについて説明した。しかしながら実施形態はこれに限定されるものではなく、例えば、正常領域及び浮腫領域に対して更に物質が追加されてもよい。一例を挙げると、散布図Ｓｃ１１は、正常領域に対応した領域Ｒ１３及び浮腫領域に対応した領域Ｒ２３に加えて、第３の物質に対応した領域Ｒ４を含むこととしてもよい。かかる場合においても、生成機能２４ｄは、図６に示した場合と同様に第２軸の設定を行なうことができる。例えば、生成機能２４ｄは、領域Ｒ１３の中心、領域Ｒ２３の中心、及び、領域Ｒ４の中心の３点を通る曲線を生成し、この曲線に沿って第２軸の設定を行なうことができる。

次に、図７について説明する。図７においては、関心領域内の各画素のプロットが領域Ｒ１４に分布しているものとして説明する。即ち、これまでは、関心領域内の各画素のプロットが散布図Ｓｃ１１において楕円形に分布するものとして説明したが、図７においては、プロットがより複雑に分布している場合について説明する。

なお、正常領域と浮腫領域とを含む関心領域が設定されている場合、散布図Ｓｃ１１には、正常領域及び浮腫領域それぞれに対応する２つの領域が現れる。即ち、散布図Ｓｃ１１には、図７に示す領域Ｒ１４の他、もう一つの領域が現れる。但し、図７においては、説明の便宜のため一方の領域を省略し、領域Ｒ１４についてのみ説明する。

例えば、生成機能２４ｄは、領域Ｒ１４に対して所定形状のモデルを当てはめる。例えば、生成機能２４ｄは、領域Ｒ１４に対して楕円形状のモデルを当てはめる。以下では、楕円形状のモデルが当てはめられた領域を、領域Ｒ１４’と記載する。

ここで、図７においては領域Ｒ１４を単一のパターンで示すが、領域Ｒ１４内におけるプロットの密度は均一ではない。そこで、生成機能２４ｄは、領域Ｒ１４内におけるプロットの密度を考慮して、モデルの当てはめを行なうこととしてもよい。例えば、生成機能２４ｄは、領域Ｒ１４を所定サイズの小領域に分割し、各小領域内におけるプロットの密度が閾値を超えるか否かを判定し、密度が閾値を超える小領域に対してモデルを当てはめる。

そして、生成機能２４ｄは、領域Ｒ１４’に基づいて第２軸の設定を行なう。例えば、生成機能２４ｄは、楕円形状の領域Ｒ１４’ の長軸に沿って軸Ａ１４を設定し、軸Ａ１４と直交するように軸Ａ２４を設定する。なお、軸Ａ１４は、第２軸及び第３軸の一例である。また、軸Ａ２４は、第２軸及び第４軸の一例である。そして、生成機能２４ｄは、軸Ａ１４及び軸Ａ２４に応じて散布図Ｓｃ１１の変換を行ない、図７に示す散布図Ｓｃ２４を生成する。また、出力機能２４ｅは、散布図Ｓｃ２４に基づいて基準物質画像Ｉｍ２４を生成し、ディスプレイ２２に表示させる。

なお、図７においては、第２軸として軸Ａ１４及び軸Ａ２４を設定する場合について説明したが、実施形態はこれに限定されるものではない。例えば、楕円形状の領域Ｒ１４’を設定した後、生成機能２４ｄは、図４Ａ〜図４Ｃに示した方法や、図６に示した方法と同様にして、第２軸の設定を行なうこともできる。

また、これまで、散布図Ｓｃ１１を取得機能２４ｃが生成するものとして説明したが、実施形態はこれに限定されるものではない。例えば、取得機能２４ｃは、Ｘ線ＣＴ装置１０等の他の装置において生成された散布図Ｓｃ１１を、ネットワークＮＷを介して取得することとしてもよい。即ち、取得機能２４ｃは、第１分布情報を自ら生成することで取得してもよいし、他の装置において生成された第１分布情報を取得してもよい。

次に、医用情報処理装置２０による処理の手順の一例を、図８を用いて説明する。図８は、第１の実施形態に係る医用情報処理装置２０の処理の一連の流れを説明するためのフローチャートである。ステップＳ１０１は、指定機能２４ｂに対応する。ステップＳ１０２は、取得機能２４ｃに対応する。ステップＳ１０３、ステップＳ１０４及びステップＳ１０７は、生成機能２４ｄに対応する。ステップＳ１０５及びステップＳ１０６は、出力機能２４ｅに対応する。

まず、処理回路２４は、被検体Ｐにおける関心領域を指定する（ステップＳ１０１）。次に、処理回路２４は、関心領域内の物質の分布を示す第１分布情報を取得する（ステップＳ１０２）。次に、処理回路２４は、第１分布情報において複数の領域を設定し、設定した複数の領域に基づいて第２軸を設定する（ステップＳ１０３）。次に、処理回路２４は、第２軸に応じて第１分布情報の変換を行なうことで、第２分布情報を生成する（ステップＳ１０４）。

次に、処理回路２４は、第２分布情報に基づく強調画像を生成し（ステップＳ１０５）、生成した強調画像をディスプレイ２２に表示させる（ステップＳ１０６）。ここで、処理回路２４は、第２軸の再設定を行なうか否かを判定する（ステップＳ１０７）。例えば、強調画像を参照してコントラストが不十分と判断した場合、ユーザは、入力インタフェース２３を介して、第２軸の再設定を行なう旨の指示を入力することができる。

第２軸の再設定を行なう場合（ステップＳ１０７肯定）、処理回路２４は、再度ステップＳ１０３に移行する。ここで、処理回路２４は、第２軸の再設定を自動で行なってもよいし、ユーザからの入力操作に基づいて第２軸の再設定を行なってもよい。例えば、図４Ａ〜図４Ｃに示した手法により第２軸の設定を行なっていた場合において、処理回路２４は、図５に示した手法により第２軸の再設定を行なうことができる。一方で、第２軸の再設定を行なわない場合（ステップＳ１０７否定）、処理回路２４は、処理を終了する。

上述したように、第１の実施形態によれば、取得機能２４ｃは、異なるＸ線エネルギーに対応する複数のデータセットに基づいて生成された、複数の指標に応じた複数の第１軸を有する座標系における物質の分布を示す第１分布情報を取得する。また、生成機能２４ｄは、第１分布情報において設定した複数の領域に基づいて、複数の第１軸とは異なる複数の第２軸を第１分布情報において設定する。また、生成機能２４ｄは、複数の第２軸に応じて第１分布情報の変換を行なうことで、複数の第２軸を有する座標系における物質の分布を示す第２分布情報を生成する。また、出力機能２４ｅは、第２分布情報に基づく強調画像を出力する。即ち、医用情報処理装置２０は、物質の分布について回帰分析を行ない、軸を再定義して、物質分離能を向上させた強調画像を出力する。

これにより、第１の実施形態に係る医用情報処理装置２０は、異なるＸ線エネルギーに対応する複数のデータセットを用いた物質弁別の処理において、物質分離能を向上させることができる。ひいては、医用情報処理装置２０は、病変等の観察を容易にし、診断能を向上させることができる。

（第２の実施形態）
さて、これまで第１の実施形態について説明したが、上述した実施形態以外にも、種々の異なる形態にて実施されてよいものである。

例えば、上述した実施形態では、生成された強調画像をディスプレイ２２に表示させるものとして説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、出力機能２４ｅは、生成された強調画像をＸ線ＣＴ装置１０等の外部装置に送信し、当該外部装置において強調画像の表示が行なわれることとしてもよい。

また、上述した実施形態では、図４Ａ〜図７等に示した方法により、第２軸を生成機能２４ｄが自動設定するものとして説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、生成機能２４ｄは、散布図Ｓｃ１１を参照したユーザからの入力操作に基づいて、第２軸の設定を行なうこととしても構わない。

また、上述した実施形態では、第１分布情報として、関心領域内の物質の分布を示す分布情報を取得するものとして説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、関心領域を指定する処理については省略し、第１分布情報として、被検体Ｐから収集されたデータセットの全領域における物質の分布を示す分布情報を取得することとしても構わない。また、この場合、処理回路２４は、指定機能２４ｂを有しないこととしても構わない。

また、上述した実施形態では、第２分布情報に基づく強調画像の例として、水成分に対応する基準物質画像について説明した。しかしながら、実施形態はこれに限定されるものではなく、カルシウム成分に対応する基準物質画像についても同様に物質分離能を向上させることが可能である。また、基準物質の選択は任意であり、水及びカルシウム以外の基準物質を選択する場合にも同様に適用が可能である。

また、上述した実施形態では、強調画像の例として基準物質画像について説明したが、実施形態はこれに限定されるものではない。例えば、出力機能２４ｅは、水成分に対応する基準物質画像とカルシウム成分に対応する基準物質画像とに基づいて重み付け計算処理を行うことにより、強調画像として、仮想単色Ｘ線画像（モノクロマティック画像）や密度画像、実効原子番号画像等を生成することも可能である。医用情報処理装置２０は、基準物質画像において物質分離能を向上させることにより、基準物質画像に基づく各種の強調画像についても物質分離能を向上させることができる。

また、上述した実施形態では、第１指標及び第２指標の例として、水の存在量に関する指標及びカルシウムの存在量に関する指標について説明した。即ち、上述した実施形態では、第１指標及び第２指標の例として、第１物質の存在量に関する指標及び第２物質の存在量に関する指標について説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、第１指標及び第２指標は、第１のＸ線エネルギーに関する指標及び第２のＸ線エネルギーに関する指標であってもよい。即ち、図３に示した散布図Ｓｃ１１の横軸及び縦軸が、それぞれ、第１のＸ線エネルギーに関する第１指標に応じた第１指標軸及び第２のＸ線エネルギーに関する第２指標に応じた第２指標軸となる場合であってもよい。

また、上述した実施形態では、デュアルエナジーを例として説明したが、３種類以上のＸ線エネルギーを用いたマルチエナジーについても同様に適用が可能である。例えば、３種類のＸ線エネルギーを用いてデータセットの収集が行われた場合、取得機能２４ｃは、第１指標に応じた第１指標軸、第２指標に応じた第２指標軸、及び、第３指標に応じた第３指標軸の３軸を有する座標系における物質の分布を示す第１分布情報を取得することができる。また、この場合、生成機能２４ｃは、３つの第２軸を第１分布情報において設定し、設定した第２軸に応じて第１分布情報の変換を行なうことで、第２分布情報を生成する。そして、出力機能２４ｅは、第２分布情報に基づく強調画像を出力する。

また、上述した実施形態では、Ｘ線ＣＴ装置１０と医用情報処理装置２０とが別体であるものとして説明したが、医用情報処理装置２０はＸ線ＣＴ装置１０に含まれることとしてもよい。例えば、Ｘ線ＣＴ装置１０におけるコンソール装置１４０が、医用情報処理装置２０としての機能を実行する場合であってもよい。

また、上述した実施形態では、異なるＸ線エネルギーに対応する複数のデータセットの例として、Ｘ線ＣＴ装置１０により収集された投影データについて説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、被検体Ｐを透過したＸ線を検出するＦＰＤ（Ｆｌａｔ Ｐａｎｅｌ Ｄｅｔｅｃｔｏｒ：平面検出器）を備えたＸ線診断装置により収集された、異なるＸ線エネルギーに対応する複数のＸ線画像についても同様に適用が可能である。

上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ：ＦＰＧＡ））等の回路を意味する。プロセッサが例えばＣＰＵである場合、プロセッサは記憶回路に保存されたプログラムを読み出し実行することで機能を実現する。一方、プロセッサが例えばＡＳＩＣである場合、記憶回路にプログラムを保存する代わりに、当該機能がプロセッサの回路内に論理回路として直接組み込まれる。なお、実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。さらに、各図における複数の構成要素を１つのプロセッサへ統合してその機能を実現するようにしてもよい。

また、図１においては、単一のメモリ２１が処理回路２４の各処理機能に対応するプログラムを記憶するものとして説明した。また、図２においては、単一のメモリ１４１が処理回路１４４の各処理機能に対応するプログラムを記憶するものとして説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、複数のメモリ２１を分散して配置し、処理回路２４は、個別のメモリ２１から対応するプログラムを読み出す構成としても構わない。同様に、複数のメモリ１４１を分散して配置し、処理回路１４４は、個別のメモリ１４１から対応するプログラムを読み出す構成としても構わない。また、メモリ２１又はメモリ１４１にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。

上述した実施形態に係る各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。即ち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行われる各処理機能は、その全部又は任意の一部が、ＣＰＵ及び当該ＣＰＵにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現されうる。

また、上述した実施形態で説明した医用情報処理方法は、予め用意された医用情報処理プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この医用情報処理プログラムは、インターネット等のネットワークを介して配布することができる。また、この医用情報処理プログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な非一過性の記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。

以上説明した少なくとも１つの実施形態によれば、物質分離能を向上させることができる。

いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行なうことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１ 医用情報処理システム
１０ Ｘ線ＣＴ装置
１１０ 架台装置
１３０ 寝台装置
１４０ コンソール装置
１４４ 処理回路
１４４ａ 制御機能
１４４ｂ 収集機能
１４４ｃ 出力機能
２０ 医用情報処理装置
２４ 処理回路
２４ａ 制御機能
２４ｂ 指定機能
２４ｃ 取得機能
２４ｄ 生成機能
２４ｅ 出力機能

Claims (12)

1. 異なるＸ線エネルギーに対応する複数のデータセットに基づいて生成された、複数の指標に応じた複数の第１軸を有する座標系における物質の分布を示す第１分布情報を取得する取得部と、
   前記第１分布情報において設定した複数の領域に基づいて前記複数の第１軸とは異なる複数の第２軸を前記第１分布情報において設定し、前記複数の第２軸に応じて前記第１分布情報の変換を行なうことで、前記複数の第２軸を有する座標系における物質の分布を示す第２分布情報を生成する生成部と、
   前記第２分布情報に基づく強調画像を出力する出力部と
   を備える、医用情報処理装置。
2. 前記取得部は、前記第１分布情報として、第１指標に応じた第１指標軸及び第２指標に応じた第２指標軸を有する座標系における物質の分布を示す分布情報を取得し、
   前記生成部は、前記複数の第２軸として、前記第１指標軸及び前記第２指標軸とは異なる第３軸及び第４軸を前記第１分布情報において設定し、前記第２分布情報として、前記第３軸及び前記第４軸を有する座標系における物質の分布を示す分布情報を生成する、請求項１に記載の医用情報処理装置。
3. 前記第１指標は、第１物質の存在量に関する指標であり、
   前記第２指標は、第２物質の存在量に関する指標である、請求項２に記載の医用情報処理装置。
4. 前記生成部は、前記第１物質の存在量を変化させた複数の対象物について生成された、前記第１指標軸及び前記第２指標軸を有する座標系における物質の分布を示す分布情報を取得し、前記複数の対象物それぞれに対応する複数の対応領域を前記分布情報において設定し、設定した複数の前記対応領域に基づいて、前記複数の第２軸に含まれる軸の１つを設定する、請求項３に記載の医用情報処理装置。
5. 前記生成部は、前記複数の領域から１つの領域を選択し、選択した領域の長軸に沿って、前記複数の第２軸に含まれる軸の１つを設定する、請求項１〜４のいずれか一項に記載の医用情報処理装置。
6. 前記生成部は、前記複数の領域から１つの領域を選択し、選択した領域の長軸に沿って、前記複数の第２軸に含まれる軸の１つを設定し、前記長軸に沿って設定した軸と直交するように、前記複数の第２軸に含まれる軸の１つを更に設定する、請求項１〜３のいずれか一項に記載の医用情報処理装置。
7. 前記生成部は、前記複数の領域を繋ぐ曲線を生成し、当該曲線に沿って、前記複数の第２軸に含まれる軸の１つを設定する、請求項１〜３のいずれか一項に記載の医用情報処理装置。
8. 前記生成部は、記複数の領域から１つの領域を選択し、選択した領域に対して所定形状のモデルを当てはめ、当てはめたモデルに基づいて、前記複数の第２軸に含まれる軸の１つを設定する、請求項１〜３のいずれか一項に記載の医用情報処理装置。
9. 前記生成部は、前記複数の領域の間の分離能が高くなるように、前記複数の第２軸を設定する、請求項１〜８のいずれか一項に記載の医用情報処理装置。
10. 被検体における関心領域を指定する指定部を更に備え、
    前記取得部は、前記第１分布情報として、前記複数の第１軸を有する座標系における、前記関心領域内の物質の分布を示す分布情報を取得する、請求項１〜９のいずれか一項に記載の医用情報処理装置。
11. 異なるＸ線エネルギーに対応する複数のデータセットに基づいて生成された、複数の指標に応じた複数の第１軸を有する座標系における物質の分布を示す第１分布情報を取得し、
    前記第１分布情報において設定した複数の領域に基づいて前記複数の第１軸とは異なる複数の第２軸を前記第１分布情報において設定し、前記複数の第２軸に応じて前記第１分布情報の変換を行なうことで、前記複数の第２軸を有する座標系における物質の分布を示す第２分布情報を生成し、
    前記第２分布情報に基づく強調画像を出力する
    ことを含む、医用情報処理方法。
12. 異なるＸ線エネルギーに対応する複数のデータセットに基づいて生成された、複数の指標に応じた複数の第１軸を有する座標系における物質の分布を示す第１分布情報を取得し、
    前記第１分布情報において設定した複数の領域に基づいて前記複数の第１軸とは異なる複数の第２軸を前記第１分布情報において設定し、前記複数の第２軸に応じて前記第１分布情報の変換を行なうことで、前記複数の第２軸を有する座標系における物質の分布を示す第２分布情報を生成し、
    前記第２分布情報に基づく強調画像を出力する
    各処理をコンピュータに実行させる医用情報処理プログラム。

Images