JP6338849B2 - 画像解析装置、及び、ｘ線診断装置 - Google Patents

Description

本発明の実施形態は、画像解析装置、及び、Ｘ線診断装置に関する。

血管内治療として、カテーテル等により血管を内側から広げる血管拡張手術と、腫瘍に栄養を送る血管を塞ぐことで腫瘍の壊死を図る血管塞栓手術とが知られている。このような血管治療の効果を確認するために、ＤＳＡ(Digital Subtraction Angiography)画像が用いられる（例えば、特許文献１参照）。

ＤＳＡ画像は、例えばＸ線診断装置により、造影剤投与の前後においてそれぞれ被検体の同一領域を時系列的に撮影することで得られる。具体的には、造影剤投与後の各時相の画像から造影剤投与前のマスク画像をそれぞれ差し引くことで得られる各時相に対応した複数の差分画像がＤＳＡ画像とされる。

従来技術では、上記のように時系列的なＤＳＡ画像を手術前、手術後のそれぞれにおいて生成し、読影者がこれらＤＳＡ画像を目視で比較することで手術の効果を確認していた。ここでの手術の効果の確認とは例えば、血管拡張手術の場合には血流改善領域の判定であり、血管塞栓手術の場合には血流閉塞領域の判定である。

特開２００６−８７６３１号公報

手術前のＤＳＡ画像の動画と、手術後のＤＳＡ画像の動画とを比較しても、血流の改善具合又は血流閉塞領域が分かりにくい場合、読影者は、手術前の複数のＤＳＡ画像と、手術後の複数のＤＳＡ画像とから任意の１フレームずつを選択し、静止画として比較していた。このように手術前と手術後のＤＳＡ画像同士を目視で比較しても、血流の僅かな差異を特定することは、必ずしも容易ではない。

このため、手術前後のＤＳＡ画像に基づいて、血管拡張手術後の血流改善領域や血管塞栓手術後の血流閉塞領域等の手術後の血流量変化領域を従来よりも容易に視認可能にする新技術が要望されていた。

以下、本発明の実施形態が取り得る態様の数例を、態様毎に説明する。
（１）一実施形態では、画像解析装置は、ＤＳＡ画像取得部と、パラメータ値取得部と、カラー割当て部と、時相画像生成部とを備える。
ＤＳＡ画像取得部は、同一の被検体に対する時系列的な複数のＤＳＡ画像の画像データを取得する。
パラメータ値取得部は、複数のＤＳＡ画像の画像データにおける、被検体の同一領域に対応する画素毎の画素値の時間変化に基づいて、画素毎にパラメータ値を取得する。
カラー割当て部は、パラメータ値に応じた色が被検体の同一領域に対応する画素毎に割り当てられように、且つ、複数の画素に有彩色が割り当てられるように、カラーマップを生成する。
時相画像生成部は、カラーマップにおける画素毎の色及びパラメータ値に応じた色変換をカラーマップの各画素に施すことで、互いに異なる複数の時相にそれぞれ対応する複数の時相画像の各カラー画像データを生成する。

（２）別の一実施形態では、Ｘ線診断装置は、Ｘ線撮影部と、上記（１）の画像解析装置とを備える。
Ｘ線撮影部は、造影剤投与の前後において被検体を透過したＸ線を検出することでＸ線画像の投影データを生成し、造影剤投与後の時系列的な複数のＸ線画像と、造影剤投与前のＸ線画像との各差分に基づいて被検体に対する時系列的な複数のＤＳＡ画像の画像データを生成する。
画像解析装置は、複数のＤＳＡ画像の画像データにおける、被検体の同一領域に対応する画素毎の画素値の時間変化に基づいて、画素毎にパラメータ値を取得し、パラメータ値に応じた色が画素毎に割り当てられようにカラーマップを生成し、カラーマップに色変換を施すことで、互いに異なる複数の時相にそれぞれ対応する複数の時相画像の各カラー画像データを生成する。

本実施形態におけるＸ線診断装置の構成の一例を示すブロック図。 パラメトリック画像の生成過程におけるパラメータ値の取得方法の一例として、造影剤濃度の時間変化の算出方法を示す模式図。 ＴＴＰをパラメータとしたパラメトリック画像の生成方法の一例を示す模式図。 時相画像生成部によるパラメトリック時相画像の生成方法の第１の例を段階的に示す模式図。 時相画像生成部によるパラメトリック時相画像の生成方法の第２の例を段階的に示す模式図。 血管拡張手術等の後における、特定処理部による血流改善領域の判定方法を段階的に示す模式図。 パラメトリック時相画像の動画表示において、血流改善領域を枠で囲うことで識別表示する一例を示す模式図。 パラメトリック時相画像の動画表示において、血流改善領域を点滅させることで識別表示する一例を示す模式図。 パラメトリック時相画像の動画表示において、特定処理部により抽出される血管壁の強調表示の一例を示す模式図 パラメトリック時相画像の動画表示において、血流改善領域のみをＤＳＡ画像によりグレースケールの動画で表示しつつ、残りの領域をパラメトリック時相画像によりカラー表示する一例を示す模式図。 血管拡張手術の前後において、本実施形態のＸ線診断装置により造影剤投与を伴った撮影及び血流改善領域の表示が行われる場合の動作の一例を示すフローチャート。 血管塞栓手術の前後において、本実施形態のＸ線診断装置により造影剤投与を伴った撮影及び血流閉塞領域の表示が行われる場合の動作の一例を示すフローチャート。

以下、本発明の実施形態を添付図面に基づいて説明する。なお、各図において同一要素には同一符号を付し、重複する説明を省略する。

＜本実施形態の構成＞
図１は、本実施形態におけるＸ線診断装置１０の構成の一例を示すブロック図である。図１に示すように、Ｘ線診断装置１０は、入力部２２と、画像解析装置２４と、システム制御部２６と、投影データ記憶部２８と、投影データ生成部３０と、Ｘ線検出器３４と、Ｃアーム３６と、天板３８と、絞り装置４０と、Ｘ線管４２と、高電圧発生器４４と、絞り制御機構４６と、天板移動機構４８と、Ｃアーム動作機構５０と、検出器移動機構５４とを有する。

Ｘ線診断装置１０の主な特徴は画像解析装置２４の機能にあるため、以下、他の構成要素の機能を簡単に説明後、画像解析装置２４の機能について詳細に説明する。
天板３８上には被検体Ｐが載置される。
Ｃアーム３６は、Ｘ線管４２、絞り装置４０、Ｘ線検出器３４を保持するアームである。Ｃアーム３６によって、Ｘ線管４２及び絞り装置４０と、Ｘ線検出器３４とは、被検体Ｐを挟んで互いに対向するように配置される。

