JP2022525674A - 単一露光多種撮像と、改善された焦点合わせと、改善された血管造影画像シーケンス表示とを備えた、患者に合わせて調整される眼科撮像システム - Google Patents

単一露光多種撮像と、改善された焦点合わせと、改善された血管造影画像シーケンス表示とを備えた、患者に合わせて調整される眼科撮像システム Download PDF

Info

Publication number
JP2022525674A
JP2022525674A JP2021556553A JP2021556553A JP2022525674A JP 2022525674 A JP2022525674 A JP 2022525674A JP 2021556553 A JP2021556553 A JP 2021556553A JP 2021556553 A JP2021556553 A JP 2021556553A JP 2022525674 A JP2022525674 A JP 2022525674A
Authority
JP
Japan
Prior art keywords
image
images
imaging
fundus
light
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
JP2021556553A
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.)
Carl Zeiss Meditec Inc
Original Assignee
Carl Zeiss Meditec Inc
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 Carl Zeiss Meditec Inc filed Critical Carl Zeiss Meditec Inc
Publication of JP2022525674A publication Critical patent/JP2022525674A/ja
Pending legal-status Critical Current

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/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/0008Apparatus for testing the eyes; Instruments for examining the eyes provided with illuminating means
    • 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/14Arrangements specially adapted for eye photography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • A61B5/14555Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases specially adapted for the eye fundus
    • 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
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/571Depth or shape recovery from multiple images from focus
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • 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/1241Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes specially adapted for observation of ocular blood flow, e.g. by fluorescein angiography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/24Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
    • 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/10024Color image
    • 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/10048Infrared image
    • 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/10064Fluorescence image
    • 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/10Image acquisition modality
    • G06T2207/10141Special mode during image acquisition
    • G06T2207/10152Varying illumination
    • 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/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20101Interactive definition of point of interest, landmark or seed
    • 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/20212Image combination
    • G06T2207/20224Image subtraction
    • 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
    • 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/30168Image quality inspection

Abstract

眼科撮像システムは、連続する走査線の差に基づく自動焦点機構を提供する。本システムはまた、眼底画像内の焦点のユーザ選択を提供する。ニューラルネットワークは、固視角度を決定するために使用され得る、FA又はICGA画像内の視神経乳頭を自動的に特定する。このシステムはまた、羞明患者及びマルチスペクトル撮像オプションに対応するように、複数の撮像モダリティのための追加の走査テーブルを提供する。

Description

