JP2016059582A - 診断支援装置、診断支援方法および診断支援プログラム - Google Patents

Abstract

【課題】ＭＲＩ画像（断層画像）を取得して乳房などの患部を診断する際の利便性を高めることができる、診断支援装置、診断支援方法および診断支援プログラムを提供する。【解決手段】通信Ｉ／Ｆ１２は、時刻Ｔ０〜Ｔｎの各々でＭＲＩ装置３０が撮像した患部のＭＲＩ画像を取得する。ＣＰＵ１４は、初期相および遅延相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表すカラーマップを通信Ｉ／Ｆ１２によって取得されたＭＲＩ画像に基づいて画面ＳＣＲに表示する。ＣＰＵ１４はポインタ移動操作に従って断層上の所望の位置を共通的に指向するポインタをカラーマップにそれぞれ多重する。ＣＰＵ１４はさらに、初期相および遅延相にそれぞれ対応しかつポインタ移動操作に従う位置における画素値の時間的変化を各々が表すグラフを通信Ｉ／Ｆ１２によって取得されたＭＲＩ画像に基づいて画面ＳＣＲに表示する。【選択図】図１

Description

この発明は、診断支援装置、診断支援方法および診断支援プログラムに関し、特に、乳房の断層を撮像する乳房ＭＲＩ(Magnetic Resonance Imaging)装置のような断層撮像装置からＭＲＩ画像（断層画像）を取得して乳房などの患部を診断する、診断支援装置、診断支援方法および診断支援プログラムに関する。

乳がんの画像診断方式には、マンモグラフィ方式、超音波方式、ＭＲＩ方式がある。このうち、マンモグラフィ方式および超音波方式は、視触診とともに乳がん検診に用いられている。これに対して、ＭＲＩ方式に関しては、検査に造影剤が使用されることもあり、乳がん検診に用いることは推奨されていない（「乳がん発症ハイリスクグループに対する乳房ＭＲＩスクリーニングに関するガイドラインVer.1.0」日本乳癌検診学会 2012）。

しかし、がんの進行度や広がりなどを調べる診断目的としては、海外では乳房ＭＲＩの有用性は広く認識されており、日本でもエビデンス不足からその有用性が認められていなかった経緯はあるものの、２０１３年に日本乳癌学会の乳癌診療ガイドラインで推奨グレードＢ（「エビデンス（有効な根拠・証拠）があり、推奨内容を日常診療で実践するように推奨する」）とされ、現在はその有用性が認識されるようになった。

また、有用性の認識が広まるとともに、米国で生まれた乳房ＭＲＩ検査のガイドライン「BI-RADS-MRI」が注目されるようになり、現在は日本版BI-RADS-MRIに則った検査方法、読影方法、レポート記述方法が多くの施設で採用されるようになってきた。

このBI-RADS-MRIが推奨する診断レポートに記載する内容の１つに「Kinetic curve assessment , Signal intensity / time curve description」という項目があり、「Initial phase（初期相）」と「Delayed phase（遅延相）」の評価を行うことが推奨されている。この評価のための検査をダイナミックＭＲマンモグラフィという。

なお、造影剤の注入開始時刻を“Ｔ０”とし、時刻Ｔ０から６０秒〜１２０秒が経過した時刻を“Ｔ１”とし、時刻Ｔ１から３００秒〜４２０秒が経過した時刻を“Ｔｎ”とすると、時刻Ｔ０から時刻Ｔ１までの時相が“初期相”に相当し、時刻Ｔ１から時刻Ｔｎまでの時相が“遅延相”に相当する。また、初期相における信号強度の変化特性は“Slow”、“Medium”、“Fast”の３つに分類され、遅延相における信号強度の変化特性は“Persistent”、“Plateau”、“Washout”の３つに分類される。

中原浩著，「ＣＡＥを用いた乳房ＭＲＩ診断」，乳癌ＢＯＯＫ２０１４，ＲａｄＦａｎ２０１４年７月臨時増刊，ｐ．７９−ｐ．８１ 株式会社ＡＺＥ著，「灌流から診るＭＲＩ乳がん解析ソフトウェア」，月刊インナービジョン２０１４年８月号，ｐ．７７

乳房ＭＲＩ検査においてBI-RADS-MRIに沿った読影を行うには、Kinetic curve assessmentにおいて、初期相および遅延相の各々における信号強度の変化特性を簡便かつ同時に評価できることが理想的である。

しかし、現状ではKinetic curve assessmentは従来のTime Intensity Curve計測で行われているため、計測に時間がかかり簡便ではない。また、断層画像を形成する画素毎の信号強度（画素値）の変化量または変化率を用いてカラーマップ化し、Time Intensity Curve計測を簡便にするアプリケーションはあるが、初期相および遅延相を同時に評価することは難しい。この結果、現状では、患部を診断する際の利便性の限界がある。

それゆえに、この発明の主たる目的は、患部を診断する際の利便性を高めることができる、診断支援装置、診断支援方法および診断支援プログラムを提供することである。

