[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP2020518915A - 自動眼底画像分析用のシステムおよび方法 - Google Patents

自動眼底画像分析用のシステムおよび方法 Download PDF

Info

Publication number
JP2020518915A
JP2020518915A JP2019559095A JP2019559095A JP2020518915A JP 2020518915 A JP2020518915 A JP 2020518915A JP 2019559095 A JP2019559095 A JP 2019559095A JP 2019559095 A JP2019559095 A JP 2019559095A JP 2020518915 A JP2020518915 A JP 2020518915A
Authority
JP
Japan
Prior art keywords
image
neural network
regions
region
classifier
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.)
Pending
Application number
JP2019559095A
Other languages
English (en)
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.)
Retinopathy Answer Ltd
retinopathy Answer Ltd Uk
Original Assignee
Retinopathy Answer Ltd
retinopathy Answer Ltd Uk
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 Retinopathy Answer Ltd, retinopathy Answer Ltd Uk filed Critical Retinopathy Answer Ltd
Publication of JP2020518915A publication Critical patent/JP2020518915A/ja
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/18Eye characteristics, e.g. of the iris
    • G06V40/197Matching; Classification
    • 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/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
    • 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
    • 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
    • A61B3/145Arrangements specially adapted for eye photography by video means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G06F18/2148Generating training patterns; Bootstrap methods, e.g. bagging or boosting characterised by the process organisation or structure, e.g. boosting cascade
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • 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
    • G06T7/0014Biomedical image inspection using an image reference approach
    • G06T7/0016Biomedical image inspection using an image reference approach involving temporal comparison
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/46Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
    • G06V10/462Salient features, e.g. scale invariant feature transforms [SIFT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/18Eye characteristics, e.g. of the iris
    • G06V40/193Preprocessing; Feature extraction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/16Ophthalmology
    • 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/20081Training; Learning
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Ophthalmology & Optometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Multimedia (AREA)
  • Data Mining & Analysis (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computational Linguistics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Image Analysis (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

病理画像の分類を行うシステムおよび方法。画像が受信されて正規化されて、正規化された画像が複数の領域に分割される。この複数の領域に対して、ニューラルネットワークを有する少なくとも1つの分類器により、疾患ベクトルが自動で判定される。判定された疾患ベクトルに基づいて、複数の領域それぞれが自動でアノテーションされる。受信された画像は、アノテーションに少なくとも基づいて、自動で評価される。ニューラルネットワークは、少なくとも1つの客観的分類基準に応じた、画像の各領域に対する専門家によるアノテーションに少なくとも基づいて訓練される。画像は眼画像、血管画像、または眼底画像であり得る。例えば疾患とは、血管疾患、血管障害、または糖尿病網膜症であり得る。【選択図】図1

Description

本発明は、ニューラルネットワーク(「NN」)またはディープニューラルネットワーク(「DNN」)を利用した自動眼底画像分析の分野に関する。
本明細書で引用する各参照文献は、全体が参照により本明細書に組み込まれる。
一般的に入手可能な各種専門的カメラにより、人の網膜写真を容易に撮影可能である。英国において、すでにそのようなカメラは大多数の眼鏡技師(特に有名チェーン店に所属)の勤務先に存在し、網膜写真が基本的な目の検査の一部に利用されることも多くなっている。
網膜画像は、多くの眼病(いわゆる一般的な網膜症)評価に利用され得、近年のデジタルカメラの解像度向上もあって、いまや網膜の状態評価に好んで使用されるようにもなっている。画像から判定できる網膜症の種類も、出血から腫瘍までと幅広い。これら疾患の多くは失明を含む深刻な状態に繋がり得るにも関わらず、早期段階では患者が気付くような兆候が発現しないことが多い。治療成功および/または完治のためには、早期処置が重要なことが多く、より高頻度の網膜画像評価を強く後押しする公衆衛生上の根拠も数多い。
網膜上の少量出血(微小動脈瘤)を伴う二型糖尿病は特に関連深い疾患である。これは処置しなければ、患者が二型糖尿病と診断される前に、急速に深刻な眼病の発症につながる可能性もある。
糖尿病網膜症(「DR」)は、英国を含む先進国において、成人の失明の主要な原因である。したがって、英国国民保険サービスは、二型糖尿病の患者に毎年無料で網膜検査を提供している。具体的には、デジタル写真が撮影されて、評価センターに送られる。そこで、熟練の評価者(人)が画像を詳細に検査し、最速で出血の兆候を発見する。この手動の処理はコストがかかるが、現状でこの種の画像の処理に対して、より信頼性の高い手法が他にはない。長年の研究にも関わらず、訓練された人間の場合と比較して、信頼性、精度で引けを取らない自動システムが実現されていない。
U.S.20170039689は、画像からDRを評価する「ディープラーニング」技術を開示する。加えてU.S.20170039412、20150110348、20150110368、20150110370、20150110372、9,008,391、9,002,085、8,885,901、8,879,813、および20170046616も参照されたい。網膜眼底画像は、血管の状態の指標を示すものであって、眼病とより全般的な血管疾患の両方の診断に有用となり得る。
「ディープラーニング」は、複数の隠れ層を含む人工NN(「ANN」)の進化版である。具体的には、データからの抽出および予測が向上している(参照文献2)。Greenspan,Hayit,Bram van Ginneken,and Ronald M.Summers.“Guest editorial deep learning in medical imaging:Overview and future promise of an exciting new technique.”IEEE Transactions on Medical Imaging 35.5(2016):1153−1159を参照されたい。畳み込みニューラルネットワーク(「CNN」)は、コンピュータビジョンに対して極めて有用なツールである。ディープCNNは、元データ(例えば画像)からの中レベルおよび高レベル抽出物を、自動で学習するように設計可能である。CNNから抽出された一般的な記述子は、自然画像内の物体認識および位置特定に有用である。
機械学習には、教師あり学習と教師なし学習の2つの体系が存在する。教師あり学習の場合、訓練用データにはラベル付けがされる。即ち、正解は外から与えられ、統計的学習基準が外部に存在する。教師なし学習の場合、訓練用データがグラウンドトゥルース(正確さや整合性をチェックするためのデータ)を参照することなく分析される。したがって、内部に基準がある。半教師あり学習には、これら両体系の要素が盛り込まれている。即ち、訓練用データセットの全てのサンプル(またはサンプルの属性)がラベル付けされるわけではない。
医療用撮像では、疾患の正確な診断および/または評価が、撮像技術および画像判断に依存する。撮像技術は近年飛躍的に向上し、装置がデータを高速かつ高解像度で取得可能となっている。したがって、本技術は画像判断処理についてのソリューションを提供する。医療用画像は多くの場合、医師により判断される。しかし、人間の画像判断は客観性、人によって見方が大幅に異なること、そして疲労のような要素により限界がある。
人間が調べた場合、人それぞれ解釈が異なるので、過誤が生じやすい。さらに、人それぞれ基準が異なる。したがって、訓練用データを得るだけでも困難が伴う。
多くの場合、診断は異常を検出するための初検処理と、経時的に得られた測定値と変化の数量化を伴う。
CNNは、医療用画像処理に利用されている。Sahinerら(参照文献4)を参照されたい。マンモグラムから、生検で特定された物体または通常の組織を含む関心領域が抽出される。CNNは1つの入力層と、2つの隠れ層と、1つの出力層とを含み、誤差逆伝播法で訓練される。CNNは、肺結節検出(参照文献5)と、マンモグラフィ上の微小石灰化(参照文献6)にも適用されている。
画像処理用の典型的なCNN構造では、一連の畳み込みフィルタ層と、一連のデータ削減またはプーリング層とが交互に配置される。畳み込みフィルタは、入力画像の小パッチに適用される。例えば、人の脳の低レベル視覚処理のように、畳み込みフィルタはより関連性の高い画像特徴を検出する。それは例えば、直線エッジを表し得る線または円(臓器検出等)、または円(結腸ポリープのような丸い物体等に対応)である。さらに、局地的または全体的形状またはテクスチャのようなより高次の特徴にも適用される。典型的なCNN出力としては、1つ以上の確率またはクラスラベルが挙げられる。訓練用データから畳み込みフィルタが学習される。これにより、特定用途向けフィルタで画像を前処理するか、計算可能な特徴を計算しない限り必要となるような、手間のかかる手作業で特徴を生成する必要性がなくなるため、好ましい。さらに、いわゆる長/短期記憶という、深層再帰型ニューラルネットワーク(「DRNN」)等のその他ネットワーク構造変数も存在する(参照文献7)。
CNNは、並列化可能なアルゴリズムであり、画像処理装置(GPU)用コンピュータチップにより、CPU処理単独の場合と比較して、処理速度を上げることができる(約40倍)。参照文献8を参照されたい。医療用画像処理において、GPUは分割、再構成、レジストレーション、および機械学習に利用できる(参照文献9,10)。
一から深層CNN(「DCNN」)を行うこと(即ちフルトレーニング)は困難である。第1に、CNNは大量のラベル付けされた訓練用データが必要である。これは、専門家によるアノテーションが高価で、疾患(例えば病変)が数限られている医療分野では、条件として困難である。第2に、DCNNの訓練には膨大な計算およびメモリ資源が必要となる。それらが揃わないと、訓練処理にとてつもない時間がかかってしまう。第3に、DCNNの訓練は、過学習や、収束という問題で、複雑となることが多い。これら問題については、ネットワークの構造または学習パラメータについて、繰り返し調整を行い、全層が十分な速度で学習していることを保証する必要があることが多い。「転移学習」や「ファインチューニング」と呼ばれるいくつからの新しい学習手法はこれら問題を解決することが確認されており、普及してきている。しかし、現状では経験則に基づいて処理が調整されている。したがって、同じような目的に対して、同じようなアルゴリズムを使用しているのも関わらず、手法、形態、および結果に大きな差異が生じている。www.kaggle.com/c/diabetic−retinopathy−detectionを参照されたい。
なお、ネットワークの訓練はアップデートできない。したがって、アルゴリズムまたは訓練用データに何かしら変化を加える場合は、ネットワーク全体を改めて最適化する必要がある。新たな技術によりこのような制限も緩和されてきている。参照により本明細書に組み込まれるU.S.Patent No.9,053,431を参照されたい。即ち、ネットワークは構築後に、ネットワークのノイズまたはエラーベクトル出力に応じて、補強され得る。
コンピュータ支援検出(CAD)は、ディープラーニングと極めて相性が良い、一般的に利用されている医療用画像分析の一分野である。CADの標準的な手法(参照文献11)では、候補となる病変が、教師あり方法またはフィルタリングや数学的モルフォロジー手法のような伝統的な画像処理技術で検出される。多くの場合、候補となる病変は分割され、大量となりがちな手動による特徴で示される。そして分類器を利用して、候補が実際の病変である確率に対して、特徴ベクトルがマッピングされる。手動で得られた特徴ではなく、単純にディープラーニングを使う方法として、病変候補を中心とした画像データ領域で動作するCNNを訓練することができる。各候補の分類には、異なる複数のCNNの組み合わせが用いられる。
分類の実現のため、数多くの研究が、教師ありCNNに対して注がれている。このネットワークは多くの用途(検出、分割、ラベル付け)に対して重要である。教師なし手法も研究されており、これは主に画像符号化、効率的な画像表示手法、教師あり手法にさらに供されるための前処理工程について有用性が示されている。制限付きボルツマンマシン(RBM)のような教師なし表示学習方法は、訓練用データから直接特徴記述を学習できるという点で、基本的なフィルタバンクの性能を凌駕し得る。RBMは、生成型学習目標で訓練される。これにより、ネットワークはラベル無しデータから標示を学習可能である。しかし、分類について最適な特徴が生成されるとは限らない。Van Tulderら(参照文献18)は、畳み込み分類用制限付きボルツマンマシンにおいて、生成型と識別型の学習目標の両方の利点を組み合わせるための研究を行った。この制限付きボルツマンマシンは、訓練用データの記述と、分類の両方に好ましいフィルタを学習する。これにより、学習目標を組み合わせることで、識別型、または生成型単独の構成を上回る性能が実現できることが示された。
CNNは、データに基づく、表現性が高い、階層状画像特徴の学習を可能とする。本稿に記載のように、このような特徴は、数多くの適用分野において極めて強力かつロバストであることが示された。このように表現豊かで、適切な分類を実現するためには、十分な訓練用データが必要となる。
十分なデータが得られない場合にも、以下の手法をとることができる。
転移学習:自然画像データセットまたは異なる医療分野により、事前に(教師あり)訓練されたCNNモデルを、対象の新たな医学的課題に利用する。一手法において、事前訓練されたCNNを入力画像に適用し、ネットワークの層から出力を抽出する。抽出された出力は、特徴として考慮され、別のパターン分類器の訓練に利用される。例えば、Barら(参照文献25,26)において、事前訓練されたCNNは、胸部病理の識別のための特徴生成に利用される。Ginnekenら(参照文献27)においては、CNNに基づく特徴と、手動で得られた特徴を統合することで、結節検出システムの性能向上が可能となっている。
ファインチューニング:対象の課題に対して、中程度の大きさのデータセットが存在しない場合、以下の手順等が推奨される。即ち、事前訓練されたCNNをネットワークの初期化に利用し、その後、対象の課題に対して新たなデータを使用して、一部(または全部)のネットワーク層に対する教師あり訓練がさらに実行される。転移学習とファインチューニングは、医療用撮像用途におけるDCNNの利用において、極めて重要な要素となる。
Shinら(参照文献17)およびTajbakhshら(参照文献28)は、ファインチューニングで事前訓練されたCNNを利用することで、最上の結果が得られることを示した。さらに、ディープファインチューニングの方が、シャロウファインチューニングよりも高性能を実現し、特定の用途領域にかかわらず、学習組が小さいほど(Tajbakhshら)、そしてあらゆるネットワーク構造において(Shinら(参照文献17))ファインチューニングの重要性が高いことを示した。GoogLeNet構造により、その他様々なより浅い構造と比較して、最先端を誇る縦隔リンパ節検出手法が実現された。Shinらを参照されたい。
グラウンドトゥルースデータは一般公開されていないため、各医療課題について収集するにはコスト面でも時間の面でも困難である。これにより、医療分野での適用が極めて困難となっている。クラウドソーシングにより、現実世界の画像に対する大規模データベースのアノテーションが可能となった。しかし、それを生物医学に活用するには、実際のアノテーション処理に対するより深い知識、即ち正確な定義が必要となる(Nguyenら(参照文献29)、McKennaら(参照文献30))。実際、専門家が行うような課題が、非専門家に外注されていることで、ユーザ間の不一致を生じる高ノイズのアノテーションに帰結してしまい得る。同様に、調整していない臨床結果データまたはユーザデータに基づくフィードバックを行うと、標準化、調整の欠如により問題のある結果となり得る。医療専門家と、非専門家とで知識を擦り合わせることは、多くの困難が伴う。即ち、情報源を如何に組み合わせるか、事前に確認された精度に応じて重みづけされた入力をどのように評価、反映するか等である。Albarqouniら(参照文献31)は、CNNに組み込まれるアグリゲーション層を組み合わせたネットワークを提案した。これにより、クラウドからの学習入力が、ネットワーク学習プロセスの一部として利用可能となる。この結果、クラウドアノテーションからのDCNN学習の機能について貴重な見識が得られた。医療分野におけるクラウドソーシングについての研究では、経験のない非専門家ユーザ群でも、医療専門家に引けを取らないという事実が示された。これは、Nguyenら(参照文献29)およびMcKennaら(参照文献30)により、放射線画像について観察されたもので、専門家の技量、知識、洞察力が結果に反映、関連しないことを示唆し得る。
Kallenbergら(参照文献32)には、マンモグラフィリスクスコアリングに対する、教師なし特徴学習が記載されている。この研究では、ラベル無しデータから、特徴階層を学習する方法が示されている。学習された特徴は、単純な分類器に入力され、2つの異なる課題、即ちi)乳房密度分割と、ii)マンモグラフィテクスチャのスコアリングに対する最先端の結果が得られる。モデル容量制御のため、疎性正規化器が導入される。これは、寿命および人口疎性の両方が組み込まれる。教師なし部分の畳み込み層は、オートエンコーダとして訓練される。教師あり部分では、(事前訓練された)重みとバイアス項が、ソフトマックス回帰により、ファインチューニングされる。
Yanら(参照文献33)では、画像分類用の多段ディープラーニングフレームワークが設計され、体の部位認識に適用される。事前訓練段階において、CNNはマルチインスタンス学習により訓練され、最も特徴的で、情報性の低いパッチがトレーニングスライスから抽出される。ブースティング段階では、事前訓練されたCNNがこれら局地的なパッチにより、画像分類用にさらにブーストされる。この方法で特筆すべきは、マルチインスタンスディープラーニングにより訓練され、最も特徴的で、情報性の低いパッチが自動で発見されたことである。したがって、手動のアノテーションは不要であった。
回帰的ネットワークは、医療用撮像分野ではあまり一般的ではない。Miaoら(参照文献34)では、リアルタイム二次元/三次元レジストレーション用のCNN回帰的手法が提示されている。元となるマッピングの回帰を簡潔化し、高強度の非線形モデリングでCNN回帰モデルを設計するため、3つのアルゴリズム式手法が提案されている。その結果により、識別ローカル(DL)手法が、その他2つの最先端加速強度式二次元/三次元レジストレーション方法よりも正確でロバストであることが確認されている。
Golkovら(参照文献35)は、拡散MRIデータ処理を、単一の最適化された工程に削減するためにDLを適用した、初期概念実証を提示している。このような変形により、高度なモデルからのスカラー尺度を、12分の1の走査時間で得られ、拡散モデルを利用しなくても異常が検出できることが実証された。拡散重みづけ信号と、微細構造組織の特性との間に有意な関係性が確認されている。Golkovら(参照文献35)は、そのような関係性が、DNNを使用することで実際に確認され得ることを実証した。ここで、モデルフィッティングで得られたスカラー尺度ではなくDWIが入力として直接使用されている。この研究では、健常組織とMS病変に対する、ボクセル探知の微細構造予測と、DWI値による、自動モデル不要の分割が示されている。拡散尖度は、12個のデータ点のみから測定され、神経突起方向分散および密度は、8個のデータ点のみから測定されている。これにより、迅速かつロバストな手順が実現され得る。これは、よりよい臨床ルーチンに寄与する。即ち、伝統的なデータ処理が、ディープラーニングで合理化できることが示される。
Kaggle(www.kaggle.org)は、眼底カラー画像からのDR検出および診断のコンペを開催した。結果、約8万枚の画像が利用可能となっている(www.kaggle.com/c/diabetic−retinopathy−detection)。多くの場合NNが利用されていた。
U.S.2014/0314288は3段階システムを開示する。このシステムは、輝度と視界の異なる眼底画像を分析し、DRの深刻度を生成する。画像の事前処理として、ヒストグラム平坦化、コントラスト強調、およびダイナミックレンジ正規化が挙げられる。第1段階として、グローバルフィルタおよび閾値技術の各種組み合わせにより、眼底画像から明るい、赤い領域が抽出される。視神経乳頭(OD)は、明るい病変に対応する構造的外観を持つ。そして、血管領域は、赤い病変に対応する画素強度特性を有する。したがって、既存の視神経乳頭検出アルゴリズムにより、視神経乳頭に対応する領域が明るい領域から除去される。さらに、単純な検出技術により、血管に対応する領域が、赤い領域から除去される。これにより、明るい候補領域を含む画像と、赤い候補領域を含む更なる画像が得られる。面積、周囲長、硬度、最小、最大、平均、標準偏差等の領域に基づく特徴が計算される。Adaboost(en.wikipedia.org/wiki/AdaBoost)により、使用される計30個の特徴が選択される。第2段階において、明るいおよび赤い候補領域が、2工程階層分類にかけられる。第1工程において、特徴に基づいて、明るいおよび赤い候補領域が非病変領域から分離される。第2工程では、分類された明るい候補領域が、さらに硬性白斑または綿花状白斑に分類される。そして分類された赤い病変領域が、さらに出血および微小動脈瘤に分類される。分類器として、GMM(en.wikipedia.org/wiki/Mixture_model#Gaussian_mixture_model)、kNN(en.wikipedia.org/wiki/K−nearest_neighbors_algorithm)、SVM(en.wikipedia.org/wiki/Support_vector_machine,en.wikipedia.org/wiki/Supervised_learning)等が利用され得る。第3段階として、画像ごとの明るいおよび赤い病変の数の合計により、DRの深刻度レベルが生成される。このようなシステムは、手動での評価を要する患者数を低減し、最もDRが深刻な患者に優先して、眼科的処置が行われることを目的としている。
U.S.8,098,907は、微小動脈瘤(MA)の自動検出を開示する。これは、血管、視神経乳頭および硬性白斑(HE)の位置のような情報をさらに考慮したものとなっている。この方法は、(i)画像を固定サイズのサブ領域に分割し、領域強調/正規化および適応分析を行うことと、(ii)血管、視神経乳頭および硬性白斑(HE)検出と、(iii)これら(i)および(ii)の結果の組み合わせにより実行される。MA検出は、トップハットフィルタリングと、局地適応的閾値処理に基づく。
U.S.2015/0104087は、自動でどちら側の眼底画像か検出すること、即ち、左目、右目か、そしてその視界の検出、さらに品質評価を開示する。この技術は、視神経乳頭および網膜黄斑の位置や、網膜血管構造の対称性の有無のような生理学的特徴を利用する。視界やどちら側かを検出することは、多チャンネルまたは単一チャンネルで実現される。このアルゴリズムは、視神経乳頭テンプレートを利用した、正規化二次元断面相関により、視神経乳頭を検出するものである。高解像度画像に対しては、血管構造密度により、どちら側か、そして視界が判断される。各画像の品質は、画像内の血管対称性を分析することで評価される。これは、二値分割画像を血管に対して取得して、特徴を抽出し、分類アルゴリズムを適用することで実現される。品質は1から5のレベルで評価される。血管の二値化については、ウェーブレットとエッジ位置特定が提案される。その後、モルフォロジー演算が実行される。特徴検出のため、画像は20個の矩形ウィンドウに分割され、ウィンドウ内の非ゼロ画素の数として、局地的血管密度(LVD)が計算される。20個のLVDが、全体的血管密度(GVD)で正規化される。GVDは、レベル5の分割二値画像の非ゼロ画素の数として計算される。そして、画像のLVDおよびGVDにより特徴ベクトルが形成される。分類に対してはSVMが提案される。画像レジストレーションおよび上記と同様の技術により、つなぎ合わされた重複画像の品質レベルを評価できる。
U.S.8,879,813、U.S.8,885,901、U.S.9,002,085、U.S.9,008,391、U.S.2015/0110348、U.S.2015/0110368、U.S.2015/0110370、U.S.2015/0110372、U.S.2017/0039412、およびU.S.2017/0039689は、診断システムの様々な態様を開示する。網膜画像は、中央値正規化により強調されてもよい。これより、画像が各画素位置にて、局地背景推定により局地的に強調される。中央値フィルタリング/収縮/膨張等により、網膜画像中のアクティブ画素が検出される。基本的に、これにより網膜ディスクが検出され、境界付近の画素が削除される。関心領域が検出される。例えば、2つのモルフォロジーフィルタリング済み画像を、2つの異なる種類およびサイズの、幾何学形状局地領域について、モルフォロジーフィルタで計算し、その間の差分をとることで、局地領域の記述子が抽出される。画像品質を、コンピュータビジョン技術で評価して、レベル付けの適切さを評価する。画像は、疾患用に自動でスクリーニングされる。患者が受診に訪れるたびに、画像に基づく病変バイオマーカーが自動分析され、受診と受診の間に画像レジストレーションが実行される。病変の変化と、解剖学的構造が計算され、統計的に数値化される。このようにして計算された統計は、病変の進行、早期発見および/または処置の効果の観察に利用可能な画像に基づくバイオマーカーとなる。
U.S.8,879,813は、網膜画像におけるアクティブ画素の自動検出を開示する。具体的には、網膜画像が評価され、異なるウィンドウサイズの2つの中央値フィルタリング済み画像が生成され、それらフィルタリング済み画像の差分画像が生成され、二値画像が生成される。
U.S.8,885,901は、網膜画像の強調を開示している。具体的には、網膜画像を取得し、中央値フィルタをかけて、元画像の値を、元画像およびフィルタリング済み画像の値に基づいて変調する。強調画像は、健康状態の検出に利用される。
U.S.9,002,085およびU.S.2015/0110372は、網膜画像内の局地領域の記述子を生成することを開示する。具体的には、網膜画像を取得し、2つのモルフォロジーフィルタリング済み画像を生成し(一方は円形/矩形ポリゴンウィンドウによるもので、他方は細長/楕円形ウィンドウによる)、フィルタリング済み画像から差分値を生成し、それを網膜画像の画素記述子に利用する。
U.S.9,008,391およびU.S.2015/0110368は、患者の網膜画像を取得することを開示する。各画像は、関心領域を含む画素サブセットをアクティブと指定し、画素レベル記述子を計算し、教師あり学習により記述子から画素単位分類を実現し、第2の記述子を計算し、教師あり学習により複数の画素について第2の分類を実現する。
U.S.2015/0110348およびU.S.2017/0039412は、網膜画像内の関心領域(ROI)の自動検出を開示する。具体的には、網膜画像を取得し、多重規模モルフォロジーフィルタバンク分析により1つ以上の所期の特性を有する領域を抽出し、二値マップを記憶する。
U.S.2015/0110370およびU.S.2017/0039689は、網膜画像の強調を開示する。具体的には、眼底画像を取得し、単一または複数の尺度で背景を推定し、第1画素位置での強度を、背景画像における同位置での強度に基づいて適応的にスケーリングする。
(参照文献)
[1]MIT Technol.Rev.,2013.Available:www.technologyreview.com/s/513696/deep−learning
[2]Y.LeCun,Y.Bengio,and G.Hinton,“Deep learning,”Nature,vol.521,No.7553,pp.436−444,2015.
[3]J.Schmidhuber,“Deep learning in neural networks:An overview,”Neural Netw.,vol.61,pp.85−117,2015.
[4]B.Sahiner et al.,“Classification of mass and normal breast tissue:A convolution neural network classifier with spatial domain and texture images,”IEEE Trans.Med.Imag.,vol.15,No.5,pp.598−610,Oct.1996.
[5]S.C.B.Lo,J.S.J.Lin,M.T.Freedman,and S.K.Mun,“Computerassisted diagnosis of lung nodule detection using artificial convolution neural−network,”Proc.SPIE Med.Imag.,Image Process.,vol.1898,pp.859−869,1993.
[6]H.−P.Chan,S.−C.Lo,B.Sahiner,K.L.Lam,and M.A.Helvie,“Computer−aided detection of mammographic microcalcifications:Pattern recognition with an artificial neural network,”Med.Phys.,vol.22,No.10,pp.1555−67,1995.
[7]S.Hochreiter and J.Schmidhuber,“Long short−term memory,”Neural Comput.,vol.9,No.8,pp.1735−1780,1997.
[8]G.E.Hinton,S.Osindero,and Y.W.Teh,“A fast learning algorithm for deep belief nets,”Neural Comput.,vol.18,No.7,pp.1527−1554,2006.
[9]D.Castano−Diez,D.Moser,A.Schoenegger,S.Pruggnaller,and A.S.Frangakis,“Performance evaluation of image processing algorithms on the GPU,”J.Struct.Biol.,vol.164,No.1,pp.153−160,2008.
[10]A.Eklund,P.Dufort,D.Forsberg,and S.M.LaConte,“Medical image processing on the GPU−Past,present and future,”Med.Image Anal.,vol.17,No.8,pp.1073−94,2013.
[11]B.van Ginneken,C.M.Schaefer−Prokop,and M.Prokop,“Computeraided diagnosis:How to move from the laboratory to the clinic,”Radiol.,vol.261,No.3,pp.719−732,2011.
[12]A.Setio et al.,“Pulmonary nodule detection in CT images using multiview convolutional networks,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1160−1169,May 2016.
[13]H.Roth et al.,“Improving computer−aided detection using convolutional neural networks and random view aggregation,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1170−1181,May 2016.
[14]Q.Dou et al.,“Automatic detection of cerebral microbleeds from MR images via 3D convolutional neural networks,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1182−1195,May 2016.
[15]K.Sirinukunwattana et al.,“Locality sensitive deep learning for detection and classification of nuclei in routine colon cancer histology images,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1196−1206,May 2016.
[16]M.Anthimopoulos,S.Christodoulidis,A.Christe,and S.Mougiakakou,“Lung pattern classification for interstitial lung diseases using a deep convolutional neural network,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1207−1216,May 2016.
[17]H.−C.Shin et al.,“Deep convolutional neural networks for computeraided detection: CNN architectures, dataset characteristics and transfer learning,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1285−1298,May 2016.
[18]G.van Tulder and M.de Bruijne,“Combining generative and discriminative representation learning in convolutional restricted Boltzmann machines,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1262−1272,May 2016.
[19]A.Depeursinge et al.,Comput.Med.Imag.Graph.,vol.36,No.3,pp.227−238,2012.
[20]F.Ghesu et al.,“Marginal space deep learning:Efficient architecture for volumetric image parsing,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1217−1228,May 2016.
[21]T.Brosch et al.,“Deep 3D convolutional encoder networks with shortcuts for multiscale feature integration applied to multiple sclerosis lesion segmentation,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1229−1239,May 2016.
[22]S.Pereira,A.Pinto,V.Alves,and C.Silva,“Brain tumor segmentation using convolutional neural networks in MRI images,” IEEE Trans.Med.Imag.,vol.35,No.5,pp.1240−1251,May 2016.
[23]P.Moeskops et al.,“Automatic segmentation of MR brain images with a convolutional neural network,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1252−1261,May 2016.
[24]M.van Grinsven,B.van Ginneken,C.Hoyng,T.Theelen,and C.Sanchez,“Fast convolutional neural network training using selective data sampling: Application to hemorrhage detection in color fundus images,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1273−1284,May 2016.
[25]Y.Bar,I.Diamant,L.Wolf,and H.Greenspan,“Deep learning with non−medical training used for chest pathology identification,”Proc.SPIE Med.Imag.Computer−Aided Diagnosis,vol.9414,2015.
[26]Y.Bar,I.Diamant,L.Wolf,S.Lieberman,E.Konen,and H.Greenspan,“Chest pathology detection using deep learning with non−medical training,”in Proc.IEEE 12th Int.Symp.Biomed.Imag.,2015,pp.294−297.
[27]B.van Ginneken,A.A.Setio,C.Jacobs,and F.Ciompi,“Off−the−shelf convolutional neural network features for pulmonary nodule detection in computed tomography scans,”in Proc.IEEE 12th Int.Symp.Biomed.Imag.,2015,pp.286−289.
[28]N.Tajbakhsh et al.,“Convolutional neural networks for medical image analysis: Full training or fine tuning?,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1299−1312,May 2016.
[29]T.B.Nguyen et al.,“Distributed human intelligence for colonic polyp classification in computer−aided detection for CT colonography,”Radiology,vol.262,No.3,pp.824−833,2012.
[30]M.T.McKenna et al.,“Strategies for improved interpretation of computer−aided detections for CT colonography utilizing distributed human intelligence,”Med.Image Anal.,No.6,pp.1280−1292,2012.
[31]S.Albarqouni,C.Baur,F.Achilles,V.Belagiannis,S.Demirci,and N.Navab, “Agg−Net:Deep learning from crowds for mitosis detection in breast cancer histology images,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1313−1321,May 2016.
[32]M.Kallenberg et al.,“Unsupervised deep learning applied to breast density segmentation and mammographic risk scoring,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1322−1331,May 2016.
[33]Z.Yan et al.,“Multi−instance deep learning: Discover discriminative local anatomies for bodypart recognition,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1332−1343,May 2016.
[34]S.Miao,Z.J.Wang,and R.Liao,“A CNN regression approach for real−time 2D/3D registration,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1352−1363,May 2016.
[35] V.Golkov et al.,“q−Space deep learning: Twelve−fold shorter and model free diffusion MRI scans,”IEEE Trans.Med.Imag.,vol.35,No.5,pp.1344−1351,May 2016.
[36]K.He,X.Zhang,S.Ren,and J.Sun,Deep residual learning for image recognition ArXiv,2015,arXiv:1512.03385,to be published
[37]H.R.Roth et al.,“A new 2.5 d representation for lymph node detection in CT,”Cancer Imag.Arch.,2015 dx.doi.org/10.7937/K9/TCIA.2015.AQIIDCNM
[38]H.R.Roth et al.,“Data from pancreas−CT,”Cancer Imag.Arch.,2016[Online].Available:dx.doi.org/10.7937/K9/TCIA.2016.tNB1kqBU
(参照)U.S.5784162;6088099;6198532;6276798;6419361;6556853;6895264;6943153;6992775;7433532;7474775;7712898;7931902;8098907;8114843;8194936;8207396;8303115;8309350;8340437;8705826;8787638;8879813;8885901;8896682;9002085;9008391;9097707;9545196;20010033364;20020052551;20020194630;20030129164;20040221855;20040254154;20060257031;20070128662;20080124344;20100061601;20100172871;20100220906;20100302507;20110097330;20110173708;20110242306;20120190094;20120257164;20120272341;20130108131;20130116670;20130137113;20130301889;20140023692;20140199277;20140276025;20150110348;20150110368;20150110370;20150110372;20150124216;20150224193;20150238767;20150379708;20160120404;20160217586;20170039412;20170039689;20170046616
本発明は、1つ以上の血管画像に基づいて、血管疾患の存在を自動で推定するシステムおよび方法を提供する。画像は好ましくは眼底画像である。したがって、血管疾患は網膜疾患、または網膜症となる。この技術の適用対象としては、特にDRのスクリーニングが挙げられる。
基本的な手段は網膜画像に限られない。ただし、眼底画像における網膜血管系の視認性は、本技術にとって特に好ましい。しかし、その他血管系および微小血管系も分析の対象となる。例えば各種条件下で、保粘膜および舌下粘膜の毛細血管、超粘膜、強膜等が撮像されてもよい。
網膜症は以下の通りに診断され得る。(H35.0)高血圧網膜症−長期高血圧による血管破裂en.wikipedia.org/wiki/Hypertensive_retinopathy;(H35.0/E10−E14)糖尿病性網膜症−糖尿病の合併症で生じる網膜の損傷。最終的に失明に至る可能性があるen.wikipedia.org/wiki/Diabetic_retinopathy;(H35.0−H35.2)網膜症−網膜の非炎症性損傷全般を指すen.wikipedia.org/wiki/Retinopathy;(H35.1)未熟児網膜症−未熟児瘢痕および網膜剥離の網膜症en.wikipedia.org/wiki/Retinopathy_of_prematurity;(H35.3)加齢黄斑変性−黄斑の光感受性細胞の機能不全で、時間が経つと機能しなくなる;(H35.3)黄斑変性−黄斑変性による中心視力の喪失(例えば、ブルズアイ黄斑症、クロロキン網膜症en.wikipedia.org/wiki/Chloroquine_retinopathy)en.wikipedia.org/wiki/Macular_degeneration;(H35.3)網膜前膜−透明な層が張り、網膜上で収縮するen.wikipedia.org/wiki/Epiretinal_membrane;(H35.4)周辺部網膜変性en.wikipedia.org/wiki/Lattice_degeneration;(H35.5)遺伝性網膜ジストロフィen.wikipedia.org/wiki/Progressive_retinal_atrophy;(H35.5)網膜色素変性症の遺伝的障害;夜盲症とその後の視野狭窄en.wikipedia.org/wiki/Retinitis_pigmentosa;(H35.6)網膜出血;(H35.7)網膜層の分離en.wikipedia.org/wiki/Retinal_haemorrhage(例えば、中心性漿液性網膜症en.wikipedia.org/wiki/Central_serous_retinopathy,網膜剥離:網膜色素上皮の剥離en.wikipedia.org/wiki/Retinal_detachment)en.wikipedia.org/wiki/Retina#Anatomy_of_vertebrate_retina;(H35.8)その他の明示された網膜障害;(H35.81)黄斑浮腫−腫れ斑による中心視のゆがみen.wikipedia.org/wiki/Macular_edema、(H35.9)網膜障害、詳細不明。世界保健機関ICD−10コード:眼及び付属器の疾患(H00−H59)。
画像は、様々な疾患に対して並行に分析されてもよい。例えば、訓練用データは、少なくとも綿花状白斑(浮腫を伴う網膜虚血の領域)、硬性白斑(脂肪沈着物)、微小動脈瘤(網膜上に小さな赤い点として現れる)、損傷した血管壁からの小さな炎上出血、血管新生(眼底での新血管形成の兆候)、糖尿病(水晶体の代謝を妨げる過剰なグルコース、およびその他の臨床徴候による早期白内障の原因ともなる)の兆候を符号化したものであってもよい。patient.info/doctor/eye−in−systemic−disease;www.ncbi.nlm.nih.gov/books/NBK221/を参照されたい。
より具体的には、本発明は、網膜眼底画像内のDR病変の自動検出および位置特定ならびに網膜眼底画像のDR自動分類のシステムおよび方法を含む。システムは、患者の眼底画像における専門家による領域の分類を示す訓練用データに基づいて、DRを認識する。
いくつかの実施形態では、コンピュータにより実施される方法は、網膜画像を分析し、当該網膜画像内で関心領域を検出することを含む。いくつかの実施形態では、方法はさらに、1つ以上の関心領域の検出に基づいて、網膜画像を分類することを含む。
本技術は、網膜画像において自動で網膜症を検出するという課題に対する新たな手法を提供する。具体的には、人の専門家が網膜症を全てマーク付け(信頼性高いと認定されたアノテーション)した何千もの網膜写真が収集される。アノテーションは、医学的分類の領域評価や、画像品質であり得る。画像そのものが光学的原則に基づくので、画像品質は少なくとも一部、教師なしNNで判定されてもよい。ただし、場合によっては、人の専門家が関わったほうが好ましくあり得る。
アノテーションは、NN(CNNまたはDNNであってもよい)の訓練に使用される。訓練用データの領域分類によって、本技術に係るネットワークは、分類されていなかった画像の各領域の分類を判定する。これは、全体的に未知の画像を分類する技術と対照的である。したがって、NNそのものがラベルそしてラベル付けされるものの両方に、基本的に依存することが明らかである。
訓練用データは、領域分類がラベルとして関連付けられた網膜画像全体であってもよい。ただしこの場合、入力データ組の自由度が増してしまい、対応するネットワーク性能向上が得られないまま、問題が複雑化してしまい得る。即ち、網膜の一領域内の網膜症は、別の領域の網膜症に極めて類似しているため、訓練は領域単位に限定されない。当然、領域単位の訓練もあり得るので、本技術はこれを無視や否定することはなく、領域毎の推定でも十分な性能が実現できるとするものである。例えば、1000×1000入力画像は106画素を含み、簡略化しないと、NNには、少なくとも106入力ニューロンが必要となる。一方、64×64画素領域でも病理特定に十分な精度が得られる上、この場合に必要な入力は少なくとも4096個で、約250分の1である。同様に、領域毎の推定により、領域間差分を処理するように設計されたネットワークの少なくとも隠れ層全体も低減可能となり得る。
処理前に、画像は色、輝度等に関して、決定論的統計またはルールベースのプロセスで正規化されてもよい。この正規化はNN内で実現されてもよいが、ネットワークとは独立的に実現されてもよい。したがって、よりシンプルなネットワークが実現される。
場合によっては、領域は分割後にまたは分割処理の一環として正規化されてもよい。これにより、統計的またはルールベース処理により分割が実現可能である。この処理は、形状または空間以外に基づく領域間境界を構築できる。前述のように、NN前の正規化により、NNをシンプルに実施できる。それでいて、NNの正規化前特徴の学習量が低減できる。正規化処理により画像から抽出され得る特徴に対する専門家の判定を利用することで、正規化されていない特徴に対するNNの特異度を上げることができる。
各網膜画像は指紋同様、極めて個別性が高い。したがって、データセット内に、画像全体が完全一致するものは含まれない。しかし、この技術が処理対象とする初期疾患は、小さく、局地的(高解像度画像において、幅数画素)という特徴がある。事実、近年ようやくデジタル写真が臨床現場でのスクリーニングに利用されたのは、疾患の超初期兆候を視認可能となる程の解像度が近年ようやく実現されたためである。したがって、実際のデータセットは実際の網膜写真から得られた小パッチ組であり、本技術のシステムは、新たな画像からの小パッチを、データセットのパッチと比較する。
入力された新たな画像には以下の高次処理が施される。
カメラから元画像が受信され、円形マスクにより、1024×1024の標準サイズにされる。
NNを実施して、領域に対して動作する分類器が実現される。好ましくは、ラベル付きデータセット画像に適用されたのと同じ前処理(領域への分割を含む)が未知の画像に適用される。正規化後の分割領域は、次に訓練用データに基づいて分類される。そして、分類されたものが、未知画像全体に反映される。これにより未知画像がアノテーションされる。その画像が診断されるか、または分類組全体に基づいて評価される。
評価は信頼性を伴うものであってもよい。即ち、第一種(偽陽性)および/または第二種(偽陰性)過誤や、その他種類の過誤をカバーするものであってもよい。評価の信頼性は、元の画像品質(輝度、焦点、コントラスト、色等)を表すものであってもよい。ここで、例えば、画像の一部の品質が悪い場合に、評価に利用可能な部分を分析してもよい。評価に利用可能な部分とは、信頼性を伴う評価である。即ち、偽陽性過誤診断の可能性および、偽陰性過誤診断の可能性を包含するものである。
任意で結果は、全体的推奨に照合されてもよい。或いは、特定の関心領域に対する付記とともに、詳細を提供してもよい。
この技術は、人の専門家が適切に評価した十分な画像を収集して、データセットを形成することを要する。偽陽性、偽陰性の結果とならないような、十分な精度および確度を持つマッチングアルゴリズムが求められる。
網膜カメラモデルは、様々な解像度の画像を生成する。典型的な網膜病変の大きさも変動し得る。数多くの実験と、専門家や英国国民保険サービス規格文書からの情報収集を経て、十分な特徴が保たれる最低解像度は、1024×1024画素であると結論づけられた。さらに、(上述の画像で)病変が特定可能な最小パッチは、32×32画素であった。
画像は、血糖値、年齢、性別、体重、病歴等の追加の医療データと併せて提供され得る。これらも分類の一部として使用される場合がある。ただし、場合によっては、画像を分類し、臨床医が他の患者データを診断および/または予後に統合できるようにすることが好ましい。実際、臨床医が分類器の出力を使用する限り、画像を分類するために画像に外因性の臨床データを使用すると、同じ要素が再び検討され得る下流側に、過剰なバイアスがかかるというリスクが生じる。
上述のような第一種および第二種過誤との間に必要な正確なバランスは、システムの用途に依存する。病院の場合、第二種過誤を最小限に抑えるのが最重要であり得る。一方で、眼鏡技師による通常のスクリーニングの場合は、第一種過誤を最小限に抑えて、顧客に無用な心配をかけないことが最重要であり得る。
なお、システムは必ずしも完全自動診断に限定されない。即ち、第二種過誤を最小限に抑える場合、偽陽性は、人の評価者のような別の手段、或いは、例えばより多くの時間および/または資源を必要とする代替または補足の自動分析によって再スクリーニングを実行してもよい。例えば、ハッシュエラーが性能品質の劣化を招く場合、データ削減量を抑えた自動分析で補ってもよい。どのようなシステムでも、本質的な性能の良し悪しにかかわらず、これら2種類の過誤間で様々なトレードオフを実現するように「調整」することができる。
このアルゴリズムは、網膜画像の疾患兆候のある領域にマーキングするのに非常に効果的であることが確認されており、画像処理/領域マッチング技術が効果的かつ効率的であることが実証された。
広視野および超広視野の網膜画像は、通常網膜眼底画像で撮像される、単一画像内の45から50度を超える網膜視野で撮影した画像である。このような画像は、特殊なカメラハードウェアを使用するか、異なる複数の視野による網膜画像を使用してモンタージュを作成することによって得られる。本明細書で説明されるシステムおよび方法は、広視野および超広視野画像に適用可能である。
フルオレセイン血管造影法では、蛍光トレーサー染料を注入し、波長490ナノメートルの光で網膜を照らすことによって放出される蛍光を測定する血管造影法が実行される。当該染料が血液中に存在するため、蛍光血管造影画像上では網膜の血管構造と病変が強調される。本明細書で説明されるシステムおよび方法は、フルオレセイン血管造影画像に適用可能である。
したがって、目的は、血管画像を分類する方法であって、血管画像を受信することと、血管画像を正規化することと、正規化血管画像を複数の領域に分割することと、NNを含む少なくとも1つの分類器により、複数の領域について自動で血管障害ベクトルを判定することと、判定された血管障害ベクトルに基づいて、自動で複数の領域それぞれをアノテーションすることと、少なくともアノテーションに基づいて、受信した血管画像を自動で評価することと、を含み、NNは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が血管画像の各領域につけたアノテーションに基づいて訓練される、方法を提供することである。
更なる目的は、血管画像を分類するシステムであって、少なくとも1つの血管画像を受信するように構成された入力と、少なくとも1つの客観的分類基準に応じて、複数の網膜画像の各領域の専門家によるアノテーションに少なくとも基づいて訓練されるNNを定義する情報を記憶するように構成されるメモリと、少なくとも1つの自動プロセッサであって、血管画像を正規化することと、正規化血管画像を複数の領域に分割することと、定義されたNNを含む少なくとも1つの分類器により、複数の領域について血管障害ベクトルを判定することと、判定された眼病ベクトルに基づいて、複数の領域それぞれをアノテーションすることと、少なくともアノテーションに基づいて、受信した血管画像を評価することとを実行するように構成された自動プロセッサと、評価を通信するように構成された出力とを備えるシステムを提供することである。
更なる目的は、少なくとも1つの自動プロセッサを制御するための指示を記憶するコンピュータ可読媒体であって、血管画像を受信するための指示と、受信した血管画像を正規化するための指示と、正規化血管画像を複数の領域に分割するための指示と、NNを有する少なくとも1つの分類器により、複数の領域について眼病ベクトルを判定するための指示と、判定された眼病ベクトルに基づいて、自動で複数の領域それぞれをアノテーションするための指示と、少なくともアノテーションに基づいて、受信した血管画像を自動で評価するための指示と、を含み、NNは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が血管画像の各領域につけたアノテーションに基づいて訓練される、媒体を提供することである。
更なる目的は、網膜症の存在について、網膜画像を分類する方法であって、眼底画像を受信することと、少なくとも三色刺激画像空間に対して眼底画像を自動で正規化することと、正規化眼画像を、少なくとも1つの分割ルールに基づいて、複数の領域に自動で分割することと、少なくとも専門家網膜アノテーションを有する眼底画像組による訓練に基づいて訓練されたNNを有する少なくとも1つの分類器により、複数の領域について網膜症ベクトルを自動で判定することと、判定した眼病ベクトルに基づいて、自動で複数の領域それぞれにアノテーションを行うことと、少なくともアノテーションに基づいて、受信した眼画像の網膜症の程度を自動で評価することと、自動で判定した網膜症の程度を出力することと、を含む方法を提供することである。
別の目的は、眼画像を分類する方法であって、眼画像を受信して正規化することと、正規化眼画像を複数の領域に分割することと、NNを有する1つの分類器により、複数の領域について自動で眼病ベクトルを判定することと、判定された眼病ベクトルに基づいて、自動で複数の領域それぞれをアノテーションすることと、少なくともアノテーションに基づいて、受信した眼画像を自動で評価することと、を含み、NNは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が眼画像の各領域につけたアノテーションに基づいて訓練される、方法を提供することである。
更なる目的は、眼画像を分類するシステムであって、少なくとも1つの眼画像を受信するように構成された入力と、少なくとも1つの客観的分類基準に応じて、複数の網膜画像の各領域の専門家によるアノテーションに少なくとも基づいて訓練されるNNを定義する情報を記憶するように構成されるメモリと、少なくとも1つの自動プロセッサであって、眼画像を正規化することと、正規化眼画像を複数の領域に分割することと、定義されたNNを有する少なくとも1つの分類器により、複数の領域について眼病ベクトルを判定することと、判定された眼病ベクトルに基づいて、複数の領域それぞれをアノテーションすることと、少なくともアノテーションに基づいて、受信した眼画像を評価することとを実行するように構成された自動プロセッサと、評価を通信するように構成された出力とを備えるシステムを提供することである。
受信した眼画像を正規化するため、自動プロセッサを制御する指示と、正規化眼画像を複数の領域に分割する指示と、NNを含む少なくとも1つの分類器により、複数の領域について眼病ベクトルを判定する指示と、判定された眼病ベクトルに基づいて、複数の領域それぞれを自動でアノテーションする指示と、少なくともアノテーションに基づいて、受信した眼画像を自動で評価する指示とを記憶するコンピュータ可読媒体であって、NNは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が眼画像の各領域につけたアノテーションに基づいて訓練される媒体を提供することを目的とする。
更なる目的は、病理学画像を分類する方法であって、画像を受信することと、画像を正規化することと、正規化画像を複数の領域に分割することと、NNを有する少なくとも1つの分類器により、複数の領域について自動で疾患ベクトルを判定することと、判定された疾患ベクトルに基づいて、自動で複数の領域それぞれをアノテーションすることと、少なくともアノテーションに基づいて、受信した画像を自動で評価することと、を含み、NNは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が画像の各領域につけたアノテーションに基づいて訓練される、方法を提供することである。
NNは、複数の隠れ層、DNN、RNN、および/またはCNNを含んでもよい。方法はさらに、血管画像とは異なるタイミングで第2の血管画像を受信することと、血管画像の経時変化に少なくとも応じた情報を出力することとをさらに含んでもよい。
NNは、網膜画像の各領域の品質に対する専門家によるアノテーションに少なくとも基づいてさらに訓練してもよい。
少なくとも1つの分類器は、例えば多クラスサポートベクターマシン分類器または勾配ブースティング分類器を含んでもよい。
少なくとも1つの分類器は、複数の各領域の眼病ベクトルまたは複数の連続領域に基づいて、各領域を分類してもよい。
分類器は、DRの兆候を持つ領域を判定する分類器を含んでもよい。
評価することはDRのレベルを評価することを含んでもよい。レベルは、レベルが正確である可能性を伴ってもよい。正確である可能性は、第一種過誤の可能性と、第二種過誤の可能性を含んでもよい。
方法は、血管画像の領域の分類を示す画像を出力することをさらに含んでもよい。
血管画像は、眼底画像を含んでもよい。アノテーションは、各眼底画像の各領域内のDRの兆候を含んでもよい。
図1は、本発明の実施形態に係るフロチャートを示す。 図2Aから2Dは網膜画像を示す。 図3は、網膜画像と、それぞれその円形領域を示す。 図4は、網膜画像と、それぞれその正方形領域を示す。 図5は、本発明の各実施形態に係る各種フロチャートを示す。 図6は、本発明の各実施形態に係る各種フロチャートを示す。 図7は、本発明の各実施形態に係る各種フロチャートを示す。 図8は、本発明の各実施形態に係る各種フロチャートを示す。 図9は、本発明の各実施形態に係る各種フロチャートを示す。 図10は、本発明の各実施形態に係る各種フロチャートを示す。 図11は、本発明の各実施形態に係る各種フロチャートを示す。 図12は、本発明の各実施形態に係る各種フロチャートを示す。 図13は、本発明の各実施形態に係る各種フロチャートを示す。 図14は、本発明の各実施形態に係る各種フロチャートを示す。 図15は、本発明の各実施形態に係る各種フロチャートを示す。 図16は、本発明の各実施形態に係る各種フロチャートを示す。 図17は、本発明の各実施形態に係る各種フロチャートを示す。 図18は、本発明の各実施形態に係る各種フロチャートを示す。 図19は、本発明の各実施形態に係る各種フロチャートを示す。 図20は、本発明の各実施形態に係る各種フロチャートを示す。 図21は、本発明の各実施形態に係る各種フロチャートを示す。 図22は、本発明の各実施形態に係る各種フロチャートを示す。 図23は、本発明の各実施形態に係る各種フロチャートを示す。
図1は本発明の実施形態を示す。網膜眼底画像Iが入力されると、ブロック100において画像正規化を実行して、正規化画像I′を生成する。この正規化工程は、Iにおいて網膜ディスクを特定し、一定のサイズにリサイズすることを含む。具体的には、p×p画素の画像において、直径p画素となるようにする。なお、一般的にp=1024となる。これを図2(A)から2(D)に示す。図2(A)と2(C)は、正規化前網膜眼底画像Iを示し、図2(B)および2(D)は対応する正規化網膜眼底画像I'を示す。網膜ディスクの自動検出は、単純な作業である。即ち、均一に比較的暗い背景において、比較的明るいディスクを検出すればよい。したがって、エッジ検出や、円形ハフ変換等、多数の既知の画像処理技術が利用できる。このような技術の詳細については、Russ,John C.The image processing handbook.CRC press,2016を参照されたい。
続いて図1のブロック110において、正規化画像I'から網膜領域が抽出される。この処理は、網膜ディスクから小さな網膜領域を抽出することを含む。本発明のいくつかの実施形態において、これら領域は網膜ディスク内全体となるが、本発明の別の実施形態では、これら領域は部分的に網膜ディスク内の一部であってもよい。具体的には、例えば網膜画素の内の少なくとも50%または75%を占めるものであってもよい。本発明のいくつかの実施形態において、これら領域は、直径q画素の円形であってもよい(q=32またはq=64)。本発明の別の実施形態においては、これら領域は大きさがq×q画素の正方形であってもよい(q=32またはq=64)。本発明の更なる実施形態においては、例えば六角形のようなさらに異なる形状の領域が利用されてもよい。
実際、同領域は任意の形状の連続的な画素領域であればよい。ただし、場合によっては非連続的な形状の領域となり得る。さらに、画素以外の領域についての画像情報も分析対象となり得る。
図3は、正規化画像I'から抽出された2つの円形網膜領域Rα=R(x,y)(中心
点:(x,y)=(579,209))、およびRβ=R(x,y)(中心点:(x,y)=(603,226))を示す。図4は、正規化画像I'から抽出された2つの正方形網膜領域Rα=R(x,y)(中心点:(x,y)=(757,327))、およびRβ=R(x,y)(中心点:(x,y)=(782,363))を示す。いずれの場合でも、領域は通常、矩形パターンで、ある程度重複して抽出される。例えば、左上隅からs個画素の移動幅でずらしながら、順次領域を抽出する。右上隅に到達すると、s個画素分下にずらして、さらに左端から右端へとs個画素の移動幅でずらしながら、順次領域を抽出する(以後同様)。通常、sは小さい値をとる。即ち、典型的には、1から数画素の間である。本発明の別の実施形態では、sの値はさらに大きくてもよく、領域同士も重複しないように抽出してもよい。いずれの場合でも、ブロック110によりn個の網膜領域R0,…,Rn-1組が生成される。これらは、まとめてRと称される。さらに正規化画像I'の(x,y)を中心点としたRk=R(x,y)とも称される。
次に、図1のブロック120の動作を考える。図1のブロック120では、網膜領域全てを個別に、そして各領域Rkについて分析する。そして、DR病変分類ベクトル
が生成される。ここで、lは、ブロック120の分類器を訓練した分類数を示し、
は領域kの分類cに対するクラス所属確率を示す。
第1の例として、本発明の一実施形態において、ブロック120における分類器は、領域をl=2クラスに分類するように訓練される。ここで、c=0は「健常」クラスラベルとなり、c=1は「病気」クラスラベルとなる。
第2の例として、本発明の別の実施形態において、ブロック120における分類器は、領域をl>2クラスに分類するように訓練されてもよい。ここで、c=0は「健常」クラスラベル、c=1は「微小動脈瘤」クラスラベル、c=2は「硬性白斑」クラスラベル、…、c=l−1は「レーザー傷」クラスラベルとなる。したがって、そのような実施形態では、網膜眼底画像に映り得る異なる様々な種類の病変、アーチファクト、傷に対して、それぞれ別個のクラスラベルが割り当てられる。
このように、本発明の異なる実施形態では、分類器が様々な数のクラスに対して訓練される。これにより、DR病変検出および/または検出されたDR病変特定性能が可能な限り引き上げられる。いずれの場合にも、本発明の様々な実施形態において、ブロック120の分類器は、
をクラス所属確率ではなく二値(0または1)で出力するように構成され得る。ここで、Dkはクラス確率分布ではなく、分類判定ベクトルとなる。
次に、図1のブロック120の構造をより詳細に検討する。本発明の一実施形態において、図1のブロック120は図5に示す構造を有する。ブロック500において、任意で網膜領域が正規化される。この任意の正規化は、例えば、平均減算(各網膜領域の各カラーチャネルの平均強度の値が0となる)、またはその他何らかの平均調整(各網膜領域の各カラーチャネルの平均強度が特定の値をとる)、または何らかの幾何正規化(例えば、回転正規化であって、各領域が、質量中心が特定の角度となるように回転される)を含む。即ち、各網膜領域Rkについて、ブロック500で正規化領域R'kが生成される。そして、ブロック510においてCNNを利用して、各領域R'kが分類され、分類ベクトルDkが生成される。CNNの訓練には、多数の訓練用サンプル、即ち訓練用ラベル付き領域が使用されているが、上述のように、本発明の様々な実施形態によると、2つのクラスのみ(例えば「健常」および「病気」)或いは2つ以上のクラス(例えば「健常」、「微小動脈瘤」、「硬性白斑」等)のみによりにより訓練される。そして、いずれの場合でも、CNNはクラス確率判定または分類判定ベクトルとしてDkを出力し得る。
図6は、図5のブロック510の構造をより詳細に示す。図6に示すように、各領域Rkは個別に処理されて、その分類ベクトルDkが生成される。図6に示すCNNは一般的な構造である。即ち、入力を学習済み畳み込みマスクにより畳み込む多くの畳み込み層と、プーリング層(即ち空間的サブサンプリング層)が交互に配置され、それに多数の全結合層が続き、領域Rkに対してDR病変分類ベクトルDkが生成される。CNNの原理(en.wikipedia.org/wiki/Convolutional_neural_network)は本明細書に記載しないが、それに関する数多くの刊行物が存在する(例えば、Krizhevsky,Alex,Ilya Sutskever,and Geoffrey E.Hinton.“Imagenet classification with deep convolutional neural networks.”Advances in neural information processing systems.2012)。
本発明の別の実施形態では、図1のブロック120は、図7に示す構造を有する。ブロック700は任意のもので、図5のブロック500と同様に動作する。即ち、各網膜領域Rkについて、ブロック700では正規化領域R'kが生成される。そして、ブロック7
10において、CNNを利用して、各領域R'kについて際立った特徴の高次ベクトルV
kが生成される。
図8は、ブロック710の構造を示す。図8に示すように、特徴ベクトルVkの生成に利用されるCNNは、図6に示すものと同一である。即ち、最後の全結合分類層の前の何れかの層からの出力として、Vkが得られる。そして、図7を再度参照すると、ブロック720において、分類器により特徴ベクトルVkが分類されて、分類ベクトルDkが生成される。ブロック720に利用される分類器は各種存在する。即ち、事前に訓練用ラベル付き領域組の特徴ベクトルで訓練された2クラスまたは多クラスサポートベクターマシン(SVM)または勾配ブースティング分類器(GBC)である。上述のように、本発明の異なる複数の実施形態では、2つまたはそれ以上のクラスで訓練された分類器が利用され、クラス確率判定または分類判定ベクトルとしてDkを出力してもよい。SVM(en.wikipedia.org/wiki/Support_ベクトル_machine)およびGBC(en.wikipedia.org/wiki/Gradient_boosting)の原理は本明細書に記載しないが、それに関する数多くの刊行物が存在する(例えば、Burges,Christopher JC.“A tutorial on support vector machines for pattern recognition.”Data mining and knowledge discovery 2.2(1998):121−167 and Schapire,Robert E.,and Yoav Freund.Boosting:Foundations and algorithms.MIT press,2012)。
本発明の別の実施形態において、図1のブロック120は、図9に示す構造を有する。ブロック900は任意のもので、図5のブロック500または図7のブロック700と同様に動作する。即ち、各網膜領域Rkについて、ブロック900では正規化領域R'kが生成される。そして、ブロック910では、線形判別分析(LDA en.wikipedia.org/wiki/Linear_discriminant_analysis)または主成分分析(PCA en.wikipedia.org/wiki/Principal_component_analysis)等の次元削減技術により、各領域R'kの次元削減表現を生成し、当該領域の特徴ベクトルVkが生成される。そしてブロック920において、分類器により各特徴ベクトルVkが分類されて、分類ベクトルDkが生成される。上述の実施形態同様、トレーニング用ラベル付き領域組の、特徴ベクトルで事前に訓練された二値または他クラスSVMまたは勾配ブースティング分類器が特に好ましい。上述のように、本発明の異なる複数の実施形態では、2つまたはそれ以上のクラスで訓練された分類器が利用され、クラス確率判定または分類判定ベクトルとしてDkを出力してもよい。
次に、図1のブロック130の動作を検討する。図1のブロック130では、網膜領域全てを個別に、そして各領域Rkについて分析する。そして、品質分類ベクトル
が生成される。ここで、mは、ブロック130の分類器を訓練した分類数を示し、
は領域kの分類zに対するクラス所属確率を示す。第1の例として、本発明の一実施形態において、ブロック130における分類器は、領域をm=2クラスに分類するように訓練されてもよい。ここで、z=0は「高品質」クラスラベルで、z=1は「低品質」クラスラベルである。第2の例として、本発明の別の実施形態では、ブロック130における分類器は、領域をm>2クラスに分類するように訓練されてもよい。ここで、z=0は「高品質」クラスラベルで、z=1は「輝度不良」クラスラベルで、z=2は「コントラスト不良」クラスラベルで、…、z=m−1は「レンズフレア」クラスラベルとなる。
したがって、そのような実施形態では、網膜眼底画像に映り得る異なる様々な種類の品質アーチファクトおよび撮像アーチファクト、それぞれ別個のクラスラベルが割り当てられる。本発明の各種実施形態では、分類器が様々な数のクラスに対して訓練される。これにより、網膜領域の品質評価性能が可能な限り引き上げられる。なお、領域品質分類器のクラス数は、DR病変分類器のクラス数選択から独立して選択される。いずれの場合にも、本発明の様々な実施形態において、ブロック130の分類器は、
をクラス所属確率ではなく二値(0または1)で出力するように構成され得る。ここで、Qkはクラス確率分布ではなく、分類判定ベクトルとなる。
次に、図1のブロック130の構造をより詳細に検討する。本発明の一実施形態において、図1のブロック130は、図10に示す構造を有する。ブロック1000において、任意で網膜領域が正規化される。この任意の正規化は、例えば強度または形状正規化を含む。即ち、各網膜領域Rkについて、ブロック1000で正規化領域R'kが生成される。そして、ブロック1010においてCNNを利用して、各領域R'kが品質に応じて分類され、品質分類ベクトルQkが生成される。CNNの訓練には、多数の訓練用サンプル、即ち訓練用ラベル付き領域が使用されているが、上述のように、本発明の様々な実施形態によると、2つのクラスのみ(例えば「高品質」および「低品質」)或いは2つ以上のクラス(例えば「高品質」、「照明不良」、「コントラスト不良」等)のみにより訓練される。そして、いずれの場合でも、CNNはクラス確率判定または分類判定ベクトルとしてQkを出力し得る。
図11は、ブロック1010のCNNの内部構造を示す。これは上述した図6のCNN構造と実質的に同じである。ただし、サイズ、幅、フィルタ数、学習済み畳み込みマスク、畳み込み層数、全結合層数等のネットワークパラメータが異なり得る。
本発明の別の実施形態では、図1のブロック130は、図12に示す構造を有する。ブロック1200は任意のもので、図10のブロック1000と同様に動作する。即ち、各網膜領域Rkについて、ブロック1200では正規化領域R'kが生成される。そして、ブロック1210において、CNNを利用して、各領域R'kについて際立った特徴の高次ベクトルVkが生成される。
図13は、このCNNの内部構造を示す。図13に示すように、特徴ベクトルVkの生成に利用されるCNNは、図11に示すものと同一である。即ち、最後の全結合分類層の前の何れかの層からの出力として、Vkが得られる。そして、図12を再度参照すると、ブロック1220において、分類器により特徴ベクトルVkが分類されて、分類ベクトルQkが生成される。この分類器は、事前に訓練用ラベル付き領域組の特徴ベクトルで訓練された2クラスまたは多クラスSVMまたはGBCである。上述のように、本発明の異なる複数の実施形態では、2つまたはそれ以上のクラスで訓練された分類器が利用され、クラス確率判定または分類判定ベクトルとしてQkを出力してもよい。
本発明の別の実施形態において、図1のブロック130は、図14に示す構造を有する。図14に示すように、各領域Rkの品質評価が、異なる複数の技術により異なる基準(輝度、コントラスト、ぼけ等を含むがこれらに限定されない)で行われる。この領域品質分析には、公知の各種画像処理技術が利用されてもよい。ブロック1400の輝度品質測定基準は、例えば
で計算される。式中、fは平均領域強度を示し、τLおよびτHはそれぞれ領域における最低および最高許容強度を示す。ブロック1410のコントラスト品質測定基準は適切な技術により計算され得る。例えば次の文献内のいずれかの技術である:Tripathi,Abhishek Kumar,Sudipta Mukhopadhyay,and Ashis Kumar Dhara.“Performance metrics for image contrast.”Image Information Processing(ICIIP),2011 International Conference on.IEEE,2011,normalized to[0,1]。同様に、ブロック1420のぼけ品質測定基準は適切な技術により計算され得る。例えば次の文献内の技術である:Marziliano,Pina,et al.“A no−reference perceptual blur metric.”Image Processing.2002.Proceedings.2002 International Conference on.Vol.3.IEEE,2002,normalized to[0,1]。
最後に、ブロック1430において、輝度、コントラスト、ぼけ等の測定基準に基づいて、「高品質」測定基準が
で計算され、最終品質評価ベクトルQkが各領域Rkについて計算される。
次に、図1のブロック140の動作について検討する。図1のブロック140では、DR病変分類ベクトルDkおよび品質分類ベクトルQkが分析され、調整済みDR病変分類ベクトルD''kが生成される。この動作は、主に3つの工程を含む。具体的には、まず(1)DR病変分類ベクトル空間確率調整。この工程は任意である。各DR病変分類ベクトルDkは網膜領域Rkに対応し、上述のように各Rk=R(x,y)は、正規化画像I'の画素(x,y)を中心とした網膜領域である。したがって、I'の画素(x,y)を中心とした網膜領域Rkに対応する各DR病変分類ベクトルDkの値を、例えば(x,y)を中心に距離d内の周辺領域の病変分類ベクトルに基づいて、調整可能である。この調整は、例えばフィルタリング、中央値フィルタリング等の形式であってもよく、網膜眼底画像に対してよりロバストで信頼性の高いDR病変分類ベクトルの作成において好ましい。次に(2)品質分類ベクトル確率調整。この工程は任意であり、工程(1)と同じ原理であるが、動作実行対象は品質分類ベクトルQとなる。そして(3)品質分類ベクトルに基づくDR病変分類ベクトル調整である。
次に、図1のブロック140の構造の詳細を検討する。図15は、図1のブロック140の構造を示す。ブロック1500において、DR病変分類ベクトル空間確率調整が実行される。即ち、I'の画素(x,y)を中心とした網膜領域Rkに対応する各DR病変分類ベクトルDkの値が、例えば(x,y)を中心に距離d内の周辺領域のDR病変分類ベクトルに基づいて、調整可能である。この調整は、例えば平均値、中央値、最小値、最大値またはその他フィルタリングであって、調整済みDR病変分類ベクトルD'kが生成される。これは次のように捉えることもできる。即ち、画像I'に対して、l個の空間確率マップが生成され得る。ここで、lはDにおけるDR分類ラベル数である。例えば、l=2個クラスの場合、c=0を「健常」クラスラベルとして、c=1を「病気」クラスラベルとして、「健常」確率マップDR0および「病気」確率マップDR1が生成され得る。そして、I'の画素(x,y)を中心とした網膜領域Rkに対応する各DR病変分類ベクトルDkについて、
となる。式中、sは領域抽出におけるサンプリング単位を示す。そして各DRcが空間処理(例えば、平均値、中央値、最小値、最大値またはその他フィルタ)、或いはその他モルフォロジー演算(膨張、収縮等)を行って、調整済み確率マップDRP'cが生成される。これを、調整済みDR病変分類ベクトルD'kにマッピングできる。ブロック1500と同様に、ブロック1510において品質分類ベクトル空間確率調整が実行される。最後に、図15のブロック1520において、品質分類ベクトルに基づいて、DR病変分類ベクトル調整が実行される。この調整は様々な方法で行われ得る。
一例として、各D'kと、l個のDRクラスの内の各クラスラベルcに対して、この調
整は以下のように実行され得る。
式中、
は領域Rkに対する「高品質」クラス確率を示す。原則として、上記関係において、領域Rkに対する「高品質」クラス確率が0.5未満であれば、領域Rkに対する全てのDR分類ベクトル確率は0に下げられる。即ち、品質が低すぎるので、当該領域に対して、システムはDR分類確率を一切出さない。
別の例として、各D'kと、l個のDRクラスの内の各クラスラベルcに対して、この
調整は以下のように実行され得る。

式中、
は領域Rkに対する「高品質」クラス確率を示す。この関係において、低品質によりDR分類確率を0に下げる切り捨ては、より漸進的に判断される。
上述の例では、各D'kは、対応するQ'kに基づく。発明の別の実施形態では、I'の画素(x,y)を中心とした網膜領域Rkに対応する各D'kは、対応するQ'kと、例えば(x,y)を中心に距離d内の周辺領域の病変分類ベクトルに基づいて、調整可能である。
次に、図1のブロック150の動作を考える。図1のブロック150では、DR病変分類ベクトルD''が分析され、最終段階網膜眼底画像DR分類ベクトルC=(C0,…,Cw-1)が生成される。ここで、wはクラスの数、Cvはクラスvのクラス所属確率を示す。本発明の一実施形態において、w=2で、v=0は「健常」クラスラベルを示し、v=1は「病気」クラスラベルを示す。第2の例として、本発明の別の実施形態において、w=4クラスで、各クラスは、対応するDRの医学的段階に対応する。いずれの場合でも、本発明の別の実施形態において、Cvは二値で0または1の値であってもよく、その場合、Cはクラス確率分布ではなく分類判定ベクトルとなる。
次に、図1のブロック150の構造をより詳細に検討する。本発明の一実施形態において、図1のブロック150は図16の構造を有する。図16のブロック1600では、例えば以下のようにDR病変分類ベクトル確率は閾値処理される。
そして、ブロック1610において、D'''は空間処理されて、D''''が生成される。このブロックは、図15のブロック1500と同様に動作する。即ち、中央値フィルタリング、膨張、収縮等の動作が行われる。本発明の別実施形態では、ブロック1610は省略されてもよいし、さらに異なる実施形態では、ブロック1600、1610の順序を逆にしてもよい。次に、ブロック1620において、D''''内の各クラスラベルlに対する領域数が数えられる。上述のように、本発明のいくつかの実施形態では、l=2クラスで、c=0は「健常」クラスラベルで、c=1は「病気」クラスラベルとなる。本発明の別の実施形態では、l>2クラスで、c=0は「健常」クラスラベルを示し、c=1「微小動脈瘤」クラスラベルを示し、c=2は「硬性白斑」クラスラベルを示し、…、c=l−1は「レーザー傷」クラスラベルを示す。したがって、ブロック1620では、ベクトルG=(G0,…,Gl-1)が生成される。各要素Gcは、対応するラベルcに対する領域数である。そして、領域数ベクトルGに基づいて、ブロック1630で最終網膜眼底画像分類判定ベクトルC=(C0,…,Cw-1)が生成される。ここで、wはクラス数で、Cvはクラスvに対する判定を示す。上述のように、本発明の一実施形態において、w=2で、v=0は「健常」クラスラベルを示し、v=1は「病気」クラスラベルを示す。本発明の別の実施形態において、w=4クラスで、各クラスは、網膜眼底画像においける異なるDR病変の種類の数に基づいて定義可能な、対応するDRの医学的段階に対応する。
本発明の別の実施形態では、図1のブロック150は図17の構造を有する。図17のブロック1700において、DR病変分類ベクトルD''は、l個の2D確率マップに変換される。lはD''におけるDR分類ラベル数を示す。DR分類ベクトルを二次元確率マップに変換する処理は、上述した図15のブロック1500の処理と実質的に同じである。二次元確率マップの作成時、ブロック1700においては、(i)確率のダイナミックレンジの変更(例えば、実際の範囲[0,1]から全範囲[0,255]に)および(i)2D確率マップのt×t画素による固定解像度へのサブサンプリング(例えば、t=56画素)が行われてもよい。したがって、ブロック1700はl個の2D確率マップP=(P0,…,Pl-1)を生成する。次に、ブロック1710において、CNNを利用して確率マップPが分類されて、最終段階網膜眼底画像DR分類ベクトルC=(C0,…,Cw-1)が生成される。ここで、wはクラスの数、Cvはクラスvのクラス所属確率を示す。CNNは事前に大量の訓練用サンプルで訓練されている。上述のように本発明の異なる実施形態では、CNNは異なる数のクラスで訓練されてもよい、例えば、w=2で、v=0は「健常」クラスラベルを示し、v=1は「病気」クラスラベルを示す。或いは、w=4クラスで、各クラスは、対応するDRの医学的段階に対応する。いずれの場合でも、本発明の別の実施形態において、Cvは0か1かの二値であり、Cは所属確率分布ではなく、分類判定ベクトルとなる。図18は、ブロック1710のCNNの内部構造を示す。これは上述した図6および図11のCNN構造と実質的に同じである。ただし、サイズ、幅、フィルタ数、学習済み畳み込みマスク、畳み込み層数、全結合層数等のネットワークパラメータが異なり得る。
本発明の別の実施形態では、図1のブロック150は図19に示す構造を有する。ブロック1900は、図17のブロック1700と同様に動作する。そして、ブロック1910において、CNNを利用して、際立った特徴の高次ベクトルVが生成される。図20はこのCNNの内部構造を示す。図20に示すように、特徴ベクトルV、の生成に利用されるCNNは、図18に示すものと同一である。即ち、最後の全結合分類層の前の何れかの層からの出力として、Vが得られる。そして、図19を再度参照すると、ブロック1920において、分類器により特徴ベクトルVが分類されて、分類ベクトルCが生成される。この分類器は、事前に訓練用サンプル組の特徴ベクトルで訓練された二値または多クラスSVMまたは勾配ブースティング分類器である。上述のように、本発明の異なる複数の実施形態では、2つまたはそれ以上のクラスで訓練された分類器が利用され、クラス確率判定または分類判定ベクトルとしてCを出力してもよい。
図21は、本発明の別の実施形態を示す。ここで、ブロック2100と2110は図1のブロック100と110と実質的に同じように動作し、画像正規化と領域抽出が実行される。
次に、図21のブロック2120の動作を考える。図21のブロック2120においては、全ての網膜領域は個別に分析され、各領域Rkについて、DR病変/品質分類合成ベクトル
が生成される。上述のように、lはDR病変クラスの数で、mは領域品質クラスの数である。即ち、ブロック2120における分類器の訓練において、各訓練用サンプルはDR病変クラスラベルと品質クラスラベルの両方が割り当てられ、分類器が領域RkについてDR病変/品質分類連結ベクトル
を生成可能になる。上述のように、ブロック2120の分類器は、クラス所属確率ではなく二値(0または1)を出力するように構成され得る
次に、図21のブロック2120の構造の詳細を検討する。図21のブロック2120は図22に示す構造を有する。ブロック2200では、図5のブロック500または図7のブロック700と実質的に同じように、網膜領域正規化が実行される。即ち、各網膜領域Rkについて、ブロック2200では正規化領域R'kが生成される。次に、ブロック2210において、CNNを利用して、各領域R'kを分類して、分類ベクトルAkが生成される。図23は、図22のブロック2210の構造をより詳細に示す。これは、上述の図6および図11のCNNの構造と実質的に同じである。ただし、サイズ、移動幅、フィルタ数、学習済み畳み込みマスク、畳み込み層数、全結合層数等のネットワークパラメータが異なり得る。
次に、図21に戻り、ブロック2130において、図1のブロック140と実質的に同じ用にして、領域DR病変確率調整が行われる。そして、ブロック2140において、図1のブロック150と実質的に同じようにして、網膜眼底画像分類が実行される。
いくつかの実施形態では、撮像処理はコンピューティングシステムによって実行される。いくつかの実施形態では、コンピューティングシステムは、1つまたは複数のコンピューティングデバイス、例えば、IBM、マッキントッシュ、マイクロソフトウィンドウズ(登録商標)、またはリナックス(登録商標)/ユニックスに対応するパーソナルコンピュータ、またはサーバまたはワークステーションを含む。一実施形態では、コンピューティングデバイスは、例えば、サーバ、ノートパソコン、スマートフォン、携帯情報端末、キオスク、またはメディアプレーヤーを含む。一実施形態では、コンピューティングデバイスは、1つまたは複数のCPUを含み、各CPUは、従来のまたは独自のマイクロプロセッサを含んでもよい。コンピューティングデバイスは、情報の一時記憶用のランダムアクセスメモリ(「RAM」)、情報の永久記憶用の1つ以上の読み取り専用メモリ(「ROM」)、ハードドライブ、ディスケット、ソリッドステートドライブ、および光学式メディアストレージデバイス等の1つ以上の大容量記憶装置等の1つ以上のメモリをさらに含む。通常、コンピューティングデバイスのモジュールは、一般的なバスシステムを使用してコンピュータに接続される。異なる様々な実施形態では、一般的なバスシステムは、例えば、周辺機器相互接続(PCI)、マイクロチャネル、小型コンピュータシステムインタフェース(SCSI)、産業標準アーキテクチャ(ISA)および拡張ISA(EISA)アーキテクチャで実現できる。加えて、コンピューティングデバイスのコンポーネントおよびモジュールに対して設けられた機能は、組み合わされてより少ないコンポーネントおよびモジュールとされ得るか、さらに分割されて追加のコンポーネントおよびモジュールに設けられてもよい。
コンピューティングデバイスは通常、Windows XP、Windows Vista、Windows 7、Windows 8、Windows Server、Embedded Windows、Unix、Linux(登録商標)、Ubuntu Linux、SunOS、Solaris、iOS、Blackberry OS、Android等のオペレーティングシステムソフトウェア、または他の対応するオペレーティングシステムによって制御および調整される。マッキントッシュシステムでは、オペレーティングシステムは、MAC OS X等の利用可能な任意のオペレーティングシステムであり得る。他の実施形態では、コンピューティングデバイスは独自のオペレーティングシステムにより制御され得る。従来のオペレーティングシステムは、実行のためにコンピュータプロセスを制御およびスケジューリングし、メモリを管理し、ファイルシステム、ネットワーク、I/Oサービスを提供し、特にグラフィカルユーザインタフェース(GUI)等のユーザインタフェースを提供する。
例示的なコンピューティングデバイスは、キーボード、マウス、タッチパッド、タッチスクリーン、およびプリンタ等の1つまたは複数の一般的に利用可なI/Oインタフェースおよびデバイスを備えてもよい。一実施形態では、I/Oインタフェースおよびデバイスは、ユーザへのデータの視覚的提示を可能にするモニターまたはタッチスクリーンモニター等の1つ以上のディスプレイデバイスを備える。より具体的には、ディスプレイデバイスは、例えば、GUI、アプリケーションソフトウェアデータ、およびマルチメディアプレゼンテーションを提示する。コンピューティングデバイスはまた、例えば、カメラ、スピーカ、ビデオカード、グラフィックスアクセラレータ、およびマイク等の1つまたは複数のマルチメディアデバイスを含んでもよい。I/Oインタフェースとデバイスは、さまざまな外部デバイスへの通信インタフェースを実現する。コンピューティングデバイスは、例えば有線、無線、または有線と無線の組み合わせの通信リンクを介して、LAN、WAN、および/またはインターネットの1つまたは複数を含むネットワークに電子的に結合される。ネットワークは、有線または無線の通信リンクを介して、さまざまなコンピューティングデバイスおよび/または他の電子デバイスと通信する。
本明細書で説明される方法およびシステムにより処理される画像は、1つ以上のデータソースからネットワークを介してコンピューティングシステムに提供され得る。データソースは、1つ以上の内部および/または外部データベース、データソース、物理データストアを含む。データソースは、上記のシステムおよび方法に応じて、撮像システムで処理されるデータを格納するデータベースを含んでもよい。或いは、データソースは、上記のシステムおよび方法により、撮像システムで処理されたデータを格納するデータベースを含んでもよい。いくつかの実施形態では、データベース、データソースの一方または両方は、Sybase、Oracle、CodeBase、MySQL、およびMicrosoft(登録商標)SQLサーバ等のリレーショナルデータベース、ならびに以下のような他の種類のデータベースを使用して実現され得る。即ち、他の種類のデータベースの例は、フラットファイルデータベース、エンティティ関係データベース、リレーショナルデータベース、オブジェクト指向データベース、NoSQLデータベース、レコードベースのデータベース等を含む。
コンピューティングシステムは、CPUによって実行される実行可能なソフトウェアコードとして大容量記憶装置に格納され得るイメージングシステムモジュールを備える。当該モジュールの例としては、ソフトウェアコンポーネント、オブジェクト指向ソフトウェアコンポーネント、クラスコンポーネントおよびタスクコンポーネント等のコンポーネント、プロセス、機能、属性、手順、サブルーチン、プログラムコードのセグメント、ドライバ、ファームウェア、マイクロコード、回路、データ、データベース、データ構造、テーブル、配列および変数が挙げられる。コンピューティングシステムは、例えば、自動化された低レベル画像処理、自動化された画像レジストレーション、自動化された画像評価、自動化されたスクリーニングを実行するため、さらに/或いは上記の新しいアーキテクチャを実装するために、イメージングシステムモジュールを実行するように構成される。
概して、本明細書で使用される「モジュール」という用語は、ハードウェアまたはファームウェアで実現されるロジック、または、例えばPython、Java(登録商標)、Lua、Cおよび/またはC++等のプログラミング言語で記述された、入口点と出口点を有し得るソフトウェア命令の集合を示す。ソフトウェアモジュールは、光学記憶媒体、フラッシュドライブ、またはその他の有形媒体等のコンピュータ可読媒体で実現されてもよい。このようなソフトウェアコードは、コンピューティングデバイスが実行されるように、コンピューティングシステム等の実行するコンピューティングデバイスのメモリデバイスに一部または全部保存され得る。前述のプロセス、方法、およびアルゴリズムのそれぞれは、1つまたは複数のコンピュータシステムまたはコンピュータハードウェアを含むコンピュータプロセッサによって実行されるコードモジュールで実現され、完全または部分的に自動化され得る。コードモジュールは、ハードドライブ、ソリッドステートメモリ、光ディスク等の任意の種類の非一時的なコンピュータ可読媒体またはコンピュータストレージデバイスに保存可能である。システムおよびモジュールは、無線式および有線/ケーブル式を含むさまざまなコンピュータ可読送媒体で、生成されたデータ信号として(例えば、搬送波または他のアナログまたはデジタル伝搬信号の一部として)送信可能なさまざまな形式をとり得る(例えば、単一または多重化されたアナログ信号の一部、または複数の個別のデジタルパケットまたはフレーム)。プロセスおよびアルゴリズムは、一部または全部、特定用途向け回路で実現できる。開示されたプロセスおよびプロセスステップの結果は、例えば揮発性または不揮発性ストレージ等の任意の種類の非一時的なコンピュータストレージに永続的にまたは別の方法で格納されてもよい。
上述のさまざまな特徴およびプロセスは、互いに独立して使用することができる。或いは、各種組み合わせ可能である。すべての可能な組み合わせおよびサブコンビネーションは、本開示の範囲内に含まれることが意図される。さらに、特定の方法またはプロセスブロックは、形態によっては省略され得る。本明細書で説明する方法およびプロセスも特定の順序に限定されず、それに関連するブロックまたは状態は、適切な他の順序で実現することができる。例えば、説明されたブロックまたは状態は、具体的に開示された順序以外の順序で実現されてもよく、または複数のブロックまたは状態が単一のブロックまたは状態に組み合わされてもよい。例示的なブロックまたは状態は、順次、並行して、または他の何らかの方法で実現されてもよい。開示された例示的な実施形態にブロックまたは状態を追加または削除することができる。本明細書で説明する例示的なシステムおよびコンポーネントは、上記とは異なる構成である。例えば、開示された例示的な実施形態に対して、要素を追加、削除、または再配置することができる。
「できる」、「可能である」、「得る」、「可能性がある」等の条件表現は概して、特に明記しない限り、または使用される文脈内で異なる理解がなされる場合を除いて、特定の実施形態が、別の実施形態に含まれない特定の特徴、要素、および/またはステップを含むことを意図している。したがって、そのような条件表現は概して、機能、要素および/またはステップが1つまたは複数の実施形態に何らかの形で必要であること、または1つまたは複数の実施形態が(ユーザ入力またはプロンプトの有無にかかわらず)これらの機能、要素、および/またはステップが特定の実施形態に含まれるまたはそれにより実行されるか否かを決定するためのロジックを必ず含むことを意味するものではない。また、「含む」という用語は、「含まれるがそれに限定されない」ことを意味し、「または」という用語は「および/または」を意味する。
本書本明細書に記載の、および/または添付の図に示されているフローまたはブロック図のプロセス記述、要素、またはブロックは、モジュール、セグメント、または特定の論理機能を実装するための1つ以上の実行可能命令を含むコードの一部またはプロセスにおける工程を表し得るものとして理解されたい。その他形態は、本明細書に記載の実施形態の範囲内に含まれる。そして当業者により要素または機能は、関連する機能に応じて削除でき、図示または説明されたのと異なる順序で実行され得る。これは実質的に同時または逆の順序で行われることを含む。
上記のすべての方法およびプロセスは、1つまたは複数の汎用コンピュータで実行されるソフトウェアコードモジュールで実施され、部分的または完全に自動化される。例えば、本明細書で説明される方法は、コンピューティングシステムおよび/または他の任意の適切なコンピューティングデバイスによって実行され得る。これらの方法は、有形のコンピュータ可読媒体から読み取られたソフトウェア命令または他の実行可能コードの実行に応じて、コンピューティングデバイス上で実行され得る。有形のコンピュータ可読媒体は、コンピュータシステムによって読み取り可能なデータを保存できるデータストレージデバイスである。コンピュータ可読媒体の例としては、リードオンリーメモリ、ランダムアクセスメモリ、その他の揮発性または不揮発性メモリデバイス、CD−ROM、磁気テープ、フラッシュドライブ、および光学データストレージデバイスが挙げられる。
上記の実施形態に対して多くの変形および修正を行うことができることが強調され、その要素は他の許容可能な例の中にあると理解されるべきである。そのような修正および変形はすべて、本開示の範囲内に含まれることが意図される。前述の説明は、特定の実施形態を詳述しているが、前述の内容が文面として如何に詳細に記載されていても、システムおよび方法は多様な実施が可能であることが理解されよう。例えば、一実施形態の特徴は、異なる実施形態の特徴と共に使用されてもよい。また、前述のように、システムおよび方法の特定の機能または態様を説明する際の特定の用語の使用は、その用語が関連するシステムおよび方法の特徴または態様の何らかの特定の特徴を含むように限定されるように、本明細書で再定義されていることを意味するものではないことに留意されたい。

Claims (46)

  1. 眼画像を分類する方法であって、
    眼画像を受信することと、
    前記眼画像を正規化することと、
    前記正規化眼画像を複数の領域に分割することと、
    ニューラルネットワークを有する少なくとも1つの分類器により、前記複数の領域について自動で眼病ベクトルを判定することと、
    前記判定された眼病ベクトルに基づいて、自動で前記複数の領域それぞれをアノテーションすることと、
    少なくとも前記アノテーションに基づいて、少なくとも2つのレベルに対して、前記受信した眼画像を自動で評価することと、を含み
    前記ニューラルネットワークは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が眼画像の各領域につけたアノテーションに基づいて訓練される、方法。
  2. 前記ニューラルネットワークは、複数の隠れ層を含む、請求項1に記載の方法。
  3. 前記ニューラルネットワークは、ディープニューラルネットワークを有する、請求項2に記載の方法。
  4. 前記ニューラルネットワークは、畳み込みニューラルネットワークを有する、請求項1に記載の方法。
  5. 前記ニューラルネットワークは再帰型ニューラルネットワークを含み、
    前記血管画像とは異なるタイミングで第2の血管画像を受信することと、
    少なくとも前記眼画像の経時変化に応じた情報を出力することと、をさらに含む、請求項1に記載の方法。
  6. 網膜画像の各領域の品質について、少なくとも専門家によるアノテーションを受信することをさらに含む、請求項1に記載の方法。
  7. 前記少なくとも1つの分類器は、多クラスサポートベクターマシン分類器を含む、請求項1に記載の方法。
  8. 前記少なくとも1つの分類器は、勾配ブースティング分類器を含む、請求項1に記載の方法。
  9. 前記少なくとも1つの分類器は、複数の各領域の前記眼病ベクトルに基づいて、各領域を分類する、請求項1に記載の方法。
  10. 前記少なくとも1つの分類器は、複数の連続した領域の前記眼病ベクトルに基づいて、各領域を分類する、請求項1に記載の方法。
  11. 前記少なくとも1つの分類器は、糖尿病網膜症の兆候を持つ領域を判定する分類器を含む、請求項1に記載の方法。
  12. 前記評価することは、糖尿病網膜症の程度を評価することを含む、請求項1に記載の方法。
  13. 前記レベルには、当該レベルが正確である可能性を伴う、請求項1に記載の方法。
  14. 前記正確である可能性は、第一種過誤の可能性と、第二種過誤の可能性を含む、請求項13に記載の方法。
  15. 前記眼画像の前記領域の分類を示す画像を出力することをさらに含む、請求項1に記載の方法。
  16. 前記眼画像は眼底画像を含み、前記アノテーションは少なくとも、各眼底画像の各領域内の糖尿病網膜症の兆候を含む、請求項1に記載の方法。
  17. 前記自動評価は、少なくとも3つのレベルに対して行われる、請求項1に記載の方法。
  18. 前記評価することは、少なくとも5つのレベルを含む、一般的な糖尿病網膜症評価基準に応じて、糖尿病網膜症の程度を評価することを含む、請求項1に記載の方法。
  19. 前記ニューラルネットワークは、網膜画像の各領域の品質に対する専門家によるアノテーションに少なくとも基づいてさらに訓練される、請求項1に記載の方法。
  20. 前記ニューラルネットワークは、前記網膜画像の各領域を、少なくともその品質およびレベルに応じて重みづけする、請求項19に記載の方法。
  21. 眼画像を分類するシステムであって、
    少なくとも1つの眼画像を受信するように構成された入力と、
    少なくとも1つの客観的分類基準に応じて、複数の網膜画像の各領域の専門家によるアノテーションに少なくとも基づいて訓練されるニューラルネットワークを定義する情報を記憶するように構成されるメモリと、
    少なくとも1つの自動プロセッサであって、
    前記眼画像を正規化することと、
    前記正規化眼画像を複数の領域に分割することと、
    前記定義されたニューラルネットワークを有する少なくとも1つの分類器により、前記複数の領域について眼病ベクトルを判定することと、
    前記判定された眼病ベクトルに基づいて、前記複数の領域それぞれをアノテーションすることと、
    少なくとも前記アノテーションに基づいて、少なくとも2つのレベルに対して、前記受信した眼画像を評価することとを実行するように構成された自動プロセッサと、
    前記評価を通信するように構成された出力とを備えるシステム。
  22. 前記ニューラルネットワークは、複数の隠れ層を含む、請求項21に記載のシステム。
  23. 前記ニューラルネットワークは、ディープニューラルネットワークを有する、請求項22に記載のシステム。
  24. 前記ニューラルネットワークは、畳み込みニューラルネットワークを有する、請求項22に記載のシステム。
  25. 前記ニューラルネットワークは再帰型ニューラルネットワークを含み、
    前記入力は前記眼画像とは異なるタイミングで第2の眼画像を受信するようにさらに構成され、
    前記少なくとも1つの自動プロセッサは、少なくとも前記眼画像の経時変化を分析するように構成される、請求項21に記載のシステム。
  26. 前記ニューラルネットワークは、網膜画像の各領域の品質について、少なくとも専門家によるアノテーションに基づいてさらに訓練される、請求項21に記載のシステム。
  27. 前記メモリは、網膜画像の各領域の品質について、少なくとも専門家によるアノテーションをさらに記憶するように構成される、請求項21に記載のシステム。
  28. 前記少なくとも1つの分類器は、多クラスサポートベクターマシン分類器を含む、請求項21に記載のシステム。
  29. 前記少なくとも1つの分類器は、勾配ブースティング分類器を含む、請求項21に記載のシステム。
  30. 前記少なくとも1つの分類器は、複数の各領域の前記眼病ベクトルに基づいて、各領域を分類する、請求項21に記載のシステム。
  31. 前記少なくとも1つの分類器は、複数の連続した領域の前記眼病ベクトルに基づいて、各領域を分類する、請求項21に記載のシステム。
  32. 前記少なくとも1つの分類器は、糖尿病網膜症の兆候を持つ領域を判定する分類器を含む、請求項21に記載のシステム。
  33. 前記評価することは、糖尿病網膜症の程度を評価することを含む、請求項21に記載のシステム。
  34. 前記レベルには、当該レベルが正確である可能性を伴う、請求項21に記載のシステム。
  35. 前記正確である可能性は、第一種過誤の可能性と、第二種過誤の可能性を含む、請求項21に記載のシステム。
  36. 前記出力は、前記眼画像の前記領域の分類を示す画像を出力するようにさらに構成される、請求項21に記載のシステム。
  37. 前記眼画像は眼底画像を含み、前記アノテーションは、少なくとも、各眼底画像の各領域内の糖尿病網膜症の兆候を含む、請求項21に記載のシステム。
  38. 前記眼画像は血管画像を含む、請求項21に記載のシステム。
  39. 前記眼病は血管障害を含む、請求項38に記載のシステム。
  40. 前記受信した画像は、少なくとも前記アノテーションに基づき、少なくとも3つのレベルに対して評価される、請求項21に記載の方法。
  41. 前記受信した眼画像に対して、少なくとも5つのレベルを含む、一般的な糖尿病網膜症評価基準に応じて、糖尿病網膜症の程度が評価される、請求項21に記載のシステム。
  42. 前記ニューラルネットワークは、網膜画像の各領域の品質に対して専門家によるアノテーションに少なくとも基づいてさらに訓練される、請求項21に記載のシステム。
  43. 前記ニューラルネットワークは、前記網膜画像の各領域を、少なくともその品質およびレベルに応じて重みづけするように構成される、請求項42に記載のシステム。
  44. 少なくとも1つの自動プロセッサを制御するための指示を記憶するコンピュータ可読媒体であって、
    眼画像を受信するための指示と、
    前記眼画像を正規化するための指示と、
    前記正規化眼画像を複数の領域に分割するための指示と、
    ニューラルネットワークを有する少なくとも1つの分類器により、前記複数の領域について眼病ベクトルを判定するための指示と、
    前記判定された眼病ベクトルに基づいて、自動で前記複数の領域それぞれをアノテーションするための指示と、
    少なくとも前記アノテーションに基づいて、少なくとも2つのレベルに対して、前記受信した眼画像を自動で評価するための指示と、を含み
    前記ニューラルネットワークは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が眼画像の各領域につけたアノテーションに基づいて訓練される、媒体。
  45. 網膜症の存在について、網膜画像を分類する方法であって、
    眼底画像を受信することと、
    少なくとも三色刺激画像空間に対して眼底画像を自動で正規化することと、
    前記正規化眼画像を、少なくとも1つの分割ルールに基づいて、複数の領域に自動で分割することと、
    少なくとも専門家による網膜アノテーションを有する眼底画像組による訓練に基づいて訓練されたニューラルネットワークを有する少なくとも1つの分類器により、前記複数の領域について網膜症ベクトルを自動で判定することと、
    前記判定した眼病ベクトルに基づいて、自動で前記複数の領域それぞれにアノテーションを行うことと、
    少なくとも前記アノテーションに基づいて、少なくとも2つの網膜症レベルに対して、前記受信した眼画像の網膜症の程度を自動で評価することと、
    前記自動で判定した網膜症の程度を出力することと、を含む方法。
  46. 病理学画像を分類する方法であって、
    画像を受信することと、
    前記画像を正規化することと、
    前記正規化画像を複数の領域に分割することと、
    ニューラルネットワークを有する少なくとも1つの分類器により、前記複数の領域について自動で疾患ベクトルを判定することと、
    前記判定された疾患ベクトルに基づいて、自動で前記複数の領域それぞれをアノテーションすることと、
    少なくとも前記アノテーションに基づいて、前記受信した画像を自動で評価することと、を含み
    前記ニューラルネットワークは、少なくとも1つの客観的な分類基準に応じて、少なくとも専門家が画像の各領域につけたアノテーションに基づいて訓練される、方法。
JP2019559095A 2017-04-27 2018-04-26 自動眼底画像分析用のシステムおよび方法 Pending JP2020518915A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201762490977P 2017-04-27 2017-04-27
US62/490,977 2017-04-27
PCT/US2018/029605 WO2018200840A1 (en) 2017-04-27 2018-04-26 System and method for automated funduscopic image analysis

Publications (1)

Publication Number Publication Date
JP2020518915A true JP2020518915A (ja) 2020-06-25

Family

ID=62186545

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019559095A Pending JP2020518915A (ja) 2017-04-27 2018-04-26 自動眼底画像分析用のシステムおよび方法

Country Status (6)

Country Link
US (2) US10719936B2 (ja)
EP (1) EP3616120B1 (ja)
JP (1) JP2020518915A (ja)
CN (1) CN111095261A (ja)
BR (1) BR112019022447A2 (ja)
WO (1) WO2018200840A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2023539483A (ja) * 2020-11-02 2023-09-14 ▲騰▼▲訊▼科技(深▲セン▼)有限公司 医用画像処理方法、装置、機器、記憶媒体及びコンピュータプログラム
WO2023228229A1 (ja) * 2022-05-23 2023-11-30 日本電気株式会社 画像処理装置、画像処理方法、およびプログラム
US12100140B2 (en) 2021-01-29 2024-09-24 Fujifilm Corporation Information processing apparatus, information processing method, and information processing program
JP7575216B2 (ja) 2020-08-25 2024-10-29 株式会社トプコン 眼科画像処理装置、眼科装置、画質評価方法及びプログラム

Families Citing this family (92)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10325176B2 (en) * 2014-08-25 2019-06-18 Agency For Science, Technology And Research Methods and systems for assessing retinal images, and obtaining information from retinal images
CN108603922A (zh) * 2015-11-29 2018-09-28 阿特瑞斯公司 自动心脏体积分割
RU2018142757A (ru) * 2016-05-13 2020-06-15 Басф Се Система и способ обнаружения болезней растений
US11195313B2 (en) * 2016-10-14 2021-12-07 International Business Machines Corporation Cross-modality neural network transform for semi-automatic medical image annotation
US10663711B2 (en) 2017-01-04 2020-05-26 Corista, LLC Virtual slide stage (VSS) method for viewing whole slide images
WO2018140596A2 (en) 2017-01-27 2018-08-02 Arterys Inc. Automated segmentation utilizing fully convolutional networks
CN108771530B (zh) * 2017-05-04 2021-03-30 深圳硅基智能科技有限公司 基于深度神经网络的眼底病变筛查系统
US11748877B2 (en) 2017-05-11 2023-09-05 The Research Foundation For The State University Of New York System and method associated with predicting segmentation quality of objects in analysis of copious image data
WO2019005722A1 (en) * 2017-06-26 2019-01-03 The Research Foundation For The State University Of New York SYSTEM, METHOD AND COMPUTER-ACCESSIBLE MEDIA FOR VIRTUAL PANCREATOGRAPHY
WO2019024568A1 (zh) * 2017-08-02 2019-02-07 上海市第六人民医院 眼底图像处理方法、装置、计算机设备和存储介质
KR101977645B1 (ko) * 2017-08-25 2019-06-12 주식회사 메디웨일 안구영상 분석방법
US11037032B2 (en) * 2017-10-06 2021-06-15 Wisconsin Alumni Research Foundation Methods, systems, and media for detecting the presence of an analyte
CN107680684B (zh) * 2017-10-12 2021-05-07 百度在线网络技术(北京)有限公司 用于获取信息的方法及装置
US10796221B2 (en) 2017-10-19 2020-10-06 General Electric Company Deep learning architecture for automated image feature extraction
US10460440B2 (en) * 2017-10-24 2019-10-29 General Electric Company Deep convolutional neural network with self-transfer learning
US10977554B2 (en) * 2017-11-14 2021-04-13 Magic Leap, Inc. Fully convolutional interest point detection and description via homographic adaptation
WO2019103912A2 (en) 2017-11-22 2019-05-31 Arterys Inc. Content based image retrieval for lesion analysis
WO2019108888A1 (en) * 2017-11-30 2019-06-06 The Research Foundation For The State University Of New York SYSTEM AND METHOD TO QUANTIFY TUMOR-INFILTRATING LYMPHOCYTES (TILs) FOR CLINICAL PATHOLOGY ANALYSIS
CN108615051B (zh) * 2018-04-13 2020-09-15 博众精工科技股份有限公司 基于深度学习的糖尿病视网膜图像分类方法及系统
US11461599B2 (en) * 2018-05-07 2022-10-04 Kennesaw State University Research And Service Foundation, Inc. Classification of images based on convolution neural networks
US10902586B2 (en) * 2018-05-08 2021-01-26 International Business Machines Corporation Automated visual recognition of a microcalcification
AU2019275232B2 (en) * 2018-05-21 2024-08-15 Corista, LLC Multi-sample whole slide image processing via multi-resolution registration
CN110648303B (zh) * 2018-06-08 2022-07-26 上海市第六人民医院 眼底图像分析方法、计算机设备和存储介质
KR102354980B1 (ko) * 2018-06-11 2022-01-24 사회복지법인 삼성생명공익재단 전안부 질환 진단 시스템 및 이를 이용한 진단 방법
CN109686423A (zh) * 2018-11-06 2019-04-26 众安信息技术服务有限公司 一种医疗影像标注方法及系统
CN109523524B (zh) * 2018-11-07 2020-07-03 电子科技大学 一种基于集成学习的眼底图像硬性渗出检测方法
WO2020107167A1 (zh) * 2018-11-26 2020-06-04 深圳先进技术研究院 乳腺密度自动分级方法及装置
CN109727256B (zh) * 2018-12-10 2020-10-27 浙江大学 一种基于玻尔兹曼和目标先验知识的图像分割识别方法
AU2019396539A1 (en) * 2018-12-12 2021-08-05 Tesseract Health, Inc. Biometric identification and health status determination
CN109602391A (zh) * 2019-01-04 2019-04-12 平安科技(深圳)有限公司 眼底出血点的自动检测方法、装置及计算机可读存储介质
CN110021009B (zh) * 2019-01-18 2023-07-21 平安科技(深圳)有限公司 一种评估眼底图像质量的方法、装置及存储介质
CN109902717A (zh) * 2019-01-23 2019-06-18 平安科技(深圳)有限公司 病灶自动识别方法、装置及计算机可读存储介质
US11068694B2 (en) * 2019-01-23 2021-07-20 Molecular Devices, Llc Image analysis system and method of using the image analysis system
CN109872322A (zh) * 2019-02-27 2019-06-11 电子科技大学 一种基于级联检测模型的核性白内障病变区域定位方法
US20220151568A1 (en) * 2019-03-13 2022-05-19 The Board Of Trustees Of The University Of Illinois Supervised machine learning based multi-task artificial intelligence classification of retinopathies
US11583188B2 (en) * 2019-03-18 2023-02-21 General Electric Company Automated detection and localization of bleeding
WO2020198380A1 (en) 2019-03-26 2020-10-01 Tempus Labs, Inc. Determining biomarkers from histopathology slide images
CN110060258B (zh) * 2019-03-27 2021-05-04 山东师范大学 基于高斯混合模型聚类的视网膜sd-oct图像分割方法和装置
US11074479B2 (en) * 2019-03-28 2021-07-27 International Business Machines Corporation Learning of detection model using loss function
EP3948773A1 (en) * 2019-03-29 2022-02-09 Ai Technologies Inc. Image-based detection of ophthalmic and systemic diseases
WO2020201942A1 (en) 2019-04-01 2020-10-08 Cathworks Ltd. Methods and apparatus for angiographic image selection
CN109977905B (zh) * 2019-04-04 2021-08-06 北京百度网讯科技有限公司 用于处理眼底图像的方法和装置
CN110136140A (zh) * 2019-04-16 2019-08-16 上海鹰瞳医疗科技有限公司 眼底图像血管影像分割方法及设备
US12050830B2 (en) * 2019-04-25 2024-07-30 Dolby Laboratories Licensing Corporation Content-aware PQ range analyzer and tone mapping in live feeds
CN110211094B (zh) * 2019-05-06 2023-05-26 平安科技(深圳)有限公司 黑眼圈智能判断方法、装置及计算机可读存储介质
WO2020243333A1 (en) 2019-05-30 2020-12-03 The Research Foundation For The State University Of New York System, method, and computer-accessible medium for generating multi-class models from single-class datasets
WO2020241794A1 (ja) * 2019-05-31 2020-12-03 株式会社ニデック 眼科画像処理装置、眼科画像処理プログラム、および眼科画像処理システム
US11915826B2 (en) * 2019-06-07 2024-02-27 Welch Allyn, Inc. Digital image screening and/or diagnosis using artificial intelligence
CA3145370A1 (en) * 2019-06-25 2020-12-30 Owkin Inc. Systems and methods for mesothelioma feature detection and enhanced prognosis or response to treatment
TWI702615B (zh) * 2019-07-26 2020-08-21 長佳智能股份有限公司 視網膜病變評估模型建立方法及系統
WO2021041068A1 (en) * 2019-08-30 2021-03-04 Google Llc Predicting disease progression from tissue images and tissue segmentation maps
KR20210035381A (ko) * 2019-09-23 2021-04-01 삼성에스디에스 주식회사 의료 진단 방법 및 장치
WO2021059165A1 (en) 2019-09-23 2021-04-01 Cathworks Ltd. Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device
KR20220077142A (ko) * 2019-10-02 2022-06-08 메모리얼 슬로안 케터링 캔서 센터 다중-클래스 이미지 분할을 위한 심층 다중-배율 네트워크들
WO2021080804A1 (en) * 2019-10-25 2021-04-29 F. Hoffmann-La Roche Ag Machine-learning techniques for prediction of future visual acuity
CN110910371A (zh) * 2019-11-22 2020-03-24 北京理工大学 基于生理指标和影像融合的肝肿瘤自动分类方法及装置
CN110992364B (zh) * 2019-12-31 2023-11-28 重庆艾可立安医疗器械有限公司 视网膜图像识别方法、装置、计算机设备和存储介质
US11734576B2 (en) * 2020-04-14 2023-08-22 International Business Machines Corporation Cooperative neural networks with spatial containment constraints
US11222201B2 (en) 2020-04-14 2022-01-11 International Business Machines Corporation Vision-based cell structure recognition using hierarchical neural networks
US11797603B2 (en) 2020-05-01 2023-10-24 Magic Leap, Inc. Image descriptor network with imposed hierarchical normalization
US11875480B1 (en) * 2020-05-15 2024-01-16 Verily Life Sciences Llc Distinguishing artifacts from pathological features in digital images
CN113706375A (zh) * 2020-05-20 2021-11-26 佛吉亚歌乐电子(厦门)有限公司 一种图像拼接装置及方法
CN111739617B (zh) * 2020-06-15 2023-07-11 中国医学科学院医学信息研究所 医疗影像人工智能质控标注方法、装置、设备及存储介质
CN111754486B (zh) * 2020-06-24 2023-08-15 北京百度网讯科技有限公司 图像处理方法、装置、电子设备及存储介质
CN115997241A (zh) 2020-08-13 2023-04-21 佩治人工智能公司 处理电子图像进行连续生物标志物预测的系统和方法
CN112001848B (zh) * 2020-09-07 2022-04-26 鹏祥智慧保安有限公司 大数据监控系统中的图像识别拼接方法及系统
US12079987B2 (en) * 2020-09-19 2024-09-03 The Cleveland Clinic Foundation Automated quality assessment of ultra-widefield angiography images
CN114254146A (zh) * 2020-09-21 2022-03-29 京东方科技集团股份有限公司 图像数据的分类方法、装置和系统
CN112184697B (zh) * 2020-10-15 2022-10-04 桂林电子科技大学 基于果蝇优化算法的糖尿病视网膜病变分级深度学习方法
CN112465839A (zh) * 2020-12-10 2021-03-09 山东承势电子科技有限公司 一种基于数据增强的眼底图像病灶分割和定量分析方法
CN112651938B (zh) * 2020-12-24 2023-12-19 平安科技(深圳)有限公司 视盘图像分类模型的训练方法、装置、设备及存储介质
CN112734697B (zh) * 2020-12-24 2024-05-28 沈阳东软智能医疗科技研究院有限公司 图像处理方法、装置、存储介质和电子设备
US12079311B2 (en) * 2021-01-08 2024-09-03 Salesforce, Inc. AI-enhanced data labeling
CN113128565B (zh) * 2021-03-25 2022-05-06 之江实验室 面向预训练标注数据不可知的图像自动标注系统和装置
CN113205523A (zh) * 2021-04-29 2021-08-03 浙江大学 多尺度表征优化的医学影像分割识别系统、终端及存储介质
US12067726B2 (en) * 2021-05-06 2024-08-20 Vestel Elektronik Sanayi Ve Ticaret A.S. Retina image annotation, and related training methods and image processing models
CN113299388B (zh) * 2021-05-12 2023-09-29 吾征智能技术(北京)有限公司 基于发热伴皮疹的跨模态医学生物特征认知疾病的系统
CN113344909B (zh) * 2021-07-01 2023-12-08 中国石油大学(北京) 一种火电锅炉透火焰高温滤镜结焦识别显示方法和装置
CN113768460B (zh) * 2021-09-10 2023-11-14 北京鹰瞳科技发展股份有限公司 一种眼底图像分析系统、方法以及电子设备
CN113850804B (zh) * 2021-11-29 2022-03-18 北京鹰瞳科技发展股份有限公司 基于生成对抗网络的视网膜图像生成系统及方法
CN114511727A (zh) * 2021-12-14 2022-05-17 广东省建筑科学研究院集团股份有限公司 基于多尺度输入神经网络的建筑风场地面粗糙度识别方法
US20230252632A1 (en) * 2022-02-10 2023-08-10 Cathworks Ltd. System and method for machine-learning based sensor analysis and vascular tree segmentation
US11877831B2 (en) * 2022-03-14 2024-01-23 O/D Vision Inc. Systems and methods for artificial intelligence based blood pressure computation based on images of the outer eye
CN114943845B (zh) * 2022-05-23 2024-08-02 天津城建大学 一种领域图片细粒度分类识别方法及系统
CN115082414B (zh) * 2022-07-08 2023-01-06 深圳市眼科医院 一种基于视觉质量分析的便携式检测方法和装置
WO2024011236A1 (en) * 2022-07-08 2024-01-11 Arcscan, Inc. Using artificial intelligence to detect and monitor glaucoma
CN115018852B (zh) * 2022-08-10 2022-12-06 四川大学 一种基于半监督学习的腹部淋巴结检测方法及装置
CN115082459A (zh) * 2022-08-18 2022-09-20 北京鹰瞳科技发展股份有限公司 对用于屈光度检测的检测模型进行训练的方法及相关产品
CN115587979B (zh) * 2022-10-10 2023-08-15 山东财经大学 基于三阶段注意力网络的糖尿病视网膜病变分级的方法
CN116168255B (zh) * 2023-04-10 2023-12-08 武汉大学人民医院(湖北省人民医院) 一种长尾分布鲁棒的视网膜oct图像分类方法
CN116612339B (zh) * 2023-07-21 2023-11-14 中国科学院宁波材料技术与工程研究所 一种核性白内障图像分级模型的构建装置及分级装置
CN117173222B (zh) * 2023-09-01 2024-07-19 深圳扬奇医芯智能科技有限公司 一种基于扩散模型的配准方法及系统

Family Cites Families (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5784162A (en) 1993-08-18 1998-07-21 Applied Spectral Imaging Ltd. Spectral bio-imaging methods for biological research, medical diagnostics and therapy
US6198532B1 (en) 1991-02-22 2001-03-06 Applied Spectral Imaging Ltd. Spectral bio-imaging of the eye
US6088099A (en) 1996-10-30 2000-07-11 Applied Spectral Imaging Ltd. Method for interferometer based spectral imaging of moving objects
US6276798B1 (en) 1998-09-29 2001-08-21 Applied Spectral Imaging, Ltd. Spectral bio-imaging of the eye
US6943153B1 (en) 1999-03-15 2005-09-13 The Regents Of The University Of California Use of recombinant gene delivery vectors for treating or preventing diseases of the eye
AU2001288372A1 (en) 2000-08-23 2002-03-04 Philadelphia Ophthalmologic Imaging Systems, Inc. System and method for tele-ophthalmology
US20030129164A1 (en) 2001-12-03 2003-07-10 Flannery John G. Expression of glial-derived neurotrophic factor for treatment of diseases of the eye
US7433532B2 (en) 2002-05-01 2008-10-07 Kestrel Corporation Max entropy optimized retinal camera
US6895264B2 (en) 2002-08-26 2005-05-17 Fovioptics Inc. Non-invasive psychophysical measurement of glucose using photodynamics
WO2004021050A2 (en) 2002-08-29 2004-03-11 Kestrel Corporation Hyperspectral imaging of the human retina
US20040221855A1 (en) 2002-10-17 2004-11-11 Paul Ashton Methods for monitoring treatment of disease
US20100172871A1 (en) 2005-02-17 2010-07-08 Flannery John G Muller Cell Specific Gene Therapy
US20110097330A1 (en) 2005-03-11 2011-04-28 Genentech, Inc. Novel Gene Disruptions, Compostitions and Methods Relating Thereto
US7474775B2 (en) 2005-03-31 2009-01-06 University Of Iowa Research Foundation Automatic detection of red lesions in digital color fundus photographs
US8098907B2 (en) 2005-07-01 2012-01-17 Siemens Corporation Method and system for local adaptive detection of microaneurysms in digital fundus images
AU2006280321A1 (en) 2005-08-15 2007-02-22 Genentech, Inc. Gene disruptions, compositions and methods relating thereto
US8114843B2 (en) 2005-11-18 2012-02-14 The Regents Of The University Of California Photoreactive regulator of protein function and methods of use thereof
WO2007118079A2 (en) 2006-04-03 2007-10-18 University Of Iowa Research Foundation Methods and systems for optic nerve head segmentation
WO2008150840A1 (en) 2007-05-29 2008-12-11 University Of Iowa Research Foundation Methods and systems for determining optimal features for classifying patterns or objects in images
US8194936B2 (en) 2008-04-25 2012-06-05 University Of Iowa Research Foundation Optimal registration of multiple deformed images using a physical model of the imaging distortion
US8705826B2 (en) 2008-05-14 2014-04-22 Agency For Science, Technology And Research Automatic cup-to-disc ratio measurement system
US8896682B2 (en) 2008-12-19 2014-11-25 The Johns Hopkins University System and method for automated detection of age related macular degeneration and other retinal abnormalities
US8303115B2 (en) 2009-05-28 2012-11-06 Annidis Health Systems Corp. Method and system for retinal health management
US8775341B1 (en) 2010-10-26 2014-07-08 Michael Lamport Commons Intelligent control with hierarchical stacked neural networks
BR112013013123A2 (pt) 2010-11-26 2018-06-19 Univ Johannesburg Witwatersrand dispositivo de entrega de fármaco
US20150379708A1 (en) 2010-12-07 2015-12-31 University Of Iowa Research Foundation Methods and systems for vessel bifurcation detection
WO2012078636A1 (en) 2010-12-07 2012-06-14 University Of Iowa Research Foundation Optimal, user-friendly, object background separation
CN103458772B (zh) * 2011-04-07 2017-10-31 香港中文大学 视网膜图像分析方法和装置
EP2554662A1 (en) 2011-08-05 2013-02-06 M Maria Pia Cosma Methods of treatment of retinal degeneration diseases
US9849034B2 (en) 2011-11-07 2017-12-26 Alcon Research, Ltd. Retinal laser surgery
WO2013165614A1 (en) 2012-05-04 2013-11-07 University Of Iowa Research Foundation Automated assessment of glaucoma loss from optical coherence tomography
US9867988B2 (en) 2012-08-29 2018-01-16 Arizona Board Of Regents On Behalf Of The University Of Arizona Apparatus for electrical stimulation of a cell and method of use
US20140276025A1 (en) 2013-03-14 2014-09-18 Carl Zeiss Meditec, Inc. Multimodal integration of ocular data acquisition and analysis
US20140314288A1 (en) 2013-04-17 2014-10-23 Keshab K. Parhi Method and apparatus to detect lesions of diabetic retinopathy in fundus images
GB201308131D0 (en) 2013-05-07 2013-06-12 London School Hygiene & Tropical Medicine Ophthalmoscope
US9905008B2 (en) 2013-10-10 2018-02-27 University Of Rochester Automated fundus image field detection and quality assessment
US9008391B1 (en) 2013-10-22 2015-04-14 Eyenuk, Inc. Systems and methods for processing retinal images for screening of diseases or abnormalities
US9836849B2 (en) 2015-01-28 2017-12-05 University Of Florida Research Foundation, Inc. Method for the autonomous image segmentation of flow systems
AU2016308097B2 (en) 2015-08-15 2018-08-02 Salesforce.Com, Inc. Three-dimensional (3D) convolution with 3D batch normalization
WO2017165538A1 (en) * 2016-03-22 2017-09-28 Uru, Inc. Apparatus, systems, and methods for integrating digital media content into other digital media content
WO2018045363A1 (en) * 2016-09-02 2018-03-08 Gargeya Rishab Screening method for automated detection of vision-degenerative diseases from color fundus images

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7575216B2 (ja) 2020-08-25 2024-10-29 株式会社トプコン 眼科画像処理装置、眼科装置、画質評価方法及びプログラム
JP2023539483A (ja) * 2020-11-02 2023-09-14 ▲騰▼▲訊▼科技(深▲セン▼)有限公司 医用画像処理方法、装置、機器、記憶媒体及びコンピュータプログラム
JP7391267B2 (ja) 2020-11-02 2023-12-04 ▲騰▼▲訊▼科技(深▲セン▼)有限公司 医用画像処理方法、装置、機器、記憶媒体及びコンピュータプログラム
US12100140B2 (en) 2021-01-29 2024-09-24 Fujifilm Corporation Information processing apparatus, information processing method, and information processing program
WO2023228229A1 (ja) * 2022-05-23 2023-11-30 日本電気株式会社 画像処理装置、画像処理方法、およびプログラム

Also Published As

Publication number Publication date
BR112019022447A2 (pt) 2020-06-09
WO2018200840A9 (en) 2018-11-29
EP3616120B1 (en) 2024-09-04
US20180315193A1 (en) 2018-11-01
US20200349710A1 (en) 2020-11-05
CN111095261A (zh) 2020-05-01
WO2018200840A1 (en) 2018-11-01
US10719936B2 (en) 2020-07-21
EP3616120A1 (en) 2020-03-04

Similar Documents

Publication Publication Date Title
US10719936B2 (en) System and method for automated funduscopic image analysis
Jiang et al. A review of deep learning-based multiple-lesion recognition from medical images: classification, detection and segmentation
Veena et al. A novel optic disc and optic cup segmentation technique to diagnose glaucoma using deep learning convolutional neural network over retinal fundus images
US20240135517A1 (en) Systems and methods for automated processing of retinal images
WO2018057714A1 (en) Systems, methods and media for automatically generating a bone age assessment from a radiograph
JP2019500110A (ja) 頭部コンピュータ断層撮影における緊急性の特徴の解釈及び定量化
Valizadeh et al. Presentation of a segmentation method for a diabetic retinopathy patient’s fundus region detection using a convolutional neural network
US11883099B2 (en) Noninvasive techniques for identifying choroidal neovascularization in retinal scans
Xiao et al. Major automatic diabetic retinopathy screening systems and related core algorithms: a review
EP4315245A2 (en) Medical image analysis using neural networks
Mathews et al. A comprehensive review on automated systems for severity grading of diabetic retinopathy and macular edema
Suman et al. Automated detection of Hypertensive Retinopathy using few-shot learning
Mookiah et al. Computer aided diagnosis of diabetic retinopathy using multi-resolution analysis and feature ranking frame work
Priya et al. Detection and grading of diabetic retinopathy in retinal images using deep intelligent systems: a comprehensive review
Sridhar et al. Lung Segment Anything Model (LuSAM): A Prompt-integrated Framework for Automated Lung Segmentation on ICU Chest X-Ray Images
US12087445B2 (en) Automatic cervical cancer diagnosis system
Daanouni et al. Automated end-to-end Architecture for Retinal Layers and Fluids Segmentation on OCT B-scans
Niemeijer Automatic detection of diabetic retinopathy in digital fundus photographs
Jain Liver Disease Detection from CT scan images using Deep Learning and Transfer Learning
Bhardwaj et al. A computational framework for diabetic retinopathy severity grading categorization using ophthalmic image processing
Mahbod Towards Improvement of Automated Segmentation and Classification of Tissues and Nuclei in Microscopic Images Using Deep Learning Approaches
Valizadeh et al. Research Article Presentation of a Segmentation Method for a Diabetic Retinopathy Patient’s Fundus Region Detection Using a Convolutional Neural Network
Devidas et al. Segmentation of Optic Disc and Exudates through Hybrid Model by Utilizing Genetic Algorithm and Root Guided Decision Tree
Vij et al. Modified deep inductive transfer learning diagnostic systems for diabetic retinopathy severity levels classification
Agurto Detection and classification of Diabetic retinopathy pathologies in fundus images