(57)【要約】      （修正有）【課題】検査部位の血流特性を評価する診断者間での評価のばらつきを抑制する画像診断支援装置を提供する。【解決手段】画像診断支援装置１０は、既定条件を満足する画素に対応する有色画素からなるカラーマップＭＰｔ０１，ＭＰｔ１２を、時刻Ｔ０〜Ｔ２にそれぞれ撮像されかつ共通の検査部位を表す断層画像ＩＭｔ０〜ＩＭｔ２に基づいてビューワ２０に表示する。表示されたカラーマップＭＰｔ０１，ＭＰｔ１２上の所望の有色画素を指定するクリック操作を受け付けると、指定された有色画素から連続する有色画素をカラーマップＭＰｔ０１，ＭＰｔ１２上で探索し、指定された有色画素および探知された有色画素の各々に対応する微小検査部位の血流特性を断層画像ＩＭｔ０〜ＩＭｔ２に基づいて解析する。その後、解析結果が特定の血流特性を示す微小検査部位が解析対象に占める割合を含む解析情報をビューワ２０に追加表示する。【選択図】図１
2017136317.tif 000002
【特許請求の範囲】
【請求項１】
既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップを複数の時刻にそれぞれ撮像されかつ共通の検査部位を表す複数の断層画像に基づいて表示する第１表示手段、
前記第１表示手段によって表示されたカラーマップ上の所望の有色画素を指定する画素指定操作を受け付ける受け付け手段、
前記画素指定操作によって指定された有色画素から連続する１または２以上の有色画素を前記第１表示手段によって表示されたカラーマップ上で探索する探索手段、
前記画素指定操作によって指定された有色画素および前記探索手段によって探知された有色画素の各々に対応する微小検査部位の血流特性を前記複数の断層画像に基づいて解析する解析手段、および
前記解析手段の解析結果が特定の血流特性を示す微小検査部位が前記解析手段の解析対象に占める割合を含む解析情報を表示する第２表示手段を備える、画像診断支援装置。
【請求項２】
前記解析手段の解析対象とされた微小検査部位に対応する有色画素の数が閾値に達した時点で前記解析手段の処理を制限する制限手段をさらに備える、請求項１記載の画像診断支援装置。
【請求項３】
前記第２表示手段によって表示される解析情報は前記解析手段の解析対象とされた微小検査部位に対応する有色画素によって形成された限定カラーマップをさらに含む、請求項１または２記載の画像診断支援装置。
【請求項４】
前記第２表示手段によって表示される解析情報は前記解析手段の解析結果が特定の血流特性と異なる血流特性を示す微小検査部位が前記解析手段の解析対象に占める割合をさらに含む、請求項１ないし３のいずれかに記載の画像診断支援装置。
【請求項５】
前記複数の時刻は前記検査部位に対する造影剤の注入前後の時刻であり、
前記特定の血流特性は早期相においてＦａｓｔを示しかつ遅延相においてＷａｓｈｏｕｔを示す特性である、請求項１ないし４のいずれかに記載の画像診断支援装置。
【請求項６】
前記カラーマップを構成する１または２以上の有色画素の各々の色は対応する画素が有する画素値に基づいて決定される、請求項１ないし５のいずれかに記載の画像診断支援装置。
【請求項７】
前記第２表示手段によって表示された解析情報を保存操作に応答して保存する保存手段をさらに備える、請求項１ないし６のいずれかに記載の画像診断支援装置。
【請求項８】
画像診断支援装置のプロセッサによって実行される画像診断支援方法であって、
既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップを複数の時刻にそれぞれ撮像されかつ共通の検査部位を表す複数の断層画像に基づいて表示する第１表示ステップ、
前記第１表示ステップによって表示されたカラーマップ上の所望の有色画素を指定する画素指定操作を受け付ける受け付けステップ、
前記画素指定操作によって指定された有色画素から連続する１または２以上の有色画素を前記第１表示ステップによって表示されたカラーマップ上で探索する探索ステップ、
前記画素指定操作によって指定された有色画素および前記探索ステップによって探知された有色画素の各々に対応する微小検査部位の血流特性を前記複数の断層画像に基づいて解析する解析ステップ、および
前記解析ステップの解析結果が特定の血流特性を示す微小検査部位が前記解析ステップの解析対象に占める割合を含む解析情報を表示する第２表示ステップを備える、画像診断支援方法。
【請求項９】
画像診断支援装置のプロセッサに、
既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップを複数の時刻にそれぞれ撮像されかつ共通の検査部位を表す複数の断層画像に基づいて表示する第１表示ステップ、
前記第１表示ステップによって表示されたカラーマップ上の所望の有色画素を指定する画素指定操作を受け付ける受け付けステップ、
前記画素指定操作によって指定された有色画素から連続する１または２以上の有色画素を前記第１表示ステップによって表示されたカラーマップ上で探索する探索ステップ、
前記画素指定操作によって指定された有色画素および前記探索ステップによって探知された有色画素の各々に対応する微小検査部位の血流特性を前記複数の断層画像に基づいて解析する解析ステップ、および
前記解析ステップの解析結果が特定の血流特性を示す微小検査部位が前記解析ステップの解析対象に占める割合を含む解析情報を表示する第２表示ステップを実行させるための、画像診断支援プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、画像診断支援装置、画像診断支援方法および画像診断支援プログラムに関し、特に、モダリティによって撮像された断層画像を参照した画像診断医の診断を支援する、画像診断支援装置、画像診断支援方法および画像診断支援プログラムに関する。
【背景技術】
【０００２】
米国で生まれた乳房ＭＲＩ検査のガイドライン「BI-RADS-MRI」が日本でも認知され、現在は日本版BI-RADS-MRIに則った検査方法，読影方法，レポート記述方法が多くの施設で採用されるようになってきた。
【０００３】
BI-RADS-MRIでは、“ＦＯＣＵＳ”と呼ばれる小造影効果が定義されている。この小造影効果が表れた病変部は、正常乳腺の一部や、線維腺腫，乳房内リンパ節など小さな良性病変を想定しており、ほとんどの場合においてカテゴリー３以下と診断される。ただし、稀に小さな浸潤癌などが含まれていることがあるため、辺縁が不整な場合や血流評価でFast/Washout を示す場合は、さらなる精査を考慮する必要がある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−１９６９９０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、小造影効果が表れた病変部は微小であるため、形態評価が難しく、Time Intensity Curveによる精度も低くなってしまう。この結果、診断者によって評価が変わる可能性が高いという問題がある。
【０００６】
これに関連して、特許文献１は、血管内に造影剤が注入された被験者を第１の時間帯に撮像して得られたフレーム群から第１の画像を生成し、当該被験者を第２の時間帯に撮像して得られたフレーム群から第２の画像を生成し、第１の画像および第２の画像において重複している血管領域内にＲＯＩを設定することを開示している。
【０００７】
より詳しくは、ＲＯＩは、実際の大動脈の大きさを超えない半径４ピクセルの領域に含まれるカウント値の合計が最も高くなる領域に設定される。なお、カウント値は、造影剤濃度と対応関係を有して各画素に割り当てられた値である。
【０００８】
しかし、特許文献１は、ＲＯＩの自動設定について開示するものの、設定されたＲＯＩの特性を如何にして定量的に評価については開示しておらず、診断者によって評価が変わる可能性がある。
【０００９】
それゆえに、この発明の主たる目的は、検査部位の血流特性を評価する診断者間での評価のばらつきを抑制することができる、画像診断支援装置、画像診断支援方法および画像診断支援プログラムを提供することである。
【課題を解決するための手段】
【００１０】
この発明に従う画像診断支援装置(10：実施例で相当する参照符号。以下同じ)は、既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップ(MPt01, MPt12)を複数の時刻にそれぞれ撮像されかつ共通の検査部位を表す複数の断層画像(IMt0~ IMt2)に基づいて表示する第１表示手段(S27~S33)、第１表示手段によって表示されたカラーマップ上の所望の有色画素を指定する画素指定操作を受け付ける受け付け手段(S15, S17)、画素指定操作によって指定された有色画素から連続する１または２以上の有色画素を第１表示手段によって表示されたカラーマップ上で探索する探索手段(S41~S45, S49~S51, S57~S63, S71~S77)、画素指定操作によって指定された有色画素および探索手段によって探知された有色画素の各々に対応する微小検査部位の血流特性を複数の断層画像に基づいて解析する解析手段(S47, S67)、および解析手段の解析結果が特定の血流特性を示す微小検査部位が解析手段の解析対象に占める割合を含む解析情報を表示する第２表示手段(S79~S121)を備える。
【００１１】
好ましくは、解析手段の解析対象とされた微小検査部位に対応する有色画素の数が閾値に達した時点で解析手段の処理を制限する制限手段(S55, S69)がさらに備えられる。
【００１２】
好ましくは、第２表示手段によって表示される解析情報は解析手段の解析対象とされた微小検査部位に対応する有色画素によって形成された限定カラーマップ(RMt01, RMt12)をさらに含む。
【００１３】
好ましくは、第２表示手段によって表示される解析情報は解析手段の解析結果が特定の血流特性と異なる血流特性を示す微小検査部位が解析手段の解析対象に占める割合をさらに含む。
【００１４】
好ましくは、複数の時刻は検査部位に対する造影剤の注入前後の時刻であり、特定の血流特性は早期相においてＦａｓｔを示しかつ遅延相においてＷａｓｈｏｕｔを示す特性である。
【００１５】
好ましくは、カラーマップを構成する１または２以上の有色画素の各々の色は対応する画素が有する画素値に基づいて決定される。
【００１６】
好ましくは、第２表示手段によって表示された解析情報を保存操作に応答して保存する保存手段(S123, S125)がさらに備えられる。
【００１７】
この発明に従う画像診断支援方法は、画像診断支援装置(10)のプロセッサ(14)によって実行される画像診断支援方法であって、既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップ(MPt01, MPt12)を複数の時刻にそれぞれ撮像されかつ共通の検査部位を表す複数の断層画像(IMt01~IMt2)に基づいて表示する第１表示ステップ(S27~S33)、第１表示ステップによって表示されたカラーマップ上の所望の有色画素を指定する画素指定操作を受け付ける受け付けステップ(S15, S17)、画素指定操作によって指定された有色画素から連続する１または２以上の有色画素を第１表示ステップによって表示されたカラーマップ上で探索する探索ステップ(S41~S45, S49~S51, S57~S63, S71~S77)、画素指定操作によって指定された有色画素および探索ステップによって探知された有色画素の各々に対応する微小検査部位の血流特性を複数の断層画像に基づいて解析する解析ステップ(S47, S67)、および解析ステップの解析結果が特定の血流特性を示す微小検査部位が解析ステップの解析対象に占める割合を含む解析情報を表示する第２表示ステップ(S79~S121)を備える。
【００１８】
この発明に従う画像診断支援プログラムは、画像診断支援装置(10)のプロセッサ(14)に、既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップ(MPt01, MPt12)を複数の時刻にそれぞれ撮像されかつ共通の検査部位を表す複数の断層画像(IMt01~IMt2)に基づいて表示する第１表示ステップ(S27~S33)、第１表示ステップによって表示されたカラーマップ上の所望の有色画素を指定する画素指定操作を受け付ける受け付けステップ(S15, S17)、画素指定操作によって指定された有色画素から連続する１または２以上の有色画素を第１表示ステップによって表示されたカラーマップ上で探索する探索ステップ(S41~S45, S49~S51, S57~S63, S71~S77)、画素指定操作によって指定された有色画素および探索ステップによって探知された有色画素の各々に対応する微小検査部位の血流特性を複数の断層画像に基づいて解析する解析ステップ(S47, S67)、および解析ステップの解析結果が特定の血流特性を示す微小検査部位が解析ステップの解析対象に占める割合を含む解析情報を表示する第２表示ステップ(S79~S121)を実行させるための、画像診断支援プログラムである。
【発明の効果】
【００１９】
探索手段は画素指定操作によって指定された有色画素に連続する１または２以上の有色画素を探索するところ、これらの有色画素の分布態様（つまり、指定された有色画素および探知された有色画素の総数ないし面積）は、検査部位の状態によって異なる。また、診断者によって関心があるのは、特定の血流特性を示す微小検査部位である。
【００２０】
これを踏まえて、指定された有色画素および探知された有色画素の各々に対応する微小検査部位の血流特性を解析し、解析結果が特定の血流特性を示す微小検査部位が解析対象に占める割合を表示するようにしている。これによって、解析結果が特定の血流特性を示す微小検査部位の大きさが正規化され、検査部位の血流特性を評価する診断者間での評価のばらつきを抑制することができる。
【００２１】
この発明の上述の目的，その他の目的，特徴および利点は、図面を参照して行う以下の実施例の詳細な説明から一層明らかとなろう。
【図面の簡単な説明】
【００２２】
【図１】この実施例に適用される医療診断支援システムの構成の一例を示すブロック図である。
【図２】（Ａ）は時刻Ｔ０にダイナミック造影ＭＲＩシーケンスで撮像された複数のダイナミック造影ＭＲＩ画像をそれぞれ収めた複数のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｂ）は時刻Ｔ１にダイナミック造影ＭＲＩシーケンスで撮像された複数のダイナミック造影ＭＲＩ画像をそれぞれ収めた複数のＤＩＣＯＭファイルの構造の一例を示す図解図であり、（Ｃ）は時刻Ｔ２にダイナミック造影ＭＲＩシーケンスで撮像された複数のダイナミック造影ＭＲＩ画像をそれぞれ収めた複数のＤＩＣＯＭファイルの構造の一例を示す図解図である。
【図３】この実施例の画像診断支援装置を構成するＤＲＡＭに設けられたビューワエリアおよびワークエリアのマッピング状態の一例を示す図解図である。
【図４】ダイナミック造影ＭＲＩ画像およびカラーマップがビューワに表示された状態の一例を示す図解図である。
【図５】カーソルがビューワに多重表示された状態の一例を示す図解図である。
【図６】（Ａ）はカラーマップ上の有色画素を解析するときの動作の一部を示す図解図であり、（Ｂ）はカラーマップ上の有色画素を解析するときの動作の他の一部を示す図解図であり、（Ｃ）はカラーマップ上の有色画素を解析するときの動作のその他の一部を示す図解図であり、（Ｄ）はカラーマップ上の有色画素を解析するときの動作のさらにその他の一部を示す図解図であり、（Ｅ）はカラーマップ上の有色画素を解析するときの動作の他の一部を示す図解図であり、（Ｆ）はカラーマップ上の有色画素を解析するときの動作のその他の一部を示す図解図である。
【図７】（Ａ）はカラーマップ上の有色画素を解析するときの動作のさらにその他の一部を示す図解図であり、（Ｂ）はカラーマップ上の有色画素を解析するときの動作の他の一部を示す図解図であり、（Ｃ）はカラーマップ上の有色画素を解析するときの動作のその他の一部を示す図解図であり、（Ｄ）はカラーマップ上の有色画素を解析するときの動作のさらにその他の一部を示す図解図であり、（Ｅ）はカラーマップ上の有色画素を解析するときの動作の他の一部を示す図解図である。
【図８】（Ａ）は早期相に対応する限定カラーマップの一例を示す図解図であり、（Ｂ）は遅延相に対応する限定カラーマップの一例を示す図解図である。
【図９】血流特性マトリクスの一例を示す図解図である。
【図１０】解析情報の一例を示す図解図である。
【図１１】解析情報がビューワに追加表示された状態の一例を示す図解図である。
【図１２】この実施例の画像診断支援装置を構成するＣＰＵの動作の一部を示すフロー図である。
【図１３】この実施例の画像診断支援装置を構成するＣＰＵの動作の他の一部を示すフロー図である。
【図１４】この実施例の画像診断支援装置を構成するＣＰＵの動作のその他の一部を示すフロー図である。
【図１５】この実施例の画像診断支援装置を構成するＣＰＵの動作のさらにその他の一部を示すフロー図である。
【図１６】この実施例の画像診断支援装置を構成するＣＰＵの動作の他の一部を示すフロー図である。
【図１７】この実施例の画像診断支援装置を構成するＣＰＵの動作のその他の一部を示すフロー図である。
【図１８】この実施例の画像診断支援装置を構成するＣＰＵの動作のさらにその他の一部を示すフロー図である。
【発明を実施するための形態】
【００２３】
図１を参照して、この実施例の医療診断支援システムは、院内ＬＡＮ５０を介して互いに接続された画像診断支援装置１０，ＭＲＩ装置（モダリティ）３０および画像サーバ４０によって構成される。画像診断支援装置１０において、バスＢＳ１には、通信Ｉ／Ｆ１２，ＣＰＵ１４，キーボード／マウス１６，ＨＤＤ１８，ビューワ２０およびＤＲＡＭ２２が接続される。
【００２４】
ＭＲＩ装置３０は、核磁気共鳴現象を利用して、検査部位である被検者の乳房の断層画像を作成する装置である。ダイナミック造影ＭＲＩシーケンスに従う乳房の撮像処理は、体内に造影剤を注入する前の時刻Ｔ０に１回実行され（これを“プレーン撮像”と呼ぶ）、造影剤の注入が開始された後の互いに異なる時刻Ｔ１およびＴ２に１回ずつ実行される。なお、時刻Ｔ１は早期相に対応する時刻であり、時刻Ｔ２は遅延相に対応する時刻である。
【００２５】
時刻Ｔ０〜Ｔ２の各々での撮像処理の結果、“水平４００画素×垂直４００画素”のマトリクスサイズと“０．８ｍｍ”の画素サイズとを有するダイナミック造影ＭＲＩ画像が、０．８ｍｍ間隔で配された１２５個の水平スライス位置の各々に対応して作成される。作成された各々のダイナミック造影ＭＲＩ画像は、ファイル情報（ファイルサイズ，マトリクスサイズ，画素サイズ，スライス間隔），患者情報，撮像時刻，スライス位置とともにＤＩＣＯＭファイルに収められる。
【００２６】
この結果、時刻Ｔ０では図２（Ａ）に示す１２５個のＤＩＣＯＭファイルが作成され、時刻Ｔ１では図２（Ｂ）に示す１２５個のＤＩＣＯＭファイルが作成され、時刻Ｔ２では図２（Ｃ）に示す１２５個のＤＩＣＯＭファイルが作成される。作成されたＤＩＣＯＭファイルはいずれも、院内ＬＡＮ５０を経て画像サーバ４０に保存される。
【００２７】
画像診断支援装置１０に設けられた通信Ｉ／Ｆ１２は、こうして画像サーバ４０に蓄積された合計３７５個のＤＩＣＯＭファイルを院内ＬＡＮ５０を通して取得する。取得されたＤＩＣＯＭファイルは、バスＢＳ１を介してＨＤＤ１８に保存される。
【００２８】
ＣＰＵ１４は、画像診断医による画像診断を支援するべく、ＨＤＤ１８に保存されたＤＩＣＯＭファイルを対象とする処理を図１２〜図１８に示すフロー図に従って実行する。なお、このフロー図に相当するプログラムは、画像診断支援プログラムとしてＨＤＤ１８に記憶される。
【００２９】
図１２を参照して、ステップＳ１ではダイナミック造影ＭＲＩ画像のスライス位置を識別する変数Ｋを初期値に設定し、ステップＳ３では検査画像をビューワ２０に表示する。ステップＳ３の処理は、具体的には図１３に示すサブルーチンに従う。
【００３０】
ステップＳ２１では、時刻Ｔ０に撮像されたＫ番目のスライス位置のダイナミック造影ＭＲＩ画像を“ＩＭｔ０”として定義し、ダイナミック造影ＭＲＩ画像ＩＭｔ０をＤＲＡＭ２２に展開する。図３を参照して、ＤＲＡＭ２２には６つの小エリアを有するビューワエリア２２ｖが設けられる。ダイナミック造影ＭＲＩ画像ＩＭｔ０は、ビューワエリア２２ｖの左上段に配された小エリアに展開される。
【００３１】
ステップＳ２３では、時刻Ｔ１に撮像されたＫ番目のスライス位置のダイナミック造影ＭＲＩ画像を“ＩＭｔ１”として定義し、ダイナミック造影ＭＲＩ画像ＩＭｔ１をビューワエリア２２ｖの中央上段に配された小エリアに展開する。
【００３２】
ステップＳ２５では、時刻Ｔ２に撮像されたＫ番目のスライス位置のダイナミック造影ＭＲＩ画像を“ＩＭｔ２”として定義し、ダイナミック造影ＭＲＩ画像ＩＭｔ２をビューワエリア２２ｖの右上段に配された小エリアに展開する。
【００３３】
こうして展開されたダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２は、ビューワ２０によって読み出される。読み出されたダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２は、図４に示す要領で画面の左上段，中央上段および右上段にそれぞれ表示される。
【００３４】
ステップＳ２７では、ビューワエリア２２ｖ上のダイナミック造影ＭＲＩ画像ＩＭｔ０およびＩＭｔ１に基づいてカラーマップＭＰｔ０１を作成する。ステップＳ２９では、ビューワエリア２２ｖ上のダイナミック造影ＭＲＩ画像ＩＭｔ１およびＩＭｔ２に基づいてカラーマップＭＰｔ１２を作成する。作成されたカラーマップＭＰｔ０１およびＭＰｔ１２は、ＤＲＡＭ２２に設けられたワークエリア２２ｗに保持される。
【００３５】
ここで、カラーマップＭＰｔ０１を構成する色（＝端的には色相）は、ダイナミック造影ＭＲＩ画像ＩＭｔ０上の画素値とダイナミック造影ＭＲＩ画像ＩＭｔ１上の対応する画素値とを比較して決定される。また、カラーマップＭＰｔ１２を構成する色は、ダイナミック造影ＭＲＩ画像ＩＭｔ１上の画素値とダイナミック造影ＭＲＩ画像ＩＭｔ２上の対応する画素値とを比較して決定される。
【００３６】
また、カラーマップＭＰｔ０１およびＭＰｔ１２のいずれについても、各画素に色を付けるか否かは、ダイナミック造影ＭＲＩ画像ＩＭｔ０の画素値を比較対象とした、ダイナミック造影ＭＲＩ画像ＩＭｔ１の画素値の増加率および増加量に基づいて判断する。具体的には、増加率が４０％以上でかつ増加量が閾値（＝増加率が３９％以下の画素値のうち最大画素値）を上回るという条件が課せられ、この条件を満足する画素に対して色が付される。
【００３７】
したがって、付される色はカラーマップＭＰｔ０１およびＭＰｔ１２の間で相違するものの、色の分布状態はカラーマップＭＰｔ０１およびＭＰｔ１２の間で一致する。
【００３８】
参考までに、有色画素（＝色が付された画素）は腫瘍の染まり具合の形状を表すところ、この形状が円形，卵型，分葉状，不整形のいずれであるかによって良悪性の診断結果が変わる。また、早期相特性（早期相における血流特性）が“Fast”に該当しかつ遅延相特性（遅延相における血流特性）が“Washout”に該当する画素に対応する微小検査部位が最も悪性であると言われている。
【００３９】
ステップＳ３１では、ワークエリア２２ｗ上のカラーマップＭＰｔ０１をビューワエリア２２ｖの左下段に配された小エリアに展開する。また、ステップＳ３３では、ワークエリア２２ｗ上のカラーマップＭＰｔ１２をビューワエリア２２ｖの中央下段に配された小エリアに展開する。
【００４０】
なお、ステップＳ３１では、たとえばダイナミック造影ＭＲＩ画像ＩＭｔ０が、ビューワエリア２２ｖの左下段の小エリアに付随的に展開される。同様に、ステップＳ３３でも、たとえばダイナミック造影ＭＲＩ画像ＩＭｔ０が、ビューワエリア２２ｖの中央下段の小エリアに付随的に展開される。カラーマップＭＰｔ０１およびＭＰｔ１２は、こうして展開されたダイナミック造影ＭＲＩ画像ＩＭｔ０に多重される。
【００４１】
ビューワエリア２２ｖに展開されたカラーマップＭＰｔ０１およびＭＰｔ１２（いずれにもダイナミック造影ＭＲＩ画像ＩＭｔ０が付随）は、ビューワ２０によって読み出され、画面の左下段および中央下段にそれぞれ表示される（図４参照）。
【００４２】
このように、ビューワ２０の画面は６つの子画面に分割され、ダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２，カラーマップＭＰｔ０１，ＭＰｔ１２はこのうちの５つの子画面にそれぞれ表示される。検査画像は、こうして表示されたダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２，カラーマップＭＰｔ０１，ＭＰｔ１２によって構成される。ステップＳ３３の処理が完了すると、上階層のルーチンに復帰する。
【００４３】
図１２に戻って、ステップＳ５では、カーソル多重表示命令をビューワ２０に与える。この命令は、互いに同期する２つのカーソルを画面上のカラーマップＭＰｔ０１，ＭＰｔ１２にそれぞれ多重表示するための命令である。ビューワ２０は、与えられた命令に従って、２つのカーソルをカラーマップＭＰｔ０１，ＭＰｔ１２にそれぞれ多重表示する（図５参照）。
【００４４】
ステップＳ７ではキーボード／マウス１６によってカーソル移動操作が行われたか否かを判別し、ステップＳ９ではキーボード／マウス１６によってページング操作が行われたか否かを判別する。
【００４５】
ステップＳ７の判別結果およびステップＳ９の判別結果のいずれもがＮＯであればステップＳ７に戻り、ステップＳ９の判別結果がＹＥＳであればステップＳ１１で変数Ｋを更新してからステップＳ３に戻り、ステップＳ７の判別結果がＹＥＳであればステップＳ１３に進む。
【００４６】
ステップＳ１３では、カーソル移動命令をビューワ２０に与える。この結果、カラーマップＭＰｔ０１，ＭＰｔ１２に多重表示された２つのカーソルが同期的に移動する。ステップＳ１５では、キーボード／マウス１６によってクリック操作が行われたか否かを判別する。また、ステップＳ１７では、カラーマップＭＰｔ０１またはＭＰｔ１２に多重表示されたカーソルが有色画素を指定しているか否かを判別する。
【００４７】
ステップＳ１５の判別結果およびステップＳ１７の判別結果のいずれか一方がＮＯであれば、そのままステップＳ７に戻る。これに対して、ステップＳ１５の判別結果およびステップＳ１７の判別結果の両方がＹＥＳであれば、ステップＳ１９で解析＆表示処理を実行してからステップＳ７に戻る。
【００４８】
ステップＳ１９の解析＆表示処理は、図１４〜図１８に示すサブルーチンに従って実行される。
【００４９】
図１４に示すステップＳ４１では、カーソルが指定する有色画素を中心とする水平３画素×垂直３画素の領域を基準領域Ｒ１として設定する（図６（Ａ），図６（Ｂ）参照）。ステップＳ４３では、基準領域Ｒ１に属する９画素のうちラスタ方向における先頭画素を対象画素ＰＩＸｒｅｆとして指定する。
【００５０】
ステップＳ４５では対象画素ＰＩＸｒｅｆが有色画素であるか否かを判別し、判別結果がＮＯであればそのままステップＳ４９に進む。これに対して判別結果がＹＥＳであれば、ステップＳ４７で対象画素ＰＩＸｒｅｆに対応する微小検査部位の早期相特性および遅延相特性を解析し、その後にステップＳ４９に進む。
【００５１】
詳しくは、ステップＳ４７の処理は、早期相特性が“Fast”，“Medium”および“Slow”のいずれに該当し、遅延相特性が “Washout”，“Plateau”および“Persistent”のいずれに該当するかを解析する処理である。早期相特性は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化に注目して解析される。また、遅延相特性は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化に注目して解析される。
【００５２】
なお、早期相特性が“Fast”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化率が１００％以上の画素である。早期相特性が“Medium”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化率が５０％〜９９％の画素である。早期相特性が“Slow”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化率が４９％以下の画素である。
【００５３】
また、遅延相特性が“Washout”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化率が−１０％以下の画素である。遅延相特性が“Plateau”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化率が−１０％〜１０％の画素である。遅延相特性が“Persistent”に該当する画素は、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化率が１０％以上の画素である。
【００５４】
ステップＳ４９では、対象画素ＰＩＸｒｅｆが基準領域Ｒ１に属する９画素のうちラスタ方向における末尾画素であるか否かを判別する。判別結果がＮＯであればステップＳ５１に進み、対象画素ＰＩＸｒｅｆを基準領域Ｒ１内でラスタ方向における次画素に更新する。更新が完了すると、ステップＳ４５に戻る。
【００５５】
ステップＳ４９の判別結果がＹＥＳであれば、基準領域Ｒ１を囲む１６個の画素の少なくとも１つが有色画素であるか否かをステップＳ５３で判別し、ステップＳ４７で解析処理を施された画素の総数が閾値ＴＨ１（＝基準領域Ｒ１を形成する画素の総数以下の値であり、たとえば“８”）以下であるか否かをステップＳ５５で判別する。
【００５６】
ステップＳ５３の判別結果またはステップＳ５５の判別結果がＮＯであれば、ステップＳ７９に進む。これに対して、ステップＳ５３の判別結果およびステップＳ５５の判別結果のいずれもがＹＥＳであれば、ステップＳ５７に進む。
【００５７】
ステップＳ５７では基準領域Ｒ１内の先頭画素を対象画素ＰＩＸｒｅｆとして設定し、ステップＳ５９では対象画素ＰＩＸｒｅｆを中心とする水平３画素×垂直３画素の領域を探索領域Ｒ２として設定する（図６（Ｃ）参照）。ステップＳ６１では、こうして設定された探索領域Ｒ２に属する９画素のうちラスタ方向における先頭画素を対象画素ＰＩＸｓｃｈとして指定する。
【００５８】
ステップＳ６３では対象画素ＰＩＸｓｃｈが有色画素であるか否かを判別し、ステップＳ６５では対象画素ＰＩＸｓｃｈが未解析の画素（ステップＳ４７の解析処理を施されていない画素）であるか否かを判別する。ステップＳ６３の判別結果またはステップＳ６５の判別結果がＮＯであればそのままステップＳ６９に進む一方、ステップＳ６３の判別結果およびステップＳ６５の判別結果のいずれもがＹＥＳであればステップＳ６７の処理を経てステップＳ６９に進む。
【００５９】
ステップＳ６７では、ステップＳ４７と同じ要領で、対象画素ＰＩＸｓｃｈに対応する微小検査部位の早期相特性および遅延相特性を解析する。つまり、ダイナミック造影ＭＲＩ画像ＩＭｔ０・ＩＭｔ１間での画素値の変化に注目して対象画素ＰＩＸｓｃｈに対応する微小検査部位の早期相における血流特性を解析し、ダイナミック造影ＭＲＩ画像ＩＭｔ１・ＩＭｔ２間での画素値の変化に注目して対象画素ＰＩＸｓｃｈに対応する微小検査部位の遅延相における血流特性を解析する。
【００６０】
ステップＳ６９ではステップＳ４７およびステップＳ６７で解析処理を施された画素の総数が閾値ＴＨ１以下であるか否かを判別し、判別結果がＮＯであればステップＳ７９に進む一方、判別結果がＹＥＳであればステップＳ７１に進む。ステップＳ７１では、対象画素ＰＩＸｓｃｈが探索領域Ｒ２に属する９画素のうちラスタ方向における末尾画素であるか否かを判別する。ステップＳ７３では、対象画素ＰＩＸｒｅｆが基準領域Ｒ１に属する９画素のうちラスタ方向における末尾画素であるか否かを判別する。
【００６１】
ステップＳ７１の判別結果がＮＯであればステップＳ７３に進み、対象画素ＰＩＸｓｃｈを探索領域Ｒ２内のラスタ方向における次画素に更新する。更新が完了すると、ステップＳ６３に戻る。ステップＳ７１の判別結果がＹＥＳでかつステップＳ７５の判別結果がＮＯであればステップＳ７７に進み、対象画素ＰＩＸｒｅｆを基準領域Ｒ１内のラスタ方向における次画素に更新する。更新が完了すると、ステップＳ５９に戻る。
【００６２】
したがって、探索領域Ｒ２の配置は、図６（Ｂ）〜図６（Ｆ），図７（Ａ）〜図７（Ｅ）の順に更新される。この結果、クリック操作によって指定された有色画素から連続する１または２以上の有色画素がカラーマップＭＰｔ０１上で探索され、探知された有色画素に対応する微小検査部位の血流特性がダイナミック造影ＭＲＩ画像ＩＭｔ０〜ＩＭｔ２に基づいて解析される。
【００６３】
ステップＳ７１の判別結果およびステップＳ７５の判別結果のいずれもがＹＥＳであれば、ステップＳ７９に進み、ステップＳ４７およびＳ６７で解析処理を施された複数の有色画素の塊を解析領域として定義する。ステップＳ８１では、定義された解析領域に外接する矩形の輪郭を有する一部のカラーマップをワークエリア２２ｗ上のカラーマップＭＰ０１ｔおよびＭＰ１２ｔから切り出す。
【００６４】
カラーマップＭＰ０１ｔから切り出された一部のカラーマップを“限定カラーマップＲＭ０１ｔ”と定義し、カラーマップＭＰ１２ｔから切り出された一部のカラーマップを“限定カラーマップＲＭ１２ｔ”と定義すると、限定カラーマップＲＭ０１ｔおよびＲＭ１２ｔはそれぞれ図８（Ａ）および図８（Ｂ）に示すように描かれ、そのいずれも、ワークエリア２２ｗに保持される。なお、解析領域は、図８（Ａ）および図８（Ｂ）の各々に示す太線で囲まれた領域に相当する。
【００６５】
ステップＳ８３ではFast/Washout画素（＝早期相特性が“Fast”に該当しかつ遅延相特性が“Washout”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ８５ではFast/Washout画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００６６】
ステップＳ８７ではFast/Plateau画素（＝早期相特性が“Fast”に該当しかつ遅延相特性が“Plateau”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ８９ではFast/Plateau画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００６７】
ステップＳ９１ではFast/Persistent画素（＝早期相特性が“Fast”に該当しかつ遅延相特性が“Persistent”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ９３ではFast/Persistent画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００６８】
ステップＳ９５ではMedium/Washout画素（＝早期相特性が“Medium”に該当しかつ遅延相特性が“Washout”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ９７ではMedium/Washout画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００６９】
ステップＳ９９ではMedium/Plateau画素（＝早期相特性が“Medium”に該当しかつ遅延相特性が“Plateau”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ１０１ではMedium/Plateau画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００７０】
ステップＳ１０３ではMedium/Persistent画素（＝早期相特性が“Medium”に該当しかつ遅延相特性が“Persistent”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ１０５ではMedium/Persistent画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００７１】
ステップＳ１０７ではSlow/Washout画素（＝早期相特性が“Slow”に該当しかつ遅延相特性が“Washout”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ１０９ではSlow/Washout画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００７２】
ステップＳ１１１ではSlow/Plateau画素（＝早期相特性が“Slow”に該当しかつ遅延相特性が“Plateau”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ１１３ではSlow/Plateau画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００７３】
ステップＳ１１５ではSlow/Persistent画素（＝早期相特性が“Slow”に該当しかつ遅延相特性が“Persistent”に該当する微小検査部位に対応する画素）の数をカウントし、ステップＳ１１７ではSlow/Persistent画素がステップＳ７９で定義された解析領域に占める割合を算出する。
【００７４】
ステップＳ１１９では、こうして算出された画素数および割合を参照して、図９に示す血流特性マトリクスを作成する。作成された血流特性マトリクスでは、横方向に早期相が割り当てられ、縦方向に遅延相が割り当てられる。また、作成された血流特性マトリクスは、ワークエリア２２ｗに保持される。
【００７５】
図９を参照して、Fast/Washout画素の割合および画素数は、Fast/Washout画素の面積とともに左上段のマトリクス要素に記載される。Fast/Plateau画素の割合および画素数は、Fast/Plateau画素の面積とともに左中段のマトリクス要素に記載される。Fast/Persistent画素の割合および画素数は、Fast/Persistent画素の面積とともに左下段のマトリクス要素に記載される。
【００７６】
Medium/Washout画素の割合および画素数は、Medium/Washout画素の面積とともに中央上段のマトリクス要素に記載される。Medium/Plateau画素の割合および画素数は、Medium/Plateau画素の面積とともに中央中段のマトリクス要素に記載される。Medium/Persistent画素の割合および画素数は、Medium/Persistent画素の面積とともに中央下段のマトリクス要素に記載される。
【００７７】
Slow/Washout画素の割合および画素数は、Slow/Washout画素の面積とともに右上段のマトリクス要素に記載される。Slow/Plateau画素の割合および画素数は、Slow/Plateau画素の面積とともに右中段のマトリクス要素に記載される。Slow/Persistent画素の割合および画素数は、Slow/Persistent画素の面積とともに右下段のマトリクス要素に記載される。
【００７８】
図１８に戻って、ステップＳ１２１では、ワークエリア２２ｗ上の限定カラーマップＲＭｔ０１，ＲＭｔ１２および血流特性マトリクスを構成要素とする解析情報（図１０参照）を、ビューワエリア２２ｖの右下段に配された小エリアに展開する。展開された解析情報はビューワ２０によって読み出され、図１１に示す要領で画面の右下段に表示される。
【００７９】
ステップＳ１２３ではキーボード／マウス１６によって保存操作が行われたか否かを判別し、判別結果がＮＯであればそのままステップＳ１２９に進む。これに対して、判別結果がＹＥＳであれば、ステップＳ１２５で解析情報の保存を通信Ｉ／Ｆ１２に命令し、ステップＳ１２７で解析情報に関連するカラーマップＭＰｔ０１およびＭＰｔ１２の保存を通信Ｉ／Ｆ１２に命令する。ステップＳ１２７の処理が完了すると、ステップＳ１２７に進む。
【００８０】
通信Ｉ／Ｆ１２は、ビューワエリア２２ｖの左下段に配された小エリアから解析情報を読み出し、ビューワエリア２２ｖの右下段および中央下段に配された２つの小エリアからカラーマップＭＰｔ０１およびＭＰｔ１２を読み出し、こうして読み出された解析情報，カラーマップＭＰｔ０１およびＭＰｔ１２を院内ＬＡＮ５０を経て画像サーバ４０に転送する。転送された解析情報，カラーマップＭＰｔ０１およびＭＰｔ１２は、今後のフォローアップ時の材料として残すために、画像サーバ４０に保存される。
【００８１】
ステップＳ１２９ではキーボード／マウス１６によって終了操作が行われたか否かを判別し、判別結果がＮＯであればステップＳ１２３に戻る一方、判別結果がＹＥＳであれば上階層のルーチンに復帰する。
【００８２】
以上の説明から分かるように、ＣＰＵ１４は、既定条件を満足する１または２以上の画素にそれぞれ対応する１または２以上の有色画素からなるカラーマップＭＰｔ０１，ＭＰｔ１２を時刻Ｔ０〜Ｔ２にそれぞれ撮像されかつ共通の検査部位を表す断層画像ＩＭｔ０〜ＩＭｔ２に基づいてビューワ２０に表示する(S27~S33)。表示されたカラーマップＭＰｔ０１，ＭＰｔ１２上の所望の有色画素を指定するクリック操作を受け付けると(S15, S17)、ＣＰＵ１４は、クリック操作によって指定された有色画素から連続する１または２以上の有色画素をカラーマップＭＰｔ０１，ＭＰｔ１２上で探索し(S41~S45, S49~S51, S57~S63, S71~S77)、指定または探知された有色画素に対応する微小検査部位の血流特性を断層画像ＩＭｔ０〜ＩＭｔ２に基づいて解析する(S47, S67)。ＣＰＵ１４はさらに、解析結果が特定の血流特性を示す微小検査部位が解析対象に占める割合を含む解析情報をビューワ２０に追加表示する(S79~S121)。
【００８３】
探索処理ではクリック操作によって指定された有色画素に連続する１または２以上の有色画素が探索されるところ、これらの有色画素の分布態様（つまり、指定された有色画素および探知された有色画素の総数ないし面積）は、検査部位の状態によって異なる。また、診断者によって関心があるのは、特定の血流特性を示す微小検査部位である。
【００８４】
これを踏まえて、指定された有色画素および探知された有色画素の各々に対応する微小検査部位の血流特性を解析し、解析結果が特定の血流特性を示す微小検査部位が解析対象に占める割合を表示するようにしている。これによって、解析結果が特定の血流特性を示す微小検査部位の大きさが正規化され、検査部位の血流特性を評価する診断者間での評価のばらつきを抑制することができる。
【符号の説明】
【００８５】
１０  …画像診断支援装置
１４  …ＣＰＵ
１６  …キーボード／マウス
２０  …ビューワ
２２  …ＤＲＡＭ
３０  …ＭＲＩ装置
４０  …画像サーバ
５０  …院内ＬＡＮ
【図１】
2017136317.tif 000003
【図２】
2017136317.tif 000004
【図３】
2017136317.tif 000005
【図４】
2017136317.tif 000006
【図５】
2017136317.tif 000007
【図６】
2017136317.tif 000008
【図７】
2017136317.tif 000009
【図８】
2017136317.tif 000010
【図９】
2017136317.tif 000011
【図１０】
2017136317.tif 000012
【図１１】
2017136317.tif 000013
【図１２】
2017136317.tif 000014
【図１３】
2017136317.tif 000015
【図１４】
2017136317.tif 000016
【図１５】
2017136317.tif 000017
【図１６】
2017136317.tif 000018
【図１７】
2017136317.tif 000019
【図１８】
2017136317.tif 000020