高電圧発生器４４は、高電圧を発生して、発生した高電圧をＸ線管４２に供給する。
Ｘ線管４２は、高電圧発生器４４から供給される高電圧を用いてＸ線を発生する。
絞り装置４０は、例えば複数の絞り羽根をスライドさせることで、発生したＸ線が被検体Ｐの撮影領域に対して選択的に照射されるように絞り込む。

絞り制御機構４６は、絞り装置４０の絞り羽根の開度を調整することで、Ｘ線の照射範囲を制御する。
Ｘ線検出器３４は、例えばマトリクス状に配列された多数のＸ線検出素子によって、被検体Ｐを透過したＸ線を電気信号に変換して蓄積し、蓄積した電気信号を投影データ生成部３０に入力する。

投影データ生成部３０は、Ｘ線検出器３４から入力される電気信号を用いて投影データを生成し、生成した投影データを投影データ記憶部２８に記憶させる。投影データは例えば、各画素が１つの画素値を有するＸ線画像の画像データであり、投影データの各画素値は、被検体の撮影領域のＸ線透過率を反映したものとなる。

ここでは一例として、投影データ生成部３０は、投影データに基づいてＤＳＡ画像の画像データも生成して投影データ記憶部２８に記憶させるものとする。
但し、画像解析装置２４のＤＳＡ画像取得部２４ａがＸ線画像の投影データを取得して、ＤＳＡ画像取得部２４ａが投影データに基づいてＤＳＡ画像の画像データを生成する構成でもよい。

入力部２２は、操作者が撮影条件や画像表示条件等の各種コマンドを入力するためのキーボード、ボタンなどを有し、入力された内容をシステム制御部２６に転送する。
システム制御部２６は、撮影条件の設定、撮影動作、表示処理においてＸ線診断装置１０全体を制御する。

次に、画像解析装置２４の各構成要素について説明する。
画像解析装置２４は、システムバスＳＢと、ＤＳＡ画像取得部２４ａと、パラメータ値取得部２４ｂと、カラー割当て部２４ｃと、時相画像生成部２４ｄと、特定処理部２４ｅと、表示部２４ｆとを有する。

システムバスＳＢは、画像解析装置２４の各構成要素を互いに電気的に接続する通信配線である。
ＤＳＡ画像取得部２４ａは、造影剤投与の前後の同一被検体に対する撮影により得られた時系列的な複数のＤＳＡ画像の画像データを投影データ記憶部２８から取得する。

パラメータ値取得部２４ｂは、各ＤＳＡ画像の画像データにおける、被検体の同一領域に対応する画素毎の画素値の時間変化に基づいて、パラメトリックイメージング用のパラメータ値を画素毎に取得する（後述の図２参照）。

上記パラメトリックイメージング(Parametric Imaging)とは、例えば、単数又は複数のパラメータをカラー又はグレースケールで画像化する処理である。パラメトリックイメージングは、広義には、投影データ生成部３０により生成されるＸ線画像の投影データも含む。Ｘ線画像の投影データの場合、各画素の画素値は、パラメータとしてのＸ線透過率の値を示すと言えるからである。

パラメトリックイメージングは、狭義には、Ｘ線画像の投影データを元にして、Ｘ線透過率以外のパラメータの値を画素毎に算出し、カラー画像にすることである。本実施形態では、狭義のパラメトリックイメージングについて、図２以降で説明する。また、以下の説明では、狭義のパラメトリックイメージングにより生成される画像を「パラメトリック画像」と称する。

カラー割当て部２４ｃは、パラメータ値に応じた色が被検体Ｐの同一領域に対応する画素毎に割り当てられように、且つ、複数の画素に有彩色が割り当てられるように、パラメトリック画像の画像データ（カラーマップ）を生成する。パラメトリック画像の画像データは、各画素が例えば赤、緑、青の３原色に対してそれぞれ画素値を有する画像データである。ここでは一例として、一連の時系列的なＤＳＡ画像のセットに対して、一のパラメトリック画像の画像データが生成されるものとする。

時相画像生成部２４ｄは、一のパラメトリック画像の画素毎に色変換を施すことで、複数の時相にそれぞれ対応する複数の「パラメトリック時相画像」のカラー画像データを生成する。パラメトリック時相画像の詳細については、後述する。また、ここでの複数の時相とは、例えば、パラメトリック画像の生成元の各ＤＳＡ画像の時相である。

特定処理部２４ｅは、手術前後の各パラメトリック画像同士の差分に基づいて、血流量変化領域を特定する。ここでの血流量変化領域とは例えば、血管拡張手術の場合には血流改善領域あり、血管塞栓手術の場合には血流閉塞領域である。

表示部（モニタ）２４ｆは、入力部２２を介して操作者からコマンドを受け付けるためのＧＵＩ(Graphical User Interface)や、Ｘ線画像などを表示する。また、本実施形態の主な特徴の１つとして、表示部２４ｆは、上記血流量変化領域が識別される複数の時相にそれぞれ対応する複数の画像の各カラー画像データを動画表示する。

以下、図２〜図６を参照しながら、パラメトリック画像の生成及び血流量変化領域の抽出に至るまでの画像解析装置２４の機能について、より詳細に説明する。
図２は、パラメトリック画像の生成過程におけるパラメータ値の取得方法の一例として、造影剤濃度の時間変化の算出方法を示す模式図である。図２の上段は、各時相のＤＳＡ画像を示し、図２の下段は、一画素に着目した造影剤濃度の時間変化の一例を示す。

例えば、Ｘ線診断装置１０の撮影により、同一の被検体Ｐの同一の撮影領域に対して、造影剤投与前の時刻ｔ＝０、造影剤投与後の時刻ｔ＝１，２，３，４，５で順に６のＸ線画像の投影データが投影データ生成部３０により生成される場合を考える。この場合、造影剤投与後の各Ｘ線画像から、ｔ＝０のＸ線画像（マスク画像）をそれぞれ差し引くことで、ｔ＝１，２，３，４，５にそれぞれ対応する５コマのＤＳＡ画像（差分画像）の画像データが得られる（図２の上段参照）。なお、図２の上段では、ｔ＝１を時相１(Time Phase 1)、ｔ＝２を時相２(Time Phase 2)としている（以下同様）。

