JP5970766B2 - 医用画像処理装置、医用画像処理方法、プログラム - Google Patents

医用画像処理装置、医用画像処理方法、プログラム Download PDF

Info

Publication number
JP5970766B2
JP5970766B2 JP2011208436A JP2011208436A JP5970766B2 JP 5970766 B2 JP5970766 B2 JP 5970766B2 JP 2011208436 A JP2011208436 A JP 2011208436A JP 2011208436 A JP2011208436 A JP 2011208436A JP 5970766 B2 JP5970766 B2 JP 5970766B2
Authority
JP
Japan
Prior art keywords
image
division
tissue
voxel
normal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2011208436A
Other languages
English (en)
Other versions
JP2013066632A (ja
JP2013066632A5 (ja
Inventor
小野　徹太郎
徹太郎 小野
智章 後藤
智章 後藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dai Nippon Printing Co Ltd
Original Assignee
Dai Nippon Printing Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dai Nippon Printing Co Ltd filed Critical Dai Nippon Printing Co Ltd
Priority to JP2011208436A priority Critical patent/JP5970766B2/ja
Priority to PCT/JP2012/073877 priority patent/WO2013047278A1/ja
Priority to CN201280046986.5A priority patent/CN103826536B/zh
Priority to US14/347,266 priority patent/US9390509B2/en
Priority to EP12837541.7A priority patent/EP2762072B1/en
Publication of JP2013066632A publication Critical patent/JP2013066632A/ja
Publication of JP2013066632A5 publication Critical patent/JP2013066632A5/ja
Application granted granted Critical
Publication of JP5970766B2 publication Critical patent/JP5970766B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30016—Brain

Description

本発明は、ＭＲＩ（Magnetic Resonance Imaging）によって撮像された脳画像を入力し、画像処理を行うことにより、疾患別の診断支援を行う医用画像処理装置等に関するものである。
本発明者らは、早期ＡＤ（Alzheimer’s Disease）診断支援システムであるＶＳＲＡＤ（登録商標）（Voxel-Based Specific Regional Analysis System for Alzheimer’s
Disease）を開発している。ＶＳＲＡＤ（登録商標）は、前駆期を含む早期アルツハイマー型認知症に特徴的にみられる海馬傍回付近（＝関心領域）の萎縮の程度をＭＲＩ画像から読み取るための画像処理・統計解析ソフトウエアである。
ＶＳＲＡＤ（登録商標）では、関心領域の萎縮の程度を自動解析し、ＭＲＩ脳画像にＺスコアマップを重畳表示する。Ｚスコアとは、被験者画像と健常者平均画像を統計比較した結果、平均値からどれだけの標準偏差分離れているかを示す値である。Ｚスコアマップは、サーモグラフィのように、Ｚスコアの分布をＺスコアの値に応じた色によって示し、被験者の脳画像に重ねて表示したものである。診断者は、Ｚスコアマップによって、萎縮の程度を視覚的に確認することができる。
ＶＳＲＡＤ（登録商標）では、ＶＢＭ（Voxel Based
Morphometry）（特許文献１参照）やＤＡＲＴＥＬ（Diffeomorphic Anatomical
Registration using Exponentiated Lie algebra）（特許文献２参照）といった技術を用いている。
ＶＢＭは、被験者の頭部を撮像して取得された脳画像を、３次元の画素であるボクセルを単位に画像処理を行う技術である。ＤＡＲＴＥＬは、それまでのＶＢＭに比べて、空間的標準化の処理の精度が優れており、アルツハイマー病の画像統計解析による診断能を向上させる技術として期待されている。ここで、空間的標準化の処理とは、個人の間に存在する脳画像の解剖学的な違いを吸収するために、脳全体の大きさに対する大局的な補正と、部分的な大きさに対する局所的な補正を行うものである。
ＭＲＩ脳画像（特に、Ｔ１強調ＭＲＩ脳画像）では、神経細胞に対応する灰色の灰白質、それより明るい神経線維に対応する白質、黒に近い脳脊髄液の３種類の組織が含まれている。ＤＡＲＴＥＬを用いることによって、被験者のＭＲＩ脳画像から組織分割により抽出される白質の異常である病変や萎縮の程度を正確に評価することが可能となる（特許文献２参照）。
特開２００５−２３７４４１号公報 ＷＯ２０１１／０４０４７３号公報
ところで、組織分割の処理は、様々な理由に起因して、灰白質組織と白質組織とが正確に分割されないことがある（以下、このような場合を「分割エラー」と称する。）。旧来のＶＢＭを用いる場合、Ｚスコアマップを確認することによって、分割エラーが生じていることを容易に判別することが出来た。例えば、Ｚスコアマップにおいて、萎縮の程度が進んでいる箇所を赤色とする場合、分割エラーが生じていると、海馬傍回付近以外の領域も含めて、脳全体が赤色になる。つまり、分割エラーに起因してＺスコアが異常な値として算出されたことが分かり、ひいては、分割エラーが生じていることを容易に判別することが出来た。
一方、ＤＡＴＥＬを用いる場合、旧来のＶＢＭに比べると、空間的標準化の処理が非常に強力になっている為、分割エラーが生じていても、それらしいＺスコアマップが表示されてしまう。つまり、分割エラーが生じていても、海馬傍回付近以外の領域も含めて、脳全体が赤色にならない場合がある。従って、ＤＡＴＥＬを用いる場合、Ｚスコアマップを確認しても、分割エラーが生じていることを判別することが難しい。そして、診断者が分割エラーに気付かずに診断を行ってしまうと、重大な誤診を招いてしまうという可能性があった。そこで、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知する仕組みが望まれている。
本発明は、前述した問題点に鑑みてなされたもので、その目的とすることは、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知することが可能な医用画像処理装置等を提供することである。
前述した目的を達成するために第１の発明は、被験者の脳画像を入力する入力手段と、前記被験者の脳画像に対して組織分割処理を行い、灰白質組織を分割して組織分割処理結果画像を作成する分割手段と、健常者の灰白質組織の画像群を保存する記憶手段と、前記分割手段により得られた前記被験者の前記組織分割処理結果画像と前記健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力する出力手段と、前記被験者の脳画像に係る前記組織分割処理結果画像のボクセル値と、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル毎の統計値と、に基づいて、前記分割手段による分割結果が正常か否かを識別する識別手段と、を具備し、前記組織分割正常画像群のボクセル毎の統計値は、前記組織分割正常画像群の各画像の対応する位置ごとのボクセル値の分布を、当該位置毎に平均値と標準偏差とで表したものであって、前記識別手段は、前記被験者の脳画像に係る前記組織分割処理結果画像のあるボクセルのボクセル値と、前記組織分割正常画像群の前記ボクセルに対応する位置の前記平均値との差に基づく値を、前記組織分割正常画像群の前記ボクセルに対応する位置の前記標準偏差で除することから定められる絶対Ｚスコアの前記組織分割処理結果画像全体に渡る平均値を、予め定義された閾値とを比較することによって、前記分割手段による分割結果が正常か否かを識別することを特徴とする医用画像処理装置である。第１の発明によって、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知することが可能となる。特に、組織分割処理による分割結果を自動的に識別するので、診断者には全く負荷がかからない。
また、これによって、組織分割処理による分割結果を精度良く識別することができる。
第１の発明は、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル値、又は対応する位置のボクセル毎の統計値、を記憶する記憶手段と、前記識別手段によって正常と識別された前記被験者に係る組織分割処理結果画像のボクセル値を前記記憶手段に反映する反映手段と、を更に具備することが望ましい。これによって、評価値（絶対Ｚスコア）の算出処理の精度が向上し、ひいては、識別処理の精度が向上する。
第１の発明は、前記識別手段によって正常ではないと識別された場合、警告情報を出力する警告手段、を更に具備することが望ましい。これによって、診断者は、分割エラーに気付かずに診断を行ってしまうことがない。
第１の発明は、前記識別手段によって正常ではないと識別された場合、前記分割手段による分割結果を表示する表示手段、を更に具備することが望ましい。これによって、診断者は、煩わしさを感じることがない。
第２の発明は、コンピュータを利用した医用画像処理方法であって、被験者の脳画像を入力する入力ステップと、前記被験者の脳画像に対して組織分割処理を行い、灰白質組織を分割して組織分割処理結果画像を作成する分割ステップと、前記分割ステップにより得られた前記被験者の前記組織分割処理結果画像と記憶手段に保存されている健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力する出力ステップと、前記被験者の脳画像に係る前記組織分割処理結果画像のボクセル値と、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル毎の統計値と、に基づいて、前記分割ステップによる分割結果が正常か否かを識別する識別ステップと、を含み、前記組織分割正常画像群のボクセル毎の統計値は、前記組織分割正常画像群の各画像の対応する位置ごとのボクセル値の分布を、当該位置毎に平均値と標準偏差とで表したものであって、前記識別ステップは、前記被験者の脳画像に係る前記組織分割処理結果画像のあるボクセルのボクセル値と、前記組織分割正常画像群の前記ボクセルに対応する位置の前記平均値との差に基づく値を、前記組織分割正常画像群の前記ボクセルに対応する位置の前記標準偏差で除することから定められる絶対Ｚスコアの前記組織分割処理結果画像全体に渡る平均値を、予め定義された閾値とを比較することによって、前記分割手段による分割結果が正常か否かを識別することを特徴とする医用画像処理方法である。第２の発明によって、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知することが可能となる。特に、組織分割処理による分割結果を自動的に識別するので、診断者には全く負荷がかからない。
第３の発明は、コンピュータを、被験者の脳画像を入力する入力手段と、前記被験者の脳画像に対して組織分割処理を行い、灰白質組織を分割して組織分割処理結果画像を作成する分割手段と、健常者の灰白質組織の画像群を保存する記憶手段と、前記分割手段により得られた前記被験者の前記組織分割処理結果画像と前記健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力する出力手段と、前記被験者の脳画像に係る前記組織分割処理結果画像のボクセル値と、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル毎の統計値と、に基づいて、前記分割手段による分割結果が正常か否かを識別する識別手段と、して機能させ、前記組織分割正常画像群のボクセル毎の統計値は、前記組織分割正常画像群の各画像の対応する位置ごとのボクセル値の分布を、当該位置毎に平均値と標準偏差とで表したものであって、前記識別手段は、前記被験者の脳画像に係る前記組織分割処理結果画像のあるボクセルのボクセル値と、前記組織分割正常画像群の前記ボクセルに対応する位置の前記平均値との差に基づく値を、前記組織分割正常画像群の前記ボクセルに対応する位置の前記標準偏差で除することから定められる絶対Ｚスコアの前記組織分割処理結果画像全体に渡る平均値を、予め定義された閾値とを比較することによって、前記分割手段による分割結果が正常か否かを識別するように機能させるためのプログラムである。第３の発明のプログラムを汎用のコンピュータにインストールすることによって、第１の発明の医用画像処理装置を得ることができるとともに、第２の発明の医用画像処理方法を実行することができる。
本発明により、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知することが可能な医用画像処理装置等を提供することができる。
医用画像処理装置の概要を示すブロック図 医用画像処理方法の事前処理を示すフローチャート 医用画像処理方法の診断支援情報作成処理を示すフローチャート 組織分割結果の正常例及び分割エラー例 医用画像処理方法の組織分割結果識別処理を示すフローチャート 分割結果の態様ごとの評価値の度数分布と閾値を示す図 分割エラーの識別性能を示すグラフ 閾値が０．８の場合の感度及び特異度を示す図
以下図面に基づいて、本発明の実施形態を詳細に説明する。図１は、医用画像処理装置の概要を示すブロック図である。医用画像処理装置１は、コンピュータに、本発明の実施の形態における医用画像処理方法を実行させるためのプログラムがインストールされることによって構築される。最初に、コンピュータのハードウエア構成について説明する。
コンピュータは、制御部、記憶部、表示部、入力部、Ｉ／Ｏ（Input／Output interface）等を備える。制御部は、ＣＰＵ（Central
Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）等から構成される。ＣＰＵは、記憶部、ＲＯＭ、記録媒体等に格納されるプログラムをＲＡＭ上のワークメモリ領域に呼び出して実行し、各装置を駆動制御し、コンピュータが行う後述する処理を実現する。ＲＯＭは、不揮発性メモリであり、コンピュータのブートプログラムやＢＩＯＳ等のプログラム、データ等を恒久的に保持している。ＲＡＭは、揮発性メモリであり、記憶部、ＲＯＭ、記録媒体等からロードしたプログラム、データ等を一時的に保持するとともに、制御部２１が各種処理を行う為に使用するワークエリアを備える。
記憶部は、例えば、ＨＤＤ（Hard Disk Drive）であり、制御部が実行するプログラム、プログラム実行に必要なデータ、ＯＳ（Operating System）等が格納される。プログラムに関しては、ＯＳに相当する制御プログラムや、後述する処理をコンピュータに実行させるためのアプリケーションプログラムが格納されている。これらの各プログラムコードは、制御部により必要に応じて読み出されてＲＡＭに移され、ＣＰＵに読み出されて各種の手段として実行される。
表示部は、液晶パネル等のディスプレイ装置、ディスプレイ装置と連携してコンピュータのビデオ機能を実現するための論理回路等（ビデオアダプタ等）を有する。入力部は、データの入力を行い、例えば、キーボード、マウス等のポインティングデバイス、テンキー等の入力装置を有する。入力部を介して、コンピュータに対して、操作指示、動作指示、データ入力等を行うことができる。尚、表示部及び入力部は、タッチパネルディスプレイのように、一体となっていても良い。
Ｉ／Ｏは、コンピュータに周辺機器（例えば、プリンタ、ネットワーク機器等）を接続させるためのポート等である。コンピュータは、Ｉ／Ｏを介して周辺機器とのデータの送受信を行う。また、コンピュータは、ネットワーク機器等を介して、医用画像の撮像装置（例えば、ＭＲＩ等）や医用画像管理サーバ等とのデータの送受信を行う。
次に、図１を参照しながら、医用画像処理装置１の機能について説明する。図１に示すように、医用画像処理装置１は、ユーザインタフェース部１０、処理部２０、及びデータベース部３０を備えている。
ユーザインタフェース部１０は、入力画像としてＭＲＩ画像を入力する画像入力機能１１と、処理部２０によって処理された結果を表示する結果出力機能１２とを有する。
処理部２０は、ユーザインタフェース部１０から入力されたＭＲＩ画像を処理する画像処理機能２１と、各種統計演算処理等を行う統計処理機能２２と、画像処理機能２１による結果が正常か否かを識別する識別処理機能２３とを有する。
データベース部３０には、処理部２０による後述する処理に使用する白質脳画像テンプレート３１、白質脳画像テンプレート３２、健常者画像データベース３３、疾患特異的ＲＯＩ３４、及び組織分割正常画像データベース３５等が保存されている。
本実施形態における医用画像処理装置１は、（１）被験者のＭＲＩ脳画像を基に診断支援情報の出力等を行う診断支援情報出力処理、及び、（２）診断支援情報出力処理に含まれる画像処理の結果が正常か否かを識別し、画像処理の結果が正常ではない場合には警告情報の出力等を行う画像処理結果識別処理を行うことによって、診断者による診断を支援する。
以下では、診断支援情報出力処理の１例として、ＷＯ２０１１／０４０４７３号公報（特許文献２）に記載の第１実施形態を説明する。本実施の形態ではフローチャートのみを図に示し、処理の流れを説明する。
図２は、医用画像処理方法の事前処理を示すフローチャートである。図２の事前処理では、図３の診断支援情報作成処理におけるステップ１３の空間的標準化において使用するテンプレートを作成する。
図２に示すように、医用画像処理装置１は、出来るだけ多くの健常者に係るＴ１強調ＭＲＩ脳画像（図２では被験者画像）を入力する。
各被験者から取得したＭＲＩ脳画像については、前処理を行っておく。具体的には、被験者の脳全体を含むように所定の厚さのスライス状に撮像した、例えば１００〜２００枚のＴ１強調ＭＲＩ画像を入力とする。また、各スライス画像におけるボクセルの各辺の長さを予め等しくなるようにスライス画像のリサンプリングを行っておく。ここで、ボクセルは、「厚さ」を持つ画像の座標単位であり、２次元画像におけるピクセルに相当する。
このような前処理を行ったＭＲＩ脳画像を入力した後、そのスライス画像の撮像方向や解像度が、予めシステムに設定されている条件に適合しているか否かをチェックする。
以上のように、ＭＲＩ脳画像が設定条件に適合して入力されていることが確認された場合には、医用画像処理装置１は、空間的位置合わせ処理を行う（ステップ１１）。
これは、入力された脳画像を、標準的な脳画像と比較する際の精度を上げるために、アフィン変換によって空間的位置と角度の補正を行っていることに相当する。
以上の空間的位置合わせが終了した後、医用画像処理装置１は、組織分割処理を行い（ステップ１２）、白質組織と灰白質組織とを各々抽出し、白質画像と灰白質画像を作成する。尚、本発明の実施の形態では、組織分割処理が白質組織と灰白質組織を分割するものとして説明するが、灰白質組織を正しく分割できれば、組織分割処理により得られた被験者の灰白質組織の画像と健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力することが可能である。従って、組織分割処理は、少なくとも灰白質組織を分割するものであれば十分である。
入力されたＴ１強調ＭＲＩ脳画像には、神経線維に対応する高い信号値を呈する白質、神経細胞に対応する中間の信号値を呈する灰白質、低い信号値を呈する脳脊髄液の３種類の組織が含まれている。そこで、医用画像処理装置１は、信号値の差に着目して白質組織と灰白質組織とを各々抽出する処理を行う。この処理については、特開２００５−２３７４４１号公報（特許文献１）に説明されている。本実施の形態では、特許文献１において行っている手法よりも抽出精度が高い統合型組織分割処理を行う。統合型組織分割処理は、標準化、組織分割、信号不均一の補正を１つのモデルに組み込んだ組織分割手法である。詳しくは、「J.Ashburner and K. J. Friston. Unified segmentation. Neuro Image.
2005; 26: 839-851.」に記載されている。統合型組織分割処理では、白質画像および灰白質画像の他に、ＭＲＩ画像の座標と標準脳の座標の対応関係を示す変換場が作成されるという特徴がある。変換場は後述する標準化において使用する。
このように、医用画像処理装置１は、予め多くの健常者のＭＲＩ脳画像から組織分割により、白質、灰白質が３次元的に抽出された白質画像、灰白質画像をそれぞれ多数のサンプルとして取得する。
以上のように、医用画像処理装置１は、多数の健常者のＭＲＩ脳画像をそれぞれ組織分割して白質が抽出された白質画像をサンプルとして作成し、作成された全てのサンプル間における空間的標準化により白質テンプレートを作成する（ステップ１３）。同様に、医用画像処理装置１は、多数の健常者のＭＲＩ脳画像をそれぞれ組織分割して灰白質が抽出された灰白質画像をサンプルとして作成し、全てのサンプル間における空間的標準化により灰白質テンプレートを作成する。
ここで実行される空間的標準化には、ダーテル（ＤＡＲＴＥＬ）・アルゴリズムが適用される。前述した通り、ＤＡＲＴＥＬは、従来のＶＢＭと比較して空間的標準化精度に優れており、アルツハイマー病の画像統計解析による診断能を向上させる技術として期待されている。また、ＤＡＲＴＥＬの空間的標準化は従来の手法よりも精密に行われるため、灰白質だけではなく、白質をも評価対象とすることができる。
ステップ１３のテンプレート作成処理では、医用画像処理装置１は、白質と灰白質それぞれについて、年齢や性別などの被験者の属性に応じた層別のテンプレートを作成し、データベース部３０に白質脳画像テンプレート３１、灰白質脳画像テンプレート３２として保存する。
以上のようにして作成された白質と灰白質のテンプレートが、年齢別、性別毎に用意されていることを前提として、医用画像処理装置１は、図３に示す診断支援情報作成処理を行う。尚、以下では、白質と灰白質のテンプレートのことをダーテルテンプレートと呼ぶことにする。
図３は、医用画像処理方法の診断支援情報作成処理を示すフローチャートである。図３に示すように、診断支援情報作成処理では、診断対象の被験者に係るＴ１強調ＭＲＩ脳画像を入力データとする。
図２の事前処理と同様、医用画像処理装置１は、被験者画像を入力し、各スライス画像におけるボクセルの各辺の長さを予め等しくなるようにスライス画像のリサンプリングを行う。そして、医用画像処理装置１は、空間的位置合わせの処理を、事前処理におけるステップ１１の場合と同様に行う（ステップ２１）。
以上の空間的位置合わせが終了した後、医用画像処理装置１は、組織分割処理を行う（ステップ２２）。この組織分割処理は、事前処理におけるステップ１２の場合と同様である。医用画像処理装置１は、白質と灰白質を抽出し、被験者の白質画像と灰白質画像をそれぞれ作成する。
以上のように作成された被験者の白質画像と灰白質画像に対して、医用画像処理装置１は、空間的標準化の処理を行う（ステップ２３）。ここで行う空間的標準化には、事前処理におけるステップ１３の場合と同様に、ダーテル・アルゴリズムが適用される。
この空間的標準化の処理は、個人の間に存在する脳画像の解剖学的な違いを吸収するために、脳全体の大きさに対する大局的な補正と、部分的な大きさに対する局所的な補正を行うものである。以下、便宜上、灰白質を中心に説明するが、白質の場合も実質的に同一の処理を行う。
ステップ２３のダーテルでの空間的標準化の処理は、次の３つのステップの処理によって構成されている。
（ステップ２３−１）初期位置決定処理
（ステップ２３−２）ダーテルテンプレートへの変換処理
（ステップ２３−３）標準脳テンプレートへの変換処理
ステップ２３−１の初期位置決定処理では、医用画像処理装置１は、前述した統合型組織分割処理によって得られる標準脳への変換場を使用して、灰白質画像、白質画像に対して、初期位置を決める処理を行う。この処理では、剛体変換を行うので、画像の形状は変わらないのが特徴である。
ステップ２３−２のダーテルテンプレートへの変換処理では、医用画像処理装置１は、ステップ２３−１が実行された画像に対して、ダーテル・アルゴリズムを用いて、ダーテルテンプレートに形状を合わせ込む。
ステップ２３−２の標準脳テンプレートへの変換処理では、医用画像処理装置１は、ステップ２３−２においてダーテルテンプレートに合わせ込まれた画像を、標準脳テンプレートに合わせ込む処理を行う。ダーテルテンプレートから標準脳テンプレートへの変換場を予め求めておき、その変換場を用いて標準脳座標系への変換を行う。
ステップ２３−２とステップ２３−３の処理では、各ボクセルの信号値の合計を保持したまま標準化を行うことによって、体積の情報が保持される為、標準化後に体積を計測することが可能となる。
ステップ２３−１では線形変換、ステップ２３−２及びステップ２３−３では線形変換と非線形変換が行われる。ステップ２３−２を例として説明すると、医用画像処理装置１は、線形変換と非線形変換を用いて、データベース部３０から読み出したステップ１２において作成された平均的な灰白質脳画像テンプレート３２との誤差の平方和が最小になるように画像処理を行う。この空間的標準化処理では、初めに線形変換による位置や大きさ、角度の大局的な補正を行い、次に非線形変換によって局所的な凹凸等の形状の補正を行う。
ここで行われる線形変換は、ステップ１１の位置合わせと同様のアフィン変換である。また、非線形変換は、ｘ方向、ｙ方向それぞれについて局所的な変位を表す変形場を推定し、この変形場によって元画像の変換を行うものである。
ステップ２３−２の処理は、ステップ１３において作成されたダーテルテンプレートを雛型として、処理対象画像を合わせ込む処理である。使用するテンプレートは、ダーテル・アルゴリズムを適用して高精度に作成されていることから、その形状がシャープである。
そのため、空間的標準化によって個々の処理対象が、個体差のない近い形に合わせ込まれるようになることから、個体間での形状は同一に近くなるが、萎縮などは局所的な密度に反映されるようになるため、空間的標準化の精度を向上できることになる。
以上のように空間的標準化を施した白質画像と灰白質画像（以下、「標準化脳画像」とも言う。）に対して、医用画像処理装置１は、画像平滑化の処理を行う（ステップ２４）。
これは、標準化脳画像のＳ／Ｎ比を向上させると共に、後に比較を行う際に標準として使用する健常者の画像群と画像のsmoothnessを等しくするための処理であり、３次元ガウシアンカーネルを使用して行う。この平滑化に使用するフィルタのＦＷＨＭ（半値幅）は８ｍｍ程度とする。
具体的には、特許文献１に説明されているように、医用画像処理装置１は、３次元的脳画像と、３次元ガウシアン関数の３次元的な畳み込み（コンボリュージョン）を行う。これは、ｘ、ｙ、ｚ各方向における１次元の畳み込みを逐次的に行うことに相当する。このように平滑化を行うことにより、空間的標準化処理において完全に一致しない個体差を低減させることができる。
ステップ２３の処理では、脳の体積の情報が保存されている。そのため、次の濃度値補正を実行する前に、白質及び灰白質の処理結果画像について、全体又は後述する関心領域（ＲＯＩ：Regions Of Interest）の積分値を体積として計測し、診断支援情報として活用しても良い。
以上のように画像平滑化を行った標準化脳画像に対して、後に比較を行う際に標準として使用する健常者の画像群におけるボクセル値の分布に合わせるために、脳全体のボクセル値を補正する濃度値補正を行う。
その後、医用画像処理装置１は、健常者群との統計的比較を行う（ステップ２５）。ここでは、以上のステップ２１〜２４の各処理を通して標準化を行った被験者の灰白質（白質）のＭＲＩ脳画像と、予め収集してデータベース部３０に健常者画像データベース３３として保存してある健常者の灰白質（白質）のＭＲＩ脳画像群との比較検定を行う。使用する健常者画像群は、被験者の年齢に近いものによって構成されていることが望ましい。
具体的には、医用画像処理装置１は、健常者画像群とボクセル単位で１：Ｎ（Ｎは健常者画像の総数）の比較検定を行い、統計的に有意な差が見られる（異常と推定される）ボクセルを検出する。
まず、全てのボクセルについて、それぞれ次式で表されるＺスコアを算出する。
Ｚ＝（μ―ｘ）／σ・・・（１）
ここで、Ｚ：Ｚスコア、μ：健常者画像群の対応するボクセル値の平均、ｘ：被験者画像のボクセル値、σ：健常者画像群の対応するボクセル値の標準偏差、である。
このように、Ｚスコアは、被験者画像のボクセル値と、健常者画像群の対応するボクセルのボクセル値平均との差を、標準偏差によってスケーリングした値であり、これは灰白質（白質）容積の相対的低下の度合を示すものである。
次に、適当な臨界値Ｚ’を定め、Ｚスコアが、
Ｚ’＜Ｚ・・・（２）
となるようなボクセルを求め、統計的に有意な差が見られるボクセルとする。臨界値には、約９５％以上の確率で異常と推定できるＺ’＝２を用いる。また、健常者よりも容積が低下している部位全てを含む臨界値の指定法として、次式も用いる。
０＜Ｚ・・・（３）
なお、ステップ２５において使用した健常者画像データベース３３は、予め別途収集した健常者の画像群のそれぞれに対して、ステップ２１〜２４の空間的位置合わせ→灰白質（白質）の組織分割処理→空間的標準化→画像平滑化等の各処理を順次行って同様に作成し、保存してあるものである。
また、この医用画像処理装置１においては、収集したこれらの健常者画像を、例えば５歳毎又は１０歳毎というように年代別に分類し、それぞれの群について算出した平均値と標準偏差を記憶部に保存しておくことにより、Ｚスコアによる検定を行うことができる。
なお、このようにＺスコアを使用する場合には、ボクセル毎に平均値と標準偏差のデータだけを持っていれば良いので、データ作成後は画像データ自体を保存しておく必要がないという利点もある。
以上のように、被験者の標準化脳画像に対して統計的比較を行った後、医用画像処理装置１は、ＲＯＩ分析を行う（ステップ２６）。
これは、脳画像を用いて異常の有無を判別する場合に、画像上に所定の大きさの関心領域（ＲＯＩ）を設定する方法であり、脳画像上において、特定の疾患に関係するとして注目されている特定部位に所定の大きさのＲＯＩを設定して比較を行うものである。
この解析方法は、特許文献１に説明されているように、統計処理により健常者と有意な差が見られた座標位置のボクセルとそのＺスコアについて、疾患に対応するＲＯＩ（疾患特異的ＲＯＩ３４）を適用することにより、罹患している度合を求めるものである。その特徴は次の２点である。
（特徴１）アルツハイマー病の疾患毎に対応する標準化された画像データとしてのＲＯＩ（疾患特異的ＲＯＩ３４）を用意しておき、被験者の症状から考えられる疾患について、被験者の脳画像データにそれぞれのＲＯＩを適用（設定）し、このＲＯＩにおけるＺスコアに基づいて最も優位性の高いものを診断結果とする。
（特徴２）ＲＯＩの部分のみのＺスコアによって疾患を判断するだけでなく、ＲＯＩを適用しない場合の脳全体のＺスコアマップと、ＲＯＩを適用した部分のみのＺスコアマップとの比較を行う。この目的は、脳全体の萎縮に対する注目部位の萎縮の度合を見ることにある。
ここでは、疾患Ａ〜Ｃの疾患別の特異的ＲＯＩが用意されている場合を例として、被験者がある疾患Ａを罹患しているか否かを判別する方法を説明する。
ステップ２５の統計的比較によって得られた被験者のＺスコアマップに対して、疾患Ａに対応するＲＯＩを用いて、（２）式及び（３）式を使って、以下の５つのパラメータを算出する。
Ｐ１＝ＲＯＩ部分において式（３）を満たすボクセルのＺスコアの合計／ＲＯＩ部分において式（３）を満たすボクセルの数
Ｐ２＝脳全体において式（２）を満たすボクセルの数／脳全体のボクセル数
Ｐ３＝ＲＯＩ部分において式（２）を満たすボクセルの数／ＲＯＩ部分のボクセル数
Ｐ４＝Ｐ３／Ｐ２
Ｐ５＝ＲＯＩ部分における全てのボクセルの中で最大のＺスコア
Ｐ１〜Ｐ５の５つのパラメータについて、予め疾患Ａを有する患者群における特性を求めておき、被験者のパラメータの値がそれに合致する場合に、被験者は疾患Ａであると判別する。
このような判別結果は、診断支援情報として出力される。また、Ｚスコアや各パラメータの値も、診断支援情報として出力される。また、サーモグラフィのように、Ｚスコアの分布をＺスコアの値に応じた色によって示し、被験者の脳画像に重ねて表示するＺスコアマップも、診断支援情報として出力される。診断支援情報は、例えば、医用画像処理装置１の表示部に表示される。
本実施形態における医用画像処理装置１は、診断支援情報を出力する前に、画像処理結果識別処理（ステップ２７）を実行する。画像処理結果識別処理については、図５を参照しながら後述する。
以上の通り、医用画像処理装置１は、診断支援情報出力処理を行う。ところで、組織分割の処理は、様々な理由に起因して、灰白質組織と白質組織とが正確に分割されない分割エラーが生じることがある。ＤＡＲＴＥＬを適用した医用画像処理装置１の処理結果では、ほとんどの場合、診断者が目視によってＺスコアマップを確認しても、分割エラーの発生の有無を判別することができない。
図４は、組織分割結果の正常例及び分割エラー例である。ＩＤ１〜ＩＤ７は、axial（横断面）スライスの画像であり、いずれも灰白質の部位を白色として示し、それ以外の部位を黒色として示している。ＩＤ１が正常例、ＩＤ２〜ＩＤ７が分割エラー例である。分割エラーの原因は、撮像時に被験者が動いてしまったり、撮像位置がずれていたり、撮像条件が誤っていたりする等、様々である。
分割エラーが生じやすい画像としては、以下のような画像が挙げられる。
・灰白質と白質のコントラストが不十分な画像
・ＳＮ比（信号対雑音比）が好ましくない画像（ノイズが目立つ画像）
・信号ムラがある画像
・アーチファクト（磁化率、体動、折り返しなどに起因する）がある画像
・撮像範囲が広すぎる画像
・白質梗塞など組織の変性が広い領域で起きている画像
ＩＤ１の画像は、灰白質と白質の分割が正常に行われており、灰白質が正しく抽出されている。
ＩＤ２の画像は、白質が多く含まれている。これは、ＩＤ１の画像と比較して観察すると理解できる。ＩＤ２の画像では、組織分割処理において、本来は白質として抽出されるべき部分が、灰白質として抽出されてしまったものである。
ＩＤ３の画像は、脳脊髄液が多く含まれている。これは、ＩＤ１の画像と比較して観察すると理解できる。ＩＤ３の画像では、組織分割処理において、本来は脳脊髄液として抽出されるべき部分が、灰白質として抽出されてしまったものである。
ＩＤ４の画像は、灰白質がほとんど出ていない。これは、ＩＤ１の画像と比較して観察すると理解できる。ＩＤ４の画像では、組織分割処理において、本来は灰白質として抽出されるべき部分が、白質や脳脊髄液として抽出されてしまったものである。
ＩＤ５の画像は、全体の形状・位置の乖離が大きい。これは、ＩＤ１の画像と比較して観察すると理解できる。ＩＤ５の画像では、撮像位置が正しくなかったものと思われる。
ＩＤ６の画像は、ノイズ、ぼやけなどがひどい。これは、ＩＤ１の画像と比較して観察すると理解できる。ＩＤ６の画像では、撮像条件が正しくなかったものと思われる。
ＩＤ７の画像は、広範囲で欠けがある。これは、ＩＤ１の画像と比較して観察すると理解できる。ＩＤ７の画像では、撮影範囲が正しくなかったか、又は一部の撮像データに欠損があったものと思われる。
本実施形態における医用画像処理装置１では、図５に示す画像処理結果識別処理によって、図４に示すような分割結果に対して、正常か否かを自動的に識別し、分割エラーが生じている或いは分割エラーが生じている可能性があれば、その旨を診断者に通知する仕組みを提供する。
図５は、医用画像処理方法の組織分割結果識別処理を示すフローチャートである。組織分割結果識別処理は、診断支援情報を出力する前に実行される（図３のステップ２７参照）。
図５に示すように、医用画像処理装置１は、診断支援情報作成処理の対象となった被験者（以下、「対象被験者」とも言う。）の灰白質画像を特定する（ステップ３１）。
ステップ３１において特定される灰白質画像は、例えば、ＭＮＩ（Montreal
Neurological Institute）座標系の原点を通るaxial（横断面）スライスの画像である。この画像は、海馬傍回付近（本実施形態における関心領域）を横断する。但し、ステップ３１において特定される灰白質画像は、この例に限定されるわけではなく、他のスライスの画像であっても良い。また、ステップ３１において特定される灰白質画像は、１枚に限らず、複数枚であっても良い。
次に、医用画像処理装置１は、ステップ３１において特定された対象被験者の灰白質画像の平滑化を行う（ステップ３２）。ここでの平滑化は、図３のステップ２４における画像平滑化の処理と同様である。
ステップ３２における平滑化は、例えば、３次元ガウシアンカーネルを使用して行う。この平滑化に使用するフィルタのＦＷＨＭ（半値幅）は８ｍｍ程度とする。
ステップ３２における平滑化は、灰白質の微細構造を一定程度除去するためと、組織分割正常画像データベース３５に登録されている画像群と画像のsmoothnessを等しくするための処理である。組織分割正常画像データベース３５には、図３のステップ２２における組織分割処理が正常に実行された灰白質画像（以下、「組織分割正常画像」とも言う。）が登録されている。医用画像処理装置１は、組織分割正常画像に対して、ステップ３２において対象被験者の灰白質画像に行った平滑化と同様の平滑化を行った後、組織分割正常画像データベース３５に登録する。
次に、医用画像処理装置１は、ステップ３２において平滑化された画像に対して、予め定義された評価関数に従って、評価値を算出する（ステップ３３）。評価関数は、以下に示す例が挙げられる。
但し、ｘ軸及びｙ軸は、axial（横断面）を構成する軸である。ｚ軸は、被験者の体軸、すなわちaxial（横断面）と直交する軸である。
input(x,y)は、ｚ軸の位置を固定し、対象被験者の灰白質画像におけるボクセル(x,y)のボクセル値である。input(x,y,z)は、対象被験者の灰白質画像におけるボクセル(x,y,z)のボクセル値である。
mean(x,y)は、ｚ軸の位置を固定し、組織分割正常画像データベース３５に登録されている画像群のボクセル(x,y)毎の平均値である。mean(x,y,z)は、組織分割正常画像データベース３５に登録されている画像群のボクセル(x,y,z)毎の平均値である。
SD(x,y)は、ｚ軸の位置を固定し、組織分割正常画像データベース３５に登録されている画像群のボクセル(x,y)毎の標準偏差である。mean(x,y,z)は、組織分割正常画像データベース３５に登録されている画像群のボクセル(x,y,z)毎の標準偏差である。
Ｘは、x軸方向のボクセル数である。Ｙは、ｙ軸方向のボクセル数である。Ｚは、ｚ軸方向のボクセル数である。
式（４）〜式（９）に示す評価関数は、いずれも、式（１）に示すＺスコアの数式に基づくものであり、かつ、値域が０以上である。そこで、式（４）〜式（９）に示す評価関数、及び、式（４）〜式（９）において分母＝１とした評価関数によって算出される評価値を、「絶対Ｚスコア」と総称する。尚、評価値は、絶対Ｚスコアに限定されず、例えば、前述したＰ１〜Ｐ５に基づく様々な変形例も考えられる。
式（４）、（６）及び（８）に基づく評価関数を用いる場合、医用画像処理装置１は、ステップ３１において、対象被験者の灰白質画像を１スライス特定する。一方、式（５）、（７）及び（９）に基づく評価関数を用いる場合、医用画像処理装置１は、ステップ３１において、Ｚ（ｚ軸方向のボクセル数）の値と一致する枚数を、対象被験者の灰白質画像として特定する。尚、対象被験者の灰白質画像を１枚だけ利用する場合であっても、図６〜図８を参照しながら後述するように、分割結果を高精度に識別することが可能である。
このように、本実施形態では、組織分割正常画像そのものを利用するのではない。そこで、組織分割正常画像データベース３５には、各組織分割正常画像のボクセル値を記憶することに代えて、組織分割正常画像群のボクセル毎の統計値（例えば、mean(x,y)、mean(x,y,z)、SD(x,y)、SD(x,y,z)等）を記憶するようにしても良い。後者の場合、組織分割正常画像そのものは記憶されないので、各被験者の医療情報を確実に保護することができる。
次に、医用画像処理装置１は、ステップ３３において算出された評価値と、予め定義された閾値とを比較し、分割結果を識別する（ステップ３４）。閾値の定義法は、図６〜図８を参照しながら後述する。
医用画像処理装置１は、分割結果が正常と識別した場合（ステップ３５の「正常」）、ユーザによってこの識別結果が正しいと確認され、入力部を介して登録する旨の指示がなされたら（ステップ３６の「はい」）、対象被験者の灰白質画像（但し、ステップ３２において平滑化された後の画像）のボクセル値を、組織分割正常画像データベース３５に登録する（ステップ３７）。つまり、医用画像処理装置１は、ステップ３４において正常と識別された被験者の脳画像に係る灰白質組織のボクセル値を組織分割正常画像データベース３５に反映する。一方、ユーザによってこの識別結果が正しくないと確認され、入力部を介して登録しない旨の指示がなされたら（ステップ３６の「いいえ」）、医用画像処理装置１は、処理を終了する。
尚、組織分割正常画像データベース３５に、組織分割正常画像群のボクセル毎の統計値を記憶するように構成する場合、医用画像処理装置１は、対象被験者の灰白質組織のボクセル値を登録することに代えて、次のような処理を行う。すなわち、医用画像処理装置１は、既に記憶されている組織分割正常画像群のボクセル毎の統計値に対して、対象被験者の灰白質画像のボクセル値を反映させる。そして、医用画像処理装置１は、対象被験者の灰白質画像のボクセル値が反映された統計値を、組織分割正常画像データベース３５に登録する。
何れの場合（対象被験者の灰白質組織のボクセル値を登録する場合、及び対象被験者の灰白質画像のボクセル値が反映された統計値を登録する場合）であっても、評価値の算出処理に用いられる組織分割正常画像のボクセル毎の統計値は、より多くの画像群に基づいて自動的に更新されることになる。従って、評価値の算出処理の精度が向上し、ひいては、識別処理の精度が向上する。
一方、医用画像処理装置１は、分割結果が正常ではないと識別した場合（ステップ３５の「正常ではない」）、分割結果が正常ではない旨を警告する（ステップ３８）。警告手段としては、例えば、分割結果が正常ではない旨のメッセージを表示部に表示したり、音によって出力したりすることが考えられる。分割結果が正常ではない旨のメッセージの１例としては、「Check segment results!」（分割結果を確認して下さい！）や、「Probable
sgment error!」（分割エラーの可能性有り！）等が挙げられる。これらのメッセージを表示する画面は、Ｚスコアマップを含む処理結果表示画面や、診断支援情報のレポート出力画面等が挙げられる。
そして、医用画像処理装置１は、分割結果を示す分割結果表示画面を表示する（ステップ３９）。ここで、分割結果表示画面には、図４に示すような灰白質画像に加えて、分割前画像、白質画像、脳脊髄液画像などが表示される。
診断支援情報作成処理の全体を通して考えると、分割結果は、診断支援情報を作成するための中間結果に過ぎない。従って、診断者は、何も問題がなければ、分割結果を確認する手順を省略したい。そこで、医用画像処理装置１は、分割結果が正常ではないと識別した場合のみ、分割結果を示す分割結果表示画面を表示する。
以上の通り、医用画像処理装置１は、画像処理結果識別処理を行う。以下では、図６〜図８を参照しながら、図５における閾値の定義法を説明する。
図６は、分割結果の態様ごとの評価値の度数分布と閾値を示す図である。図７は、分割エラーの識別性能を示すグラフである。図８は、閾値が０．８の場合の感度及び特異度を示す図である。
図６〜図８は、図４に示す正常例及び分割エラー例に従う分割結果を、それぞれ複数用意し、閾値の定義を行ったものである。図８に示すように、各分割結果のデータ数は以下の通りである。
・識別ＩＤ＝１（分割結果は正常）の例：８０５件
・識別ＩＤ＝２（白質が多く含まれる）の例：４件
・識別ＩＤ＝３（脳脊髄液が多く含まれる）の例：１３件
・識別ＩＤ＝４（灰白質がほとんど出ない）の例：０件
・識別ＩＤ＝５（全体の形状・位置の乖離が大きい）の例：８件
・識別ＩＤ＝６（ノイズ・ぼやけがひどい）の例：４件
・識別ＩＤ＝７（広範囲で欠けがある）の例：０件
尚、識別ＩＤ＝４（灰白質がほとんど出ない）の例や識別ＩＤ＝７（広範囲で欠けがある）の例は、Ｚスコアマップを目視にて確認することによって、分割エラーが生じているか否かを容易に判別することができたことから、本実施形態では考慮に入れていない。
図６の度数分布図は、横軸が分割結果の態様を識別するＩＤ、縦軸が図５のステップ３３において算出される評価値である。図６に示す評価値は、式（６）の評価関数によって求めたものである。
図６では、分割結果が正常な例のプロットを実線の矩形によって囲み、分割結果がエラーの例のプロットを点線の矩形によって囲んでいる。図６では、閾値＝０．８によって識別し、閾値以上であれば識別結果をエラーとし、閾値未満であれば識別結果を正常としている。
図６において閾値を０．８と定義した理由は、図７によって説明する。図７の識別性能のグラフは、横軸が評価値（或いは閾値）、縦軸が感度・特異度である。
感度とは、ある検査について「陽性（本実施形態では「分割結果が正常ではない」）と判定されるべきものを正しく陽性と判定する確率」として定義される値である。感度が高い（高感度である）とは、「陽性と判定されるべきものを正しく陽性と判定する可能性が高い」、或いは「陽性と判定されるべきものを間違って陰性と判定する可能性が低い」という意味である。
特異度は、感度と対となる表現であり、図７に示すように、感度との兼ね合いによって決まる。特異度とは、ある検査について「陰性（本実施形態では「分割結果が正常」）と判定されるべきものを正しく陰性と判定する確率」として定義される値である。特異度が高い（高特異度である）とは、「陰性と判定されるべきものを正しく陰性と判定する可能性が高い」、或いは「陰性と判定されるべきものを間違って陽性と判定する可能性が低い」という意味である。
本実施形態では、診断者が分割エラーに気付かずに診断を行ってしまうと、重大な誤診を招いてしまうという可能性を排除するために、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知することが目的である。そこで、閾値は高感度となるように定義する、つまり、「分割結果が正常ではない」と判定されるべきものを正しく「分割結果が正常ではない」と判定するような値とすることが望ましい。
図７のグラフから、感度が１００％となる中で、最も高特異度となる値を閾値として定義すると、閾値＝０．８となる。
図８では、閾値＝０．８としたときの感度及び特異度を示している。閾値＝０．８であれば、識別ＩＤ＝２〜７の分割エラー例の全てを分割エラーと識別し、感度＝１００％となっている。また、識別ＩＤ＝１の正常例についても、ほとんどが正常と識別し、高特異度（特異度＝９９．５％）となっている。
以上、本実施の形態における医用画像処理装置１によって、組織分割の処理において分割エラーが生じている或いは分割エラーが生じている可能性があることを診断者に通知することが可能となる。
特に、医用画像処理装置１は、組織分割処理による分割結果を自動的に識別するので、診断者には全く負荷がかからない。
また、医用画像処理装置１は、組織分割処理による分割結果を正常ではないと識別した場合、診断支援情報とともに、分割エラーが生じている或いは分割エラーが生じている可能性がある旨の警告情報を出力するので、診断者は、分割エラーに気付かずに診断を行ってしまうことがない。
また、医用画像処理装置１は、組織分割処理による分割結果を正常ではないと識別した場合にのみ、分割結果を示す分割結果表示画面を表示するので、診断者は、煩わしさを感じることがない。
また、医用画像処理装置１は、組織分割処理が正常に行われた脳画像群に係る灰白質組織のボクセル値又はボクセル毎の統計値を組織分割正常画像データベース３４に記憶し、更に、正常と識別された対象被験者の脳画像に係る灰白質組織のボクセル値を組織分割正常画像データベース３４に反映するので、評価値の算出処理の精度が向上し、ひいては、識別処理の精度が向上する。
尚、前述の説明では、組織分割正常画像データベース３５には、被験者の属性や撮像条件等を考慮せずに全てを纏めて記憶するものとしたが、組織分割正常画像データベース３５には、被験者の属性や撮像条件毎に分けて記憶するようにしても良い。被験者の属性は、例えば、性別や年齢等である。撮像条件は、ＭＲＩ機種名、磁場強度、ボクセル数、ＦＯＶ（Field Of View）の大きさ、スライス厚、撮像時間、ＴＲ（Repetition
Time）、ＴＥ（Echo Time）、ＦＡ（Flip
Angle）等である。この場合、医用画像処理装置１は、ボクセル毎の統計値を、被験者の属性や撮像条件毎に分けて算出する。そして、医用画像処理装置１は、図５のステップ３３における評価値を算出する際、対象被験者の属性や、対象被験者の脳画像を撮像したときの撮像条件に応じたボクセル毎の統計値を評価関数に代入する。
以上、添付図面を参照しながら、本発明に係る医用画像処理装置等の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
１………医用画像処理装置
１０………ユーザインタフェース部
１１………画像入力機能
１２………結果出力機能
２０………処理部
２１………画像処理機能
２２………統計処理機能
２３………識別処理機能
３０………データベース部
３１………白質脳画像テンプレート
３２………灰白質脳画像テンプレート
３３………健常者画像データベース
３４………疾患特異的ＲＯＩ
３５………組織分割正常画像データベース

Claims (6)

  1. 被験者の脳画像を入力する入力手段と、
    前記被験者の脳画像に対して組織分割処理を行い、灰白質組織を分割して組織分割処理結果画像を作成する分割手段と、
    健常者の灰白質組織の画像群を保存する記憶手段と、
    前記分割手段により得られた前記被験者の前記組織分割処理結果画像と前記健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力する出力手段と、
    前記被験者の脳画像に係る前記組織分割処理結果画像のボクセル値と、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル毎の統計値と、に基づいて、前記分割手段による分割結果が正常か否かを識別する識別手段と、
    を具備し、
    前記組織分割正常画像群のボクセル毎の統計値は、前記組織分割正常画像群の各画像の対応する位置ごとのボクセル値の分布を、当該位置毎に平均値と標準偏差とで表したものであって、
    前記識別手段は、前記被験者の脳画像に係る前記組織分割処理結果画像のあるボクセルのボクセル値と、前記組織分割正常画像群の前記ボクセルに対応する位置の前記平均値との差に基づく値を、前記組織分割正常画像群の前記ボクセルに対応する位置の前記標準偏差で除することから定められる絶対Ｚスコアの前記組織分割処理結果画像全体に渡る平均値を、予め定義された閾値とを比較することによって、前記分割手段による分割結果が正常か否かを識別することを特徴とする医用画像処理装置。
  2. 前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル値、又は対応する位置のボクセル毎の統計値、を記憶する記憶手段と、
    前記識別手段によって正常と識別された前記被験者に係る組織分割処理結果画像のボクセル値を前記記憶手段に反映する反映手段と、
    を更に具備することを特徴とする請求項１に記載の医用画像処理装置。
  3. 前記識別手段によって正常ではないと識別された場合、警告情報を出力する警告手段、
    を更に具備することを特徴とする請求項１又は請求項２に記載の医用画像処理装置。
  4. 前記識別手段によって正常ではないと識別された場合、前記分割手段による分割結果を表示する表示手段、
    を更に具備することを特徴とする請求項１乃至請求項３のいずれかに記載の医用画像処理装置。
  5. コンピュータを利用した医用画像処理方法であって、
    被験者の脳画像を入力する入力ステップと、
    前記被験者の脳画像に対して組織分割処理を行い、灰白質組織を分割して組織分割処理結果画像を作成する分割ステップと、
    前記分割ステップにより得られた前記被験者の前記組織分割処理結果画像と記憶手段に保存されている健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力する出力ステップと、
    前記被験者の脳画像に係る前記組織分割処理結果画像のボクセル値と、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル毎の統計値と、に基づいて、前記分割ステップによる分割結果が正常か否かを識別する識別ステップと、
    を含み、
    前記組織分割正常画像群のボクセル毎の統計値は、前記組織分割正常画像群の各画像の対応する位置ごとのボクセル値の分布を、当該位置毎に平均値と標準偏差とで表したものであって、
    前記識別ステップは、前記被験者の脳画像に係る前記組織分割処理結果画像のあるボクセルのボクセル値と、前記組織分割正常画像群の前記ボクセルに対応する位置の前記平均値との差に基づく値を、前記組織分割正常画像群の前記ボクセルに対応する位置の前記標準偏差で除することから定められる絶対Ｚスコアの前記組織分割処理結果画像全体に渡る平均値を、予め定義された閾値とを比較することによって、前記分割手段による分割結果が正常か否かを識別することを特徴とする医用画像処理方法。
  6. コンピュータを、
    被験者の脳画像を入力する入力手段と、
    前記被験者の脳画像に対して組織分割処理を行い、灰白質組織を分割して組織分割処理結果画像を作成する分割手段と、
    健常者の灰白質組織の画像群を保存する記憶手段と、
    前記分割手段により得られた前記被験者の前記組織分割処理結果画像と前記健常者の灰白質組織の画像群との統計的比較に基づいて診断支援情報を出力する出力手段と、
    前記被験者の脳画像に係る前記組織分割処理結果画像のボクセル値と、前記組織分割処理が正常に行われた脳画像群である組織分割正常画像群のボクセル毎の統計値と、に基づいて、前記分割手段による分割結果が正常か否かを識別する識別手段と、
    して機能させ、
    前記組織分割正常画像群のボクセル毎の統計値は、前記組織分割正常画像群の各画像の対応する位置ごとのボクセル値の分布を、当該位置毎に平均値と標準偏差とで表したものであって、
    前記識別手段は、前記被験者の脳画像に係る前記組織分割処理結果画像のあるボクセルのボクセル値と、前記組織分割正常画像群の前記ボクセルに対応する位置の前記平均値との差に基づく値を、前記組織分割正常画像群の前記ボクセルに対応する位置の前記標準偏差で除することから定められる絶対Ｚスコアの前記組織分割処理結果画像全体に渡る平均値を、予め定義された閾値とを比較することによって、前記分割手段による分割結果が正常か否かを識別するように機能させるためのプログラム。
JP2011208436A 2011-09-26 2011-09-26 医用画像処理装置、医用画像処理方法、プログラム Active JP5970766B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2011208436A JP5970766B2 (ja) 2011-09-26 2011-09-26 医用画像処理装置、医用画像処理方法、プログラム
PCT/JP2012/073877 WO2013047278A1 (ja) 2011-09-26 2012-09-19 医用画像処理装置、医用画像処理方法、プログラム
CN201280046986.5A CN103826536B (zh) 2011-09-26 2012-09-19 医用图像处理装置、医用图像处理方法
US14/347,266 US9390509B2 (en) 2011-09-26 2012-09-19 Medical image processing device, medical image processing method, program
EP12837541.7A EP2762072B1 (en) 2011-09-26 2012-09-19 Medical image processing device, medical image processing method, program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011208436A JP5970766B2 (ja) 2011-09-26 2011-09-26 医用画像処理装置、医用画像処理方法、プログラム

Publications (3)

Publication Number Publication Date
JP2013066632A JP2013066632A (ja) 2013-04-18
JP2013066632A5 JP2013066632A5 (ja) 2014-08-21
JP5970766B2 true JP5970766B2 (ja) 2016-08-17

Family

ID=47995310

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011208436A Active JP5970766B2 (ja) 2011-09-26 2011-09-26 医用画像処理装置、医用画像処理方法、プログラム

Country Status (5)

Country Link
US (1) US9390509B2 (ja)
EP (1) EP2762072B1 (ja)
JP (1) JP5970766B2 (ja)
CN (1) CN103826536B (ja)
WO (1) WO2013047278A1 (ja)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11216498B2 (en) 2005-10-26 2022-01-04 Cortica, Ltd. System and method for generating signatures to three-dimensional multimedia data elements
US11361014B2 (en) 2005-10-26 2022-06-14 Cortica Ltd. System and method for completing a user profile
US9646138B2 (en) * 2011-12-30 2017-05-09 The Johns Hopkins University Bioimaging grid
JP6402437B2 (ja) * 2013-09-20 2018-10-10 大日本印刷株式会社 培地情報登録システム、検査システム、衛生管理システム、及び、プログラム
JP2016064004A (ja) * 2014-09-25 2016-04-28 大日本印刷株式会社 医用画像表示処理方法、医用画像表示処理装置およびプログラム
KR101665032B1 (ko) * 2014-11-13 2016-10-11 삼성전자 주식회사 자기 공명 영상 장치 및 자기 공명 영상의 처리 방법
CN105045220B (zh) * 2015-05-08 2018-03-23 上海质晟生物科技有限公司 一种用于实验室诊断领域或工业生产领域的基于多变量z分数质量控制图的质量控制方法
JP2018517207A (ja) * 2015-05-18 2018-06-28 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 自己認識画像セグメント化方法及びシステム
US11195043B2 (en) 2015-12-15 2021-12-07 Cortica, Ltd. System and method for determining common patterns in multimedia content elements based on key points
EP3248543B8 (en) 2016-05-24 2024-02-21 Siemens Healthineers AG Apparatus and method for visualizing tissue macro- and microstructure for pathology evaluation in magnetic resonance imaging
CN106897993B (zh) * 2017-01-12 2019-07-26 华东师范大学 基于定量磁化率成像人脑灰质核团概率图谱的构建方法
CN108877922A (zh) * 2017-05-15 2018-11-23 沈渊瑶 病变程度判断系统及其方法
US10140421B1 (en) 2017-05-25 2018-11-27 Enlitic, Inc. Medical scan annotator system
JP7140606B2 (ja) * 2018-08-29 2022-09-21 富士フイルムヘルスケア株式会社 画像処理装置、画像処理方法、画像処理プログラム及び磁気共鳴イメージング装置
US11700356B2 (en) 2018-10-26 2023-07-11 AutoBrains Technologies Ltd. Control transfer of a vehicle
US10748038B1 (en) 2019-03-31 2020-08-18 Cortica Ltd. Efficient calculation of a robust signature of a media unit
US11145059B2 (en) 2018-11-21 2021-10-12 Enlitic, Inc. Medical scan viewing system with enhanced training and methods for use therewith
US11457871B2 (en) 2018-11-21 2022-10-04 Enlitic, Inc. Medical scan artifact detection system and methods for use therewith
US11011257B2 (en) 2018-11-21 2021-05-18 Enlitic, Inc. Multi-label heat map display system
US11282198B2 (en) 2018-11-21 2022-03-22 Enlitic, Inc. Heat map generating system and methods for use therewith
JP7313165B2 (ja) * 2019-03-18 2023-07-24 株式会社Erisa アルツハイマー病生存分析装置及びアルツハイマー病生存分析プログラム
EP3712844B1 (en) * 2019-03-20 2023-08-09 Siemens Healthcare GmbH Method and system for fast assessment of brain change normality
US20220237785A1 (en) * 2019-05-23 2022-07-28 H. Lee Moffitt Cancer Center And Research Institute, Inc. System and method of evaluating neural networks to segment medical images
TWI697686B (zh) 2019-06-20 2020-07-01 臺北榮民總醫院 基於磁振造影分析腦組織成分的系統與方法
US11462315B2 (en) 2019-11-26 2022-10-04 Enlitic, Inc. Medical scan co-registration and methods for use therewith
CN112651924B (zh) * 2020-12-04 2024-03-26 深圳博脑医疗科技有限公司 一种数据生成装置、方法、终端和存储介质
US11669678B2 (en) 2021-02-11 2023-06-06 Enlitic, Inc. System with report analysis and methods for use therewith
JP2022135392A (ja) * 2021-03-05 2022-09-15 コニカミノルタ株式会社 医用情報管理装置、医用情報管理方法および医用情報管理プログラム
DE102021209169A1 (de) 2021-08-20 2023-02-23 Siemens Healthcare Gmbh Validierung KI-basierter Ergebnisdaten

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4022587B2 (ja) * 2004-02-23 2007-12-19 国立精神・神経センター総長 脳疾患の診断支援方法及び装置
JP4025823B2 (ja) * 2004-02-24 2007-12-26 国立精神・神経センター総長 脳疾患の診断支援方法及び装置
US8218835B2 (en) * 2005-08-22 2012-07-10 Dai Nippon Printing Co., Ltd. Method for assisting in diagnosis of cerebral diseases and apparatus thereof
US9101282B2 (en) * 2005-09-22 2015-08-11 Brainlab Ag Brain tissue classification
WO2007114238A1 (ja) * 2006-03-30 2007-10-11 National University Corporation Shizuoka University 脳萎縮判定装置、脳萎縮判定方法及び脳萎縮判定プログラム
CN101479762A (zh) * 2006-06-23 2009-07-08 皇家飞利浦电子股份有限公司 用于在包含图像值的图像中确定阈值的方法、系统和计算机程序
GB2451416B (en) * 2006-09-07 2010-05-19 Siemens Molecular Imaging Ltd ROI-based assessment of abnormality using transformation invariant features
JP5098393B2 (ja) * 2007-03-28 2012-12-12 大日本印刷株式会社 関心領域決定装置
WO2009108135A1 (en) * 2008-02-29 2009-09-03 Agency For Science, Technology And Research A method and system for anatomy structure segmentation and modeling in an image
JP2010046103A (ja) * 2008-08-19 2010-03-04 Osaka Univ 診断画像生成システム、診断画像生成プログラム、及び診断画像生成方法
JP5532730B2 (ja) * 2009-08-05 2014-06-25 大日本印刷株式会社 Ｃｎｒ測定装置，方法及びコンピュータプログラム
JP5229175B2 (ja) * 2009-09-25 2013-07-03 大日本印刷株式会社 医用画像表示処理方法、装置およびプログラム
WO2011040473A1 (ja) * 2009-09-29 2011-04-07 大日本印刷株式会社 医用画像処理方法、装置およびプログラム
CN102858240B (zh) * 2010-04-22 2015-07-29 株式会社日立医疗器械 磁共振成像装置
EP2647335B1 (en) * 2010-12-02 2020-06-03 Dai Nippon Printing Co., Ltd. Medical image processing device

Also Published As

Publication number Publication date
JP2013066632A (ja) 2013-04-18
EP2762072A1 (en) 2014-08-06
CN103826536A (zh) 2014-05-28
US9390509B2 (en) 2016-07-12
WO2013047278A1 (ja) 2013-04-04
EP2762072B1 (en) 2018-09-05
CN103826536B (zh) 2016-11-16
EP2762072A4 (en) 2015-03-11
US20140341471A1 (en) 2014-11-20

Similar Documents

Publication Publication Date Title
JP5970766B2 (ja) 医用画像処理装置、医用画像処理方法、プログラム
US10517504B2 (en) Systems and methods for automated voxelation of regions of interest for magnetic resonance spectroscopy
JP4025823B2 (ja) 脳疾患の診断支援方法及び装置
JP5699936B2 (ja) 医用画像処理装置の作動方法、装置およびプログラム
JP5468905B2 (ja) 神経変性疾患の診断を支援するためのツール
US8050734B2 (en) Method and system for performing patient specific analysis of disease relevant changes of a disease in an anatomical structure
US9361686B2 (en) Method and apparatus for the assessment of medical images
WO2007023522A1 (ja) 脳疾患の診断支援方法及び装置
Kim et al. Automatic hippocampal segmentation in temporal lobe epilepsy: impact of developmental abnormalities
JP5098393B2 (ja) 関心領域決定装置
Dwyer et al. Improved longitudinal gray and white matter atrophy assessment via application of a 4-dimensional hidden Markov random field model
US20200265276A1 (en) Copd classification with machine-trained abnormality detection
JP6273291B2 (ja) 画像処理装置および方法
JP2012045387A (ja) 局所的臨床特徴を分析して視覚化するためのシステムおよび方法
EP2747658B1 (en) Method to compute and present brain amyloid in gray matter
US20110064289A1 (en) Systems and Methods for Multilevel Nodule Attachment Classification in 3D CT Lung Images
CN112700451A (zh) 用于自动分割3d医学图像的方法、系统和计算机可读介质
JP4022587B2 (ja) 脳疾患の診断支援方法及び装置
US9020215B2 (en) Systems and methods for detecting and visualizing correspondence corridors on two-dimensional and volumetric medical images
Wittens et al. Inter-and intra-scanner variability of automated brain volumetry on three magnetic resonance imaging systems in Alzheimer’s disease and controls
JP5524589B2 (ja) 診断支援システム、方法及びコンピュータプログラム
Saurabh A robust and sensitive voxel-based method for measuring cortical thickness change using fuzzy correspondence

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140707

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140707

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150728

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150924

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20160119

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160419

A911 Transfer of reconsideration by examiner before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20160426

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160614

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160627

R150 Certificate of patent or registration of utility model

Ref document number: 5970766

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150