この発明に従う診断支援装置(10：実施例で相当する参照符号。以下同じ)は、複数の時刻(T0~Tn)の各々で断層撮像装置(30)が撮像した患部の断層画像を取得する取得手段(12)、複数の時相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表す複数のカラーマップ(MP01, MP1n)を取得手段によって取得された断層画像に基づいて表示するカラーマップ表示手段(S5, S9, S21, S31)、ポインタ移動操作に従って断層上の所望の位置を共通的に指向する複数のポインタ(PT01, PT1n)をカラーマップ表示手段によって表示された複数のカラーマップにそれぞれ多重するポインタ多重手段(S33, S43, S45)、および複数の時相にそれぞれ対応しかつポインタ移動操作に従う位置における画素値の時間的変化を各々が表す複数のグラフ(G01, G1n)を取得手段によって取得された断層画像に基づいて表示するグラフ表示手段(S15, S35~S37, S63~S65)を備える。

好ましくは、グラフの横軸および縦軸には時間および画素値がそれぞれ割り当てられ、

画素値に応じて異なる色を示すカラーバーをグラフの縦軸に沿って表示するカラーバー表示手段(S17)がさらに備えられる。

好ましくは、複数の時相をそれぞれ代表する複数の代表画像(IM01, IM1n)を取得手段によって取得された断層画像に基づいて表示する代表画像表示手段(S3, S7, S19, S29)がさらに備えられ、カラーマップ表示手段は代表画像表示手段によって表示された複数の代表画像に複数のカラーマップをそれぞれ多重するカラーマップ多重手段(S21, S31)を含む。

好ましくは、複数の時相は初期相および遅延相を含み、複数のグラフは、初期相に対応しかつ初期値からの画素値の変化を表す第１グラフ(G01)、および遅延相に対応しかつ初期値からの画素値の変化を表す第２グラフ(G1n)を含む。

さらに好ましくは、断層撮像装置は患部への造影剤の注入を想定した装置であり、複数の時刻は患部に造影剤を注入する前の時刻を基準時刻(T0)として含み、グラフ表示手段は、基準時刻と造影剤の注入を開始した時刻(T0’)との差分を補償する補償手段(S63)、および補償手段による補償を踏まえて第１グラフを作成する第１グラフ作成手段(S65)を含む。

好ましくは、取得手段が複数の時刻の各々に対応して取得する断層画像は複数のスライス位置の各々の断層を表し、カラーマップ表示手段およびグラフ表示手段は複数のスライス位置のうちの共通のスライス位置を対象として処理を実行する。

さらに好ましくは、カラーマップ表示手段およびグラフ表示手段が対象とするスライス位置をスライス位置変更操作に従って変更するスライス位置変更手段(S23~S25, S39~S41)がさらに備えられる。

この発明に従う診断支援方法は、複数の時刻(T0~Tn)の各々で断層撮像装置(30)が撮像した患部の断層画像を取得する取得手段(12)を備える診断支援装置(10)によって実行される診断支援方法であって、複数の時相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表す複数のカラーマップ(MP01, MP1n)を取得手段によって取得された断層画像に基づいて表示するカラーマップ表示ステップ(S5, S9, S21, S31)、ポインタ移動操作に従って断層上の所望の位置を共通的に指向する複数のポインタ(PT01, PT1n)をカラーマップ表示ステップによって表示された複数のカラーマップにそれぞれ多重するポインタ多重ステップ(S33, S43, S45)、および複数の時相にそれぞれ対応しかつポインタ移動操作に従う位置における画素値の時間的変化を各々が表す複数のグラフ(G01, G1n)を取得手段によって取得された断層画像に基づいて表示するグラフ表示ステップ(S15, S35~S37, S63~S65)を備える。

この発明に従う診断支援プログラムは、複数の時刻(T0~Tn)の各々で断層撮像装置(30)が撮像した患部の断層画像を取得する取得手段(12)を備える診断支援装置(10)のプロセッサ(14)に、複数の時相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表す複数のカラーマップ(MP01, MP1n)を取得手段によって取得された断層画像に基づいて表示するカラーマップ表示ステップ(S5, S9, S21, S31)、ポインタ移動操作に従って断層上の所望の位置を共通的に指向する複数のポインタ(PT01, PT1n)をカラーマップ表示ステップによって表示された複数のカラーマップにそれぞれ多重するポインタ多重ステップ(S33, S43, S45)、および複数の時相にそれぞれ対応しかつポインタ移動操作に従う位置における画素値の時間的変化を各々が表す複数のグラフ(G01, G1n)を取得手段によって取得された断層画像に基づいて表示するグラフ表示ステップ(S15, S35~S37, S63~S65)を実行させるための、診断支援プログラムである。

患部を診断するには複数の時相における血流状態を把握する必要があるところ、複数のカラーマップおよび複数のグラフはいずれも、複数の時相にそれぞれ対応して患部の血流状態を表す。具体的には、カラーマップは画素値に応じて異なる平面的色分布を表し、グラフは断層上の或る位置における画素値の時間的変化を表す。これによって、複数の時相の各々における血流状態を複数の視点で把握することができる。

さらに、グラフを表示する際に注目される位置は、複数のカラーマップに共通的に多重表示された複数のポインタによって指向され、かつ各ポインタの指向位置はポインタ移動操作に従って変更される。この結果、グラフの表示態様は、ポインタ移動操作に従って変化する。

ポインタの指向位置およびグラフの表示態様の動的な変化は、カラーマップおよびグラフが互いに異なる視点で血流状態を表すことと相俟って、患部を診断する際の利便性の向上に繋がる。