ここで、パラメータ値取得部２４ｂは、５コマのＤＳＡ画像を通して、同位置の画素毎に画素値の時相変化（ｔ＝１〜５まで）を算出することで、画素毎の造影剤濃度の時間変化を算出する。図２の下段は、各ＤＳＡ画像（この例では、画素数５×５）の左下の一画素に着目した造影剤濃度の時間変化の一例を示す。図２の下段において、縦軸は造影剤濃度(Intensity of Contrast Medium)を示し、横軸は時相（経過時刻ｔ）を示す。

より詳細には、造影剤は、Ｘ線吸収率が体組織よりも高いので、造影剤濃度が高い位置に対応するＸ線検出素子の受線量は少なく、造影剤はＸ線画像において周囲よりも暗く投影される。ＤＳＡ画像の各画素値は、（造影剤投与前の）マスク画像における同位置の画素値との差分であるから、同位置の一画素に着目し、この画素の画素値の時相変化に符号反転等の適切な処理を施せば、造影剤濃度の時間変化と等価になる。

パラメータ値取得部２４ｂは、上記のように算出される造影剤濃度の時間変化に基づいて、各画素に対して例えば以下の３つのパラメータ値の少なくとも１つを取得する。

第１には、ＴＴＰ(Time To Peak)であり、これは、どの時相で造影剤濃度がピークに達するかを示す。上記の左下の一画素の例では、時相３において造影剤濃度がピークに達するため、３というパラメータ値が与えられる。

第２には、ＰＨ(Peak Height)であり、これは、造影剤濃度のピーク値を示す。上記の左下の一画素の例では、時相３における縦軸の値に該当する。なお、縦軸の値の範囲は、例えばダイナミックレンジに合わせて任意に設定することができる。

第３には、ＡＵＣ(Area Under Curve)であり、これは、ＤＳＡ画像の最初の時相から最後の時相までの造影剤濃度の時間積分値を示す。上記の左下の一画素の例では、図２の下段において、右下がりの破線の斜線で示した領域の面積に相当する。

パラメータ値取得部２４ｂは、残りの全ての画素についても、同様に造影剤濃度の時間変化を算出することで、上記３つのパラメータ値の少なくとも１つを算出する。

図３は、ＴＴＰをパラメータとしたパラメトリック画像の生成方法の一例を示す模式図である。図３の上段は、例えば画素数５×５として、ＤＳＡ画像の複数コマを通して同位置となる画素毎に算出されたＴＴＰの値の一例を示す。

図３の中段は、カラー割当て部２４ｃに記憶されたＴＴＰ用のカラーテーブルの一例を示す。図３では一例として、ＴＴＰの値が１の画素には赤色(RED)が割り当てられ、ＴＴＰの値が２の画素には黄色(YELLOW)が割り当てられ、ＴＴＰの値が３の画素には緑色(GREEN)が割り当てられ、ＴＴＰの値が４の画素には青色(BLUE)が割り当てられ、ＴＴＰの値が５の画素には紫色(PURPLE)が割り当てられる。

カラーの割り当て方法は任意であるが、有彩色が含まれるように割り当てることが望ましい。グレーススケール表示では、血管部分を周囲から識別しづらいからである。上記のカラー割り当ての例では、赤、緑、青の３原色が用いられるが、例えば赤、緑、青の３原色の各値の内、１つの値が常に０となるようにしてもよいし、２つの値が常に０となるようにしてもよい。

ここでは便宜上、グレースケールの横方向のバーでカラーテーブルを示すが、実際には有彩色によるカラーバーとしてカラーテーブルを記憶してもよい。或いは、カラー割当て部２４ｃは、所定のビット数表示における赤、緑、青の３原色の各値のセットが各々のＴＴＰの値に対して示されるテーブルデータとして、カラーテーブルを記憶してもよい。例えば８ビットの場合、ＴＴＰ＝１に対しては、赤、緑、青の各値として(255, 0, 0)が割り当てられる。

図３の下段は、図３の上段の各画素のＴＴＰの値と、図３の中段のカラーテーブルとに従って規定される各画素の色を示す。即ち、各画素の色を図３下段に示すカラーで表示したものがＴＴＰのパラメトリック画像である。

なお、カラー割当て部２４ｃは、ＡＵＣ用のカラーテーブル、ＰＨ用のカラーテーブル等の他のパラメータ用のカラーテーブルも記憶している。上記の例では、ＴＴＰをパラメータとしてパラメトリック画像を生成しているが、前述のＡＵＣやＰＨをパラメータとして、パラメトリック画像を生成してもよい。また、パラメータの種類については、ＴＴＰ、ＡＵＣ、ＰＨに限定されるものではなく、例えば、造影剤濃度の時間変化における、ピークまで立ち上がる傾きの平均値をパラメータとしてもよい。

以上の図２、図３のように、パラメータ値取得部２４ｂ、カラー割当て部２４ｃにより、時系列的な複数のＤＳＡ画像に対して一のパラメトリック画像が生成される。
次に、本実施形態の特徴の１つとして、時相画像生成部２４ｄは、一のパラメトリック画像から、複数のＤＳＡ画像にそれぞれ対応する時相の複数画像のカラー画像データを生成する。以下の説明では、一のパラメトリック画像から生成される時系列的なカラー画像の各々を「パラメトリック時相画像」と称する。

図４は、時相画像生成部２４ｄによるパラメトリック時相画像の生成方法の第１の例を段階的に示す模式図である。
図４の最上段は、図２で示した５フレームの時系列的なＤＳＡ画像である。図４の上段の一例では、各ＤＳＡ画像を通して左下の一画素の画素値(Pixel Value)は、時相１から時相５まで順に０，１７９，２５５，１２７，２６となっており、この中の最大値は時相３の２５５である。

これらの画素値を最大値に基づいて正規化すると、時相１から順に０％、７５％，１００％，５０％，１０％となる（図４の上から２番目参照）。ここで、前述の例では、ＴＴＰをパラメータとするパラメトリック画像において左下の画素に割り当たられた色は、緑(GREEN)である（図４の上から３番目参照）。

ここでは一例として、時相画像生成部２４ｄは、パラメトリック画像の左下の画素の色に対して、ＤＳＡ画像における同位置の画素の正規化された各画素値に応じた色変換を施す。

即ち、パラメトリック画像の左下の画素の色に対し、彩度が０％（時相１のＤＳＡ画像における左下の画素の正規化画素値）になるように公知の色変換処理が施され、この色変換後の色が、時相１のパラメトリック時相画像の左下の画素の色となる。

同様に、パラメトリック画像の左下の画素の色に対し、彩度が７５％（時相２のＤＳＡ画像における当該画素の正規化画素値）になるように色変換が施され、この色変換後の色が、時相２のパラメトリック時相画像の左下の画素の色となる。

