JP2004174217A

JP2004174217A - 画像処理装置、画像処理方法、及び当該画像処理をコンピュータに実行させるプログラムを格納する記録媒体

Info

Publication number: JP2004174217A
Application number: JP2003182671A
Authority: JP
Inventors: Sakura Kadowaki; さくら門脇; Kazuto Sato; 和人佐藤; Hirokazu Madokoro; 洋和間所
Original assignee: Japan Science and Technology Corp
Current assignee: Japan Science and Technology Agency
Priority date: 2002-10-01
Filing date: 2003-06-26
Publication date: 2004-06-24

Abstract

【課題】例えば医用画像において、診断等に必要な関心領域を解剖学的基準に基づいて自動抽出する画像処理装置等を提供すること。
【解決手段】学習近傍領域を可変とした自己組織化マップ（ＳＯＭ）による輪郭抽出処理によって取得された脳表面輪郭線画像と、ＳＯＭによる組織の自動分類処理によって取得された組織カラーマップと、から得られる解剖学的構造情報（脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点等）に基づいて、組織カラーマップに対して、幾何学的処理を施し例えば関心領域としての前頭葉を設定する。設定された関心領域は、当該関心領域のみを切り出して表示する形態、当該関心領域と他の領域とを色別表示する形態、又は当該関心領域を組織別に色別表示する形態等によって表示され、医師等に提供される。
【選択図】図１

Description

