(57)【要約】      （修正有）【課題】得られた機能画像から疾患が各ステージに該当する可能性を表示でき、ユーザによるステージ判定を支援することができるステージ判定支援システムを提供する。【解決手段】ステージ判定支援システムは、記憶手段、第１判定手段、第２判定手段、算出手段及び表示手段を備えている。前記記憶手段は、被検体の複数種の機能画像を記憶する。前記第１判定手段は、前記各機能画像に基づいて、疾患に対する複数種の解析指標の各々の度合を判定する。前記第２判定手段は、前記第１判定手段による各々の判定結果と所定のステージ判定ルールとに基づいて、解析指標毎に、当該判定結果に該当するステージを判定する。前記算出手段は、前記第２判定手段による各々の判定結果に基づいて、前記ステージ毎に、前記疾患が当該ステージに該当する可能性を示す値を算出する。前記表示手段は、前記算出された値を前記ステージ毎に表示する。【選択図】図１
2014104097.tif 000002
【特許請求の範囲】
【請求項１】
重症度を示す複数のステージのうち、いずれのステージに被検体の疾患が該当するかの判定を支援するためのステージ判定支援システムであって、
前記被検体の複数種の機能画像を記憶する記憶手段と、
前記各機能画像に基づいて、前記被検体の正常とみなされる領域に対する複数種の解析指標の各々に比較して、前記疾患に対する前記複数種の解析指標の各々の度合を判定する第１判定手段と、
前記第１判定手段による各々の判定結果と所定のステージ判定ルールとに基づいて、前記解析指標毎に、当該判定結果に該当するステージを判定する第２判定手段と、
前記第２判定手段による各々の判定結果に基づいて、前記ステージ毎に、前記疾患が当該ステージに該当する可能性を示す値を算出する算出手段と、
前記算出された値を前記ステージ毎に表示する表示手段と
を備えたことを特徴とするステージ判定支援システム。
【請求項２】
請求項１に記載のステージ判定支援システムにおいて、
前記疾患は脳虚血疾患であることを特徴とするステージ判定支援システム。
【請求項３】
請求項２に記載のステージ判定支援システムにおいて、
前記複数種の解析指標は、脳血液量、脳血流量、脳酸素摂取率及び脳酸素代謝量であり、
前記ステージ判定ルールは、前記脳血液量、前記脳血流量、前記脳酸素摂取率及び前記脳酸素代謝量に基づいて、前記脳虚血疾患に対するパワーズのステージ分類を実行するための判定ルールであることを特徴とするステージ判定支援システム。
【請求項４】
請求項１乃至請求項３のいずれか１項に記載のステージ判定支援システムにおいて、
前記算出手段は、前記ステージ毎に、前記第２判定手段による各々の判定結果を重み付け加算することにより、前記値を算出することを特徴とするステージ判定支援システム。
【請求項５】
請求項１乃至請求項４のいずれか１項に記載のステージ判定支援システムにおいて、
前記第１判定手段は、前記各機能画像にユーザが設定した関心領域、前記各機能画像が表す正常領域との画素値の比較によって異常と判定された領域、前記各機能画像が表す組織全体を複数のセグメントに分割した領域、又は前記各機能画像上の各画素単位、を判定対象領域として、前記度合の判定を実行することを特徴とするステージ判定支援システム。
【請求項６】
請求項１乃至請求項４のいずれか１項に記載のステージ判定支援システムにおいて、
前記第１判定手段は、前記各機能画像において判定対象領域と前記正常とみなされる領域との画素値の比較により、前記度合の判定を実行することを特徴とするステージ判定支援システム。
【請求項７】
請求項１乃至請求項６のいずれか１項に記載のステージ判定支援システムにおいて、
前記正常とみなされる領域は、前記各機能画像において判定対象領域と左右対称の関係にある対側領域であることを特徴とするステージ判定支援システム。
【請求項８】
請求項１乃至請求項６のいずれか１項に記載のステージ判定支援システムにおいて、
前記正常とみなされる領域は、小脳の一部領域であることを特徴とするステージ判定支援システム。
【請求項９】
請求項１乃至請求項８のいずれか１項に記載のステージ判定支援システムにおいて、
前記各機能画像の位置合わせを実行する位置合わせ手段を更に備え、
前記第１判定手段は、前記位置合わせが実行された各機能画像に基づいて、前記度合の判定を実行することを特徴とするステージ判定支援システム。
【請求項１０】
請求項１乃至請求項９のいずれか１項に記載のステージ判定支援システムにおいて、
前記表示手段は、前記算出された値に基づいて、数値表示、グラフ表示又はカラースケール表示を実行することを特徴とするステージ判定支援システム。
【請求項１１】
請求項１乃至請求項４のいずれか１項に記載のステージ判定支援システムにおいて、
前記第１判定手段は、前記各機能画像上の各画素単位を対象として、前記度合の判定を実行し、
前記第２判定手段は、前記各画素単位を対象として、前記ステージの判定を実行し、
前記算出手段は、前記各画素単位を対象として、前記ステージ毎に前記値を算出し、
前記表示手段は、前記各画素単位を対象として、前記ステージ毎に前記算出された値のうち、最も高い値に該当するステージを表すように色分け表示を実行することを特徴とするステージ判定支援システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、ステージ判定支援システムに関する。
【背景技術】
【０００２】
近年、各種の機能画像を観察することにより、心臓や脳などの疾患の判定が可能となってきている。例えば、急性期脳梗塞、一過性虚血発作、慢性期頭頚部血管狭窄などの脳虚血疾患を診断する場合、各種の医用画像診断装置により得られた２次元又は３次元画像としての灌流（Perfusion）画像に基づいて、脳循環代謝量が評価される。
【０００３】
脳循環代謝量の具体的な解析指標としては、脳血流量ＣＢＦ［ml/100g/min］、脳血液量ＣＢＶ［ml/100g］、脳酸素摂取率ＯＥＦ［％］及び脳酸素代謝量ＣＭＲＯ２［ml/100g/min］等の種類がある。ここで、ＣＢＦは、Cerebral Blood Flowの略語であり、ＣＢＶは、Cerebral Blood Volumeの略語である。ＯＥＦは、Cerebral Oxygen Extraction Fractionの略語であり、ＣＭＲＯ２は、Cerebral Metabolic Rate of Oxygenの略語である。
【０００４】
医師は、通常、虚血部位におけるこれらの解析指標の大小を複数の機能画像上で観察することにより、脳循環代謝の重症度を総合的に判断し、治療方針を決定する。
【０００５】
各解析指標のうち、脳血流量ＣＢＦ及び脳血液量ＣＢＶの大小は、灌流画像上で観察される。なお、灌流画像は、Ｘ線ＣＴ（computed tomography）装置、ＭＲＩ（magnetic resonance imaging）装置、ＳＰＥＣＴ（single photon emission CT）装置又はＰＥＴ（positron emission tomography）装置等から得られる。
【０００６】
また、各解析指標のうち、脳酸素摂取率ＯＥＦ及び脳酸素代謝量ＣＭＲＯ２の大小は、特にＰＥＴ装置から得られるＰＥＴ画像上で観察される。
【０００７】
各解析指標の大小の観察後、医師が脳循環代謝の重症度を判断する基準となる１つの指標として、パワーズ（Powers）のステージ分類と呼ばれる代表的な重症度判定の考え方がある。パワーズのステージ分類は、上述した各種の解析指標の大小の組合せにより、重症度のステージを３段階に分ける指標である。この指標は、具体的には、上述した４つの解析指標（ＣＢＦ、ＣＢＶ、ＯＥＦ、ＣＭＲＯ２）の大小により、図１１に示すように、脳循環代謝の重症度をステージＩ、ＩＩ、ＩＩＩに分類して示している。一般に、ステージＩＩＩは梗塞（脳細胞壊死）のため治療不能な状態とみなされる。ステージＩＩは虚血性ペナンブラ（可逆性組織障害）と呼ばれる状態に相当し、治療による血流回復により正常化可能な状態とみなされる。虚血性ペナンブラ領域の残っている割合が治療方針の決定に重要である。
【０００８】
具体的なステージ分類の方法を図１２により説明する。図１２に示す例は、患側（脳の虚血が起きている側）の脳血流量ＣＢＦが減少し、脳血液量ＣＢＶが増加し、脳酸素摂取率ＯＥＦが増加し、脳酸素代謝量ＣＭＲＯ２が減少しているので、ステージＩＩＩに相当し、脳細胞壊死状態とみなされる。
【先行技術文献】
【非特許文献】
【０００９】
【非特許文献１】Powers, et al. Stroke 1985, 16:361-376
【発明の概要】
【発明が解決しようとする課題】
【００１０】
しかしながら、以上のようなステージ分類の方法は、４つの解析指標の大小の組合せがどのステージに相当するのかよく考えないと判定できないため、判定に時間がかかったり、誤判定したりする可能性がある。
【００１１】
また、実際の臨床においては、どのステージにも該当しない組合せが起こり得るため、そのような場合、どのような診断結果を下せば良いのか悩むことになる。
【００１２】
さらに、全ての機能画像が得られていない場合や、逆に１つの解析指標に異なる医用画像診断装置から得られた複数の機能画像が得られている場合、過不足な情報に対して主観的に判定しがちになり、適切な診断結果を下せない可能性がある。
【００１３】
目的は、得られた機能画像から疾患が各ステージに該当する可能性を表示でき、ユーザによるステージ判定を支援し得るステージ判定支援システムを提供することである。
【課題を解決するための手段】
【００１４】
実施形態のステージ判定支援システムは、重症度を示す複数のステージのうち、いずれのステージに被検体の疾患が該当するかの判定を支援するためのシステムである。
【００１５】
前記ステージ判定支援システムは、記憶手段、第１判定手段、第２判定手段、算出手段及び表示手段を備えている。
【００１６】
前記記憶手段は、前記被検体の複数種の機能画像を記憶する。
【００１７】
前記第１判定手段は、前記各機能画像に基づいて、前記被検体の正常とみなされる領域に対する複数種の解析指標の各々に比較して、前記疾患に対する前記複数種の解析指標の各々の度合を判定する。
【００１８】
前記第２判定手段は、前記第１判定手段による各々の判定結果と所定のステージ判定ルールとに基づいて、前記解析指標毎に、当該判定結果に該当するステージを判定する。
【００１９】
前記算出手段は、前記第２判定手段による各々の判定結果に基づいて、前記ステージ毎に、前記疾患が当該ステージに該当する可能性を示す値を算出する。
【００２０】
前記表示手段は、前記算出された値を前記ステージ毎に表示する。
【図面の簡単な説明】
【００２１】
【図１】第１の実施形態に係るステージ判定支援システム及びその周辺構成を示す模式図である。
【図２】同実施形態における動作を説明するためのフローチャートである。
【図３】同実施形態における脳を分割したセグメントの一例を示す模式図である。
【図４】一般的な虚血性ペナンブラと中心部梗塞領域を説明するための模式図である。
【図５】同実施形態におけるステージ分類上の該当範囲の決定を説明するための模式図である。
【図６】同実施形態におけるステージ可能性の計算方法を説明するための模式図である。
【図７】同実施形態におけるステージ可能性の他の計算方法を説明するための模式図である。
【図８】第２の実施形態における動作を説明するためのフローチャートである。
【図９】同実施形態における各ステージの可能性が高い領域の表示画像の一例を示す模式図である。
【図１０】同実施形態における各ステージの可能性が高い領域を色分け表示した画像の一例を示す模式図である。
【図１１】一般的なパワーズのステージ分類を示す模式図である。
【図１２】一般的な内頚動脈狭窄症における脳循環代謝量の一例を示す図である。
【発明を実施するための形態】
【００２２】
以下、各実施形態に係るステージ判定支援システムについて図面を用いて説明する。以下のステージ判定支援システムは、それぞれハードウェア構成、又はハードウェア資源とソフトウェアとの組合せ構成のいずれでも実施可能となっている。組合せ構成のソフトウェアとしては、予めネットワーク又は記憶媒体からコンピュータにインストールされ、ステージ判定支援システムの各機能を当該コンピュータに実現させるためのステージ判定支援プログラムが用いられる。
【００２３】
＜第１の実施形態＞
図１は第１の実施形態に係るステージ判定支援システム及びその周辺構成を示す模式図である。ステージ判定支援システムは、重症度を示す複数のステージのうち、いずれのステージに被検体の疾患が該当するかの判定を支援するためのシステムである。ステージ判定支援システムは、例えば、医用画像ワークステーション１０、入力装置（例えば、キーボード又はマウス等）１４及びモニタ装置１５によって構成されている。医用画像ワークステーション１０の内部にはハードディスク１１、コンピュータ１２及びネットワークＩ／Ｆ（例えば、Ethernet（登録商標） Cardなど）１３がある。なお、ステージ判定支援システムにおける「システム」の用語は「装置」と読み替えてもよく、入力装置１４及びモニタ装置１５における「装置」の用語は「部」と読み替えてもよい。また、ステージ判定支援システムは、各装置２０〜５０からネットワークを介して機能画像１１ｇを収集する画像サーバ装置に接続され、当該画像サーバ装置から機能画像１１ｇを読出可能なクライアント装置として実現してもよい。
【００２４】
ハードディスク１１にはステージ判定支援用ソフトウェア（ステージ判定支援プログラム）１１ｓと、被検体（患者）の複数種の機能画像１１ｇが保管されている。コンピュータ１２は、メインメモリ１２ｍ及びプロセッサ１２ｐを備えている。
【００２５】
プロセッサ１２ｐは、ハードディスク１１内のステージ判定支援用ソフトウェア１１ｓと各機能画像１１ｇをメインメモリ（記憶手段）１２ｍに読み出してステージ判定支援用ソフトウェア１１ｓの処理を実行する。
【００２６】
プロセッサ１２ｐは、各機能画像１１ｇをネットワークＩ／Ｆ１３によりネットワークを介してＸ線ＣＴ装置２０、ＭＲＩ装置３０、ＳＰＥＣＴ装置４０及びＰＥＴ装置５０などから取得し、ハードディスク１１に保管しておく。
【００２７】
但し、プロセッサ１２ｐは、各機能画像１１ｇをハードディスク１１に保管せずにメインメモリ１２ｍに転送し、ハードディスク１１から読み出したステージ判定支援用ソフトウェア１１ｓを実行して、メインメモリ１２ｍ内の各機能画像１１ｇを解析しても良い。
【００２８】
なお、ステージ判定支援用ソフトウェア１１ｓは、プロセッサ１２ｐに実行され、コンピュータ１２に第１判定機能、第２判定機能、算出機能及び表示機能を実現させるためのステージ判定支援プログラムを含んでいる。すなわち、ステージ判定支援用ソフトウェア１１ｓ及びプロセッサ１２ｐは、第１判定手段、第２判定手段及び算出手段を構成している。ステージ判定支援用ソフトウェア１１ｓ、プロセッサ１２ｐ及びモニタ装置１５は、表示手段を構成している。
【００２９】
ここで、第１判定機能は、メインメモリ１２ｍ内の各機能画像１１ｇに基づいて、被検体の正常とみなされる領域に対する複数種の解析指標の各々に比較して、疾患に対する複数種の解析指標の各々の度合を判定する機能である。疾患としては、例えば、脳虚血疾患が適用可能となっている。複数種の解析指標としては、例えば、脳血液量ＣＢＶ、脳血流量ＣＢＦ、脳酸素摂取率ＯＥＦ及び脳酸素代謝量ＣＭＲＯ２が適用可能となっている。また、第１判定機能としては、例えば以下の（ａ），（ｂ），（ｃ）又は（ｄ）を判定対象領域として、度合の判定を実行する機能としてもよい。なお、本実施形態では（ａ）を判定対象領域とした場合を例に挙げて述べている。
【００３０】
（ａ）各機能画像１１ｇにユーザが設定した関心領域ＲＯＩ。
【００３１】
（ｂ）各機能画像１１ｇが表す正常領域との画素値の比較によって異常と判定された領域。
【００３２】
（ｃ）各機能画像１１ｇが表す組織全体を複数のセグメントに分割した領域。
【００３３】
（ｄ）各機能画像１１ｇ上の各画素単位。
【００３４】
また、第１判定機能としては、各機能画像１１ｇにおいて当該判定対象領域（上記（ａ），（ｂ），（ｃ）又は（ｄ））と正常とみなされる領域との画素値の比較により、度合の判定を実行する機能としてもよい。ここで、正常とみなされる領域としては、例えば、判定対象領域と左右対称の関係にある対側領域としてもよく、又は、（梗塞しにくく、虚血に陥る可能性が低い）小脳の一部領域としてもよい。
【００３５】
また、ステージ判定支援プログラムが各機能画像１１ｇの位置合わせを実行する位置合わせ機能を更にコンピュータ１２に実現させる場合には、第１判定機能は、当該位置合わせが実行された各機能画像１１ｇに基づいて、度合の判定を実行する機能としてもよい。
【００３６】
第２判定機能は、第１判定機能による各々の判定結果と所定のステージ判定ルールとに基づいて、解析指標毎に、当該判定結果に該当するステージを判定する機能である。ステージ判定ルールとしては、例えば、脳血液量ＣＢＶ、脳血流量ＣＢＦ、脳酸素摂取率ＯＥＦ及び脳酸素代謝量ＣＭＲＯ２に基づいて、脳虚血疾患に対するパワーズのステージ分類を実行するための判定ルールが適用可能となっている。
【００３７】
算出機能は、第２判定機能による各々の判定結果に基づいて、ステージ毎に、疾患が当該ステージに該当する可能性を示す値を算出する機能である。この算出機能としては、例えば、ステージ毎に、第２判定機能による各々の判定結果を重み付け加算することにより、可能性を示す値を算出する機能としてもよい。
【００３８】
表示機能は、当該算出された値をステージ毎にモニタ装置１５に表示する機能である。表示機能は、例えば、当該算出された値に基づいて、数値表示、グラフ表示又はカラースケール表示をモニタ装置１５により実行する機能としてもよい。
【００３９】
次に、以上のように構成されたステージ判定支援システムの動作について図２のフローチャートを用いて説明する。なお、ステージ判定支援システムの医用画像ワークステーション１０は、ユーザによる入力装置１４の操作により、プロセッサ１２ｐがハードディスク１１内のステージ判定支援用ソフトウェア１１ｓをメインメモリ１２ｍに読み出して実行しているものとする。
【００４０】
ステップＳＴ１において、医用画像ワークステーション１０は、ユーザによる入力装置１４の操作により、プロセッサ１２ｐが、ユーザが指定した脳循環代謝量の各解析指標を担当する各機能画像１１ｇをハードディスク１１からメインメモリ１２ｍに読み出す。さらに、プロセッサ１２ｐは、各機能画像１１ｇをモニタ装置１５に送出し、モニタ装置１５は、各機能画像１１ｇを表示する。
【００４１】
続いて、プロセッサ１２ｐは、解剖上の形態に基づいて各機能画像１１ｇの位置合わせを実行する。位置合わせの方法としては、公知の代表的な標準脳との位置合わせ法である３Ｄ−ＳＳＰ（stereotactic surface projection）法又はＳＰＭ（statistical parametric mapping）法などが使用可能となっている。
【００４２】
３Ｄ−ＳＳＰ法は、例えば、文献［１］に記載されている。
【００４３】
［１］S Minoshima, RA Koeppe, KA Frey, et al., “Anatomic standardization: Linear scaling and nonlinear warping of functional brain images”, J Nucl Med, 35（9）, 1528-1537,（1994）.
ＳＰＭ法は、例えば、文献［２］及び［３］に記載されている。
【００４４】
［２］KJ Friston, J Ashburner, CD Frith, et al., “Spatial registration and normalization of images”, Human Brain Mapping, 2, 165-189,（1995）.
［３］J Ashburner, and KJ Friston, “Nonlinear spatial normalization using basis functions”, Hum Brain Mapp, 7（4）, 254-266,（1999）.
１つの解析指標を担当する機能画像１１ｇは１つでもよく、複数でもよく、０個でもよい。具体的には、脳血流量ＣＢＦ及び脳血液量ＣＢＶを担当する機能画像１１ｇとしては、主にＸ線ＣＴ装置２０、ＭＲＩ装置３０、ＳＰＥＣＴ装置４０又はＰＥＴ装置５０から得られる灌流画像を用いる。また、脳酸素摂取率ＯＥＦ及び脳酸素代謝量ＣＭＲＯ２を担当する機能画像１１ｇとしてはＰＥＴ画像を用いる。また、図１１から分かるように、脳酸素代謝量ＣＭＲＯ２の減少が脳細胞壊死と同等な状態とみなせるため、機能画像１１ｇの代わりにＭＲＩの拡散強調画像ＤＷＩ（Diffusion Weighted Image）の高信号値を解析指標として代用してもよい。
【００４５】
なお、同一の医用画像診断装置２０〜５０で１回の撮像で得られたデータを用いて解析して得られた複数の機能画像１１ｇ同士であれば、位置ズレが発生しないため、本位置合わせ処理は不要である。
【００４６】
ステップＳＴ２において、プロセッサ１２ｐは、ユーザが指定した全ての機能画像１１ｇをモニタ装置１５の画面上に表示する。プロセッサ１２ｐは、ユーザによる入力装置１４の操作により、表示された任意の機能画像１１ｇ上の、対側と比較して低信号又は高信号で虚血状態と思われる領域に関心領域ＲＯＩ（Region of Interest）を設定する。
【００４７】
なお、関心領域ＲＯＩの形状は楕円等でも良いし、任意形状でも良い。また、プロセッサ１２ｐは、内部処理により信号値の左右差が大きい領域を検出してその領域に自動的に関心領域ＲＯＩを設定しても良い。或いは、図３に一例を示すように、予め脳全体（又は脳の外周側領域のみ）を複数の細かいセグメントＬ１，Ｌ２，…，Ｒ１，Ｒ２，…に分割しておき、機能画像１１ｇ上で複数のセグメントＬ１，Ｌ２，…，Ｒ１，Ｒ２，…に個別に対応した複数の領域に複数の関心領域ＲＯＩを設定しても良い。なお、各セグメントＬ１，Ｌ２，…，Ｒ１，Ｒ２，…は、左右対称に分割しておくことが、関心領域ＲＯＩを左右対称に設定する観点から好ましい。
【００４８】
また、ユーザによる入力装置１４の操作に応じて、複数の関心領域ＲＯＩを設定することも可能とする。ユーザが設定した複数の関心領域ＲＯＩが重なっている場合は、大きい方の関心領域ＲＯＩは、小さい方の関心領域ＲＯＩの領域を除いた領域とする。これは、図４に示すように、脳虚血領域ｒ１における中心部梗塞（脳細胞壊死）領域（ischemic core）ｒ２及びそれ以外の虚血性ペナンブラ領域（ischemic penumbra）ｒ３といった各領域ｒ２，ｒ３にそれぞれステージ判定を行いたい場合を想定している。
【００４９】
ステップＳＴ３において、プロセッサ１２ｐは、患側（脳の虚血側）に設定された関心領域ＲＯＩに対応する関心領域ＲＯＩを健側（脳の正常側）に左右対称に設定し、各機能画像１１ｇに対して左右の関心領域ＲＯＩの平均値を比較する。プロセッサ１２ｐは、健側と比較して患側の関心領域ＲＯＩ平均値が大きいか小さいか同程度かの度合を判定する。大、小、同程度の判定の閾値は、健側との比が１２０％以上は大、７０％未満は小、７０％以上〜１２０％未満は同程度などとして判定する。また、両側が虚血に陥っている場合も考えられる。両側が虚血状態の場合にも判定可能とするため、左右比ではなく、虚血に陥る可能性が低い小脳などに正常領域の関心領域ＲＯＩを設定して虚血領域との比較をしても良い。
【００５０】
ステップＳＴ４において、プロセッサ１２ｐは、各機能画像１１ｇについて、健側との比較による判定結果（大、小、同程度）及びパワーズのステージ分類上のステージ判定ルールに基づいて、解析指標毎に、当該判定結果に該当するステージを判定する。例えば図５に示す場合、脳血液量ＣＢＶは大と判定されたため、ステージＩ、ＩＩ、ＩＩＩの全てのステージに該当することが判定され、該当範囲がステージＩ〜ＩＩＩと決定される。脳血流量ＣＢＦは小と判定されたため、ステージＩＩ、ＩＩＩに該当することが判定され、該当範囲がステージＩＩ〜ＩＩＩと決定される。脳酸素摂取率ＯＥＦは大と判定されたため、ステージＩＩ、ＩＩＩに該当することが判定され、該当範囲がステージＩＩ〜ＩＩＩと決定される。脳酸素代謝量ＣＭＲＯ２は同程度と判定されたため、ステージＩ、ＩＩに該当することが判定され、該当範囲がステージＩ〜ＩＩと決定される。
【００５１】
ステップＳＴ５において、プロセッサ１２ｐは、ユーザによる入力装置１４の操作により、予め各解析指標の重み係数を設定する。解析指標毎の重み係数は、ユーザの考える重要度に応じた割合に設定される。解析指標毎の重み係数の合計は、この例では１とするが、これに限らず、例えば１００でもよい（予め百分率で重み付け加算する場合）。
【００５２】
プロセッサ１２ｐは、ステージ毎に、各解析指標が該当か非該当の判定結果に応じて１か０を設定し、これら各設定値（１又は０）に重み係数（０．３、０．３、０．１又は０．３）をかけて加算した重み付け加算値を各ステージの可能性を示す値として求める。例えば図６に示す場合、ステージ１の重み付け加算値は、０．３×１＋０．３×０＋０．１×０＋０．３×１＝０．６と計算され、１００を乗じて百分率で表すと、可能性を示す値が６０％となる。ステージＩＩ、ＩＩＩについても同様の計算を行い、それぞれ１００％、７０％となる。なお、このような重み付け加算を行わない場合には、全ての重み係数を１に設定すればよい。
【００５３】
ここで、対象の機能画像１１ｇが存在しない解析指標は、各ステージＩ〜ＩＩＩの可能性を０としてもよく、どのステージＩ〜ＩＩＩの可能性も均等にあるとして１／３として計算してもよい。また、１つの解析指標に異なる医用画像診断装置から得られた複数の機能画像１１ｇが存在する場合は、単純に複数の結果を加算すれば良い。このように、複数の機能画像１１ｇの結果を加算することにより、その解析指標のステージ可能性の確からしさの信頼性が増す効果が得られる。
【００５４】
ステップＳＴ６において、プロセッサ１２ｐは、算出した値をステージ毎にモニタ装置１５に表示する。表示方法としては、図６に示すようにパワーズのステージ分類図とともに数値で表示しても良いし、棒グラフ、円グラフ又はカラースケール画像など別の表現方法で表示しても良い。
【００５５】
上述したように本実施形態によれば、各機能画像１１ｇに基づいて解析指標の各々の度合を判定し、各々の判定結果と所定のステージ判定ルールとに基づいて、解析指標毎に、当該判定結果に該当するステージを判定し、各々の判定結果に基づいて、ステージ毎に、該当する可能性を示す値を算出して表示する構成により、得られた機能画像１１ｇから疾患が各ステージに該当する可能性を表示でき、ユーザによるステージ判定を支援することができる。
【００５６】
また、これにより、ひと目でどのステージの可能性がどの程度なのかが分かるので、ユーザが画像を見て頭で考えて判定する必要がなくなり、即座に客観的にステージ分類を判定できるようになり、誤判定もなくなる。また、解析指標の組合せがどのステージにも該当しない場合でも、ユーザが各ステージの可能性の高さを客観的に見て、ステージ分類を判定することができる。さらに、複数モダリティの機能画像１１ｇが得られている場合には、可能性を示す値の信頼性が上がるため、判定の精度を更に向上させることができる。
【００５７】
なお、本実施形態は図７に示すように変形してもよい。図７に示す例の場合、ステップＳＴ４において健側との比率により大、小、同程度の度合を判定するのではなく、健側との比率自体を直接ステージ判定に用いる。例えば、比率の最大値（例、＋１００％）と最小値（例、−１００％）をユーザが予め設定しておく。プロセッサ１２ｐは、比率の最大値又は最小値を限界値として健側との比率を示す値と、各解析指標のステージに該当する当該比率を示す値の範囲とに基づいて、各解析指標に対してステージを判定する。
【００５８】
図７に示す例では、脳血液量ＣＢＶにおいて、ステージＩに該当する比率を示す値の範囲が０％〜＋７４％であり、ステージＩＩ，ＩＩＩに該当する比率を示す値の範囲が＋７５％である。
【００５９】
脳血流量ＣＢＦにおいては、ステージＩに該当する比率を示す値の範囲が０％であり、ステージＩＩに該当する比率を示す値の範囲が−１％〜−５０％であり、ステージＩＩＩに該当する比率を示す値の範囲が−５１％〜−１００％である。
【００６０】
脳酸素摂取率ＯＥＦにおいては、ステージＩに該当する比率を示す値の範囲が０％であり、ステージＩＩに該当する比率を示す値の範囲が＋１％〜＋９９％であり、ステージＩＩＩに該当する比率を示す値の範囲が＋１００％である。
【００６１】
脳酸素代謝量ＣＭＲＯ２においては、ステージＩ，ＩＩに該当する比率を示す値の範囲が０％であり、ステージＩＩＩに該当する比率を示す値の範囲が−１％〜−１００％である。
【００６２】
従って、図７に示す例の場合、脳血液量ＣＢＶは、健側との比率を示す値が＋７０％であるため、０％〜＋７４％の範囲に該当するステージＩと判定される。脳血流量ＣＢＦ、脳酸素摂取率ＯＥＦ及び脳酸素代謝量ＣＭＲＯ２についても同様の方法で、それぞれステージＩＩ、ステージＩＩ及びステージＩ，ＩＩと判定される。各ステージの可能性を示す値の計算方法及び表示方法は、ステップＳＴ５〜ＳＴ６に記載の方法と同じである。
【００６３】
＜第２の実施形態＞
次に、第２の実施形態に係るステージ判定支援システムについて図１を用いて説明する。なお、前述した部分と同一部分についてはその詳しい説明を省略し、ここでは異なる部分について主に述べる。
【００６４】
第２の実施形態は、第１の実施形態の変形例であり、関心領域ＲＯＩを用いる第１の実施形態とは異なり、各機能画像１１ｇ上の各画素単位を対象として、前述した第１判定機能、第２判定機能、算出機能及び表示機能を実行する形態である。
【００６５】
すなわち、第１判定機能は、前述した第１判定機能において、各機能画像１１ｇ上の各画素単位を対象として、度合の判定を実行する。
【００６６】
第２判定機能は、前述した第２判定機能において、各画素単位を対象として、ステージの判定を実行する。
【００６７】
算出機能は、前述した算出機能において、各画素単位を対象として、ステージ毎に値を算出する。
【００６８】
表示機能は、前述した表示機能において、各画素単位を対象として、ステージ毎に当該算出された値のうち、最も高い値に該当するステージを表すように色分け表示を実行する機能となっている。
【００６９】
他の構成は、第１の実施形態と同様である。
【００７０】
次に、以上のように構成されたステージ判定支援システムの動作について図８のフローチャートを用いて説明する。
【００７１】
始めに、開始時からステップＳＴ１の各機能画像１１ｇの位置合わせまでの動作は、第１の実施形態と同様に実行される。また、ステップＳＴ１の完了後、前述したステップＳＴ２（関心領域ＲＯＩの設定）は行わない。
【００７２】
ステップＳＴ３ａにおいて、プロセッサ１２ｐは、前述したステップＳＴ３に相当する処理（各機能画像１１ｇに対する健側との比較による度合の判定）を左右対称の全ての画素毎に行う。
【００７３】
ステップＳＴ４ａにおいて、プロセッサ１２ｐは、前述したステップＳＴ４に相当する処理（度合の判定結果に基づくステージの判定）を左右対称の全ての画素毎に行う。
【００７４】
ステップＳＴ５ａにおいて、プロセッサ１２ｐは、前述したステップＳＴ５に相当する処理（ステージの判定結果に基づく、各ステージの可能性を示す重み付け加算値の計算）を左右対称の全ての画素毎に行う。
【００７５】
ステップＳＴ６ａにおいて、プロセッサ１２ｐは、画素毎にどのステージの可能性を示す値が最も高い値かを判定し、判定結果が同じ値の画素を３つの集合領域に領域分割（クラスタリング）する。
【００７６】
ステップＳＴ６ｂにおいて、プロセッサ１２ｐは、可能性が最も高いステージを表す領域分割結果をモニタ装置１５により色分け表示する。例えば、プロセッサ１２ｐは、図９に示すように、中央の赤色領域（ステージＩＩＩ）の周囲を黄土色領域（ステージＩＩ）及び青色領域（ステージＩ）が順に囲んだ形状の領域分割結果をステップＳＴ６ａで得たとする。ステップＳＴ６ｂでは、プロセッサ１２ｐは、この領域分割結果を各ステージの可能性が高い領域として機能画像１１ｇ上に重ねた画像を作成し、この画像をモニタ装置１５により表示する。
【００７７】
ステップＳＴ６ｂの表示方法の変形例としては、ステージ毎に左右対称の全ての画素の可能性を示す値をカラースケール等の方法により任意の機能画像１１ｇ上に色分け表示した画像を表示する方法がある。図１０（ａ）〜図１０（ｃ）に示す例では、ステージＩの機能画像１１ｇ上には、青色（０％）、黄色（５０％）、茶色（７０％）及び赤色（１００％）で各画素を色分け表示した画像が表示されている。同様に、ステージＩＩの機能画像１１ｇ上には、青色（０％）、黄色（５０％）、茶色（７０％）及び紫色（９０％）で各画素を色分け表示した画像が表示されている。また、ステージＩＩＩの機能画像１１ｇ上に、緑色（３０％）、黄色（５０％）、赤色（１００％）及び紫色（９０％）で各画素を色分け表示した画像が表示されている。なお、色分けした画像は、白黒の二色図面の便宜上、各色の境界に破線を描画したが、実際には破線が無く、虹色のように不明瞭な境界を有して色が変化している。また、図１０（ａ）〜図１０（ｃ）に示すように、３つのステージのカラースケール画像を同時に表示しても良いし、各ステージのカラースケール画像を切り替えて表示してもよい。
【００７８】
上述したように本実施形態によれば、関心領域ＲＯＩに代えて、各機能画像１１ｇ上の各画素単位を対象として、ステージ判定を実行する構成により、第１の実施形態の効果に加え、前述した関心領域ＲＯＩよりも細かい画素単位で各ステージの可能性を示す値を表示することができる。
【００７９】
以上説明した少なくとも一つの実施形態によれば、各機能画像１１ｇに基づいて解析指標の各々の度合を判定し、各々の判定結果と所定のステージ判定ルールとに基づいて、解析指標毎に、当該判定結果に該当するステージを判定し、各々の判定結果に基づいて、ステージ毎に、該当する可能性を示す値を算出して表示する構成により、得られた機能画像１１ｇから疾患が各ステージに該当する可能性を表示でき、ユーザによるステージ判定を支援することができる。
【００８０】
なお、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００８１】
１０…医用画像ワークステーション、１１…ハードディスク、１１ｓ…ステージ判定支援用ソフトウェア、１１ｇ…機能画像、１２…コンピュータ、１２ｍ…メインメモリ、１２ｐ…プロセッサ、１３…ネットワークＩ／Ｆ、１４…入力装置、１５…モニタ装置、２０…Ｘ線ＣＴ装置、３０…ＭＲＩ装置、４０…ＳＰＥＣＴ装置、５０…ＰＥＴ装置。
【図１】
2014104097.tif 000003
【図２】
2014104097.tif 000004
【図３】
2014104097.tif 000005
【図４】
2014104097.tif 000006
【図５】
2014104097.tif 000007
【図６】
2014104097.tif 000008
【図７】
2014104097.tif 000009
【図８】
2014104097.tif 000010
【図９】
2014104097.tif 000011
【図１０】
2014104097.tif 000012
【図１１】
2014104097.tif 000013
【図１２】
2014104097.tif 000014