同様に、パラメトリック画像の左下の画素の色に対し、彩度が１００％（時相３のＤＳＡ画像における当該画素の正規化画素値）になるように色変換が施され、この色変換後の色が、時相３のパラメトリック時相画像の左下の画素の色となる。

同様に、パラメトリック画像の左下の画素の色に対し、彩度が５０％（時相４のＤＳＡ画像における当該画素の正規化画素値）になるように公知の色変換処理が施され、色変換後の色が、時相４のパラメトリック時相画像の左下の画素の色となる。

同様に、パラメトリック画像の左下の画素の色に対し、彩度が１０％（時相５のＤＳＡ画像における当該画素の正規化画素値）になるように色変換が施され、この色変換後の色が、時相５のパラメトリック時相画像の左下の画素の色となる。

なお、上記の例では、パラメトリック画像の左下の画素の色は、８ビットでの３原色の各値が赤0, 緑255, 青0で示され、元々彩度１００％であるため、時相３のパラメトリック時相画像の左下の画素は、パラメトリック画像の左下の画素と同じ色になる。

これにより、時相１〜時相５の各パラメトリック時相画像における、左下の画素の色が決まる。他の画素についても、上記同様の色変換処理を施すことで、時相１〜時相５の各パラメトリック時相画像の全画素の色が決まる。時相画像生成部２４ｄは、このようにして、時相１〜時相５の各パラメトリック時相画像のカラー画像データを生成する。

なお、時相画像生成部２４ｄは、彩度の代わりに、明度を上記同様に変えることで、時相１〜時相５の各パラメトリック時相画像を生成してもよい。即ち、パラメトリック画像の左下の画素の色に対し、彩度が０％，７５％，１００％，５０％，１０％になるように色変換がそれぞれ施され、この色変換後の各色が、時相１〜時相５の各パラメトリック時相画像の左下の画素の色となる。

また、時相画像生成部２４ｄは、彩度の代わりに、透過度を上記同様に変えることで、時相１〜時相５の各パラメトリック時相画像を生成してもよい。

以上、ＴＴＰをパラメータ値とするパラメトリック画像（カラーマップ）から、各ＤＳＡ画像における同位置の画素の正規化画素値に応じて彩度、明度、透過度のいずれか１つが変化するように色変換を施すことで、時相１〜時相５の各パラメトリック時相画像を生成する例を述べた。本発明の実施形態は、かかる態様に限定されるものではない。

時系列的な複数のＤＳＡ画像における同位置の画素値の時間変化から、前述のＡＵＣをパラメータ値として得られるパラメトリック画像（カラーマップ）を用いて、上記同様に時相１〜時相５の各パラメトリック時相画像を生成してもよい。

或いは、時系列的な複数のＤＳＡ画像における同位置の画素値の時間変化から、前述のＰＨをパラメータ値として得られるパラメトリック画像（カラーマップ）を用いて、上記同様に時相１〜時相５の各パラメトリック時相画像を生成してもよい。

何をパラメータとしてパラメトリック画像を生成するかで３通りあり、彩度、明度、透過度のどれが変わるように色変換するかで３通りある場合、図４の手法では、少なくとも、３×３＝９通りのパラメトリック時相画像の生成方法がある。

但し、前述のように、何をパラメータとしてパラメトリック画像を生成するかは、ＴＴＰ、ＡＵＣ、ＰＨの他にも、造影剤濃度の時間変化曲線における立ち上がりの傾き等がある。また、彩度、明度、透過度のいずれか２つが変化するように、各ＤＳＡ画像における同位置の画素の正規化画素値に応じて色変換を実行してもよい。或いは、彩度、明度、透過度の３つが変化するように、各ＤＳＡ画像における同位置の画素の正規化画素値に応じて色変換を実行してもよい。

これらを加味すると、図４の手法では、パラメトリック時相画像の生成方法は、２０通り以上ある。
このように生成される各パラメトリック時相画像を時相１から時相５まで時系列順に表示することで、動画像のように表示することができる。

図５は、時相画像生成部２４ｄによるパラメトリック時相画像の生成方法の第２の例を段階的に示す模式図である。図４と同様に、図５の上段に示す５フレームの時系列的な各ＤＳＡ画像を通して左下の一画素に着目して説明する。
左下の画素は、ＴＴＰが３であるから、図５の中段に示すように、時相３においてのみハイコントラスト画素となるように表示される。ここでのハイコントラスト画素とは、ある時相のパラメトリック時相画像内において、明度を上げて表示される画素である。

即ち、時相１、２、４、５の４つのパラメトリック時相画像において、左下の画素は、ハイコントラスト画素には該当しないから、明度を落として表示される。より詳細には、時相１、２、４、５の各パラメトリック時相画像における左下の画素の色は、パラメトリック画像（カラーマップ）で割り当てられた緑色に対して明度が下がるように色変換処理が施された色となる。明度の下げ方については、一律的に同じ明度としてもよいし、正規化された画素値（図４の上から２段目に参照）が低い程、明度が下がるようにしてもよい。

一方、時相３のパラメトリック時相画像における左下の画素の色は、パラメトリック画像（カラーマップ）で割り当てられた緑色に対して明度が上がるように色変換処理が施された色となる。これにより、時相１〜５の各パラメトリック時相画像における左下の画素の色がそれぞれ決まる。図５の下段に示すように、時相画像生成部２４ｄは、他の画素についても、同様の色変換処理を施すことで、時相１〜５の各パラメトリック時相画像の全画素の色を算出する。

このようにハイコントラスト画素は明るく、それ以外の画素は暗くして時系列順にパラメトリック時相画像を表示すれば、動画像のように見せることができる。即ち、ＴＴＰを反映した色合いとして算出される時系列順のパラメトリック時相画像は、造影剤の到達時間情報を反映した動画像となる。

なお、以上の説明では簡単化のため、ＤＳＡ画像、パラメトリック画像、パラメトリック時相画像の画素数が２５の例を述べたが、これらの画素数は２６以上でもよいし、２４以下でもよい。また、ＤＳＡ画像、パラメトリック画像、パラメトリック時相画像の３者で縦横の画素数が同一である必要はない。ＤＳＡ画像の画素数が極めて多い場合、パラメトリック画像、及び、パラメトリック時相画像の画素数は、例えば、ＤＳＡ画像の画素数を縦方向に半分、横方向に半分に間引いてもよい。
また、以上の説明では簡単化のため、５時相の例で説明したが、時相は６以上でもよい。