この発明の上述の目的，その他の目的，特徴および利点は、図面を参照して行う以下の実施例の詳細な説明から一層明らかとなろう。

この実施例に適用される医療診断支援システムの構成の一例を示すブロック図である。 （Ａ）は時刻Ｔ０に収集された２００スライスのＭＲＩ画像をそれぞれ収めた２００個のＤＩＣＯＭ(Digital Imaging and COmmunication in Medicine)ファイルの構造の一例を示す図解図であり、（Ｂ）は時刻Ｔ１に収集された２００スライスのＭＲＩ画像をそれぞれ収めた２００個のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｃ）は時刻Ｔ２に収集された２００スライスのＭＲＩ画像をそれぞれ収めた２００個のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｄ）は時刻Ｔ３に収集された２００スライスのＭＲＩ画像をそれぞれ収めた２００個のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｅ）は時刻Ｔｎに収集された２００スライスのＭＲＩ画像をそれぞれ収めた２００個のＤＩＣＯＭファイルの構造の一例を示す図解図である。 （Ａ）は先頭スライス位置のＭＲＩ画像に補完処理を施して作成された初期相を代表する補完画像の一例を示す図解図であり、（Ｂ）は先頭スライス位置のＭＲＩ画像に補完処理を施して作成された遅延相を代表する補完画像の一例を示す図解図である。 （Ａ）は図３（Ａ）に示す補完画像に基づいて作成されたカラーマップの一例を示す図解図であり、（Ｂ）は図３（Ｂ）に示す補完画像に基づいて作成されたカラーマップの一例を示す図解図である。 この実施例の診断支援装置に設けられたディスプレイ装置に表示される診断画像の一例を示す図解図である。 （Ａ）は所望のスライス位置のＭＲＩ画像に補完処理を施して作成された初期相を代表する補完画像の一例を示す図解図であり、（Ｂ）は所望のスライス位置のＭＲＩ画像に補完処理を施して作成された遅延相を代表する補完画像の一例を示す図解図である。 （Ａ）は図６（Ａ）に示す補完画像に基づいて作成されたカラーマップの一例を示す図解図であり、（Ｂ）は図６（Ｂ）に示す補完画像に基づいて作成されたカラーマップの一例を示す図解図である。 この実施例の診断支援装置に設けられたディスプレイ装置に表示される診断画像の他の一例を示す図解図である。 ポインタによって指定された画素の値の初期相における時間的変化の一例を示すグラフである。 ポインタによって指定された画素の値の初期相および遅延相における時間的変化の一例を示すグラフである。 この実施例の診断支援装置に設けられたＣＰＵの動作の一部を示すフロー図である。 この実施例の診断支援装置に設けられたＣＰＵの動作の他の一部を示すフロー図である。 この実施例の診断支援装置に設けられたＣＰＵの動作のその他の一部を示すフロー図である。 他の実施例においてディスプレイ装置に表示されるグラフ表示メニューの設定の一例を示す図解図である。 図１４に示す設定に従うグラフの表示態様の一例を示す図解図である。 他の実施例においてディスプレイ装置に表示されるグラフ表示メニューの設定の他の一例を示す図解図である。 図１６に示す設定に従うグラフの表示態様の一例を示す図解図である。 他の実施例の診断支援装置に設けられたＣＰＵの動作の一部を示すフロー図である。

［第１実施例］
図１を参照して、この実施例の医療診断支援システムは、通信ネットワーク４０を介して互いに接続された診断支援装置１０およびＭＲＩ装置３０によって構成される。診断支援装置１０において、バスＢＳ１には、通信Ｉ／Ｆ１２，ＣＰＵ１４，キーボード１６，マウス１８，ディスプレイ装置２０，ＤＲＡＭ２２，および主記憶装置２４が接続される。ディスプレイ装置２０には、画面ＳＣＲが設けられる。

ＭＲＩ装置３０は、核磁気共鳴現象を利用して患者の乳房のＭＲＩ画像を作成する装置である。撮像は、造影剤を乳房に注入する前に１回実行され（これを“プレーン撮像”と呼ぶ）、造影剤の注入が開始された後は６０秒に１回の割合で合計４回実行される。造影剤の注入開始時刻を“Ｔ０”と定義すると、その後に６０秒毎に到来する４つの撮像時刻はそれぞれ“Ｔ１”，“Ｔ２”，“Ｔ３”，“Ｔｎ”と定義される。

なお、各回の撮像は数秒を要し、時刻Ｔ０〜Ｔｎの各々は厳密には撮像開始時刻を示す。また、この実施例では、造影剤の注入開始時刻がプレーン撮像時刻と一致することを前提としており、時刻Ｔ０はプレーン撮像時刻でもある。さらに、時刻Ｔ０から時刻Ｔ１までの時相は“初期相”と呼ばれ、時刻Ｔ１から時刻Ｔｎまでの時相は“遅延相”と呼ばれる。

ＭＲＩ装置３０では、互いに異なる２００個のスライス位置Ｐ０〜Ｐ１９９にそれぞれ対応する２００フレームのＭＲＩ画像Ｆ０〜Ｆ１９９が、時刻Ｔ０〜Ｔｎの各々で作成される。作成されたＭＲＩ画像をなす各画素の値（＝信号強度）は、撮像時にスライス位置の断層を流れる造影剤の濃度つまりは血液濃度を反映する。