【０００１】
【発明の属する技術分野】
本発明は、例えば磁気共鳴映像装置（Magnetic Resonance Imaging device）等の医用画像診断装置に組み込まれる、又は医療用ワークステーション等によって提供される画像処理装置等に関する。
【０００２】
【従来の技術】
医用画像は、被検体についての多くの臨床情報を画像により提供するものであり、疾病の診断、治療や手術計画等を初めとする多くの医療行為において重要な役割を果たしている。医用画像は、医療用画像機器（例えば、超音波診断装置、Ｘ線ＣＴ装置、核医学診断装置、磁気共鳴診断装置等）によって収集されたデータに、当該医療用画像機器に内蔵される、又は医療用ワークステーション等によって実現される画像処理装置により、所定の画像処理を施すことで生成される。
【０００３】
この医用画像を利用した診断では、被検体の所定部位に関する画像において、診断部位を含む関心領域（Region of Interest ：以下「ＲＯＩ」とも称する。）に注目する。医用画像を利用した脳萎縮の診断を例にとれば、医師は、一般に頭部ＭＲ画像のＴ２強調水平断を対象とし、側脳室前角が現れるスライス画像によって画像診断を行う。医用画像における関心領域の抽出法は、診断部位によって様々である。例えば、上記脳萎縮を診断する場合では、医師がＭＲ画像上に現れる脳の解剖学的構造を目視で確認しながら、マニュアル操作により、脳萎縮の診断に必要な関心領域を選択している。
【０００４】
しかしながら、脳構造は個人により微妙な差異があり、特に白質と灰白質との境界などは、脳の解剖図譜に基づいて作成された標準脳地図を参照し、関心領域をマニュアル操作にて抽出するのは困難である。特に、前頭葉萎縮の程度を定量的に解析するためには、個人の脳構造に合わせて注目する関心領域を選択する必要がある。しかし、従来のマニュアル操作による抽出法では、その抽出根拠は医師毎に異なるため、客観的な定量的評価を行うことはできない。この事情は、診断部位が脳である場合に限らず、客観的基準に基づく関心領域の抽出法が確立されていない診断部位についても同様である。
【０００５】
【発明が解決しようとする課題】
本発明は、上記事情を鑑みてなされたもので、例えば医用画像において、診断等に必要な関心領域を解剖学的基準に基づいて自動抽出する画像処理装置、画像処理方法、及び当該画像処理をコンピュータに実行させるプログラムを格納する記録媒体を提供することを目的としている。
【０００６】
【課題を解決するための手段】
本発明は、上記目的を達成するため、次のような手段を講じている。
【０００７】
本発明の第１の視点は、被検体の所定部位の輪郭線に関する第１の画像と、前記所定部位の組織分布を表す第２の画像と、を記憶する記憶手段と、前記第１の画像及び前記第２の画像から得られる解剖学的構造情報に基づいて、第３の画像に幾何学的処理を施し、当該第３の画像に診断対象を含む関心領域を設定する関心領域設定手段と、前記関心領域を所定の形態にて表示する表示手段とを具備することを特徴とする画像処理装置である。
【０００８】
本発明の第２の視点は、第１の視点に係る装置において、前記解剖学的構造情報は、脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点のうちの少なくとも一つを含むこと、を特徴とするものである。
【０００９】
本発明の第３の視点は、第１又は第２の視点に係る装置において、前記第１の画像は、学習近傍領域を可変とした自己組織化マップによる輪郭抽出処理によって取得された画像であり、前記第２の画像は、自己組織化マップによる組織の自動分類処理によって取得された画像であることを特徴とするものである。
【００１０】
本発明の第４の視点は、第１乃至３のうちいずれかの視点に係る装置において、前記表示手段は、前記関心領域のみを切り出して表示する第１の形態、前記関心領域と他の領域とを色別表示する第２の形態、又は前記関心領域を組織別に色別表示する第３の形態の何れかによって前記関心領域を表示することを特徴とするものである。
【００１１】
本発明の第５の視点は、被検体の所定部位の輪郭線に関する第１の画像と、前記所定部位の組織分布を表す第２の画像と、から得られる解剖学的構造情報に基づいて、第３の画像に幾何学的処理を施し、当該第３の画像に診断対象を含む関心領域を設定する関心領域設定ステップと、前記関心領域を所定の形態にて表示する表示ステップと、を具備することを特徴とする画像処理方法である。
【００１２】
本発明の第６の視点は、第５の視点に係る方法において、前記所定部位は脳であり、前記解剖学的構造情報は、脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点のうちの少なくとも一つを含むこと、を特徴とするものである。
【００１３】
本発明の第７の視点は、第５又は第６の視点に係る方法において、前記第１の画像は、学習近傍領域を可変とした自己組織化マップによる輪郭抽出処理によって取得された画像であり、前記第２の画像は、自己組織化マップによる組織の自動分類処理によって取得された画像であること、を特徴とするものである。
【００１４】
本発明の第８の視点は、第５乃至７のうちいずれかの視点に係る方法において、前記所定の形態は、前記関心領域のみを切り出して表示する第１の形態、前記関心領域と他の領域とを色別表示する第２の形態、又は前記関心領域を組織別に色別表示する第３の形態の何れかであることを特徴とするものである。
【００１５】
本発明の第９の視点は、コンピュータに、被検体の所定部位の輪郭線に関する第１の画像と、前記所定部位の組織分布を表す第２の画像と、から得られる解剖学的構造情報に基づいて、第３の画像に幾何学的処理を施し、当該第３の画像に診断対象を含む関心領域を設定する関心領域設定ステップと、前記関心領域を所定の形態にて表示する表示ステップと、を実行させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体である。
【００１６】
本発明の第１０の視点は、第９の視点に係る記憶媒体において、前記所定部位は脳であり、
前記解剖学的構造情報は、脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点のうちの少なくとも一つを含むこと、を特徴とするものである。
【００１７】
本発明の第１１の視点は、第９又は第１０の視点に係る記憶媒体において、前記第１の画像は、学習近傍領域を可変とした自己組織化マップによる輪郭抽出処理によって取得された画像であり、前記第２の画像は、自己組織化マップによる組織の自動分類処理によって取得された画像であること、を特徴とするものである。
【００１８】
本発明の第１２の視点は、第９乃至１１のうちいずれかの視点に係る記憶媒体方法において、前記所定の形態は、前記関心領域のみを切り出して表示する第１の形態、前記関心領域と他の領域とを色別表示する第２の形態、又は前記関心領域を組織別に色別表示する第３の形態の何れかであることを特徴とするものである。
【００１９】
このような構成によれば、例えば医用画像において、診断等に必要な関心領域を解剖学的基準に基づいて自動抽出する画像処理装置、画像処理方法、及び当該画像処理をコンピュータに実行させるプログラムを格納する記録媒体として実現することができる。
【００２０】
【発明の実施の形態】
以下、本実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００２１】
本実施形態に係る画像処理装置は、医用画像を用いた画像診断において、被検体の所定部位を含む画像を対象に、解剖学的構造情報に基づいて、診断部位を含む関心領域の自動設定を行うものである。以下においては、説明を具体的にするため、脳萎縮の診断において、所定の被検体に関する頭部ＭＲ画像のＴ_２強調水平断において、診断対象である前頭葉を関心領域とし、これを解剖学的構造情報に基づいて画像上に自動設定する場合を例とする。
【００２２】
しかしながら、本発明の技術的思想は、これに限定されない。例えば、組織分類の対象とする画像は、Ｘ線ＣＴ画像等の他の断層画像でもよい。また、本実施形態で採用する解剖学的構造情報、及びこれに基づく幾何学的処理は一例であり、他の解剖学的構造情報及び幾何学的手法を採用してもよい。さらに、本技術的思想は、脳以外の診断においても有益であることは、言うまでもない。
【００２３】
図１は、本実施形態に係る画像処理装置１０のブロック構成図を示している。同図に示すように、画像処理装置１０は、ＭＲＩデータベース１１、組織カラーマップ記憶部１２、構造情報抽出フェーズ１３、関心領域設定フェーズ１５、関心領域記憶部１７、インタフェース部１９、脳表面輪郭線画像記憶部２１を具備している。
【００２４】
組織カラーマップ記憶部１２は、診断に供するＭＲ画像に対して、組織の自動分類法を適用することで得られる組織カラーマップを記憶している。ここで、自動分類法とは、診断に供するＭＲ画像の画像特性を自己組織化マップ（Self-Organizing Map：以下「ＳＯＭ」と称する。）により自己写像することで、画像を組織毎に分類するものである。この自動分類法については、後で詳しく説明する。
【００２５】
ＭＲＩデータベース１１は、被検体（患者）毎のＭＲ画像を格納している。
【００２６】
脳表面輪郭線画像記憶部２１は、脳の組織カラーマップに対して、輪郭抽出処理を行うことで得られる脳表面輪郭線画像を記憶している。ここで、輪郭抽出処理とは、脳のおおまかな輪郭線上の点の位置をＳＯＭにより自己写像することで、より正確で客観性の高い脳表面の輪郭線を生成するものである。この輪郭抽出処理についても、後で詳しく説明する。
【００２７】
構造情報抽出フェーズ１３は、処理対象領域設定部１３０、大脳縦裂連結線抽出部１３１、脳脊髄液細線化部１３２、側脳室端点抽出部１３３、外側溝端点抽出部１３４を有している。処理対象領域設定部１３０は、組織カラーマップ記憶部１２が記憶する組織カラーマップと、脳表面輪郭線画像記憶部２１が記憶する脳表面輪郭線画像とに基づいて、処理対象領域を設定する。この処理対象領域は、組織カラーマップに対して設定される。しかしながら、これに限定する趣旨ではなく、例えば当該組織カラーマップ作成の原画像（今の場合、ＭＲ画像）に対して処理対象領域を設定し、これに対して関心領域の自動抽出を行う構成であってもよい。大脳縦裂連結線抽出部１３１は、脳表面輪郭線から上大脳縦裂特徴点と下大脳縦裂特徴点とを抽出し、双方を線分で結んで大脳縦裂連結線を定義する。脳脊髄液細線化部１３２は、大脳縦裂連結線の二等分線を脳の仮想中心とし、当該仮想中心を通る水平線分より上に存在する脳脊髄液を細線化する。側脳室端点抽出部１３３は、上記仮想中心から所定の半径Ｒ内に存在する線分を追跡し、解剖学的構造情報と類似する端点を抽出して側脳室端点と定義する。
【００２８】
外側溝端点抽出部１３４は、上記側脳室端点からユークリッド距離が最小になる脳脊髄液細線を検索してその一端を抽出し、外側溝端点と定義する。
【００２９】
関心領域設定フェーズ１５は、脳表面接点抽出部１５０、大脳縦裂連結線接点抽出部１５１、関心領域設定・表示部１５２を有している。脳表面接点抽出部１５０は、外側溝端点から脳表面輪郭線に水平線を引き、その接点を抽出する。大脳縦裂連結線接点抽出部１５１は、側脳室端点から大脳縦裂連結線に垂線を引き、その接点を抽出する。関心領域設定・表示部１５２は、解剖学的構造情報に基づく一連の幾何学的処理によって得られた情報に基づいて、関心領域を設定し、これをインタフェース部１９に出力する。
【００３０】
関心領域記憶部１７は、取得された関心領域の位置に関する情報や当該関心領域のみを切り出した画像等を記憶する。
【００３１】
インタフェース部１９は、表示部、操作部等からなり、設定された関心領域を表示する。
【００３２】
なお、図１においては、ＭＲＩデータベース１１、組織カラーマップ記憶部１２、脳表面輪郭線画像記憶部２１をそれぞれ個別に表現した。しかしながら、単一のメモリによって各構成を実現する構成であってもよい。また、ＭＲＩデータベース１１等を本画像処理装置１０と別体として設け、ネットワークを介してＭＲ画像等を取得する構成であってもよい。
【００３３】
（自己組織化マップ：Self-Organizing Map ）
次に、自己組織化マップ（ＳＯＭ）について説明する。ＳＯＭは、入力データの位相を保存し、トポロジカルなマッピングを行うものであり、その学習過程は、明示的な教師を必要としない。
【００３４】
一般のＳＯＭは、図２に示すように、入力ユニットを含む入力層と、出力ユニットを含むマッピング層との二層からなる。典型的なＳＯＭの学習アルゴリズムは、次の様である。
【００３５】
（１）ｗ_ｉｊ（１≦ｉ≦ｎ）を時刻ｔにおける入力ユニットｉから出力ユニットｊへの重み係数とする。ユニットの重み係数を乱数で初期化し、ノードｊの近傍の初期範囲を大きく設定する。
【００３６】
（２）ｘ_ｉ（１≦ｉ≦ｎ）を時刻ｔにおけるノードｉへの入力とする。
【００３７】
（３）入力データと出力ノードｊとのユークリッド距離ｄ_ｊを次の式（１）にて計算する。
【００３８】
【数１】

