(57)【要約】      （修正有）【課題】所望の画素値変化率を示す画素を診断対象画像上で容易に発見することができる画像診断支援装置等を提供する。【解決手段】ＣＰＵは、時刻Ｔ０〜Ｔ２にそれぞれ撮像されたダイナミック造影ＭＲＩ画像に基づいて、診断対象画像をビューワ２０に表示する。ＣＰＵはまた、早期相および遅延相に対応する２つの座標軸の各々に画素値変化率が割り振られたＳＣＰ用グラフをビューワ２０に表示する。ＣＰＵはさらに、診断対象画像上の複数の画素にそれぞれ対応しかつ早期相および遅延相の各々における画素値変化率を示す複数の変化率セットを検出し、検出した複数の変化率セットにそれぞれ対応する複数のドットをＳＣＰ用グラフに描画する。所望のドットがキーボード／マウスによってクリックされると、ＣＰＵは、クリックされたドットに対応する画素を診断対象画像上でガイドする。【選択図】図１１
2017153852.tif 000002
【特許請求の範囲】
【請求項１】
複数の時相を定義する複数の時刻にそれぞれ撮像されかつ検査部位の共通の断層を表す複数の断層画像に基づいて診断対象画像を表示する第１表示手段、
前記複数の時相にそれぞれ対応する複数の座標軸の各々に画素値変化率が割り振られた変化率グラフを表示する第２表示手段、
前記第１表示手段によって表示された診断対象画像上の複数の画素にそれぞれ対応しかつ前記複数の時相の各々における画素値変化率を示す複数の変化率セットを検出する検出手段、
前記検出手段によって検出された複数の変化率セットにそれぞれ対応する複数の指標を前記第２表示手段によって表示された変化率グラフに描画する指標描画手段、
前記指標描画手段によって描画された複数の指標のいずれか１つを選択する指標選択操作を受け付ける受け付け手段、および
前記指標選択操作によって選択された指標に対応する画素を前記第１表示手段によって表示された診断対象画像上でガイドするガイド手段を備える、画像診断支援装置。
【請求項２】
前記第１表示手段によって表示された診断対象画像上で解析領域を定義する定義手段をさらに備え、
前記検出手段は前記定義手段によって定義された解析領域に属する画素を対象として処理を実行する、請求項１記載の画像診断支援装置。
【請求項３】
前記第１表示手段によって表示された診断対象画像に第１カーソルを多重表示する第１多重表示手段をさらに備え、
前記定義手段は前記第１多重表示手段によって多重表示された第１カーソルが指示する画素を含む一部の領域を前記解析領域として定義する、請求項２記載の画像診断支援装置。
【請求項４】
前記診断対象画像は既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップを含み、
前記定義手段は前記第１カーソルが指示する画素が有色画素であるときに処理を実行し、
前記解析領域は前記第１カーソルが指示する有色画素から連続する複数の有色画素によって形成された領域である、請求項３記載の画像診断支援装置。
【請求項５】
前記第１多重表示手段によって多重表示された第１カーソルが指示する画素に対応する指標を前記指標描画手段によって描画された複数の指標の中から特定する指標特定手段、および
前記指標特定手段によって特定された指標を指示する第２カーソルを前記第２表示手段によって表示された変化率グラフに多重表示する第２多重表示手段をさらに備える、請求項３または４記載の画像診断支援装置。
【請求項６】
前記第２多重表示手段によって多重表示された第２カーソルをカーソル移動操作に従って移動させるカーソル移動手段をさらに備え、
前記指標選択操作によって選択される指標は前記第２多重表示手段によって多重表示された第２カーソルが前記指標選択操作の時点で指示する指標に相当する、請求項５記載の画像診断支援装置。
【請求項７】
時間および画素値係数がそれぞれ割り振られた２つの座標軸を有する画素値係数グラフを表示する第３表示手段、および
前記第２多重表示手段によって多重表示された第２カーソルが指示する指標に対応する画素値の時間的な変化特性を前記第３表示手段によって表示された画素値係数グラフに描画する変化特性描画手段をさらに備える、請求項５または６記載の画像診断支援装置。
【請求項８】
前記複数の時刻は、前記検査部位に対する造影剤の注入前後の時刻であり、かつ前記造影剤が注入された後に到来する基準時刻を含み、
前記複数の時相は、前記基準時刻よりも前の第１時相、および前記基準時刻よりも後の第２時相を含む、請求項１ないし７のいずれかに記載の画像診断支援装置。
【請求項９】
画像診断支援装置のプロセッサによって実行される画像診断支援方法であって、
複数の時相を定義する複数の時刻にそれぞれ撮像されかつ検査部位の共通の断層を表す複数の断層画像に基づいて診断対象画像を表示する第１表示ステップ、
前記複数の時相にそれぞれ対応する複数の座標軸の各々に画素値変化率が割り振られた変化率グラフを表示する第２表示ステップ、
前記第１表示ステップによって表示された診断対象画像上の複数の画素にそれぞれ対応しかつ前記複数の時相の各々における画素値変化率を示す複数の変化率セットを検出する検出ステップ、
前記検出ステップによって検出された複数の変化率セットにそれぞれ対応する複数の指標を前記第２表示ステップによって表示された変化率グラフに描画する指標描画ステップ、
前記指標描画ステップによって描画された複数の指標のいずれか１つを選択する指標選択操作を受け付ける受け付けステップ、および
前記指標選択操作によって選択された指標に対応する画素を前記第１表示ステップによって表示された診断対象画像上でガイドするガイドステップを備える、画像診断支援方法。
【請求項１０】
画像診断支援装置のプロセッサに、
複数の時相を定義する複数の時刻にそれぞれ撮像されかつ検査部位の共通の断層を表す複数の断層画像に基づいて診断対象画像を表示する第１表示ステップ、
前記複数の時相にそれぞれ対応する複数の座標軸の各々に画素値変化率が割り振られた変化率グラフを表示する第２表示ステップ、
前記第１表示ステップによって表示された診断対象画像上の複数の画素にそれぞれ対応しかつ前記複数の時相の各々における画素値変化率を示す複数の変化率セットを検出する検出ステップ、
前記検出ステップによって検出された複数の変化率セットにそれぞれ対応する複数の指標を前記第２表示ステップによって表示された変化率グラフに描画する指標描画ステップ、
前記指標描画ステップによって描画された複数の指標のいずれか１つを選択する指標選択操作を受け付ける受け付けステップ、および
前記指標選択操作によって選択された指標に対応する画素を前記第１表示ステップによって表示された診断対象画像上でガイドするガイドステップを実行させるための、画像診断支援プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、画像診断支援装置、画像診断支援方法および画像診断支援プログラムに関し、特に、モダリティによって撮像された断層画像に基づく診断医の画像診断を支援する、画像診断支援装置、画像診断支援方法および画像診断支援プログラムに関する。
【背景技術】
【０００２】
CT、MRIなどで実施されるDynamic Studyでは、病変部内部に設定した関心領域における造影剤の継時的変化を計測・評価することが、病変内部の性状を把握する上で非常に重要である。この評価には従来からTime Intensity Curve（以下、TICと呼ぶ。）計測が用いられており、TICの軌跡に注目してその評価を行ってきた。近年のモダリティ装置の性能向上に伴い、造影剤の継時的変化を微小関心領域（数画素）で計測し病変内部の性状をより細かく把握するようになり、さらに最近では画素単位で性状を把握しその結果を病理画像と対比する試みが始まっている。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】中原浩著、「CAEを用いた乳房MRI診断」、乳癌BOOK２０１４、RadFan２０１４年７月臨時増刊、p.79-p.81
【非特許文献２】株式会社AZE著、「灌流から診るMRI乳がん解析ソフトウェア」、月間インナービジョン２０１４年８月号、p.77
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、病変内部の造影剤の継時的変化を画素単位でTIC計測しその血流特性を表示する機能はあるが、逆に計測・評価後の血流特性から診断画像上の画素を特定できる機能がないため、血流特性と診断画像上の画素を１対１に対応付けて把握することが難しく、改善が求められていた。
【０００５】
それゆえに、この発明の主たる目的は、特定の血流特性を示す微小検査部位を断層画像上で容易に発見することができる、画像診断支援装置、画像診断支援方法および画像診断支援プログラムを提供することである。
【課題を解決するための手段】
【０００６】
この発明に従う画像診断支援装置(10：実施例で相当する参照符号。以下同じ)は、複数の時相を定義する複数の時刻にそれぞれ撮像されかつ検査部位の共通の断層を表す複数の断層画像(IMt0~IMt2)に基づいて診断対象画像を表示する第１表示手段(S21~S33)、複数の時相にそれぞれ対応する複数の座標軸の各々に画素値変化率が割り振られた変化率グラフを表示する第２表示手段(S35)、第１表示手段によって表示された診断対象画像上の複数の画素にそれぞれ対応しかつ複数の時相の各々における画素値変化率を示す複数の変化率セットを検出する検出手段(S65~S69, S73~S75)、検出手段によって検出された複数の変化率セットにそれぞれ対応する複数の指標を第２表示手段によって表示された変化率グラフに描画する指標描画手段(S71)、指標描画手段によって描画された複数の指標のいずれか１つを選択する指標選択操作を受け付ける受け付け手段(S89)、および指標選択操作によって選択された指標に対応する画素を第１表示手段によって表示された診断対象画像上でガイドするガイド手段(S103)を備える。
【０００７】
好ましくは、第１表示手段によって表示された診断対象画像上で解析領域を定義する定義手段(S59~S63)がさらに備えられ、検出手段は定義手段によって定義された解析領域に属する画素を対象として処理を実行する。
【０００８】
さらに好ましくは、第１表示手段によって表示された診断対象画像に第１カーソル(CS1)を多重表示する第１多重表示手段(S41)がさらに備えられ、定義手段は第１多重表示手段によって多重表示された第１カーソルが指示する画素を含む一部の領域を解析領域として定義する。
【０００９】
或る局面では、診断対象画像は既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップ(MPt012)を含み、定義手段は第１カーソルが指示する画素が有色画素であるときに処理を実行し、解析領域は第１カーソルが指示する有色画素から連続する複数の有色画素によって形成された領域である。
【００１０】
他の局面では、第１多重表示手段によって多重表示された第１カーソルが指示する画素に対応する指標を指標描画手段によって描画された複数の指標の中から特定する指標特定手段(S79)、および指標特定手段によって特定された指標を指示する第２カーソル(CS2)を第２表示手段によって表示された変化率グラフに多重表示する第２多重表示手段(S81)がさらに備えられる。
【００１１】
或る例では、第２多重表示手段によって多重表示された第２カーソルをカーソル移動操作に従って移動させるカーソル移動手段(S87, S91)がさらに備えられ、指標選択操作によって選択される指標は第２多重表示手段によって多重表示された第２カーソルが指標選択操作の時点で指示する指標に相当する。
【００１２】
他の例では、時間および画素値係数がそれぞれ割り振られた２つの座標軸を有する画素値係数グラフを表示する第３表示手段(S37, S39)、および第２多重表示手段によって多重表示された第２カーソルが指示する指標に対応する画素値の時間的な変化特性を第３表示手段によって表示された画素値係数グラフに描画する変化特性描画手段(S85, S95, S107)がさらに備えられる。
【００１３】
好ましくは、複数の時刻は、検査部位に対する造影剤の注入前後の時刻であり、かつ造影剤が注入された後に到来する基準時刻を含み、複数の時相は、基準時刻よりも前の第１時相、および基準時刻よりも後の第２時相を含む。
【００１４】
この発明に従う画像診断支援方法は、診断支援装置(10)のプロセッサ(14)によって実行される診断支援方法であって、複数の時相を定義する複数の時刻にそれぞれ撮像されかつ検査部位の共通の断層を表す複数の断層画像(IMt0~IMt2)に基づいて診断対象画像を表示する第１表示ステップ(S21~S33)、複数の時相にそれぞれ対応する複数の座標軸の各々に画素値変化率が割り振られた変化率グラフを表示する第２表示ステップ(S35)、第１表示ステップによって表示された診断対象画像上の複数の画素にそれぞれ対応しかつ複数の時相の各々における画素値変化率を示す複数の変化率セットを検出する検出ステップ(S65~S69, S73~S75)、検出ステップによって検出された複数の変化率セットにそれぞれ対応する複数の指標を第２表示ステップによって表示された変化率グラフに描画する指標描画ステップ(S71)、指標描画ステップによって描画された複数の指標のいずれか１つを選択する指標選択操作を受け付ける受け付けステップ(S89)、および指標選択操作によって選択された指標に対応する画素を第１表示ステップによって表示された診断対象画像上でガイドするガイドステップ(S103)を備える。
【００１５】
この発明に従う画像診断支援プログラムは、診断支援装置(10)のプロセッサ(14)に、複数の時相を定義する複数の時刻にそれぞれ撮像されかつ検査部位の共通の断層を表す複数の断層画像(IMt0~IMt2)に基づいて診断対象画像を表示する第１表示ステップ(S21~S33)、複数の時相にそれぞれ対応する複数の座標軸の各々に画素値変化率が割り振られた変化率グラフを表示する第２表示ステップ(S35)、第１表示ステップによって表示された診断対象画像上の複数の画素にそれぞれ対応しかつ複数の時相の各々における画素値変化率を示す複数の変化率セットを検出する検出ステップ(S65~S69, S73~S75)、検出ステップによって検出された複数の変化率セットにそれぞれ対応する複数の指標を第２表示ステップによって表示された変化率グラフに描画する指標描画ステップ(S71)、指標描画ステップによって描画された複数の指標のいずれか１つを選択する指標選択操作を受け付ける受け付けステップ(S89)、および指標選択操作によって選択された指標に対応する画素を第１表示ステップによって表示された診断対象画像上でガイドするガイドステップ(S103)を実行させるための、診断支援プログラムである。
【発明の効果】
【００１６】
変化率グラフを構成する複数の座標軸は複数の時相にそれぞれ対応し、各座標軸には画素値変化率が割り振られる。また、変化率グラフに描画される複数の指標は診断対象画像上の複数の画素にそれぞれ対応し、各指標は複数の時相における対象画素の画素値変化率を示す。変化率グラフに描画された複数の指標のいずれか１つが指標選択操作によって選択されると、選択された指標に対応する画素が診断対象画像上でポイントされる。これによって、所望の画素値変化率を示す画素（＝特定の血流特性を示す微小検査部位）を診断対象画像上で容易に発見することができる。
【００１７】
この発明の上述の目的，その他の目的，特徴および利点は、図面を参照して行う以下の実施例の詳細な説明から一層明らかとなろう。
【図面の簡単な説明】
【００１８】
【図１】この実施例に適用される医療診断支援システムの構成の一例を示すブロック図である。
【図２】（Ａ）は時刻Ｔ０にダイナミック造影ＭＲＩシーケンスで撮像された複数のダイナミック造影ＭＲＩ画像をそれぞれ収めた複数のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｂ）は時刻Ｔ１にダイナミック造影ＭＲＩシーケンスで撮像された複数のダイナミック造影ＭＲＩ画像をそれぞれ収めた複数のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｃ）は時刻Ｔ２にダイナミック造影ＭＲＩシーケンスで撮像された複数のダイナミック造影ＭＲＩ画像をそれぞれ収めた複数のＤＩＣＯＭファイルの構造の一例を示す図解図である。
【図３】この実施例の画像診断支援装置を構成するＤＲＡＭに設けられたビューワエリアおよびワークエリアのマッピング状態の一例を示す図解図である。
【図４】ＳＣＰ用グラフの一例を示す図解図である。
【図５】カラーマップとこれに多重されたカーソル，ＳＣＰ用グラフ，ＫＣ用グラフ，ＴＩＣ用グラフがビューワに表示された状態の一例を示す図解図である。
【図６】（Ａ）は拡大操作に応答してカラーマップを拡大した状態の一例を示す図解図であり、（Ｂ）はパン操作に応答してカラーマップをパンした状態の一例を示す図解図であり、（Ｃ）はカーソル移動操作に応答してカーソルを移動させた状態の一例を示す図解図であり、（Ｄ）はクリック操作に応答して解析領域が定義された状態の一例を示す図解図である。
【図７】ＳＣＰ用グラフに複数のドットが描画された状態の一例を示す図解図である。
【図８】カラーマップとこれに多重されたカーソル，ＳＣＰ用グラフとこれに描画された複数のドット，ＫＣ用グラフ，ＴＩＣ用グラフがビューワに表示された状態の一例を示す図解図である。
【図９】カラーマップ，ＳＣＰ用グラフとこれに描画ないし多重された複数のドットおよびカーソル，ＫＣ用グラフとこれに描画されたKinetic Curve，ＴＩＣ用グラフがビューワに表示された状態の一例を示す図解図である。
【図１０】カラーマップ，ＳＣＰ用グラフとこれに描画ないし多重された複数のドットおよびカーソル，ＫＣ用グラフとこれに描画されたKinetic Curve，ＴＩＣ用グラフがビューワに表示された状態の他の一例を示す図解図である。
【図１１】カラーマップとこれに多重されたＲＯＩガイド，ＳＣＰ用グラフとこれに描画ないし多重された複数のドット，カーソルおよびＲＯＩガイド，ＫＣ用グラフとこれに描画されたKinetic Curve，ＴＩＣ用グラフとこれに描画されたTime Intensity Curveがビューワに表示された状態の一例を示す図解図である。
【図１２】検査レポートの一例を示す図解図である。
【図１３】この実施例の画像診断支援装置を構成するＣＰＵの動作の一部を示すフロー図である。
【図１４】この実施例の画像診断支援装置を構成するＣＰＵの動作の他の一部を示すフロー図である。
【図１５】この実施例の画像診断支援装置を構成するＣＰＵの動作のその他の一部を示すフロー図である。
【図１６】この実施例の画像診断支援装置を構成するＣＰＵの動作のさらにその他の一部を示すフロー図である。
【図１７】この実施例の画像診断支援装置を構成するＣＰＵの動作の他の一部を示すフロー図である。
【図１８】この実施例の画像診断支援装置を構成するＣＰＵの動作のその他の一部を示すフロー図である。
【図１９】ＳＣＰ用グラフの変形例を示す図解図である。
【図２０】ＳＣＰ用グラフの他の変形例を示す図解図である。
【発明を実施するための形態】
【００１９】
図１を参照して、この実施例の医療診断支援システムは、院内ＬＡＮ５０を介して互いに接続された画像診断支援装置１０，ＭＲＩ装置（モダリティ）３０および画像サーバ４０によって構成される。画像診断支援装置１０において、バスＢＳ１には、通信Ｉ／Ｆ１２，ＣＰＵ１４，キーボード／マウス１６，ＨＤＤ１８，ビューワ２０およびＤＲＡＭ２２が接続される。
【００２０】
ＭＲＩ装置３０は、核磁気共鳴現象を利用して、検査部位である被検者の乳房の断層画像を作成する装置である。ダイナミック造影ＭＲＩシーケンスに従う乳房の撮像処理は、体内に造影剤を注入する前の時刻Ｔ０に１回実行され（これを“プレーン撮像”と呼ぶ）、造影剤の注入が開始された後の互いに異なる時刻Ｔ１およびＴ２に１回ずつ実行される。
【００２１】
このうち、時刻Ｔ１を基準時刻とすると、時刻Ｔ０から時刻Ｔ１までの時相が早期相（第１時相）に相当し、時刻Ｔ１から時刻Ｔ２までの時相が遅延相（第２時相）に相当する。
【００２２】
時刻Ｔ０〜Ｔ２の各々での撮像処理の結果、“水平４００画素×垂直４００画素”のマトリクスサイズと“０．８ｍｍ”の画素サイズとを有するダイナミック造影ＭＲＩ画像が、０．８ｍｍ間隔で配された１２５個の水平スライス位置の各々に対応して作成される。作成された各々のダイナミック造影ＭＲＩ画像は、ファイル情報（ファイルサイズ，マトリクスサイズ，画素サイズ，スライス間隔），患者情報，撮像時刻，スライス位置とともにＤＩＣＯＭファイルに収められる。
【００２３】
この結果、時刻Ｔ０では図２（Ａ）に示す１２５個のＤＩＣＯＭファイルが作成され、時刻Ｔ１では図２（Ｂ）に示す１２５個のＤＩＣＯＭファイルが作成され、時刻Ｔ２では図２（Ｃ）に示す１２５個のＤＩＣＯＭファイルが作成される。作成されたＤＩＣＯＭファイルはいずれも、院内ＬＡＮ５０を経て画像サーバ４０に保存される。
【００２４】
画像診断支援装置１０に設けられた通信Ｉ／Ｆ１２は、こうして画像サーバ４０に蓄積された合計３７５個のＤＩＣＯＭファイルを院内ＬＡＮ５０を通して取得する。取得されたＤＩＣＯＭファイルは、バスＢＳ１を介してＨＤＤ１８に保存される。
【００２５】
ＣＰＵ１４は、画像診断医による画像診断を支援するべく、ＨＤＤ１８に保存されたＤＩＣＯＭファイルを対象とする処理を図１３〜図１８に示すフロー図に従って実行する。なお、このフロー図に相当するプログラムは、画像診断支援プログラムとしてＨＤＤ１８に記憶される。
【００２６】
図１３を参照して、ステップＳ１ではダイナミック造影ＭＲＩ画像のスライス位置を識別する変数Ｋを初期値に設定し、ステップＳ３では検査画像をビューワ２０に表示する。ステップＳ３の処理は、具体的には図１４に示すサブルーチンに従う。
【００２７】
ステップＳ２１では時刻Ｔ０に撮像されたＫ番目のスライス位置のダイナミック造影ＭＲＩ画像を取得し、ステップＳ２３では時刻Ｔ１に撮像されたＫ番目のスライス位置のダイナミック造影ＭＲＩ画像を取得し、ステップＳ２５では時刻Ｔ２に撮像されたＫ番目のスライス位置のダイナミック造影ＭＲＩ画像を取得する。
【００２８】
以下では、ステップＳ２１で取得したダイナミック造影ＭＲＩ画像を“ＩＭｔ０”と定義し、ステップＳ２３で取得したダイナミック造影ＭＲＩ画像を“ＩＭｔ１”と定義し、ステップＳ２５で取得したダイナミック造影ＭＲＩ画像を“ＩＭｔ２”と定義する。
【００２９】
また、図３に示すように、ＤＲＡＭ２２には、ワークエリア２２ｗが設けられる。ステップＳ２１〜Ｓ２５で取得したダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２は、ＤＲＡＭ２２のワークエリア２２ｗに保持される。
【００３０】
ステップＳ２７では、ワークエリア２２ｗ上のダイナミック造影ＭＲＩ画像ＩＭｔ０およびＩＭｔ１に基づいてカラーマップＭＰｔ０１を作成する。ステップＳ２９では、ワークエリア２２ｗ上のダイナミック造影ＭＲＩ画像ＩＭｔ１およびＩＭｔ２に基づいてカラーマップＭＰｔ１２を作成する。作成されたカラーマップＭＰｔ０１およびＭＰｔ１２もまた、ワークエリア２２ｗに保持される。
【００３１】
ここで、カラーマップＭＰｔ０１およびＭＰｔ１２はいずれも、画素値に応じて異なる平面的色分布を表すマップである。カラーマップＭＰｔ０１を構成する色（＝端的には色相）は、ダイナミック造影ＭＲＩ画像ＩＭｔ０上の画素値とダイナミック造影ＭＲＩ画像ＩＭｔ１上の対応する画素値とを比較して決定される。同様に、カラーマップＭＰｔ１２を構成する色は、ダイナミック造影ＭＲＩ画像ＩＭｔ１上の画素値とダイナミック造影ＭＲＩ画像ＩＭｔ２上の対応する画素値とを比較して決定される。
【００３２】
また、カラーマップＭＰｔ０１およびＭＰｔ１２のいずれについても、各画素に色を付けるか否かは、ダイナミック造影ＭＲＩ画像ＩＭｔ０の画素値を比較対象とした、ダイナミック造影ＭＲＩ画像ＩＭｔ１の画素値の増加率および増加量に基づいて判断する。具体的には、増加率が４０％以上でかつ増加量が閾値（＝増加率が３９％以下の画素値のうち最大画素値）を上回るという条件が課せられ、この条件を満足する画素に対して色が付される。
【００３３】
したがって、付される色はカラーマップＭＰｔ０１およびＭＰｔ１２の間で相違するものの、色の分布状態はカラーマップＭＰｔ０１およびＭＰｔ１２の間で一致する。
【００３４】
参考までに、有色画素（＝色が付された画素）は腫瘍の染まり具合の形状を表すところ、この形状が円形，卵型，分葉状，不整形のいずれであるかによって良悪性の診断結果が変わる。また、早期相特性（早期相における血流特性）が“Fast”に該当しかつ遅延相特性（遅延相における血流特性）が“Washout”に該当する画素に対応する微小検査部位が最も悪性であると言われている。
【００３５】
ステップＳ３１では、ステップＳ２７およびＳ２９で作成されたカラーマップＭＰｔ０１およびＭＰｔ１２に基づいてカラーマップＭＰｔ０１２を作成する。カラーマップＭＰｔ０１２もまた、画素値に応じて異なる平面的色分布を表すマップである。カラーマップＭＰｔ０１２の詳細は後述する。
【００３６】
ステップＳ３３では、ダイナミック造影ＭＲＩ画像ＩＭｔ０にカラーマップＭＰｔ０１２を多重し、これによって得られた合成画像を“診断対象画像”としてＤＲＡＭ２２上のビューワエリア２２ｖに展開する。図３に示すように、ビューワエリア２２ｖは、４つの小エリアを有する。診断対象画像は、４つの小エリアのうち左上に配された小エリアに展開される。
【００３７】
なお、ステップＳ３３では、後述するステップＳ４７またはＳ５３の処理に備えるべく、診断対象画像の表示倍率を“１．０”に設定し、診断対象画像の切り出し位置を“中央”に設定する。
【００３８】
ステップＳ３５では、ビューワエリア２２ｖの左下に配された小エリアにＳＣＰ用グラフ（変化率グラフ）を展開する。ステップＳ３７では、ビューワエリア２２ｖの右上に配された小エリアにＫＣ用グラフを展開する。ステップＳ３９では、ビューワエリア２２ｖの右下に配された小エリアにＴＩＣ用グラフを展開する。
【００３９】
ここで、“ＳＣＰ”はScatter Plotの略であり、ＳＣＰ用グラフは図４に示すように構成される。図４によれば、早期相における画素値変化率が縦軸に割り振られ、遅延相における画素値変化率が横軸に割り振られる。また、縦軸には早期相特性を識別する“Fast”，“Medium”，“Slow”の文字列が付され、横軸には遅延相特性を識別する“Persistent”，“Plateau”，“Washout”の文字列が付される。
【００４０】
さらに、縦軸はステップＳ２７で作成したカラーマップＭＰｔ０１に対応し、横軸はステップＳ２９で作成したカラーマップＭＰｔ１２に対応する。つまり、ＳＣＰ用グラフ上のカラーは、縦軸近傍ではカラーマップＭＰｔ０１と同等のカラーとなり、横軸近傍ではカラーマップＭＰｔ１２と同等のカラーとなる。また、縦軸および横軸から離れた場所のカラーは、その座標に応じてカラーマップＭＰｔ０１およびＭＰｔ１２に重み付けされたカラーとなる。
【００４１】
早期相特性が“Fast”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化率が１００％以上の画素である。早期相特性が“Medium”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化率が５０％〜９９％の画素である。早期相特性が“Slow”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化率が４９％以下の画素である。
【００４２】
遅延相特性が“Washout”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化率が−１０％以下の画素である。遅延相特性が“Plateau”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化率が−１０％〜１０％の画素である。遅延相特性が“Persistent”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化率が１０％以上の画素である。
【００４３】
また、“ＫＣ”はKinetic Curveであり、“ＴＩＣ”はTime Intensity Curveの略である。ＫＣ用グラフおよびＴＩＣ用グラフのいずれも、BI-RADS 規格に沿って作成され、横軸に時間が割り振られる一方、縦軸に画素値係数が割り振られる。より詳しくは、ＫＣ用グラフの縦軸には画素値変化率が割り振られ、ＴＩＣ用グラフの縦軸には画素値が割り振られる。
【００４４】
なお、ＫＣ用グラフに描画されるKinetic Curveは、対象画素の画素値変化率の時間的な変化特性を示し、ＴＩＣ用グラフに描画されるTime Intensity Curveは、対象画素の画素値の時間的変化特性を示す。
【００４５】
ビューワエリア２２ｖに展開された診断対象画像，ＳＣＰ用グラフ，ＫＣ用フラグおよびＴＩＣ用フラグは、ビューワ２０によって読み出され、図５に示す要領で画面に表示される。検査画像は、こうして表示された診断対象画像，ＳＣＰ用グラフ，ＫＣ用フラグおよびＴＩＣ用フラグによって構成される。ステップＳ３９の処理が完了すると、上階層のルーチンに復帰する。
【００４６】
図１３に戻って、ステップＳ５ではキーボード／マウス１６によってページング操作が行われたか否かを判別し、ステップＳ７ではキーボード／マウス１６によって解析開始操作が行われたか否かを判別する。
【００４７】
ステップＳ５の判別結果およびステップＳ７の判別結果のいずれもがＮＯであれば、ステップＳ５に戻る。ステップＳ５の判別結果がＹＥＳであれば、ステップＳ９で変数Ｋを更新し、その後にステップＳ３に戻る。ステップＳ７の判別結果がＹＥＳであれば、ステップＳ１１で解析＆表示処理を実行し、その後にステップＳ５に戻る。
【００４８】
ステップＳ１１の解析＆表示処理は、図１５〜図１８に示すサブルーチンに従って実行される。
【００４９】
図１５に示すステップＳ４１では、カーソル多重表示命令をビューワ２０に与える。ビューワ２０は、与えられた命令に従って、カーソルＣＳ１を診断対象画像に多重表示する（図５参照）。
【００５０】
ステップＳ４３ではキーボード／マウス１６によって拡大／縮小操作が行われたか否かを判別し、ステップＳ４５では拡大／縮小操作によって指定された表示倍率が可変範囲に属するか否かを判別する。ステップＳ４３の判別結果およびステップＳ４５の判別結果のいずれか一方がＮＯであればそのままステップＳ４９に進み、ステップＳ４３の判別結果およびステップＳ４５の判別結果の両方がＹＥＳであればステップＳ４７で以下の処理を実行してからステップＳ４９に進む。
【００５１】
つまり、ステップＳ４７では、拡大／縮小操作に従って表示倍率を更新し、更新された表示倍率および現時点の切り出し位置の設定に従う一部の画像をワークエリア２２ｗ上の診断対象画像から切り出す。
【００５２】
ステップＳ４７ではまた、切り出された画像のサイズを調整し、調整後のサイズを有する画像を更新後の診断対象画像としてビューワエリア２２ｖの左上の小エリアに上書きする。したがって、図５に示す表示状態で拡大操作が行われると、ビューワ２０の左上に表示された診断対象画像は図６（Ａ）に示すように更新される。
【００５３】
ステップＳ４９ではキーボード／マウス１６によってパン／チルト操作が行われたか否かを判別し、ステップＳ５１ではパン／チルト操作によって指定された切り出し位置が切り出し可能範囲に属するか否かを判別する。ステップＳ５１の判別では、現時点の表示倍率も参照する。
【００５４】
ステップＳ４９の判別結果およびステップＳ５１の判別結果のいずれか一方がＮＯであればそのままステップＳ５５に進み、ステップＳ４９の判別結果およびステップＳ５１の判別結果の両方がＹＥＳであればステップＳ５３で以下の処理を実行してからステップＳ５５に進む。
【００５５】
つまり、ステップＳ５３では、パン／チルト操作に従って切り出し位置を更新し、更新された切り出し位置および現時点の表示倍率の設定に従う一部の画像をワークエリア２２ｗ上の診断対象画像から切り出す。
【００５６】
ステップＳ５３ではまた、切り出された画像のサイズを調整し、調整後のサイズを有する画像を更新後の診断対象画像としてビューワエリア２２ｖの左上の小エリアに上書きする。したがって、図６（Ａ）に示す表示状態でパン操作が行われると、ビューワ２０の左上に表示された診断対象画像は図６（Ｂ）に示すように更新される。
【００５７】
ステップＳ５５ではキーボード／マウス１６によってカーソル移動操作が行われたか否かを判別し、判別結果がＮＯであればそのままステップＳ５９に進む一方、判別結果がＹＥＳであればステップＳ５７でカーソル移動命令をビューワ２０に与えてからステップＳ４７に進む。したがって、図６（Ｂ）に示す表示状態でカーソル移動操作が行われた場合、カーソルＣＳ１はたとえば図６（Ｃ）に示す位置に移動する。
【００５８】
ステップＳ５９では、キーボード／マウス１６によってクリック操作が行われたか否かを判別する。また、ステップＳ６１では、クリックされた画素（＝カーソルＣＳ１が指示する画素）が有色画素であるか否かを判別する。ステップＳ５９の判別結果およびステップＳ６１の判別結果のいずれか一方がＮＯであればステップＳ４３に戻り、ステップＳ５９の判別結果およびステップＳ６１の判別結果の両方がＹＥＳであればステップＳ６３に進む。
【００５９】
ステップＳ６３では、カーソルＣＳ１が指示する有色画素から平面的に連続する複数の有色画素によって形成された領域を“解析領域”として定義する。ステップＳ６３ではまた、定義された解析領域の輪郭に沿うように、ビューワエリア２２ｖ上の診断対象画像に線を描画する。描画された線は、ビューワ２０の表示にも反映される（図６（Ｄ）参照）。
【００６０】
ステップＳ６５では、定義された解析領域に属する複数の有色画素の中からラスタ方向における先頭の有色画素を特定し、特定した有色画素に対応する画素をダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２の各々の上で指定する。
【００６１】
ステップＳ６７では、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での指定画素値の変化率を早期相における画素値変化率として検出する。ステップＳ６９では、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での指定画素値の変化率を遅延相における画素値変化率として検出する。以下では、こうして検出された２つの画素値変化率を“変化率セット”と定義する。
【００６２】
ステップＳ７１では、ステップＳ６７およびＳ６９の処理によって得られた変化率セットに対応する座標をＳＣＰ用グラフ上で検出し、検出した座標にドット（指標）を描画する。ステップＳ７３では、指定画素が解析領域に属する複数の有色画素のうちラスタ方向における末尾の有色画素に対応する画素であるか否かを判別する。判別結果がＮＯであればステップＳ７５に進み、指定画素を解析領域においてラスタ方向の次画素に更新する。更新が完了すると、ステップＳ６７に戻る。
【００６３】
この結果、ステップＳ７３の判別結果がＹＥＳとなるまでステップＳ６７〜Ｓ６９およびＳ７５の処理が繰り返され、ＳＣＰ用グラフには図７に示す要領で複数のドットが描画される。つまり、ステップＳ５９〜Ｓ７３の処理が実行されることで、ビューワ２０の表示は図８に遷移する。
【００６４】
ステップＳ７３の判別結果がＹＥＳであればステップＳ７７に進み、カーソルＣＳ１の非表示をビューワ２０に命令する。この結果、カーソルＣＳ１が診断対象画像上から消失する。
【００６５】
ステップＳ７９では、クリック画素（＝消失前にカーソルＣＳ１が指示していた有色画素）に対応するドットをＳＣＰ用グラフに描画された複数のドットの中から特定する。ステップＳ８１では、特定したドットを指示する位置へのカーソルＣＳ２の表示をビューワ２０に命令する。ビューワ２０は、命令に従ってカーソルＣＳ２をＳＣＰ用グラフに多重表示する（図９の左下参照）。
【００６６】
ステップＳ８３では、カーソルＣＳ２が指示するドットに対応する２つの画素値変化率（＝変化率セット）を検出し、検出した２つの画素値変化率を示す文字をビューワエリア２２ｖの左下の小エリアおよび右上の小エリアに描画する。この結果、描画された文字が図９に示す要領でＳＣＰ用グラフおよびＫＣ用グラフに多重表示される。図９によれば、カーソルＣＳ２が指示するドットについては、早期相における画素値変化率が“９８％”であり、遅延相における画素値変化率は“−３％”である。
【００６７】
ステップＳ８５では、ステップＳ８３で検出された２つの画素値変化率を参照してKinetic Curveを作成し、作成されたKinetic Curveをビューワエリア２２ｖの右上の小エリアに描画する。この結果、描画されたKinetic Curveが図９に示す要領でＫＣ用グラフに多重表示される。
【００６８】
ステップＳ８７ではキーボード／マウス１６によってカーソル移動操作が行われたか否かを判別し、ステップＳ８９ではキーボード／マウス１６によってクリック操作（指標選択操作）が行われたか否かを判別する。ステップＳ８７の判別結果およびステップＳ８９の判別結果のいずれもがＮＯであればステップＳ８７に戻り、ステップＳ８７の判別結果がＹＥＳであればステップＳ９１に進み、ステップＳ８９の判別結果がＹＥＳであればステップＳ９９に進む。
【００６９】
ステップＳ９１では、カーソルＣＳ２の移動をビューワ２０に命令する。この結果、カーソルＣＳ２が命令に従う方向に移動する。ステップＳ９３では、移動後のカーソルＣＳ２がＳＣＰ用グラフ上の複数のドットのいずれか１つを指示しているか否かを判別する。判別結果がＹＥＳであればステップＳ９５に進み、判別結果がＮＯであればステップＳ９７に進む。
【００７０】
ステップＳ９５では、ＳＣＰ用グラフおよびＫＣ用グラフの各々に明示された２つの画素値変化率をカーソルＣＳ２が指示するドットに対応する２つの画素値変化率に更新し、更新された２つの画素値変化率に従ってＫＣ用グラフ上のKinetic Curveを更新する。ビューワ２０の表示は、図９から図１０に遷移する。ステップＳ９５の処理が完了すると、ステップＳ８７に戻る。
【００７１】
ステップＳ９７では、ＳＣＰ用グラフおよびＫＣ用グラフの各々に明示された２つの画素値変化率とＫＣグラフに描画されたKinetic Curveとを非表示とする。ステップＳ９７の処理が完了すると、ステップＳ８７に戻る。
【００７２】
図１８に示すステップＳ９９では、カーソルＣＳ２がＳＣＰ用グラフ上の複数のドットのいずれか１つを指示しているか否かを判別する。判別結果がＮＯであればステップＳ８７に戻り、判別結果がＹＥＳであればステップＳ１０１に進む。
【００７３】
ステップＳ１０１では、カーソルＣＳ２が指示するドットを強調するＲＯＩガイドをビューワエリア２２ｖの左下の小エリアに描画する。描画されたＲＯＩガイドは、図１１に示す要領でビューワ２０上のＳＣＰグラフに多重表示される。図１１によれば、ＳＣＰ用グラフ上のＲＯＩガイドは、カーソルＣＳ２が指示するドットから引き出された直線と“Ｒ１”の文字列とによって構成される。
【００７４】
ステップＳ１０３では、カーソルＣＳ２が指示するドットに対応する有色画素を診断対象画像上で検出し、検出した有色画素を中心とする５つの有色画素の領域を“関心領域”として定義する。ステップＳ１０３では、こうして定義された関心領域を強調するＲＯＩガイドをビューワエリア２２ｖの左上の小エリアに描画する。
【００７５】
描画されたＲＯＩガイドは、図１１に示す要領でビューワ２０上の診断対象画像に多重表示される。図１１によれば、診断対象画像上のＲＯＩガイドは、関心領域を囲む円とこの円から引き出された直線と“Ｒ１”の文字列によって構成される。
【００７６】
ステップＳ１０７では、関心領域に属する５つの有色画素にそれぞれ対応する５つの変化率セットを検出し、検出した５つの変化率セットを参照してTime Intensity Curveを作成し、作成したTime Intensity Curveをビューワエリア２２ｖの右下の小エリアに描画する。描画されたTime Intensity Curveは、図１１に示す要領でビューワ２０上のＴＩＣ用グラフに多重表示される。
【００７７】
ステップＳ１０９では、ＫＣ用グラフとこれに多重表示された文字およびKinetic Curveとを非表示とするべく、ビューワエリア２２ｖの右上の小エリアをクリアする。ステップＳ１１１では、カーソルＣＳ２が指示するドットに対応する２つの画素値変化率（早期相：１６７％，遅延相：−２３％）と、関心領域を構成する画素数（＝５画素）と、関心領域の中心座標（２６６，１４６）とを表す解析情報をビューワエリア２２ｖの右上の小エリアに描画する。解析情報は、ビューワ２０の右上に図１２に示す要領で表示される。
【００７８】
ステップＳ１１３ではキーボード／マウス１６によって復帰操作が行われたか否かを判別し、ステップＳ１１５ではキーボード／マウス１６によって保存操作が行われたか否かを判別する。ステップＳ１１３の判別結果およびステップＳ１１５の判別結果のいずれもがＮＯであればステップＳ１１３に戻り、ステップＳ１１３の判別結果がＹＥＳであればステップＳ１１９に進み、ステップＳ１１５の判別結果がＹＥＳであればステップＳ１１７に進む。
【００７９】
ステップＳ１１７では、ビューワ２０の表示に相当する検査レポートをワークエリア２２ｗ上で作成し、作成した検査レポートの保存を通信Ｉ／Ｆ１２に命令する。通信Ｉ／Ｆ１２は、ワークエリア２２ｗから検査レポートを読み出し、読み出された検査レポートを院内ＬＡＮ５０を介して画像サーバ４０に転送する。転送された検査レポートは、対象患者に関連付けられた状態で画像サーバ４０に保存される。ステップＳ１１７の処理が完了すると、ステップＳ１１９に進む。
【００８０】
ステップＳ１１９では、解析情報が表示される前の画面に復帰するべく、上述したステップＳ３７およびＳ８５の処理と同様の処理を実行する。この結果、ビューワ２０の表示は図１２から図１１に遷移する。ステップＳ１１９の処理が完了すると、ステップＳ８７に戻る。
【００８１】
以上の説明から分かるように、ＣＰＵ１４は、早期相および遅延相を定義する時刻Ｔ０〜Ｔ２にそれぞれ撮像されかつ検査部位の共通の断層を表すダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２に基づいて、診断対象画像をビューワ２０に表示する(S21~S33)。ＣＰＵ１４はまた、早期相および遅延相にそれぞれ対応する２つの座標軸の各々に画素値変化率が割り振られたＳＣＰ用グラフをビューワ２０に表示する(S35)。
【００８２】
ＣＰＵ１４はさらに、診断対象画像上の複数の画素にそれぞれ対応しかつ早期相および遅延相の各々における画素値変化率を示す複数の変化率セットを検出し(S65~S69, S73~S75)、検出した複数の変化率セットにそれぞれ対応する複数のドットをＳＣＰ用グラフに描画する(S71)。こうして描画された複数のドットのいずれか１つがキーボード／マウス１６によってクリックされると、ＣＰＵ１４は、クリックされたドットに対応する画素を診断対象画像上でガイドする(S103)。
【００８３】
これによって、所望の画素値変化率を示す画素（＝特定の血流特性を示す微小検査部位）を診断対象画像上で容易に発見することができる。
【００８４】
なお、この実施例では、図４に示すグラフをＳＣＰ用グラフとして採用しているが、これに代えて図１９に示すグラフまたは図２０に示すグラフをＳＣＰ用グラフとして採用するようにしてもよい。
【００８５】
図１９によれば、早期相特性を識別する“Fast”，“Medium”，“Slow”の文字および遅延相特性を識別する“Persistent”，“Plateau”，“Washout”の文字の各々に、画素値変化率の範囲を示す数字が付される。
【００８６】
また、図２０によれば、早期相特性を識別する“Fast”，“Medium”，“Slow”の文字に代えて“多”，“中”，“少”の文字が付され、遅延相特性を識別する“Persistent”，“Plateau”，“Washout”の文字に代えて“増加”，“変化なし”，“減少”の文字が付される。
【符号の説明】
【００８７】
１０  …画像診断支援装置
１４  …ＣＰＵ
１６  …キーボード／マウス
２０  …ビューワ
２２  …ＤＲＡＭ
３０  …ＭＲＩ装置
４０  …画像サーバ
５０  …院内ＬＡＮ
【図１】
2017153852.tif 000003
【図２】
2017153852.tif 000004
【図３】
2017153852.tif 000005
【図４】
2017153852.tif 000006
【図５】
2017153852.tif 000007
【図６】
2017153852.tif 000008
【図７】
2017153852.tif 000009
【図８】
2017153852.tif 000010
【図９】
2017153852.tif 000011
【図１０】
2017153852.tif 000012
【図１１】
2017153852.tif 000013
【図１２】
2017153852.tif 000014
【図１３】
2017153852.tif 000015
【図１４】
2017153852.tif 000016
【図１５】
2017153852.tif 000017
【図１６】
2017153852.tif 000018
【図１７】
2017153852.tif 000019
【図１８】
2017153852.tif 000020
【図１９】
2017153852.tif 000021
【図２０】
2017153852.tif 000022