ＭＲＩ装置３０は、こうして作成されたＭＲＩ画像をファイルサイズ，患者情報，撮像時刻，スライス位置とともにＤＩＣＯＭファイルに収める。なお、撮像時刻は、ＤＩＣＯＭ規格で規定されている収集時刻（０００８，００３２）のタグで管理される。

この結果、時刻Ｔ０では図２（Ａ）に示す２００個のＤＩＣＯＭファイルが作成され、時刻Ｔ１では図２（Ｂ）に示す２００個のＤＩＣＯＭファイルが作成され、時刻Ｔ２では図２（Ｃ）に示す２００個のＤＩＣＯＭファイルが作成される。また、時刻Ｔ３では図２（Ｄ）に示す２００個のＤＩＣＯＭファイルが作成され、時刻Ｔｎでは図２（Ｅ）に示す２００個のＤＩＣＯＭファイルが作成される。

診断支援装置１０に設けられた通信Ｉ／Ｆ１２は、こうして作成された合計１０００個のＤＩＣＯＭファイルを通信ネットワーク４０を通して診断支援装置１０から取得する。取得されたＤＩＣＯＭファイルは、バスＢＳ１を介して主記憶装置２４に保存される。

なお、ＤＩＣＯＭファイルは、図示しない外部サーバによって管理するようにしてもよい。この場合、診断支援装置１０は、所望のＤＩＣＯＭファイルを所望のタイミングで外部サーバから取得する。

ＣＰＵ１４は、主記憶装置２４に保存されたＤＩＣＯＭファイルに基づいて診断画像をディスプレイ装置２０の画面ＳＣＲに表示するべく、図１１〜図１３に示すフロー図に従う処理を実行する。なお、このフロー図に相当するプログラムは、診断支援プログラムとして主記憶装置２４に記憶される。

まず、スライス位置を識別する変数ＫをステップＳ１で“０”に設定する。ステップＳ３では、時刻Ｔ０およびＴ１の各々におけるスライス位置ＰＫの断層を表すＭＲＩ画像を主記憶装置２４からＤＲＡＭ２２に転送し、転送された２フレームのＭＲＩ画像に補完処理を施して補完画像ＩＭ０１を作成する。ステップＳ５では、ステップＳ３でＤＲＡＭ２２に転送された２フレームのＭＲＩ画像に基づいて、スライス位置ＰＫに対応するカラーマップＭＰ０１を作成する。

ステップＳ７では、時刻Ｔ１〜Ｔｎの各々におけるスライス位置ＰＫの断層を表すＭＲＩ画像を主記憶装置２４からＤＲＡＭ２２に転送し、転送された４フレームのＭＲＩ画像に補完処理を施して補完画像ＩＭ１ｎを作成する。ステップＳ９では、ステップＳ７でＤＲＡＭ２２に転送された４フレームのＭＲＩ画像に基づいて、スライス位置ＰＫに対応するカラーマップＭＰ１ｎを作成する。

補完画像ＩＭ０１は初期相を代表する代表画像に相当し、補完画像ＩＭ１ｎは遅延相を代表する代表画像に相当する。また、カラーマップＭＰ０１は初期相に対応し、補完画像ＩＭ０１をなす各画素の値（信号強度）に対応する平面的色分布を表す。カラーマップＭＰ１ｎは遅延相に対応し、補完画像ＩＭ１ｎをなす各画素の値（信号強度）に対応する平面的色分布を表す。

ステップＳ１１では変数Ｋが“１９９”に達したか否かを判別し、判別結果がＮＯであればステップＳ１３で変数ＫをインクリメントしてからステップＳ３に戻る一方、判別結果がＹＥＳであればステップＳ１５に進む。したがって、ステップＳ３〜Ｓ９の処理は２００回繰り返され、この結果、補完画像ＩＭ０１およびＩＭ１ｎとカラーマップＭＰ０１およびＭＰ１ｎとがスライス位置Ｐ０〜Ｐ１９９の各々に対応して作成される。

ステップＳ１５では、グラフの枠をＤＲＡＭ２２に描画し、この枠の表示をディスプレイ装置２０に要求する。また、ステップＳ１７では、カラーバーをＤＲＡＭ２２に描画し、このカラーバーの表示をディスプレイ装置２０に要求する。

ディスプレイ装置２０は、グラフの枠およびカラーバーをＤＲＡＭ２２から読み出し、読み出された枠およびカラーバーを図５に示す要領で画面ＳＣＲの上段に表示する。図５によれば、グラフの枠は画面ＳＣＲの上段中央に表示され、カラーバーはグラフの枠の左側に縦軸に沿って表示される。

グラフの枠の横軸および縦軸には、時間および画素値がそれぞれ割り当てられる。詳しくは、横軸には、上述した時刻Ｔ０を起点として時刻Ｔｎまでの時間が割り当てられる。また、縦軸には、時刻Ｔ０における対象スライス位置のＭＲＩ画像をなす対象画素の値を初期値として、３００％までの値が割り当てられる。また、縦軸に沿って表示されるカラーバーは、画素値に応じて異なる色を示す。