本発明は、概して、眼科撮像システムの分野を対象とする。より具体的には、眼科撮像システムのユーザ操作を容易にするための技術に関する。
検眼鏡、光干渉断層撮像法(OCT)、及び他の眼科撮像システムを含む、様々なタイプの眼科検査システムがある。眼科撮像の一例は、スリット走査式眼底撮像又は幅広線眼底撮像であり(例えば、以下の明細書の全ての内容が、参照により本明細書に組み込まれる、(特許文献1)、(特許文献2)、(特許文献3)、(特許文献4)、及び(特許文献5)を参照)、この眼底撮像は、人間の網膜の高解像度生体内撮像を達成するための有望な技術である。撮像手法は、共焦点撮像システムと広視野撮像システムとのハイブリッドである。走査中に網膜の狭く細長い部分を照明することにより、その照明は視認経路の外に留まり、従来の眼底カメラで使用される環状リング照明と比較して網膜のより大部分のより明瞭な視認を可能にする。
良好な画像を得るためには、照明ストリップが良好に焦点合わせされ、弱まることなく眼の瞳孔を通過して眼の眼底に達することが望ましい。これには、システムの慎重な焦点合わせと、眼と眼科撮像システムとの位置合わせとが必要となる。眼科撮像を更に複雑にするのは、患者が撮像中に静止したままでいることが困難であり得ることである。このことは、複数の異なるタイプの画像が必要となる場合、又は患者が高強度照明に対する嫌悪感、例えば羞明を有する場合に、特に問題となる可能性がある。結果として、そのようなシステムを使用するときに高度な能力を達成するのに、概して、かなりの訓練が必要となる。
米国特許第4170398号明細書 米国特許第4732466号明細書 国際公開第2012059236号 米国特許出願公開第2014/0232987号明細書 米国特許出願公開第2015/0131050号明細書
本発明の目的は、眼科撮像又は検査システムの焦点合わせを容易にするためのツールを提供することである。
本発明の別の目的は、眼科画像の撮像を高速化するための様々な方法を提供することである。
本発明の更なる目的は、眼科画像の取り込み中に、患者の不安又は不快感を軽減するための方法を提供することである。
上記目的は、改善された焦点及び機能を有するシステム/方法で達成される。まず、連続線走査間の差に基づく焦点機構は、焦点ずれ測定値を決定する目的で幅広線を細線に効果的に変換する。これらの細線は更に、網膜のトポロジを決定するために使用され得る。
好ましい実施形態では、システムオペレータには、予め指定された複数の焦点補助位置を有するプレビュー画面が提示される。プレビュー画面上に画像が表示され得、且つシステムオペレータは、焦点をよりしっかりと合わせるべくプレビュー画面上の任意の点を自由に選択し得る。次いで、プレビュー画面は、焦点補助位置からの焦点情報を組み合わせることにより、選択点において画像に焦点を合わせる。
いくつかの実施形態では、焦点は更に、患者の固視角度に基づいて調整される。固視角度は、患者の視神経乳頭(ONH)の位置を特定することにより決定され得る。本明細書では、患者の真の固視角度を確認できるように、赤外線(IR)プレビュー画像内のONHを特定するための技術が提供される。
ONHの特定には付加的用途があるので、任意のタイプの画像内のONHを特定することが有益であることに留意されたい。本明細書では、FA(蛍光眼底血管造影)及び/又はICGA(インドシアニングリーン血管造影)画像内及び/又は他の任意の撮像モダリティにおいてONHを見つけ出す方法が提示される。システムは、FA又はICGA画像を取り込む前にIR(赤外線)プレビュー画像を連続的に取り込み得るので、システムはIRプレビュー画像内のONHの位置を特定し、次いで、位置の特定されたONHを取り込まれたFA又はICGA画像に転写する。
これを達成するために、システムは、赤外線画像内のONHの位置を特定する方法を学習するために訓練用カラー画像の大規模なライブラリを利用するために移行訓練を使用する、Uネット(Uネットの概要は、ロンネベルガー(Ronneberger)他著、「Uネット:医用画像セグメンテーションのための畳み込みネットワーク(U-Net:Convolutional Networks for Biomedical Image Segmentation)」コンピュータビジョンとパターン認識(Computer Vision and Pattern Recognition)、2015年、アーカイブ(arXiv):1505.04597[cs.CV]に見られ得る)を提供する。まず、Uネットは、カラー画像のみを使用して訓練される。その後、Uネットの所定数の初期層が固定され、残りの層は、より小さい訓練セットの赤外線画像で訓練される。赤外線画像内のONHの位置特定には、依然として問題がある可能性があり、ONHの位置特定を最適化するための追加のステップが提供される。
FA及びICGA検査を更に改善するために、本システムは、FA又はICGA画像シーケンスの簡略化した取り込みを提供する。本システムは、FA又はICGA画像シーケンスを取り込んでいる間に、オペレータが画像の輝度を調節する必要がないように、非常に高いダイナミックレンジを提供する。生画像が記憶され、シーケンスにおける画像間の相対的輝度情報が生データから決定され、例えばプロットとして、ユーザに正確に提供され得る。加えて、オペレータには、取得された画像シーケンス間の相対的な輝度を維持しながら視認のための記憶された画像を高輝度化するためのオプションが与えられる。
患者に更に対応するように、本システムは、患者に合わせて調整された撮像を提供する。このシステムは、羞明の指標がないか患者の眼の物理的特徴(及び/又は患者の医療記録)を調べる。患者が羞明である候補であるとの疑いがあるか又は以前に羞明と診断されている場合、システムオペレータは、この事実への注意を喚起される。次いで、システムオペレータは、患者の羞明の程度に基づいてシステムに調光撮像スキームを使用させることを選択し得る。
本システムの他の改良形態は、単一の取り込み指示で複数の撮像モダリティを行うための様々な機構を含む。これにより、検査に必要な時間が短縮されるだけでなく、追加の診断に使用され得る追加の情報が提供される。例えば、システムは、全てのカラー画像と共に、IR画像を自動的に取り込み得、このIR画像は、一部の組織の視認性を高めることができる。
本発明のその他の目的及び達成事項は、本発明のより十分な理解と共に、添付の図面と併せて解釈される以下の説明と特許請求の範囲を参照することにより明らかとなり、理解されるであろう。
本明細書で開示される実施形態は例にすぎず、本開示の範囲はそれらに限定されない。1つの請求カテゴリ、例えばシステムにおいて記載される何れの実施形態の特徴も、他の請求カテプリ、例えば方法においても特許請求できる。付属の請求項中の従属性又は後方参照は、形式的な理由のためにのみ選択されている。しかしながら、それ以前の請求項への慎重な後方参照から得られる何れの主題もまた特許請求でき、それによって請求項及びその特徴のあらゆる組合せが開示され、付属の特許請求の範囲の中で選択された従属性に関係なく、特許請求できる。
特許又は出願ファイルは、カラーで作成された少なくとも1つの図面を含む。カラー図面を伴う本特許又は特許出願公開の複写は、請求および必要な手数料の納付に応じて特許商標庁により提供される。
走査される対象物上に走査線が生成され得るときの走査線の簡略化したパターンを図示する。 走査線間のステップサイズが走査線の幅よりも小さい、重複線走査シーケンスを図示する。 2Bは、図2Aの連続する走査線の2つの連続して取り込まれた画像の差し引き効果を図示し、2Cが、2つの走査線からの単一の細線の作成と、線幅を定義する、2次モーメントσの特定とを図示する。 3Aは、重心の変位が眼の深さ(又は高さ)の変化(例えば、逆畳み込み測定値)にどのように対応し得るかを図示し、3Bは、重心の変位が眼の深さ(又は高さ)の変化(例えば、逆畳み込み測定値)にどのように対応し得るかを図示する。 4Aは、複数の焦点補助位置が眼底画像に重ね合わされた、電子ディスプレイ上のプレビューウィンドウ内のプレビュー眼底画像を図示し、4Bは、焦点を合わせるべきプレビュー画像の任意の点をユーザがどのように自由に選択し得るかの例を提供する。 血管造影検査の異なる段階で撮影された2つの血管画像を提供する。 本発明に従って動作中にFA/ICGA画像内のONHを見つけ出すためのセグメンテーション及び位置特定フレームワークを提供する。 7Aは、図6のシステムを使用して特定されたIRプレビュー画像内のONH領域の位置特定を図示し、7Bは、図6のIRプレビュー画像を取り込んだ直後に取り込まれたFA画像を図示する。 血管造影検査の一部として取り込まれた一連の画像内の個々の画像を処理するためのフレームワークを提供する。 各個々の画像がそのそれぞれの個々の白色レベルに従って調節されたFAタイムラプス画像シーケンスを示す。 輝度調節がどの画像にも適用されていない生FA画像シーケンス(図9の画像シーケンスに対応する)を示す。 全ての画像が同じ全体的スケーリング係数(例えば、255/wmax)により高輝度化された、全体的に調節されたFA画像シーケンス(図10の画像シーケンスに対応する)を示す。 同じ全体的スケーリング係数により全ての画像を高輝度化することに等しい、圧縮パラメータ∝がゼロに設定される効果を図示する。 各画像の記憶された白色レベルに従って各画像を個々に自動的に高輝度化することに等しい、圧縮パラメータ∝が1に設定される効果を図示する。 圧縮パラメータ∝の修正が、∝の異なる値について画像シーケンスの輝度にどのように影響を与えるかを図示する。 圧縮パラメータ∝の修正が、∝の異なる値について画像シーケンスの輝度にどのように影響を与えるかを図示する。 圧縮パラメータ∝の修正が、∝の異なる値について画像シーケンスの輝度にどのように影響を与えるかを図示する。 圧縮パラメータ∝の修正が、∝の異なる値について画像シーケンスの輝度にどのように影響を与えるかを図示する。 眼科撮像装置に好適なグラフィカルユーザインタフェース(GUI)の例を提供する。 患者の相対的な光感度に応じて減光量を調整(例えば、調節)するための例示的な工程を提供する。 走査テーブル(例えば、画像取得に使用される強度レベル)が瞳孔の大きさ及び/又は虹彩色に基づいてどのように選択され得るか及び/又は修正され得るかを図示する。 マルチスペクトルオプションが選択された単一の取り込み指示入力に応答して取り込まれ得る、複数の波長にわたる例示的な一連の複数の画像を図示する。 単一の取り込み指示に応答してFA+ICGA機能を提供する2つの従前の方法と本発明とを比較する。 4つのチャネルで構成されたカラー眼底画像を示す。 赤色、緑色、青色、及び赤外チャネルを別個に視認するためのオプションを提供する「チャネル分割」レビュー画面の例を示す。 酸素化及び脱酸素化ヘモグロビン吸光スペクトルを示す。 酸素飽和度マップの例を提供する。 MPOD測定の原理(左)及びMPODプロファイル(右)を図示する。 8つの光源(赤色、緑色、青色、及びIRを各2つずつ)を備えたライトボックス設計を図示する。 眼底を撮影するためのスリットスキャニング眼科システムの例を図解する。 本発明と使用するのに適した眼の3D画像データを収集するために使用される一般型周波数領域光干渉断層撮影システムを図解する。 en face脈管画像の例を示す。 多層パーセプトロン(MLP)ニューラルネットワークの例を図解する。 入力層、隠れ層、及び出力層からなる単純化されたニューラルネットワークを示す。 例示的な畳み込みニューラルネットワークアーキテクチャを図解する。 例示的なU-Netアーキテクチャを図解する。 例示的なコンピュータシステム(又は計算装置又はコンピュータ装置)を図解する。
眼底撮像システム及び光干渉断層撮像法(OCT)撮像システムのセクションで後述するように、様々なタイプの眼科撮像システムがある。本発明の態様は、そのような眼科撮像システムのいずれか又は全てに適用され得る。概して、本発明は、眼科撮像システムの動作及びユーザインタフェースに様々な拡張をもたらす。
焦点及び深さの分析のための細線走査
本発明の一態様は、焦点合わせ用途(例えば、自動焦点)及び逆畳み込み(deconvolution)用途(例えば、トポグラフィ)のための画像測定値を決定する改善された方法を提供する。特定の例として、本発明の向上した焦点合わせ(及び逆畳み込み)技術及び用途は、組み合わせてサンプルの合成画像を構築できる一連の画像セグメントを作成するためにサンプルを横切って走査される、幅広線などの、線状光ビームを使用する眼科撮像システムに適用されるものとして説明されているが、本発明は、他のタイプの眼科撮像システムにも適用され得ることを理解されたい。
例えば、図1は、走査される対象物上に走査線(例えば、スリット又は幅広線)が生成され得るときの走査線の簡略化したパターンを図示する。本例では、走査線は、垂直方向の走査パターンであるV走査で複数の走査線L1~Liを生成するように垂直方向に走査(例えば、横断)される。そのような走査パターンは、線走査眼底撮像システム(又はOCTに基づくシステム)により使用され得、概して、走査線は、走査線(L1~Li)に直交する(例えば、Y軸に沿った)焦点外光のある程度の共焦点抑制(confocal suppression)を維持し得るが、線に沿った(例えば、X軸に沿った)共焦点抑制を欠くことがある。走査線はまた、撮像を改善するために使用され得る。例えば、照明ストリップのエッジの鮮明度は、検出器による取得中に照明が大きく移動しなかった場合(典型的には、走査ビームが、段階的に走査され、取得中に相対的に静止している場合)、線走査システムのために最適化された焦点を見つけ出すために使用され得る。線走査眼底撮像システムに好適な焦点合わせ方法は、全体が参照により本明細書に組み込まれる、国際公開第2018178269号パンフレットに開示されている。別の改善例として、眼の焦点外領域からの、背景光レベル、例えば、迷光レベルを評価するために、走査線により直接照明されない網膜上の位置が検出され(例えば、画像が取り込まれ)得、次いで、この背景レベルは、取り込まれた線画像から差し引かれ得る。線走査撮像装置はまた、瞳分割と組み合わせられている(例えば、全体が参照により本明細書に組み込まれる、ミュラー(Muller)らの米国特許第8488895号明細書を参照)。
要約すると、線走査撮像システムは、眼の網膜上の個々の走査線の画像を記録することにより眼底画像を構築するが、記録された走査線は、より多くの情報を抽出するために使用することができる。例えば、線幅は、焦点の直接の測定値である。この情報は、システムの(例えば、自動)焦点合わせのために又は逆畳み込みアルゴリズムへの追加の入力として使用することができる。また、検出器(例えば、カメラ)上の線位置(重心)は、サンプルの高さ(例えば、特定の点における網膜の深さ)の直接の測定値である。このトポグラフィ情報は、例えば、血管、血管交差、又は網膜上の他の構造(体積を測定すべき腫瘍など)を分析するために使用することができる。
線走査の極めて重要な技術的詳細の1つは、線幅である。線があまりに幅広である場合には、不可能ではないにしても、反射率を位置及び線幅と区別することが困難であり、線の位置特定の不確実性が増し、例えば、幅広線は、広い領域にわたる情報を自動的に平均化し得る。幅広線を使用して画像を走査することには利点があるが、以下に説明するように、幅広線は、位置特定又は幅分析に理想的ではない場合がある。本発明は、走査シーケンスを利用して有効線幅を狭める方法を提供し、この方法は、焦点及び/又はトポグラフィ分析のために使用することができる。
図2Aは、走査線間のステップサイズが走査線の幅よりも小さい、重複線走査シーケンスを図示する。例えば、第1の走査線l1と第2の走査線l2との間の走査ステップサイズS1は、第2の走査線l2が、第1の走査線l1により定義された走査領域に重なるように、走査線l1(及び走査線l2)の幅Wよりも小さい。図2Bは、図2Aの連続する走査線(例えば、隣り合う走査線)の2つの連続して取り込まれた画像(例えば、走査画像)の差し引き効果(例えば、l1から差し引かれたl2)を図示する。隣り合う走査線は、変位幅S1を有する二重線D1/D2をもたらす小さな変位S1がある複数の矩形としてモデル化され得る。ここで、2つの線D1/D2の幅S1及び位置に基づいて、2つの線D1/D2(一方は正であり、他方は負である)を分析することができる。これにより、小さなステップサイズS1の幅広線(例えば、l1及び/又はl2)の走査画像が、同じステップサイズの(二重)細線D1/D2の走査画像に効果的に変換される。
取り込まれる走査画像の線位置を決定するために、取り込まれる走査画像の(光)強度の重心(例えば、標的領域)が算出され得る。これは、非常に高速なアルゴリズムであるが、ノイズの影響を受けやすいことがある。ロバスト性を改善するために、走査画像が、前景及び背景(例えば、線及びノイズ)セグメントにセグメント化され得、且つ前景セグメントについて重心が算出され得る。重心の位置は、深さ感度を高める、サブピクセル精度で算出することができる。
図3A及び図3Bは、重心の変位が眼の深さ(又は高さ)の変化(例えば、逆畳み込み測定値)にどのように対応し得るかを図示する。(予測位置に対する重心の)変位は、網膜トポグラフィの直接の測定値であり、較正システムで実際の高さ(計測)に変換することができる。すなわち、(予測位置からの)線位置の位置ずれは、走査対象物の高さの変化又は凸部(例えば、互いに交差する2つの血管などに起因する眼の深さの変化)によるものであり得る。図3Aは、深さに不規則性のない平面対象物の走査などからの、走査線l3のずれのない走査画像を示す。図3Bは、走査対象物がその表面に不規則性(例えば、凸部)を有することに起因して凸部B1を有する同じ走査線l3の走査画像を示す。取り込まれた走査画像の長さに沿った重心の位置の変化に追従することにより決定できる、凸部B1の高さΔzは、走査対象物の表面における差(凸部)の高さに相当する(例えば、高さの測定値である)。この手法は、カメラ及び照明光学系の較正から利益を得ることがある。
線位置の位置特定と同様に、線幅は、強度分布の2次モーメント(注記:重心は1次モーメントに対応する)を使用して前景セグメントから決定され得る。幅は、少なくとも走査線に直交する方向における、画像の非鮮明度(例えば、焦点ずれ)の直接の測定値である。この幅測定値は、図2Bに図示する作成された細線のいずれかから決定され得る。
例えば、二重細線D1/D2は、分析のために正の線(例えば、D1)又は負の線(例えば、D2)の一方を選択することにより、単一の線に変換され得る。例えば、(l1-l2)の最大値を選択して残りをゼロ(例えば、max[(l1-l2),0])に設定することにより細線D1が単独で定義され得るか、又は(l1-l2)の最小値を選択して残りをゼロ(min[(l1-l2),0])に設定することにより細線D2が単独で定義され得る。図2Cは、max[(l1’-l2’),0]を定義することによる2つの走査線l1’及びl2’からの単一の細線D1’の作成を図示し、且つD1’の線幅とその位置とを定義する、2次モーメントσを特定する。代替的に、又は第2のモーメントσの使用に加えて、線の位置及び幅は、観察データ(例えば、D1’)に機能的形状(例えば、ガウス分布)をフィッティングさせることにより定義され得る。次いで、このフィッティングのパラメータにより位置と幅が得られる。この手法の利点は、走査線の形状(例えば矩形照明)に関する予備知識を考慮に入れることができることである。加えて、ノイズ感度を下げるためにロバストフィッティング法を使用することができる。
このようにして、眼底撮像システムでの幅広線(図2Aに示すものなど)は、小ステップ走査画像の差として細線に(例えば、デジタル処理で)変換され得る。次に、幅(焦点ずれ)及び位置(対象物の深さ)に関して細線を分析することができる。
線幅及び線の変位は、焦点ずれの直接の測定値であり、網膜の深さ、すなわち、網膜の高さに等しいとみなされ得る。これらの測定値は、例えば、自動焦点を提供するために、撮像システム内で使用され得る。図2を参照して説明するような、より高い解像度では、この情報を他の目的にも使用することが可能である(例えば、深さは、血管をセグメント化する、光ディスクをセグメント化する、腫瘍の体積を推定するなどのために使用することができる)。基本的に、本方法は、解像度の向上をもたらす。
この高解像度化は、以下の撮像モデルを使用して説明され得る。
Figure 2022525674000002
ここで、画像は観測された強度であり、pillumは照明の点広がり関数であり、pdetectは検出の点広がり関数であり、且つ
Figure 2022525674000003
は、照明の矩形ストライプである。アスタリスク記号「*」は、畳み込みを表し、ドット「・」は、乗算を表す。このモデルは、照明において線形であり、以下を意味する。
Figure 2022525674000004
解像度を向上させるために、
Figure 2022525674000005
という項が除去され得る。ストライプを差し引く基本概念は、ストライプをδピーク:
Figure 2022525674000006
に変換することである。簡単な差し引きにより2つのデルタが得られるが、上記で説明したように、これは、最小-最大演算により容易に対処される。これにより、撮像システムに入力される照明の幅が最小限に抑えられるが、解像度に対する制限係数として照明のぼやけが依然として残る(δ*pillum=pillum)。
代替案として、重みを使用して正弦照明が概算され得る。
Figure 2022525674000007
これは、(複素)正弦波が畳み込みの固有関数であることによるものである(例えば、何かで畳み込まれた正弦/余弦は、スケーリングされた正弦/余弦である)。ここで、照明のぼやけを取り除くために、pillumを逆畳み込みする、アルファ(α)スケーリング係数が除去され得る。強度の正規化などにより、アルファ(α)スケーリング係数を取り除いた後に、正弦波を(逆の重みで)再結合して、照明としてデルタ又は矩形のいずれかに戻すことができる。
本技術を使用して、(単方向)焦点ずれマップは、走査画像に沿った複数の点における焦点(及び/又は重心位置)を決定することにより構築され得、次いで、焦点ずれマップは、自動焦点ツールとして又は逆畳み込みアルゴリズムを改良するために使用され得る。
代替的に、単一の走査線を分析する代わりに、3D走査で使用されるパターンなどの構造化光パターンが生成され得る。照明の矩形形状がぼやけている(すなわち、よりガウスベル(又は凸部)形状である)ので、交互に代わる符号(例えば、正符号及び負符号)と組み合わせて、正弦照明をリサンプリングすることができる。そのような3つのパターンは、深さ情報を抽出するための位相シフトのために使用することができる。別の例は、相関関係により矩形パターンから近似させることができる三角形の照明パターンの使用である。
上記の例では、照明のぼやけを取り除く(又は軽減する)ために、振幅(例えば、強度)の測定値であり得る、アルファ(α)スケーリング係数を除去したが、アルファ(α)スケーリング係数は、照明の焦点ずれに伴って変化し(例えば、焦点ずれを表し、又は焦点ずれに関連し)、それゆえ、システムの焦点ずれの直接の測定値を提供し得ることが分かっている。例えば、焦点ずれ測定値は、正弦波状照明の振幅に基づいて決定され得る。更に、アルファ(α)スケーリング係数は周波数に依存することが分かっている。結果として、異なる周波数の正弦波を使用することにより、システムの完全な(又は実質的に完全な)点広がり関数を表すために使用され得る一連のアルファ値が得られる。振幅は使用される周波数に依存するので、振幅情報を使用して、システム光学系、眼、及びストライプ幅を含み得る、疑似変調伝達関数(MTF)も生成され得る。これらのアルファ(α)値は、自動焦点合わせの目的及び/又は深さ測定の目的(例えば、表面トポグラフィ/トポロジ)で使用され得る「アルファマップ」を構成し得る。例えば、異なる周波数は異なるアルファ(α)値(例えば、異なる振幅)をもたらすので、周波数を掃引することにより、トポグラフィマップが取得され得る。トポグラフィ情報は、走査ごとに振幅データの高速フーリエ変換(FFT)を行うことにより決定され得る。したがって、網膜のトポグラフィ及び/又は焦点ずれマップは、画像ストライプ(又は走査)の線形結合を使用して決定され得る。
バランスの取れた多点自動焦点
(眼を撮像/走査する前の予備画像又はプレビュー画像上などの)眼底画像上の複数の点又は位置における上記で説明した焦点ずれを決定する方法のいずれか(又は焦点ずれを測定する他の既知の方法)は、眼底の使用可能な画像を取り込むユーザの(例えば、人間オペレータの)能力を向上させるために、グラフィカルユーザインタフェース(GIU)と共に使用され得る。
広視野眼底カメラの焦点深さは、網膜深さの変化よりも小さい可能性があり、画像の周辺の焦点ずれを起こす原因となる。これは、近視網膜の画像に非常に頻繁に見られる。また、網膜剥離及び腫瘍の場合、ユーザは、腫瘍若しくは隆起した網膜の又は周囲の網膜の焦点画像を取得するのに苦労する可能性がある。眼科撮像システムで使用される従来の自動焦点方法では、(眼底)画像は、画像の中心に焦点が合わせられるが、この方法は、臨床医が網膜の(例えば、中心から外れた)特定の部分に焦点を合わせたいと望む又は全体的に最も良く焦点の合った画像を取得したいと単に望む場合には役に立たない。本発明はこの問題に対処する。
典型的には、自動焦点は、概して、多数の異なる技術を使用してカメラの視野(FOV)の中心において提供される。この手法の限界は、純粋に画像の中心に基づいて焦点が選択され、したがって、視野の中心以外の領域には使用することができないことである。この手法はまた、撮像システムが画像のより多くに焦点を合わせることができることを防止する。理解されるように、眼の後部は平坦ではなく、そのため、カメラのFOV内の眼の眼底の異なる部分は、それらに焦点を合わせるのに異なる焦点設定を必要とし得る。
本発明は、各部分が焦点補助位置を定義する、眼の異なる部分における焦点を決定し得る。次いで、焦点補助位置は、新たな画像を取り込むのに最適な焦点を決定するために使用され得る。本発明は、例えば、上記で且つ/又は国際公開第2018178269号パンフレットで説明されているように、複数の焦点補助/機構のいずれかを使用して眼内の複数の別個の位置(焦点補助位置)における焦点読み取り値を取得し得る(例えば、複数の別個の位置(焦点補助位置)についての焦点ずれ測定値を決定し得る)。各位置における焦点は、患者の固視方向を考慮して、様々な焦点ずれ状態(眼の乱視の領域(不規則な形状の角膜)、飛蚊症、眼の曲率などの)を補償するように調節され得る。追加的に、システムは、ユーザ(撮像システムオペレータ)からの入力命令に従ってどのように画像に焦点を合わせるかを調節し得る。例えば、ユーザがどのように画像に焦点を合わせたいとを望むか(例えば、バランスの取れた焦点合わせ又はある点での焦点合わせ)に応じて、複数の焦点補助位置の各々に対する重みが算出され得る。任意選択で、集められた情報は、単一の焦点読み取り値を算出するために組み合わせられ得る。このようにして、本発明は、眼の任意のユーザ指定部分に焦点を合わせるか、又は撮像システムのFOVのより広い領域に焦点を合わせ得る。
図4Aは、複数の焦点補助位置f~fが眼底画像11に重ね合わされた、電子ディスプレイ上のプレビューウィンドウ(又は表示領域)15内のプレビュー眼底画像11を図示する。例示の目的で7つの焦点補助位置が示されていることと、任意の数の複数の焦点補助位置が提供され得ることを理解されたい。更に、焦点補助位置f~fは、予め指定された(任意選択で固定された)位置に位置決めされ且つ眼底画像11又は表示領域15の予め定められた領域(固定又は可変)にわたり得ることを理解されたい。任意選択で、焦点補助位置f~fは、眼底画像11上に表示され(例えば、眼底画像11に重ね合わされ)得るか、又は表示領域15の視界から隠され得る。すなわち、ユーザは、任意の数の焦点補助位置f~fを視認して選択することが許可され得るか、又は焦点補助位置f~fは、ユーザの視界から隠され得る。ユーザは、好ましくは、画像を取り込むために焦点を合わせるべき眼底画像11内の1つ又は複数の領域(又は点)を選択する。ユーザは、任意の既知の入力装置又は機構(例えば、キーボード、ポインタ、ローラ、タッチスクリーン、音声、入力スクリプトなど)を用いて、この選択を入力し得る。
動作中には、以下に説明するように、撮像システムのカメラの視認経路(例えば、集光経路)に対してオフセットされた、患者の瞳孔内に入る光路を通して、網膜(眼底)上の各焦点補助位置が、(例えば、特別な走査線又は他の形状の光ビームであり得る、それぞれの焦点補助照明により)照明される。患者の眼の屈折力の異常は、眼から外に出る視認経路とは異なるように、眼内に入る焦点補助照明を屈曲させ、結果として、屈折異常に依存する各焦点補助位置の位置(例えば、予測位置と異なる位置)がもたらされ、網膜上の焦点補助照明の位置を鮮明に撮像するのに必要な焦点設定に対するシステムのカメラの焦点ずれ測定値に変換される。
各焦点補助f~fは、他の方向への屈折とは異なる可能性がある、瞳孔内のオフセット方向への屈折成分を感知し、次いで、特に乱視がある場合に、各焦点補助位置での最良の画像鮮明度のための焦点を決定する。人間の平均的な眼は、数ジオプトリの乱視を有し、調節可能な非点収差のない眼底カメラ(又は他の眼科撮像システム)は、典型的には、カメラの中心を見るように眼を回転させたときに平均的な眼の乱視を矯正するだけである。任意の所与の固視標においてなど、カメラに対して眼を回転させたときに、本システムは、眼の乱視と眼底カメラの非点収差補正との不整合を考慮に入れることができる。これは、以下の式を使用して行われ得る。
Figure 2022525674000008
ここで、fは、所与の焦点補助位置における焦点であり、
Hは、水平方向視線配向係数であり、
Fixationは、固視標のX位置であり、
focus aidは、焦点補助のX位置であり、
Vは、上下方向視線配向係数であり、
Fixationは、固視標のY位置であり、
focus aid=焦点補助のY位置である。
したがって、ユーザは、所与の焦点補助位置又はプレビュー眼底画像11上の任意の位置を選択することなどにより、システムのFOV内の任意の所与の点/領域での焦点合わせのためにシステムを最適化することを選択し得る。焦点合わせのためにユーザが選択した点/領域が、プレビュー画像11上に表示される又はされないことがある、特定の焦点補助位置(f~f)に対応していない場合に、ユーザが選択した点に最適な焦点は、ユーザが選択した点の位置に対して重み付けされた焦点補助(f~f)の適切な組み合わせにより決定され得る。代替的に、ユーザは、焦点補助位置の各々に対して所定の重みを適用し得る、プレビュー画像にわたる所定の焦点のバランスを使用することを選択し得る。例えば、バランスの取れた焦点は、以下のように算出され得る。
Figure 2022525674000009
ここで:w=焦点位置における重みであり、
=ある位置における非点収差補正焦点である。
図4Bは、焦点を合わせるべきプレビュー画像11の任意の点をユーザがどのように自由に選択し得るかの例を提供する。本例では、ユーザが入力した焦点は、星により特定される。ユーザは、入力装置(例えば、クリック点)を用いて、又はタッチスクリーンを用いて、所望の焦点を入力し得る。この場合、ユーザが入力した点に対する最適な焦点は、ユーザが入力した点に対する焦点補助位置の重み付けされた組み合わせとして決定され得る。例えば、入力されたクリック点からそれぞれの焦点補助位置f~fまでの個々の距離d~dが決定される。次いで、クリック点(星印)に対する焦点は、以下のように決定され得る。
Figure 2022525674000010
ここで、dは、所望の焦点位置と所与の焦点補助位置との間の距離である。
要約すると、本手法は、眼上の複数の別個の位置(例えば、f~f)での焦点読み取りを行う。次に、各位置での焦点読み取りは、固視角度を考慮して調整される。重みは、ユーザが入力した焦点に対する各焦点位置について算出される(重みは、所与の別個の焦点位置からユーザが選択した焦点までの距離に依存し得る)。最後に、画像取り込みのための焦点調節は、別個の各焦点位置の重み付けされた焦点に基づいて、ユーザが選択した焦点について算出され得る。
システムオペレータは患者に固視点を提供し得るが、患者は、提供された固視点に焦点を維持できないことがあることに留意されたい。すなわち、患者の真の固視点は、システムが提供する固視点と異なることがある。それゆえ、本システムは更に、上述のように、最適な焦点設定を決定するときに患者の真の固視方向を決定するための機構を提供する。例えば、システムは、赤外線撮像を使用して、プレビュー画像11であり得る、検査画像を取得し得る。次いで、システムは、検査(例えば、プレビュー)画像内の視神経乳頭を特定し、特定された視神経乳頭の位置に基づいて患者の視線角度を決定し得る。視神経の特定は、通常は未発達の視神経に関連付けられる、視神経低形成を患者が患っている場合、複雑になることがある。それゆえ、システムは更に、中心窩又は特有の血管系などの、視線角度を示す他の目印を探し得る。次いで、患者の視線は、IRプレビュー画像内の特定された目印の位置と、患者の眼に対するシステムの既知の撮像設定(例えば、照明及び視野角)とに基づいて決定され得る。IRプレビュー画像の使用は、眼がIR波長に敏感ではないので有利であり、IRプレビュー画像は、システムの位置合わせの目的などで、眼の邪魔をせずに連続的に取り込まれ得る。特定の目印は、特定の撮像処理アルゴリズム又は機械学習技術(例えば、深層学習)を使用して特定され得る。可視光画像内の目印は、蛍光眼底血管造影(FA)又はインドシアニングリーン血管造影(ICGA)撮像など、特定のタイプの画像を撮影した場合に、見つけ出すのがより困難であり得るので、このタスクでは、赤外線(プレビュー)画像が、可視光(プレビュー)画像よりも優先される。しかしながら、IRプレビュー画像から決定された視神経乳頭の特定された位置は、取り込まれるFA及び/又はICGA画像に時間的に実質的に近接して(例えば、取り込まれるFA及び/又はICGA画像の直前又はそれらと同時に)IRプレビュー画像が撮影される場合に、取り込まれたFA及び/又はICGA画像に直接適用され得る。理解されるように、眼の固視角を決定するためのこのIR撮像の使用では、OCTなどの、他の眼科撮像システムと共に使用され得る。更に、IR画像は、本明細書では固視検出のために使用され、したがって眼の眼底(又は後部)を撮像し、眼の瞳孔(又は前部)を撮像する、視線追跡には使用されないことに留意されたい。
広視野眼底画像内の視神経乳頭の中心を見つけ出すためのUネット
本明細書では、深層学習を使用して、広視野眼底画像内などの、視神経乳頭の中心を特定する例が提供される。上述のように、視神経乳頭の中心は、患者の視線方向を特定するために使用され得る。
視神経乳頭(ONH)は、眼底画像内に観察される最も目立つ目印の1つであり、眼底画像処理ではONHを自動的に見つけ出す機能が非常に望ましい。ONH位置特定アルゴリズムは、病変を探すことを開始するために標準化オフセット又は関心領域(ROI)を適用する基準を提供することにより中心窩又は他の目印を見つけ出すのに役立つことができる。ONH位置特定で利用可能な文献のほとんどは、カラー(例えば、可視光)眼底画像内のONHの位置を特定することに対処する。しかしながら、いくつかのタイプの眼底画像と、カラー画像に適用可能なONH位置特定技術があり、その技術の多くは、他のタイプの画像にはあまり適していない。例えば、血管造影検査の様々な段階で色素が眼血管系を通って進むにつれてONHの外観が変化するので、蛍光眼底血管造影(FA)画像及びインドシアニングリーン血管造影(ICGA)血管造影画像内の視神経の位置を特定することが非常に困難である。本明細書では、FA及びICGA画像の全ての段階でONHを確実に見つけ出すために赤外線(IR)プレビュー画像に適用されるUネットベースの深層学習アーキテクチャが提供される。
従来のONH位置特定アルゴリズムでは、手作りの特徴/フィルタが使用される。例えば、シンタナヨシン(Sinthanayothin)他著、「デジタルカラー眼底画像からの視神経乳頭、中心窩、及び網膜血管の自動位置特定(Automatic Localisation of the Optic Disk, Fovea, and Retinal Blood Vessels from Digital Colour Fundus Images)」、英国眼科学会誌(Br.J.Ophthalmol.)、第83巻、第8号、p.902~910、1999年では、分散フィルタを使用して、視神経乳頭における血管の存在に起因する強度の急速な変化を利用する。しかしながら、この手法は、カラー眼底画像に対してのみ有効であり、小さなデータセットで検査された。また、この手法は、白色病変及び浮き出た脈絡膜血管を有する眼底では失敗する可能性が高いことも示されている。アキタ(Akita)他著、「眼底画像を理解するコンピュータによる方法(A Computer Method of Understanding Ocular Fundus Images)」、パターン認識(Pattern Recognition)、第15巻、第6号、1982年、431~443頁で提供されている別の手法は、血管追跡を使用してONHの位置を特定する。しかしながら、このアルゴリズムは、血管セグメンテーションアルゴリズムの良好な適用に依存している。アルゴリズムが血管系に依存しているので、異常な血管系は偽検出につながることがある。メンデル(Mendels)他著(「動的輪郭を使用した網膜画像内の視神経乳頭境界の特定(Identification of the Optic Disk Boundary in Retinal Images Using Active Contours)」、アイルランドマシンビジョン及び画像処理会議の会議録(Proc.Irish Machine Vision Image Processing Conf.)、1999年9月、p.103~115)では、ONHの検出のために動的輪郭を使用する。この手法の欠点は、画質の変化に対してロバストではない、モルフォロジフィルタリングにより生成される初期輪郭への依存性にある。アルゴリズムが動的輪郭に基づくので、時間効率も高くない。スカール(Sekar)他著、「網膜眼底画像内の視神経乳頭及び中心窩の自動位置特定」、2008年の第16回ヨーロッパ信号処理会議(2008 16th European Signal Processing Conference)、ローザンヌ(Lausanne)、2008年、p.1~5において提供されている更に別の手法は、モルフォロジ演算とその後のハフ変換とを使用して、ONHを見つけ出す。このアルゴリズムでは、ハフ変換は、所与の限界半径を有する円のフィッティングを行うために使用される。アルゴリズムは、グレイレベル強度の変化が大きい画像領域を見つけ出し、その領域をONHを見つけ出すための基準として使用する。このことは通常の画像(例えば、健康な眼の画像)では問題ではない場合があるが、病変のある画像では、視神経乳頭が、強度変化の大きい唯一の構造ではない可能性があり、このことがエラーにつながる場合がある。ジュ(Zhu)他著、「円のハフ変換を使用した網膜の眼底画像内の視神経乳頭の検出(Detection of the Optic Nerve Head in Fundus Images of the Retina Using the Hough Transform for Circles)」、デジタル画像ジャーナル(J Digit Imaging)、2009年、第23巻(第3号)、p.332~41では、ONHを見つけ出すためにエッジ情報及びハフ変換が使用される。しかしながら、このアルゴリズムは、良好なエッジ情報(画質)を抽出してある半径の円に一致させることができるかどうかに依存するので、変化する画質及びONH関連の病変に対してロバストではない可能性が高い。ランガイヤン(Rangayyan)他著、「ガボールフィルタと位相ポートレート解析とを用いた網膜の眼底画像内の視神経乳頭の検出(Detection of the Optic Nerve Head in Fundus Images of the Retina with Gabor Filters and Phase Portrait Analysis)」、デジタル画像学会誌(Journal of digital imaging)、第23巻(第4号)、p.438~53で説明されている、別の手法は、血管の検出にガボールフィルタを使用し、位相ポートレートモデリングを適用してONHの位置を特定しようと試みる。このアルゴリズムは、具体的には、網膜血管の集合点としてのONHの特性に基づいている。血管又は集合中心点が画質又は病変により不明瞭である場合、このアルゴリズムの性能が影響を受ける。畳み込みニューラルネットワーク(CNN)を使用して視神経の位置を特定するために深層学習アルゴリズムを使用するいくつかの手法も報告されている。CNNの説明が以下に提供される。CNNを使用して視神経を見つけ出す例は、D.ニウ(D.Niu)他著、「眼底画像における深層学習に基づく視神経乳頭の自動位置特定(Automatic Localization of Optic Disc Based on Deep Learning in Fundus Images)」、2017年IEEE第2回信号及び画像処理に関する国際会議(2017 IEEE 2nd International Conference on Signal and Image Processing)(ICSIP)、シンガポール(Singapore)、2017年、p.208~212、ゴンザレス・ヘルナンデス・D.(Gonzalez-Hernandez D.)他著、「健常者及び緑内障患者におけるヘモグロビン含有量を測定するための深層学習に基づく視神経乳頭のセグメンテーション(Segmentation of the Optic Nerve Head Based on Deep Learning to Determine its Hemoglobin Content in Normal and Glaucomatous Subjects)」、臨床的&実験的眼科学会誌(J Clin Exp Opthamol)、第9巻:760、並びにH.S.アルガムディ(H.S.Alghamdi)他著、「眼底画像内の視神経乳頭の異常の自動検出:深層学習法(Automatic Optic Disc Abnormality Detection in Fundus Images:A Deep Learning Approach)」、眼科医用画像解析国際ワークショップ(the Ophthalmic Medical Image Analysis International Workshop)(MICCAI、2016年と同時にギリシャのアテネで開催されたOMIA 2016年)の議事録、アイオワリサーチオンライン(Iowa Research Online)、p.17~24、2016年10月において説明されている。上記の公開された著作物の全ては、それらの全体が参照により本明細書に組み込まれる。しかしながら、これらの公表された著作物の全てにおいて、視神経乳頭がカラー眼底画像内に検出されることに留意されたい。すなわち、これらの手法の全ては、特徴を使用するか、又は(例えば、ONHの)カラー眼底画像を使用して(深層学習)アルゴリズムを訓練する。
本発明のONH位置特定手法は、FA及びICGA画像内などの、カラー画像以外の画像タイプ(例えば、撮像モダリティ)でONHの位置を特定するのに好適であるという点で、上記の例とは異なる。
FA及びICGA血管造影検査では、光反応性色素(例えば、FAでは蛍光色素、及びICGAではインドシアニングリーン色素)を被検者の血流に注入した後に、一連のタイムラプス画像が取り込まれる。高コントラストの画像は、色素を励起するために選択された特定の光周波数を使用して取り込まれる。色素が眼内を流れるときに、眼の様々な部分を明るく輝かせ(例えば、蛍光させ)、眼を通る色素の進行ひいては血流を識別することを可能にする。典型的には、色素が眼内に入っていない(又はほとんど入っていない)間は暗色であり、より多くの色素が眼内に入って眼内を流れるにつれて輝度が増し、色素が眼から出ると再び暗くなることにより、取り込まれた一連の画像が特徴付けられるように、より多くの色素を含む部分がより明るく輝く。結果として、血管造影検査の様々な段階で色素が血管を通って進むにつれてONHの外観が経時的に変化し、この変化は、概して、蛍光眼底血管造影(FA)及びインドシアニングリーン血管造影(ICGA)画像内のONHの位置を特定することを非常に困難にする。例示の目的で、図5は、血管造影検査の異なる段階(又は血管造影検査中の異なる時点)で撮像された、2つの血管画像A1及びA2を提供する。図示のように、画像A1内のONHは、画像A2内のONHとは非常に異なって見え、このことは、血管造影検査の全ての段階に好適なフィルタの作製を複雑にする。したがって、ONHの位置特定で利用可能な文献のほとんどは、概して、それ自体をカラー眼底画像に限定しており、FA又はICGA画像内のONHの位置を特定する難題を避けている。しかしながら、本明細書では、FA及びICGA画像又は任意の他のタイプの画像/走査の全ての段階でONHを確実に見つけ出すためにIRプレビュー画像に適用されるUネットベースの深層学習アーキテクチャが提示される。
概して、眼底/網膜のIRプレビュー画像は、検査のために走査(撮像)を行う前に収集される。上記で説明したように、IRプレビュー画像は、焦点を調節するために使用され得、以下に説明するように、IRプレビュー画像は、患者に対するシステムの正確な位置合わせを確実にし、眼の正しい領域が撮像されることを確実にするために使用され得る。したがって、IRプレビュー画像は、FA及び/又はICGA画像の取り込み直前に収集/記録される。本明細書では、FA又はICGA画像(又は他の任意のタイプの画像)を取り込む前にIRプレビュー画像内に直ぐに見つけ出されるONHの位置は、取り込まれたFA又はICGA画像内においても同じである(変わらないままである)ことが述べられている。すなわち、IRプレビュー画像内のOHNの位置は、FA及びICGA画像内のONHの位置特定のために使用される。本発明は、FA及び/又はICGA画像の取り込みを開始する直前に撮影された(IR画像を使用して訓練されたアルゴリズムを使用して)網膜のIRプレビュー画像内のONHを特定して、IRプレビュー画像からONHの位置をFA及び/又はICGA画像に転写することにより、FA及び/又はICGA画像内のOHNの位置を間接的に特定する。
以前の従来の手法は、カラー眼底画像内のONHの位置を特定することに限定されていた。そして、カラー眼底画像内のONHを特定することは、FA及びICGA画像のような動的撮像モダリティ(撮像タイプ)でONHの位置を特定することと比較して、(かなりの訓練画像が与えられた場合には)比較的簡単である。FA及び/又はICGA撮像の様々な段階の各々について、かなりの訓練画像セット(例えば、1000枚を超える画像)を収集することは、特に光反応性色素の注入を伴うので、非常に困難である。本発明は、これらの困難の一部を回避する。
本発明は、機械学習モデルを訓練するために大量の取得するのが困難な血管造影画像(例えば、FA及び/又はICGA画像)を必要とせずに、ONH(又は他の任意の所定の構成又は構造)の位置をどうにか特定する。本発明の一実施形態では、移行訓練は、IR画像を含むように機械学習モデルの動作を拡張することによりカラー画像内のONHの位置を特定するように訓練される既存の機械学習モデルに基づくように使用され得る。移行訓練は、第1のタスクのために開発された機械学習モデルが、第2のタスクに関して別の機械モデルを訓練するための開始点として再使用される、機械学習方法である。これらの問題に関するニューラルネットワークモデルの開発に必要な膨大な計算及び時間資源を考慮して、予め訓練されたモデルが開始点として使用される。例えば、本機械学習モデルは、まず、カラー画像(又はカラー画像から抽出された個々のカラーチャネル)を使用して訓練し、次いで、IR画像に対する移行訓練がなされる。移行訓練の使用は、本発明の別の新規な特徴である。例えば、本アルゴリズムは、2つのステップで訓練され得る、すなわち、第1に、カラー画像及び/又はカラー画像から抽出されたカラーチャネル(例えば、約2000枚の訓練画像)を使用して訓練され、第2に、IR画像(例えば、約1000枚の画像)に対する移行訓練がなされる。すなわち、アルゴリズムは、最初にカラー画像及び/又は個々のカラーチャネル画像(例えば、モノクロ画像)で訓練され、次いでIR画像での訓練を終了し得る。この理由から、この2ステップ手法は、よくあることであるが、多数のカラー画像がすぐに利用可能であり(例えば、2000枚の画像は、ニューラルネットワークの全てのパラメータを訓練するのに十分な量である)、且つ訓練に利用可能なIR画像の数が少ない場合であり得る。例えば、カラー画像から抽出された赤色チャネル画像はまず、第1のステップでネットワークを訓練するために使用され得、IRプレビュー画像は、同じネットワーク上の機械モデルを移行訓練するために使用され得る。この2ステップ手法では、例えば、カラー画像を使用してニューラルネットワークを訓練した後に、所与の数の初期層が固定され(例えば、重み/パラメータが、その後の訓練で変化しないように固定され)、ニューラルネットワークが再び、IRプレビュー画像を使用して残りの固定されていない層のみで訓練される(又は所与の数の最後の、例えば2つの層のみでの訓練)。この第2のステップでの重み/パラメータ訓練は、第1のステップでの重み/パラメータ訓練よりも制限される。しかしながら、IRプレビュー画像の利用可用性は、概して、カラー画像の利用可能性よりも低いので、ネットワーク内の全ての重み/パラメータを訓練するのに十分なIRプレビュー画像がないことがある。この理由から、一般的に利用可能なモダリティ画像(例えば、カラー画像)を使用して機械モデルを訓練し、いくつかの層を固定し、利用可能なIRプレビュー画像を使用してニューラルネットワークのごく少数の層を移行訓練することが有益である。
代替的に、十分な数のIR(プレビュー)画像(例えば、5000枚を超える画像)が訓練に利用可能である場合には、機械モデルは、IR画像のみで、又は1回(例えば、1つの訓練ステップ)でIR画像とカラー(又はカラーチャネル)画像との組み合わせで訓練され得る。しかしながら、カラー眼底画像が概してIR画像よりも広く利用可能であるので、移行訓練はなお、パフォーマンスを向上させるために使用され得る。実験は、1ステップ工程と比較して、本ステップ工程を使用したときの精度の測定可能な向上を示した。
図6は、本発明に従って動作中に(配備システムにおいて)FA/ICGA画像内のONHを見つけ出すためのセグメンテーション及び位置特定フレームワークを提供する。システムは、ブロックB3においてFA及び/又はICGA画像を取り込むことが後続する、IR(例えば、プレビュー)画像を取り込むことにより、ブロックB1で開始し得る。代替的に、ブロックB3では、IR画像と異なる別のタイプの画像(例えば、FAF、色など)は、FA及び/又はICGA画像の代わりに(又はそれに加えて)撮影され得、且つ本フレームワークは、同様に、他のタイプの画像内のONHを特定するのに効果的である。ブロックB1からの取り込まれたIR画像は、処理のために、既に訓練されたUネット、深層学習機械モデルに(入力検査画像として)送られる(ブロックB7)前に、アンチエイリアシングフィルタ及びダウンサンプリングを適用し得る、任意選択の前処理ブロックB5に送られ得る。
好適なUネットアーキテクチャの例が以下で提供される。ブロックB7のUネットは、検査(IR)画像を受信して検査画像内のONHを特定するように訓練される。好ましくは、Uネットは、上述のように、2ステップ工程で、IR眼底画像(例えば、IRプレビュー画像)の訓練セット及びカラー眼底/チャネル抽出カラー眼底のセットを使用して訓練される。Uネットは、まず、カラー画像及び/又は視神経乳頭セグメント(例えば、手動マーキングにより又はカラー画像内のOHNを特定する自動アルゴリズムを用いて境界が定められたONHを有する画像セグメント)を使用してカラー画像から抽出されたチャネル-チャネルを使用して訓練され得る。この第1の訓練ステップは、Uネットの全ての層に初期重み/パラメータを割り当てる。結果として得られる重みは、手動マーキング又は自動アルゴリズムから生成された視神経乳頭セグメントを用いたIRプレビュー画像を使用して移行訓練される。すなわち、Uネットは、初期重み/パラメータで固定された初期層のいくつかと、IRプレビュー画像を使用して訓練された残りの固定されていない層(例えば、最後の2つの層)とを用いて訓練され得る。したがって、Uネットは、検査(IR)画像を受信して、検査画像内において特定された多数の可能性のあるONHセグメントを出力する。
ブロックB9は、Uネットにより特定された全ての可能性のあるONHセグメンテーションを受信して、サイズ基準に基づいて候補ONHセグメントを選択する。この基準は、(例えば、面積又は直径が)第1の閾値よりも大きく且つ第2の閾値よりも小さい受信した全てのONHセグメントを除外することを含み得る。許容可能なサイズ範囲は、基準範囲に従って決定され得る。
次いで、選択された候補ONHセグメントがブロックB11に送信され、B11では、候補ONHセグメントの1つが真のONHセグメントとして選択される。この選択は、候補ONHセグメントの形状に基づき得る。例えば、最も丸みを帯びた形状を有する候補ONHセグメントが、真のONHセグメントとして選択され得る。この真のONHセグメントは、各候補ONHセグメントの短軸と長軸(半径)との差を決定することなどにより、各候補ONH領域の真円度を測ることにより決定され得る。短軸と長軸との差が最も小さいONHセグメントが、真のONHセグメントとして選択され得る。複数のONHセグメントが同じ真円度を有すると分かった場合には、これらのONHセグメントの中で最も大きいONHセグメントが真のONH領域と指定され得る。
このように特定された真のONHセグメントを用いて、次のステップは、IR(プレビュー)画像内の真のONHセグメントの位置を決定することである。IRプレビュー画像内の候補ONHセグメントの位置特定は、候補ONHセグメントの重心を見つけ出すことにより行われ得る(ブロックB13)。
最後に、ブロックB13では、特定されたONHセグメンテーション及び位置特定は、ブロックB1の取り込まれたIRプレビュー画像並びに/或いはブロックB3の取り込まれたFA及び/又はICGA画像に適用され得る。
図7Aは、図6のシステムを使用して特定されたIRプレビュー画像IR-1内のONH領域ONH-1の位置特定を図示する。図7Bは、図6のIRプレビュー画像IR-1を取り込んだ直後に取り込まれたFA画像FA-1(例えば、血管造影画像)を図示し、例えば、画像FA-1は、眼によりONHの実質的な移動がなされないように十分に短いと判断された、画像IR-1の取り込みの後に続く時間期間内に取り込まれる。図7Bは、画像IR-1からFA画像FA-1上に転写された、位置が特定されたONH領域ONH-1を示す。このようにして、血管造影検査(FA又はICGA)の任意の段階でのONHは、素早く且つ効率的に特定され得る。任意選択で、撮像システムにより支援される場合、追加のIRプレビュー画像は、ONHの位置が、連続的に取り込まれたIR画像から、必要に応じて、連続的に更新され得るように、血管造影検査の複数の(又は全ての)段階を通して連続的に取り込まれ得る。
次いで、取り込まれたFA又はICGA画像シーケンスは、特定されたONH領域が強調表示されるか否かにかかわらず表示され得る。それに関係なく、取り込まれた一連のFA又はICGA画像の表示に伴う困難がある。
血管造影画像シーケンスの自動調節
血管造影撮像、蛍光眼底血管造影(FA)及びインドシアニングリーン血管造影(ICGA)の両方は、注入された色素が眼を通り抜けるときに暗色から明色までの広範囲の光強度を取り込む必要がある。画像シーケンスのこの強度の変化を取り込むようにカメラを設定することは、困難である可能性があり、結果として、いくつかの画像が飽和する(例えば、非常に多くの信号利得が過度に明るい又は色あせた画像を生じさせる)又は過度に薄暗くなる。それゆえ、システムオペレータは典型的には、FA又はICGA検査のための画像シーケンスが取り込まれているときに、カメラの取り込み設定を手動で調節する。例えば、システムオペレータは、現在取り込まれた画像を監視して、画像のシーケンスにおける後続の画像の利得を増減するようにカメラを調節し得る。すなわち、ほとんどの眼底カメラは、オペレータが血管造影中に、照明の輝度、センサ利得、又は開口サイズを調節して、信号レベルの多様性に対処し、表示のための良好な輝度を有する画像(すなわち、飽和していない又は過度に薄暗くない画像)を取り込もうとすることを必要とする。これにはいくつかの明確な欠点がある。第1の欠点は、カメラの設定を調節する(例えば、輝度を制御する)必要性によりオペレータの注意がそらされることである。第2の欠点は、オペレータの過誤又は不作為が、画像の露光不足又は露光過多につながる可能性があることである。第3の欠点は、結果として得られる血管造影画像シーケンスが互いに異なるように高輝度化され得ることであり、(血流を示し且つ疾患の評価に関連し得る)時間の経過に伴う画像強度の真の変化が部分的又は完全に不明瞭であることを意味する。いくつかの眼科撮像システムは、システムオペレータからカメラを調節する負担を軽減するために、自動輝度制御(又は感度制御又は自動利得制御)を提供し得る。このタイプの自動制御は、第1及び第2の欠点を克服するのに役立つ可能性があるが、概して、第3の欠点に対処しない。実際に、自動輝度制御は典型的には個々の画像ごとに適用されるので、取り込まれたシーケンスにおける全ての画像は、同様に明るい画像を提供するように個々に調節され得る。この自動輝度制御は、シーケンスにおける画像間の輝度の自然に発生する相対的な変化を更に隠すので、第3の欠点を悪化させる。
本発明は、血管造影検査中に(異なる画像間で)信号強度の真の比例した変化を維持しながら、画像輝度の自動調節を提供する。追加的に、本発明はまた、特定の眼疾患に臨床的関連性を有し得る、時間に対する光強度レベルの傾向を提供し得る。
第1に、本発明の眼科撮像システムは、好ましくは、高ダイナミックレンジ(光)センサを使用して、血管造影検査(例えば、FA又はICGA)において画像シーケンスを取り込む。本発明の新規な特徴は、システム/ハードウェアの調節なしに(例えば、照明輝度、センサ利得、及び/又は開口サイズに対する調節なしに)血管造影検査中に画像シーケンスにおける全ての画像の取り込み中の露光不足及び飽和を回避するのに十分に広いダイナミックレンジを有するセンサの使用である。システムは、(例えば、将来の処理のために)取り込まれた生画像を記憶するが、画像の生データを変更せずに、光学表示のために各画像を任意選択で調節し得る。本システムは、画像の領域を調べて、ディスプレイの最大限の活用のために画像を最適に調節する「画像高輝度化係数」を決定し得る。システムは、一連の画像内の画像間の相対的な輝度が維持されるという制約の中で、ディスプレイの最大限の活用のために一連の画像における全ての画像を調節する「シーケンス高輝度化係数」を決定するために、一連の血管造影の一環として取得された全ての画像の「画像高輝度化係数」を更に調べ得る。この場合、各画像の個々の「画像高輝度化係数」ではなく、同じ「シーケンス高輝度化係数」が、表示のための各画像に適用され得る。「シーケンス高輝度化係数」又は生データを使用して、本システムは更に、色素注入からの経過時間に対する画像強度レベルの傾向を決定して記録し得る。
本システムは、血管造影のために最適化される。本システムは、血管造影に必要な蛍光範囲に対応するのに十分なダイナミックレンジを提供する。結果として、画像が露光過多になるリスクがごく僅かであり、これにより、血管造影中に照度又は検出器感度を常に調節する必要性がない。それゆえ、オペレータは、画像取り込みシーケンスの輝度又は感度を監視又は制御する必要性に気を取られることは決してない。光センサに対する調節がなされないので、取り込まれた画像シーケンスは、医師により検討される場合、多くの疾患/病状において重要である、血管造影を通して取得された画像強度の真の相対的な変化を忠実に表す。例えば、静脈又は動脈の閉塞は、網膜を通る色素の通過を遅らせることがあり、経時的な画像強度の低下率を低下させる可能性があることを意味する。画像間の蛍光性の相対的な差も、炎症状態を診断するために重要である。対照的に、従来の眼科撮像システムは、データ取得中に輝度又は感度設定を手動又は自動で調節することに起因して画像強度の真の相対的な変化を不明瞭にし、画像を検討する医師が存在しない可能性も高かった。
血管造影画像の強度は、色素の通過段階、色素の投与量、網膜の色素沈着、眼球媒体の状態、及び疾患状態によって大きく異なる。ほとんどの眼底撮像システムは、ダイナミックレンジが制限されており、それゆえ、露光不足又は露光過多を回避するために血管造影中に照明輝度又は検出器感度設定の調節を必要とし得る。血管造影中の蛍光信号の自然な変化(通常は初期段階で最も強く、後期段階へと徐々に減少する)は、眼の健康を評価するのに臨床医にとって有用であり得る、網膜を通る血液の通過動力学に関する情報を含む。
図8は、血管造影検査の一部として取り込まれた一連の画像内の個々の画像を処理するためのフレームワークを提供する。取得中に、取得ウィンドウ内の「自動調節」チェックボックスをオンにするなどにより、システムオペレータにより有効にされた場合、取り込まれた画像は、(例えば、血管造影検査の経過時間を表す画像のグループ又はシーケンスの一部として表示するためではなく)個々の表示のために(グレイレベル範囲の点で)最適化される。この最適化により、専門家が画質を適切に評価するのを妨げる可能性がある、個々に表示される画像が非常に暗くなること又は過飽和になることが防止される。まず、画像20が取り込まれる。画像20の取り込みは、複数の走査線(例えば、ストライプデータ)22を取り込むことと、ストライプデータを合成して(組み合わせて)(未調節)生画像20を作成すること24を含み得る。各画像20について、取得された生画素強度が記憶される。画像を表示するのに最適化された白色レベルが、生画像20に対して決定される。白色レベルは、取得された画像データ20の一部26又は全てを調べて、99.9パーセンタイル(又は他の適切に選択されたパーセンタイル)ピクセル強度を算出することにより決定され得る。決定された最適な白色レベルは、生画像20に関連付けられる画像メタデータ28に記憶される。メタデータ28には、(光)センサの電子ノイズをマスキングするレベルで予め設定され得、且つ構成ファイルにより提供され得る、画像を表示するための黒色レベルも含まれる。画像メタデータ28は、個々の画像20についての自動高輝度化情報を(例えば、画像20の白色レベルに従って)効果的に提供する。個々の画像20の「画像高輝度化係数」は、画像20の個々に決定された白色レベル及び/又は黒色レベルに基づき得る。
表示のために個々の画像20が選択されたときに、個々の画像20は、それに関連するメタデータ28に応じて(例えば、記憶された白色レベル及び/又は黒色レベルに応じて)画像調節30(例えば、自動高輝度化)を受ける。次いで、調節された画像は、電子ディスプレイ34上のレビュー画面32に出力され得る。任意選択で、調節された画像はまた、更なる処理のためにレビューソフトウェア36に送信され得る。
所与の血管造影検査のために全ての撮像が完了した(例えば、全てのタイムラプス画像が取り込まれた)後に、タイムラプス画像はグループとして表示され得る。例えば、グループ画像は、レビュー画面32に表示され得る。各画像の対応するメタデータは個々の表示調節のための情報を提供するが、個々に最適化された表示設定に従って各画像を表示することにより、画像の多くは、時間の経過に伴う注入された色素の移動による視覚的変化がある画像と別の画像とで容易に識別できないような、同様の外観を有する。例えば、図9は、各個々の画像がそのそれぞれの個々の白色レベルに従って調節された(例えば、高輝度化された)FAタイムラプス画像シーケンス40を示す。後期画像は、真の蛍光信号が時間の経過に伴って減少するにもかかわらず、初期画像と同様に明るく見え、それゆえ、後期画像は初期画像よりも暗く見えるべきであることに留意されたい。
色素の移行に起因するある画像から別の画像への視覚的変化は、画像のそれぞれの生データでは利用できるが、画像のそれぞれの表示最適化形式では利用できない。各画像の生データ形式に従って各画像をタイムラプスシーケンスで表示することは、色素移行情報を保持するが、表示最適化が適用されないので、各画像の詳細を識別することが困難であり得る。例えば、図10は、任意の画像(例えば、各画像の生信号レベルを表示する)に輝度調節が適用されていない、生FA画像シーケンス42(図9の画像シーケンスに対応する)を示す。結果として、画像の多くは薄暗く、評価するのが困難である。
グループ表示の場合、本発明の目的は、色素移行に起因する全体的な輝度の点で個々の画像間の相対的な差(例えば、視覚的変化)を保持しながら、可能な限り詳細を表示するように、タイムラプスシーケンスにおける全ての画像をグループとして調節することである。(図9に図示するように)シーケンスにおける全ての画像をそれぞれの白色レベルに従って個々に自動的に高輝度化するのではなく、本発明の一実施形態は、シーケンス内の全ての画像の白色レベル(例えば、各画像の対応するメタデータに記憶された白色レベル)を検査して、それらの白色レベルの中で最大白色レベルを見つけ出し、この最大白色レベルは、本明細書では「wmax」呼ばれる。全体的スケーリング係数(例えば、シーケンスにおける全ての画像に及ぶスケーリング係数)は、利用可能な画像強度レベルの最大数をwmax(例えば、255/wmax)で割ったものとして定義される。次に、この同じ全体的スケーリング係数は、画像間の輝度の相対的な差が保持されるように、シーケンスにおける全ての画像に適用される。図11は、全ての画像が同じ全体的スケーリング係数(例えば、255/wmax)により高輝度化された、全体的に調節されたFA画像シーケンス44(図10の画像シーケンスに対応する)を示す。図示のように、画像は、真の蛍光信号を示す相対的な強度変化を保持するが、図10の生FA画像シーケンス42とは異なり、図11の画像は、視認可能な細部を最適化するために高輝度化されている。
代替実施形態では、シーケンスにおける全ての画像に同じ全体的スケーリング係数を適用するのではなく、各画像は、各画像のそれぞれ記憶された白色レベルに基づいて個々に調節され得るが、各画像の輝度設定はそのグループのwmaxを考慮して調節される。このようにして、シーケンス(例えば、グループ)における画像の相対的な輝度は、画像が概してより明るく且つより解釈しやすいが、血管造影中の蛍光信号の変化のある程度の類似性を依然として保持するように圧縮される。この圧縮高輝度化技術は、カスタマイズされた輝度スケーリング係数をシーケンスにおける個々の画像に適用することにより実施することができ、この輝度スケーリング係数は、次のように計算され得る。
Figure 2022525674000011
ここで、B(T)は、血管造影検査のタイムラプスシーケンス中に時間Tで取得された画像に適用されたカスタマイズされた輝度スケーリングであり、w(T)は、時間Tで取得された画像の白色レベルであり、∝は、所定の範囲内(例えば、[0,1]の範囲内)で選択された圧縮パラメータであり、且つwmaxは、例えば上記で説明したような、全体的スケーリング係数である。圧縮パラメータ∝=0を設定することは、図11に図示するように、体的スケーリング係数wmaxを(修正なしに)適用することに等しく、それに対して、∝=1を設定することは、図9に図示するように、個々に最適化された表示設定に従って(修正なしに)画像を自動的に高輝度化することに等しい。
圧縮パラメータ∝は、個々の血管造影検査中に固定値に設定され得(例えば、好ましくは圧縮パラメータ∝は画像シーケンスの取得中に変更されない)、且つ理想的には臨床医(例えば、システムオペレータ)により選択された値に設定され得る。すなわち、システムオペレータからの入力信号は、圧縮パラメータ∝の値を調節することにより、シーケンスにおける画像間の相対的な輝度圧縮の量を修正し得る。例えば、システムオペレータは、輝度スライドバーを用いて、増分/減分ボタンを用いて、及び/又は∝の数値を直接入力することにより、圧縮パラメータ∝を調節し得る。
図12A~図12Fは、この圧縮高輝度化技術が、圧縮パラメータ∝の異なる値に対するタイムラプスシーケンスにおける個々の画像の輝度にどのように影響を及ぼすかのいくつかの例を提供する。図12A~図12Fの各々では、結果として得られる高輝度化された画像シーケンスが左側に示され、血管造影の経過時間に対する輝度のプロットが右側に示されている。図12Aでは、圧縮パラメータ∝がゼロに設定され、このことは、図11を参照して上記で説明したように、同じ全体的スケーリング係数、例えば、この場合は「wmax」により、全ての画像を高輝度化することに等しい。比較の目的で、取り込まれた生データも記憶されるので、複数のプロットが示される。プロット50は、取り込まれた一連の画像内の時間の経過に伴う(例えば、ゼロから10分までの)生画像間の相対的な輝度を図示する。プロット52は、概して画像に同様の輝度をもたらす、各画像の記憶された白色レベルに従って各画像が自動的に輝度調節されたときの、画像間の相対的な輝度を図示する。プロット54は、同じwmax値に従って全ての画像が輝度調節されたときの、画像間の相対的な輝度に対応する。図示のように、wmaxプロット54の相対的な差は、生画像プロット50の相対的な差に従う。
図12Bでは、圧縮パラメータ∝が1に設定され、このことは、図9を参照して上述したように、各画像の記憶された白色レベルに従って各画像を個々に自動的に高輝度化することに等しい。この場合も、プロット50は、生画像間の相対的な輝度を図示し、プロット54は、同じ輝度係数wmaxを使用した画像間の相対的な輝度を図示し、且つプロット52は、時間の経過を伴う画像間の相対的な輝度を示し、各画像は、画像自体の個々の白色レベルに基づいて高輝度化される。プロット52により示すように、画像のそれぞれの白色レベル(のみ)に応じて各画像を自動的に高輝度化することにより、時間の経過に伴う画像間の相対的な輝度の差が隠される。
図12C~図12Fは、圧縮パラメータ∝の修正が、∝の異なる値に対する画像シーケンスの輝度にどのように影響を及ぼすかを図示している。図12C~図12Fでは、圧縮パラメータ∝は、図12Cの∝=0.2から図12Fの∝=0.8まで、0.2の増分で増加する。図12C~図12Fの各々は、タイムラプスシーケンスにおける画像(それぞれ56C~56F)間の結果として得られる相対的な輝度のプロットを提供する。比較の目的で、図12C~図12Fの各々は、プロット50(生画像)、プロット52(記憶された白色レベルに基づいて最大限に自動的に高輝度化された画像)、プロット54(wmaxに基づいて高輝度化された画像)を示す。プロット56C~56Fが示すように、∝を増加させることは、生画像プロット50と同様の、画像間の相対的な差を維持しながら、生画像プロット50と比較して、表示された画像シーケンスの輝度を増加させる効果を有する。
グラフィカルユーザインタフェース
上述のように、本システムは、拡張された焦点合わせ機能を提供する。ユーザが焦点合わせのために多数の所定の点/領域の中から選択することを可能にすること又はユーザが焦点合わせのために任意のランダム点を自由に選択することを可能にすることに加えて、本発明は、複数の撮像タイプ(撮像モダリティ)を取り込むための且つ/又は特定の患者に合わせて特定の撮像タイプを修正するための追加の特徴を提供する。例示として、図13は、眼科撮像装置に好適なグラフィカルユーザインタフェース(GUI)の例を提供し、これらの追加の特徴のいくつかを強調する。GUIは、患者の画像/走査を取得する(例えば、取り込む)ために使用され得る取得(インタフェース)画面21を提供する。任意選択で、現在の患者を特定する患者データ23を取得画面21の上部に表示され得る。取得画面21は、複数の表示域(例えば、情報を視認するための表示画面上の囲み領域)、並びに取得パラメータの設定と画像の取り込み及び検討とデータの分析とを通じて機器オペレータをガイドするための表示要素及びアイコンを提供し得る。例えば、複数の焦点補助位置f~fがプレビュー(眼底)画像(例えば、IRプレビュー画像)11に重ね合わされたプレビューウィンドウ15は、表示域内に提供され得る。上記で説明するように、システムオペレータは、ユーザ入力装置(例えば、電子ポインタ装置、タッチスクリーンなど)を用いて焦点を合わせるためのプレビュー画像11上の任意の点(例えば、標的焦点領域)を選択し得る。
所望であれば、(ステップ調節)焦点ボタン33又は焦点スライドバー35などの、追加の焦点オプションが提供され得る。例示的な実施態様では、プレビューウィンドウ15の標的焦点領域、及び焦点ボタン33、又は焦点スライドバー35は、一方の使用により他方が無効にされるように、相互に排他的であり得る。例えば、焦点ボタン33又は焦点スライドバー35の使用により、プレビューウィンドウ15上の以前に手動で選択された任意の標的焦点領域が無効にされ、焦点ボタン33及び/又は焦点スライドバー35に従って眼底画像11全体に従来の全体的焦点調節が適用され得る。同様に、レビュー画面15上で標的焦点領域を割り当てる(例えば選択する)ことにより、焦点ボタン33又は焦点スライドバー35の任意の焦点設定が無効にされ得る。代替的に、プレビューウィンドウ15の標的焦点領域及び焦点ボタン33又はスライドバー35は協働して機能し得る。例えば、ユーザが、焦点を合わせる眼底画像11内の標的焦点領域を選択した場合、本システムは、上述したように、それに応じて、選択された焦点領域における焦点を調節し得る。しかしながら、ユーザが、本システムにより自動的に提供される焦点領域以外の選択された標的焦点領域の焦点設定を更に調節したいと望む場合、ユーザは、焦点ボタン33及び/又はスライドバー35を使用して、選択された標的焦点領域における焦点を更に調節し得る。したがって、本システムは、眼底画像11内(例えば、プレビュー画面15内)の選択された任意の標的焦点領域に対して手動焦点を提供し得る。
取得画面21は、眼科撮像システムと患者との間の位置合わせを容易にするために使用され得る、異なる視野角からの眼の瞳孔のライブ画像(例えば、赤外線画像)の1つ又は複数の瞳孔ストリーム39を表示し得る。例えば、半透明帯域47又は十字線49などの、オーバーレイガイドは、良好な画像取得のために患者の瞳孔が位置決めされるべきである、標的位置又は標的位置範囲を明らかにするために、ライブストリームに重ね合わせることができる。オーバーレイガイドは、撮像のための標的瞳孔位置と患者の瞳孔中心の現在位置との間の現在のオフセットに関してシステムオペレータに知らせ、その結果、是正措置が講じられ得る。
オペレータが取り込み設定に満足した時点で、取り込み作動入力、例えば、取り込みボタン37を用いて画像が取り込まれ得る。
画像が取り込まれるときに、画像は、取得画面21の取り込みビンセクション41にサムネイル形式で表示され得る。フィルタオプションのドロップダウンリストを提供し得る、表示フィルタ43(例えば、左右性、走査(又は画像)タイプなど)は、取り込みビン41内に表示されるサムネイルを絞り込む(選択又は制限する)ために使用され得る。複数の画像の取得を必要とする撮像モードが選択される場合、(例えば、2つの取り込まれた画像を合成し得る、超広視野オプション(UWF)、又は互いにオフセットされた所定数の取り込まれた画像を合成し得る、自動合成)、取り込みビン41には、取り込むべき画像のプレースホルダを示すプレースホルダアイコン45(例えば、破線円)が予め用意され得る。必要な各画像が取り込まれるときに、プレースホルダアイコン45は、撮像モードの実行状況の表示をユーザに提供するために、取り込まれた画像のサムネイルに置き換えられ得る。
しかしながら、理解されるように、システムオペレータには、画像を取り込む(例えば、走査動作を作動させる)前に、複数のオプション設定が提供され得る。例えば、追加の表示要素及びアイコンは、システムオペレータ(ユーザ)が、取得すべき画像(走査)のタイプを選択して、適切な焦点と、眼科撮像システムに対する患者の適切な位置合わせとを確保することを可能にするために提供される。例えば、検査される眼の左右性は、例えば、左右性アイコン25により、表示/選択され得、右眼又は左眼のどちらの眼(例えば、右眼(オキュラスデクスタ(oculus dexter))OD)又は左眼(オキュラスシニスタ(oculus sinister))OS)が検査/撮像されているかを強調表示する。
様々なユーザ選択可能な走査(例えば、撮像)オプションは、取得画面21上に表示され得る。走査オプションは、取得すべき1つ又は複数の画像の撮像FOVをユーザが選択するためのFOVボタン27を含み得る。FOVオプションボタンは、標準的な単一画像操作のための広視野オプション(WF)、2つの画像を取り込んで合成する超広視野オプション(UWF)、所定数(例えば、4つの)画像を収集して合成する予め設定された合成シーケンスを提供する「自動合成」オプション、及び/又はユーザがユーザ指定の合成シーケンスを提示することを可能にするユーザ定義可能な「合成」オプションを含み得る。加えて、チェックボックスは、ユーザがステレオ撮像を行いたいと望む及び/又は外部固視標を使用したいと望むかどうかを示すために提供され得る。外部固視標を選択することにより、撮像中にシステム内の内部固視標が無効にされ得る。
眼底撮像システムは、複数の走査タイプ(例えば、撮像モダリティ)を支援し得る。走査タイプの例としては、フルカラーで眼底ビューを提供し、且つ患者が不快に感じ得る明るい光を必要とする、カラー撮像が挙げられ得る。別のタイプの画像はIR画像であり、IR画像は、患者には見えず、したがって、より快適である。IR画像は、概してモノクロ画像であり、カラー画像よりも細部を特定するのに役立つ可能性がある。いくつかの種類の組織は、特定の波長、例えば500~800nmで蛍光を発するように作られ得る自然光感知分子を有し、これらの組織は、カラー又はIR周波数で容易に定義されないことがある病変の可能性がある領域/構造を特定するために使用され得る。標的型の組織を蛍光させるために選択された特定の波長を使用して眼を撮像することは、眼底自発蛍光撮像(FAF)と呼ばれる。異なる周波数でのFAF撮像は、異なる種類の組織を自動的に蛍光させ、それにより、異なる種類の病変に対する異なる診断ツールを提供し得る。他の撮像モダリティは、上記で説明したように、一連の画像内の血流を追跡するために、注入された色素を利用する、FA及びICGAを含む。
本例では、取得画面21は、複数の撮像モダリティの中から選択するための走査型セクション29を提供する。選択可能な走査タイプの例としては、カラー(例えば、赤色、緑色、及び青色の光成分/チャネルなどの、可視光を使用したトゥルーカラー撮像)、IR(非可視赤外光を使用した、例えばIRレーザの使用による撮像)、FAF-緑色(緑色励起による眼底自発蛍光)、及びFAF-青色(青色励起による眼底自発蛍光)、FAF-NIA(近赤外光を用いた眼底自発蛍光)、FA(例えば、蛍光色素を患者の血流に注入することによる、蛍光眼底血管造影)、FA-NIA(近赤外光を用いた蛍光眼底血管造影)及びICGA(例えば、インドシアニングリーン色素を患者の血流に注入することによる、インドシアニングリーン血管造影)が挙げられ得る。
取得画面21は更に、拡張した眼(散瞳モード)、拡張していない眼(非散瞳モード)、及び/又は光感度を有し得る眼(例えば、羞明を患う患者)のためなどの、様々な画像取り込み設定のための入力を提供し得る。ユーザは、それぞれのモードのためにボタン31を選択することにより、散瞳(Myd)撮像モードと非散瞳(Non-Myd)撮像モードのいずれかを選択し得る。任意選択で、システムは、眼の状態に基づいてどのモードが適切であるかを自動的に選択(例えば、自動選択)し得る。
オペレータは、画像取り込みシーケンスを開始する前に、複数の羞明関連設定を利用し得る。任意選択で、システムは、以前に羞明を患っていると又は光に敏感である候補であると診断されているとして患者を特定し得る/患者にフラグを立て得る。いずれの場合も、システムは、オペレータにその事実を警告し、それに応じて撮像シーケンスが調節されることを提案し得る。例えば、患者が光に敏感である候補である場合、警告アイコン/ボタン61は、オペレータの注意を引くように、点灯され、色を変化させ、且つ/又は点滅し得る。次いで、オペレータは、ボタンを選択し且つ/又は光強度修正オプション63のリストから手動で選択することにより既定の明所視設定を選択し得、それらオプション63の各々は、適用された光強度を画像取得中に調節し得る。患者の医学的条件(例えば、羞明)に応じた画像取得設定の調整(調節)は、患者に合わせて調整された診断の一例である。
患者に合わせて調整された診断
眼科撮像システムは、眼底画像を取得するために、走査パターンなどでの、閃光を使用し得る。例えば、カール ツァイス メディテック インコーポレイテッド社(Carl Zeiss Meditec Inc.)のクラルス500(商標)(CLARUS 500(商標))は、閃光を使用して高解像度の広視野画像を取り込み得る。FAF-青色モード撮像のために使用される閃光などの、いくつかのタイプの閃光は、明るいものとすることができ、視界不良、弱視、及び/又は色覚低下を含む、患者の短期視覚への影響を及ぼすことができる。羞明、例えば、露光による眼の不快感又は疼痛の症状がある患者に対しては、眼科撮像システムの光源は、かなり影響力が大きく、疼痛、片頭痛、吐き気、流涙、及び他の症状を含む、短期的影響を誘発する可能性がある。
典型的には、既知の光感度を有する患者を撮像するときに、システムオペレータは、患者の露光量を低減するためのステップを行い得る。例えば、システムオペレータは、患者の眼の散大を防止し且つ暗室設定で検査(例えば、撮像)を行うことにより、光感度を有する患者に対して診断を実行(例えば、患者の眼を撮像)し得る。可能であれば、システムオペレータは、システムの光強度を低減し、より小さいFOVの撮像を使用しようと試み得るが、これは、画質の低下をもたらす可能性がある。
本明細書では、光感度を有する患者に対する診断モード(例えば、撮像モード)を有する眼科撮像システム及び/又は光に敏感である候補であり得る患者を特定することが可能なシステムが提示される。
本システムは、眼科撮像システムの光出力を低減する、光に敏感な患者のために「調整された」(例えば、羞明モード又は光敏感モード)撮像モードを提供し得るが、診断のための(例えば、最低限)実現可能な画質をもたらす。図14は、患者の相対的な光感度に応じて減光量を調整(例えば、調節)するための例示的な工程を提供する。この工程は、(例えば、患者の医療記録で)光に敏感であるとして特定された患者に対する任意選択の自動化モードと、画像取得画面内(例えば、図13のボタン61及び/又は63)などの、グラフィカルユーザインタフェース(GUI)内でシステムオペレータに手動選択/選択解除オプションを提供することを含み得る。したがって、工程は、患者の記録を確認するなどにより、光敏感モードを自動又は手動で設定すること(ステップ60)により開始し得る。任意選択で、患者の光感度レベルの評価は、虹彩サイズ、虹彩色、及び/又は明/暗網膜検出などの、追加の情報/パラメータを提供すること(ステップ62)により調節/決定され得る。
走査に基づく眼科撮像システムは、眼底上の所与の位置に適用される走査線の光周波数、強度、及び/又は持続時間を指定する走査テーブルを使用し得る。典型的には、1つの走査テーブルは、撮像タイプごとに構築/定義される。しかしながら、本羞明モードの実施態様は、所与の撮像タイプについての現在の/標準的な光設定(例えば、既定の光設定)の一部を使用して1つ又は複数の走査テーブルを生成することを含み得る。走査テーブルは、光強度の低減と持続時間の短縮(例えば、より短時間の閃光)の両方を提供し得る。走査テーブルはまた、患者の瞳孔サイズ及び/又は虹彩色に基づいた固有の走査設定(例えば、追加的な強度及び持続時間の調節)を含み得る。各撮像モード/タイプ(カラー、FAF-B、FAF-G、IGC、FA)での最小露光量に基づく最低限実現可能な画質は、臨床検査に基づいて予め決定され得る。例えば、光感度を有する患者はまず、患者の光感度量に基づいてシステムに代替の走査テーブルの1つを自動的に選択させ得ること(ステップ64)により、光敏感モードを使用して検査(撮像)され得る。代替的に、ユーザは、所定の代替の走査テーブルを選択し得る。次いで、システムは、選択された走査テーブルを使用して患者を走査(撮像)し(ステップ66)、その後、取得画像の評価が行われる(ステップ68)。この評価は、取得画像の画質測定値を決定するなどにより、システムにより自動的に行われ得るか、又はシステムオペレータにより行われ得る。結果として得られた画像が十分な画質を有しない場合には、追加的に、画像は、必要に応じて、通常モード強度まで強度をステップ関数的に増加させて撮像され得る(ステップ70)。この段階的増加の手動選択は、例えば、図13の入力63を用いて行われ得る。
本システムは、瞳孔サイズ、虹彩色の自動検出、自動明暗網膜検出(例えば、網膜のIRプレビュー画像を使用する)などの、任意選択の要素/パラメータ(例えば、任意選択のステップ62で示す)を含み得る。淡色の網膜は、暗色の網膜よりも重度の羞明になる傾向がある。直接的なテーブル相関関係は、羞明指標と所与の画像取得動作のための対応する光/持続時間設定(例えば、走査テーブル選択)との異なる組み合わせの間で作成され得る。代替的に、画質に影響を及ぼすパラメータの組み合わせに基づく最適(最小)電力設定のために数学的方程式が作成され得る。追加的に、本システムは更に、DNAデータを統合して、可能な光感度を示し得る1つのマーカ又はマーカの集合に基づいて診断装置設定を自動的に調節し得る。これは更に、最適な設定のためのマーカセットの数学的評価に拡張することができる。
図15は、走査テーブル(例えば、画像取得に使用される強度レベル)が瞳孔の大きさ及び/又は虹彩色に基づいてどのように選択され得るか及び/又は修正され得るかを図示する。任意選択で、システムの既定の強度及び持続時間設定は、淡色の眼に関連付けられる最低限の設定に設定され得、且つ強度(及び/又は持続時間)は、羞明ではない眼では増加され得る。更に任意選択で、システムの既定の強度及び持続時間設定は、淡色の眼に関連付けられる最低限の設定に設定され得、且つ取り込まれた画質が最小値以下である場合には、強度及び/又は持続時間は、最低限の画質の画像が得られるまで現在の画像取り込みシーケンス内の追加の画像を取り込むために、所定の段階的増分で自動的に増加され得る。
本発明の患者に合わせて調整された診断は、OCT撮像のためにドライアイモードを追加してなど、患者の既知の状態又は疾患に基づく複数の動作モードを用いて既存の眼科撮像システムを増強し得る。本手法は、ストレス又は不安などの不快感が、閃光又は輝度強度により誘発され得る、精神的健康状態を有する患者にも適用され得る。例えば、本発明の患者に合わせて調整された診断は、患者の不快感及び有害な症状を最小限に抑えるために、心的外傷後ストレス障害(PTSD)を持つ患者に適用され得る。
図13に戻ると、走査型セクション29により示すように、本システムは、1つのボタンを用いて複数の画像タイプ(異なる撮像モダリティ)を取得するための機能を更に提供し得る。以下により十分に説明するように、これは、単一の画像取り込みシーケンス内の複数の撮像モダリティを支援する走査テーブルを提供することにより達成され得る。例えば、マルチスペクトルオプション/ボタン63は、波長のスペクトルにわたって複数の画像を取り込むために提供され得る。図16は、マルチスペクトルオプション63が選択された単一の取り込み指示入力に応答して取り込まれ得る、複数の波長にわたる例示的な一連の複数の画像を図示する。別の例は、以下により十分に説明するように、異なる選択波長で複数の画像を取り込むことにより血液中の酸素の測定値を決定する酸素測定オプション/ボタン65であり得る。また、単一の画像取り込みボタンに応答して単一の画像取り込みシーケンスでカラー画像とIR画像の両方を取り込む、カラー+IRオプション/ボタン67が提供される。別のオプションは、単一の画像取り込み指示に応答して、FA画像とICGA画像の両方又は画像のシーケンスを撮影するためのFA+ICGAオプション/ボタン69であり得る。本明細書では、FA+ICGAオプションを実行するための好ましい方法が提供される。
単一の「取り込み」指示に応答したFA画像及びICGA画像の連続取得
蛍光眼底血管造影(FA)とインドシアニングリーン血管造影(ICGA)の両方は、単一の検査の一環として行うことができる。蛍光色素の直前又は直後にICG(インドシアニングリーン)色素を注入することができ、又は混合ボーラスが使用され得る。FA画像は、青色/緑色照明源を使用して取得され得、且つICGA画像は、近赤外照明源を使用して取得される。両方のモダリティを1つの検査に組み込むことにより、患者が椅子に座る時間が短縮される。しかしながら、血管造影の異なる段階を通してFAとICGAの両方の十分な回数の取り込みを行うことは、患者(特により重度の羞明の患者)とオペレータの両方にとって負担が大きくなることがある。この負担は、単一の露光シーケンスでFA画像とICGA画像の両方を取得することが可能である撮像装置を使用することにより軽減することができる。
FA画像とICGA画像の両方を同時に取得することは、それぞれのモダリティに対する異なる光源及び光学フィルタ要件のため、システム設計における課題を提示する。眼球運動のリスクがあるので、医師が眼の状態を評価するのに十分な画質を達成しながら、短時間(約100ms)で露光を完了する必要がある。
網膜に可視光を送る任意の機器に関しては、患者の快適さが、コンプライアンスを確保する上で重要な考慮事項であり、羞明の症状がある患者は、撮像中に過度に動く若しくはまばたきをすることがあり、又は検査の完了を拒否することさえある。J.M.ストリンガム(J.M.Stringham)他著、「羞明についての作用スペクトル(Action spectrum for photophobia)」、米国光学会誌A(JOSA A)、2003年、第20巻(第10号)、p.1852~1858で説明されているように、羞明感は波長が短くなるにつれて増大する。FAに典型的に使用される光の短波長(470~510nm)は、ICGAに使用される近赤外光よりもはるかに不快である。更に、知覚される閃光持続時間が閃光持続時間の関数として増加することが示されている(オサカ エヌ.(Osaka N.)著、「周辺視野における閃光持続時間の関数としての知覚輝度(Perceived brightness as a function of flash duration in the peripheral visual field)」知覚&精神物理学学会誌(Perception & Psychophysics)、1977年、第22巻(第1号)、p.63~69)。それゆえ、患者の快適さを助長するために、(閃光エネルギーを最小限に抑えることのみならず)痛みを伴う可能性のある露光の持続時間を最小限に抑えることが望ましい。
図17は、単一の取り込み指示に応答してFA+ICGA機能を提供する2つの従前の方法と本発明とを比較する。第1の従前の手法82は、青色及び赤外線レーザ源の両方を使用して、眼に光を同時に送り、2つの別々の光検出器を用いて、放射された蛍光を同時に集光する。第2の従前の手法84は、照明及び検出をインターリーブする(又は交互に行う)。すなわち、レーザ源は、FA画像の1つの線が走査され、次いでICGAの1つの線が走査され、次いでFAの第2の線が走査されるなどのように、交互に入れ替えられる。これらの2つの手法は、眼球運動が生じても、取得されたFA及びICGA画像が共に正確に登録されるという利点を有する。この利点は、画像後処理及び視覚化ツールに対して利益をもたらし得る。
しかしながら、これらの2つの従来の手法にも欠点がある。全く同時の手法82は、比較的複雑なハードウェアを必要とし、例えば、別々の検出器及びフィルタがFA及びICGA収集経路の各々に必要とされる。インターリーブ手法84は、(ICGA源よりも患者にとってはるかに厄介である)青色/緑色FA照明源が、露光(すなわち、長い閃光)の全期間にわたって持続すると患者により知覚されるので、患者の快適さの点で欠点を有し得る。この手法は、羞明反応(例えば、まばたき、縮瞳、眼球運動、ベル現象)を引き起こし得る。
現時点で好ましい手法86は、ユーザからの単一の「取り込み」指示に応答して、FA及びICGA画像を順番に取り込む。すなわち、単一の画像取り込みシーケンスで(例えば、1回の閃光時間内に)、完全なICGA画像(又は画像シーケンス)が取り込まれ、その直後に、完全なFA画像(又は画像シーケンス)が取り込まれる。
好ましい実施形態では、必要な赤外線源は患者にとって不快ではないので、ICGA画像が最初に走査される。ICGA画像取得走査が完了すると、FA走査が開始される。これは、両方の光源を用いて一度に眼を照明する(82)、又は両方の画像を構築するために複数回のFA及びICGA取得をインターリーブする(84)、FA+ICGAの同時取り込みのための2つの従前の手法とは異なる。
この連続的FA+ICGA手法86は、(インターリーブFA+ICGA手法84と比較して)知覚可能な可視光露光の時間間隔を短縮するので、患者の羞明不快感をあまり引き起こさず、まばたき及び/又は眼球運動アーチファクトの起こる確率を低減し得る。すなわち、可視光(例えば、FAの場合)の開始から終了までの時間は、同じエネルギーを送るインターリーブされる閃光の場合よりもはるかに短く、それゆえ、患者に羞明が誘発され得る間隔が最小限に抑えられ得る。これは、血管造影検査からの患者の快適さ及び画質の改善につながることがある。本手法86はまた、第1の手法82の複雑なハードウェア要件を回避する。実際には、本手法86は、適切に定義された走査テーブルを用いて実施され得る。
このFA+ICGAに加えて、本システムは、図13を参照して示すように、追加のマルチスペクトル眼底撮像モード63~67を提供し得る。
スリット走査式検眼鏡のための単一取り込みマルチスペクトル撮像
マルチスペクトル眼底撮像は、複数の波長帯の光(例えば、可視及び非可視反射率並びに自発蛍光)を用いて取得された画像から得られた情報を統合することにより、網膜の構造及び機能の可視性及び識別の向上をもたらすことができる。
トゥルーカラー眼底撮像は、網膜検査の標準である。近赤外光は、より長い波長は血液及びメラニンにより吸収され難いので、網膜のより深部に入り込むことができる。したがって、赤外線眼底撮像は、網膜下の特徴(網膜色素上皮(RPE)下の小さなドルーゼンなど)のより良好な視認性を実現することにより、トゥルーカラー撮像に相補的な情報を提供できることがある。
血液酸素化は、マルチスペクトル撮像手法を使用して網膜内において非侵襲的に調べることができる。J.M.ビーチ(Beach)他著、「2波長撮像による網膜血管の酸素測定:色素沈着の補正及び影響(Oximetry of retinal vessels by dual-wavelength imaging:calibration and influence of pigmentation)」、臨床生理学会誌(Journal of Applied Physiology)、第86巻、第2号(1999年)、p.748~758で説明されているように、2つ以上の励起波長(血液吸収が酸素化に反応しない励起波長)は、血管内径及び網膜色素沈着に対して好適な補正が得られた場合に、血液酸素飽和度を推測するのに使用できる、反射率画像を取り込むために使用され得る。
眼底自発蛍光(FAF)撮像では、取得された画像情報は、励起及び集光波長帯のスペクトル成分に応じて異なる。これらの相違は、網膜の状態に関する臨床情報を抽出するために使用することができる。別々の波長帯でFAFを集光することは、M.ハマー著、「加齢黄斑変性症及び糖尿病性網膜症におけるカラー自発蛍光撮像(Color Autofluorescence Imaging in Age-Related Macular Degeneration and Diabetic Retinopathy)」、調査眼科&視覚科学誌(Investigative Ophthalmology&Visual Science)、第49巻、第13号(2008年)、p.4207~4207で説明されているように、例えば、網膜内の蛍光体を区別するために使用することができる。同様に、FAF励起波長帯を変化させることは、デロリ・フランソワ・C.(Delori,Francois C.)他著、「自発蛍光分光法により測定された黄斑色素密度:反射率測定法及び交照法との比較(Macular pigment density measured by autofluorescence spectrometry:comparison with reflectometry and heterochromatic flicker photometry)」、米国光学会誌A(JOSA A)、第18巻、第6号(2001年)、p.1212~1230において説明されているように、黄斑色素の光学密度を推定するのに使用することができる。
本明細書に記載される全ての参照文献は、それらの全体が参照により組み込まれる。
ほとんどの市販の眼底撮像システムでは、マルチスペクトルデータを取得することは、典型的には、オペレータにとって時間がかかり且つ患者にとって不快である、複数の取り込みを必要とする。本システムは、単一の取り込みプロトコルのみを使用して、マルチスペクトル眼底撮像を行うために、広視野スリット走査式眼底撮像を拡張する経路を提供する。
いくつかの従来の眼底撮像システムは、複数の波長での撮像を提供するが、これらは、典型的には、各モードが専用の取り込み入力を必要とする、個々の走査モードとして提供される。したがって、例えば、RGBカラー画像と赤外線画像の両方を取得するには、2つの別々の走査が必要であり、結果として得られる画像は、(走査間の不可避の眼球運動/再位置合わせに起因して)登録されない可能性が高い。従来の酸素測定手法は、同時2波長酸素測定法を用いる。カラーセンサのベイヤー(Bayer)フィルタは、波長チャネルを分離するために使用され、分析ソフトウェアは、血管ヘモグロビン酸素飽和度の算出を行う。学術研究では、黄斑色素評価のために2つの別々のFAF撮像モダリティを用いる従来のシステムが使用されている。CLARUS(商標)500はまた、2つの異なるFAF励起/検出帯域を提供する。しかしながら、各モードは、別々の専用の取り込みを必要とする。
本明細書では、同時のカラー撮像と赤外線撮像のための技術が提示される。
ある用途では、順序付けられたマルチスペクトル走査は、追加の走査としてIR画像を取得する必要性を排除する。代わりに、完全に登録された赤外線画像は、全てのトゥルーカラー走査で自動的に提供される。機器ソフトウェアのレビュー画面における4帯域チャネル分割機能により、追加の画像情報をユーザに利用できるようにし得る。例えば、図18は、4つのチャネルで構成されたカラー(例えば、トゥルーカラー)眼底画像88-Cを示す。すなわち、このカラー画像88-Cは、典型的な赤色チャネル88-R、緑色チャネル88-G、及び青色チャネル88-Bに加えて、赤外チャネル88-IRを追加的に有する。
本発明の眼科撮像システムは、順序付けられた光源及びモノクロカメラを使用して網膜を撮像する。この撮像は、スリット取得ごとにハードウェアに特定の光源のスイッチを入れるように命令する走査テーブルを使用して達成される。この方法は、異なるタイプの光の走査及び順序付けの両方に関して、大きな柔軟性を与える。
走査テーブル指示の順序は、取得されたスリット画像を別々の画像に仕上げることができるように、スリット画像を個々の色によりグループ化するために使用され、次いで、それらの別々の画像は、トゥルーカラー画像を形成する処理において互いに組み合わされる。走査テーブルは、特定の走査位置、励起波長(例えば、チャネル)などを特定するためにシステムにより使用されるパラメータの横列及びテーブルに配置され得る。例えば、トゥルーカラーモードのための典型的な走査テーブルでは、テーブルにおける各横列はスリット取得を指示し、「チャネル」縦列は、その特定の取得のために光の色を追跡し、後に画像カラーチャネルを別々に組み合わせるために使用される。この場合、チャネルは、赤色、緑色、及び青色のLEDの取得を表すようにR、G、及びBであり得る。赤外線画像の取り込みを含むようにカラー画像の取り込みを拡張するために、走査テーブルは、(走査テーブルの「チャネル」縦列に「I」と表される)追加のIR取得を組み込むように修正され得、その結果、カラー画像の「チャネル」縦列におけるR、G、及びBの各シーケンスは、R、G、B、及びIのシーケンスに拡張する。これにより、4チャネルのマルチスペクトル画像が効果的に提供される。したがって、全てのトゥルーカラー画像はまた、走査の一部分として完全に登録されたIR画像を提供する。IR光は患者にとって不快ではないので、閃光輝度に大きな影響を及ぼさないことに留意されたい。複数の励起を伴うFAFの場合、緑色及び青色照明は、ロングパスバリアフィルタ(例えば、送られたFAF信号を収集するために使用される>650nm)を用いて、走査テーブルにおいて同様に順序付けられる。
全てのカラーチャネル成分画像は、単一の閃光シーケンスで取得されるので、(走査中に眼球運動がなければ)完全に登録される。各チャネルビューの選択は、「チャネル分割」レビュー画面において提供することができる。例えば、図19は、オペレータが赤色、緑色、青色、及び赤外チャネルを別個に視認することを可能にするシステムオペレータによる使用のためのオプション92を提供する「チャネル分割」レビュー画面90の例を示す。別の実施形態では、IR画像情報は、カラー画像と組み合わせて合成画像を形成することができる。例えば、カラー画像の赤色チャネルは、「深さが強調された」トゥルーカラー眼底画像を提供するために、IR画像情報と線形的に結合させることができる。
別個の取り込みモードとしてIR撮像を提供する従前のシステムとは異なり、本技術は、全てのトゥルーカラー取り込みの一環としてIR画像を自動的に提供し得る。これには、椅子に座る追加の時間を必要とせず、患者にとって認識できる更なる不快感がなく、結果として得られたIR画像がカラー画像に完全に登録されるという利点がある。この特徴は、認識できるコストを払わずに、IRにより与えられた追加情報が網膜専門家に提供されるので、網膜専門家の間でのIR撮像への関心を高めるのに役立つことがある。
本明細書では、網膜血管の酸素測定のための技術も提示される。血液酸素飽和度は、マルチスペクトル反射率撮像法を使用して網膜内において非侵襲的に調べることができる。基本的な手法では、2つの励起波長のみが必要とされる。1つの励起波長は、ヘモグロビン(血液中の酸素が結合する)による波長の吸収が酸素飽和(等吸収)度に依存しないように選択される。他の励起波長は、酸素ヘモグロビンと脱酸素ヘモグロビンとの間に吸収における有意差を示すはずである。
図20は、酸素化及び脱酸素化ヘモグロビン吸光スペクトルを示す。好適な等吸収点は、390、422、452、500、530、545、570、584、及び797nmに見られる。用いられる波長は、網膜内での光搬送(例えば、散乱)の波長依存性の差を最小限に抑えるために、できるだけ近い波長である必要がある。酸素測定の現時点で好ましい手法は、広視野を提供し、(例えば、従来の方法により使用されるベイヤー(Bayer)フィルタ分離ではなく)順序付けられた照明に加えて、スリット走査技術を使用する。
本システムは、スリットホイールに専用の励起フィルタを用いて酸素測定を達成し得る。フィルタは、2波長の光を通過させるための5~10nmのノッチを備えた、デュアルバンドパスフィルタ又は積層フィルタであり得る。波長ノッチ中心は、波長ノッチ中心が2波長酸素測定に好適であり且つ波長ノッチ中心をライトボックス内の別々のLED(例えば、570nm及び615nm)により提供できるように選択され得る。次に、酸素測定を実施する方法は、以下のようなものであり得る。
a)単一の走査パターンで緑色光と赤色光を順序付け、次いで、単一の取り込みで570nmと615nmの2つの画像(これらを登録する必要はない)を収集する。
b)血管をセグメント化する。
c)ビーチ(Beach)他により用いられた方法を適用する[1](血管直径と網膜色素沈着の変化のある程度の補正を伴う)。
d)酸素飽和度マップの表示。図21は、酸素飽和度マップの例を提供する。
本システムは更に、黄斑色素光学密度を測定するために使用され得る。この場合、2波長励起での自発蛍光画像は、青色及び緑色励起を順序付けることにより、単一の走査で得られる。取得データは、黄斑色素光学密度(MPOD)の客観的な測定を可能にする。MPODは、加齢黄斑変性症(AMD)を発症する患者のリスクを判定するのに臨床的に有用である可能性がある。AMDは、西欧諸国では失明の主な原因である。AMDは、網膜の黄斑領域に影響を及ぼして、中心視の高解像度の色覚の喪失をもたらす、進行性の不治の疾患である。この疾患は、主に50歳を超える人々に影響を与える。発症率は、年齢プロファイル、食習慣、及び眼のより軽度の平均色素沈着に起因して、一般的に西欧諸国では高くなっている。MPODの測定は、特にスクリーニングの設定において、医師にいくつかの利点をもたらす可能性がある。MPOD測定は、医療提供者が栄養補助食品(その多くは市販されている)を提供して患者の食の健康にその範囲を拡大する可能性を開く。現時点で好ましい方法は、単一の画像取り込みを使用して、黄斑色素光学密度の客観的な測定を可能にするために、複数の励起帯域を使用する順序付けられたFAF取り込みを使用し得る。
好ましい実施形態では、2波長励起での自発蛍光画像は、FAF-緑色に通常使用される>650nmのバリアフィルタを用いて青色及び緑色励起を順序付けることにより、単一の走査で得られる。ODスクリーニング設定において小瞳孔撮像を可能にするために、非散瞳モードを使用することができ、閃光に反応する患者の快適さは、走査を黄斑とその周辺領域とに走査を限定することにより最適化することができる。取得データは、黄斑色素光学密度の客観的な測定を可能にするはずである。
図22は、MPOD測定の原理(左)及びMPODプロファイル(右)を図示する。FAF-緑色画像は、FAF-緑色とFAF-青色との間の正規化された差の(所定の色、例えば黄色の)オーバーレイを有する。この正規化された差は、MPODプロファイルを推測するために使用することができる。
本明細書では、マルチスペクトル撮像に関するいくつかの考慮事項が提供される。様々な波長の追加の光源は、単一露光の、順序付けられた、マルチスペクトル撮像のためにライトボックス設計に組み込むことができる。これにより、(多くの露光ではなく)単一露光と完全に登録された画像出力とを伴うマルチスペクトル撮像機能を提供することができる。
図23は、8つの光源(赤色、緑色、青色、及びIRを各2つずつ)を備えたライトボックス設計を図示する。典型的には、(より良好な画像SNR(信号対雑音比)を取得し、反射をより効果的に抑制するために)より大きな瞳孔に可能であれば各色の両方の光源を使用することに対していくつかの利点があるが、各光源の1つのみが、非散瞳モードでの撮像に必要である。しかしながら、ライトボックス設計は、追加の光源を組み込むように修正され得る。この設計は、図16に図示するように、500~940nmの範囲にわたる波長で、LED(発光ダイオード)光源とフィルタとの複数の組み合わせを使用して、撮像をもたらす、図13のマルチスペクトル機能63を提供することができる。
図13に戻ると、システムオペレータからの単一の制御入力に応答して2つ以上のタイプの画像を撮影する複数の多機能(組み合わせ)キー63~69が提供される。これらの組み合わせは、予め決定されており、固定多機能キー63~69として提示されるが、任意選択で、オペレータは、単一の取り込み指示に応答して異なる撮像モダリティの複数の画像を取り込むカスタム多機能キーを定義し得る。例えば、オペレータには、追加のタイプの画像を単一の画像タイプの選択肢にリンクするための設定ウィンドウが提供され得る。異なるタイプの画像が異なるタイプのフィルタを必要とし得ることと、システムが一度に1つのフィルタの使用のみを支援する場合に、これにより、互いにリンクされ得る画像のタイプの数が制限され得る(例えば、同様のフィルタを必要とする画像のタイプのみが互いにリンクされ得る)ことに留意されたい。システムは、それに応じて利用可能な選択肢を制限し得る。代替的に、システムが、各々がそれ自体のフィルタを有する複数の検出器を有する場合に、光スプリッタは、それぞれのフィルタを介して異なる波長帯の光を異なる検出器に向けるために使用され得る。これにより、単一の取り込み指示に応答してリンクされ且つ撮影され得る画像のタイプの数が増加する。
以下に、本発明に好適な様々なハードウェア及びアーキテクチャの説明が提供される。
眼底撮像システム
眼底を撮像するために使用される撮像システムの2つのカテゴリは、投光照明撮像システム(又は投光照明撮像装置)及び走査照明撮像システム(又は走査撮像装置)である。投光照明撮像装置は、閃光ランプを用いるなどにより、被検査物の対象となる視野(FOV)全体を同時に光であふれさせ、フルフレームカメラ(例えば、全体として、所望のFOVを取り込むのに十分なサイズの2次元(2D)光センサアレイを有するカメラ)を用いて被検査物(例えば、眼底)のフルフレーム画像を取り込む。例えば、投光照明眼底撮像装置は、眼の眼底を光であふれさせ、カメラの単一の画像取り込みシーケンスで眼底のフルフレーム画像を取り込む。走査撮像装置は、対象物、例えば眼を横切って走査される走査ビームを提供し、走査ビームは、走査ビームが対象物を横切って走査され、所望のFOVの合成画像を作成するように再構成され、例えば合成され得る一連の画像セグメントを作成するときに、異なる走査位置で結像される。走査ビームは、点、線、又はスリット若しくは幅広線などの2次元領域とすることができる。
図24は、眼水晶体(又は水晶体)CLとは反対側の眼Eの内表面であり且つ網膜、視神経乳頭、黄斑、中心窩、及び後極を含み得る、眼底Fを撮像するためのスリット走査式眼科システムSLO-1の例を図示する。本例では、撮像システムは、走査線ビームSBは、眼底F全体にわたって走査されるように、眼Eの光学構成要素(角膜Crn、虹彩Irs、瞳孔Ppl、及び水晶体を含む)を横断する、いわゆる「スキャン-デスキャン」構成である。投光眼底撮像装置の場合には、スキャナは不要であり、光は一度に所望の視野(FOV)全体に照射される。他の走査構成は当技術分野で知られており、特定の走査構成は本発明にとって重要ではない。図示のように、撮像システムは、1つ又は複数の光源LtSrc、好ましくは、エタンデュが好適に調節された多色LEDシステム又はレーザシステムを含む。任意選択のスリットSlt(調節可能又は静止)は、光源LtSrcの前に位置決めされ、走査線ビームSBの幅を調節するために使用され得る。追加的に、スリットSltは、撮像中に静止したままであり得るか、又は特定の走査に関して若しくは反射抑制に使用される走査中に異なる共焦点レベル及び異なる用途を可能にするために異なる幅に調整され得る。任意選択の対物レンズObjLは、スリットSltの前に配置することができる。対物レンズObjLは、限定されるものではないが、屈折、回折、反射、又はハイブリッドレンズ/システムを含む最先端のレンズのいずれか1つとすることができる。スリットSltからの光は、瞳分割ミラーSMを通過して、スキャナLnScnへ導かれる。走査面と瞳面とをできるだけ近づけて、システムの口径食を低減することが望ましい。2つの構成要素の画像間の光学距離を操作するために、オプションの光学系DLが含まれ得る。瞳分割ミラーSMは、光源LtSrcからスキャナLnScnに照明ビームを通過させ、スキャナLnScnからの検出ビーム(例えば、眼Eから戻る反射光)をカメラCmrに向けて反射し得る。瞳分割ミラーSMのタスクは、照明ビームと検出ビームとを分割し、システム反射の抑制を補助することである。スキャナLnScnは、回転ガルボスキャナ又は他のタイプのスキャナ(例えば、圧電又はボイスコイル、微小電気機械システム(MEMS)スキャナ、電気光学偏向器、及び/又は回転ポリゴンスキャナ)とすることができる。瞳分割がスキャナLnScnの前に行われるか又は後に行われるかに応じて、1つのスキャナが照明経路内にあり、別々のスキャナが検出経路内にある2つのステップに走査を分割することができる。特定の瞳分割配置は、米国特許第9456746号明細書において詳細に説明されており、その全体は参照により本明細書に組み込まれる。
スキャナLnScnから、照明ビームは、眼Eの瞳孔がシステムの画像瞳に結像されることを可能にする1つ又は複数の光学系、この場合は、走査レンズSL及び眼科用又は接眼レンズOLを通過する。概して、走査レンズSLは、複数の走査角度(入射角)のいずれかでスキャナLnScnから走査照明ビームを受け取り、略平面焦点面(例えば、コリメートされた光路)を有する走査線ビームSBを生成する。眼科用レンズOLは、走査線ビームSBを眼Eの眼底F(又は網膜)に集束させ、眼底を撮像し得る。このようにして、走査線ビームSBは、眼底Fを横切って移動する横断走査線を作成する。これらの光学系の可能な構成の1つは、2つのレンズ間の距離がほぼテレセントリックな中間眼底画像(4-f構成)を生成するように選択されるケプラー型望遠鏡である。眼科用レンズOLは、単レンズ、色消しレンズ、又は異なるレンズの配置とすることができる。全てのレンズは、当業者に知られているように、屈折性、回折性、反射性、又はハイブリッドとすることができる。眼科用レンズOL、走査レンズSLの焦点距離、瞳分割ミラーSM及びスキャナLnScnのサイズ及び/又は形状は、所望の視野(FOV)に応じて異なる可能性があり、そのため、例えば、視野に応じて、光学系のフリップ、電動ホイール、又は着脱可能な光学要素を使用することにより、複数の構成要素をビーム経路の内外に切り替えることができる配置を想定することができる。視野の変化により瞳孔上に異なるビームサイズがもたらされるので、FOVの変更に合わせて瞳分割をも変更することができる。例えば、45°~60°の視野は、眼底カメラについての典型的な又は標準的なFOVである。60°~120°以上のより高視野、例えば広視野FOVも実現可能であり得る。広視野FOVは、幅広線眼底撮像装置(BLFI)と光干渉断層撮像法(OCT)などの別の撮像モダリティとの組み合わせに望ましい場合がある。視野の上限は、人間の眼の周りの生理学的条件と組み合わせたアクセス可能な作動距離により決定され得る。典型的な人間の網膜は、FOVが水平140°及び垂直80°~100°であるため、システム上で可能な限り高いFVOに関して非対称の視野を有することが望ましいことがある。
走査線ビームSBは、眼Eの瞳孔Pplを通過して、網膜又は眼底、すなわち表面Fへ導かれる。スキャナLnScn1は、眼Eの横方向位置の範囲が照明されるように、網膜又は眼底F上の光の位置を調整する。反射光又は散乱光(又は蛍光撮像の場合は放射光)は、照明と同様の経路に沿って導かれ、カメラCmrへの検出経路上に集光ビームCBを定義する。
本発明の例示的なスリット走査式眼科システムSLO-1の「スキャン-デスキャン」構成では、眼Eから戻る光は、瞳分割ミラーSMへの途中でスキャナLnScnにより「デスキャン」される。すなわち、スキャナLnScnは、瞳分割ミラーSMからの照明ビームSBを走査して、眼Eを横切る走査照明ビームSBを定義するが、スキャナLnScnは、同じ走査位置で眼Eからの戻り光も受け取るので、戻り光をデスキャンして(例えば、走査動作をキャンセルして)、スキャナLnScnから瞳分割ミラーSMへの非走査(例えば、定常又は静止)集光ビームを定義するという効果があり、瞳分割ミラーSMは、集光ビームをカメラCmrに向けて折り返す。瞳分割ミラーSMにおいて、反射光(又は蛍光撮像の場合には放射光)は、画像を取り込むために光センサを有するデジタルカメラであり得る、カメラCmrへ導かれる、検出経路上への照明光から分離される。撮像(例えば、対物レンズ)レンズImgLは、眼底がカメラCmrに結像されるように検出経路内に位置決めされ得る。対物レンズObjLの場合のように、撮像レンズImgLは、当技術分野で知られている任意のタイプのレンズ(例えば、屈折、回折、反射又はハイブリッドレンズ)であり得る。追加の動作の詳細、特に、画像におけるアーチファクトを低減する方法が、国際公開第2016/124644号パンフレットで説明されており、その内容の全体は参照により本明細書に組み込まれる。カメラCmrは、受信した画像を取り込み、例えば、画像ファイルを作成し、この画像ファイルは、1つ又は複数の(電子)プロセッサ又は計算装置(例えば、図31に示すコンピュータシステム)により更に処理することができる。したがって、集光ビーム(走査線ビームSBの全ての走査位置から戻る)は、カメラCmrにより収集され、フルフレーム画像Imgは、個々に取り込まれた集光ビームの合成から、モンタージュなどにより、構築され得る。しかしながら、照明ビームが眼Eを横切って走査され、集光ビームがカメラの光センサアレイを横切って走査されるものを含む、他の走査構成も想定される。参照により本明細書に組み込まれる、国際公開第2012/059236号パンフレット及び米国特許出願公開第2015/0131050号明細書は、戻り光がカメラの光センサアレイを横切って掃引される設計、戻り光がカメラの光センサアレイを横切って掃引されない設計などの様々な設計を含む、スリット走査式検眼鏡のいくつかの実施形態を説明している。
本例では、カメラCmrは、プロセッサ(例えば、処理モジュール)Procとディスプレイ(例えば、表示モジュール、コンピュータ画面、電子画面など)Dsplとに接続され、プロセッサとディスプレイの両方は、画像システム自体の一部とすることができ、又は、無線ネットワークを含むケーブル又はコンピュータネットワークを介してデータがカメラCmrからコンピュータシステムに渡されるコンピュータシステムなどの、別個の専用の処理及び/又は表示ユニットの一部であり得る。ディスプレイ及びプロセッサは、一体型ユニットとすることができる。ディスプレイは、従来の電子ディスプレイ/画面とするか又はタッチスクリーン型とすることができ、機器オペレータ又はユーザに情報を表示し且つ機器オペレータ又はユーザから情報を受信するためのユーザインタフェースを含むことができる。ユーザは、限定されるものではないが、マウス、ノブ、ボタン、ポインタ、及びタッチスクリーンを含む、当技術分野で知られているような任意のタイプのユーザ入力装置を使用してディスプレイと対話することができる。
撮像の実行中に患者の視線が固定されたままであることが望ましい場合がある。視線の固定を達成する1つの方法は、患者に凝視するように指示できる固視標を提供することである。固視標は、眼のどの領域が撮像されるかに応じて、機器の内部又は外部とすることができる。内部固視標の一実施形態が図24に示されている。撮像のために使用される一次光源LtSrcに加えて、1つ又は複数のLEDなどの、任意選択の第2の光源FxLtSrcは、レンズFxL、走査要素FxScn及び反射器/ミラーFxMを使用して光パターンが網膜に結像されるように位置決めすることができる。固視スキャナFxScnは、光パターンの位置を移動させることができ、反射器FxMは、光パターンを固視スキャナFxScnから眼Eの眼底Fへ導く。好ましくは、固視スキャナFxScnは、所望の固視位置に応じて網膜/眼底上の光パターンを移動させることができるように固視スキャナFxScnがシステムの瞳面に位置するように配置される。
スリット走査式検眼システムは、用いられる光源及び波長選択フィルタリング要素に応じて、異なる撮像モードで動作することが可能である。トゥルーカラー反射率撮像(手持ち式又は細隙灯検眼鏡を使用して眼を検査するときに臨床医により観察される撮像と同様の撮像)は、一連の有色LED(赤色、青色、緑色)を用いて眼を撮像する場合に達成することができる。各カラーの画像は、各走査位置で各LEDをオンにした状態で段階的に構築することができ、又は各カラー画像を別個に完全に撮影することができる。3つのカラー画像は、トゥルーカラー画像を表示するために組み合わせることができ、又は網膜の異なる特徴を強調するために個々に表示することができる。赤色チャネルは、脈絡膜を最も良く強調し、緑色チャネルは網膜を強調し、青色チャネルは網膜前層を強調する。追加的に、特定の周波数の光(例えば、個々の有色LED又はレーザ)は、眼内の異なる蛍光体(例えば、自発蛍光)を励起するために使用することができ、結果として得られる蛍光は、励起波長をフィルタで除去することにより検出できる。
眼底撮像システムはまた、赤外線レーザ(又は他の赤外線光源)を使用するなどにより、赤外線(IR)反射率画像を提供することができる。赤外線(IR)モードは、眼がIR波長に敏感ではないという点で有利である。この赤外線(IR)モードは、機器の位置合わせ中にユーザを補助するために(例えば、プレビュー/位置合わせモードで)眼の邪魔をせずにユーザが画像を連続的に撮影することを可能にし得る。また、IR波長は、組織を通る透過性が高く、且つ脈絡膜構造の視覚化の改善をもたらし得る。加えて、蛍光眼底血管造影(FA)及びインドシアニングリーン血管造影(ICG)撮像は、蛍光色素が被検者の血流に注入された後に画像を収集することにより達成することができる。
光干渉断層撮影システム
眼底写真、眼底自発蛍光(FAF)、蛍光眼底造影法(FA)に加えて、眼科画像はその他の撮像モダリティ、例えば光干渉断層撮影法(OCT)、OCT眼底造影法(OCTA)、及び/又は眼エコー写真によって製作されてもよい。本発明又は本発明の少なくとも一部は、当業界で理解されるように、若干の改良を加えて、これらのその他の眼科撮像モダリティにも応用されてもよい。より具体的には、本発明はまた、OCT及び/又はOCTA画像を生成するOCT/OCTAシステムにより生成された眼科画像にも応用されてよい。例えば、本発明は、en face OCT/OCTA画像にも応用されてよい。眼底画像の例は米国特許第8967806号明細書及び同第8998411号明細書において提供されており、OCTシステムの例は米国特許第6741359号明細書及び同第9706915号明細書において提供され、OCTA撮像システムの例は米国特許第9700206号明細書及び同第9759544号明細書に見られるかもしれず、これらすべての内容の全体を参照によって本願に援用する。万全を期して、例示的なOCT/OCTAシステムの例を本願で提供する。
図25は、本発明との使用に適した眼の3D画像データ収集用の一般型周波数領域光干渉断層撮影(FD-OCT)システムを図解する。FD-OCTシステムOCT_1は、光源LtSrc1を含む。典型的な光源には、時間コヒーレンス長が短い広帯域光源、又は掃引レーザ源が含まれるがこれらに限定されない。光源LtScr1からの光のビームは、典型的に光ファイバFbr1によってサンプル、例えば眼Eを照明するように誘導され、典型的なサンプルは人間の眼内組織である。光源LrSrc1は、スペクトル領域OCT(SD-OCT)の場合の短い時間コヒーレンス長の広帯域光源か、掃引光源OCT(SS-OCT)の場合の波長調整可能レーザ源の何れかとすることができる。光は、典型的には光ファイバFbr1の出力とサンプルEとの間のスキャナScnr1でスキャンされ、それによって光のビーム(破線Bm)はサンプルの画像撮影対象領域を横方向に(x及びyに)スキャンされる。フルフィールドOCTの場合、スキャナは不要であり、光は一度に所望の視野(FOV)全体に当てられる。サンプルから散乱した光は、典型的に照明用の光を案内するために使用されるものと同じ光ファイバFbr1へと集光される。同じ光源LtSrc1から派生する参照光は別の経路に沿って移動し、この場合、これには光ファイバFbr2及び調整可能な光学遅延を有する逆反射板RR1が含まれる。当業者であればわかるように、透過性参照経路も使用でき、調整可能遅延はサンプル又は干渉計の参照アームの中に設置できる。集光されたサンプル光は、典型的にファイバカプラCplr1において参照光と結合され、OCT光検出器Dtctr1(例えば、光検出器アレイ、デジタルカメラ等)内の光干渉を形成する。1つのファイバポートが検出器Dtctr1に到達するように示されているが、当業者であればわかるように、干渉信号のバランス又はアンバランス検出のために様々な設計の干渉計を使用できる。検出器Dtctr1からの出力は、プロセッサCmp1(例えば、コンピューティングデバイス)に供給され、それが観察された干渉をサンプルの深さ情報へと変換する。深さ情報は、プロセッサCmp1に関連付けられるメモリ内に保存され、及び/又はディスプレイ(例えば、コンピュータ/電子ディスプレイ/スクリーン)Scn1に表示されてよい。処理及び保存機能はOCT機器内に配置されてよく、又は機能は収集されたデータが転送される外部処理ユニット(例えば、図20に示されるコンピュータシステム)上で実行されてもよい。このユニットは、データ処理専用とすることも、又はごく一般的で、OCTデバイス装置に専用ではないその他のタスクを実行することもできる。プロセッサCmp1は例えば、ホストプロセッサに供給される前に、又は並行してデータ処理ステップの一部又は全部を実行するフィールドプログラマブルゲートアレイ(FPGA)、デジタル信号プロセッサ(DSP)、特定用途集積回路(ASIC)、グラフィクス処理ユニット(GPU)、システムオンチップ(SoC)、中央処理ユニット(CPU)、汎用グラフィクス処理ユニット(GPGPU)、又はそれらの組合せを含んでいてよい。
干渉計内のサンプルアームと参照アームは、バルク光学系、ファイバ光学系、又はハイブリッドバルク光学システムで構成でき、また、当業者の間で知られているように、マイケルソン、マッハ・ツェンダ、又は共通光路系設計等、異なるアーキテクチャを有することができる。光ビームとは、本明細書において使用されるかぎり、慎重に方向付けられるあらゆる光路と解釈されるべきである。ビームを機械的にスキャンする代わりに、光の場が網膜の1次元又は2次元エリアを照明して、OCTデータを生成できる(例えば、米国特許第9332902号明細書、ディー.ヒルマン(D.Hillmann)他著、「ホロスコピ-ホログラフィック光干渉断層撮影(Holoscopy-holographic optical coherence tomography)」オプティクスレターズ(Optics Letters)、第36巻(13)、p.2290、2011年、ワイ.ナカムラ(Y.Nakamura)他著、「ラインフィールドスペクトルドメイン光干渉断層撮影法による高速3次元ヒト網膜撮像(High-Speed three dimensional human retinal imaging by line field spectral domain optical coherence tomography)」、オプティクスエクスプレス(Optics Express)、第15巻(12)、p.7103、2007年、ブラスコヴィッチ(Blazkiewicz)他著、「フルフィールドフーリエドメイン光干渉断層撮影法の信号対ノイズ比の研究(Signal-to-noise ratio study of full-field Fourier-domain optical coherence tomography)」、アプライド・オプティクス(Applied Optics)、第44巻(36)、p.7722(2005年)参照)。時間領域システムでは、参照アームは干渉を生じさせるために調整可能な光学遅延を有する必要がある。バランス検出システムは典型的にTD-OCT及びSS-OCTシステムで使用され、分光計はSD-OCTシステムのための検出ポートで使用される。本明細書に記載の発明は、何れの種類のOCTシステムにも応用できる。本発明の様々な態様は、何れの種類のOCTシステムにも、又はその他の種類の眼科診断システム及び/又は、眼底撮像システム、視野試験装置、及び走査型レーザ偏光計を含むがこれらに限定されない複数の眼科診断システムにも適用できる。
フーリエドメイン光干渉断層撮影法(FD-OCT)において、各測定値は実数値スペクトル制御干渉図形(Sj(k))である。実数値スペクトルデータには典型的に、背景除去、分散補正等を含む幾つかの後処理ステップが行われる。処理された干渉図形のフーリエ変換によって、複素OCT信号出力Aj(z)=|Aj|eiφが得られる。この複素OCT信号の絶対値、|Aj|から、異なる経路長での散乱強度、したがってサンプル内の深さ(z-方向)に関する散乱のプロファイルが明らかとなる。同様に、位相φjもまた、複素OCT信号から抽出できる。深さに関する手散乱のプロファイルは、軸方向スキャン(A-スキャン)と呼ばれる。サンプル内の隣接する位置において測定されたA-スキャンの集合により、サンプルの断面画像(断層画像又はB-スキャン)が生成される。サンプル上の横方向の異なる位置で収集されたBスキャンの集合が、データボリューム又はキューブを構成する。特定のデータボリュームについて、速い軸とは1つのB-スキャンに沿ったスキャン方向を指し、遅い軸とは、それに沿って複数のB-スキャンが収集される軸を指す。「クラスタスキャン」という用語は、血流を識別するために使用されてよいモーションコントラストを解析するために、同じ(又は実質的に同じ)位置(又は領域)での反復的取得により生成されるデータの1つのユニット又はブロックを指してよい。クラスタスキャンは、サンプル上のほぼ同じ位置において比較的短い時間間隔で収集された複数のA-スキャン又はB-スキャンで構成できる。クラスタスキャンのスキャンは同じ領域のものであるため、静止構造はクラスタスキャン中のスキャン間で比較的変化しないままであるのに対し、所定の基準を満たすスキャン間のモーションコントラストは血液流として識別されてよい。B-スキャンを生成するための様々な方法が当業界で知られており、これには、水平又はx方向に沿ったもの、垂直又はy方向に沿ったもの、x及びyの対角線に沿ったもの、又は円形若しくは螺旋パターンのものが含まれるがこれらに限定されない。B-スキャンは、x-z次元内であってよいが、z次元を含む何れの断面画像であってもよい。
OCT血管造影法又は関数型OCTにおいて、解析アルゴリズムは、動き又は流れを解析するために、サンプル上の同じ、又はほぼ同じサンプル位置において異なる時間に収集された(例えば、クラスタスキャン)OCTデータに適用されてよい(例えば、米国特許出願公開第2005/0171438号明細書、同第2012/0307014号明細書、同第2010/0027857号明細書、同第2012/0277579号明細書、及び米国特許第6549801号明細書を参照されたく、これらの全ての全体を参照によって本願に援用する)。OCTシステムでは、血流を識別するために多くのOCT血管造影法処理アルゴリズム(例えば、モーションコントラストアルゴリズム)のうちの何れの1つを使用してもよい。例えば、モーションコントラストアルゴリズムは、画像データから導出される強度情報(強度に基づくアルゴリズム)、画像データからの位相情報(位相に基づくアルゴリズム)、又は複素画像データ(複素に基づくアルゴリズム)に適用できる。en face画像は3D OCTデータの2D投射である(例えば、個々のA-スキャンの各々の強度を平均することにより、これによって、各A-スキャンが2D投射内のピクセルを画定する)。同様に、en face脈管画像は、モーションコントラスト信号を表示する画像であり、その中で深さに対応するデータディメンション(例えば、A-スキャンに沿ったz方向)は、典型的にはデータの全部又は隔離部分を加算又は集積することによって、1つの代表値(例えば、2D投射画像内のピクセル)として表示される(例えば、米国特許第7301644号明細書を参照されたく、その全体を参照によって本願に援用する)。血管造影機能を提供するOCTシステムは、OCT血管造影(OCTA)システムと呼ばれてよい。
図26は、en face脈管構造画像の例を示す。データを処理し、当業界で知られるモーションコントラスト法の何れかを用いてモーションコントラストをハイライトした後に、網膜の内境界膜(ILM:internal limiting membrane)の表面からのある組織深さに対応するピクセル範囲を加算して、その脈管構造のen face(例えば、正面図)画像が生成されてよい。
ニューラルネットワーク
前述のように、本発明はニューラルネットワーク(NN)機械学習(ML)モデルを使用してよい。万全を期して、本明細書ではニューラルネットワークについて概説する。発明は、下記のニューラルネットワークアーキテクチャの何れも、単独でも組み合わせても使用してよい。ニューラルネットワーク、又はニューラルネットは、相互接続されたニューロンの(ノードを介した)ネットワークであり、各ニューロンはネットワーク内のノードを表す。ニューロンの集合は層状に配置されてよく、1つの層の出力は多層パーセプトロン(MLP)配置の中の次の層へと順方向に供給される。MLPは、入力データの集合を出力データの集合にマッピングするフィードフォワードニューラルネットワークと理解されてよい。
図27は、多層パーセプトロン(MLP)ニューラルネットワークの例を図解する。その構造は、複数の隠れ(例えば内側)層HL1~HLnを含んでいてよく、これは入力層InL(入力(又はベクトル入力)の集合in_1~in_3を受け取る)を出力層OutLにマッピングし、それが出力(又はベクトル出力)の集合、例えばout_1及びout_2を生成する。各層は、何れの数のノードを有していてもよく、これらはここでは説明のために各層内の円として示されている。この例では、第一の隠れ層HL1は2つのノードを有し、隠れ層HL2、HL3、及びHLnは各々3つのノードを有する。一般に、MLPが深いほど(例えば、MLP内の隠れ層の数が多いほど)、その学習容量は大きい。入力層InLは、ベクトル入力(説明のために、in_1、in_2、及びin_3からなる3次元ベクトルとして示されている)を受け取り、受け取ったベクトル入力を隠れ層のシーケンス内の第一の隠れ層HL1に供給してよい。出力層OutLは、多層モデル内の最後の隠れ層、例えばHLnからの出力を受け取り、ベクトル出力結果(説明のためにout_1及びout_2からなる2次元ベクトルとして示されている)を生成する。
典型的に、各ニューロン(すなわちノード)は1つの出力を生成し、それがその直後の層のニューロンへと順方向に供給される。しかし、隠れ層内の各ニューロンは、入力層から、又はその直前の隠れ層内のニューロンの出力から、複数の入力を受け取るかもしれない。一般に、各ノードはその入力に関数を適用して、そのノードのための出力を生成してよい。隠れ層(例えば、学習層)内のノードは、それぞれの入力に同じ関数を適用して、それぞれの出力を生成してよい。しかしながら、幾つかのノード、例えば入力層InL内のノードは1つの入力しか受け取らず、受動的であってよく、これは、それらが単純にその1つの入力の値をその出力へと中継することを意味し、例えばこれらはその入力のコピーをその出力に提供し、これは説明のために入力層InLのノード内の破線矢印によって示されている。
説明を目的として、図28は、入力層InL’、隠れ層HL1’、及び出力層OutL’からなる単純化されたニューラルネットワークを示す。入力層InL’は、2つの入力ノードi1及びi2を有するように示されており、これらはそれぞれ入力Input_1及びInput_2を受け取る(例えば、層InL’の入力ノードは、2次元の入力ベクトルを受け取る)。入力層InL’は、2つのノードh1及びh2を有する1つの隠れ層HL1’へと順方向に供給し、それが今度は、2つのノードo1及びo2の出力層OutL’に順方向に供給する。ニューロン間の相互接続、又はリンクは(説明のために実線の矢印で示されている)は重みw1~w8を有する。典型的に、入力層を除き、ノード(ニューロン)は入力としてその直前の層のノードの出力を受け取るかもしれない。各ノードは、その入力の各々に各入力の対応する相互接続重みを乗じ、その入力の積を加算し、その特定のノードに関連付けられるかもしれない他の重み又はバイアス(例えば、それぞれノードh1、h2、o1、及びo2に対応するノード重みw9、w10、w11、w12)により定義される定数を加算し(又は、それを乗じ)、その後、その結果に非線形関数又は対数関数を適用することによってその出力を計算してよい。非線形関数は、活性化関数又は伝達関数と呼ばれてよい。複数の活性化関数が当業界で知られており、特定の活性化関数の選択はこの説明には重要ではない。しかしながら、留意すべき点として、MLモデルの演算、ニューラルネットの挙動は重みの値に依存し、これはニューラルネットワークがある入力のための所望の出力を提供するように学習されてよい。
ニューラルネットは、訓練、又は学習段階中に、ある入力にとって望ましい出力を実現するための適当な重み値を学習する(例えば、それを特定するように訓練される)。ニューラルネットが訓練される前に、各重みは個々に初期の(例えば、ランダムな、任意選択によりゼロ以外の)値、例えば乱数シードに割り当てられてもよい。初期重みを割り当てる様々な方法が当業界で知られている。すると、重みは、ある訓練ベクトル入力について、ニューラルネットワークが所望の(所定の)訓練ベクトル出力に近い出力を生成するように訓練される(最適化される)。例えば、重みはバックプロパゲーションと呼ばれる方法によって、何千回もの繰返しサイクルで徐々に調整されてよい。バックプロパゲーションの各サイクルで、訓練入力(例えば、ベクトル入力又は訓練入力画像/サンプル)はニューラルネットワークを通じてフォワードパスが行われて、その実際の出力(例えば、ベクトル出力)が提供される。その後、各出力ニューロン、又は出力ノードのエラーが、実際のニューロンの出力及びそのニューロンのための教師値訓練出力(例えば、現在の訓練入力画像/サンプルに対応する訓練出力画像/サンプル)に基づいて計算される。すると、それはニューラルネットワークを通じて逆方向に(出力層から入力層へと逆方向に)伝搬し、各重みが全体的エラーに対してどの程度の影響を有するかに基づいて重みが更新され、それによってニューラルネットワークの出力は所望の訓練出力に近付く。このサイクルはその後、ニューラルネットワークの実際の出力がその訓練入力のための所望の訓練出力の容認可能なエラー範囲内になるまで繰り返される。理解されるように、各訓練入力は、所望のエラー範囲を実現するまでに多くのバックプロパゲーションイテレーションを必要とするかもしれない。典型的に、エポックは全ての訓練サンプルの1つのバックプロパゲーションイテレーション(例えば、1回のフォワードパスと1回のバックワードパス)を指し、ニューラルネットワークの訓練には多くのエポックが必要かもしれない。一般に、訓練セットが大きいほど、訓練されるMLモデルのパフォーマンスは向上するため、各種のデータ拡張方法が、訓練セットのサイズを大きくするために使用されてよい。例えば、訓練セットが対応する訓練入力画像と訓練出力画像のペアを含む場合、訓練画像は複数の対応する画像セグメント(又はパッチ)に分割されてよい。訓練入力画像及び訓練出力画像からの対応するパッチがペアにされて、1つの入力/出力画像ペアから複数の訓練パッチペアが画定されてよく、それによって訓練セットが拡張される。しかしながら、大きい訓練セットを訓練することによって、コンピュータリソース、例えばメモリ及びデータ処理リソースへの要求が高まる。演算要求は、大きい訓練セットを複数のミニバッチに分割することによって軽減されるかもしれず、このミニバッチのサイズは1回のフォワード/バックワードパスにおける訓練サンプルの数が決まる。この場合、そして1つのエポックは複数のミニバッチを含んでいてよい。他の問題は、NNが訓練セットを過剰適合して、特定の入力から異なる入力へと一般化するその能力が減少する可能性である。過剰適合の問題は、ニューラルネットワークのアンサンブルを作るか、又は訓練中にニューラルネットワーク内のノードをランダムにドロップアウトすることによって軽減されるかもしれず、これはドロップされたリードをニューラルネットワークから有効に除去する。インバースドロップアウト等、各種のドロップアウト調整方法が当業界で知られている。
留意すべき点として、訓練済みのNN機械モデルの演算は、演算/解析ステップの単純なアルゴリズムではない。実際に、訓練済みのNN機械モデルが入力を受け取ると、その入力は従来の意味では解析されない。むしろ、入力の主旨や性質(例えば、ライブ画像/スキャンを画定するベクトル、又は人口構造的説明又は活動の記録等のその他何れかのエンティティを画定するベクトル)に関係なく、入力は、訓練済みニューラルネットワークの同じアーキテクチャ構築(例えば、同じノード/層配置、訓練済み重み及びバイアス値、所定の畳み込み/逆畳み込み演算、活性化関数、プーリング演算等)の対象となり、訓練済みネットワークのアーキテクチャ構築がその出力をどのように生成するかは明らかでないかもしれない。さらに、訓練された重みとバイアスの値は、決定的ではなく、そのニューラルネットワークに付与される訓練のための時間の量(例えば、訓練におけるエポック数)、訓練開始前の重みのランダムな開始値、NNがそこで訓練されるマシンのコンピュータアーキテクチャ、訓練サンプルの選択、複数のミニバッチ間の訓練サンプルの分布、活性化関数の選択、重みを変更するエラー関数の選択、さらには訓練が1つのマシン(例えば、第一のコンピュータアーキテクチャを有する)で中断され、他のマシン(例えば、異なるコンピュータアーキテクチャを有する)で完了したか等、多くの要素に依存する。ポイントは、訓練済みのMLモデルが特定の出力になぜ到達したかの理由は明白でなく、MLモデルがその出力の基礎とする要素を特定しようとする多くの研究が現在行われている、ということである。したがって、ライブデータに対するニューラルネットワークの処理は、単純なステップのアルゴリズムまで減少させることはできない。むしろ、その演算は、その訓練アーキテクチャ、訓練サンプルセット、訓練シーケンス、及びMLモデルの訓練における様々な状況に依存する。
概略すると、NN機械学習モデルの構成は、学習(又は訓練)ステージと分類(又は演算)ステージを含んでいてよい。学習ステージでは、ニューラルネットワークは特定の目的のために訓練されてよく、また訓練例の集合が提供されてよく、これには訓練(サンプル)入力及び訓練(サンプル)出力が含まれ、任意選択により、訓練の進行状況を試験するためのバリデーションの例の集合が含まれる。この学習プロセス中、ニューラルネットワーク内のノード及びノード相互接続に関係付けられる各種の重みが徐々に調整されて、ニューラルネットワークの実際の出力と所望の訓練出力との間のエラーが縮小される。このようにして、多層フィードフォワードニューラルネットワーク(例えば前述のもの)は、何れの測定可能関数を何れの所望の精度までも概算できるかもしれない。学習ステージの結果として得られるのは、学習した(例えば、訓練済みの)(ニューラルネットワーク)機械学習(ML)である。演算ステージで、試験入力(又はライブ入力)の集合が学習済み(訓練済み)MLモデルに提供されてよく、それが学習したことを応用して、試験入力に基づいて出力予測を生成するかもしれない。
図26及び図27の通常のニューラルネットワークと同様に、畳み込みニューラルネットワーク(CNN)もまた、学習可能な重みとバイアスを有するニューロンで構成される。各ニューロンは入力を受け取り、演算(例えば、ドット積)を行い、任意選択によってそれに非線形変換が続く。しかしながら、CNNは、一方の端(例えば入力端)で生の画像ピクセルを受け取り、反対の端(例えば、出力端)で分類(又はクラス)のスコアを提供する。CNNは入力として画像を予想するため、これらはボリューム(例えば、画像のピククセル高さと幅及び、画像の深さ、例えば赤、緑、及び青の3色で定義されるRGB深さ等の色深さ)を扱うように最適化される。例えば、CNNの層は、3次元で配置されるニューロンのために最適化されてよい。CNN層内のニューロンは、完全に接続されたNNのニューロンの全部ではなく、その前の層の小さい領域に接続されてもよい。CNNの最終的な出力層は、フル画像を深さの次元に沿って配置される1つのベクトル(分類)に縮小するかもしれない。
図29は、例示的な畳み込みニューラルネットワークアーキテクチャを提供する。畳み込みニューラルネットワークは、2つ又はそれ以上の層(例えば、層1~層N)の連続として定義されてよく、層は(画像)畳み込みステップ、(結果の)加重和ステップ、及び非線形関数ステップを含んでいてよい。畳み込みはその入力データについて、例えばその入力データにわたる移動ウィンドウ上のフィルタ(又はカーネル)を適用して特徴マップを生成することによって行われてよい。各層及び層の構成要素は、異なる所定のフィルタ(フィルタバンクからのもの)、重み(又は重み付けパラメータ)、及び/又は関数パラメータを有していてよい。この例において、入力データは、あるピクセル高さと幅の画像であり、この画像の生のピクセル値であってもよい。この例において、入力画像は3つの色チャネルRGB(赤、緑、青)の深さを有するように描かれている。任意選択により、入力画像には様々な前処理が行われてよく、前処理の結果が生の画像データの代わりに、又はそれに加えて入力されてもよい。画像処理の幾つかの例には、網膜血管マップセグメンテーション、色空間変換、適応型ヒストグラム均一化、接続構成要素生成等が含まれていてよい。ある層内で、ドット積がある重みとそれらが入力ボリューム内で接続された小さい領域との間で計算されてよい。CNNを構成するための多くの方法が当業界で知られているが、例として、層はゼロにおけるmax(0,x)閾値等、要素ごと活性化関数を適用するために構成されてもよい。プーリング関数は、ボリュームをダウンサンプルするために(例えばx-y方向に沿って)行われてもよい。完全に接続された層は、分類出力を特定し、画像認識及び分類に有益であることが判明している1次元出力ベクトルを生成するために使用されてよい。しかしながら、画像セグメンテーションのためには、CNNは各ピクセルを分類する必要がある。各CNN層は入力画像の解像度を低下させる傾向があるため、画像をその当初の解像度へとアップサンプルするための別のステージが必要である。これは、転置畳み込み(又は逆畳み込み)ステージTCの適用によって実現されてよく、これは典型的に、何れの所定の補間方法も使用せず、その代わりに学習可能パラメータを有する。
畳み込みニューラルネットワークは、コンピュータビジョンの多くの問題にうまく適用されている。前述のように、CNNを訓練するには一般に、大きな訓練データセットが必要である。U-NetアーキテクチャはCNNに基づいており、一般に従来のCNNより小さい訓練データセットで訓練できる。
図30は、例示的なU-Netアーキテクチャを図解する。この例示的なU-Netは、入力モジュール(又は入力層若しくはステージ)を含み、これは何れかのサイズ(例えば、128×128ピクセルのサイズ)の入力U-in(例えば、入力画像又は画像パッチ)を受け取る。入力画像は、眼底画像、OCT/OCTA en face、B-スキャン画像等であってよい。しかしながら、理解すべき点として、入力は何れの大きさ及び次元のものであってもよい。例えば、入力画像はRGBカラー画像、モノクロ画像、体積画像等であってよい。入力画像は一連の処理層を経て、その各々は例示的な大きさで図解されているが、これらの大きさは説明を目的としているにすぎず、例えば画像のサイズ、畳み込みフィルタ、及び/又はプーリングステージに依存するであろう。このアーキテクチャは、収束経路(4つの符号化モジュールを含む)とそれに続く拡張経路(4つの復号モジュールを含む)、及び対応するモジュール/ステージ間にあり、収束経路内の1つの符号化モジュールの出力をコピーして、それを拡張経路内の対応する復号モジュールの入力に結合する4つのコピー・アンド・クロップリンク(例えば、CC1~CC4)からなる。その結果、特徴的なU字型となり、そこからこのアーキテクチャが名付られている。収束経路はエンコーダと同様であり、その基本機能はコンパクトの特徴マップを介してコンテキストを捕捉することである。この例において、収束経路内の各符号化モジュールは2つの畳み込みニューラルネットワーク層を含み、それに続いて1つの最大プーリング層(例えば、ダウンサンプリング層)があってよい。例えば、入力画像U_inは2つの畳み込み層を経るが、各々が32の特徴マップを有する。特徴マップの数は各プーリングにおいて、第一のブロックの32の特徴マップから始まり、第二のブロックで64等々と、2倍になる。それゆえ、収束経路は、各々が畳み込みステージを提供する複数の符号化モジュール(又はステージ)と、それに続く活性化関数(例えば、正規化線形ユニットReLU、又はシグモイド層)及び最大プーリング演算からなる畳み込みネットワークを形成する。拡張経路はデコーダと同様であり、その機能は、収縮ステージで行われたダウンサンプリング及び何れの最大プーリングにもかかわらず、局所化を提供し、空間情報を保持することである。収束経路では、空間情報が縮小され、特徴情報は増大される。拡張経路は、複数の復号モジュールを含み、各復号モジュールは、その現在の値を対応する符号化モジュールの出力と結合する。すなわち、特徴及び空間情報は拡張経路においてアップコンボリューション(例えば、アップサンプリング又は転置畳み込み、すなわち逆畳み込み)と収束経路からの高解像度特徴との結合(例えば、CC1~CC4を介する)の連続を通じて組み合わされる。それゆえ、逆畳み込み層の出力は、収束経路からの対応する(任意選択によりクロップされた)特徴マップと、それに続いて2つの畳み込み層及び活性化関数(任意選択によるバッチ正規化)に結合される。拡張経路内の最後のモジュールからの出力は、分類器ブロック等、他の処理/訓練ブロック又は層に供給されてよく、これはU-Netアーキテクチャと共に訓練されてもよい。
収束経路と拡張経路との間のモジュール/ステージ(BN)は、「ボトルネック」と呼ばれてもよい。ボトルネックBNは、2つの畳み込み層(バッチ正規化と任意選択によるドロップアウトを伴う)からなっていてよい。
コンピューティングデバイス/システム
図31は、例示的なコンピュータシステム(又はコンピューティングデバイス又はコンピュータデバイス)を図解する。幾つかの実施形態において、1つ又は複数のコンピュータシステムは本明細書において記載又は図解された機能を提供し、及び/又は本明細書において記載又は図解された1つ又は複数の方法の1つ又は複数のステップを実行してよい。コンピュータシステムは、何れの適当な物理的形態をとってもよい。例えば、コンピュータシステムは、埋込みコンピュータシステム、システムオンチップ(SOC)、又はシングルボードコンピュータシステム(SBC)(例えば、コンピュータ・オン・モジュール(COM)又はシステム・オン・モジュール(SOM)等)、デスクトップコンピュータシステム、ラップトップ若しくはノートブックコンピュータシステム、コンピュータシステムのメッシュ、携帯電話、携帯型情報端末(PDA)、サーバ、タブレットコンピュータシステム、拡張/仮想現実装置、又はこれらのうちの2つ以上の組合せであってよい。適当であれば、コンピュータシステムはクラウド内にあってよく、これは1つ又は複数のクラウドコンポーネントを1つ又は複数のネットワーク内に含んでいてよい。
幾つかの実施形態において、コンピュータシステムはプロセッサCpnt1、メモリCpnt2、ストレージCpnt3、入力/出力(I/O)インタフェースCpnt4、通信インタフェースCpnt5、及びバスCpnt6を含んでいてよい。コンピュータシステムは、任意選択により、ディスプレイCpnt7、例えばコンピュータモニタ又はスクリーンも含んでいてよい。
プロセッサCpnt1は、コンピュータプログラムを構成するもの等、命令を実行するためのハードウェアを含む。例えば、プロセッサCpnt1は、中央処理ユニット(CPU)又は汎用コンピューティング・オン・グラフィクス処理ユニット(GPGPU)であってもよい。プロセッサCpnt1は、命令を内部レジスタ、内部キャッシュ、メモリCpnt2、又はストレージCpnt3から読み出し(又はフェッチし)、この命令を復号して実行し、1つ又は複数の結果を内部レジスタ、内部キャッシュ、メモリCpnt2、又はストレージCpnt3に書き込んでよい。特定の実施形態において、プロセッサCpnt1は、データ、命令、又はアドレスのための1つ又は複数の内部キャッシュを含んでいてよい。プロセッサCpnt1は、1つ又は複数の命令キャッシュ、1つ又は複数のデータキャッシュを、例えばデータテーブルを保持するために含んでいてよい。命令キャッシュ内の命令は、メモリCpnt2又はストレージCpnt3内の命令のコピーであってもよく、命令キャッシュはプロセッサCpnt1によるこれらの命令の読出しをスピードアップするかもしれない。プロセッサCpnt1は、何れの適当な数の内部レジスタを含んでいてもよく、1つ又は複数の算術論理演算ユニット(ALU:arithmetic logic units)を含んでいてよい。プロセッサCpnt1は、マルチコアプロセッサであるか、又は1つ若しくは複数のプロセッサCpnt1を含んでいてよい。本開示は特定のプロセッサを説明し、図解しているが、本開示は何れの適当なプロセッサも想定している。
メモリCpnt2は、処理を実行し、又は処理中に中間データを保持するプロセッサCpnt1のための命令を保存するメインメモリを含んでいてよい。例えば、コンピュータシステムは、命令又はデータ(例えば、データテーブル)をストレージCpnt3から、又は他のソース(例えば、他のコンピュータシステム)からメモリCpnt2にロードしてもよい。プロセッサCpnt1は、メモリCpnt2からの命令とデータを1つ又は複数の内部レジスタ又は内部キャッシュにロードしてもよい。命令を実行するために、プロセッサCpnt1は内部レジスタ又は内部キャッシュから命令を読み出して復号してもよい。命令の実行中又はその後に、プロセッサCpnt1は1つ又は複数の結果(これは、中間結果でも最終結果でもよい)を内部レジスタ、内部キャッシュ、メモリCpnt2、又はストレージCpnt3に書き込んでよい。バスCpnt6は、1つ又は複数のメモリバス(これは各々、アズレスバスとデータバスを含んでいてよい)を含んでいてよく、プロセッサCpnt1をメモリCpnt2及び/又はストレージCpnt3に連結してよい。任意選択により、1つ又は複数のメモリ管理ユニット(MMU)は、プロセッサCpnt1とメモリCpnt2との間のデータ伝送を容易にする。メモリCpnt2(これは、高速揮発性メモリであってもよい)には、ランダムアクセスメモリ(RAM)、例えばダイナミックRAM(DRAM)又はスタティックRAM(SRAM)が含まれていてよい。ストレージCpnt3には、データ又は命令のための長期又は大容量メストレージを含んでいてよい。ストレージCpnt3はコンピュータシステムに内蔵されても外付けでもよく、ディスクドライブ(例えば、ハードディスクドライブHDD、又はソリッドステートドライブSSD)、フラッシュメモリ、ROM、EPROM、光ディスク、磁気光ディスク、磁気テープ、ユニバーサルシリアルバス(USB)-アクセス可能ドライブ、又はその他の種類の不揮発性メモリのうちの1つ又は複数を含んでいてよい。
I/OインタフェースCpnt4は、ソフトウェア、ハードウェア、又はそれら両方の組合せであってよく、I/Oデバイスと通信するための1つ又は複数のインタフェース(例えば、シリアル又はパラレル通信ポート)を含んでいてよく、これはヒト(例えば、ユーザ)との通信を可能にしてもよい。例えば、I/Oデバイスとしては、キーボード、キーパッド、マイクロフォン、モニタ、マウス、プリンタ、スキャナ、スピーカ、スチールカメラ、スタイラス、テーブル、タッチスクリーン、トラックボール、ビデオカメラ、他の適当なI/Oデバイス、又はこれら2つ以上の組合せが含まれていてよい。
通信インタフェースCpnt5は、他のシステム又はネットワークと通信するためのネットワークインタフェースを提供してもよい。通信インタフェースCpnt5は、Bluetooth(登録商標)インタフェース又はその他の種類のパケットベース通信を含んでいてよい。例えば、通信インタフェースCpnt5は、ネットワークインタフェースコントローラ(NIC)及び/又は、無線ネットワークとの通信のための無線NIC若しくは無線アダプタを含んでいてよい。通信インタフェースCpnt5は、WI-FIネットワーク、アドホックネットワーク、パーソナルエリアネットワーク(PAN)、無線PAN(例えば、Bluetooth WPAN)、ローカルエリアネットワーク(LAN)、ワイドエリアネットワーク(WAN)、メトロポリタンエリアネットワーク(MAN)、携帯電話ネットワーク(例えば、汎欧州デジタル移動電話方式(Global System for Mobile Communications)(GSM(登録商標))ネットワーク等)、インタネット、又はこれらの2つ以上の組合せとの通信を提供してよい。
バスCpnt6は、コンピューティングシステムの上述のコンポーネント間の通信リンクを提供してよい。例えば、バスCpnt6は、アクセラレーテッド・グラフィックス・ポート(Accelerated Graphics Port)(AGP)若しくはその他のグラフィクスバス、拡張業界標準(Enhanced Industry Standard)アーキテクチャ(EISA)バス、フロントサイドバス(FSB)、ハイパートランスポート(HyperTransport)(HT)インタコネクト、業界標準アーキテクチャ(ISA)バス、インフィニバンド(InfiniBand)バス、low-pin-count(LPC)バス、メモリバス、マイクロチャネルアーキテクチャ(MCA)バス、ペリフェラル・コンポーネント・インターコネクト(Peripheral Component Interconnect)(PCI)バス、PCI-Express(PCIe)バス、シリアル・アドバンスト・テクノロジ・アタッチメント(serial advanced technology attachment)(SATA)バス、ビデオ・エレクトロニクス・スタンダーズ・アソシエーション・ローカル(Video Electronics Standards Association local)(VLB)バス、若しくはその他の適当なバス、又はこれらの2つ以上の組合せを含んでいてよい。
本開示は、特定の数の特定のコンポーネントを特定の配置で有する特定のコンピュータシステムを説明し、図解しているが、本開示は何れの適当な数の何れの適当なコンポーネントを何れの適当な配置で有する何れの適当なコンピュータシステムも想定している。
本明細書において、コンピュータ可読非一時的記憶媒体は、1つ又は複数の半導体ベース又はその他の集積回路(IC)(例えば、フィールドプログラマブルゲートアレイ(FPGA)若しくは特定用途IC(ASIC))、ハードディスクドライブ(HDD)、ハイブリッドハードドライブ(HHD)、光ディスク、光ディスクドライブ(ODD)、磁気光ディスク、磁気光ドライブ、フロッピディスケット、フロッピディスクドライブ(FDD)、磁気テープ、ソリッドステートドライブ(SSD)、RAM-ドライブ、SECURE DIGITALカード若しくはドライブ、その他のあらゆる適当なコンピュータ可読非一時的記憶媒体、又は適当であればこれらの2つ以上あらゆる適当な組合せを含んでいてよい。コンピュータ可読非一時的記憶媒体は、揮発性、不揮発性、又は適当であれば揮発性と不揮発性の組合せであってよい。
本発明は幾つかの具体的な実施形態と共に説明されているが、当業者にとっては明白であるように、上記の説明を参照すれば多くのその他の代替案、改良、及び変形型が明らかである。それゆえ、本明細書に記載の発明は、付属の特許請求の範囲の主旨と範囲に含まれるかもしれないあらゆるこのような代替案、改良、応用、及び変形型の全てを包含することが意図されている。

