(57)【要約】      （修正有）【課題】術前の画像と術後の画像とを正確に比較することを可能にするＸ線診断装置及び画像処理装置を提供する。【解決手段】Ｘ線診断装置１は、プロファイル計測部２５と、補正ファクター同定部２６と、ＤＳＡ画像補正部２７と、制御部３４とを備える。プロファイル計測部は、被検体の頭部に対して略同一方向から撮影された撮影時期の異なる２枚の差分画像における血管を含む略同一位置にそれぞれ設定された関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する。補正ファクター同定部は、プロファイル計測部によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する。ＤＳＡ画像補正部は、同定された補正ファクターに基づいて、２枚の差分画像うちの少なくとも一方を補正する。制御部は、ＤＳＡ画像補正部によって少なくとも一方が補正された２枚の差分画像に基づいた情報を表示部にて表示するように制御する。【選択図】図１
2014193320.tif 000002
【特許請求の範囲】
【請求項１】
被検体の頭部に対して略同一方向から撮影された撮影時期の異なる２枚の差分画像における血管を含む略同一位置にそれぞれ設定された関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する計測部と、
前記計測部によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する同定部と、
前記同定部によって同定された補正ファクターに基づいて、前記２枚の差分画像のうちの少なくとも一方を補正する補正部と、
前記補正部によって少なくとも一方が補正された２枚の差分画像に基づいた情報を所定の表示部にて表示するように制御する制御部と、
を備える、Ｘ線診断装置。
【請求項２】
前記関心領域は、前記血管である主幹動脈または治療に影響が無い毛細血管を含む領域である、請求項１記載のＸ線診断装置。
【請求項３】
前記２枚の差分画像のうち少なくとも一方の画像における前記主幹動脈の位置を算出する算出部をさらに有し、
前記関心領域は、前記算出部で算出された前記主幹動脈を含む領域である、請求項２記載のＸ線診断装置。
【請求項４】
前記表示部は、前記少なくとも一方が補正された２枚の差分画像を差分した結果に基づくパーフュージョン画像を表示する、請求項１記載のＸ線診断装置。
【請求項５】
前記同定部は、前記２枚の差分画像における２つのプロファイルが略一致するように、造影剤到達までの時間のずれ、前記被検体の拍出量に関するゲイン及び血流速度のうち、少なくとも一つを同定する、請求項１記載のＸ線診断装置。
【請求項６】
前記２枚の差分画像は、撮影時期が治療前及び治療後に得られた画像である、請求項１記載のＸ線診断装置。
【請求項７】
前記差分画像の下端から所定の範囲内に含まれる血管のうち、血管径が最大となる血管を抽出し、抽出した血管において、所定の領域ごとに濃度値の部分積分値のピーク値を算出し、算出したピーク値が最低となる領域を前記関心領域として設定する設定部をさらに備える、請求項１記載のＸ線診断装置。
【請求項８】
前記設定部は、前記差分画像の生成に用いられたマスク画像の濃度値に基づいて、前記所定の領域ごとのＸ線透過度を算出し、算出したＸ線透過度が低い領域を前記関心領域の対象外とする、請求項７記載のＸ線診断装置。
【請求項９】
前記制御部は、補正後の差分画像に基づいて算出された部位ごとのプロファイルを元に、血流量差及び差異度を算出する、請求項１記載のＸ線診断装置。
【請求項１０】
前記制御部は、Ｘ線撮影条件、撮影角度、ＦＯＶ、ＳＩＤ、コリメータの位置、補償フィルターの位置、Ｘ線焦点サイズ、線質調整フィルター、造影剤注入までの時間及び造影条件のうち、少なくとも１つを一致させて前記撮影時期の異なる２枚の差分画像を撮影するように制御する、請求項１記載のＸ線診断装置。
【請求項１１】
被検体の頭部に対して略同一方向から造影剤を注入しながら撮影された撮影時期の異なる２枚のＸ線動画像を元に、造影剤の影響のある複数の画像と造影剤の影響の殆どない画像とを引き算し差分画像を算出する差分画像算出部と、
前記２枚のＸ線動画像の中の造影剤の影響が殆どない画像を元に位置ずれを同定する位置ずれ同定部と、
前記同定された位置ずれ情報を元に２枚のＸ線差分画像の内、少なくとも一方を位置ずれ補正する位置ずれ補正部と、
前記少なくとも一方が補正された２枚のＸ線差分画像に含まれる血管を含む略同一位置にそれぞれ設定された関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する計測部と、
前記計測部によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する同定部と、
前記同定部によって同定された補正ファクターに基づいて、前記２枚の差分画像のうちの少なくとも一方を補正する補正部と、
前記補正部によって少なくとも一方が補正された２枚の差分画像に基づいた情報を所定の表示部にて表示するように制御する制御部と、
を備える、Ｘ線診断装置。
【請求項１２】
被検体の頭部に対して略同一方向から撮影された撮影時期の異なる２枚の差分画像における血管を含む略同一位置にそれぞれ設定された関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する計測部と、
前記計測部によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する同定部と、
前記同定部によって同定された補正ファクターに基づいて、前記２枚の差分画像のうちの少なくとも一方を補正する補正部と、
前記補正部によって少なくとも一方が補正された２枚の差分画像に基づいた情報を所定の表示部にて表示するように制御する制御部と、
を備える、画像処理装置。
【請求項１３】
被検体の頭部に対して略同一方向から造影剤を注入しながら撮影された撮影時期の異なる２枚のＸ線動画像を元に、造影剤の影響のある複数の画像と造影剤の影響の殆どない画像とを引き算し差分画像を算出する差分画像算出部と、
前記２枚のＸ線動画像の中の造影剤の影響が殆どない画像を元に位置ずれを同定する位置ずれ同定部と、
前記同定された位置ずれ情報を元に２枚のＸ線差分画像の内、少なくとも一方を位置ずれ補正する位置ずれ補正部と、
前記少なくとも一方が補正された２枚のＸ線差分画像に含まれる血管を含む略同一位置にそれぞれ設定された関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する計測部と、
前記計測部によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する同定部と、
前記同定部によって同定された補正ファクターに基づいて、前記２枚の差分画像のうちの少なくとも一方を補正する補正部と、
前記補正部によって少なくとも一方が補正された２枚の差分画像に基づいた情報を所定の表示部にて表示するように制御する制御部と、
を備える、画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施の形態は、Ｘ線診断装置及び画像処理装置に関する。
【背景技術】
【０００２】
従来、Ｘ線ＣＴ（Computed  Tomography）装置による脳梗塞の診断においては、ブレインパーフュージョン（Brain  Perfusion）解析が知られている。ブレインパーフュージョン解析においては、造影剤を注入して撮影した脳の断層像から血液の還流状態を表す画像であるブレインパーフュージョン画像が脳梗塞の診断に用いられる。
【０００３】
また、近年、Ｘ線ＣＴ装置によって撮像された頭部断面画像に対して、頭部の左右半球を２つに分割する境界線を手動又は自動で設定し、境界線で分割された一方の画像を反転させたうえで他方の画像に重ね、重ね合わせた画像の差分画像を作成する解析方法が知られている。この解析方法によれば、差分画像において病変領域が浮かび上がるので左右の比較読影が容易になり、脳梗塞などの血流異常をチェックすることが可能となる。しかしながら、上述した従来技術においては、術前の画像と術後の画像との正確な比較ができない場合があった。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−１５３８７０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明が解決しようとする課題は、術前の画像と術後の画像とを正確に比較することを可能にするＸ線診断装置及び画像処理装置を提供することである。
【課題を解決するための手段】
【０００６】
実施の形態のＸ線診断装置は、計測部と、同定部と、補正部と、表示制御部とを備える。計測部は、被検体の頭部に対して略同一方向から撮影された撮影時期の異なる２枚の差分画像における血管を含む略同一位置にそれぞれ設定された関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する。同定部は、前記計測部によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する。補正部は、前記同定部によって同定された補正ファクターに基づいて、前記２枚の差分画像うちの少なくとも一方を補正する。表示制御部は、前記補正部によって少なくとも一方が補正された２枚の差分画像に基づいた情報を所定の表示部にて表示するように制御する。
【図面の簡単な説明】
【０００７】
【図１】図１は、第１の実施形態に係るＸ線診断装置の構成の一例を示す図である。
【図２】図２は、第１の実施形態に係るＲＯＩ設定部の構成の一例を示す図である。
【図３】図３は、第１の実施形態に係るＸ線診断装置による処理の手順を示すフローチャートである。
【図４】図４は、第１の実施形態に係るＲＯＩ設定部によって設定されたＲＯＩの一例を示す図である。
【図５】図５は、第１の実施形態に係るＲＯＩ設定部による処理の手順を示すフローチャートである。
【図６】図６は、第１の実施形態に係る血管探索部による処理の一例を説明するための図である。
【図７Ａ】図７Ａは、第１の実施形態に係る最低ピーク領域同定部による部分領域ごとの積分値の算出処理を説明するための図である。
【図７Ｂ】図７Ｂは、第１の実施形態に係る最低ピーク領域同定部による部分領域ごとの積分値の算出処理を説明するための図である。
【図７Ｃ】図７Ｃは、第１の実施形態に係る最低ピーク領域同定部による部分領域ごとの積分値の算出処理を説明するための図である。
【図８】図８は、第１の実施形態に係る最低ピーク領域同定部による最低ピークの抽出処理を説明するための図である。
【図９】図９は、第１の実施形態に係るプロファイル計測部による処理の一例を説明するための図である。
【図１０】図１０は、第１の実施形態に係る表示部にて表示される血流チェック画像の一例を示す図である。
【図１１Ａ】図１１Ａは、第１の実施形態に係る表示部にて表示される解析結果の一例を示す図である。
【図１１Ｂ】図１１Ｂは、第１の実施形態に係る表示部にて表示される解析結果の一例を示す図である。
【図１１Ｃ】図１１Ｃは、第１の実施形態に係る表示部にて表示される解析結果の一例を示す図である。
【図１２】図１２は、第１の実施形態に係る表示部にて表示される表示情報の一例を示す図である。
【発明を実施するための形態】
【０００８】
以下に、本発明に係るＸ線診断装置及び画像処理装置の実施形態を図面に基づいて詳細に説明する。なお、以下に示す実施形態によって本発明が限定されるものではない。
【０００９】
（第１の実施形態）
図１は、第１の実施形態に係るＸ線診断装置１の構成の一例を示す図である。図１に示すように、第１の実施形態に係るＸ線診断装置１は、Ｘ線撮影機構１０と、画像処理装置１００とを有する。Ｘ線撮影機構１０は、Ｘ線管球１１と、検出器（ＦＰＤ（Flat  Panel  Detector））１２と、Ｃ型アーム１３と、寝台１４とを有する。Ｃ型アーム１３は、Ｘ線管球１１及び検出器１２を支持し、土台（図示を省略）に設けられたモータでプロペラのように被検体Ｐの周りを高速回転する。  
【００１０】
画像処理装置１００は、図１に示すように、Ａ／Ｄ（Analog／Digital）変換部２１と、位置ずれ同定部２２と、位置ずれ補正部２３と、ＲＯＩ（Region  of  Interest）設定部２４と、プロファイル計測部２５と、補正ファクター同定部２６と、ＤＳＡ（Digital  Subtraction  Angiography）画像補正部２７と、第１サブトラクション部２８と、第２サブトラクション部２９と、フィルタリング部３０と、アフィン変換部３１と、ＬＵＴ（Look  Up  Table）３２と、画像メモリ３３と、制御部３４と、表示部４０とを有する。また、画像処理装置１００は、図示していないが、例えば、マウスやキーボード、トラックボール、ポインティングデバイスなど、Ｘ線診断装置１に対する各種操作を操作者から受け付ける入力部を有する。
【００１１】
表示部４０は、画像処理装置１００によって処理された各種画像や、ＧＵＩ（Graphical  User  Interface）などの各種情報を表示する。例えば、表示部４０は、ＣＲＴ（Cathode  Ray  Tube）モニタや液晶モニタなどである。Ａ／Ｄ変換部２１は、検出器１２に接続され、検出器１２から入力されたアナログ信号をディジタル信号に変換し、変換したディジタル信号をＸ線収集画像として画像メモリ３３に格納する。画像メモリ３３は、Ｘ線収集画像を記憶する。
【００１２】
位置ずれ同定部２２は、異なるタイミングで撮影した２枚のＤＳＡ画像の位置ずれを同定する。位置ずれ補正部２３は、位置ずれ同定部２２によって同定された位置ずれを元に、２枚のＤＳＡ画像の位置ずれを補正する。ＲＯＩ設定部２４は、ＤＳＡ画像上にＲＯＩを設定する。具体的には、ＲＯＩ設定部２４は、ＤＳＡ画像の下端から所定の範囲内に含まれる血管のうち、血管径が最大となる血管を抽出し、抽出した血管において、所定の領域ごとに濃度値の部分積分値のピーク値を算出し、算出したピーク値が最低となる領域をＲＯＩとして設定する。
【００１３】
図２は、第１の実施形態に係るＲＯＩ設定部２４の構成の一例を示す図である。図２に示すように、ＲＯＩ設定部２４は、血管探索部２４ａと、最大血管同定部２４ｂと、最低ピーク領域同定部２４ｃとを有する。血管探索部２４ａは、ＤＳＡ画像の下端から一定範囲離れた横軸上でＤＳＡ画像を元に血管を探索する。最大血管同定部２４ｂは、血管探索部２４ａによって探索された血管から血管径の最も大きい血管を同定する。最低ピーク領域同定部２４ｃは、最大血管同定部２４ｂによって同定された最大血管を下方向にトレースし、濃度値の部分積分値のピークが最も低い部位を同定する。
【００１４】
図１に戻って、プロファイル計測部２５は、ＤＳＡ画像に設定されたＲＯＩ上でプロファイルを計測する。具体的には、プロファイル計測部２５は、被検体の頭部に対して略同一方向から撮影された時相（撮影時期）の異なる２枚のＤＳＡ画像における血管を含む略同一位置にそれぞれ設定されたＲＯＩにおける造影剤濃度に関するプロファイルをそれぞれ計測する。例えば、プロファイル計測部２５は、血管である主幹動脈または治療に影響が無い毛細血管を含む関心領域における造影剤濃度に関するプロファイルをそれぞれ計測する。領域補正ファクター同定部２６は、２つのプロファイルが略一致するように補正ファクターを同定する。具体的には、補正ファクター同定部２６は、プロファイル計測部２５によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する。例えば、補正ファクター同定部２６は、２枚の差分画像における２つのプロファイルが略一致するように、造影剤到達までの時間のずれ、被検体の拍出量に関するゲイン及び血流速度のうち、少なくとも一つを同定する。
【００１５】
ＤＳＡ画像補正部２７は、補正ファクター同定部２６によって同定された補正ファクターを元に２枚のＤＳＡ画像のうちの少なくとも一方のＤＳＡ画像を補正する。第１サブトラクション部２８は、造影剤投与前の画像と造影剤投与後の画像とをサブトラクションする。第２サブトラクション部２９は、ＤＳＡ画像補正部２７によって補正されたＤＳＡ画像と、もう一方のＤＳＡ画像とをサブトラクションする。フィルタリング部３０は、高周波強調フィルタリングなどを行う。アフィン変換部３１は、画像の拡大や縮小、移動などを行う。ＬＵＴ３２は、諧調変換を行う。
【００１６】
制御部３４は、Ｘ線診断装置１全体を制御する。具体的には、制御部３４は、Ｘ線収集画像の収集、表示画像の生成、表示部４０における表示画像の表示などに係る各種処理を制御する。例えば、制御部３４は、ＤＳＡ画像補正部２７によって少なくとも一方が補正された２枚のＤＳＡ画像に基づいた情報（例えば、差分情報など）を表示部４０にて表示するように制御する。
【００１７】
第１の実施形態に係るＸ線診断装置１は、上述した構成により、従来技術では困難であった術前の画像と術後の画像との正確な比較を可能にする。
【００１８】
例えば、Ｘ線診断装置を用いた頭部の治療には、カテーテルを狭窄部まで挿入し、カテーテルの周りに設置されているバルーンを膨らませて狭窄部を広げる治療がある。かかる治療は、インターベンションと呼ばれる。インターベンションでは、バルーンを広げた際に狭窄部の一部の小片が抹消に流れ、その小片によって脳の毛細血管に梗塞が生じてしまうこともある。
【００１９】
そこで、インターベンションの術後、ブレインパーフュージョン画像により梗塞が生じているか否かの診断が行われる。また、従来技術においては、例えば、左右の脳の血流状態を比較読影することにより、血流異常のチェックが行われる。しかしながら、上述した従来技術では、術前の状態が反映されないため、治療前に既に発生していた梗塞や、多少の濃染の傾向などを誤認してしまう場合がある。例えば、慢性的に梗塞が生じており、その梗塞に対する治療の必要がない患者であっても、既に発生していた梗塞を治療によって新たに発生した梗塞として診断されてしまい、無駄な血栓溶解療法が施されてしまう。
【００２０】
そこで、本実施形態に係るＸ線診断装置１は、以下、詳細に説明する処理により、インターベンションや、血栓溶解療法などの手技の術前の画像と、術後の画像とを正確に比較することを可能にする。図３は、第１の実施形態に係るＸ線診断装置１による処理の手順を示すフローチャートである。
【００２１】
図３に示すように、Ｘ線診断装置１においては、まず、インターベンションや、血栓溶解療法などの手技開始前に術前のＤＳＡ画像を収集する（ステップＳ１０１）。具体的には、Ｘ線診断装置１においては、Ｃアーム１３を任意の方向にセットし、造影前に数フレーム分のマスク画像を撮影し、その後、血管に造影剤が流れている間に連続的にコントラスト画像を撮影する。造影前の数フレーム分のマスク画像と複数のコントラスト画像は、Ａ／Ｄ変換部２１によってディジタル信号に変換され、画像メモリ３３に記憶される。
【００２２】
第１サブトラクション部２８は、画像メモリ３３によって記憶された数フレーム分のマスク画像を読み出し、読み出したマスク画像を加算平均してノイズの少ない平均マスク画像を生成する。そして、第１サブトラクション部２８は、複数のコントラスト画像から平均マスク画像をそれぞれサブトラクション（Ｌｏｇサブ）することで、術前のＤＳＡ画像を生成する。
【００２３】
第１サブトラクション部２８によって術前のＤＳＡ画像が生成されると、制御部３４は、生成された術前のＤＳＡ画像を表示部４０にて動画表示させる。なお、生成された術前のＤＳＡ画像は、画像メモリ３３によって記憶される。
【００２４】
その後、例えば、インターベンションや、血栓溶解療法などの手技が行われ、治療が完了すると、Ｘ線診断装置１は、術前と同一のデータ収集プログラムを指定することで術前と同様のデータ収集を行って、術後のＤＳＡ画像を収集する（ステップＳ１０２）。
【００２５】
そして、制御部３４は、生成された術後のＤＳＡ画像を表示部４０にて動画表示させる。なお、生成された術後のＤＳＡ画像は、画像メモリ３３によって記憶される。
【００２６】
そして、図示しないＧＵＩを介してユーザから術前のＤＳＡ画像と術後のＤＳＡ画像による血流比較の指示を受付けると、位置ずれ同定部２２は、画像メモリ３３から術前のＤＳＡ画像及び術後のＤＳＡ画像それぞれの平均マスク画像を読み出し、位置ずれ量を同定する（ステップＳ１０３）。例えば、位置ずれ同定部２２は、以下に示す式（１）により、術前のＤＳＡ画像のマスク画像（Ｍpre（ｉ，ｊ））と、術後のＤＳＡ画像のマスク画像（Ｍpost（ｉ，ｊ））との位置ずれ量を同定する。ここで、式（１）におけるＣＲ（Δｉ，Δｊ）は、（Ｍpre（ｉ，ｊ））と（Ｍpost（ｉ，ｊ））との位置ずれ量を示す。また、式（１）におけるＮは、画像サイズを示す。
【００２７】
【数１】
2014193320.tif 000003
【００２８】
位置ずれ同定部２２は、逐次近似アルゴリズム（algorithm）で（Δｉ，Δｊ）を変化させながら、ＣＲ（Δｉ，Δｊ）が最小になる（Δｉ，Δｊ）を探索する。そして、位置ずれ同定部２２は、探索した（Δｉ，Δｊ）を位置ずれ補正部２３に送る。なお、ここでは説明を簡単にするため、２次元の方向への位置ずれを検出する場合について説明するが、さらに回転方向の位置ずれを検出するのが望ましい。すなわち、位置ずれ同定部２２は、（Δｉ，Δｊ，Δθ）を同定することが望ましい。
【００２９】
続いて、位置ずれ補正部２３は、位置ずれ同定部２２から位置ずれデータを受付けると、受付けた位置ずれデータ「（Δｉ，Δｊ）又は（Δｉ，Δｊ，Δθ）」を用いて、術前のＤＳＡ画像又は術後のＤＳＡ画像のうちの一方を補正する（ステップＳ１０４）。
【００３０】
そして、ＲＯＩ設定部２４は、ＤＳＡ画像にＲＯＩを設定する（ステップＳ１０５）。具体的には、ＲＯＩ設定部２４は、位置ずれが補正された術前のＤＳＡ画像及び術後のＤＳＡ画像のうち、一方のＤＳＡ画像を用いて、術前のＤＳＡ画像及び術後のＤＳＡ画像の略同一位置にＲＯＩを設定する。例えば、ＲＯＩ設定部２４は、２枚の差分画像のうち少なくとも一方の画像における主幹動脈の位置を算出し、算出した主幹動脈を含む領域をＲＯＩとして設定する。図４は、第１の実施形態に係るＲＯＩ設定部２４によって設定されたＲＯＩの一例を示す図である。
【００３１】
図４に示すように、ＲＯＩ設定部２４は、術前ＤＳＡ画像及び術後ＤＳＡ画像の略同一位置にＲＯＩ５０及びＲＯＩ５１をそれぞれ設定する。ここで、以下、ＲＯＩ設定部２４による処理の詳細について説明する。図５は、第１の実施形態に係るＲＯＩ設定部２４による処理の手順を示すフローチャートである。なお、図５に示す処理は、図３のステップＳ１０５における処理に相当する。
【００３２】
図５に示すように、ＲＯＩ設定部２４においては、血管探索部２４ａは、位置ずれが補正された術前のＤＳＡ画像及び術後のＤＳＡ画像のうち、一方のＤＳＡ画像を用いて、ＤＳＡ画像の下端から所定の範囲内に含まれる血管を探索する（ステップＳ２０１）。図６は、第１の実施形態に係る血管探索部２４ａによる処理の一例を説明するための図である。ここで、図６の（Ａ）においては、画像サイズ「Ｎ×Ｎ」のＤＳＡ画像内の血管を探索する場合について示す。
【００３３】
例えば、血管探索部２４ａは、図６の（Ａ）に示すように、ＤＳＡ画像の下端からＮ／４の位置を血管探索の初期位置とする。そして、血管探索部２４ａは、図６の（Ｂ）に示すように、ＤＳＡ画像のＮ／４の位置の横軸上のＤＳＡ値のプロファイルを測定する。その後、血管探索部２４ａは、測定したプロファイルにおいて、ＤＳＡ値が所定の閾値「Ｔｈ」を越えたピークを血管として抽出する。例えば、血管探索部２４ａは、図６に示すように、「Ｔｈ」を越えたピークに対応する血管６０及び血管６１を抽出する。そして、血管探索部２４ａは、プロファイルと血管６０及び血管６１の位置情報を最大血管同定部２４ｂに送信する。
【００３４】
最大血管同定部２４ｂは、血管探索部２４ａによって抽出された血管において血管径が最大となる血管を同定する（ステップＳ２０２）。例えば、最大血管同定部２４ｂは、図６の（Ｂ）に示すプロファイルの結果から血管径（ピークの幅）を計測し、血管径が最大である血管６１を最大血管として同定する。そして、最大血管同定部２４ｂは、計測した血管径の情報及び血管６１の位置情報を最低ピーク領域同定部２４ｃに送信する。上述したように、血管探索部２４ａが、ＤＳＡ画像の下端からＮ／４の位置を血管探索の初期位置として血管を探索し、最大血管同定部２４ｂが、血管径が最大となる血管を同定することで、ＲＯＩを設定する血管として主幹動脈を抽出することができる。
【００３５】
最低ピーク領域同定部２４ｃは、最大血管同定部２４ｂによって同定された血管の部分領域ごとに、時間ごとの濃度値の積分値を算出する（ステップＳ２０３）。具体的には、最低ピーク領域同定部２４ｃは、ＤＳＡ画像の縦軸のＮ／４の位置から下に向かって血管中心をトレースしながら、部分領域の濃度値の積分値を時間ごとに算出する。図７Ａ〜図７Ｃは、第１の実施形態に係る最低ピーク領域同定部２４ｃによる部分領域ごとの積分値の算出処理を説明するための図である。なお、図７Ｂにおいては、図７Ａの領域５３を拡大した図１を示す。
【００３６】
例えば、最低ピーク領域同定部２４ｃは、図７Ａに示すように、ＤＳＡ画像の縦軸の下端からＮ／４の位置から下に向かって血管６１の血管中心をトレースしながら部分領域ごとの濃度値の積分値を算出する。ここで、最低ピーク領域同定部２４ｃは、例えば、図７Ｂに示すように、横「１．５Ｌ（＝１．５×血管径）」×縦「３ピクセル」の矩形を部分領域として設定し、各部分領域ごとに、積分値を算出する。なお、部分領域のサイズは、ユーザによって任意に設定することができる。
【００３７】
最低ピーク領域同定部２４ｃは、造影剤が流れている間のＤＳＡ画像それぞれについて、各部分領域ごとに濃度値の積分値を算出する。すなわち、最低ピーク領域同定部２４ｃは、部分領域ごとに、図７Ｃに示すような積分値の時間変化を算出する。
【００３８】
そして、最低ピーク領域同定部２４ｃは、算出した時間ごとの積分値のピーク値を部分領域ごとに抽出する（ステップＳ２０４）。例えば、最低ピーク領域同定部２４ｃは、図７Ｃに示すように、部分領域ごとに積分値がピークとなる代表値を抽出する。
【００３９】
その後、最低ピーク領域同定部２４ｃは、最低ピーク値を示す部分領域をＲＯＩに設定する（ステップＳ２０５）。具体的には、最低ピーク領域同定部２４ｃは、各部分領域の代表値を比較して、最低値を示す部分領域をＲＯＩとして設定する。図８は、第１の実施形態に係る最低ピーク領域同定部２４ｃによる最低ピークの抽出処理を説明するための図である。図８においては、横軸に図７Ａに示す「０」から「Ｎ／４」を示し、縦軸に積分値を示す。すなわち、図８においては、各部分領域の代表値をプロットしたグラフを示す。
【００４０】
例えば、最低ピーク領域同定部２４ｃは、図８に示すように、各部分領域の代表値をグラフ上にプロットして、最低ピークを同定する。そして、最低ピーク領域同定部２４ｃは、同定した最低ピークに対応する部分領域をＲＯＩとして設定する。
【００４１】
上述したように、最低ピーク領域同定部２４ｃは、各部分領域の代表値の中で最低値を示す部分領域をＲＯＩとして設定する。すなわち、最低ピーク領域同定部２４ｃは、積分値の時間変化が最も少ない部分領域をＲＯＩとして設定する。ＲＯＩを設定する主幹動脈は、血管径がほとんど変化しないが、実際の主幹動脈における濃度値は大きく変化する。これは、主に、血管の走行と、ビームハードニングに起因する。
【００４２】
例えば、血管の走行がＸ線の照射方向に対して９０度をなしている場合には、Ｘ線が曝射される造影剤の量は血管径の厚み分のみとなり、濃度値は最低となる。しかしながら、血管の走行がＸ線の照射方向と平行である場合には、Ｘ線が曝射される造影剤の量は血管が走行している奥行き分となり、濃度値が最大となる。
【００４３】
従って、血管の走行がＸ線の照射と平行となっている場合、血流を正しく評価できないため、血管の走行がＸ線の照射に対して９０度に近い角度をなしている領域をＲＯＩとして設定することが望ましい。そこで、最低ピーク領域同定部２４ｃは、濃度値の積分値が最低となる部分領域をＲＯＩとして設定することで、血管の走行がＸ線の照射に対して最も９０度に近い領域を抽出する。
【００４４】
また、例えば、厚い骨と重なった血管領域は、ビームハードニングによって濃度値が低下する。そこで、ＲＯＩ設定部２４が、ＤＳＡ画像の生成に用いられたマスク画像の濃度値に基づいて、部分領域ごとのＸ線透過度を算出し、算出したＸ線透過度が低い領域をＲＯＩの対象外とするようにしてもよい。
【００４５】
以上、ＲＯＩ設定部２４による処理について説明した。図３に戻って、上述したようにＤＳＡ画像にＲＯＩが設定されると、ＲＯＩの情報がプロファイル計測部２５に送られる。プロファイル計測部２５は、ＲＯＩ設定部２４によって設定されたＲＯＩ内のプロファイルを計測する（ステップＳ１０６）。
【００４６】
具体的には、プロファイル計測部２５は、位置ずれ補正されたＤＳＡ画像と、もう一方のＤＳＡ画像とを画像メモリ３３から読み出し、読み出したＤＳＡ画像それぞれについて、ＲＯＩ内のＤＳＡ画素値の平均値又は総和値のプロファイルを計測する。図９は、第１の実施形態に係るプロファイル計測部２５による処理の一例を説明するための図である。
【００４７】
例えば、プロファイル計測部２５は、図９に示すように、術前ＤＳＡ画像のプロファイル「ｆ（ｔ）」と、術後ＤＳＡ画像のプロファイル「ｇ（ｔ）」とを計測する。ここで、図９に示す「ΔＴ」は、ＲＯＩに造影剤が到達するまでの遅れ時間を示す。プロファイル計測部２５によって計測されたプロファイルの情報は、補正ファクター同定部２６に送られる。
【００４８】
補正ファクター同定部２６は、ＲＯＩ内のプロファイルの情報を取得すると、術前のＤＳＡ画像及び術後のＤＳＡ画像のＲＯＩ内のプロファイルが略一致する補正ファクターを同定する（ステップＳ１０７）。例えば、補正ファクター同定部２６は、以下に示す式（２）により、術前ＤＳＡ画像のプロファイル「ｆ（ｔ）」と、術後ＤＳＡ画像のプロファイル「ｇ（ｔ）」とを略一致させる補正ファクターを同定する。
【００４９】
【数２】
2014193320.tif 000004
【００５０】
補正ファクター同定部２６は、逐次近似アルゴリズム（algorithm）で、ゲイン「α」、造影剤到達までの遅れ時間「Δｔ」及び血流速度「Ｔ」を変化させながら、「Ｅ」が最小になる「α」、「Δｔ」及び「Ｔ」を探索する。すなわち、補正ファクター同定部２６は、術前、術後の患者の心臓の拍出量の違いを補正するためのゲイン「α」と、術前、術後の造影剤投与時のカテーテルの位置ずれによる造影剤到達までの遅れ時間を補正するための「Δｔ」と、術前、術後の患者の心拍数の違いを補正するための血流速度「Ｔ」を用いて、２つのプロファイルを略一致させる。そして、補正ファクター同定部２６は、探索した「α」、「Δｔ」及び「Ｔ」をＤＳＡ画像補正部２７に送る。
【００５１】
なお、上述した例では、「α」、「Δｔ」及び「Ｔ」のすべての補正ファクターを求める場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、患者の心拍が安定しており、血流速度がほぼ一定と仮定できる場合には、「Ｔ＝１」とする場合であってもよい。このようにすることで、補正ファクターの同定を高速に実行することが可能となる。
【００５２】
ＤＳＡ画像補正部２７は、補正ファクター同定部２６から受付けた補正ファクター「α」、「Δｔ」及び「Ｔ」を用いて、術前のＤＳＡ画像又は術後のＤＳＡ画像を補正する（ステップＳ１０８）。ここで、ＤＳＡ画像の補正は、一般的に、術前のＤＳＡ画像を補正するが、本実施形態においては、どちらのＤＳＡ画像を補正する場合であってもよい。ＤＳＡ画像補正部２７によって補正されたＤＳＡ画像は、画像メモリ３３に記憶される。
【００５３】
なお、ここでは、補正ファクターは術前・術後の血流状態をチェックするＤＳＡ画像を用いて同定された。しかし、ほぼ同時期に撮影された治療などの影響のないＤＳＡ画像を用いて補正ファクターを同定しても良い。具体的には、例えば、右内頸動脈の血管内治療を行った場合、右内頸動脈の術前・術後のＤＳＡ画像を撮影するタイミングの前後で左内頸動脈の術前・術後のＤＳＡ画像を撮影する。この画像には治療などの影響がないと考えられるため、補正ファクターを安定して同定できると言うメリットがある。
【００５４】
ＤＳＡ画像補正部２７によってＤＳＡ画像が補正されると、第２サブトラクション部２９は、補正されたＤＳＡ画像ともう一方のＤＳＡ画像とを画像メモリ３３から読み出し、サブトラクションを実行する（ステップＳ１０９）。例えば、第２サブトラクション部２９は、以下に示す式（３）により、補正されたＤＳＡ画像ともう一方のＤＳＡ画像とをサブトラクションする。
【００５５】
【数３】
2014193320.tif 000005
【００５６】
ここで、式（３）においては、Ｃ（ｉ，ｊ）が血流チェック画像を示す。また、式（３）においては、ｆ(i,j)（ｔ）が術前のＤＳＡ画像を示し、ｇ(i,j)（ｔ）が術後のＤＳＡ画像を示す。また、式（３）においては、術後のＤＳＡ画像が補正された場合の式を示す。第２サブトラクション部２９は、式（３）に示すように、補正された術後のＤＳＡ画像から術前のＤＳＡ画像を差分することで、血流チェック画像を生成する。第２サブトラクション部２９によって生成された血流チェック画像（サブトラクション画像）は、画像メモリ３３に記憶される。
【００５７】
そして、第２サブトラクション部２９によって生成された血流チェック画像は、画像メモリ３３に記憶されると同時に制御部３４に送られ、制御部３４によって表示部４０にてカラーで表示される（ステップＳ１１０）。例えば、表示部４０は、少なくとも一方が補正された２枚のＤＳＡ画像を差分した結果に基づくパーフュージョン画像を表示する。図１０は、第１の実施形態に係る表示部４０にて表示される血流チェック画像の一例を示す図である。
【００５８】
例えば、制御部３４は、図１０に示すように、画素値をカラーに変換して（例えば、血流が十分増加している領域（図１０の領域Ｒ１など）を赤、血流が十分低下している領域（図１０の領域Ｒ２など）を青、中間の領域（図１０の領域Ｒ３など）を黄色に変換して）、変換した画像を表示部４０にて表示させる。なお、表示される血流チェック画像は、白黒の画像であってもよい。
【００５９】
上述したように、第１の実施形態に係るＸ線診断装置１は、術前及び術後の少なくとも一方が補正された２枚のＤＳＡ画像を差分した血流チェック画像を時系列に沿って連続して表示することでカラー又は白黒の動画像を表示する。ここで、第１の実施形態に係るＸ線診断装置１は、上述した動画像の他にも種々の解析結果を表示することができる。例えば、Ｘ線診断装置１は、ＤＳＡ画像の所定の領域ごとに、術前と術後の血流の状態の変化を示すグラフを表示することも可能である。
【００６０】
図１１Ａ〜図１１Ｃは、第１の実施形態に係る表示部にて表示される解析結果の一例を示す図である。図１１Ａ〜図１１Ｃにおいては、図中の（Ａ）に所定の領域における術前及び術後のプロファイルを示し、図中の（Ｂ）に術前のプロファイルと術後のプロファイルとの差分情報を示す。なお、図１１Ａ〜図１１Ｃにおいては、術後のプロファイル「ｇ（ｔ）」から術前のプロファイル「ｆ（ｔ）」を差分した場合の差分情報を示す。
【００６１】
例えば、制御部３４は、主幹動脈を含む領域或いは治療の影響のない毛細血管を含む領域によって同定された補正ファクターによって術前及び術後の少なくとも一方が補正された２枚のＤＳＡ画像の所定の領域について、術前及び術後のＤＳＡ画像のプロファイルと、差分情報を表示部４０にて表示させる。例えば、制御部３４は、図１１Ａの（Ａ）に示すように、術前のＤＳＡ画像における所定の領域のプロファイル「ｆ（ｔ）」と、術後のＤＳＡ画像における同一の領域のプロファイル「ｇ（ｔ）」とを表示部４０にて表示する。さらに、制御部３４は、図１１Ａの（Ｂ）に示すように、図１１Ａの（Ａ）に示す術後のプロファイル「ｇ（ｔ）」から術前のプロファイル「ｆ（ｔ）」を差分したグラフ７１を表示部４０にて表示する。
【００６２】
ここで、図示したプロファイル及びグラフを表示する所定の領域については、操作者によって任意に指定することができる。また、所定のサイズの領域で画像全体を網羅的に表示する場合であってもよい。これによって操作者は、プロファイルやグラフが表示された領域が術前と術後でどのように変化しているのか詳細に観察することが可能になる。例えば、図１１Ａに示すようなプロファイルやグラフが表示された場合には、操作者は、対応する領域において、術後に血流の状態が改善していることを判断することができる。
【００６３】
同様に、操作者は、図１１Ｂ及び図１１Ｃに示すプロファイル及びグラフによって、対応する領域が術前と術後でどのように変化しているのか容易に判断することが可能になる。例えば、表示部４０に図１１Ｂに示すプロファイルやグラフが表示された場合には、操作者は、術後の血流に若干の遅れが生じていると判断することができ、軽い狭窄が生じている恐れを疑うことができる。また、例えば、表示部４０に図１１Ｃに示すプロファイルやグラフが表示された場合には、操作者は、術後の血流が悪くなっていると判断することができ、梗塞が生じている恐れを疑うことができる。
【００６４】
なお、図１１Ａ〜図１１Ｃにおいては、術後のプロファイル「ｇ（ｔ）」から術前のプロファイル「ｆ（ｔ）」を差分したグラフを示す場合について説明したが、実施形態はこれに限定されるものではない。すなわち、術前のプロファイル「ｆ（ｔ）」から術後のプロファイル「ｇ（ｔ）」を差分したグラフを示す場合であってもよい。
【００６５】
また、第１の実施形態に係るＸ線診断装置１は、血流チェック画像と、血流チェック画像の領域に対応するプロファイルなどを表示することも可能である。図１２は、第１の実施形態に係る表示部にて表示される表示情報の一例を示す図である。図１２に示すように、操作者が入力部を介して血流チェック画像上で点や領域を指定すると、制御部３４は、指定された点や領域の術前及び術後のプロファイルを表示させる。なお、図１２においてはプロファイルのみ示しているが、差分のグラフを表示することも可能である。
【００６６】
このように、第１の実施形態に係るＸ線診断装置１においては、血流チェック画像に加えて、所定の領域ごとのプロファイルや差分のグラフなどを表示することができる。さらに、第１の実施形態に係るＸ線診断装置１では、所定の領域の血流量の差を示す血流量差や、術前と術後の血流の違いを示す差異度を表示することも可能である。以下、式（４）〜式（７）を用いて血流量差及び差異度について説明する。なお、式（４）〜式（７）においては、上述した補正ファクター「α」、「Δｔ」及び「Ｔ」のうち、「α」及び「Δｔ」によって補正する場合を一例に挙げて説明する。
【００６７】
すなわち、補正ファクター同定部２６は、プロファイル計測部２５によって計測された術前のＤＳＡ画像のＲＯＩ内のプロファイル「ｆ（ｔ）」と、術後のＤＳＡ画像のＲＯＩ内のプロファイル「ｇ（ｔ）」とが、式（４）となるような「α」及び「Δｔ」を同定する。なお、「α」は、術前、術後の患者の心臓の拍出量の違いを補正するためのゲインであり、「Δｔ」は、術前、術後の造影剤投与時のカテーテルの位置ずれによる造影剤到達までの遅れ時間を補正するためのものである。
【００６８】
【数４】
2014193320.tif 000006
【００６９】
そして、第２サブトラクション部２９は、以下に示す式（５）により、補正されたＤＳＡ画像ともう一方のＤＳＡ画像とをサブトラクションする。ここで、式（５）においては、ＣＯＭＰ（ｔ）がＤＳＡ画像の差分を示す。また、式（５）においては、ｆ(i,j)（ｔ）が術前のＤＳＡ画像を示し、ｇ(i,j)（ｔ）が術後のＤＳＡ画像を示す。また、式（５）においては、術後のＤＳＡ画像が補正された場合の式を示す。
【００７０】
【数５】
2014193320.tif 000007
【００７１】
式（５）に示すように、第２サブトラクション部２９は、補正された術後のＤＳＡ画像と術前のＤＳＡ画像とを画素ごとにサブトラクションする。制御部３４は、第２サブトラクション部２９によって算出されたＣＯＭＰ（ｔ）を用いて、以下に示す式（６）により、血流量差を算出する。すなわち、制御部３４は、補正後の差分画像に基づいて算出された部位ごとのプロファイルを元に、血流量差及び差異度を算出する。ここで、式（６）においては、ＤｉｆＶが血流量差を示す。また、式（６）においては、Ｎが画像サイズを示す。また、式（６）においては、ＣＯＭＰ（ｎ）が画素ごとの差分を示す。
【００７２】
【数６】
2014193320.tif 000008
【００７３】
式（６）に示すように、制御部３４は、画素ごとの差分の平均を血流量の差として算出して、表示部４０に表示する。また、制御部３４は、第２サブトラクション部２９によって算出されたＣＯＭＰ（ｔ）を用いて、以下に示す式（７）により、差異度を算出する。ここで、式（７）においては、ＤｉｆＥが差異度を示す。また、式（７）においては、Ｎが画像サイズを示す。また、式（７）においては、ＣＯＭＰ（ｎ）が画素ごとの差分を示す。
【００７４】
【数７】
2014193320.tif 000009
【００７５】
式（７）に示すように、制御部３４は、画素ごとの差分の２乗の平均を血流の差異度として算出して、表示部４０に表示する。このように、第１の実施形態に係るＸ線診断装置１は、血流チェック画像や、プロファイル及び差分のグラフに加えて、血流量差及び差異度を算出して表示することができる。なお、血流量差及び差異度については、操作者によって状況に応じて使い分けられる。
【００７６】
ここで、上述した血流量差及び差異度の式はあくまでも一例であり、その他の式で算出される場合であってもよい。例えば、術前、術後のプロファイルによって規格化する場合であってもよい。以下、術前、術後のプロファイルによって規格化する場合について、式（８）〜式（１３）を用いて説明する。なお、式（８）〜式（１３）においては、血流量差ＤｉｆＶ及び差異度ＤｉｆＥの右辺のみを示す。
【００７７】
例えば、制御部３４は、以下に示す式（８）により、術前のプロファイルによって規格化した血流量差を算出し、以下に示す式（９）により、術前のプロファイルによって規格化した差異度を算出する。ここで、式（８）及び（９）においては、「ｆ（ｎ）」が術前のプロファイルを示す。ここで、式（８）及び（９）における「β」は、分母が０になることを避けるためのもので、「ｆ（ｎ）」の血管領域の値より十分小さい値が設定される。
【００７８】
【数８】
2014193320.tif 000010
【００７９】
【数９】
2014193320.tif 000011
【００８０】
また、例えば、制御部３４は、以下に示す式（１０）により、術後のプロファイルによって規格化した血流量差を算出し、以下に示す式（１１）により、術後のプロファイルによって規格化した差異度を算出する。ここで、式（１０）及び（１１）においては、「ｇ（ｎ）」が術前のプロファイルを示し、術後の画像が補正ファクターによって補正されている場合について示す。ここで、式（１０）及び（１１）における「β」は、分母が０になることを避けるためのもので、「ｆ（ｎ）」の血管領域の値より十分小さい値が設定される。
【００８１】
【数１０】
2014193320.tif 000012
【００８２】
【数１１】
2014193320.tif 000013
【００８３】
また、例えば、制御部３４は、以下に示す式（１２）により、術前及び術後のプロファイルの平均によって規格化した血流量差を算出し、以下に示す式（１３）により、術前及び術後のプロファイルの平均によって規格化した差異度を算出する。ここで、式（１２）及び（１３）においては、「ｆ（ｎ）」が術前のプロファイルを示し、「ｇ（ｎ）」が術前のプロファイルを示し、術後の画像が補正ファクターによって補正されている場合について示す。ここで、式（１２）及び（１３）における「β」は、分母が０になることを避けるためのもので、「ｆ（ｎ）」の血管領域の値より十分小さい値が設定される。
【００８４】
【数１２】
2014193320.tif 000014
【００８５】
【数１３】
2014193320.tif 000015
【００８６】
上述した例では、補正ファクター「α」及び「Δｔ」によって補正する場合を一例に挙げて説明したが、実施形態はこれに限定されるものではない。すなわち、補正ファクター「α」、「Δｔ」及び「Ｔ」を用いて血流量差及び差異度が算出される場合であってもよい。かかる場合には、補正ファクター同定部２６は、プロファイル計測部２５によって計測された術前のＤＳＡ画像のＲＯＩ内のプロファイル「ｆ（ｔ）」と、術後のＤＳＡ画像のＲＯＩ内のプロファイル「ｇ（ｔ）」とが、以下に示す式（１４）となるような「α」、「Δｔ」及び「Ｔ」を同定する。なお、「α」は、術前、術後の患者の心臓の拍出量の違いを補正するためのゲインであり、「Δｔ」は、術前、術後の造影剤投与時のカテーテルの位置ずれによる造影剤到達までの遅れ時間を補正するためのものである。また、「Ｔ」は、術前、術後の患者の心拍数の違いを補正するための血流速度を示す。
【００８７】
【数１４】
2014193320.tif 000016
【００８８】
式（１４）に示すように、第１の実施形態に係るＸ線診断装置１は、術前、術後の患者の心拍数の違いを補正するための血流速度を補正して術前、術後の血流量差及び差異度を算出することができる。これにより、術前、術後で心拍数が大きく異なる場合であっても、それらを補正したうえで解析することを可能にする。
【００８９】
上述したように、第１の実施形態によれば、プロファイル計測部２５は、被検体の頭部に対して略同一方向から撮影された時相（撮影時期）の異なる２枚のＤＳＡ画像における血管を含む略同一位置にそれぞれ設定されたＲＯＩにおける造影剤濃度に関するプロファイルをそれぞれ計測する。そして、補正ファクター同定部２６は、プロファイル計測部２５によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する。ＤＳＡ画像補正部２７は、補正ファクター同定部２６によって同定された補正ファクターに基づいて、２枚のＤＳＡ画像のうちの少なくとも一方を補正する。そして、制御部３４は、ＤＳＡ画像補正部２７によって少なくとも一方が補正された２枚のＤＳＡ画像に基づいた情報を表示部４０にて表示するように制御する。従って、第１の実施形態に係るＸ線診断装置１は、術前のＤＳＡ画像及び術後のＤＳＡ画像における造影剤による染まりの効果を略一致させた後の差分画像を生成して表示することができ、術前の画像と術後の画像とを正確に比較することを可能にする。
【００９０】
また、第１の実施形態によれば、第１サブトラクション部２８は、被検体の頭部に対して略同一方向から造影剤を注入しながら撮影された時相（撮影時期）の異なる２枚のＸ線動画像を元に、造影剤の影響のある複数の画像と造影剤の影響の殆どない画像とを引き算しＤＳＡ画像を算出する。位置ずれ同定部２２は、２枚のＸ線動画像の中の造影剤の影響が殆どない画像を元に位置ずれを同定する。位置ずれ補正部２３は、同定された位置ずれ情報を元に２枚のＤＳＡ画像の内、少なくとも一方を位置ずれ補正する。プロファイル計測部２５は、少なくとも一方が補正された２枚のＤＳＡ画像における血管を含む略同一位置にそれぞれ設定されたＲＯＩのプロファイルをそれぞれ計測する。補正ファクター同定部２６は、プロファイル計測部２５によってそれぞれ計測された２つのプロファイルが略一致するように補正ファクターを同定する。ＤＳＡ画像補正部２７は、補正ファクター同定部２６によって同定された補正ファクターに基づいて、２枚のＤＳＡ画像うちの少なくとも一方を補正する。制御部３４は、ＤＳＡ画像補正部２７によって少なくとも一方が補正された２枚のＤＳＡ画像に基づいた情報を表示部４０にて表示するように制御する。従って、第１の実施形態に係るＸ線診断装置１は、術前のＤＳＡ画像及び術後のＤＳＡ画像における造影剤による染まりの効果を略一致させた後の差分画像を生成して表示することができ、術前の画像と術後の画像とを正確に比較することを可能にする。
【００９１】
また、第１の実施形態によれば、ＲＯＩは、血管である主幹動脈または治療に影響が無い毛細血管を含む領域である。従って、第１の実施形態に係るＸ線診断装置１は、画像内の種々の領域を用いて補正ファクターを同定することを可能にする。
【００９２】
また、第１の実施形態によれば、ＲＯＩ設定部２４は、２枚のＤＳＡ画像のうち少なくとも一方の画像における主幹動脈の位置を算出し、算出した主幹動脈を含む領域をＲＯＩとして設定する。従って、第１の実施形態に係るＸ線診断装置１は、治療の影響がなく、造影剤濃度の変化が顕著に表れる主幹動脈を含む領域を用いることができ、より正確な補正ファクターを同定することを可能にする。
【００９３】
また、第１の実施形態によれば、表示部４０は、少なくとも一方が補正された２枚のＤＳＡ画像を差分した結果に基づくパーフュージョン画像を表示する。従って、第１の実施形態に係るＸ線診断装置１は、より正確な解析結果を表示することを可能にする。
【００９４】
また、第１の実施形態によれば、補正ファクター同定部２６は、２枚のＤＳＡ画像における２つのプロファイルが略一致するように、造影剤到達までの時間のずれ、ゲイン及び血流速度のうち、少なくとも一つを同定する。従って、第１の実施形態に係るＸ線診断装置１は、造影剤の効果が大きく変化する補正ファクターを用いてＤＳＡ画像を補正することができ、術前のＤＳＡ画像と、術後のＤＳＡ画像の比較をより正確に行うことを可能にする。
【００９５】
また、第１の実施形態によれば、ＲＯＩ設定部２４は、ＤＳＡ画像の下端から所定の範囲内に含まれる血管のうち、血管径が最大となる血管を抽出し、抽出した血管において、部分領域ごとに濃度値の積分値のピーク値を算出し、算出したピーク値が最低となる部分領域をＲＯＩとして設定する。従って、第１の実施形態に係るＸ線診断装置１は、主幹動脈上の造影剤の効果が高い領域をＲＯＩとして設定することができ、正確な比較読影を可能にする。
【００９６】
また、第１の実施形態によれば、ＲＯＩ設定部２４は、ＤＳＡ画像の生成に用いられたマスク画像の濃度値に基づいて、部分領域ごとのＸ線透過度を算出し、算出したＸ線透過度が低い部分領域をＲＯＩの対象外とする。従って、第１の実施形態に係るＸ線診断装置１は、ＲＯＩの設定に係る工程を高速化することを可能にする。
【００９７】
また、第１の実施形態によれば、血流チェック画像（サブトラクション画像）は動画像として作成されているが、動画像の内、ピクセル毎に最大誤差を同定し、最大誤差画像を静止画像として表示しても良い。
【００９８】
（第２の実施形態）
さて、これまで第１の実施形態について説明したが、上述した第１の実施形態以外にも、種々の異なる形態にて実施されてよいものである。
【００９９】
（変形例）
上述した第１の実施形態では、術前のＤＳＡ画像を収集したときと同一のデータ収集プログラムを用いて術後のＤＳＡ画像を収集する場合について説明した。しかしながら、実施形態は、これに限定されるものではなく、例えば、造影剤の染まりの効果を変化させるファクターを任意に選択して、選択したファクターを術前と術後で一致させるようにしてもよい。
【０１００】
かかる場合には、例えば、制御部３４が、Ｘ線撮影条件、撮影角度、ＦＯＶ（Field  Of  View）、ＳＩＤ（Source  Image  Distance）、コリメータの位置、補償フィルターの位置、Ｘ線焦点サイズ、線質調整フィルター、造影剤注入までの時間及び造影条件のうち、少なくとも１つについて、術前のＤＳＡ画像の撮影時と一致させて、術後のＤＳＡ画像を撮影するように制御する。
【０１０１】
上述した第１の実施形態では、ＲＯＩ設定部２４が、ＤＳＡ画像に自動的にＲＯＩを設定する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、ユーザによって設定される場合であってもよい。
【０１０２】
かかる場合には、例えば、ＲＯＩ設定部２４が、ユーザに対してＲＯＩを設定するように要求する。一例を挙げると、ＲＯＩ設定部２４は、表示部４０にＲＯＩの設定を促すメッセージを表示させる。ユーザは、図示しないＧＵＩを介して、ＲＯＩを設定する。このとき、ＲＯＩ設定部２４は、ユーザに対して、造影剤を投与するカテーテルに近く、患部から遠く、さらに血管が検出器面に対して平行になるようにＲＯＩを設定させる。
【０１０３】
上述した第１の実施形態では、Ｘ線診断装置１が血流チェック画像を生成する場合について説明したが、上述した処理は、ワークステーションなどの画像処理装置によって実行される場合であってもよい。かかる場合には、例えば、ネットワークを介してＸ線診断装置や、画像保管装置などと接続されたワークステーションが、Ｘ線診断装置や、画像保管装置などから画像データを取得する。そして、ワークステーションは、取得した画像データを用いて上述した処理を実行する。
【０１０４】
以上説明したとおり、第１〜２の実施形態によれば、本実施形態のＸ線診断装置及び画像処理装置は、術前の画像と術後の画像とを正確に比較することを可能にする。
【０１０５】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１０６】
１  Ｘ線診断装置
２２  位置ずれ同定部
２３  位置ずれ補正部
２４  ＲＯＩ設定部
２４ａ  血管探索部
２４ｂ  最大血管同定部
２４ｃ  最低ピーク領域同定部
２５  プロファイル計測部
２６  補正ファクター同定部
２７  ＤＳＡ画像補正部
２８  第１サブトラクション部
２９  第２サブトラクション部
３４  制御部
１００  画像処理装置
【図１】
2014193320.tif 000017
【図２】
2014193320.tif 000018
【図３】
2014193320.tif 000019
【図４】
2014193320.tif 000020
【図５】
2014193320.tif 000021
【図６】
2014193320.tif 000022
【図７Ａ】
2014193320.tif 000023
【図７Ｂ】
2014193320.tif 000024
【図７Ｃ】
2014193320.tif 000025
【図８】
2014193320.tif 000026
【図９】
2014193320.tif 000027
【図１０】
2014193320.tif 000028
【図１１Ａ】
2014193320.tif 000029
【図１１Ｂ】
2014193320.tif 000030
【図１１Ｃ】
2014193320.tif 000031
【図１２】
2014193320.tif 000032