ステップＳ１９では、先頭スライス位置（＝スライス位置Ｐ０）に対応する補完画像ＩＭ０１およびＩＭ１ｎの表示をディスプレイ装置２０に要求する。また、ステップＳ２１では、先頭スライス位置に対応するカラーマップＭＰ０１およびＭＰ１ｎの多重表示をディスプレイ装置２０に要求する。

ディスプレイ装置２０は、対象の補完画像ＩＭ０１およびＩＭ１ｎをＤＲＡＭ２２から読み出し、読み出された補完画像ＩＭ０１およびＩＭ１ｎを画面ＳＣＲに表示する。ディスプレイ装置２０はまた、対象のカラーマップＭＰ０１およびＭＰ１ｎをＤＲＡＭ２２から読み出し、読み出されたカラーマップＭＰ０１およびＭＰ１ｎを画面ＳＣＲに表示された補完画像ＩＭ０１およびＩＭ１ｎにそれぞれ多重する。

先頭スライス位置の補完画像ＩＭ０１およびＩＭ１ｎが図３（Ａ）および図３（Ｂ）に示す画像をなし、先頭スライス位置のカラーマップＭＰ０１およびＭＰ１ｎが図４（Ａ）および図４（Ｂ）に示すマップをなす場合、補完画像ＩＭ０１〜ＩＭ１ｎ，カラーマップＭＰ０１〜ＭＰ１ｎは、図５に示す要領で画面ＳＣＲに表示される。つまり、補完画像ＩＭ０１およびＩＭ１ｎは画面ＳＣＲの下段の位置に左右並んで表示され、カラーマップＭＰ０１およびＭＰ１ｎは補完画像ＩＭ０１およびＩＭ１ｎにそれぞれ多重される。

なお、初期相に対応するカラーマップＭＰ０１の色は、カラーバーの色と同期する。この結果、カラーマップＭＰ０１をカラーバーと照らし合わせることで、初期相における各画素値の変化率の把握が容易になる。

ステップＳ２３では、キーボード１６および／またはマウス１８によってスライス位置変更操作が行われたか否かを繰り返し判別する。判別結果がＮＯからＹＥＳに更新されると、ステップＳ２５で変数Ｋの値を変更し、ステップＳ２７でポインタＰＴ０１およびＰＴ１ｎの表示位置を初期化する。初期化によって、ポインタＰＴ０１およびＰＴ１ｎは断層上の初期位置を指向する。

ステップＳ２９では、スライス位置ＰＫに対応する補完画像ＩＭ０１およびＩＭ１ｎの表示をディスプレイ装置２０に要求する。また、ステップＳ３１では、スライス位置ＰＫに対応するカラーマップＭＰ０１およびＭＰ１ｎの多重表示をディスプレイ装置２０に要求する。

ディスプレイ装置２０は、対象の補完画像ＩＭ０１およびＩＭ１ｎをＤＲＡＭ２２から読み出し、読み出された補完画像ＩＭ０１およびＩＭ１ｎを画面ＳＣＲに表示する。ディスプレイ装置２０はまた、対象のカラーマップＭＰ０１およびＭＰ１ｎをＤＲＡＭ２２から読み出し、読み出されたカラーマップＭＰ０１およびＭＰ１ｎを画面ＳＣＲに表示された補完画像ＩＭ０１およびＩＭ１ｎにそれぞれ多重する。

ステップＳ３３では、ポインタＰＴ０１およびＰＴ１ｎの各々を表すキャラクタ画像（以下、単に「ポインタＰＴ０１」または「ポインタＰＴ１ｎ」と呼ぶ。）をＤＲＡＭ２２に展開し、ポインタＰＴ０１およびＰＴ１ｎの表示をディスプレイ装置２０に要求する。ディスプレイ装置２０は、ポインタＰＴ０１およびＰＴ１ｎをＤＲＡＭ２２から読み出し、読み出されたポインタＰＴ０１およびＰＴ１ｎを画面ＳＣＲ上のカラーマップＭＰ０１およびＭＰ１ｎに多重表示する。

ポインタＰＴ０１およびＰＴ１ｎは断層上の共通の位置を指向し、かつ指向位置はステップＳ２７または後述するステップＳ４５の処理に従う。したがって、カラーマップＭＰ０１およびＭＰ１ｎに共通の座標を割り当てた場合、ポインタＰＴ０１の多重位置を示す座標値は、ポインタＰＴ１ｎの多重位置を示す座標値と一致する。

スライス位置ＰＫの補完画像ＩＭ０１およびＩＭ１ｎが図６（Ａ）および図６（Ｂ）に示す画像をなし、スライス位置ＰＫのカラーマップＭＰ０１およびＭＰ１ｎが図７（Ａ）および図７（Ｂ）に示すマップをなす場合、補完画像ＩＭ０１〜ＩＭ１ｎ，カラーマップＭＰ０１〜ＭＰ１ｎは、図８に示す要領で画面ＳＣＲに代替表示される。ポインタＰＴ０１およびＰＴ１ｎは、断層上の共通の位置を指向するように、カラーマップＭＰ０１およびＭＰ１ｎにそれぞれ多重される。

ステップＳ３５では、時刻Ｔ０〜Ｔｎの各々におけるスライス位置ＰＫの断層を表すＭＲＩ画像を主記憶装置２４からＤＲＡＭ２２に転送し、転送された６フレームのＭＲＩ画像の各々からポインタＰＴ０１またはＰＴ１ｎが指向する画素の値を抽出する。ステップＳ３５ではまた、時刻Ｔ０およびＴ１にそれぞれ対応する２つの画素値に基づいてグラフＧ０１を作成するとともに、時刻Ｔ１〜Ｔｎにそれぞれ対応する４つの画素値に基づいてグラフＧ１ｎを作成する。