Claims (32)

  1. 眼科撮像システムの焦点ずれ測定値を決定するための方法であって、
    眼の網膜の第1の領域を照明することと、
    前記第1の領域から検出器上に戻る光の第1の画像の収集することと、
    前記眼の網膜の第2の領域を照明することであって、前記第2の領域は前記第1の領域に部分的に重なる、照明することと、
    前記第2の領域から前記検出器上に戻る光の第2の画像の収集することと、
    前記第1の画像と前記第2の画像との差を求めることにより第3の画像を定義することであって、前記第3の画像は線形領域を含む、定義することと、
    前記線形領域の幅測定値を決定することと、
    前記幅測定値に基づいて前記焦点ずれ測定値を決定することと
    を含む、方法。
  2. 前記幅測定値は、前記第3の画像の強度分布の2次モーメントに基づく、請求項1に記載の方法。
  3. 眼科撮像システムを使用して眼のトポグラフィを特徴付ける方法であって、
    前記眼の網膜の領域を照明することと、
    網膜の前記領域から検出器上に戻る光を集光することと、
    前記検出器上の予め決定された位置に対する前記検出器上の集光された光の位置のオフセットを決定することであって、前記オフセットは網膜の高さに対応する、決定することと、
    決定された前記オフセット又は前記オフセットの更なる分析を記憶又は表示することと
    を含む、方法。
  4. 被照明領域は線形領域であり、
    前記集光された光は、集光された長尺領域であり、
    前記オフセットを決定することは、前記長尺領域の長さに沿った強度重心を前記検出器上の予め決定された線形位置と比較してトポロジ線を定義することを含む、請求項3に記載の方法。
  5. 網膜の領域に及ぶ複数の前記トポロジ線を収集することと、
    前記複数のトポロジ線に基づいて網膜の前記領域のトポロジを定義することと
    を更に含む、請求項4に記載の方法。
  6. 眼科撮像システムの焦点ずれ測定値を決定するための方法であって、
    眼の網膜の第1の領域に正弦波状照明を照射することと、
    前記第1の領域から検出器上に戻る光の第1の画像の収集することと、
    前記第1の画像に沿った複数の点における焦点を決定することにより焦点ずれマップを決定することと、
    前記焦点ずれマップに基づいて前記焦点ずれ測定値を決定することと
    を含む、方法。
  7. 前記第1の画像は、強度を表すアルファ(α)スケーリング係数により特徴付けられ、
    前記第1の画像に沿った複数の点における焦点を決定することは、前記アルファ(α)スケーリング係数を除去することと、強度重心の位置を決定することとを含む、請求項6に記載の方法。
  8. 前記アルファ(α)スケーリング係数は、強度正規化により除去される、請求項7に記載の方法。
  9. 照射された正弦波状照明は、複数の異なる周波数にわたって掃引され、
    前記第1の画像は、強度を表すアルファ(α)スケーリング係数により特徴付けられ、前記アルファ(α)スケーリング係数は、掃引された周波数に依存し、
    前記アルファ(α)スケーリング係数を収集してアルファマップが作成され、
    前記アルファマップに基づいて前記焦点ずれマップが生成される、請求項6に記載の方法。
  10. 前記アルファ(α)スケーリング係数に基づいて前記システムの点広がり関数を決定することを含む、請求項9に記載の方法。
  11. 標的眼底画像内の視神経乳頭(ONH)の位置を特定する方法であって、
    標的画像の取り込み直前又は直後に第1の撮像モダリティの第1の眼底画像を取り込むことであって、前記標的画像は第2の撮像モダリティによるものである、取り込むことと、
    前記第1の眼底画像内の複数の可能性のあるONH領域を特定する深層学習ニューラルネットワークに前記第1の眼底画像を送ることと、
    所定のサイズ範囲内にある、ニューラルネットワークにより特定される、任意の可能性のあるONH領域を候補ONH領域として特定することと、
    最も均一に丸みを帯びた前記候補ONH領域を前記第1の眼底画像内の真のONH領域として指定することと、
    前記第1の眼底画像上の前記真のONH領域に対応する前記標的眼底画像上の領域を前記標的眼底画像内の視神経乳頭として特定することと
    を含む、方法。
  12. 前記標的眼底画像は、蛍光眼底血管造影(FA)画像又はインドシアニングリーン血管造影(ICGA)画像の一方であり、
    前記第1の眼底画像は、人間の眼では実質的に知覚できない光を使用して取り込まれる、請求項11に記載の方法。
  13. 前記第1の眼底画像は赤外線画像であり、
    前記深層学習ニューラルネットワークは、カラー画像と赤外線画像との組み合わせで訓練されたUネットである、請求項11に記載の方法。
  14. 前記Uネットは、カラー画像のみを使用して前記Uネットが訓練される第1のステップと、前記Uネットの所定数の初期ネットワーク層が固定され且つ赤外線画像のみを使用して残りのネットワーク層が訓練される第2のステップとを含む、2ステップ工程を使用して訓練される、請求項13に記載の方法。
  15. 血管造影画像を表示する方法であって、
    前記血管造影画像のシーケンスを取得することと、
    取得された各血管造影画像についての輝度測定値を計算することと、
    計算された最も高い輝度測定値に基づいて、取得された各血管造影画像についての輝度スケーリング係数を決定することと、
    各血管造影画像の決定された輝度スケーリング係数に従って、取得された各血管造影画像を表示のために修正することと、
    修正された画像の前記シーケンスを表示することと
    を含む、方法。
  16. 各血管造影画像の前記輝度スケーリング係数は、表示のために修正される前と修正された後で取得されたシーケンスの前記血管造影画像間の相対的な輝度が実質的に同じままであるように決定される、請求項15に記載の方法。
  17. 決定された前記輝度スケーリング係数は、血管造影画像の前記シーケンスにおける全ての画像に対して同じである、請求項15に記載の方法。
  18. 各血管造影画像についての決定された前記輝度スケーリング係数は、前記各血管造影画像の個々の白色レベル、計算された最も高い輝度測定値、及び輝度圧縮係数に基づく、請求項15に記載の方法。
  19. 血管造影画像の前記シーケンスは、血管造影画像の前記シーケンスにおける任意の画像よりも大きい撮像ダイナミックレンジを有する眼科撮像装置を使用して取得される、請求項15に記載の方法。
  20. 血管造影画像の前記シーケンスは、前記眼科撮像装置の利得を調節せずに取得される、請求項19に記載の方法。
  21. 前記シーケンスにおける各血管造影画像が取得されたときに、
    現在取得された血管造影画像の生データを記憶することと、
    前記現在取得された血管造影画像の白色レベルを決定して、前記現在取得された血管造影画像に関連付けられるメタデータに前記白色レベルを記憶することと、
    前記現在取得された血管造影画像の白色レベルに基づいて、前記現在取得された血管造影画像に初期輝度係数を適用して、仮画像を定義することと、
    前記仮画像を表示することと
    を含む、請求項15に記載の方法。
  22. 血管造影画像の前記シーケンス全体にわたる強度の時間経過を記憶又は表示することを更に含む、請求項15に記載の方法。
  23. 眼科撮像装置のための焦点合わせ方法であって、
    眼の検査画像を取り込むことと、
    電子ディスプレイ上に前記検査画像を表示して、前記電子ディスプレイ内の標的領域を示すユーザ入力を受け取ることと、
    前記検査画像内の複数の基準領域の各々についての個々の焦点ずれ測定値を決定することと、
    前記標的領域に対する前記基準領域の位置に基づいて基準焦点ずれ測定値を組み合わせることにより前記標的領域についての標的焦点ずれ測定値を定めることと
    を含む、焦点合わせ方法。
  24. 前記複数の基準領域は、前記電子ディスプレイ上で特定され、
    前記標的領域は、前記複数の基準領域の中から選択される、請求項23に記載の焦点合わせ方法。
  25. 前記標的焦点ずれ測定値は、前記標的領域と選択された基準領域との間の距離により重み付けされた組み合わせに基づく、請求項23に記載の焦点合わせ方法。
  26. 各基準領域についての前記焦点ずれ測定値は、前記眼の固視角度に少なくとも部分的に基づく、請求項23に記載の焦点合わせ方法。
  27. 前記検査画像は、赤外線プレビュー画像であり、
    前記固視角度は、前記赤外線プレビュー画像内の前記眼の視神経乳頭を特定することにより決定される、請求項26に記載の焦点合わせ方法。
  28. 眼底画像を取り込む方法であって、
    画面上にグラフィカルユーザインタフェース(GUI)を表示することであって、前記GUIは、複数の撮像型オプションを提供し、各撮像型オプションは、撮像パラメータの個々の既定のセットに関連付けられる、表示することと、
    電子プロセッサを使用して、患者が所定の医学的条件を満たすことに応答して、拡張機能オプションを選択的に設定することと、
    拡張機能が設定されたときに画像取り込み入力指示に応答して、
    現在選択された前記撮像型オプションの第1の撮像型オプションを特定することと、
    第1の撮像型オプションに関連付けられる撮像パラメータの前記既定のセットを前記拡張機能オプションにより指定された調整された撮像パラメータのセットに置き換えることと、
    画像取り込みシーケンスを実行して、撮像パラメータの調整された前記セットを使用して眼底の画像を取り込むことと
    を含む、方法。
  29. 前記医学的条件は、羞明、心的外傷後ストレス障害、及び所定のDNAマーカのうちの1つ又は複数である、請求項28に記載の方法。
  30. 撮像パラメータの前記既定のセットは、画像取り込み動作中に照射される光の強度を制御する光強度パラメータ及び前記画像取り込み動作中に光が照射される持続時間を制御する光持続時間パラメータのうちの1つ又は複数を含み、
    調整された前記撮像パラメータは、前記光強度パラメータ及び前記光持続時間パラメータのうちの一方又は両方を低下させる、請求項28に記載の方法。
  31. 調整された撮像パラメータは、眼の虹彩サイズ、虹彩色、及び網膜明度レベルのうちの1つ又は複数を含む、前記眼の生理学的特性に基づく、請求項28に記載の方法。
  32. 眼底画像を取り込む方法であって、
    スリット走査式検眼鏡を使用して、眼底画像取り込みシーケンスを開始することであって、
    眼底上の複数の走査位置を指定する走査テーブルを評価することであって、各走査位置は複数の別個の光波範囲に関連付けられる、評価することと、
    前記走査テーブルにより指定される眼底の走査シーケンスを開始することであって、
    (a)少なくとも1つの光源を使用して、現在の走査位置に関連付けられる別個の各光波範囲を順番に適用することと、
    (b)少なくとも1つの集光器を使用して、順番に適用される光波範囲ごとに前記現在の走査位置からの戻り光を集光することと
    を含む、開始することと、
    各光波範囲のそれぞれ集光された戻り光に基づいて光波範囲ごとに別々の眼底画像を生成することと、
    個々の眼底画像の所定の領域内における前記別個の光波範囲の少なくとも2つからの集光された戻り光のレシオメトリックを決定することと、
    前記眼底画像を処理又は表示することと
    を含む、開始することを含み、
    前記レシオメトリックは、黄斑色素密度及び酸素化のうちの少なくとも1つの測定値である、方法。