図６は、血管拡張手術等の後における、特定処理部２４ｅによる血流改善領域の判定方法を段階的に示す模式図である。前述の例と同様に、手術前において、被検体Ｐの同一の撮影領域から５フレームの手術前ＤＳＡ画像１、２、３、４、５が生成され、手術後において同様に、手術前と同一の被検体Ｐの同一の撮影領域から５フレームの手術後ＤＳＡ画像１’、２’、３’、４’、５’が生成された場合を考える。

この場合、パラメータ値取得部２４ｂ及びカラー割当て部２４ｃは、前述と同様にして、手術前ＤＳＡ画像１〜５からＰＨをパラメータとして、手術前パラメトリック画像を１つ生成する（図６の上段の左側参照）。同様に、パラメータ値取得部２４ｂ及びカラー割当て部２４ｃは、手術後ＤＳＡ画像１’〜５’からＰＨをパラメータとして、手術後のパラメトリック画像を１つ生成する（図６の上段の右側参照）。

次に、特定処理部２４ｅは、手術前後で同位置に該当する画素同士で、手術後のパラメトリック画像の各画素値から、手術前のパラメトリック画像の各画素値を差し引くことで、差分画像を生成する。特定処理部２４ｅは、この差分画像において、画素値が閾値以上の画素領域を血流改善領域として抽出する（図６の中段参照）。

なお、ＰＨの代わりに、ＡＵＣをパラメータとして手術前、手術後の各パラメトリック画像を生成し、同様に差分画像に閾値処理を施すことでも、血流改善領域を抽出（算出）できる。

さらに、特定処理部２４ｅは、血流改善領域内の血管の輪郭（血管壁）が強調された血管壁強調画像データを生成する。具体的には例えば、特定処理部２４ｅは、上記の差分画像において抽出した血流改善領域に対してラプラシアンフィルターなどの輪郭強調処理を適用することで、血管壁強調画像データを生成できる（図６の下段参照）。

一方、血管塞栓手術等の後における血流閉塞領域は、差分画像生成時の引かれる側を逆にすることを除き、同様に判定される。即ち、特定処理部２４ｅは、手術前のパラメトリック画像の各画素値から、手術後パラメトリック画像の各画素値を差し引くことで差分画像を生成し、この差分画像において画素値が閾値以上の画素領域を血流閉塞領域として算出する。

＜本実施形態の表示方式＞
図７は、パラメトリック時相画像の動画表示において、血流改善領域を枠で囲うことで識別表示する一例を示す模式図である。ここでは一例として、表示部２４ｆの画面左側において手術前のＤＳＡ画像の各時相にそれぞれ対応する各パラメトリック時相画像が時系列順に表示され、画面右側において手術後のＤＳＡ画像の各時相にそれぞれ対応する各パラメトリック時相画像が時系列順に表示される。この点、以下の図８〜図１０も同様である。

図７に示すように、表示部２４ｆは、手術前、手術後の双方のパラメトリック時相画像の動画表示において、特定処理部２４ｅにより抽出された血流改善領域を長方形の枠で囲われた関心領域(ROI: Region of Interest)３００として識別表示することができる。

なお、図７では煩雑となるので時相３までを図示しているが、実際には時相４以降も表示される（以下の図８〜図１０も同様）。また、対比のために、手術前のパラメトリック時相画像も動画表示しているが、血流改善領域を識別しつつ手術後のパラメトリック時相画像のみを表示してもよい（以下の図８〜図１０も同様）。

また、図７では、手術前の動画像に対しても、手術後の血流改善領域と同一の画素領域を枠で囲って識別表示しているが、次の図８のように、手術後の動画像においてのみ、血流改善領域を識別表示してもよい。

図８は、パラメトリック時相画像の動画表示において、血流改善領域を点滅させることで識別表示する一例を示す模式図である。図８の例では、表示部２４ｆは、血流改善領域を点滅表示させることで、関心領域３００として識別表示する。点滅表示(BLINKING DISPLAY)に際しては、例えば、関心領域３００の画素の明度を断続的に上げ下げすればよい。或いは、時相数が十分多い場合には、図８の例のように、偶数の時相でのみ、関心領域３００の明度を下げて表示すればよい。

図９は、パラメトリック時相画像の動画表示において、特定処理部２４ｅにより抽出される血管壁の強調表示の一例を示す模式図である。図９に示すように、表示部２４ｆは、手術後の各パラメトリック時相画像において、枠で囲う等によって関心領域（血流改善領域）３００を識別表示しつつ、特定処理部２４ｅにより抽出された血管壁を関心領域３００内で例えば太線により強調表示することができる。

図１０は、パラメトリック時相画像の動画表示において、血流改善領域のみをＤＳＡ画像によりグレースケールの動画で表示しつつ、残りの領域をパラメトリック時相画像によりカラー表示する一例を示す模式図である。上述のように、複数のＤＳＡ画像の各時相に対応して複数のパラメトリック時相画像が生成される場合、表示部２４ｆは、パラメトリック時相画像の動画表示において、関心領域（血流改善領域）３００のみを同時相のＤＳＡ画像の同一領域に置換して表示できる。

より詳細には、パラメトリック時相画像では各画素が例えば赤、青、緑の３原色の画素値をそれぞれ有するが、ＤＳＡ画像の各画素は一の画素値のみを有する。ＤＳＡ画像の各画素値は、造影剤投与後のある時相の画素値と、マスク画像の画素値とのＸ線透過率の差分に相当するからである。従って、ＤＳＡ画像の画素値により表示される関心領域３００は、有彩色のデータがないので、グレースケール表示となる。但し、残りの領域がカラー表示されるため、グレースケールで表示される関心領域３００のみが識別的表示態様となる。

なお、図７〜図１０では便宜上、グレースケールの模式図で示したが、各パラメトリック時相画像はカラーの画像データであるため、実際には図１０の関心領域３００を除いてカラー表示される。
また、図７〜図１０では血流改善領域を識別表示する例を述べたが、血管塞栓手術後における血流閉塞領域を識別表示する場合も、図７〜図１０と同様の表示態様を適用することができる。

＜本実施形態の動作説明＞
図１１は、血管拡張手術の前後において、本実施形態のＸ線診断装置１０により造影剤投与を伴った撮影及び血流改善領域の表示が行われる場合の動作の一例を示すフローチャートである。以下、前述の各図を適宜参照しながら、図１１に示すフローチャートのステップ番号に従って、Ｘ線診断装置１０の動作を説明する。

