(57)【要約】【課題】異なる疾患の比較に適した診断支援装置等を提供する。特に、複数の疾患が想定される場合において、特定の疾患を精度よく判別する。【解決手段】本発明に係る診断支援装置は、想定される複数の疾患から特定の疾患を判定する診断支援装置であって、脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出手段と、算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定手段と、算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定手段と、を備えることを特徴とする。【選択図】図１
2019136519.tif 000002
【特許請求の範囲】
【請求項１】
想定される複数の疾患から特定の疾患を判定する診断支援装置であって、
脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出手段と、
算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定手段と、
算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定手段と、
を備えることを特徴とする診断支援装置。
【請求項２】
前記算出手段は、前記萎縮の評価値として、前記特定の疾患とは異なる他の疾患の関心領域における前記萎縮度に対する前記特定の疾患の関心領域における前記萎縮度の比である萎縮比を複数の組織で算出する
ことを特徴とする請求項１に記載の診断支援装置。
【請求項３】
前記算出手段は、前記萎縮比を算出する際に、前記他の疾患の関心領域における萎縮度として、当該他の疾患の関心領域において萎縮の傾向が強く現れる組織の萎縮度を用いる
ことを特徴とする請求項２に記載の診断支援装置。
【請求項４】
前記特定の疾患はレビー小体型認知症であり、前記他の疾患はアルツハイマー型認知症であり、
前記設定手段は、各組織の萎縮比がともに０．２以下の閾値より大きいことを、レビー小体型認知症を判定するための条件として設定する
ことを特徴とする請求項２または請求項３に記載の診断支援装置。
【請求項５】
コンピュータが、想定される複数の疾患から特定の疾患を判定する診断支援方法であって、
前記コンピュータが、
脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出ステップと、
算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定ステップと、
算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定ステップと、
を実行することを特徴とする診断支援方法。
【請求項６】
コンピュータを、想定される複数の疾患から特定の疾患を判定する診断支援装置として機能させるプログラムであって、
前記コンピュータを、
脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出手段、
算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定手段、
算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定手段、
として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＭＲＩ等により撮像された脳画像に基づき脳疾患の診断支援を行う診断支援技術に関し、特に、複数の疾患が想定される場合に適した診断支援を行う技術に関する。
【背景技術】
【０００２】
近年、ＳＰＥＣＴ（Ｓｉｎｇｌｅ  Ｐｈｏｔｏｎ  Ｅｍｉｓｓｉｏｎ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）やＰＥＴ（Ｐｏｓｉｔｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇａｒｐｈｙ）等の核医学検査や、ＣＴ（Ｃｏｍｐｕｔｅｒｉｚｅｄ  Ｔｏｍｏｇｒａｐｈｙ）やＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｃｅ  Ｉｍａｇｉｎｇ）によって脳の状態に関する情報が取得可能になってきている。
【０００３】
特に、脳の組織の萎縮に関しては、ＭＲＩ画像によって特定部位の容積を求め、その相対的な大きさを比較して異常の有無を判別できる。例えば、特許文献１では、アルツハイマー型認知症の診断支援を行うシステムが開示されており、ＭＲＩ画像を用いて内側側頭部の萎縮を定量的に評価することで、アルツハイマー型認知症の診断支援を行うことを可能としている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許４０２５８２３号
【特許文献２】特開２０１３−６６６３２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来の診断支援システム等においては、特定の疾患を対象とした有効な診断支援情報は得られるものの、同時に異なる疾患が想定される場合に、これら疾患を比較する有効な診断支援情報を提供するまでには至っていない。
【０００６】
本発明は、上述した課題に鑑みてなされたものであり、その目的は、異なる疾患の比較に適した診断支援装置等を提供することにある。
【課題を解決するための手段】
【０００７】
前述した目的を達成するための第１の発明は、想定される複数の疾患から特定の疾患を判定する診断支援装置であって、脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出手段と、算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定手段と、算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定手段と、を備えることを特徴とする診断支援装置である。
【０００８】
前述した目的を達成するための第２の発明は、コンピュータが、想定される複数の疾患から特定の疾患を判定する診断支援方法であって、前記コンピュータが、脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出ステップと、算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定ステップと、算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定ステップと、を実行することを特徴とする診断支援方法である。
【０００９】
前述した目的を達成するための第３の発明は、コンピュータを、想定される複数の疾患から特定の疾患を判定する診断支援装置として機能させるプログラムであって、前記コンピュータを、脳画像から、特定の疾患の関心領域における脳の萎縮の評価値を、複数の組織で算出する算出手段、算出した各組織の萎縮の評価値のそれぞれに前記特定の疾患を判定するための条件を設定する設定手段、算出した各組織の萎縮の評価値のそれぞれが対応する前記条件を満たす場合に、前記特定の疾患であると判定する判定手段、として機能させることを特徴とするプログラムである。
【発明の効果】
【００１０】
本発明により、異なる疾患の比較に適した診断支援装置等を提供することができる。特に、複数の疾患が想定される場合において、特定の疾患を精度よく判別することができる。
【図面の簡単な説明】
【００１１】
【図１】本実施形態に係る診断支援装置の機能を示すブロック図
【図２】本実施形態に係る診断支援装置の処理の手順を示すフローチャート
【図３】萎縮スコアの算出処理の手順を示すフローチャート
【図４】診断支援情報等の表示の一例を示す図
【図５】スライス画像を拡大した図
【図６】組織毎の萎縮比を組み合わせて疾患の判別を行う例を示す図
【図７】ＡＤとＤＬＢの萎縮比のプロット図
【発明を実施するための形態】
【００１２】
以下図面に基づいて、本発明の実施形態を詳細に説明する。
【００１３】
図１は、診断支援装置１の機能を示すブロック図である。診断支援装置１は、ユーザインターフェース部１０、処理部２０、データベース部３０を有している。
【００１４】
ユーザインターフェース部１０は、主に、ＭＲＩ画像の入力を受付ける画像入力機能１１と、処理部２０で処理された結果を表示する表示機能１３とを有する。
【００１５】
処理部２０は、主に、ユーザインターフェース部１０から入力されたＭＲＩ画像を処理する画像処理機能２１と、Ｚスコア等の各種統計指標を算出する統計処理機能２３と、比較する各疾患に特異的な部位（関心領域）を特定する部位特定機能２５、萎縮度を算出する萎縮度算出機能２７、萎縮比を算出する萎縮比算出機能２９等を有する。
【００１６】
また、データベース部３０には、処理部２０の処理に供する灰白質脳画像テンプレート３１、白質脳画像テンプレート３３、健常者画像データベース３５、関心領域ＲＯＩ３７等が保存されている。
【００１７】
上記の灰白質脳画像テンプレート３１、白質脳画像テンプレート３３は、灰白質と白質それぞれについて作成されたものが予め前記データベース部３０に保存されているものとする。各テンプレートは年齢や性別などの被験者の属性に応じて層別に作成されていてもよい。
尚、本実施形態においては、上記テンプレートを作成する際の解剖学的標準化の手法としてＤＡＲＴＥＬ（Ｄｉｆｆｅｏｍｏｒｐｈｉｃ  Ａｎａｔｏｍｉｃａｌ  Ｒｅｇｉｓｔｒａｔｉｏｎ  Ｔｈｒｏｕｇｈ  Ｅｘｐｏｎｅｎｔｉａｔｅｄ  Ｌｉｅ  ａｌｇｅｂｒａ）を採用する。ＤＡＲＴＥＬを用いたテンプレート作成の処理は、特許文献１と同様であるため、記載を省略する。
【００１８】
［診断支援装置１の処理］
図２は、本実施形態に係る診断支援装置１の処理を示すフローチャートである。尚、この処理はコンピュータからなる処理部２０においてプログラムにより実行可能なものである。
【００１９】
ステップＳ１において、診断支援装置１（画像入力機能１１）は、被験者のＭＲＩ脳画像の入力を受付ける。
【００２０】
ステップＳ２において、診断支援装置１は、ステップＳ１において入力された被験者のＭＲＩ脳画像に基づき脳の萎縮の程度を表す「萎縮スコア」を算出する。
【００２１】
＜萎縮スコア算出処理＞
ここで、図３のフローチャートを参照しながら、上記ステップＳ２における萎縮スコアの算出処理を説明する。
【００２２】
（画像再構成）
診断支援装置１は、入力された被験者のＭＲＩ脳画像に対し「画像再構成」を行う（ステップＳ２１）。
【００２３】
画像再構成は、最初に、入力された被験者のＭＲＩ脳画像を、脳全体を含むように所定の厚さのスライス状に撮像した、例えば１００〜２００枚のＴ１強調ＭＲＩ画像に変換する。この際、各スライス画像におけるボクセルの各辺の長さが予め等しくなるようにスライス画像のリサンプリングを行う。
【００２４】
そして、上記処理を施した被験者のＭＲＩ脳画像に対し、標準的な脳画像と空間的な位置合わせを行う。具体的には、被験者のＭＲＩ脳画像に対して、線形変換（アフィン変換）、トリミング等を施し、標準的な脳画像と位置、角度、サイズ等を合わせる。これにより、ＭＲＩ撮影時の被験者の頭の位置のずれ等が画像上で補正され、標準的な脳画像と比較する際の精度が向上する。
【００２５】
（組織分割）
ステップＳ２１の画像再構成がなされた後、診断支援装置１は「組織分割」を行い、灰白質と白質を抽出した灰白質脳画像と白質脳画像を作成する（ステップＳ２２）。
【００２６】
前述のＴ１強調ＭＲＩ脳画像には、神経線維に対応する高い信号値を呈する白質、神経細胞に対応する中間の信号値を呈する灰白質、低い信号値を呈する脳脊髄液の３種類の組織が含まれているため、この信号値の差に着目して灰白質と白質とをそれぞれ抽出する処理を行う。この処理は、特許文献１や特許文献１より抽出精度を高めた特許文献２に記載されている処理と同様であるため、記載を省略する。
【００２７】
（解剖学的標準化）
そして、診断支援装置１は、ステップＳ２２において作成された灰白質脳画像及び白質脳画像に対し「解剖学的標準化」を行う（ステップＳ２３）。
【００２８】
解剖学的標準化とは、標準脳画像へのボクセルの位置合わせを行うものである。本実施形態では、ＤＡＲＴＥＬによる解剖学的標準化を実行する。ＤＡＲＴＥＬの処理については特許文献１と同様であるため、記載を省略する。
【００２９】
ＤＡＲＴＥＬによる解剖学的標準化を施した灰白質脳画像と白質脳画像に対して、Ｓ／Ｎ比の向上等を目的に、画像平滑化の処理を行う。このように画像平滑化を行うことにより、解剖学的標準化処理で完全に一致しない個体差を低減させることができる。こちらも具体的な処理については、特許文献１と同様である。
【００３０】
また、その後、比較対象となる健常者の画像群におけるボクセル値の分布と合わせるために、脳全体のボクセル値を補正する濃度値補正を行う。
【００３１】
（比較）
ステップＳ２４において、診断支援装置１は、健常者のＭＲＩ画像との比較を行い、被験者の脳の萎縮の程度を示す「萎縮スコア」を算出する。本実施形態では、萎縮スコアとして統計的指標である「Ｚスコア」を用いる。
【００３２】
具体的には、上記したステップＳ２３により解剖学的標準化、画像平滑化等を施した被験者の灰白質脳画像及び白質脳画像と、予め収集して前記データベース部３０の健常者画像データベース３５に保存してある健常者の灰白質及び白質のＭＲＩ脳画像群との統計的比較を行い、ＭＲＩ脳画像の全ボクセル又は特定領域のボクセルについて灰白質及び白質のＺスコアを次のように算出する。以降、灰白質のＺスコアをＺ[灰白質]、白質のＺスコアをＺ[白質]と表す。
【００３３】
【数１】
2019136519.tif 000003
【００３４】
【数２】
2019136519.tif 000004
【００３５】
上式に示すように、Ｚスコアは、被験者画像のボクセル値と、健常者画像群の対応するボクセル値の平均との差を、標準偏差によってスケーリングした値であり、灰白質及び白質の容積の相対的低下の度合を示すものである。
【００３６】
尚、Ｚスコアに限らず、被験者画像と健常者画像とのボクセル値の大小が判断できるその他の指標を萎縮の程度を示す萎縮スコアとして用いてもよい（例えばｔスコア等）。
【００３７】
また、ステップＳ２４において使用した健常者画像データベース３５に保存される健常者の灰白質及び白質のＭＲＩ脳画像は、予め収集した健常者の画像群のそれぞれに対しステップＳ２１〜ステップＳ２３の「画像再構成」→「組織分割」→「解剖学的標準化」及び画像平滑化等の同様の処理を順次適用して作成されたものである。尚、この処理はコンピュータからなる処理部２０においてプログラムにより実行可能なものである。
【００３８】
以上説明した処理により、被験者のＭＲＩ脳画像から萎縮スコア（本実施形態ではＺスコア）が算出される。
【００３９】
図２のフローチャートに戻る。ステップＳ３において、診断支援装置１は、比較する各疾患に特異的な脳の部位（関心領域）を特定する。これは主に、処理部２０の部位特定機能２５により実現される。
【００４０】
例えば、診断支援装置１は、統計的処理に基づいて各疾患に関連する関心領域を特定する。具体的には、ある疾患に対応する関心領域を特定する場合、その疾患を持つ患者のＭＲＩ画像群（疾患画像群）と、それ以外の人の画像群（非疾患者画像群）とについて、ボクセル単位で２群間の有意差を統計的に検定する２標本t検定を行い、有意差が認められたボクセルを、その疾患における特徴的なボクセルとみなし、その座標の集合をその疾患に対応する関心領域（Ｒｅｇｉｏｎ  Ｏｆ  Ｉｎｔｅｒｅｓｔ：ＲＯＩ）として特定する。
また、特許５０９８３９３号記載のように有意水準と経験則の両方を考慮してＲＯＩを特定してもよい。
また、疾患画像（群）のみからＲＯＩを特定してもよい。例えば、疾患画像（群）について、脳全体における萎縮の大きさに相関して萎縮が大きい部位をＲＯＩとして特定してもよい。
その他、診断者等の主観により手動でＲＯＩを特定してもよい。
【００４１】
以降、本実施形態においては異なる疾患Ａと疾患Ｂを想定し、疾患Ａの関心領域ＲＯＩＡ、疾患Ｂの関心領域ＲＯＩＢが、ステップＳ３において特定されたものとして説明を行う。
【００４２】
（比較表示）
ステップＳ４において、診断支援装置１は、ステップＳ３において特定した各部位の診断支援情報等を比較表示する。
【００４３】
ここで、ステップＳ４において表示される「萎縮度」及び「萎縮比」について説明する。これらの指標は、主に、処理部２０の主に萎縮度算出機能２５、萎縮比算出機能２７により算出される。
【００４４】
＜萎縮度＞
診断支援装置１は、「関心領域内」における萎縮の程度を示す「萎縮度」を算出する。また、萎縮度を灰白質及び白質の組織毎に算出することで、各疾患に関連する部位の萎縮の程度を組織毎に定量的に評価することができるようにする。
【００４５】
具体的には、関心領域ＲＯＩＡ内の「灰白質」の萎縮度（数式３）、及び「白質」の萎縮度（数式４）は次のようにＺスコアから算出される。
【００４６】
【数３】
2019136519.tif 000005
【００４７】
【数４】
2019136519.tif 000006
【００４８】
また、関心領域ＲＯＩＢ内の「灰白質」の萎縮度（数式５）、及び「白質」の萎縮度（数式６）は次のように算出される。
【００４９】
【数５】
2019136519.tif 000007
【００５０】
【数６】
2019136519.tif 000008
【００５１】
尚、本実施形態では、萎縮度として関心領域内の「正のＺスコアの平均値」を採用しているが、これに限らず、任意に閾値を定めた「閾値を超えるＺスコアの平均値」または単に「Ｚスコアの平均値」を採用してもよい。また、ＲＯＩ内の総ボクセル数に対してＺスコアが閾値を超えるボクセルが占める割合を採用してもよい。
【００５２】
＜萎縮比＞
診断支援装置１は、上記した萎縮度に基づき更に「萎縮比」を算出する。ここで、「萎縮比」とは異なる疾患が想定されるときにある疾患を基準として他の疾患の特徴がどれだけ大きいかを示す指標のことを言う。前述した萎縮度は、疾患毎の関心領域における萎縮の程度を個別的に把握できるものの、各疾患の関連性を一意に把握できる指標ではないため、各疾患を鑑別支援する指標としては十分ではない。そこで、本実施形態では、更に、上記算出した各疾患の萎縮度の比である「萎縮比」を定義し、これを各疾患の鑑別支援の指標として用いる。
例えば、疾患Ａ、疾患Ｂが想定される場合、疾患Ａを基準とした疾患Ｂの萎縮比は次のように算出される。
【００５３】
【数７】
2019136519.tif 000009
【００５４】
【数８】
2019136519.tif 000010
【００５５】
萎縮比も萎縮度と同様に、灰白質及び白質について算出する。尚、基準となる疾患Ａの萎縮度（数式の分母）については、疾患Ａの患者においてＲＯＩＡ内で萎縮の傾向が強く現れる組織（「灰白質」又は「白質」のいずれか）を選択することが望ましい。
上式により、萎縮比の値が小さいと疾患Ａの傾向が強く、萎縮比の値が大きいと疾患Ｂの傾向が強いものと判断でき、各疾患の鑑別を支援する指標となりうる。例えば、適切な閾値を設定し、数式（７）（８）の萎縮比が閾値より小さい場合には疾患Ａ、萎縮比が閾値より大きい場合には疾患Ｂ、といった具合に各疾患を判別可能である。
【００５６】
尚、本実施形態においては、灰白質と白質の萎縮度及び萎縮比を個別的に算出しているが、灰白質と白質を合成して１つの萎縮度及び萎縮比を算出してもよい。これにより、両方の組織の萎縮を１つの指標を用いて評価できる。これは、例えば、ある疾患において灰白質及び白質の両方の組織が萎縮するような場合、若しくは、灰白質及び白質のいずれの組織が萎縮するかを判断できないような場合に有効な指標となりうる。
但し、この場合、健常者データベース３５には予め灰白質脳画像と白質脳画像を合成した健常者画像群が予め用意されており、これら画像群と被験者の灰白質脳画像及び白質脳画像を合成した被験者画像との比較により、Ｚスコアを算出し、当該Ｚスコアから萎縮度及び萎縮比を算出する。
【００５７】
例えば、Ｚスコアは下記のように算出される。
【００５８】
【数９】
2019136519.tif 000011
【００５９】
また、疾患Ａの関心領域ＲＯＩＡ、疾患Ｂの関心領域ＲＯＩＢにおける萎縮度は次のように算出される。
【００６０】
【数１０】
2019136519.tif 000012
【００６１】
【数１１】
2019136519.tif 000013
【００６２】
さらに、上記した萎縮度から例えば疾患Ａを基準とした疾患Ｂの萎縮比は次のように算出される。
【００６３】
【数１２】
2019136519.tif 000014
【００６４】
尚、基準となる疾患Ａの萎縮度（数式の分母）については、疾患Ａの患者においてＲＯＩＡ内で萎縮の傾向が強く現れる組織の形態（「灰白質及び白質」、「灰白質」、又は「白質」のいずれか）を選択することが望ましい。
【００６５】
以上説明した種々の指標（Ｚスコア、萎縮度、萎縮比）は、ユーザインターフェース部１０（表示機能１３）に表示される。
図４は、診断支援装置１のユーザインターフェース部１０の表示例を示す。
【００６６】
図４の表示エリア４１には、脳のスライス画像が所定間隔おきに並べて表示される。そして、スライス画像上に「灰白質」のＺスコア（数式１）の分布（Ｚスコアマップ）が重ねて表示され、また、スライス画像上に疾患Ａの関心領域と疾患Ｂの関心領域が表示される。
【００６７】
図５は、一枚のスライス画像を拡大した図である。スライス画像上に、Ｚスコアマップ５ａが表示され、また、スライス画像面での疾患Ａの関心領域５ｂが実線、疾患Ｂの関心領域５ｃが破線で表示される。これにより、スライス画像全体での萎縮の程度を把握できるとともに、スライス画像中の着目部位（関心領域５ｂ、５ｃ）の萎縮の程度を把握することができる。
【００６８】
図４の表示エリア４２には、表示エリア４１と同様に脳のスライス画像が所定間隔おきに並べて表示される。但し、表示エリア４２においては、スライス画像上に「白質」のＺスコア（数式２）の分布が重ねて表示される。
【００６９】
このように表示エリア４１及び表示エリア４２において、組織毎（灰白質、白質）にＺスコアの分布を表示することで、組織毎の萎縮の違いを把握することができる。
【００７０】
尚、Ｚスコアの表示の手段は様々であり、例えば、Ｚスコアの値に応じて濃淡を変えて表示してもよいし、等高線を用いて表示してもよい。また、各疾患における関心領域の表示の手段も様々であり、例えば、疾患ごとに関心領域を色分けして表示するようにしてもよい。
【００７１】
図４の表示エリア４３には、疾患Ａの関心領域（図４では部位Ａと示す）の「灰白質」の萎縮度（数式３）、及び「白質」の萎縮度（数式４）が数値表示され、また、疾患Ｂの関心領域（図４では部位Ｂと示す）の「灰白質」の萎縮度（数式５）、及び「白質」の萎縮度（数式６）が数値表示される。
【００７２】
そして、図４の表示エリア４４には、疾患Ａを基準とした疾患Ｂの「灰白質」の萎縮比（数式７）、及び「白質」の萎縮比（数式８）が数値表示される。
【００７３】
以上、本実施形態においては、診断支援装置１によって、各疾患の関心領域が特定され、特定した各疾患の関心領域に関する種々の指標を用いて各疾患を比較表示する。これにより、異なる疾患を比較又は鑑別支援する有効な診断支援情報が診断者に提供される。
【００７４】
[実施例]
ここでは、一実施例として、アルツハイマー型認知症（以下、「ＡＤ」と呼ぶ）とレビー小体型認知症（以下、「ＤＬＢ」と呼ぶ）の２つの疾患を対象として鑑別支援の可能性を検証した。
【００７５】
このうち、ＡＤにおいては、内側側頭部の灰白質に強い萎縮が観測されるため、ＭＲＩ画像を用いて内側側頭部の萎縮を定量的に評価することで、ＡＤの診断支援が可能であることは既に知られている。
【００７６】
一方、ＤＬＢに関しては、ＭＲＩにおける疾患特異性について未だエビデンスが少ないのが現状である。しかしながら、近年の研究により、ＤＬＢは中脳（背側）の灰白質が萎縮するとの報告がなされている（  Whitwell, Jennifer L．et al．”Focal atrophy in dementia with Lewy bodies on MRI: a distinct
pattern from Alzeimer’s disease．” Brain(2007)  ）。また、別の研究によれば、ＤＬＢは中脳（背側）・橋（背側）・小脳の白質が萎縮するとの報告がなされている（  Nakatsuka, et al．
“Discrimination of dementia with Lewy
bodies from Alzheimer’s disease using voxel-based
morphometry of white matter by statistical parametric mapping 8 plus
diffeomorphic anatomic registration through exponentiated Lie algebra．
”Neuroradiology(2013)  ）。これら先行研究の知見によれば、ＤＬＢにおいては後部脳幹付近に特異な傾向があるものと推測される。
【００７７】
実際に、診断支援装置１の部位特定機能２５によりＡＤ及びＤＬＢに特異的な部位（萎縮が大きい部位）を特定した結果、ＡＤの場合は「内側側頭部」付近が特異的な部位として現れ、ＤＬＢの場合は「後部脳幹」付近が特異的な部位として現れた。そこで、本実施例では、これら部位をＡＤとＤＬＢの関心領域に設定した。ここで、ＡＤの関心領域（内側側頭部付近）をＲＯＩＡ、ＤＬＢの関心領域（後部脳幹付近）をＲＯＩＢとする。
【００７８】
次に、診断支援装置１の萎縮度算出機能２７により算出する萎縮度を検討する。ＡＤにおいては、前述したように内側側頭部付近の灰白質の萎縮が大きいため、数式３に基づく「灰白質」の萎縮度を採用した。ＤＬＢにおいては、前述した先行研究によれば、後部脳幹付近の灰白質、白質のいずれか又はその両方に影響が出る可能性があるため、組織毎に評価することとした。そこで、ＤＬＢの萎縮度として、数式５に基づく「灰白質」の萎縮度、及び数式６に基づく「白質」の萎縮度の両方を利用することとした。
【００７９】
そして、診断支援装置１の萎縮比算出機能２９により算出する萎縮比として、数式７に基づく灰白質の萎縮比（以降、「τ１」と示す）、及び数式８に基づく白質の萎縮比（以降、「τ２」と示す）を採用した。ここで、数式７、８の萎縮比の分母に相当するＲＯＩＡ（内側側頭部付近）内のＡＤの萎縮度は、ＡＤ患者に萎縮の傾向が大きく表れる「灰白質」の萎縮度を用いた。
【００８０】
図６は、上記した灰白質の萎縮比τ１及び白質の萎縮比τ２によりＡＤとＤＬＢを判別する一例を示す。
図６に示すように、萎縮比τ１に対して閾値α１、萎縮比τ２に対して閾値α２を設定し、τ１＞α１及びτ２＞α２の両方を満たす場合、ＤＬＢの疑いがあると判別できる。それ以外の場合は、ＡＤの疑いがあると判別できる。
【００８１】
図７は、実際にＡＤの患者とＤＬＢの患者を対象として上記した萎縮比τ１及び萎縮比τ２を診断支援装置１により算出した値を、模式的にプロットした図を
示す。灰色のドットがＡＤと診断された患者を示し、黒色のドットがＤＬＢと診断された患者を示す。本実施形態では閾値α１、α２はともに０．２に設定した。
図７に示すようにτ１＞α１及びτ２＞α２を満たすエリアにＤＬＢ患者が多く分布し、それ以外のエリアにＡＤ患者が多く分布しているのが分かり、良好な判別結果が得られた。これにより「萎縮比」がＡＤとＤＬＢを対象とした鑑別支援の指標として有効であることが確認された。
【００８２】
以上、添付図面を参照しながら、本発明の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【００８３】
１・・・・診断支援装置
１０・・・ユーザインターフェース部
１１・・・画像入力機能
１３・・・表示機能
２０・・・処理部
２１・・・画像処理機能
２３・・・統計処理機能
２５・・・部位特定機能
２７・・・萎縮度算出機能
２９・・・萎縮比算出機能
３０・・・データベース部
３１・・・灰白質脳画像テンプレート
３３・・・白質脳画像テンプレート
３５・・・健常者データベース
３７・・・関心領域ＲＯＩ
【図１】
2019136519.tif 000015
【図２】
2019136519.tif 000016
【図３】
2019136519.tif 000017
【図４】
2019136519.tif 000018
【図５】
2019136519.tif 000019
【図６】
2019136519.tif 000020
【図７】
2019136519.tif 000021