(57)【要約】【課題】判別器の学習装置、方法およびプログラム、梗塞領域を判別する判別器、並びに梗塞領域判別装置、方法およびプログラムにおいて、ＣＴ画像を用いて脳梗塞が発症してからの経過時間に応じて梗塞領域を迅速に判別できるようにする。【解決手段】画像取得部２１が、脳梗塞を発症している被検体の脳のＣＴ画像および１以上のＭＲＩ画像を取得し、梗塞領域抽出部２２が、ＭＲＩ画像から発症からの経過時間に応じた梗塞領域を抽出する。位置合わせ部２３が、ＣＴ画像とＭＲＩ画像との位置合わせを行い、梗塞領域特定部２５が、位置合わせ結果に基づいて、ＣＴ画像Ｂｃ０における発症からの経過時間に応じた梗塞領域を特定する。学習部２５が、ＣＴ画像Ｂｃ０において特定された発症からの経過時間に応じた梗塞領域を教師データとして、判別対象のＣＴ画像における、発症からの経過時間に応じた梗塞領域を判別する判別器２６の学習を行う。【選択図】図２
2020032045.tif 000002
【特許請求の範囲】
【請求項１】
脳梗塞を発症している被検体の脳のＣＴ画像および該ＣＴ画像を取得するＣＴ装置とは異なるモダリティにより取得された、前記被検体の少なくとも１つの脳画像を取得する画像取得部と、
少なくとも前記脳画像に基づいて、前記脳画像から、発症からの経過時間に応じた梗塞領域を抽出する梗塞領域抽出部と、
前記ＣＴ画像と前記脳画像との位置合わせを行う位置合わせ部と、
前記位置合わせの結果に基づいて、前記ＣＴ画像における発症からの経過時間に応じた前記梗塞領域を特定する梗塞領域特定部と、
前記ＣＴ画像において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う学習部とを備えた判別器の学習装置。
【請求項２】
前記脳画像はＭＲＩ画像である請求項１に記載の学習装置。
【請求項３】
前記ＭＲＩ画像は拡散強調画像、ＦＬＡＩＲ画像およびＴ２強調画像の少なくとも１つである請求項２に記載の判別器の学習装置。
【請求項４】
請求項１から３のいずれか１項に記載の判別器の学習装置により学習がなされた判別器。
【請求項５】
判別対象のＣＴ画像を取得する画像取得部と、
前記判別対象のＣＴ画像における梗塞領域を発症からの経過時間に応じて判別する、請求項４に記載の判別器とを備えた梗塞領域判別装置。
【請求項６】
前記判別器による前記梗塞領域の判別結果を表示部に表示する表示制御部をさらに備えた請求項５に記載の梗塞領域判別装置。
【請求項７】
脳梗塞を発症している被検体の脳のＣＴ画像および該ＣＴ画像を取得するＣＴ装置とは異なるモダリティにより取得された、前記被検体の少なくとも１つの脳画像を取得し、
少なくとも前記脳画像に基づいて、前記脳画像から、発症からの経過時間に応じた梗塞領域を抽出し、
前記ＣＴ画像と前記脳画像との位置合わせを行い、
前記位置合わせの結果に基づいて、前記ＣＴ画像における発症からの経過時間に応じた前記梗塞領域を特定し、
前記ＣＴ画像において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う判別器の学習方法。
【請求項８】
判別対象のＣＴ画像を取得し、
請求項４に記載の判別器により、前記判別対象のＣＴ画像における梗塞領域を発症からの経過時間に応じて判別する梗塞領域判別方法。
【請求項９】
脳梗塞を発症している被検体の脳のＣＴ画像および該ＣＴ画像を取得するＣＴ装置とは異なるモダリティにより取得された、前記被検体の少なくとも１つの脳画像を取得する手順と、
少なくとも前記脳画像に基づいて、前記脳画像から、発症からの経過時間に応じた梗塞領域を抽出する手順と、
前記ＣＴ画像と前記脳画像との位置合わせを行う手順と、
前記位置合わせの結果に基づいて、前記ＣＴ画像における発症からの経過時間に応じた前記梗塞領域を特定する手順と、
前記ＣＴ画像において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う手順とをコンピュータに実行させる判別器の学習プログラム。
【請求項１０】
判別対象のＣＴ画像を取得する手順と、
請求項４に記載の判別器により、前記判別対象のＣＴ画像における梗塞領域を発症からの経過時間に応じて判別する手順とをコンピュータに実行させる梗塞領域判別プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、脳画像における梗塞領域を発症からの経過時間に応じて判別する判別器の学習装置、方法およびプログラム、梗塞領域を判別する判別器、並びに梗塞領域判別装置、方法およびプログラムに関する。
【背景技術】
【０００２】
近年、ＣＴ（Computed Tomography）装置およびＭＲＩ（Magnetic Resonance Imaging）装置等の医療機器の進歩により、より質の高い高解像度の医用画像を用いての画像診断が可能となってきている。とくに、対象部位を脳とした場合においては、ＣＴ画像およびＭＲＩ画像等を用いた画像診断により、脳梗塞および脳出血等の脳の血管障害を起こしている領域を特定することができる。このため、画像診断を支援するための各種手法が提案されている。
【０００３】
例えば、特許文献１においては、超音波画像、ＣＴ画像およびＭＲＩ画像等の被検体の生体組織を撮像した複数の医用画像を、予め定められた所定の種類毎に分類するパターンを求めて記憶する学習演算を行い、入力画像を学習データに基づいて認識して表示する手法が提案されている。また、特許文献２においては、ＭＲＩの拡散強調画像（ＤＷＩ：Diffusion Weighted Image）に含まれる脳梗塞部位を検出し、拡散強調画像の異常部位と健常者の拡散強調画像とから、両者の解剖学的位置合わせに必要な位置変換データを取得し、患者の脳の各組織位置が健常者の脳の各組織位置に合致するように、位置変換データを用いて、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置で撮影されたＳＰＥＣＴ画像を変換し、ＳＰＥＣＴ画像上で脳梗塞部位を判別する手法が提案されている。また、特許文献３においては、ＣＴ画像とＭＲＩ画像とを位置合わせして表示して、診断を行う手法が提案されている。
【０００４】
ところで、脳梗塞患者に対しては、抗血小板薬を使った血栓溶解療法が行われる。しかしながら、血栓溶解療法は、脳梗塞の未発症が確認された時刻から４．５時間以内が適用対象であり、時間の経過により梗塞範囲が広がるほど、治療後に出血が生じるリスクが高くなることが知られている。このため、血栓溶解療法の適否を判断するためには、医用画像を用いて梗塞範囲を迅速かつ適切に判別することが必要である。ここで、拡散強調画像においては、梗塞領域は他の領域とは異なる信号値を有するものとなる。とくに急性期の脳梗塞の領域は、拡散強調画像においては他の領域との信号値の相違が顕著となる。このため、梗塞領域の確認には拡散強調画像が使用されることが多い。さらに、血栓溶解療法が適用対象となる時間を経過してしまった場合には、脳梗塞に対して血栓溶解療法とは異なる治療が必要である。したがって、ＭＲＩの拡散強調画像のみならず、Ｔ１強調画像、Ｔ２強調画像またはＦＬＡＩＲ（FLuid-Attenuated Inversion Recovery）画像等を用いて、超急性期、急性期、亜急性期および慢性期（陳旧性）といった、脳梗塞が発症してからの経過時間を認識することも脳梗塞を治療する上で重要である。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】再公表２００９−１４８０４１号公報
【特許文献２】特開２０１３−１６５７６５号公報
【特許文献３】特開２０１４−５１８５１６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
一方、脳の診断においては、脳梗塞を確認する前に、脳内の出血の有無を確認することが多い。脳内の出血はＣＴ画像において明確に確認することができるため、脳疾患が疑われる患者に対しては、まずＣＴ画像を用いた診断が行われる。しかしながら、ＣＴ画像においては、急性期の脳梗塞の領域と他の領域との信号値の相違はそれほど大きくないため、ＣＴ画像を用いて急性期の梗塞を特定することは、困難であることが多い。このため、ＣＴ画像を用いた診断後にＭＲＩ画像を取得して脳梗塞が発症しているか否かの診断が行われることとなる。
【０００７】
しかしながら、ＣＴ画像を用いた診断後にＭＲＩ画像を取得して脳梗塞が発症しているか否かの診断を行っていたのでは、脳梗塞が発症してからの経過時間が長くなるため、血栓溶解療法による治療後の出血のリスクが高くなる可能性がある。さらに、上述したような脳梗塞が発症してからの経過時間を認識するためには複数のＭＲＩ画像の取得が必要となるため、患者の負担が大きくなる。
【０００８】
本開示は上記事情に鑑みなされたものであり、ＣＴ画像を用いて脳梗塞が発症してからの経過時間に応じて梗塞領域を迅速に判別できるようにすることを目的とする。
【課題を解決するための手段】
【０００９】
本開示による判別器の学習装置は、脳梗塞を発症している被検体の脳のＣＴ画像およびＣＴ画像を取得するＣＴ装置とは異なるモダリティにより取得された、被検体の少なくとも１つの脳画像を取得する画像取得部と、
少なくとも脳画像に基づいて、脳画像から、発症からの経過時間に応じた梗塞領域を抽出する梗塞領域抽出部と、
ＣＴ画像と脳画像との位置合わせを行う位置合わせ部と、
位置合わせの結果に基づいて、ＣＴ画像における発症からの経過時間に応じた梗塞領域を特定する梗塞領域特定部と、
ＣＴ画像において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う学習部とを備える。
【００１０】
「少なくとも脳画像に基づいて」とは、脳画像以外の他の画像、例えばＣＴ画像にも基づいてもよいことを意味する。
【００１１】
なお、本開示による判別器の学習装置においては、脳画像はＭＲＩ画像であってもよい。
【００１２】
また、本開示による判別器の学習装置においては、ＭＲＩ画像は拡散強調画像、ＦＬＡＩＲ画像およびＴ２強調画像の少なくとも１つであってもよい。
【００１３】
本開示による判別器は、本開示による判別器の学習装置により学習がなされたものである。
【００１４】
本開示による梗塞領域判別装置は、判別対象のＣＴ画像を取得する画像取得部と、
判別対象のＣＴ画像における梗塞領域を発症からの経過時間に応じて判別する、本開示による判別器とを備える。
【００１５】
なお、本開示による梗塞領域判別装置においては、判別器による梗塞領域の判別結果を表示部に表示する表示制御部をさらに備えるものであってもよい。
【００１６】
本開示による判別器の学習方法は、脳梗塞を発症している被検体の脳のＣＴ画像およびＣＴ画像を取得するＣＴ装置とは異なるモダリティにより取得された、被検体の少なくとも１つの脳画像を取得し、
少なくとも脳画像に基づいて、脳画像から、発症からの経過時間に応じた梗塞領域を抽出し、
ＣＴ画像と脳画像との位置合わせを行い、
位置合わせの結果に基づいて、ＣＴ画像における発症からの経過時間に応じた梗塞領域を特定し、
ＣＴ画像において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う。
【００１７】
本開示による梗塞領域判別方法は、判別対象のＣＴ画像を取得し、
本開示による判別器により、判別対象のＣＴ画像における梗塞領域を発症からの経過時間に応じて判別する。
【００１８】
なお、本開示による判別器の学習方法および本開示による梗塞領域判別方法をコンピュータに実行させるためのプログラムとして提供してもよい。
【００１９】
本開示による他の判別器の学習装置は、コンピュータに実行させるための命令を記憶するメモリと、
記憶された命令を実行するよう構成されたプロセッサとを備え、プロセッサは、
脳梗塞を発症している被検体の脳のＣＴ画像およびＣＴ画像を取得するＣＴ装置とは異なるモダリティにより取得された、被検体の少なくとも１つの脳画像を取得し、
少なくとも脳画像に基づいて、脳画像から、発症からの経過時間に応じた梗塞領域を抽出し、
ＣＴ画像と脳画像との位置合わせを行い、
位置合わせの結果に基づいて、ＣＴ画像における発症からの経過時間に応じた梗塞領域を特定し、
ＣＴ画像において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う処理を実行する。
【００２０】
本開示による他の梗塞領域判別装置は、コンピュータに実行させるための命令を記憶するメモリと、
記憶された命令を実行するよう構成されたプロセッサとを備え、プロセッサは、
判別対象のＣＴ画像を取得し、
本開示による判別器により、判別対象のＣＴ画像における梗塞領域を発症からの経過時間に応じて判別する処理を実行する。
【発明の効果】
【００２１】
本開示によれば、ＣＴ画像のみを用いて、発症からの経過時間に応じた脳梗塞の診断を迅速に行うことができる。
【図面の簡単な説明】
【００２２】
【図１】本開示の実施形態による学習装置、判別器および梗塞領域判別装置を適用した、診断支援システムの概要を示すハードウェア構成図
【図２】本実施形態による梗塞領域判別装置の概略構成を示す図
【図３】拡散強調画像およびＦＬＡＩＲ画像の例を示す図
【図４】ＣＴ画像、拡散強調画像およびＦＬＡＩＲ画像の信号値と梗塞領域の発症からの経過時間との関係を示す図
【図５】ＣＴ画像と拡散強調画像との位置合わせを説明するための図
【図６】ＣＴ画像における発症からの経過時間に応じた梗塞領域の特定を説明するための図
【図７】判別結果の表示の例を示す図
【図８】本実施形態において学習時に行われる処理を示すフローチャート
【図９】本実施形態において発症からの経過時間に応じた梗塞領域の判別時に行われる処理を示すフローチャート
【発明を実施するための形態】
【００２３】
以下、図面を参照して本開示の実施形態について説明する。図１は、本開示の実施形態による判別器の学習装置、判別器および梗塞領域判別装置を適用した、診断支援システムの概要を示すハードウェア構成図である。図１に示すように、診断支援システムでは、本実施形態による梗塞領域判別装置１、３次元画像撮影装置２、および画像保管サーバ３が、ネットワーク４を経由して通信可能な状態で接続されている。なお、梗塞領域判別装置１には、本実施形態による学習装置および判別器が内包される。
【００２４】
３次元画像撮影装置２は、被検体の診断対象となる部位を撮影することにより、その部位を表す３次元画像を医用画像として生成する装置である。３次元画像撮影装置２により生成された医用画像は画像保管サーバ３に送信され、保存される。なお、本実施形態においては、被検体である患者の診断対象部位は脳であり、３次元画像撮影装置２はＣＴ装置２ＡおよびＭＲＩ装置２Ｂである。そして、ＣＴ装置２Ａにおいて、被検体の脳を含む３次元のＣＴ画像Ｂｃ０を生成し、ＭＲＩ装置２Ｂにおいて、被検体の脳を含む少なくとも１つの３次元のＭＲＩ画像を生成する。なお、本実施形態においては、ＭＲＩ画像は、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の２種類の画像が生成されるものとする。また、本実施形態においては、ＣＴ画像Ｂｃ０は、造影剤を使用しないで撮影を行うことにより取得される非造影ＣＴ画像とするが、造影剤を使用して撮影を行うことにより取得した造影ＣＴ画像を用いてもよい。なお、ＭＲＩ装置２ＢがＣＴ装置２Ａとは異なるモダリティであり、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２が異なるモダリティにより取得された少なくとも１つの脳画像である。なお、以降の説明において、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２をＭＲＩ画像Ｂｍ０で代表させる場合があるものとする。
【００２５】
画像保管サーバ３は、各種データを保存して管理するコンピュータであり、大容量外部記憶装置およびデータベース管理用ソフトウェアを備えている。画像保管サーバ３は、有線あるいは無線のネットワーク４を介して他の装置と通信を行い、画像データ等を送受信する。具体的には３次元画像撮影装置２で生成されたＣＴ画像およびＭＲＩ画像の画像データを含む各種データをネットワーク経由で取得し、大容量外部記憶装置等の記録媒体に保存して管理する。なお、画像データの格納形式およびネットワーク４経由での各装置間の通信は、ＤＩＣＯＭ（Digital Imaging and Communication in Medicine）等のプロトコルに基づいている。
【００２６】
梗塞領域判別装置１は、１台のコンピュータに、本開示の学習プログラムおよび梗塞領域判別プログラムをインストールしたものである。コンピュータは、診断を行う医師が直接操作するワークステーションまたはパーソナルコンピュータでもよく、それらとネットワークを介して接続されたサーバコンピュータでもよい。学習プログラムおよび梗塞領域判別プログラムは、ＤＶＤ（Digital Versatile Disc）あるいはＣＤ−ＲＯＭ（Compact Disc Read Only Memory）等の記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、ネットワークに接続されたサーバコンピュータの記憶装置、もしくはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じて医師が使用するコンピュータにダウンロードされ、インストールされる。
【００２７】
図２は、コンピュータに学習プログラムおよび梗塞領域判別プログラムをインストールすることにより実現される、本実施形態による梗塞領域判別装置の概略構成を示す図である。図２に示すように、梗塞領域判別装置１は、標準的なワークステーションの構成として、ＣＰＵ（Central Processing Unit）１１、メモリ１２およびストレージ１３を備えている。また、梗塞領域判別装置１には、液晶ディスプレイ等の表示部１４、並びにキーボードおよびマウス等の入力部１５が接続されている。
【００２８】
ストレージ１３には、ハードディスクドライブ等からなり、ネットワーク４を経由して画像保管サーバ３から取得した、被検体の医用画像、並びに処理に必要な情報を含む各種情報が記憶されている。
【００２９】
また、メモリ１２には、学習プログラムおよび梗塞領域判別プログラムが記憶されている。学習プログラムは、ＣＰＵ１１に実行させる処理として、脳梗塞を発症している被検体の脳のＣＴ画像Ｂｃ０、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２を取得する画像取得処理、少なくとも拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２に基づいて、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２から発症からの経過時間に応じた梗塞領域を抽出する梗塞領域抽出処理、ＣＴ画像Ｂｃ０と拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方との位置合わせを行う位置合わせ処理、位置合わせの結果に基づいて、ＣＴ画像Ｂｃ０における梗塞領域を発症からの経過時間に応じて特定する梗塞領域特定処理、ＣＴ画像Ｂｃ０において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像Ｂｃ１における発症からの経過時間に応じた梗塞領域を判別する判別器の学習を行う学習処理を規定する。また、梗塞領域判別プログラムは、ＣＰＵ１１に実行させる処理として、判別対象のＣＴ画像Ｂｃ１を取得する画像取得処理、判別対象のＣＴ画像Ｂｃ１における梗塞領域を発症からの経過時間に応じて判別する判別処理、および判別結果を表示部１４に表示する表示制御処理を規定する。
【００３０】
そして、ＣＰＵ１１がプログラムに従いこれらの処理を実行することで、コンピュータは、画像取得部２１、梗塞領域抽出部２２、位置合わせ部２３、梗塞領域特定部２４、学習部２５、判別器２６および表示制御部２７として機能する。ここで、画像取得部２１、梗塞領域抽出部２２、位置合わせ部２３、梗塞領域特定部２４および学習部２５が、本実施形態の判別器の学習装置を構成する。また、判別器２６および表示制御部２７が、本実施形態の梗塞領域判別装置を構成する。
【００３１】
画像取得部２１は、判別器２６の学習のために、脳梗塞を発症している被検体の脳のＣＴ画像Ｂｃ０、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２を画像保管サーバ３から取得する。また、発症からの経過時間に応じた梗塞領域の判別のために、梗塞領域の判別対象のＣＴ画像Ｂｃ１を画像保管サーバ３から取得する。なお、ＣＴ画像Ｂｃ０、ＣＴ画像Ｂｃ１、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２が既にストレージ１３に記憶されている場合には、画像取得部２１は、ストレージ１３からＣＴ画像Ｂｃ０、ＣＴ画像Ｂｃ１、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２を取得するようにしてもよい。また、画像取得部２１は、後述する判別器２６の学習のために、多数の被検体についてのＣＴ画像Ｂｃ０、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２を取得する。
【００３２】
ここで、梗塞領域の発症からの経過時間について説明する。脳梗塞の病期は、発症後２４時間以内の超急性期、発症後１〜７日の急性期、発症後１〜４週間の亜急性期、それ以降の慢性期(陳旧性)がある。本実施形態においては、亜急性期と慢性期とをまとめて亜急性期等とし、発症からの経過時間として、超急性期、急性期および亜急性期等を判別するものとする。
【００３３】
梗塞領域抽出部２２は、少なくとも拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２から脳の梗塞領域を抽出する。なお、本実施形態においては、ＣＴ画像Ｂｃ０をも用いて梗塞領域を抽出する。図３は拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の例を示す図である。なお、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２は３次元画像であるが、ここでは説明のため、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の１つの断層面における２次元の断層画像を用いて説明する。なお、図３に示すように、拡散強調画像Ｂｍ１は、頭蓋骨が消去された脳実質のみを含む画像となる。一方、ＦＬＡＩＲ画像Ｂｍ２は頭蓋骨も含んでいる。また、図３に示す拡散強調画像Ｂｍ１は超急性期および急性期の梗塞領域を含み、ＦＬＡＩＲ画像Ｂｍ２は急性期の梗塞領域を含むものとする。図３に示すように、拡散強調画像Ｂｍ１においては、超急性期および急性期の梗塞領域は他の領域と比較して高い信号値（低い濃度）を有するものとなる。また、ＦＬＡＩＲ画像Ｂｍ２においては、急性期の梗塞領域は他の領域と比較して高い信号値（低い濃度）を有するものとなる。
【００３４】
なお、図３には現れていないが、ＦＬＡＩＲ画像Ｂｍ２においては、超急性期の梗塞領域は他の領域と同程度か低い信号値（高い濃度）を有するものとなる。また、拡散強調画像Ｂｍ１においては、亜急性期等の梗塞領域は他の領域と同程度か低い信号値（高い濃度）を有するものとなる。また、ＣＴ画像Ｂｃ０においては、亜急性期等の梗塞領域は他の領域よりも明らかに低い信号値（高い濃度）を有するものとなる。
【００３５】
図４は、ＣＴ画像、拡散強調画像およびＦＬＡＩＲ画像の信号値と梗塞領域の発症からの経過時間との関係を示す図である。梗塞領域抽出部２２は、ＣＴ画像Ｂｃ０、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の信号値に基づいて、図４に示す関係を参照して、発症からの経過時間に応じた梗塞領域を抽出する。すなわち、梗塞領域抽出部２２は、拡散強調画像Ｂｍ１において高信号かつＦＬＡＩＲ画像Ｂｍ２において等信号となる領域を、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方から超急性期の梗塞領域として抽出する。また、梗塞領域抽出部２２は、拡散強調画像Ｂｍ１において高信号かつＦＬＡＩＲ画像Ｂｍ２において高信号となる領域を、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方から急性期の梗塞領域として抽出する。また、梗塞領域抽出部２２は、拡散強調画像Ｂｍ１において低信号かつＣＴ画像Ｂｃ０において低信号となる領域を、拡散強調画像Ｂｍ１から亜急性期等の梗塞領域として抽出する。なお、図４に示す関係はストレージ１３に保存しておけばよい。
【００３６】
なお、梗塞領域抽出部２２は、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方から発症からの経過時間に応じた梗塞領域を抽出するように学習がなされた判別器を用いて、発症からの経過時間に応じた梗塞領域を抽出するようにしてもよい。これにより、図３に示す梗塞領域Ａ１および梗塞領域Ａ２の少なくとも一方が抽出される。
【００３７】
位置合わせ部２３は、ＣＴ画像Ｂｃ０と拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方との位置合わせを行う。図５は、ＣＴ画像Ｂｃ０と拡散強調画像Ｂｍ１との位置合わせを説明するための図である。なお、ＣＴ画像Ｂｃ０および拡散強調画像Ｂｍ１はともに３次元画像であるが、ここでは説明のため、ＣＴ画像Ｂｃ０および拡散強調画像Ｂｍ１の対応する１つの断層面における２次元の断層画像を用いて説明する。図５に示すように、同一被検体においては、脳の形状は略同一となる。一方、拡散強調画像Ｂｍ１においては、超急性期および急性期の梗塞領域が他の領域と比較して高い信号値（低い濃度）を有するものとなる。これに対して、ＣＴ画像Ｂｃ０においては、超急性期および急性期の梗塞領域と他の領域とは信号値の差異が拡散強調画像Ｂｍ１ほど大きくない。なお、ＣＴ画像Ｂｃ０は、拡散強調画像Ｂｍ１とは異なり、頭蓋骨および脳実質を含んでいる。このため、位置合わせ部２３は、ＣＴ画像Ｂｃ０から脳実質の領域を脳領域として抽出し、抽出した脳領域と拡散強調画像Ｂｍ１との位置合わせを行う。
【００３８】
本実施形態において、位置合わせ部２３は、非剛体位置合わせにより、ＣＴ画像Ｂｃ０および拡散強調画像Ｂｍ１のいずれか一方を他方に対して位置合わせする。本実施形態においては、ＣＴ画像Ｂｃ０を拡散強調画像Ｂｍ１に対して位置合わせしているが、拡散強調画像Ｂｍ１をＣＴ画像Ｂｃ０に対して位置合わせしてもよい。
【００３９】
なお、非剛体位置合わせとしては、例えばＢスプラインおよびシンプレートスプライン等の関数を用いて、ＣＴ画像Ｂｃ０における特徴点をＭＲＩ画像Ｂｍ０における特徴点に対応する対応点に非線形に変換することによる手法を用いることができるが、これに限定されるものではない。
【００４０】
なお、ＦＬＡＩＲ画像Ｂｍ２を用いる場合、位置合わせ部２３は、拡散強調画像Ｂｍ１の場合と同様に、ＣＴ画像Ｂｃ０とＦＬＡＩＲ画像Ｂｍ２との位置合せを行う。
【００４１】
梗塞領域特定部２４は、位置合わせ部２３による位置合わせの結果に基づいて、ＣＴ画像Ｂｃ０における発症からの経過時間に応じた梗塞領域を特定する。図６はＣＴ画像Ｂｃ０における発症からの経過時間に応じた梗塞領域の特定を説明するための図である。図６に示すように、梗塞領域特定部２４は、ＭＲＩ画像Ｂｍ０において抽出した梗塞領域Ａ１またはＦＬＡＩＲ画像Ｂｍ２において抽出した梗塞領域Ａ２に対応するＣＴ画像Ｂｃ０のボクセルの領域を梗塞領域Ａ３に特定する。なお、特定した梗塞領域Ａ３には、発症からの経過時間に応じたラベルが付与される。図６においては、「急性期」のラベルが付与された状態を示している。
【００４２】
学習部２５は、ＣＴ画像Ｂｃ０において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像における発症からの経過時間に応じた梗塞領域を判別する判別器２６の学習を行う。本実施形態においては、判別器２６は、判別対象となるＣＴ画像Ｂｃ１が入力されると、ＣＴ画像Ｂｃ１における発症からの経過時間に応じた梗塞領域を判別する。具体的には、判別器２６は、判別対象となるＣＴ画像Ｂｃ１の各ボクセルを、梗塞領域であるか梗塞領域以外の領域であるか、さらには梗塞領域の場合は発症からの経過時間に応じた１以上の梗塞領域のいずれであるかを判別する。発症からの経過時間に応じた１以上の梗塞領域のいずれであるか判別するとは、梗塞領域が、超急性期、急性期および亜急性期等のいずれであるかを判別することを意味する。このため、学習部２５は、複数の被検体のＣＴ画像Ｂｃ０において特定した梗塞領域Ａ３から、あらかじめ定められたサイズ（例えば３×３等）の領域内の特徴量を取得し、取得した特徴量を判別器２６に入力し、発症からの経過時間に応じた梗塞領域の判別結果を出力するように、判別器２６の学習、すなわち機械学習を行う。
【００４３】
このように学習が行われることにより、ＣＴ画像Ｂｃ１が入力されると、ＣＴ画像Ｂｃ１のボクセルを、発症からの経過時間に応じた１以上の梗塞領域と梗塞領域外の領域とに分類して、判別対象のＣＴ画像Ｂｃ１における発症からの経過時間に応じた梗塞領域を判別する判別器２６が生成される。
【００４４】
表示制御部２７は、判別対象となるＣＴ画像Ｂｃ１の判別器２６による判別結果を表示部１４に表示する。図７は判別結果の表示の例を示す図である。なお、図７においては、判別対象となるＣＴ画像Ｂｃ１の１つの断面における断層画像を示している。図７に示すように、表示部１４に表示された判別結果においては、判別対象のＣＴ画像Ｂｃ１において、梗塞領域Ａ４が特定されている。また、梗塞領域Ａ４の発症からの経過時間である「急性期」のラベルＬ０が表示されている。
【００４５】
なお、機械学習の手法としては、サポートベクタマシン（ＳＶＭ(Support Vector Machine)）、ディープニューラルネットワーク（ＤＮＮ(Deep Neural Network)）、畳み込みニューラルネットワーク（ＣＮＮ(Convolutional Neural Network)）、およびリカレントニューラルネットワーク（ＲＮＮ(Recurrent Neural Network)）等を用いることができる。
【００４６】
次いで、本実施形態において行われる処理について説明する。図８は本実施形態において学習時に行われる処理を示すフローチャートである。まず、画像取得部２１が、脳梗塞を発症している被検体の脳のＣＴ画像Ｂｃ０およびＭＲＩ画像Ｂｍ０（すなわち拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２）を取得し（ステップＳＴ１）、梗塞領域抽出部２２が、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２から発症からの経過時間に応じた梗塞領域Ａ１，Ａ２をそれぞれ抽出する（ステップＳＴ２）。次いで、位置合わせ部２３が、ＣＴ画像Ｂｃ０と拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方との位置合わせを行い（ステップＳＴ３）、梗塞領域特定部２４が、位置合わせの結果に基づいて、ＣＴ画像Ｂｃ０における発症からの経過時間に応じた梗塞領域Ａ３を特定する（ステップＳＴ４）。そして、学習部２５が、ＣＴ画像Ｂｃ０において特定された、発症からの経過時間に応じた梗塞領域を教師データとして、入力されたＣＴ画像Ｂｃ１における梗塞領域を判別する判別器２６を学習し（ステップＳＴ６）、処理を終了する。
【００４７】
図９は本実施形態において発症からの経過時間に応じた梗塞領域の判別時に行われる処理を示すフローチャートである。まず、画像取得部２１が、判別対象となるＣＴ画像Ｂｃ１を取得し（ステップＳＴ１１）、判別器２６が判別対象のＣＴ画像における、発症からの経過時間に応じた梗塞領域を判別する（ステップＳＴ１２）。そして、表示制御部２７が、判別結果を表示部１４に表示し（ステップＳＴ１３）、処理を終了する。
【００４８】
このように、本実施形態においては、脳梗塞を発症している被検体の脳のＣＴ画像Ｂｃ０、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２を取得し、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２から発症からの経過時間に応じた梗塞領域を抽出し、ＣＴ画像Ｂｃ０と拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方との位置合わせを行う。そして、位置合わせの結果に基づいて、ＣＴ画像Ｂｃ０における発症からの経過時間に応じた梗塞領域Ａ３を特定し、特定された梗塞領域Ａ１を教師データとして、入力されたＣＴ画像Ｂｃ１における、発症からの経過時間に応じた梗塞領域を判別する判別器２６の学習を行うようにした。
【００４９】
ここで、ＣＴ画像においては、出血領域が他の領域とは大きく異なる信号値を有するため、ＣＴ画像において出血領域を特定することは容易である。しかしながら、ＣＴ画像においては、梗塞領域は他の領域と異なる信号値を有するものの、その差は出血領域と他の領域との差よりも小さい。一方、ＭＲＩ画像においては、梗塞領域が他の領域とは大きく異なる信号値を有する。特に、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２においては、発症からの経過時間に応じて信号値が異なるため、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２、さらにはＣＴ画像Ｂｃ０をも用いることにより、梗塞領域の発症からの経過時間を特定することができる。このため、脳梗塞を発症している被検体の脳の拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２の少なくとも一方とＣＴ画像Ｂｃ０とを位置合わせすれば、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２における梗塞領域Ａ１，Ａ２に基づいて、ＣＴ画像Ｂｃ０における発症からの経過時間に応じた梗塞領域Ａ３を特定することができる。したがって、特定された発症からの経過時間に応じた梗塞領域Ａ３を教師データとして判別器２６の学習を行うことにより、学習された判別器２６によって、判別対象となるＣＴ画像Ｂｃ１における、発症からの経過時間に応じた梗塞領域を判別できることとなる。これにより、ＣＴ画像のみを用いて、脳の出血領域のみならず、発症からの経過時間に応じた梗塞領域を判別することができる。よって、本実施形態によれば、ＣＴ画像のみを用いて、発症からの経過時間に応じた脳梗塞の診断を迅速に行うことができる。
【００５０】
なお、上記実施形態においては、ＭＲＩ画像Ｂｍ０として拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２を用いているが、拡散強調画像Ｂｍ１およびＦＬＡＩＲ画像Ｂｍ２以外のＭＲＩ画像を用いてもよい。例えば、Ｔ１強調画像およびＴ２強調画像等を用いてもよい。また、拡散強調画像Ｂｍ１、ＦＬＡＩＲ画像Ｂｍ２、Ｔ１強調画像およびＴ２強調画像等から選択された１以上の画像を用いてもよい。これらのＭＲＩ画像の組み合わせにより、亜急性期の梗塞領域と慢性期の梗塞領域とを特定することも可能となる。具体的には、拡散強調画像Ｂｍ１において等〜低信号であり、Ｔ２強調画像において高信号の場合、亜急性期の梗塞領域であると特定することができる。また、拡散強調画像Ｂｍ１において低信号であり、Ｔ２強調画像において高信号の場合、慢性期の梗塞領域であると特定することができる。したがって、さらにこのような関係を用いて判別器２６の学習を行うことにより、超急性期、急性期、亜急性期および慢性期の梗塞領域を判別することが可能となる。
【００５１】
また、上記実施形態においては、梗塞領域の発症からの経過時間を特定するためにＣＴ画像Ｂｃ０を用いているが、ＭＲＩ画像のみを用いてもよい。
【００５２】
また、上記実施形態においては、判別器２６の学習に用いるＣＴ画像Ｂｃ０として、非造影ＣＴ画像または造影ＣＴ画像を用いているが、造影ＣＴ画像および非造影ＣＴ画像の双方を判別器２６の学習に用いるようにしてもよい。このように学習された判別器２６を用いることにより、判別対象のＣＴ画像が造影ＣＴ画像および非造影ＣＴ画像のいずれであっても、梗塞領域を判別できることとなる。
【００５３】
また、上記実施形態においては、梗塞領域判別装置１が学習装置を内包しているが、これに限定されるものではない。すなわち、診断支援システムにおいて、梗塞領域判別装置１とは別個に、画像取得部２１、梗塞領域抽出部２２、位置合わせ部２３、梗塞領域特定部２４、および学習部２５を備え、判別器２６の学習を行う学習装置を備えるようにしてもよい。この場合、梗塞領域判別装置１は、画像取得部２１、判別器２６および表示制御部２７のみを備えるものとなる。
【００５４】
また、上記実施形態において、例えば、画像取得部２１、梗塞領域抽出部２２、位置合わせ部２３、梗塞領域特定部２４、学習部２５、判別器２６および表示制御部２７といった各種の処理を実行する処理部（processing unit）のハードウェア的な構造としては、次に示す各種のプロセッサ（processor）を用いることができる。上記各種のプロセッサには、上述したように、ソフトウェア（プログラム）を実行して各種の処理部として機能する汎用的なプロセッサであるＣＰＵに加えて、ＦＰＧＡ（Field  Programmable Gate Array）等の製造後に回路構成を変更可能なプロセッサであるプログラマブルロジックデバイス（Programmable Logic Device :PLD）、ＡＳＩＣ（Application Specific Integrated Circuit）等の特定の処理を実行させるために専用に設計された回路構成を有するプロセッサである専用電気回路等が含まれる。
【００５５】
１つの処理部は、これらの各種のプロセッサのうちの１つで構成されてもよいし、同種または異種の２つ以上のプロセッサの組み合わせ（例えば、複数のＦＰＧＡの組み合わせまたはＣＰＵとＦＰＧＡとの組み合わせ）で構成されてもよい。また、複数の処理部を１つのプロセッサで構成してもよい。
【００５６】
複数の処理部を１つのプロセッサで構成する例としては、第１に、クライアントおよびサーバ等のコンピュータに代表されるように、１つ以上のＣＰＵとソフトウェアとの組み合わせで１つのプロセッサを構成し、このプロセッサが複数の処理部として機能する形態がある。第２に、システムオンチップ（System On Chip:SoC）等に代表されるように、複数の処理部を含むシステム全体の機能を１つのＩＣ（Integrated Circuit）チップで実現するプロセッサを使用する形態がある。このように、各種の処理部は、ハードウェア的な構造として、上記各種のプロセッサの１つ以上を用いて構成される。
【００５７】
さらに、これらの各種のプロセッサのハードウェア的な構造としては、より具体的には、半導体素子等の回路素子を組み合わせた電気回路（circuitry）を用いることができる。
【符号の説明】
【００５８】
１    梗塞領域判別装置
２    ３次元画像撮影装置
３    画像保管サーバ
４    ネットワーク
１１    ＣＰＵ
１２    メモリ
１３    ストレージ
１４    ディスプレイ
１５    入力部
２１    画像取得部
２２    梗塞領域抽出部
２３    位置合わせ部
２４    梗塞領域特定部
２５    学習部
２６    判別器
２７    表示制御部
Ａ１，Ａ２，Ａ３    梗塞領域
Ｂｃ０    ＣＴ画像
Ｂｃ１    判別対象のＣＴ画像
Ｂｍ１    拡散強調画像
Ｂｍ２    ＦＬＡＩＲ画像
【図１】
2020032045.tif 000003
【図２】
2020032045.tif 000004
【図３】
2020032045.tif 000005
【図４】
2020032045.tif 000006
【図５】
2020032045.tif 000007
【図６】
2020032045.tif 000008
【図７】
2020032045.tif 000009
【図８】
2020032045.tif 000010
【図９】
2020032045.tif 000011