［ステップＳ１］血管拡張手術前において、被検体Ｐの同一領域に対して、公知の動作により時系列的な複数のＸ線画像の投影データが生成される。具体的には、システム制御部２６（図１参照）は、入力部２２を介して入力された撮影領域、管電流、管電圧、Ｘ線パルス幅等の一部の撮影条件に基づいて、造影剤投与の前後における各時相での全撮影条件を設定する。高電圧発生器４４は、システム制御部２６の制御に従って高電圧をＸ線管４２に供給し、Ｘ線管４２はＸ線を発生し、絞り装置４０によって被検体Ｐに対するＸ線照射範囲が制御される。

Ｘ線検出器３４は、被検体Ｐを透過したＸ線を検出して電気信号に変換し、これを投影データ生成部３０に入力する。投影データ生成部３０は、入力された電気信号からＸ線画像の投影データを生成して、投影データ記憶部２８に記憶させる。このようにして、造影剤投与の前後において、各画素の輝度がＸ線検出器３４のＸ線検出素子（図示せず）毎の受線量に応じた輝度となるように、被検体Ｐの同一領域に対して多時相の投影データが生成される。

なお、造影剤投与後のＸ線画像は複数の時相に対応するように複数必要であるが、造影剤投与前のＸ線画像は１つのみでもよい。
さらに、投影データ生成部３０は、投影データに基づいて、被検体の同一領域に対する時系列的な複数のＤＳＡ画像の画像データを生成し、これらＤＳＡ画像の画像データを投影データ記憶部２８に記憶させる。この後、ステップＳ２に進む。

［ステップＳ２］画像解析装置２４のＤＳＡ画像取得部２４ａは、上記の時系列的な複数のＤＳＡ画像の画像データを投影データ記憶部２８から取得する。パラメータ値取得部２４ｂは、時系列的な複数のＤＳＡ画像の画像データから、パラメータ値を同位置の画素毎に算出する（図２参照）。なお、ＴＴＰ，ＰＨ，ＡＵＣ等の内のどれをパラメータ値として用いるかについては、ステップＳ１で入力部２２を介して選択される構成としてもよいし、パラメータ値取得部２４ｂが自動決定する構成でもよい。

次に、カラー割当て部２４ｃは、前述の手順により、時系列的な複数のＤＳＡ画像に対して一のパラメトリック画像の画像データを生成する（図３参照）。

次に、時相画像生成部２４ｄは、前述の手順により、時系列的な複数のＤＳＡ画像にそれぞれ対応する時相の複数のパラメトリック時相画像の画像データを生成する。この生成方法については、図４で説明した第１の方法でも、図５で説明した第２の方法でもよい。この後、ステップＳ３に進む。

［ステップＳ３］血管拡張手術後において、ステップＳ１と同様にして、ステップＳ１と同一の被検体Ｐの同一の撮影領域に対して時系列的なＸ線画像が生成される。これにより、造影剤投与前後において、被検体Ｐの同一領域に対して多時相の投影データが生成される。さらに、投影データ生成部３０は、前述同様に、被検体の同一領域に対する時系列的な複数のＤＳＡ画像の画像データを生成し、これら画像データを投影データ記憶部２８に記憶させる。この後、ステップＳ４に進む。

［ステップＳ４］ステップＳ３と同様にして、画像解析装置２４により、手術後における、時系列的な複数のＤＳＡ画像にそれぞれ対応する時相の複数のパラメトリック時相画像の画像データが生成される。この後、ステップＳ５に進む。

［ステップＳ５］画像解析装置２４の特定処理部２４ｆは、手術後のパラメトリック画像の各画素値から、手術前パラメトリック画像の各画素値を差し引くことで差分画像を生成し、差分画像において、画素値が閾値以上の画素領域を血流改善領域として抽出する（図６参照）。また、特定処理部２４ｅは、前述のように、血流改善領域内の血管壁が強調された血管壁強調画像データを生成する。この後、ステップＳ６に進む。

［ステップＳ６］画像解析装置２４の表示部２４ｆは、ステップＳ５で抽出された血流改善領域を識別態様にしつつ、手術後のパラメトリック時相画像を時系列的に表示する。即ち、動画表示される。このとき、図７〜図１０のように、対比のために、手術前のパラメトリック時相画像を画面片側に、手術後のパラメトリック時相画像を画面の反対側に表示してもよい。

また、関心領域（血流改善領域）３００を識別表示については、図７のように長方形の枠で囲ってもよいし、図８のように点滅表示させてもよいし、図９のように血管壁を関心領域３００内で強調表示してもよい。或いは、図１０のように、血流改善領域のみをＤＳＡ画像によりグレースケールの動画で表示しつつ、残りの領域をパラメトリック時相画像によりカラー表示してもよい。以上が図１１のフローの説明である。

図１２は、血管塞栓手術の前後において、本実施形態のＸ線診断装置１０により造影剤投与を伴った撮影及び血流閉塞領域の表示が行われる場合の動作の一例を示すフローチャートである。前述の各図を適宜参照しながら、図１２に示すフローチャートのステップ番号に従って、Ｘ線診断装置１０の動作を説明する。

［ステップＳ１１］血管塞栓手術前において、図１１のステップＳ１と同様に、被検体Ｐの同一の撮影領域に対して造影剤投与前後の時系列的な複数のＸ線画像の投影データ及びＤＳＡ画像の画像データが生成される。この後、ステップＳ１２に進む。

［ステップＳ１２］図１１のステップＳ２と同様にして、画像解析装置２４は、手術前の各時相のパラメトリック時相画像の画像データを生成する。この後、ステップＳ１３に進む。

［ステップＳ１３］血管塞栓手術後において、前述同様に、ステップＳ１１と同一の被検体Ｐの同一の撮影領域に対して造影剤投与前後の時系列的な複数のＸ線画像の投影データ及びＤＳＡ画像の画像データが生成される。この後、ステップＳ１４に進む。

［ステップＳ１４］前述同様に、画像解析装置２４により、手術後における各時相のパラメトリック時相画像の画像データが生成される。この後、ステップＳ１５に進む。

［ステップＳ１５］特定処理部２４ｆは、手術前のパラメトリック画像の各画素値から、手術後のパラメトリック画像の各画素値を差し引くことで差分画像を生成し、差分画像において、画素値が閾値以上の画素領域を血流閉塞領域として抽出する。この後、ステップＳ１６に進む。なお、図１２のフローでは、手術の効果は血流閉塞なので、抽出された関心領域（血流閉塞領域）に対する血管壁強調の処理は実行されない。