JP2021556553A 2019-03-20 2020-03-18 単一露光多種撮像と、改善された焦点合わせと、改善された血管造影画像シーケンス表示とを備えた、患者に合わせて調整される眼科撮像システム Pending JP2022525674A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201962821283P 2019-03-20 2019-03-20
US62/821,283 2019-03-20
PCT/EP2020/057522 WO2020188007A1 (en) 2019-03-20 2020-03-18 A patient tuned ophthalmic imaging system with single exposure multi-type imaging, improved focusing, and improved angiography image sequence display

Publications (1)

Publication Number Publication Date
JP2022525674A true JP2022525674A (ja) 2022-05-18

Family

ID=69941371

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021556553A Pending JP2022525674A (ja) 2019-03-20 2020-03-18 単一露光多種撮像と、改善された焦点合わせと、改善された血管造影画像シーケンス表示とを備えた、患者に合わせて調整される眼科撮像システム

Country Status (5)

Country Link
US (1) US20220160228A1 (ja)
EP (1) EP3941333A1 (ja)
JP (1) JP2022525674A (ja)
CN (1) CN113613545A (ja)
WO (1) WO2020188007A1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022087469A1 (en) * 2020-10-23 2022-04-28 Genentech, Inc. Multimodal geographic atrophy lesion segmentation
CN112669256B (zh) * 2020-11-23 2024-02-20 西安电子科技大学 一种基于迁移学习的医学图像分割与显示方法
EP4256513A1 (en) * 2020-12-03 2023-10-11 Genentech, Inc. Multimodal prediction of geographic atrophy growth rate
WO2022190113A1 (en) * 2021-03-10 2022-09-15 Forus Health Pvt. Ltd. A portable retinal imaging system
US20230015951A1 (en) * 2021-07-12 2023-01-19 Welch Allyn, Inc. Multispectral fundus imaging
CN113724214B (zh) * 2021-08-23 2024-02-23 唯智医疗科技(佛山)有限公司 一种基于神经网络的图像处理方法及装置
CO2022011610A1 (es) * 2022-08-17 2024-02-26 Univ Simon Bolivar Método implementado por aprendizaje de máquina para la detección de vasos sanguíneos a través de clustering de píxeles
CN116416417B (zh) * 2023-06-09 2023-09-05 远洋装饰工程股份有限公司 基于增强现实的卫生间隐藏式排水结构辅助检修系统
CN117058132B (zh) * 2023-10-11 2024-01-23 天津大学 基于神经网络的文物照明视觉舒适度量化评估方法及系统
CN117351009B (zh) * 2023-12-04 2024-02-23 江苏富翰医疗产业发展有限公司 基于多光谱眼底图像生成血氧饱和度数据的方法及系统

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4170398A (en) 1978-05-03 1979-10-09 Koester Charles J Scanning microscopic apparatus with three synchronously rotating reflecting surfaces
US4732466A (en) 1985-04-04 1988-03-22 Humphrey Instruments, Inc. Fundus camera
US6549801B1 (en) 1998-06-11 2003-04-15 The Regents Of The University Of California Phase-resolved optical coherence tomography and optical doppler tomography for imaging fluid flow in tissue with fast scanning speed and high velocity sensitivity
US6741359B2 (en) 2002-05-22 2004-05-25 Carl Zeiss Meditec, Inc. Optical coherence tomography optical scanner
US7359062B2 (en) 2003-12-09 2008-04-15 The Regents Of The University Of California High speed spectral domain functional optical coherence tomography and optical doppler tomography for in vivo blood flow dynamics and tissue structure
US7301644B2 (en) 2004-12-02 2007-11-27 University Of Miami Enhanced optical coherence tomography for anatomical mapping
US7365856B2 (en) 2005-01-21 2008-04-29 Carl Zeiss Meditec, Inc. Method of motion correction in optical coherence tomography imaging
US8488895B2 (en) 2006-05-31 2013-07-16 Indiana University Research And Technology Corp. Laser scanning digital camera with pupil periphery illumination and potential for multiply scattered light imaging
AU2007300310B2 (en) 2006-09-26 2011-12-08 Oregon Health & Science University In vivo structural and flow imaging
JP5341386B2 (ja) * 2008-04-24 2013-11-13 キヤノン株式会社 眼科撮影装置
JP5581690B2 (ja) * 2009-12-25 2014-09-03 ソニー株式会社 厚み情報取得装置、厚み情報取得方法、厚み情報取得プログラム及び顕微鏡
CN103002794B (zh) 2010-02-08 2016-08-17 奥勒冈保健科学大学 用于超高灵敏光学微血管显影术之设备及方法
DE102010050693A1 (de) 2010-11-06 2012-05-10 Carl Zeiss Meditec Ag Funduskamera mit streifenförmiger Pupillenteilung und Verfahren zur Aufzeichnung von Fundusaufnahmen
US8433393B2 (en) 2011-07-07 2013-04-30 Carl Zeiss Meditec, Inc. Inter-frame complex OCT data analysis techniques
US8998411B2 (en) 2011-07-08 2015-04-07 Carl Zeiss Meditec, Inc. Light field camera for fundus photography
DE102011053880B4 (de) 2011-09-23 2023-11-09 Carl Zeiss Ag Vorrichtung und Verfahren zum Abbilden eines Augenhintergrunds
US9332902B2 (en) 2012-01-20 2016-05-10 Carl Zeiss Meditec, Inc. Line-field holoscopy
US9456746B2 (en) 2013-03-15 2016-10-04 Carl Zeiss Meditec, Inc. Systems and methods for broad line fundus imaging
US9759544B2 (en) 2014-08-08 2017-09-12 Carl Zeiss Meditec, Inc. Methods of reducing motion artifacts for optical coherence tomography angiography
US9700206B2 (en) 2015-02-05 2017-07-11 Carl Zeiss Meditec, Inc. Acquistion and analysis techniques for improved outcomes in optical coherence tomography angiography
JP6822965B2 (ja) 2015-02-05 2021-01-27 カール ツアイス メディテック アクチエンゲゼルシャフト ブロードライン眼底イメージングにおいて散乱光を低減するための方法および装置
WO2018178269A1 (en) 2017-03-31 2018-10-04 Carl Zeiss Meditec, Inc. Systems and methods for characterizing refraction with ophthalmic imaging systems

