JP7078948B2 - 眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体 - Google Patents

眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体 Download PDF

Info

Publication number
JP7078948B2
JP7078948B2 JP2017124825A JP2017124825A JP7078948B2 JP 7078948 B2 JP7078948 B2 JP 7078948B2 JP 2017124825 A JP2017124825 A JP 2017124825A JP 2017124825 A JP2017124825 A JP 2017124825A JP 7078948 B2 JP7078948 B2 JP 7078948B2
Authority
JP
Japan
Prior art keywords
data
information processing
papilla
processing system
learning
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
JP2017124825A
Other languages
English (en)
Other versions
JP2019005319A (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.)
Tohoku University NUC
Topcon Corp
RIKEN Institute of Physical and Chemical Research
Original Assignee
Tohoku University NUC
Topcon Corp
RIKEN Institute of Physical and Chemical Research
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 Tohoku University NUC, Topcon Corp, RIKEN Institute of Physical and Chemical Research filed Critical Tohoku University NUC
Priority to JP2017124825A priority Critical patent/JP7078948B2/ja
Priority to US15/925,922 priority patent/US10765316B2/en
Publication of JP2019005319A publication Critical patent/JP2019005319A/ja
Application granted granted Critical
Publication of JP7078948B2 publication Critical patent/JP7078948B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/14Arrangements specially adapted for eye photography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0016Operational features thereof
    • A61B3/0025Operational features thereof characterised by electronic signal processing, e.g. eye models
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/102Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • A61B3/1225Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation
    • A61B3/1233Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation for measuring blood flow, e.g. at the retina
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10101Optical tomography; Optical coherence tomography [OCT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30041Eye; Retina; Ophthalmic

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Ophthalmology & Optometry (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Signal Processing (AREA)
  • Hematology (AREA)
  • Eye Examination Apparatus (AREA)

Description

特許法第30条第2項適用 生体医用画像研究会 第4回若手発表会、開催日 平成29年3月25日 生体医用画像研究会 第4回若手発表会 プログラム・抄録集、発行日 平成29年3月 電子情報通信学会、開催日 平成29年5月26日 電子情報通信学会 プログラム・講演論文詳細、ウェブサイトの掲載日 平成29年3月
本発明は、眼科情報を処理するためのシステム、方法、プログラム、及び記録媒体に関する。
主要な眼疾患の一つに緑内障がある。緑内障は、視神経損傷及び視力障害を引き起こす一群の進行性疾患である。緑内障は多因子疾患であり、病態の細分化や治療の個別化が望まれている。そのため、緑内障の診断には、眼圧検査、隅角検査、眼底検査(視神経乳頭陥凹の度合い、網膜神経線維層の欠損、乳頭出血の有無、乳頭周囲網脈絡膜萎縮の有無)、視野検査、画像解析(乳頭形状解析、層厚解析、篩状板解析)など、多様な検査が利用される。
緑内障の病態分類として視神経乳頭の形状に応じた分類(乳頭形状分類)があり、その代表的なものにニコレラ(Nicolela)らによる分類がある。ニコレラ分類では、視神経乳頭の形状に応じた4つの類型(局所虚血型、近視型、加齢性硬化型、全体的拡大型)が定められている。局所虚血型(Focal Ischemia、 FI)は、リムの一部にノッチが存在し、視神経線維層の局所的な欠損が見られる類型であり、女性に多く、片頭痛や発作を伴うことが多い。近視型(Myopic、MY)は、視神経乳頭が傾斜しており、耳側に陥凹を伴う三日月型の乳頭周囲網脈絡膜萎縮(PPA)が見られる類型であり、若年層に多く、近視を伴うことが多い。加齢性硬化型(Senile Sclerotic、SS)は、視神経乳頭が円形であり、陥凹が浅く、乳頭の周辺にハロ(Halo)が見られる類型であり、高齢者に多く、心血管障害を伴うことが多い。全体的拡大型(Generalized Enlargement、GE)は、大きく深い円形の陥凹を呈する類型であり、眼圧が高いことが多い。
乳頭形状分類は、典型的には、眼底カメラや走査型レーザー検眼鏡(SLO)により取得された眼底像を医師が読影することで行われる。つまり、分類は医師の主観によって行われる。一方、緑内障の有病率が比較的高いことや、早期発見が重要であることを鑑み、分類の自動化が望まれている。
しかし、眼底カメラやSLOで得られる眼底像は、撮影条件などに画質が大きく左右されるため、乳頭形状パラメータの定量化が難しく、分類の自動化を阻む要因の一つになっている。
また、従来の分類では眼底像のみ(つまり乳頭形状パラメータのみ)を参照しているため、良好な精度及び確度で分類を行うことが難しい。一方、分類のために参照可能な眼科情報には多種多様なものがあるが、それら全てを考慮すると処理が複雑かつ膨大となり、処理システムの規模やリソースも極めて大きくなるため、実用化が困難な場合も想定される。
特開平10-155744号公報 特開2000-245700号公報 特開2005-253796号公報 特開2006-263127号公報 特開2007-252707号公報 特開2008-220617号公報 特開2008-237839号公報 特開2008-237840号公報 特開2011-83555号公報 特開2016-2380号公報 特開2016-187475号公報
Nicolela MT, Drance SM, "Various glaucomatous optic nerve appearances: clinical correlations", Ophthalmology, Volume 103, Issue 4, April 1996, Pages 640-649
本発明の目的の一つは、良好な精度及び確度での乳頭形状分類を自動で行うための技術を提供することにある。
本発明の目的の他の一つは、大規模なシステムを用いることなく、良好な精度及び確度での乳頭形状分類を自動で行うための技術を提供することにある。
例示的な実施形態に係る眼科情報処理システムは、被検者の視神経乳頭の形態データと当該被検者の背景データとを受け付ける受付部と、前記形態データと同種のデータ群及び前記背景データと同種のデータ群に基づき機械学習を実行することにより生成された学習情報を予め記憶する記憶部と、前記学習情報と前記形態データと前記背景データとに基づいて、前記機械学習に用いられた人工知能エンジンと同種の人工知能エンジンを適用することにより、視神経乳頭形状分類を実行するプロセッサとを含む。
例示的な実施形態によれば、良好な精度及び確度での乳頭形状分類を自動で行うことが可能である。
例示的な実施形態によれば、大規模なシステムを用いることなく、良好な精度及び確度での乳頭形状分類を自動で行うことが可能である。
例示的な実施形態に係る眼科情報処理システムの構成の一例を表す概略図である。 例示的な実施形態に係る眼科情報処理システムにおいて適用可能なパラメータのリストの一例である。 例示的な実施形態に係る眼科情報処理システムを用いて実施可能な使用形態の一例を表すフローチャートである。 例示的な実施形態に係る眼科情報処理システムにより表示可能な分類結果の一例である。
例示的な実施形態に係る眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体について以下に説明する。例示的な眼科情報処理方法は、例示的な眼科情報処理システムによって実現される。例示的な眼科情報処理システムは、相互に通信が可能な2以上の装置(例えば、1以上のコンピュータ、1以上の記憶装置など)を含んでよい。或いは、例示的な眼科情報処理システムは、単一の装置(例えば、記憶装置が搭載されたコンピュータ)であってもよい。
例示的な眼科情報処理方法を実現するためのハードウェア及びソフトウェアは、以下に例示する眼科情報処理システムには限定されず、その実現に寄与する任意のハードウェアと任意のソフトウェアとの組み合わせが含まれていてよい。典型的な例として、眼科情報処理システムは、人工知能エンジンとして機能するハードウェア及びソフトウェアを含んでいてよい。
例示的な眼科情報処理方法は、医学的な検査を行うために用いられ、より具体的には、眼科情報を処理して視神経乳頭の形状を自動で分類するために用いられる。乳頭形状分類には、被検者の眼底の形態データと、当該被検者の背景データとが用いられる。
眼底形態データは、眼底の形態(形状、サイズ等)を表すデータである。眼底形態データは、多種多様な眼科情報から選択された種別のパラメータに関する。パラメータの選択は、典型的には、人工知能等のコンピュータを利用して実行される。或いは、医師等の専門家がパラメータの選択を行ってもよい。また、コンピュータにより選択されたパラメータと専門家により選択されたパラメータとの双方を考慮してもよい。眼底形態データの種別の具体例やパラメータの選択処理の具体例については後述する。
背景データは、被検者に関するデータであって、典型的には、医学的検査以外の方法で取得されたデータを含む。また、背景データは、電子カルテ等に記録(又は付帯)された情報を含んでよい。例えば、電子カルテには、医学的検査以外の方法で取得された患者情報(年齢、性別など)に加え、医学的検査により取得された検査データも記録される。検査データとしては、視力値、屈折力値、等価球面度数、眼圧値、角膜厚値、視野検査結果(Mean Defect値、Pattern Standard Deviation値など)などがある。背景データは、各種情報から選択された種別のパラメータに関する。眼底形態データの場合と同様に、パラメータの選択は、例えば、コンピュータ及び専門家の少なくとも一方によって実行される。背景データの種別の具体例やパラメータの選択処理の具体例については後述する。
例示的なプログラムは、例示的な眼科情報処理システム等のコンピュータに、例示的な眼科情報処理方法を実行させる。また、例示的な記録媒体は、コンピュータ読み取り可能な記録媒体であって、例示的なプログラムを記録する。例示的な記録媒体は、非一時的な(non-transitory)記録媒体である。例示的な記録媒体は、磁気、光、光磁気、半導体などを利用した電子媒体であってよい。典型的には、例示的な記録媒体は、磁気テープ、磁気ディスク、光ディスク、光磁気ディスク、フラッシュメモリ、ソリッドステートドライブなどである。
〈例示的な眼科情報処理システム〉
眼科情報処理システムについて説明する。図1に示す例示的な眼科情報処理システム1は、被検者の眼底形態データ及び背景データに基づいて、当該被験者の視神経乳頭形状分類を実行する。視神経乳頭形状分類は、典型的にはニコレラ分類であるが、他の分類法でもよい。
眼科情報処理システム1は、制御部10と、記憶部20と、データ処理部30と、データ受付部40と、通信部50とを含む。ユーザインターフェイス(UI)100は、眼科情報処理システム1に含まれてもよいし、含まれなくてもよい。眼科情報処理システム1に含まれる要素は、単一の装置又は2以上の装置として構成される。例えば、眼科情報処理システム1は、これら要素の全てを備えた単一のコンピュータを含む。
眼科情報処理システム1が2以上の装置を含む場合の例として、制御部10を含むコンピュータと、記憶部20を含むコンピュータと、データ処理部30を含むコンピュータとを、それぞれ個別に設けることができる。或いは、眼科情報処理システム1は、制御部10、記憶部20及びデータ処理部30のうちのいずれか2つを備えたコンピュータと、他の1つを備えたコンピュータとを含む。異なるコンピュータの間の通信形態は、有線通信及び/又は無線通信を含んでよく、専用回線及び/又は公衆回線を含んでよく、ローカルエリアネットワーク(LAN)、ワイドエリアネットワーク(WAN)、近距離通信及びインターネットのうちの少なくとも1つを含んでよい。
本実施形態において、各種の処理はコンピュータにより実行される。コンピュータは、人工知能エンジンとして機能するよう構成されてもよい。コンピュータは、1以上のプロセッサを含む。プロセッサは、例えば、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、ASIC(Application Specific Integrated Circuit)、プログラマブル論理デバイス(例えば、SPLD(Simple Programmable Logic Device)、CPLD(Complex Programmable Logic Device)、FPGA(Field Programmable Gate Array))等の回路を意味する。プロセッサは、例えば、記憶装置(記憶回路)に格納されているプログラムを読み出して実行することにより所望の機能を実現する。
プロセッサは、記憶装置及び出力装置を制御することができる。記憶装置は、コンピュータの内部に含まれてもよいし、コンピュータの外部に設けられてもよい。出力装置についても同様である。
眼科情報処理システム1は、OCTシステム200、カルテ入力装置250、及びデータベース(DB)300の少なくともいずれかを含んでよい。OCTシステム200及びカルテ入力装置250は、病院内LAN等のネットワークを介してデータベース(DB)300に接続されている。
〈OCTシステム200〉
OCTシステム200は、光コヒーレンストモグラフィ(OCT)を使用して被検眼(眼底)のデータを収集するデータ収集機能を備えてよい。データ収集機能は、例えば、公知のフーリエドメインOCT(周波数ドメインOCT)の手法によって実現される。データ収集機能を実現する構成(OCT装置)は、従来のフーリエドメインOCTと同様に、光源と、干渉光学系と、光検出器と、OCTデータ構築プロセッサとを含む。
スペクトラルドメインOCTの場合、OCT装置は、以下の要素を含む:低コヒーレンス光源;低コヒーレンス光源から出力された低コヒーレンス光を測定光と参照光とに分割し、被検眼からの測定光の戻り光と参照光路を経由した参照光とを重ね合わせる干渉光学系;干渉計により生成された干渉光を空間分割(スペクトル分解)してスペクトルデータを検出する分光器;分光器からの出力信号にフーリエ変換等の信号処理を適用してOCT画像データを構築するOCTデータ構築プロセッサ。
スウェプトソースOCTの場合、OCT装置は、以下の要素を含む:波長可変光源;波長可変光源から時分割で出力される光を測定光と参照光とに分割し、被検眼からの測定光の戻り光と参照光路を経由した参照光とを重ね合わせる干渉光学系;干渉計により生成された干渉光を検出する光検出器(バランスドフォトダイオード等);光検出器からの出力信号を収集するデータ収集システム(DAQ又はDAS);データ収集システムにより収集されたスペクトルデータにフーリエ変換等の信号処理を適用してOCT画像データを構築するOCTデータ構築プロセッサ。
OCTシステム200は、OCT画像データを解析する解析機能を備えてよい。解析機能は、OCTを使用して収集された眼底のデータ(OCT画像データ)を解析して眼底形態データを生成する解析プロセッサによって実現される。解析プロセッサは、例えば、コンピュータ又は専門家により事前に選択された1以上の種別のパラメータの一部又は全部について、眼底形態データを生成する。解析プロセッサは、コンピュータ又は専門家により事前に選択された種別以外のパラメータについて、眼底形態データを生成するように構成されてもよい。解析プロセッサが実行する解析処理は、眼科分野(より広くは医学分野、更に広くは任意分野)において使用される公知の処理を含む。
〈カルテ入力装置250〉
カルテ入力装置250は、電子カルテの作成に用いられるコンピュータである。典型的には、カルテ入力装置250は、医師が使用するコンピュータ端末、入力データからカルテ記載事項を読み取るプロセッサ、用紙等に記入された情報を読み取るデータスキャナなどであってよい。
〈データベース300〉
データベース300は各種データを管理する。データベース300は、例えば、病院情報システム(HIS)と、医療用画像管理システム(PACS)とを含む。病院情報システムには、電子カルテシステムが含まれる。OCTシステム200又は他のデータ取得装置(検査装置、撮影装置等)により取得されたデータは、電子カルテシステム又は医療用画像管理システムに保存される。カルテ入力装置250等のコンピュータを用いて作成された電子カルテは、電子カルテシステムに保存される。
データベース300は、健康診断や検診等のスクリーニングによって取得されたデータを保存するスクリーニングデータベースを含んでもよい。
〈制御部10〉
制御部10は各種の制御を実行する。例えば、制御部10は、眼科情報処理システム1の各要素の制御や、2以上の要素の連係的制御を実行する。また、制御部10は、眼科情報処理システム1の外部の装置の制御を行うことができる。例えば、ユーザインターフェイス100が眼科情報処理システム1に含まれない場合、制御部10はユーザインターフェイス100の制御を行うことができる。制御部10はプロセッサを含む。
〈出力制御部11〉
出力制御部11は、各種の出力装置を制御する。出力装置は、例えば、表示装置、音声出力装置、印刷装置、通信機器、及びデータライタのうちの少なくとも1つを含む。すなわち、出力制御部11は、次の制御のうちの少なくともいずれかを実行する:情報を表示するための表示装置の制御;音声情報を出力するための音声出力装置の制御;用紙に情報を印刷するための印刷装置の制御;外部装置に情報を送信させるための通信機器の制御;記録媒体に情報を記録するためのデータライタの制御。
図1に示す表示部101は表示装置に相当する。また、通信部50は通信機器に相当する。図示は省略するが、音声出力装置、印刷装置、及びデータライタのうちの1つ以上が設けられていてよい。なお、他の種別の出力装置が設けられていてもよい。
〈記憶部20〉
記憶部20は各種のデータを記憶する。記憶部20に記憶されるデータの例として、データ受付部40により外部から受け付けられたデータがある。典型的には、記憶部20は、眼底形態データと背景データとを記憶する。記憶部20は、例えば、半導体記憶装置、磁気記憶装置、光記憶装置、及び光磁気記憶装置のうちの少なくとも1つを含む。
〈学習情報21〉
記憶部20には学習情報21が予め記憶されている。学習情報21は、視神経乳頭形状分類において参照される。分類処理は、任意のプロセッサにより実行され、典型的には人工知能技術を利用して実行される。また、分類処理と同種の人工知能技術が学習においても適用される。
学習は、任意の機械学習技術を利用して実行される。機械学習では、公知技術と同様に、準備されたサンプルデータセットを解析することにより規則・ルール・知識表現・判断基準などを抽出し、アルゴリズムを発展させる。機械学習のアルゴリズムは任意であり、典型的には教師あり学習を適用することができるが、教師なし学習、半教師あり学習、強化学習、トランスダクション、マルチタスク学習などの技術を適用してもよい。
本実施形態では、データ処理部30に設けられた学習器32によって機械学習が実行される。他方、眼科情報処理システム1の外部に設けられたコンピュータ(学習器)によって機械学習を行うことができる。学習器32(又は外部学習器)は、所定の人工知能技術を利用して所定のサンプルデータセットに基づく機械学習を実行することにより分類モデルを生成する。生成された分類モデルは、学習情報21に含まれる。
本実施形態では、ニューラルネットワーク(NN)、勾配ブースティング決定木(GBDT)、サポートベクターマシン(SVM)、及び、ベイズ(例えばナイーブベイズ(NB))のうちのいずれかの人工知能技術が利用されるが、他の人工知能技術を利用することも可能である。本実施形態では、利用される人工知能技術に応じた種別のパラメータが学習に用いられる。
視神経乳頭形状分類にニューラルネットワークが利用される場合、陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、及び、上方乳頭周囲網膜神経線維平均層厚(これらのうちの少なくとも1つのパラメータでもよい)が、眼底形態データのパラメータとして分類のために参照される。更に、年齢及び等価球面度数(これらのうちの少なくとも1つのパラメータでもよい)が、背景データのパラメータとして分類のために参照される。また、分類のためのパラメータと同じ種別のパラメータが学習に用いられる。なお、等価球面度数の代わりに、これと同等の情報や類似の情報を用いることも可能である。例えば、被検眼の眼軸長や収差から等価球面度数を代替可能なパラメータを求め、これを背景データに含めることが可能である。
視神経乳頭形状分類に勾配ブースティング決定木が利用される場合、水平乳頭傾斜度、鼻側リム乳頭面積比、陥凹平均深さ、陥凹乳頭面積比、及び、上鼻側乳頭周囲網膜神経線維平均層厚(これらのうちの少なくとも1つのパラメータでもよい)が、眼底形態データのパラメータとして分類のために参照される。更に、年齢及び等価球面度数(これらのうちの少なくとも1つのパラメータでもよい)が、背景データのパラメータとして分類のために参照される。また、分類のためのパラメータと同じ種別のパラメータが学習に用いられる。
視神経乳頭形状分類にサポートベクターマシンが利用される場合、乳頭水平径、陥凹平均深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、及び、乳頭網膜色素上皮高低差(これらのうちの少なくとも1つのパラメータでもよい)が、眼底形態データのパラメータとして分類のために参照される。更に、年齢及び等価球面度数(これらのうちの少なくとも1つのパラメータでもよい)が、背景データのパラメータとして分類のために参照される。また、分類のためのパラメータと同じ種別のパラメータが学習に用いられる。
視神経乳頭形状分類にナイーブベイズが利用される場合、上方乳頭周囲網膜神経線維平均層厚、鼻側リム乳頭面積比、水平乳頭傾斜度、及び、陥凹平均深さ(これらのうちの少なくとも1つのパラメータでもよい)が、眼底形態データのパラメータとして分類のために参照される。更に、年齢及び等価球面度数(これらのうちの少なくとも1つのパラメータでもよい)が、背景データのパラメータとして分類のために参照される。また、分類のためのパラメータと同じ種別のパラメータが学習に用いられる。
人工知能技術の種別に応じたパラメータ(眼底形態データ、背景データ)の選択について説明する。発明者らは、OCT画像の読影・解析、カラー眼底像、及び患者背景データに基づいて医師たちが診断した広義開放隅角緑内障の患者のN個の眼(サンプル眼)に関するサンプルデータセットに基づき、ニューラルネットワーク学習器、勾配ブースティング決定木学習器、サポートベクターマシン学習器、及び、ナイーブベイズ学習器のそれぞれを用いて機械学習を実施した。
サンプルデータセットには、少なくとも、3次元OCTデータを解析して得られた眼底形態データが含まれる。例えば、3次元OCTデータに対してセグメンテーションが適用され、その結果から乳頭外縁と陥凹外縁とが検出され、これら検出結果から眼底形態データが生成される。このような眼底形態データとして、視神経乳頭形状パラメータや、乳頭周囲網膜神経線維層厚などがある。更に、サンプルデータセットは、レーザー散乱を利用した血流画像化装置であるレーザースペックルフローグラフィ(LSFG)により得られた血流動態データを含んでいてよい。
発明者らは、スウェプトソースOCTを利用して、22種類の視神経乳頭形状パラメータと、26種類の乳頭周囲網膜神経線維層厚(cpRNFLT)パラメータとを準備した。更に、レーザースペックルフローグラフィを利用して、36種類の血流パラメータを準備した。また、患者のカルテから7種類の背景パラメータを準備した。
発明者らが考慮したパラメータ(特徴量)のリストを図2に示す。このリストにおいて、No.1~7は7種類の背景パラメータであり、No.8~29は22種類の視神経乳頭形状パラメータであり、No.30~55は26種類の乳頭周囲網膜神経線維層厚パラメータであり、No.56~91は36種類の血流パラメータである。
更に、発明者らは、N個のサンプル眼をランダムにN1個のトレーニングデータとN2個のテストデータとに分け、トレーニングデータに基づき特徴選択及び機械学習によって分類モデルを作成し、その性能をテストデータに基づき確認した。
発明者らは、N個のサンプル眼についての医師たちによるニコレラ分類の結果を教師データとした。更に、ニューラルネットワーク、勾配ブースティング決定木、サポートベクターマシン、及びナイーブベイズのそれぞれについて教師あり学習を実施して分類モデルを構築し、その性能を評価した。このとき、91個のパラメータ(特徴量)のうちから最適なパラメータを選択するために特徴選択を行った。
なお、勾配ブースティング決定木では、特徴選択の手法として埋め込み特徴選択法が適用され、性能評価視標として10分割交差検証が適用された。また、ニューラルネットワーク、サポートベクターマシン、及びナイーブベイズでは、フィルター型特徴選択法の一つであるmRMR(minimum Redundancy Maximum Relevance)と、ラッパー型特徴選択法の一つである遺伝的アルゴリズム(GA)とを組み合わせて特徴選択を行い、10分割交差検証により性能評価を行った。
以上のような分類評価の結果、ニューラルネットワーク、勾配ブースティング決定木、サポートベクターマシン、及びナイーブベイズの全てにおいて、80パーセント以上の正解率を得た。なお、これら4つのうち最も正解率が高かったのはニューラルネットワークであった。
更に、最適なパラメータ種別として、ニューラルネットワークでは、陥凹面積(No.9)、陥凹平均深さ(No.18)、陥凹最大深さ(No.17)、鼻側リム乳頭面積比(No.22)、水平乳頭傾斜度(No.26)、上鼻側乳頭周囲網膜神経線維平均層厚(No.38)、上方乳頭周囲網膜神経線維平均層厚(No.32)、年齢(No.2)、及び等価球面度数(No.3)が選択された。
勾配ブースティング決定木では、最適なパラメータ種別として、水平乳頭傾斜度(No.26)、鼻側リム乳頭面積比(No.22)、陥凹平均深さ(No.18)、陥凹乳頭面積比(No.15)、上鼻側乳頭周囲網膜神経線維平均層厚(No.38)、年齢(No.2)、及び等価球面度数(No.3)が選択された。
サポートベクターマシンでは、最適なパラメータ種別として、乳頭水平径(No.12)、陥凹平均深さ(No.18)、鼻側リム乳頭面積比(No.22)、水平乳頭傾斜度(No.26)、上鼻側乳頭周囲網膜神経線維平均層厚(No.38)、上方乳頭周囲網膜神経線維平均層厚(No.32)、及び、乳頭網膜色素上皮高低差(No.28)、年齢(No.2)、及び等価球面度数(No.3)が選択された。
ナイーブベイズでは、最適なパラメータ種別として、上方乳頭周囲網膜神経線維平均層厚(No.32)、鼻側リム乳頭面積比(No.22)、水平乳頭傾斜度(No.26)、陥凹平均深さ(No.18)、年齢(No.2)、及び等価球面度数(No.3)が選択された。
すなわち、ニューラルネットワーク、勾配ブースティング決定木、サポートベクターマシン、及びナイーブベイズのいずれの場合においても、OCTに基づく幾つかの眼底形態データと、2つの背景データ(年齢、等価球面度数)とが、最適なパラメータ種別として選択された。このような研究結果を利用して眼科情報処理システム1が構築される。
〈データ受付部40〉
データ受付部40は、被検者の眼底形態データと背景データとを受け付ける。本実施形態では、例えば、眼底形態データは、OCTシステム200により生成されてデータベース300(電子カルテシステム、PACS等)に保存され、且つ、背景データは、カルテ入力装置250を用いて生成されてデータベース300(電子カルテシステム)に保存される。
データ受付部40は、眼底形態データ及び背景データをデータベース300から受け付ける。例えば、医師等が1以上の患者を指定したとき、指定された各患者のデータが、データ受付部40を介して眼科情報処理システム1に入力される。或いは、患者ID等が眼科情報処理システム1に入力されたとき、制御部10は、その患者ID又はそれに関連付けられた識別情報を検索クエリとして、データ受付部40を介してデータベース300からデータを取得する。
〈通信部50〉
通信部50は、他のコンピュータ(外部コンピュータ400)にデータを送信する処理と、外部コンピュータ400からデータを受信する処理とを行う。通信部50は、外部コンピュータ400との間の通信方式に応じた公知の通信機器を含む。なお、データ受付部40と通信部50とが一体であってもよい。
〈データ処理部30〉
データ処理部30は、各種のデータ処理を実行する。データ処理部30は、分類器31と、学習器32とを含む。
〈分類器31〉
分類器31は、データ受付部40により受け付けられた眼底形態データと背景データとに基づいて、視神経乳頭形状分類を実行するプロセッサの例である。本実施形態では、分類器31は、人工知能エンジンを含み、記憶部20に記憶されている学習情報21と、データ受付部40により受け付けられた眼底形態データ及び背景データとに基づいて、視神経乳頭形状分類を実行する。
ニューラルネットワークを利用した人工知能エンジンが適用される場合、例えば、眼底形態データとして、陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、及び、上方乳頭周囲網膜神経線維平均層厚が使用され、且つ、背景データとして年齢及び等価球面度数が使用される。また、学習情報21は、これら眼底形態データと同種のデータ群とこれら背景データと同種のデータ群とに基づく学習により生成された分類モデルを予め含んでいる。この分類モデルの生成は、ニューラルネットワークを利用した学習器32により実行される。ニューラルネットワークを利用した分類器31(ニューラルネットワーク分類器)は、この分類モデルと、被検者の眼底形態データ及び背景データとに基づいて、ニコレラ分類の4つの類型のそれぞれについて確信度を求める。なお、4つの類型の全てについて確信度を求める代わりに、いずれか1つ、2つ、又は3つの類型についてのみ確信度を求めるようにしてもよい。
勾配ブースティング決定木を利用した人工知能エンジンが適用される場合、例えば、眼底形態データとして、水平乳頭傾斜度、鼻側リム乳頭面積比、陥凹平均深さ、陥凹乳頭面積比、及び、上鼻側乳頭周囲網膜神経線維平均層厚が使用され、且つ、背景データとして年齢及び等価球面度数が使用される。また、学習情報21は、これら眼底形態データと同種のデータ群とこれら背景データと同種のデータ群とに基づく学習により生成された分類モデルを予め含んでいる。この分類モデルの生成は、勾配ブースティング決定木を利用した学習器32により実行される。勾配ブースティング決定木を利用した分類器31(勾配ブースティング決定木分類器)は、この分類モデルと、被検者の眼底形態データ及び背景データとに基づいて、ニコレラ分類の4つの類型のそれぞれについて確信度を求める。なお、いずれか1つ、2つ、又は3つの類型についてのみ確信度を求めるようにしてもよい。
サポートベクターマシンを利用した人工知能エンジンが適用される場合、例えば、眼底形態データとして、乳頭水平径、陥凹平均深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、及び、乳頭網膜色素上皮高低差が使用され、且つ、背景データとして年齢及び等価球面度数が使用される。また、学習情報21は、これら眼底形態データと同種のデータ群とこれら背景データと同種のデータ群とに基づく学習により生成された分類モデルを予め含んでいる。この分類モデルの生成は、サポートベクターマシンを利用した学習器32により実行される。サポートベクターマシンを利用した分類器31(サポートベクターマシン分類器)は、この分類モデルと、被検者の眼底形態データ及び背景データとに基づいて、ニコレラ分類の4つの類型のそれぞれについて確信度を求める。なお、いずれか1つ、2つ、又は3つの類型についてのみ確信度を求めるようにしてもよい。
ナイーブベイズを利用した人工知能エンジンが適用される場合、例えば、眼底形態データとして、上方乳頭周囲網膜神経線維平均層厚、鼻側リム乳頭面積比、水平乳頭傾斜度、及び、陥凹平均深さが使用され、且つ、背景データとして年齢及び等価球面度数が使用される。また、学習情報21は、これら眼底形態データと同種のデータ群とこれら背景データと同種のデータ群とに基づく学習により生成された分類モデルを予め含んでいる。この分類モデルの生成は、ナイーブベイズを利用した学習器32により実行される。ナイーブベイズを利用した分類器31(ナイーブベイズ分類器)は、この分類モデルと、被検者の眼底形態データ及び背景データとに基づいて、ニコレラ分類の4つの類型のそれぞれについて確信度を求める。なお、いずれか1つ、2つ、又は3つの類型についてのみ確信度を求めるようにしてもよい。
〈学習器32〉
学習器32は、眼底形態データと同種のデータ群と背景データと同種のデータ群とに基づき機械学習を実行することにより学習情報21を生成する。
ニューラルネットワークを利用した人工知能エンジンが適用される場合、学習器32は、例えば、陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、年齢、及び等価球面度数に基づき、ニューラルネットワークを利用した機械学習を行うことにより、学習情報21(分類モデル)を生成する。ニューラルネットワーク分類器31は、この分類モデルに基づきニコレラ分類を行う。
勾配ブースティング決定木を利用した人工知能エンジンが適用される場合、学習器32は、例えば、水平乳頭傾斜度、鼻側リム乳頭面積比、陥凹平均深さ、陥凹乳頭面積比、上鼻側乳頭周囲網膜神経線維平均層厚、年齢、及び等価球面度数に基づき、勾配ブースティング決定木を利用した機械学習を行うことにより、学習情報21(分類モデル)を生成する。勾配ブースティング決定木分類器31は、この分類モデルに基づきニコレラ分類を行う。
サポートベクターマシンを利用した人工知能エンジンが適用される場合、学習器32は、例えば、乳頭水平径、陥凹平均深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、乳頭網膜色素上皮高低差、年齢、及び等価球面度数に基づき、サポートベクターマシンを利用した機械学習を行うことにより、学習情報21(分類モデル)を生成する。サポートベクターマシン分類器31は、この分類モデルに基づきニコレラ分類を行う。
ナイーブベイズを利用した人工知能エンジンが適用される場合、学習器32は、例えば、上方乳頭周囲網膜神経線維平均層厚、鼻側リム乳頭面積比、水平乳頭傾斜度、陥凹平均深さ、年齢、及び等価球面度数に基づき、ナイーブベイズを利用した機械学習を行うことにより、学習情報21(分類モデル)を生成する。ナイーブベイズ分類器31は、この分類モデルに基づきニコレラ分類を行う。
〈ユーザインターフェイス100〉
ユーザインターフェイス100は、表示部101と操作部102とを含む。表示部101は、フラットパネルディスプレイ等の表示デバイスを含む。操作部102は、マウス、キーボード、トラックパッド、ボタン、キー、ジョイスティック、操作パネル等の操作デバイスを含む。
表示部101と操作部102は、それぞれ個別のデバイスとして構成される必要はない。たとえばタッチパネルのように、表示機能と操作機能とが一体化されたデバイスを用いることも可能である。その場合、操作部102は、このタッチパネルとコンピュータプログラムとを含んで構成される。操作部102に対する操作内容は、電気信号として制御部10に入力される。また、表示部101に表示されたグラフィカルユーザインターフェイス(GUI)と、操作部102とを用いて、操作や情報入力を行うようにしてもよい。
〈使用形態〉
眼科情報処理システム1によって実施可能な使用形態について説明する。使用形態の例を図3に示す。この例示的な使用形態では、学習情報21の生成に続いて、特定の被検者に関する視神経乳頭形状分類(ニコレラ分類)を実行しているが、これらを連続的に行う必要はない。例えば、学習情報21を所定期間毎に更新するように構成することが可能である。
(S1:サンプルデータセットを受け付ける)
まず、眼科情報処理システム1のデータ受付部40が、データベース300に保存されているサンプルデータセット(眼底形態データ群及び患者背景データ)を、ネットワークを通じて受け付ける。
例えば、制御部10は、サンプルデータセットの送信要求や、送信対象のデータの種別を示す情報などを、データベース300に送信するように、データ受付部40を制御する。具体例として、ニューラルネットワークが採用されている場合、制御部10は、サンプルデータセットの送信要求(対象疾患名を含む)と、ニューラルネットワークに対応するデータ種別(陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、年齢、及び等価球面度数)を示す情報とを、データベース300に送信するように、データ受付部40を制御する。
データベース300は、例えば、電子カルテ等を参照することで対象疾患名(例えば、開放隅角緑内障)に該当する被検者(患者)を特定し、特定された被検者の電子カルテから年齢及び等価球面度数を取得し、且つ、特定された被検者の陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、及び上方乳頭周囲網膜神経線維平均層厚を取得する。
なお、対象となるデータがデータベース300に保存されていない場合、そのデータを生成することができる。この処理は、眼科情報処理システム1以外の装置が行ってもよいし、眼科情報処理システム1が行ってもよい。
例えば、眼底形態データのいずれかが存在しない場合、データベース300(PACS等)に保存されている当該被検者のOCTデータを、解析要求(解析種別情報を含む)とともにOCTシステム200に送る。OCTシステム200は、解析要求に対応する解析処理をこのOCTデータに適用することで目的のデータを生成し、データベース300に送る。データベース300は、このデータを眼底形態データとして眼科情報処理システム1に送ることができる。他の例として、電子カルテに等価球面度数が記録されていない場合、電子カルテに記録されている球面度数と乱視度数から等価球面度数を算出することができる。
対象データを眼科情報処理システム1が生成する場合、データベース300は、その対象データを求めるために必要なデータ(予め決定される)を眼科情報処理システム1に送る。データ処理部30は、データベース300から入力されたデータに基づいて対象データを生成する。
(S2:機械学習を実行する)
学習器32(又は、眼科情報処理システム1の外部に設けられた学習器)が、ステップS1で受け付けられたサンプルデータセットに基づき機械学習を実行することにより、学習情報21(分類モデル)を生成する。
(S3:学習情報を保存する)
制御部10は、ステップS2で生成された学習情報21を記憶部20に保存する。なお、過去に生成された学習情報が記憶部20に保存されている場合、過去の学習情報を新たな学習情報で更新することができる。或いは、生成日時情報とともに学習情報を蓄積するように構成してもよい。
(S4:被検者を指定する)
医師等のユーザ、眼科情報処理システム1、又は外部コンピュータ400が、1以上の被検者を指定する。
ユーザが被検者を指定する場合、眼科情報処理システム1又は外部コンピュータ400は、被検者指定用の画面を表示デバイスに表示する。ユーザは、操作デバイスを用いて1以上の被検者を指定する。このとき、例えば、検診を受けたグループをまとめて指定できるようにしてもよい。また、特定の疾患の疑いが有る被検者を(まとめて)指定できるようにしてもよい。
眼科情報処理システム1又は外部コンピュータ400が被検者を指定する場合、日時、特定疾患名、特定疾患の疑いの有無、スクリーニングの種別などに基づいて、被検者を指定することが可能である。
(S5:眼底形態データと背景データを受け付ける)
ステップS4で指定された1以上の被検者のそれぞれについて、眼科情報処理システム1のデータ受付部40は、データベース300に保存されている眼底形態データ群及び患者背景データを、ネットワークを通じて受け付ける。
ここで、送信要求や、送信対象のデータの種別を示す情報や、データベース300の動作や、眼科情報処理システム1の動作については、ステップS1と同様であってよい。
(S6:各類型の確信度を求める)
分類器31は、ステップS5で受け付けられた眼底形態データ群及び患者背景データと、ステップS3で記憶部20に保存された学習情報21とに基づいて、各被検者の各被検眼に関するニコレラ分類を実行する。それにより、例えば、ニコレラ分類の4つの類型のそれぞれについて確信度を求める。
(S7:分類結果を表示する)
出力制御部11は、ステップS6で求められた各類型の確信度を表示部101に表示する。
確信度の表示態様の例を図4に示す。本例の表示態様では、ニコレラ分類の4類型の確信度(Confidence)が、その値の大きい順に上から列挙される。図4に示す例は、上から順に、加齢性硬化型(SS)の確信度「Confidence(SS)」、局所虚血型(FI)の確信度「Confidence(FI)」、近視型(MY)の確信度「Confidence(MY)」、全体的拡大型(GE)の確信度「Confidence(GE)」が配置されたテーブル情報を含む。更に、確信度の数値とともに、その大きさを示す棒グラフ(点線で示す)が提示されている。
このような分類結果の表示態様によれば、ニコレラ分類の各類型の確信度を把握することができ、更に、複数の類型の大きさの相対的な差を容易に把握することができる。よって、被検眼の視神経乳頭の形状が、ニコレラ分類のどの類型に属するか容易に知ることができる。また、眼科情報処理システム1により行われた分類がどの程度の信頼性を有するか容易に知ることができる。
(S8:分類結果を保存する)
制御部10は、ステップS6で求められた各類型の確信度や、図7で表示された情報を、記憶部20に保存することができる。また、制御部10は、このような分類結果を、データベース300や外部コンピュータ400に送信することができる。データベース300や外部コンピュータ400は、眼科情報処理システム1から受信した分類結果を保存することができる。
〈作用・効果〉
以上に説明した例示的な実施形態の作用及び効果について説明する。
例示的な実施形態に係る眼科情報処理システムは、受付部と、プロセッサとを含む。受付部は、被検者の眼底の形態データと当該被検者の背景データとを受け付ける。プロセッサは、受付部により受け付けられた形態データ及び背景データに基づいて、視神経乳頭形状分類を実行する。
上記の眼科情報処理システム1では、データ受付部40が受付部に相当し、分類器31がプロセッサに相当する。また、眼科情報処理システム1は、視神経乳頭形状分類としてニコレラ分類を実施する。
例示的な実施形態に係る眼科情報処理システムは、記憶部を更に含んでいてよい。本例では、視神経乳頭形状分類に用いられる形態データと同種のデータ群及び背景データと同種のデータ群に基づく学習によって、学習情報が生成される。記憶部は、この学習情報を予め記憶する。プロセッサは、記憶部に記憶されている学習情報と、受付部により受け付けられた形態データ及び背景データに基づいて、視神経乳頭形状分類を実行する。
上記の眼科情報処理システム1では、記憶部20が記憶部に相当し、学習情報21が学習情報に相当する。
例示的な実施形態に係る眼科情報処理システムにおいて、形態データは、陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、及び、上方乳頭周囲網膜神経線維平均層厚のうちの少なくとも1つを含んでよい。また、背景データは、年齢及び等価球面度数の少なくとも一方を含んでよい。更に、学習情報は、ニューラルネットワークを用いた学習によって生成された分類モデルを含む。加えて、プロセッサは、ニューラルネットワーク分類器を含む。
上記の実施形態においてニューラルネットワークが適用される場合、眼科情報処理システム1における典型的な形態データは、陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、及び、上方乳頭周囲網膜神経線維平均層厚の全てを含む。更に、典型的な背景データは、年齢及び等価球面度数の双方を含む。また、分類器31がニューラルネットワーク分類器に相当する。
例示的な実施形態に係る眼科情報処理システムにおいて、形態データは、水平乳頭傾斜度、鼻側リム乳頭面積比、陥凹平均深さ、陥凹乳頭面積比、及び、上鼻側乳頭周囲網膜神経線維平均層厚のうちの少なくとも1つを含む。また、背景データは、年齢及び等価球面度数の少なくとも一方を含む。更に、学習情報は、勾配ブースティング決定木を用いた学習によって生成された分類モデルを含む。加えて、プロセッサは、勾配ブースティング決定木分類器を含む。
上記の実施形態において勾配ブースティング決定木が適用される場合、眼科情報処理システム1における典型的な形態データは、水平乳頭傾斜度、鼻側リム乳頭面積比、陥凹平均深さ、陥凹乳頭面積比、及び、上鼻側乳頭周囲網膜神経線維平均層厚の全てを含む。更に、典型的な背景データは、年齢及び等価球面度数の双方を含む。また、分類器31が勾配ブースティング決定木分類器に相当する。
例示的な実施形態に係る眼科情報処理システムにおいて、形態データは、乳頭水平径、陥凹平均深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、及び、乳頭網膜色素上皮高低差のうちの少なくとも1つを含む。また、背景データは、年齢及び等価球面度数の少なくとも一方を含む。更に、学習情報は、サポートベクターマシンを用いた学習によって生成された分類モデルを含む。加えて、プロセッサは、サポートベクターマシン分類器を含む。
上記の実施形態においてサポートベクターマシンが適用される場合、眼科情報処理システム1における典型的な形態データは、乳頭水平径、陥凹平均深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、及び、乳頭網膜色素上皮高低差の全てを含む。更に、典型的な背景データは、年齢及び等価球面度数の双方を含む。また、分類器31がサポートベクターマシン分類器に相当する。
例示的な実施形態に係る眼科情報処理システムにおいて、形態データは、上方乳頭周囲網膜神経線維平均層厚、鼻側リム乳頭面積比、水平乳頭傾斜度、及び、陥凹平均深さのうちの少なくとも1つを含む。また、背景データは、年齢及び等価球面度数の少なくとも一方を含む。更に、学習情報は、ベイズ(の定理)(例えばナイーブベイズ)を用いた学習によって生成された分類モデルを含む。加えて、プロセッサは、ベイズ分類器を含む。
上記の実施形態においてベイズが適用される場合、眼科情報処理システム1における典型的な形態データは、上方乳頭周囲網膜神経線維平均層厚、鼻側リム乳頭面積比、水平乳頭傾斜度、及び、陥凹平均深さの全てを含む。更に、典型的な背景データは、年齢及び等価球面度数の双方を含む。また、分類器31がベイズ分類器に相当する。
例示的な実施形態に係る眼科情報処理システムは、学習器を更に含んでよい。学習器は、視神経乳頭形状分類に用いられる形態データと同種のデータ群及び背景データと同種のデータ群に基づく学習を実行することで学習情報を生成する。
上記の眼科情報処理システム1では、学習器32が学習器に相当する。
例示的な実施形態に係る眼科情報処理システムにおいて、プロセッサは、視神経乳頭形状分類により、視神経乳頭形状の1以上の類型のそれぞれについての確信度を求める。更に、例示的な実施形態に係る眼科情報処理システムは、視神経乳頭形状の1以上の類型の確信度を表示デバイスに表示させる表示制御部を更に含んでよい。
上記の眼科情報処理システム1の典型的な使用形態では、ニコレラ分類の4類型の全てについて確信度を求めている。更に、上記の眼科情報処理システム1では、出力制御部11が表示制御部に相当し、表示部101が表示デバイスに相当する。なお、表示デバイスは眼科情報処理システムに含まれてもよいし、含まれなくてもよい。
例示的な実施形態に係る眼科情報処理システムにおいて、視神経乳頭形状分類に用いられる形態データの少なくとも一部は、光コヒーレンストモグラフィを使用して収集された眼底のデータを解析することにより生成されてよい。
上記の眼科情報処理システム1では、眼底の3次元領域がOCTスキャンされる。解析では、このOCTスキャンにより収集された3次元OCTデータにセグメンテーションを適用することで、眼底の層構造が特定される。この層構造に基づいて、乳頭領域(disc)、陥凹領域(cup)、リム領域(rim)、乳頭周囲網膜神経線維層などが特定される。そして、乳頭領域の面積や径、陥凹領域の面積や径、リム領域の面積や径などが算出される。また、乳頭領域、陥凹領域、及びリム領域のうちのいずれか2つの面積比や径比が算出される。更に、陥凹領域の深さ、乳頭領域の傾きなどが算出される。また、複数のセクターにおける乳頭周囲網膜神経線維層厚が算出され、更に、これらセクターにおける層厚値に基づき特徴量が算出される。
例示的な実施形態に係る眼科情報処理システムは、解析プロセッサを更に含んでよい。解析プロセッサは、光コヒーレンストモグラフィを使用して収集された眼底のデータを解析して形態データを生成する。
上記の眼科情報処理システム1では、OCTシステム200(解析機能)、データ処理部30、及び外部コンピュータ400のうちの少なくとも1つが、解析プロセッサとして機能する。
例示的な実施形態に係る眼科情報処理システムは、光コヒーレンストモグラフィを使用して眼底のデータを収集する光コヒーレンストモグラフィ装置を更に含んでよい。
上記の眼科情報処理システム1では、OCTシステム200(データ収集機能)が光コヒーレンストモグラフィ装置として機能する。
例示的な実施形態に係る眼科情報処理システムにおいて、受付部は、背景データの少なくとも一部を電子カルテシステムから直接的に又は間接的に受け付けるように構成されてよい。
上記の眼科情報処理システム1では、データベース300に電子カルテシステムが含まれ、年齢及び等価球面度数が電子カルテシステムから取得される。
例示的な実施形態に係る眼科情報処理方法は、コンピュータを用いて実現される。コンピュータは、被検者の眼底の形態データと当該被検者の背景データとを受け付け、受け付けられた形態データと背景データとに基づいて視神経乳頭形状分類を実行する。
例示的な実施形態に係るプログラムは、このような例示的な実施形態に係る眼科情報処理方法をコンピュータに実行させるプログラムである。
例示的な実施形態に係る記録媒体は、このような例示的な実施形態に係るプログラムを記録したコンピュータ読み取り可能な記録媒体である。
以上のように構成された例示的な実施形態によれば、医師の主観に頼る従来の乳頭形状分類に代わる、自動的な乳頭形状分類が可能となる。それにより、例えば、有病率が比較的高い疾患や早期発見が重要な疾患のスクリーニングを広く提供することが可能になる。また、緑内障などの疾患の診療において、乳頭形状に応じた治療方針の決定を容易に行うことができる。
更に、例示的な実施形態によれば、光コヒーレンストモグラフィを用いて収集されたデータを解析して乳頭形状パラメータを定量化することが可能である。したがって、眼底カメラやSLOで得られた眼底像を参照する従来のケースと比較して、乳頭形状パラメータの定量化が容易であり、乳頭形状分類の自動化がより促進される。
また、例示的な実施形態によれば、乳頭形状分類に対する貢献度が高いパラメータを選択的に参照することができるので、無駄に多くのデータを参照することなく良好な精度及び確度で分類を行うことが可能である。
このように、例示的な実施形態によれば、良好な精度及び確度での乳頭形状分類を自動で行うことが可能である。
また、前述した特徴選択を実行可能に構成された例示的な実施形態によれば、大規模なシステムを用いることなく、良好な精度及び確度での乳頭形状分類を自動で行うことが可能である。
実施形態に係る眼科情報処理方法に含まれる工程や、眼科情報処理システムに含まれる要素(構成、動作等)は、上記の例示に限定されるものではない。
以上に説明した実施形態は本発明の一例に過ぎない。本発明を実施しようとする者は、本発明の要旨の範囲内における変形(省略、置換、付加等)を任意に施すことが可能である。
例示的な変形として、前述したような特徴選択を実行しない構成を適用することができる。例えば、ランダムフォレスト等の学習アルゴリズムを用いることにより、特徴選択を行うことなく視神経乳頭形状分類を行うように構成することが可能である。
他の例示的な変形として、眼底カメラやSLOを用いて取得された眼底像を利用することも可能である。例えば、眼底カメラ及び/又はSLOを使用して取得された眼底像を解析することにより形態データの少なくとも一部を生成することが可能である。眼底カメラ又はSLOを用いて取得された眼底像から生成可能な形態データの例として、乳頭面積、陥凹面積、リム面積、乳頭径、陥凹径、リム径などがある。
1 眼科情報処理システム
10 制御部
11 出力制御部
20 記憶部
21 学習情報
30 データ処理部
31 分類器
32 学習器
40 データ受付部
50 通信部
100 ユーザインターフェイス
101 表示部
102 操作部
200 OCTシステム
250 カルテ入力装置
300 データベース
400 外部コンピュータ

Claims (14)

  1. 被検者の視神経乳頭の形態データと当該被検者の背景データとを受け付ける受付部と、
    前記形態データと同種のデータ群及び前記背景データと同種のデータ群に基づき機械学習を実行することにより生成された学習情報を予め記憶する記憶部と、
    前記学習情報と前記形態データと前記背景データとに基づいて、前記機械学習に用いられた人工知能エンジンと同種の人工知能エンジンを適用することにより、視神経乳頭形状分類を実行するプロセッサと
    を含む眼科情報処理システム。
  2. 前記形態データは、陥凹面積、陥凹平均深さ、陥凹最大深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、及び、上方乳頭周囲網膜神経線維平均層厚のうちの少なくとも1つを含み、
    前記背景データは、年齢及び等価球面度数の少なくとも一方を含み、
    前記学習情報は、ニューラルネットワークを用いた学習によって生成された分類モデルを含み、
    前記プロセッサは、ニューラルネットワーク分類器を含む
    ことを特徴とする請求項に記載の眼科情報処理システム。
  3. 前記形態データは、水平乳頭傾斜度、鼻側リム乳頭面積比、陥凹平均深さ、陥凹乳頭面積比、及び、上鼻側乳頭周囲網膜神経線維平均層厚のうちの少なくとも1つを含み、
    前記背景データは、年齢及び等価球面度数の少なくとも一方を含み、
    前記学習情報は、勾配ブースティング決定木を用いた学習によって生成された分類モデルを含み、
    前記プロセッサは、勾配ブースティング決定木分類器を含む
    ことを特徴とする請求項に記載の眼科情報処理システム。
  4. 前記形態データは、乳頭水平径、陥凹平均深さ、鼻側リム乳頭面積比、水平乳頭傾斜度、上鼻側乳頭周囲網膜神経線維平均層厚、上方乳頭周囲網膜神経線維平均層厚、及び、乳頭網膜色素上皮高低差のうちの少なくとも1つを含み、
    前記背景データは、年齢及び等価球面度数の少なくとも一方を含み、
    前記学習情報は、サポートベクターマシンを用いた学習によって生成された分類モデルを含み、
    前記プロセッサは、サポートベクターマシン分類器を含む
    ことを特徴とする請求項に記載の眼科情報処理システム。
  5. 前記形態データは、上方乳頭周囲網膜神経線維平均層厚、鼻側リム乳頭面積比、水平乳頭傾斜度、及び、陥凹平均深さのうちの少なくとも1つを含み、
    前記背景データは、年齢及び等価球面度数の少なくとも一方を含み、
    前記学習情報は、ベイズを用いた学習によって生成された分類モデルを含み、
    前記プロセッサは、ベイズ分類器を含む
    ことを特徴とする請求項に記載の眼科情報処理システム。
  6. 前記形態データと同種のデータ群及び前記背景データと同種のデータ群に基づき前記機械学習を実行することにより前記学習情報を生成する学習器を更に含む
    ことを特徴とする請求項1~5のいずれかに記載の眼科情報処理システム。
  7. 前記プロセッサは、前記視神経乳頭形状分類により、視神経乳頭形状の1以上の類型のそれぞれについての確信度を求め、
    前記1以上の類型の前記確信度を表示デバイスに表示させる表示制御部を更に含む
    ことを特徴とする請求項1~のいずれかに記載の眼科情報処理システム。
  8. 前記形態データの少なくとも一部は、光コヒーレンストモグラフィを使用して収集された眼底のデータを解析することにより生成される
    ことを特徴とする請求項1~のいずれかに記載の眼科情報処理システム。
  9. 光コヒーレンストモグラフィを使用して収集された眼底のデータを解析して形態データを生成する解析プロセッサを更に含む
    ことを特徴とする請求項に記載の眼科情報処理システム。
  10. 光コヒーレンストモグラフィを使用して眼底のデータを収集する光コヒーレンストモグラフィ装置を更に含む
    ことを特徴とする請求項に記載の眼科情報処理システム。
  11. 前記受付部は、前記背景データの少なくとも一部を電子カルテシステムから直接的に又は間接的に受け付ける
    ことを特徴とする請求項1~10のいずれかに記載の眼科情報処理システム。
  12. コンピュータを用いた眼科情報処理方法であって、
    前記コンピュータが、
    被検者の視神経乳頭の形態データと当該被検者の背景データとを受け付け、
    前記形態データと同種のデータ群及び前記背景データと同種のデータ群に基づき機械学習を実行することにより生成された学習情報を取得し、
    前記学習情報と前記形態データと前記背景データとに基づいて、前記機械学習に用いられた人工知能エンジンと同種の人工知能エンジンを適用することにより、視神経乳頭形状分類を実行する
    ことを特徴とする眼科情報処理方法。
  13. 請求項12に記載の眼科情報処理方法をコンピュータに実行させるプログラム。
  14. 請求項13に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
JP2017124825A 2017-06-27 2017-06-27 眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体 Active JP7078948B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2017124825A JP7078948B2 (ja) 2017-06-27 2017-06-27 眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体
US15/925,922 US10765316B2 (en) 2017-06-27 2018-03-20 Ophthalmic information processing system, ophthalmic information processing method, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017124825A JP7078948B2 (ja) 2017-06-27 2017-06-27 眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体

Publications (2)

Publication Number Publication Date
JP2019005319A JP2019005319A (ja) 2019-01-17
JP7078948B2 true JP7078948B2 (ja) 2022-06-01

Family

ID=64691631

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017124825A Active JP7078948B2 (ja) 2017-06-27 2017-06-27 眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体

Country Status (2)

Country Link
US (1) US10765316B2 (ja)
JP (1) JP7078948B2 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101977645B1 (ko) * 2017-08-25 2019-06-12 주식회사 메디웨일 안구영상 분석방법
CN111757705B (zh) * 2018-03-05 2023-07-28 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法
JP2020103579A (ja) * 2018-12-27 2020-07-09 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
US11191492B2 (en) * 2019-01-18 2021-12-07 International Business Machines Corporation Early detection and management of eye diseases by forecasting changes in retinal structures and visual function
WO2020202680A1 (ja) * 2019-03-29 2020-10-08 キヤノン株式会社 情報処理装置及び情報処理方法
KR102282334B1 (ko) * 2019-08-08 2021-07-28 인천대학교 산학협력단 시신경 유두 분류 방법
JP2021145805A (ja) * 2020-03-18 2021-09-27 株式会社トプコン 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム
JP7483247B2 (ja) * 2020-04-13 2024-05-15 株式会社トーメーコーポレーション 検査対象分類システム、検査対象分類方法および検査対象分類プログラム
KR102482680B1 (ko) * 2020-06-03 2022-12-28 고려대학교 산학협력단 안저영상을 기초로 생체 계측치를 예측하는 장치 및 방법
KR102596340B1 (ko) * 2020-08-25 2023-11-02 사회복지법인 삼성생명공익재단 안질환 시야 손상 예후 예측 시스템 및 방법과 이를 위한 컴퓨터 프로그램
EP4241652A4 (en) 2020-11-04 2024-10-16 Topcon Corp OPHTHALMIC INFORMATION PROCESSING DEVICE, OPHTHALMIC DEVICE, OPHTHALMIC INFORMATION PROCESSING METHOD, AND PROGRAM
JP2022104239A (ja) 2020-12-28 2022-07-08 株式会社トプコン 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム
WO2023017875A1 (ko) * 2021-08-12 2023-02-16 고려대학교 산학협력단 안저영상을 기초로 생체 계측치를 예측하는 장치 및 방법
KR102704938B1 (ko) * 2021-12-29 2024-09-09 단국대학교 산학협력단 녹내장 진단 장치, 시스템 및 방법
WO2023195554A1 (ko) * 2022-04-06 2023-10-12 경상국립대학교병원 녹내장 위험도 진단 장치 및 방법

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100277691A1 (en) 2009-04-30 2010-11-04 University Of Southern California Methods for Diagnosing Glaucoma Utilizing Combinations of FD-OCT Measurements from Three Anatomical Regions of the Eye
JP2010282366A (ja) 2009-06-03 2010-12-16 Canon Inc 診断支援システム及びその診断支援方法、情報処理装置
JP2016002380A (ja) 2014-06-18 2016-01-12 キヤノン株式会社 画像処理装置、その作動方法及びプログラム
US20170112372A1 (en) 2015-10-23 2017-04-27 International Business Machines Corporation Automatically detecting eye type in retinal fundus images

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3585331B2 (ja) 1996-12-03 2004-11-04 株式会社ニデック 眼底立体画像の解析方法
JP2000245700A (ja) 1999-03-01 2000-09-12 Nidek Co Ltd 眼底計測装置及び眼底計測プログラムを記録した記録媒体
US7500751B2 (en) 2004-03-12 2009-03-10 Yokohama Tlo Company Ltd. Ocular fundus portion analyzer and ocular fundus portion analyzing method
JP2006263127A (ja) 2005-03-24 2006-10-05 Gifu Univ 眼底画像診断支援システム、及び眼底画像診断支援プログラム
JP5278984B2 (ja) 2006-03-24 2013-09-04 興和株式会社 画像解析装置及び画像解析プログラム
JP4915737B2 (ja) 2007-03-13 2012-04-11 興和株式会社 画像解析システム、及び画像解析プログラム
JP2008237840A (ja) 2007-03-29 2008-10-09 Gifu Univ 画像解析システム、及び画像解析プログラム
JP2008237839A (ja) 2007-03-29 2008-10-09 Gifu Univ 画像解析システム、及び画像解析プログラム
US9700210B2 (en) * 2009-03-02 2017-07-11 Canon Kabushiki Kaisha Image processing apparatus and method for controlling the same
JP5701494B2 (ja) 2009-10-19 2015-04-15 株式会社ニデック 眼底画像処理装置、及び眼底画像処理方法
US10149610B2 (en) * 2014-04-25 2018-12-11 Carl Zeiss Meditec, Inc. Methods and systems for automatic detection and classification of ocular inflammation
JP6568375B2 (ja) 2015-03-30 2019-08-28 キヤノン株式会社 眼科情報処理システム、画像処理装置、および画像処理方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100277691A1 (en) 2009-04-30 2010-11-04 University Of Southern California Methods for Diagnosing Glaucoma Utilizing Combinations of FD-OCT Measurements from Three Anatomical Regions of the Eye
JP2010282366A (ja) 2009-06-03 2010-12-16 Canon Inc 診断支援システム及びその診断支援方法、情報処理装置
JP2016002380A (ja) 2014-06-18 2016-01-12 キヤノン株式会社 画像処理装置、その作動方法及びプログラム
US20170112372A1 (en) 2015-10-23 2017-04-27 International Business Machines Corporation Automatically detecting eye type in retinal fundus images

Also Published As

Publication number Publication date
US10765316B2 (en) 2020-09-08
JP2019005319A (ja) 2019-01-17
US20180368679A1 (en) 2018-12-27

Similar Documents

Publication Publication Date Title
JP7078948B2 (ja) 眼科情報処理システム、眼科情報処理方法、プログラム、及び記録媒体
Mursch-Edlmayr et al. Artificial intelligence algorithms to diagnose glaucoma and detect glaucoma progression: translation to clinical practice
Sarhan et al. Machine learning techniques for ophthalmic data processing: a review
Zhang et al. A survey on computer aided diagnosis for ocular diseases
Vessani et al. Comparison of quantitative imaging devices and subjective optic nerve head assessment by general ophthalmologists to differentiate normal from glaucomatous eyes
US9357911B2 (en) Integration and fusion of data from diagnostic measurements for glaucoma detection and progression analysis
Zhu et al. Predicting visual function from the measurements of retinal nerve fiber layer structure
Huang et al. Diagnostic power of optic disc morphology, peripapillary retinal nerve fiber layer thickness, and macular inner retinal layer thickness in glaucoma diagnosis with fourier-domain optical coherence tomography
US20220400943A1 (en) Machine learning methods for creating structure-derived visual field priors
Bowd et al. Relevance vector machine and support vector machine classifier analysis of scanning laser polarimetry retinal nerve fiber layer measurements
Russakoff et al. A 3D deep learning system for detecting referable glaucoma using full OCT macular cube scans
US20110190657A1 (en) Glaucoma combinatorial analysis
Kauppi Eye fundus image analysis for automatic detection of diabetic retinopathy
Zhang et al. The application of artificial intelligence in glaucoma diagnosis and prediction
US20230140881A1 (en) Oct en face pathology segmentation using channel-coded slabs
Panda et al. Describing the structural phenotype of the glaucomatous optic nerve head using artificial intelligence
Saleh et al. The role of medical image modalities and AI in the early detection, diagnosis and grading of retinal diseases: a survey
WO2016185737A1 (ja) 眼科検査支援システム
Kapoor et al. The role of artificial intelligence in the diagnosis and management of glaucoma
Chen et al. Applications of artificial intelligence and deep learning in glaucoma
JP6499918B2 (ja) 眼科検査支援システム及び眼科検査支援サーバ
JP7457315B2 (ja) 情報処理装置、推論モデル生成方法、情報処理方法、及びプログラム
Huang et al. Rule extraction for glaucoma detection with summary data from StratusOCT
EP4241652A1 (en) Ophthalmic information processing device, ophthalmic device, ophthalmic information processing method, and program
Pahlitzsch et al. Significance of the disc damage likelihood scale objectively measured by a non-mydriatic fundus camera in preperimetric glaucoma

Legal Events

Date Code Title Description
A80 Written request to apply exceptions to lack of novelty of invention

Free format text: JAPANESE INTERMEDIATE CODE: A80

Effective date: 20170726

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200427

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200610

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210331

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210406

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210603

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211116

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211214

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: 20220419

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220512

R150 Certificate of patent or registration of utility model

Ref document number: 7078948

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150