グラフＧ０１は、ポインタＰＴ０１が指向する画素値の初期相における時間的変化を表す。また、グラフＧ１ｎは、ポインタＰＴ１ｎが指向する画素値の遅延相における時間的変化を表す。また、作成されたグラフＧ０１およびＧ１ｎは、ＤＲＡＭ２２に描画される。

BI-RADS-MRIによれば、初期相に対応するグラフＧ０１の特性は“Slow”，“Medium”および“Fast”の３つに分類され、遅延相に対応するグラフＧ１ｎの特性は“Persistent”，“Plateau”および“Washout”の３つに分類される。

図９〜図１０においては、説明の便宜上、初期相において“Slow”の特性を表すグラフを“Ｇ０１ｓｌｗ”と定義し、初期相において“Medium”の特性を表すグラフを“Ｇ０１ｍｄ”と定義し、初期相において“Fast”の特性を表すグラフを“Ｇ０１ｆｓｔ”と定義する。また、遅延相において“Persistent”の特性を表すグラフを“Ｇ１ｎｐｓｔ”と定義し、遅延相において“Plateau”の特性を表すグラフを“Ｇ１ｎｐｌｔ”と定義し、遅延相において“Washout ”の特性を表すグラフを“Ｇ１ｎｗｓｈ”と定義する。

さらに、図１０では、初期相においてグラフＧ０１ｆｓｔが描画される場合を想定して、グラフＧ１ｎｐｓｔ，Ｇ１ｎｐｌｔおよびＧ１ｎｗｓｈをグラフＧ０１ｆｓｔに連続して描画している。このとき、グラフＧ１ｎｐｓｔ，Ｇ１ｎｐｌｔおよびＧ１ｎｗｓｈは、初期値からの画素値の増減率を示す。

ステップＳ３７では、ＤＲＡＭ２２に描画されたグラフＧ０１およびＧ１ｎの表示をディスプレイ装置２０に要求する。ディスプレイ装置２０は、グラフＧ０１およびＧ１ｎをＤＲＡＭ２２から読み出し、読み出されたグラフＧ０１およびＧ１ｎをステップＳ１５で表示された枠に合わせて画面ＳＣＲに表示する。グラフＧ０１およびＧ１ｎは、たとえば図８に示す要領で表示される。

ステップＳ３９ではキーボード１６および／またはマウス１８によってスライス位置変更操作が行われたか否かを判別し、ステップＳ４３ではマウス１８によってポインタ移動操作が行われたか否かを判別する。ステップＳ３９の判別結果がＹＥＳであれば、ステップＳ４１で変数Ｋの値を変更し、その後にステップＳ２９に戻る。また、ステップＳ４５の判別結果がＹＥＳであれば、ステップＳ４５でポインタＰＴ０１およびＰＴ１ｎの位置を共通的に変更し、その後にステップＳ３３に戻る。

したがって、補完画像ＩＭ０１〜ＩＭ１ｎおよびマラーマップＭＰ０１〜ＭＰ１ｎの表示態様はスライス位置変更操作に従って動的に変化し、ポインタＰＴ０１〜ＰＴ１ｎおよびグラフＧ０１〜Ｇ１ｎの表示態様はポインタ移動操作に従って動的に変化する。

以上の説明から分かるように、通信Ｉ／Ｆ１２は、時刻Ｔ０〜Ｔｎの各々でＭＲＩ装置３０が撮像した患部（乳房）のＭＲＩ画像を取得する。ＣＰＵ１４は、初期相および遅延相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表すカラーマップＭＰ０１およびＭＰ１ｎを通信Ｉ／Ｆ１２によって取得されたＭＲＩ画像に基づいて画面ＳＣＲに表示する(S5, S9, S21, S31)。ＣＰＵ１４はまた、ポインタ移動操作に従って断層上の所望の位置を共通的に指向するポインタＰＴ０１およびＰＴ１ｎを画面ＳＣＲに表示されたカラーマップＭＰ０１およびＭＰ１ｎにそれぞれ多重する(S33, S43, S45)。ＣＰＵ１４はさらに、初期相および遅延相にそれぞれ対応しかつポインタ移動操作に従う位置における画素値の時間的変化を各々が表すグラフＧ０１およびＧ１ｎを通信Ｉ／Ｆ１２によって取得されたＭＲＩ画像に基づいて画面ＳＣＲに表示する(S15, S35~S37)。

患部を診断するには初期相および遅延相における血流状態を把握する必要があるところ、カラーマップＭＰ０１〜ＭＰ１ｎおよびグラフＧ０１〜Ｇ１ｎはいずれも、初期相および遅延相にそれぞれ対応して患部の血流状態を表す。具体的には、カラーマップＭＰ０１〜ＭＰ１ｎは画素値に応じて異なる平面的色分布を表し、グラフＧ０１〜Ｇ１ｎは断層上の或る位置における画素値の時間的変化を表す。これによって、初期相および遅延相の各々における血流状態を複数の視点で把握することができる。