【００３９】
（４）ユークリッド距離ｄ_ｊが最小となる出力ユニットを検索する。
【００４０】
（５）Ｎ_ｃ（ｔ）で定義される近傍に含まれるユニットへの重み係数を、次の式（２）にて更新する。
ｗ_ｉｊ（ｔ＋１）＝ｗ_ｉｊ（ｔ）＋α（ｔ）（ｘ_ｉ（ｔ）−ｗ_ｉｊ（ｔ））（２）
ここで、α（ｔ）は学習率係数（０＜α＜１）、Ｎ_ｃ（ｔ）は近傍領域のサイズであり、時間とともに減少させる。
【００４１】
（６）上記（２）〜（５）の処理を繰り返す。
【００４２】
次に、上記ＳＯＭ処理を利用した組織の自動分類法による、組織カラーマップの作成処理について説明する。
【００４３】
すなわち、図３に示すように、組織分類の対象となるＭＲ画像（以下、「分類対象画像」と称する。）に対して、ランダムに抽出される複数の局所ブロックに関する４つの特徴量（注目画素の輝度値、局所ブロックを構成する各画素の輝度値の平均値、注目画素の輝度値と局所ブロックを構成する画素の輝度値の最大値との差、注目画素の輝度値と局所ブロックを構成する画素の輝度値の最小値との差）を入力層における入力とし、所定の重み係数によって患者の所定部位の組織分布を分類対象画像の特徴空間における画像の位相特性として自己写像する。この自己写像により、入力としての各局所ブロックは、ＳＯＭ処理により背景（ＢＧ）、脳脊髄液（ＣＳＦ）、白質（ＷＭ）、灰白質（ＧＭ）の４つの出力ユニットＡ、Ｂ、Ｃ及びＤのいずれかに写像される。
【００４４】
すなわち、いまｘ１＝注目画素の輝度値、ｘ２＝局所ブロックを構成する各画素の輝度値の平均値、ｘ３＝注目画素の輝度値と局所ブロックを構成する画素の輝度値の最大値との差、ｘ４＝注目画素の輝度値と局所ブロックを構成する画素の輝度値の最小値との差、とする。各局所ブロックにおける各特徴量ｘ_ｉ（ただし、ｉは１≦ｉ≦４を満たす整数）を入力し、式（１）に基づいて、初期設定の重み係数ｗ_ｉｊによってユークリッド距離ｄ_ｊを算出する。なお、今の場合、ｊ＝１が出力ユニットＣＳＦに対応、ｊ＝２が出力ユニットＧＭに対応、ｊ＝３が出力ユニットＷＭ、ｊ＝４が出力ユニットＢＧに対応するものとして、ｊは１≦ｊ≦４を満たす整数とする。
【００４５】
入力は、算出されたユークリッド距離ｄ_１、ｄ_２、ｄ_３、ｄ_４のうち、最小のものに対応する出力ユニットに写像される。なお、当該自己写像においてユークリッド距離が最小になる出力ユニットは、勝者ユニットと呼ばれる。
【００４６】
また、この抽出された局所ブロックに対するＳＯＭ処理は、式（２）によって逐次重み係数を更新しながら、複数回（例えば１００００００回）実行される。この様に、重み係数を更新しながら、ＳＯＭ処理を複数回繰り返し学習する。学習が進むに従って、各局所ブロックは、所定の出力ユニットに収束するようになる。
【００４７】
このＳＯＭ処理の結果と、分類対象画像において生成される、各画素を注目画素とする複数の局所ブロック（３×３画素）の４つの特徴量とに基づいて、分類対象画像の各画素を背景（ＢＧ）、脳脊髄液（ＣＳＦ）、白質（ＷＭ）、灰白質（ＧＭ）のいずれかに組織分類し、各画素に所定の色を割り当てることで、組織カラーマップを生成することができる。
【００４８】
次に、学習近傍領域を可変としたＳＯＭ処理による輪郭抽出処理について説明する。輪郭抽出処理とは、当該患者の脳の輪郭線を代表する所定の代表点を学習データとするＳＯＭ処理を繰り返し実行することで、より客観性の高い代表点を取得するものである。本実施形態においては、初期輪郭線と上矢状静脈洞のエッジ部分とを合成して生成される脳の大まかな輪郭線を構成する複数の点（画素）の一部又は全部を、学習データとしての代表点とする。
【００４９】
また、以下に説明する例は、二種類のＳＯＭ処理を実行しその結果を合成することで、対象の複雑さに応じた輪郭抽出を実現するものである。二種類のＳＯＭ処理の一方は「初期近傍領域が３／４のＳＯＭ処理」と呼ばれ、他方は「初期近傍領域が１／８のＳＯＭ処理」と呼ばれる。
【００５０】
初期近傍領域が３／４のＳＯＭ処理では、図５に示すように、学習データ上の所定の各画素に関するＸ座標及びＹ座標を入力層における入力とし、所定の各重み係数によって、マッピング層における各出力ユニットに写像し、その中から最適な出力ユニット（勝者ユニット）を選択するＳＯＭ処理が繰り返し実行される。その際、各ＳＯＭ処理においては、学習近傍領域内に存在する出力ユニットに対応する各重み係数の更新が実行される。ここで、学習近傍領域とは、得られた勝者ユニットを少なくとも含む領域であり、学習回数の進行に従って収束するものである。また、初期近傍領域とは、学習近傍領域の初期設定を意味する。本ＳＯＭ処理では、初期近傍領域を、出力ユニット総数の３／４として設定し、所定の直線又は曲線に従って収束させる。
【００５１】
定量的な例により、上記初期近傍領域が３／４のＳＯＭ処理の内容をより詳しく説明すると、次のようである。すなわち、学習データ上の所定画素の座標（ｘ_ｉ，ｙ_ｉ）を入力とし、当該所定画素に対して式（１）に基づいてユークリッド距離ｄ_ｊを算出する。ここで、ｉは１≦ｉ≦最大学習回数（１０００００回）を満たす整数、ｊは１≦ｊ≦マッピング層の出力ユニット総数（１０００ユニット）を満たす整数である。算出されたユークリッド距離ｄ_ｊが最小となるユニットを勝者ユニットとし、勝者ユニット及び当該勝者ユニットを中心とした近傍領域内の出力ユニット（近傍ユニット）に対して式（２）に基づいて重み係数を更新する。今の場合、初期近傍領域が３／４であるから、当該領域内に存在する出力ユニットの総数は７５０ユニットとなる。
【００５２】
このように学習近傍領域を収束させつつ、当該領域内に存在する出力ユニットに対応する重み係数を式（２）によって逐次更新しながら、学習データ上の各画素についてのＳＯＭ処理を複数回（例えば１０００００回）実行する。最終的に、学習データ上の各画素の座標は、それぞれ所定の出力ユニットに収束するようになる。
【００５３】
また、初期近傍領域が１／８のＳＯＭ処理では、図５に示すように、初期近傍領域を最初に出力ユニット総数の１／８として設定し、所定の直線又は曲線に従って収束させながら、同様の処理をおこなう。
【００５４】
なお、本実施形態では、図５に示すように、各ＳＯＭ処理における複数の出力ユニットは、第１出力ユニット（最初の出力ユニット）と第１０００出力ユニット（最後の出力ユニット）とが連結された閉空間を構成している。これにより、第１０００出力ユニットは、第１出力ユニットの近傍ユニットとして判断され、式（２）による重み係数の更新に影響することになる。この様に閉空間の出力ユニットを採用するのは、抽出しようとする輪郭が閉曲線であるからである。これにより、第１出力ユニットと第１０００出力ユニットとの相関度合いが高くなり、より精度の高い結果を得ることができる。
【００５５】
また、処理の収束性の観点から、学習回数が進行するに従って、各ＳＯＭ処理における近傍領域のサイズ及び学習率係数を線形で減少させることが好ましい。
【００５６】
次に、上記のように構成した画像処理装置１０によって実行される関心領域の自動設定処理について、図６、図７を参照しながら説明する。なお、以下においては、説明を解りやすくするために、主に脳の右半球に関する処理を例として説明するが、脳の左半球に関する処理も同様の内容である。また、右半球に関する処理と左半球に関する処理とは、並行して実行されるものとする。
【００５７】
図６、図７は、画像処理装置１０が実行する関心領域の自動設定処理の流れを示したフローチャートである。図６において、まず、前段の処理として、既述の自動分類法により組織カラーマップを作成し、組織カラーマップ記憶部１２に格納する（ステップＳ１）。同様に、前段の処理として、既述の輪郭抽出法により、脳表面の輪郭を抽出する（ステップＳ２）。
【００５８】
次に、処理対象領域設定部１３０は、脳表面輪郭線内に存在する組織カラーマップを処理対象領域として設定する。大脳縦裂連結線抽出部１３１は、脳表面輪郭線に基づいて、図８に示すように、上記処理対象領域に上大脳縦裂特徴点ａと下大脳縦裂特徴点ｂとを抽出し、双方を線分Ｌで連結して大脳縦裂連結線を定義する（ステップＳ３）。
【００５９】
次に、脳脊髄液細線化部１３２は、図８に示すように大脳縦裂連結線の２等分点を脳の仮想中心ｓと定義し、当該仮想中心ｓを通る水平線分より上の脳脊髄液を細線化する（ステップＳ４）。以下、細線化された脳脊髄液を「脳脊髄液細線」と呼ぶ。
【００６０】
次に、側脳室端点抽出部１３３は、上記仮想中心ｓから半径Ｒ（ここでは、Ｒ＝線分ａｓの１／３）内に存在する脳脊髄液細線の線分を追跡し、図８に示すように解剖学的構造情報と類似する端点を側脳室端点ｃとする（ステップＳ５）。
【００６１】
次に、外側溝端点抽出部１３４は、上記側脳室端点ｃからユークリッド距離が最小になる右半球での脳脊髄液細線の端点を検索し（ステップＳ６、ステップＳ７）、検出した場合には、当該検出点を右半球での外側溝端点（単に、「右端点」とも言う。）として設定する（ステップＳ８）。一方、ユークリッド距離が最小になる右半球での脳脊髄液細線の端点を検出しない場合には、側脳室細線化端点の座標を擬似的な右半球での外側溝端点（単に、「右疑似端点」とも言う。）として設定する（ステップＳ８´）。
【００６２】
同様に、外側溝端点抽出部１３４は、上記側脳室端点ｃからユークリッド距離が最小になる左半球での脳脊髄液細線の端点を検索し（ステップＳ９）、検出した場合には、当該検出点を左半球での外側溝端点（単に、「左端点」とも言う。）として設定する（ステップＳ１０）。一方、ユークリッド距離が最小になる左半球での脳脊髄液細線の端点を検出しない場合には、側脳室細線化端点の座標を擬似的な左半球での外側溝端点（単に、「左疑似端点」とも言う。）として設定する（ステップＳ１０´）。
【００６３】
次に、外側溝端点抽出部１３４は、ステップＳ７、Ｓ９での判断の結果に基づいて、左右の外側溝端点の最終決定を行う。すなわち、外側溝端点抽出部１３４は、ステップＳ７において右端点有りと判断し、且つステップＳ９において左端点有りと判断した場合には、各端点をそれぞれの半球での外側溝端点ｄと決定する（ステップ１２）。また、外側溝端点抽出部１３４は、ステップＳ７において右端点無しと判断し、且つステップＳ９において左端点無しと判断した場合には、各疑似端点を各半球での外側溝端点ｄと決定する。さらに、外側溝端点抽出部１３４は、ステップＳ７、Ｓ９のいずれか一方においてのみ右端点又は左端点有りと判断した場合には、検出した一方の端点座標を大脳縦裂連結線に対して対象移動し、検出できなかった半球での外側溝端点ｄとして決定する（ステップＳ１２´）。
【００６４】
次に、脳表面接点抽出部１５０は、図８に示すように、外側溝端点ｄから脳表面輪郭に水平線を引き、その足を点ｅと定義する（ステップＳ１３）。また、大脳縦裂連結線接点抽出部１５１は、図８に示すように、側脳室端点ｃから大脳縦裂連結線に垂線を引き、その足を点ｆと定義する（ステップＳ１４）。関心領域設定・表示部１５２は、図８に示すように点ａ，ｆ，ｃ，ｄ，ｅを結んだ線分内を右前頭葉関心領域として設定する。こうして設定された関心領域の位置情報は、関心領域記憶部１７に記憶される。また、設定された関心領域は、所定の形態にてインタフェース１９の表示部に表示される。以下、関心領域の表示形態について、幾つかの例を挙げる。
【００６５】
図９（ａ）は頭部ＭＲ画像を示している。図９（ｂ）は、図９（ａ）に示した頭部ＭＲ画像から切り出された左前頭葉関心領域を示している。図９（ｃ）は、図９（ａ）に示した頭部ＭＲ画像から切り出された右前頭葉関心領域を示している。図９（ｂ）、（ｃ）に示すように、関心領域以外をマスクした形態にて表示することで、観察しやすい診断画像を提供することができる。
【００６６】
また、図１０（ａ）、（ｂ）、（ｃ）、（ｄ）は、頭部ＭＲ画像に設定された関心領域のカラー表示を示している。すなわち、図１０（ａ）では、左右前頭葉関心領域の脳脊髄液（ＣＳＦ）のみをカラー表示している。図１０（ｂ）では、左右前頭葉関心領域の白質のみをカラー表示している。図１０（ｃ）では、左右前頭葉関心領域の灰白質のみをカラー表示している。図１０（ｄ）では、左右前頭葉関心領域の脳脊髄液（ＣＳＦ）、白質、灰白質を色別表示し、前頭葉萎縮が判断しやすいように表示している。これらのカラー表示は、組織カラーマップを利用することで実現することができる。
【００６７】
以上述べた構成によれば、以下の効果を得ることができる。
【００６８】
本画像処理装置においては、解剖学的構造情報（例えば、外側溝、側脳室前角等）に基づいて、例えば前頭葉等の関心領域が自動的に抽出される。従って、操作者の主観を必要とせず、より客観性の高い関心領域の抽出を迅速かつ簡便に実施することができる。その結果、質の高い画像診断に寄与することができる。
【００６９】
また、本画像処理装置は、個々の患者毎の解剖学的構造情報を利用して、前頭葉等の関心領域を自動抽出している。従って、対象とする画像の種類、対象となる被検体が変わっても、同様により客観性の高い関心領域の抽出を迅速かつ簡便に実施することができる。さらに、画像上に肉眼で解剖学的構造情報が確認できるか否かに関係なく、関心領域を自動抽出することができる。その結果、質の高い画像診断に寄与することができる。
【００７０】
また、本画像処理装置は、抽出した関心領域を種々の形態にて表示することができる。従って、医師等の観察者は、自身が観察しやすい所望の形態を自由に選択することができ、診断に有益な医用画像を使い勝手よく提供することができる。
【００７１】
以上、本発明を実施形態に基づき説明したが、本発明の思想の範疇において、当業者であれば、各種の変更例及び修正例に想到し得るものであり、それら変形例及び修正例についても本発明の範囲に属するものと了解される。例えば以下に示すように、その要旨を変更しない範囲で種々変形可能である。
【００７２】
本実施形態で説明した関心領域の自動設定処理は、当該方法をコンピュータに実行させるプログラムによって、パーソナルコンピュータやワークステーション等によっても実現することができる。この場合、コンピュータに当該方法を実行させることのできるプログラムは、磁気ディスク（フロッピーディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することができる。
【００７３】
また、各実施形態は可能な限り適宜組み合わせて実施してもよく、その場合組合わせた効果が得られる。さらに、上記実施形態には種々の段階の発明が含まれており、開示される複数の構成要件における適宜な組合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題が解決でき、発明の効果の欄で述べられている効果の少なくとも１つが得られる場合には、この構成要件が削除された構成が発明として抽出され得る。
【００７４】
【発明の効果】
以上本発明によれば、例えば医用画像において、診断等に必要な関心領域を解剖学的基準に基づいて自動抽出する画像処理装置、画像処理方法、及び当該画像処理をコンピュータに実行させるプログラムを格納する記録媒体を実現できる。
【図面の簡単な説明】
【図１】図１は、本実施形態に係る画像処理装置１０のブロック構成図を示している。
【図２】図２は、自己組織化マップの概念を説明するための図である。
【図３】図３は、組織の自動分類法に使用される局所ブロックを説明するための図である。
【図４】図４は、組織の自動分類法において実行される自己組織化マップの概念を説明するための図である。
【図５】図５は、輪郭抽出法において実行される自己組織化マップの概念を説明するための図である。
【図６】図６は、画像処理装置１０が実行する関心領域の自動設定処理の流れを示したフローチャートである。
【図７】図７は、画像処理装置１０が実行する関心領域の自動設定処理の流れを示したフローチャートである。
【図８】図８は、関心領域の自動設定処理を説明するための図である。
【図９】図９（ａ）は頭部ＭＲ画像を示している。図９（ｂ）は、図９（ａ）に示した頭部ＭＲ画像から切り出された左前頭葉関心領域を示している。図９（ｃ）は、図９（ａ）に示した頭部ＭＲ画像から切り出された右前頭葉関心領域を示している。
【図１０】図１０（ａ）、（ｂ）、（ｃ）、（ｄ）は、頭部ＭＲ画像に設定された関心領域のカラー表示を示している。
【符号の説明】
１０…画像処理装置
１１…ＭＲＩデータベース
１２…組織カラーマップ記憶部
１３…構造情報抽出フェーズ
１５…関心領域設定フェーズ
１７…関心領域記憶部
１９…インタフェース部
２１…脳表面輪郭線画像記憶部
１３０…処理対象領域設定部
１３１…大脳縦裂連結線抽出部
１３２…脳脊髄液細線化部
１３３…側脳室端点抽出部
１３４…外側溝端点抽出部
１５０…脳表面接点抽出部
１５１…大脳縦裂連結線接点抽出部
１５２…関心領域設定・表示部