［ステップＳ１６］表示部２４ｆは、前述の図７、図８、図１０と同様にして、ステップＳ１５で抽出された血流閉塞領域を識別態様にしつつ、手術後のパラメトリック時相画像を時系列的に表示する。
以上が本実施形態のＸ線診断装置１０の動作説明である。

＜本実施形態の効果＞
従来技術では、手術前又は手術後にパラメトリック画像が１つ生成され、パラメトリック画像は静止画表示されていた。静止画表示では、手術の効果である血流改善領域又は血流改善領域の目視による判別は、必ずしも容易ではなかった。

そこで本実施形態では、１のパラメトリック画像から、画素値の時間変化に応じて、複数の時相のＤＳＡ画像にそれぞれ対応する複数の時相のパラメトリック時相画像がそれぞれ生成される。そして、手術前、手術後のパラメトリック時相画像がそれぞれ動画表示されるので、血流量変化領域を従来よりも目視で判別し易くなる。

さらに本実施形態では、手術前のパラメトリック画像と、手術後のパラメトリック画像との差分画像に基づいて血流改善領域又は血流閉塞領域が自動抽出される。この抽出された血流改善領域又は血流閉塞領域は、パラメトリック時相画像の動画表示において、例えば図７〜図１０のような手法で識別表示される。従って、読影者は、手術の効果である血流改善領域又は血流改善領域を目視で容易に判別できる。

以上の画期的且つ新規な構成により、手術前と手術後との差異が僅かであっても、上記識別表示の効果により、読影者の見落としが少なくなる。従って、パラメトリック画像の静止画表示や、ＤＳＡ画像の動画表示、ＤＳＡ画像の各時相の静止画表示等を切り替えながら、僅かな差異を発見するといった従来の作業負担はなくなる。この結果、ユーザの利便性は大いに向上する。

＜本実施形態の補足事項＞
［１］Ｘ線診断装置１０内の一構成要素として画像解析装置２４が配置される例を述べた。本発明は、かかる実施形態に限定されるものではない。画像解析装置２４は、Ｘ線ＣＴ装置(X-Ray Computed Tomography Apparatus)や磁気共鳴イメージング装置などの他のモダリティ内に配置されてもよい。

或いは、画像解析装置２４が画像保管通信システム(PACS: Picture Archiving and Communication Systems)に接続された画像解析サーバとして機能する態様でもよい。上記画像保管通信システムは、Ｘ線診断装置などのモダリティにより得られた画像データを画像保管サーバに保存し、参照端末において必要な画像を呼び出して表示するシステムである。この場合、画像解析装置２４は、画像保管通信システムを介して、同一の被検体Ｐに対する時系列的な複数のＤＳＡ画像を取得し、上記実施形態と同様に血流量変化領域を識別態様にしつつ、パラメトリック時相画像を動画表示してもよい。

［２］上記実施形態では、取得される全てのＤＳＡ画像の各時相にそれぞれ対応するパラメトリック時相画像を生成し、生成した全てのパラメトリック時相画像を時系列順に表示する例を述べた。ＤＳＡ画像の全時相にそれぞれ対応するようにパラメトリック時相画像を生成する必要はない。

取得されるＤＳＡ画像の数が多い場合、即ち、ＤＳＡ画像の時相数が多い場合、動画として自然に見える程度に、パラメトリック時相画像の時相数を間引いてもよい。例えば１００時相のＤＳＡ画像を取得した場合、偶数の時相のＤＳＡ画像のみを用いて５０のパラメトリック時相画像を生成し、時系列順に表示してもよい。

［３］上記実施形態では、造影剤投与前後の撮影で得られる時系列的な２次元の画像データに基づいて時系列的な複数のＤＳＡ画像が生成され、これらＤＳＡ画像にそれぞれ対応する複数のパラメトリック画像が生成されて２次元的に動画表示される例を述べた。本発明の実施形態は、かかる態様に限定されるものではない。

例えば、Ｘ線ＣＴ装置での撮影や磁気共鳴イメージング装置での撮像のように、ボリュームデータが時系列的に得られる場合、時系列的な複数のパラメトリック画像（の視差画像）を裸眼式３次元表示装置で動画表示してもよい。この場合、表示部２４ｆは裸眼式３次元表示装置として構成され、各撮影時刻に対応する一のパラメトリック時相画像に対して複数の視差画像の画像データが生成される。

具体的には、ボリュームデータの例えば中央の断面の各画素値に基づいて上記のように時系列的なパラメトリック時相画像１〜５を生成する（前述のように、フレーム数は５に限定されるものではない）。そして、ボリュームデータの全断面の画素値に基づいて奥行き情報を算出し、奥行き情報に基づいて、複数の時相のパラメトリック時相画像のそれぞれに対して複数の視差画像を生成する。

各時相（時刻）に対して複数の視差画像が生成されれば、裸眼式３次元表示装置により、パラメトリック時相画像を立体的に動画表示することができる。視差画像の生成方法や、裸眼式の３次元画像表示装置については、例えば特開２００７−９４０２２号公報などに記載の従来技術を用いればよい。

［４］上記実施形態では、手術によって血流が改善又は閉塞した領域が自動抽出され、パラメトリック時相画像の動画表示において血流改善領域又は血流閉塞領域が識別的態様で表示される例を述べた。本発明の実施形態は、かかる態様に限定されるものではない。血流改善領域又は血流閉塞領域の抽出処理は必須ではなく、単にパラメトリック時相画像を動画表示してもよい。

動画表示に際しては、例えば、各時相のＤＳＡ画像と共に表示してもよい。具体的には例えば、表示部２４ｆの画面の右側ではＤＳＡ画像を時系列順に表示し、画面左側では、画面右側に表示されているＤＳＡ画像に対応する時相のパラメトリック時相画像を追従して表示してもよい。画面左側におけるパラメトリック時相画像の動画表示に際しては、手術前の画像との差分によって上記実施形態のように抽出される血流改善領域又は血流閉塞領域を識別表示することが好ましい。

或いは、手術後のＤＳＡ画像、パラメトリック画像、パラメトリック時相画像を共に表示してもよい。例えば、画面上段において上記同様に手術後の各ＤＳＡ画像を時系列順に表示し、画面中段において上記同様にパラメトリック時相画像を時系列順に表示し、画面下段において手術後の一のパラメトリック画像を静止画表示してもよい。

［５］請求項の用語と実施形態との対応関係を説明する。なお、以下に示す対応関係は、参考のために示した一解釈であり、本発明を限定するものではない。
カラー割当て部２４ｃにより生成されるパラメトリック画像の画像データは、請求項記載のカラーマップの一例である。