さらに、グラフＧ０１〜Ｇ１ｎを表示する際に注目される位置は、カラーマップＭＰ０１〜ＭＰ１ｎに共通的に多重表示されたポインタＰＴ０１〜ＰＴ１ｎによって指向され、かつポインタＰＴ０１〜ＰＴ１ｎの指向位置はポインタ移動操作に従って変更される。この結果、グラフＧ０１〜Ｇ１ｎの表示態様は、ポインタ移動操作に従って変化する。

ポインタＰＴ０１〜ＰＴ１ｎの指向位置およびグラフＧ０１〜Ｇ１ｎの表示態様の動的な変化は、カラーマップＭＰ０１〜ＭＰ１ｎおよびグラフＧ０１〜Ｇ１ｎが互いに異なる視点で血流状態を表すことと相俟って、患部を診断する際の利便性の向上に繋がる。

［第２実施例］
BI-RADS-MRIによれば、遅延相に対応するグラフＧ１ｎの特性は“Persistent”，“Plateau”および“Washout”に分類され、信号強度ないし画素値の変化の態様が注目されるのに対して、初期相に対応するグラフＧ０１の特性は“Slow”，“Medium”および“Fast”に分類され、信号強度ないし画素値の変化の速度が注目される。したがって、BI-RADS-MRIに沿った読影を初期相について行う上では、グラフＧ０１の傾きが重要になる。

他方、プレーン撮像の終了から造影剤の注入開始までには、実際には、準備や被検者の状態などにより十秒〜数十秒のタイムラグが発生する。にも関わらず、第１実施例では、造影剤の注入開始時刻がプレーン撮像時刻と一致することを前提としている。この結果、グラフＧ０１の傾きがこのタイムラグに起因して減少し、初期相における画素値の変化の速度を適切に評価することが難しくなる。

そこで、この実施例では、プレーン撮像の終了から造影剤の注入開始までのタイムラグを排除したグラフＧ０１の描画を可能とするべく、図１８に示すフロー図に従う処理を第１実施例に対して追加的に実行するようにしている。なお、図１８に示すステップＳ５１へは、図１３に示すステップＳ４３の判別結果がＮＯであるときに移行する。

ステップＳ５１では、キーボード１６および／またはマウス１８によってグラフ設定操作が行われたか否かを判別する。判別結果がＮＯであればステップＳ３９に戻り、判別結果がＹＥＳであればステップＳ５３に進む。

ステップＳ５３では、グラフ設定メニューを表すメニュー画像（以下、単に「グラフ設定メニュー」と呼ぶ。）をＤＲＡＭ２２に展開し、グラフ設定メニューの表示をディスプレイ装置２０に要求する。ディスプレイ装置２０は、グラフ設定メニューをＤＲＡＭ２２から読み出し、読み出されたグラフ設定メニューを画面ＳＣＲ上の診断画像に多重表示する。

図１４または図１６を参照して、グラフ設定メニューでは、プレーン撮像時刻（＝Ｔ０）を横軸の起点とするか、造影剤の注入開始時刻（＝Ｔ０´）を横軸の起点とするかを選択できる。前者を“項目１”とし、後者を“項目２”とすると、図１４は項目１にチェックマークが表示された状態を表し、図１６は項目２にチェックマークが表示された状態を示す。また、時刻Ｔ０´は時刻Ｔ１を基準として任意に設定可能であり、図１４および図１６では時刻Ｔ１の６０秒前が時刻Ｔ０´として設定されている。

図１８に戻って、ステップＳ５５ではチェック項目変更操作が行われたか否かを判別し、ステップＳ５７では秒数変更操作が行われたか否かを判別する。ステップＳ５５の判別結果がＹＥＳであれば、ステップＳ５９でチェックマークの表示対象を項目１および項目２の間で変更する。ステップＳ５７の判別結果がＹＥＳであれば、ステップＳ６１で変更後の秒数を項目２の所定の位置に表示する。

ステップＳ５９またはＳ６１の処理が完了するとステップＳ６３に進み、グラフ設定メニューの設定に従ってグラフの枠の横軸を調整する。したがって、グラフ設定メニューの項目２が選択された場合には、プレーン撮像時刻と造影剤の注入開始時刻との差分がステップＳ６３の処理によって補償される。

ステップＳ６５では、調整された横軸を有する枠に合わせてグラフＧ０１およびＧ１ｎを再描画する。項目１がチェックされた場合、画面ＳＣＲの上段には、図１５に示す要領でグラフの枠およびグラフＧ０１〜Ｇ１ｎが表示される。これに対して、項目２がチェックされた場合、画面ＳＣＲの上段には、図１７に示す要領でグラフの枠およびグラフＧ０１〜Ｇ１ｎが表示される。

この実施例では、造影剤の注入開始時刻（＝Ｔ０´）がプレーン撮像時刻（＝Ｔ０）よりも遅れる実情を想定している。したがって、図１５では、時刻Ｔ０´と時刻Ｔ０との誤差の分だけ、グラフＧ０１の傾きが減少する。これに対して、図１６では、造影剤の注入開始時刻（＝Ｔ０´）が横軸の起点となるため、グラフＧ０１の傾きは現実を反映する。