Also Published As

Publication number Publication date
US20220160228A1 (en) 2022-05-26
WO2020188007A1 (en) 2020-09-24
EP3941333A1 (en) 2022-01-26
CN113613545A (zh) 2021-11-05

Similar Documents

Publication Publication Date Title
JP2022525674A (ja) 単一露光多種撮像と、改善された焦点合わせと、改善された血管造影画像シーケンス表示とを備えた、患者に合わせて調整される眼科撮像システム
US11935241B2 (en) Image processing apparatus, image processing method and computer-readable medium for improving image quality
US20200394789A1 (en) Oct-based retinal artery/vein classification
US20220084210A1 (en) Segmentation and classification of geographic atrophy patterns in patients with age related macular degeneration in widefield autofluorescence images
US20220400943A1 (en) Machine learning methods for creating structure-derived visual field priors
US10102621B2 (en) Apparatus, method, and program for processing image
US20210398259A1 (en) Image processing apparatus, image processing method, and storage medium
US20220151483A1 (en) Ophthalmic apparatus, method for controlling ophthalmic apparatus, and computer-readable medium
Trucco et al. Computational retinal image analysis: tools, applications and perspectives
JP7348374B2 (ja) 眼科情報処理装置、眼科撮影装置、眼科情報処理方法、及びプログラム
JP2021037239A (ja) 領域分類方法
JP2022155690A (ja) 画像処理装置、画像処理方法、及びプログラム
JP7362403B2 (ja) 画像処理装置および画像処理方法
JP2020039430A (ja) 画像処理装置、画像処理方法及びプログラム
US20230140881A1 (en) Oct en face pathology segmentation using channel-coded slabs
WO2022145129A1 (ja) 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム
US20220390369A1 (en) Systems And Methods For Imaging And Characterizing Objects Including The Eye Using Non-Uniform Or Speckle Illumination Patterns
JP7344847B2 (ja) 画像処理装置、画像処理方法、及びプログラム
JP7478216B2 (ja) 眼科装置、眼科装置の制御方法、及びプログラム
US20230143051A1 (en) Real-time ir fundus image tracking in the presence of artifacts using a reference landmark

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230308

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240220