Claims

被検体の所定部位の輪郭線に関する第１の画像と、前記所定部位の組織分布を表す第２の画像と、を記憶する記憶手段と、
前記第１の画像及び前記第２の画像から得られる解剖学的構造情報に基づいて、第３の画像に幾何学的処理を施し、当該第３の画像に診断対象を含む関心領域を設定する関心領域設定手段と、
前記関心領域を所定の形態にて表示する表示手段と、
を具備することを特徴とする画像処理装置。
前記所定部位は脳であり、
前記解剖学的構造情報は、脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点のうちの少なくとも一つを含むこと、
を特徴とする請求項１記載の画像処理装置。
前記第１の画像は、学習近傍領域を可変とした自己組織化マップによる輪郭抽出処理によって取得された画像であり、
前記第２の画像は、自己組織化マップによる組織の自動分類処理によって取得された画像であること、
を特徴とする請求項１又は２記載の画像処理装置。
前記表示手段は、前記関心領域のみを切り出して表示する第１の形態、前記関心領域と他の領域とを色別表示する第２の形態、又は前記関心領域を組織別に色別表示する第３の形態の何れかによって前記関心領域を表示することを特徴とする請求項１乃至３のうちいずれか一項記載の画像処理装置。
被検体の所定部位の輪郭線に関する第１の画像と、前記所定部位の組織分布を表す第２の画像と、から得られる解剖学的構造情報に基づいて、第３の画像に幾何学的処理を施し、当該第３の画像に診断対象を含む関心領域を設定する関心領域設定ステップと、
前記関心領域を所定の形態にて表示する表示ステップと、
を具備することを特徴とする画像処理方法。
前記所定部位は脳であり、
前記解剖学的構造情報は、脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点のうちの少なくとも一つを含むこと、
を特徴とする請求項５記載の画像処理方法。
前記第１の画像は、学習近傍領域を可変とした自己組織化マップによる輪郭抽出処理によって取得された画像であり、
前記第２の画像は、自己組織化マップによる組織の自動分類処理によって取得された画像であること、
を特徴とする請求項５又は６記載の画像処理方法。
前記所定の形態は、前記関心領域のみを切り出して表示する第１の形態、前記関心領域と他の領域とを色別表示する第２の形態、又は前記関心領域を組織別に色別表示する第３の形態の何れかであることを特徴とする請求項５乃至７のうちいずれか一項記載の画像処理方法。
コンピュータに、
被検体の所定部位の輪郭線に関する第１の画像と、前記所定部位の組織分布を表す第２の画像と、から得られる解剖学的構造情報に基づいて、第３の画像に幾何学的処理を施し、当該第３の画像に診断対象を含む関心領域を設定する関心領域設定ステップと、
前記関心領域を所定の形態にて表示する表示ステップと、
を実行させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
前記所定部位は脳であり、
前記解剖学的構造情報は、脳表面輪郭線、上大脳縦裂特徴点と下大脳縦裂特徴点とを連結する大脳縦裂連結線、脳脊髄液を細線化した脳脊髄液細線、側脳室端点、外側溝端点のうちの少なくとも一つを含むこと、
を特徴とする請求項９記載のコンピュータ読み取り可能な記録媒体。
前記第１の画像は、学習近傍領域を可変とした自己組織化マップによる輪郭抽出処理によって取得された画像であり、
前記第２の画像は、自己組織化マップによる組織の自動分類処理によって取得された画像であること、
を特徴とする請求項９又は１０記載のコンピュータ読み取り可能な記録媒体。
前記所定の形態は、前記関心領域のみを切り出して表示する第１の形態、前記関心領域と他の領域とを色別表示する第２の形態、又は前記関心領域を組織別に色別表示する第３の形態の何れかであることを特徴とする請求項９乃至１１のうちいずれか一項記載のコンピュータ読み取り可能な記録媒体。