ステップＳ６５の処理が完了するか、またはステップＳ５５およびＳ５７の判別結果がいずれもＮＯであれば、キーボード１６および／またはマウス１８によって終了操作が行われたか否かをステップＳ６７で判別する。判別結果がＮＯであればステップＳ５５に戻り、判別結果がＹＥＳであればステップＳ３９に戻る。したがって、終了操作が行われない限り、グラフの枠およびグラフＧ０１〜Ｇ１ｎの表示態様は、チェック項目変更操作または秒数変更操作が行われる毎に動的に変化する。

この実施例によれば、グラフ設定メニューを操作することで、プレーン撮像時刻と造影剤の注入開始時刻との差分が補償された状態でグラフの枠およびグラフＧ０１〜Ｇ１ｎが画面ＳＣＲに表示される。これによって、BI-RADS-MRIに沿った読影が可能となる。

１０ …診断支援装置
１４ …ＣＰＵ
２０ …ディスプレイ装置
２４ …主記憶装置
３０ …ＭＲＩ装置
４０ …通信ネットワーク

Claims (9)

1. 複数の時刻の各々で断層撮像装置が撮像した患部の断層画像を取得する取得手段、
   複数の時相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表す複数のカラーマップを前記取得手段によって取得された断層画像に基づいて表示するカラーマップ表示手段、
   ポインタ移動操作に従って前記断層上の所望の位置を共通的に指向する複数のポインタを前記カラーマップ表示手段によって表示された複数のカラーマップにそれぞれ多重するポインタ多重手段、および
   前記複数の時相にそれぞれ対応しかつ前記ポインタ移動操作に従う位置における画素値の時間的変化を各々が表す複数のグラフを前記取得手段によって取得された断層画像に基づいて表示するグラフ表示手段を備える、診断支援装置。
2. 前記グラフの横軸および縦軸には時間および前記画素値がそれぞれ割り当てられ、
   前記画素値に応じて異なる色を示すカラーバーを前記グラフの縦軸に沿って表示するカラーバー表示手段をさらに備える、請求項１記載の診断画像処理装置。
3. 前記複数の時相をそれぞれ代表する複数の代表画像を前記取得手段によって取得された断層画像に基づいて表示する代表画像表示手段をさらに備え、
   前記カラーマップ表示手段は前記代表画像表示手段によって表示された複数の代表画像に前記複数のカラーマップをそれぞれ多重するカラーマップ多重手段を含む、請求項１または２記載の診断支援装置。
4. 前記複数の時相は初期相および遅延相を含み、
   前記複数のグラフは、前記初期相に対応しかつ初期値からの前記画素値の変化を表す第１グラフ、および前記遅延相に対応しかつ前記初期値からの前記画素値の変化を表す第２グラフを含む、請求項１ないし３のいずれかに記載の診断支援装置。
5. 前記断層撮像装置は前記患部への造影剤の注入を想定した装置であり、
   前記複数の時刻は前記患部に前記造影剤を注入する前の時刻を基準時刻として含み、
   前記グラフ表示手段は、前記基準時刻と前記造影剤の注入を開始した時刻との差分を補償する補償手段、および前記補償手段による補償を踏まえて前記第１グラフを作成する第１グラフ作成手段を含む、請求項４記載の診断支援装置。
6. 前記取得手段が前記複数の時刻の各々に対応して取得する断層画像は複数のスライス位置の各々の断層を表し、
   前記カラーマップ表示手段および前記グラフ表示手段は前記複数のスライス位置のうちの共通のスライス位置を対象として処理を実行する、請求項１ないし５のいずれかに記載の診断支援装置。
7. 前記カラーマップ表示手段および前記グラフ表示手段が対象とするスライス位置をスライス位置変更操作に従って変更するスライス位置変更手段をさらに備える、請求項６記載の診断支援装置。
8. 複数の時刻の各々で断層撮像装置が撮像した患部の断層画像を取得する取得手段を備える診断支援装置によって実行される診断支援方法であって、
   複数の時相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表す複数のカラーマップを前記取得手段によって取得された断層画像に基づいて表示するカラーマップ表示ステップ、
   ポインタ移動操作に従って前記断層上の所望の位置を共通的に指向する複数のポインタを前記カラーマップ表示ステップによって表示された複数のカラーマップにそれぞれ多重するポインタ多重ステップ、および
   前記複数の時相にそれぞれ対応しかつ前記ポインタ移動操作に従う位置における画素値の時間的変化を各々が表す複数のグラフを前記取得手段によって取得された断層画像に基づいて表示するグラフ表示ステップを備える、診断支援方法。
9. 複数の時刻の各々で断層撮像装置が撮像した患部の断層画像を取得する取得手段を備える診断支援装置のプロセッサに、
   複数の時相にそれぞれ対応しかつ画素値に応じて異なる平面的色分布を各々が表す複数のカラーマップを前記取得手段によって取得された断層画像に基づいて表示するカラーマップ表示ステップ、
   ポインタ移動操作に従って前記断層上の所望の位置を共通的に指向する複数のポインタを前記カラーマップ表示ステップによって表示された複数のカラーマップにそれぞれ多重するポインタ多重ステップ、および
   前記複数の時相にそれぞれ対応しかつ前記ポインタ移動操作に従う位置における画素値の時間的変化を各々が表す複数のグラフを前記取得手段によって取得された断層画像に基づいて表示するグラフ表示ステップを実行させるための、診断支援プログラム。
   
   

Images