時相画像生成部２４ｄにより生成されるパラメトリック時相画像の画像データは、請求項記載の時相画像のカラー画像データの一例である。

特定処理部２４ｅにより生成される手術前のパラメトリック画像と、手術後のパラメトリック画像との差分画像は、請求項記載の差分マップの一例である。
特定処理部２４ｅにより抽出される血流改善領域及び血流閉塞領域は、請求項記載の血流量変化領域の一例である。

被検体Ｐに対するＸ線撮影を実行することでＸ線画像の投影データ及びＤＳＡ画像の画像データを生成するシステム制御部２６、投影データ記憶部２８、投影データ生成部３０、Ｘ線検出器３４、Ｃアーム３６、天板３８、絞り装置４０、Ｘ線管４２、高電圧発生器４４、絞り制御機構４６、天板移動機構４８、Ｃアーム動作機構５０、検出器移動機構５４は、請求項記載のＸ線撮影部の一例である。

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

１０：Ｘ線診断装置，２４：画像解析装置，
２４ａ：ＤＳＡ画像取得部，２４ｂ：パラメータ値取得部，２４ｃ：カラー割当て部，
２４ｄ：時相画像生成部，２４ｅ：特定処理部，２４ｆ：表示部

Claims (14)

1. 同一の被検体に対する時系列的な複数のＤＳＡ画像の画像データを取得するＤＳＡ画像取得部と、
   前記複数のＤＳＡ画像の画像データにおける、前記被検体の同一領域に対応する画素毎の画素値の時間変化に基づいて、画素毎にパラメータ値を取得するパラメータ値取得部と、
   前記パラメータ値に応じた色が前記被検体の同一領域に対応する画素毎に割り当てられように、且つ、複数の画素に有彩色が割り当てられるように、一のカラー画像の画像データを生成するカラー割当て部と、
   前記一のカラー画像の画像データの各画素の色を変換することで、互いに異なる複数の時相にそれぞれ対応する複数の時相画像のカラー画像データを生成する時相画像生成部と
   を備えることを特徴とする画像解析装置。
2. 請求項１記載の画像解析装置において、
   前記複数の時相画像間の差分画像に基づいて、血流量変化領域を特定する特定処理部をさらに備える
   ことを特徴とする画像解析装置。
3. 請求項２記載の画像解析装置において、
   前記ＤＳＡ画像取得部は、造影剤投与の前後にそれぞれ撮影された前記時系列的な複数のＤＳＡ画像の画像データを取得し、
   前記特定処理部は、造影剤投与前に撮影された前記複数のＤＳＡ画像に基づいて生成される前記一のカラー画像の画像データと、造影剤投与後に撮影された前記複数のＤＳＡ画像に基づいて生成部により生成される前記一のカラー画像の画像データとの間の差分マップを算出し、前記差分マップに各画素に対して閾値処理を施すことで前記血流量変化領域を特定する
   ことを特徴とする画像解析装置。
4. 請求項２又は請求項３に記載の画像解析装置において、
   前記時相画像生成部は、各画素の明度、彩度、透過度の少なくとも１つが画像間で変わるように色変換を実行することで、前記複数のＤＳＡ画像に時系列的にそれぞれ対応するように前記複数の時相画像のカラー画像データを生成する
   ことを特徴とする画像解析装置。
5. 請求項２乃至請求項４のいずれか１項に記載の画像解析装置において、
   前記複数の時相画像をカラー動画として表示する表示部をさらに備える
   ことを特徴とする画像解析装置。
6. 請求項５記載の画像解析装置において、
   前記表示部は、前記複数の時相画像の表示に際して、前記血流量変化領域を枠で囲って識別表示する
   ことを特徴とする画像解析装置。
7. 請求項５記載の画像解析装置において、
   前記表示部は、前記複数の時相画像の表示に際して、前記血流量変化領域を点滅させることで識別表示する
   ことを特徴とする画像解析装置。
8. 請求項５記載の画像解析装置において、
   前記表示部は、前記複数の時相画像の表示に際して、前記血流量変化領域のみをグレースケールで表示すると共に残りの領域をカラー表示することで、前記血流量変化領域を識別表示する
   ことを特徴とする画像解析装置。
9. 請求項５記載の画像解析装置において、
   前記特定処理部は、前記血流量変化領域の血管壁が強調された血管壁強調画像データを生成し、
   前記表示部は、前記複数の時相画像の表示に際して、前記血管壁強調画像データに基づいて血管壁を識別表示する
   ことを特徴とする画像解析装置。
10. 造影剤投与の前後において被検体を透過したＸ線を検出することでＸ線画像の投影データを生成し、造影剤投与後の時系列的な複数のＸ線画像と、造影剤投与前のＸ線画像との各差分に基づいて前記被検体に対する時系列的な複数のＤＳＡ画像の画像データを生成するＸ線撮影部と、
    前記複数のＤＳＡ画像の画像データにおける、前記被検体の同一領域に対応する画素毎の画素値の時間変化に基づいて、画素毎にパラメータ値を取得し、前記パラメータ値に応じた色が画素毎に割り当てられように一のカラー画像の画像データを生成し、前記一のカラー画像の画像データの各画素の色を変換することで、互いに異なる複数の時相にそれぞれ対応する複数の時相画像の各カラー画像データを生成する請求項１乃至請求項９のいずれか１項に記載の画像解析装置と
    を備えることを特徴とするＸ線診断装置。
11. 前記時相画像生成部は、前記一のカラー画像の画像データの各画素の色を、前記複数の時相にそれぞれ対応する前記複数のＤＳＡ画像の画像データを用いて変換することにより、前記複数の時相画像のカラー画像データを生成する、
    請求項１乃至請求項９のいずれか１項に記載の画像解析装置。
12. 前記時相画像生成部は、前記一のカラー画像の画像データの各画素の色を、前記複数の時相にそれぞれ対応する複数の前記パラメータ値を用いて変換することにより、前記複数の時相画像のカラー画像データを生成する、
    請求項１乃至請求項９のいずれか１項に記載の画像解析装置。
13. 前記時相画像生成部は、各画素の明度、彩度、透過度の少なくとも１つが画像間で変わるように色変換を実行することで、前記複数のＤＳＡ画像に時系列的にそれぞれ対応するように前記複数の時相画像のカラー画像データを生成する
    請求項１に記載の画像解析装置。
14. 前記複数の時相画像をカラー動画として表示する表示部をさらに備える
    請求項１に記載の画像解析装置。