CN101778593B - 用于分析光学相干断层照相图像的方法 - Google Patents

用于分析光学相干断层照相图像的方法 Download PDF

Info

Publication number
CN101778593B
CN101778593B CN2008801033823A CN200880103382A CN101778593B CN 101778593 B CN101778593 B CN 101778593B CN 2008801033823 A CN2008801033823 A CN 2008801033823A CN 200880103382 A CN200880103382 A CN 200880103382A CN 101778593 B CN101778593 B CN 101778593B
Authority
CN
China
Prior art keywords
pattern
oct
flv
thickness
deviation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN2008801033823A
Other languages
English (en)
Other versions
CN101778593A (zh
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.)
University of Southern California USC
Original Assignee
University of Southern California USC
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 University of Southern California USC filed Critical University of Southern California USC
Publication of CN101778593A publication Critical patent/CN101778593A/zh
Application granted granted Critical
Publication of CN101778593B publication Critical patent/CN101778593B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/02—Subjective types, i.e. testing apparatus requiring the active assistance of the patient
    • A61B3/024—Subjective types, i.e. testing apparatus requiring the active assistance of the patient for determining the visual field, e.g. perimeter types
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062—Arrangements for scanning
    • A61B5/0066—Optical coherence imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30041—Eye; Retina; Ophthalmic
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00—Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03—Recognition of patterns in medical or anatomical images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00—Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03—Recognition of patterns in medical or anatomical images
    • G06V2201/031—Recognition of patterns in medical or anatomical images of internal organs

Landscapes

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

Abstract

本发明公开了用于分析视网膜断层照相图以检测诸如青光眼、视神经炎、前部缺血性视神经病的视神经疾病的图案的方法。成图的区域包括以凹窝为中心的黄斑，以及以视神经乳头为中心的区域。所分析的视网膜层包括神经纤维、神经节细胞、内丛状层和内核层以及它们的组合。还可以分析总视网膜厚度。对这些图应用图案分析，以创建用于对青光眼和视神经炎进行诊断和发展分析的单个参数。

Description

用于分析光学相干断层照相图像的方法
本发明至少部分是在美国国家卫生研究院（NIH）授权码R01EY013516的支持下做出的。美国政府具有本发明的某些权利。
相关申请的相交引用
本申请要求2007年6月15日提交的、申请号60/944,449、发明名称为“通过光学相干分析的用于诊断视神经疾病的视网膜地形图图案分析”的美国临时申请的优先权。在此要求35USC§119（e）下的美国临时申请的优先权。通过引用将上述在线申请并入到本文中。
技术领域
本发明属于眼科领域，更具体地说，本发明属于用于获取并且分析光学相干断层照相图像以检测视神经疾病的方法。
背景技术
所有的视神经病变都会在根本上影响视网膜的内层。尤其，神经纤维层（NFL）、神经节细胞层（GCL）和内丛状层（IPL）是最受影响的。相反，内核层（INL）是受到影响较小的，而外层是不受影响的【1】。由于视网膜中三个内层各自包含视网膜神经节细胞的不同部分（NFL包含轴突，GCL包含细胞体，IPL包含树突），因此将该局部解剖变异考虑进去的诊断方法通常具有更好的诊断特异性。例如，尽管对视网膜总厚度的测量会提供一般的诊断信息，但是聚焦在视神经乳头（ONH，还称为视神经盘）附近区域周围的测量将提供更多的诊断信息，因为NFL在这个区域是最厚的。在黄斑（凹窝周围的区域）中，所有三个内层都提供诊断信息，因此，最好测量包括NFL，GCL和IPL的组合内视网膜层。
虽然这个原理在理论上是很简单的，但是其在实践上却是不容易的。以青光眼的诊断为例。青光眼视神经病变是由于眼科解剖学中的几处步进性交替而引起的：视网膜神经节细胞（RGC）的丢失，视网膜神经纤维层（NFL）的变薄以及视神经盘成杯状。因此，顺理成章，这些解剖学变化可以用作青光眼的诊断指示符。不幸的是，实践中，在诊断方法中不能容易地运用该知识。在传统的裂隙灯眼科检查中无法看到RGC丢失。同样，在临床检查中很难检测到NFL束缺陷。尽管无赤光眼底照相术能够检测到视网膜的血管系统以及神经纤维中的变化，但是该技术很少用在临床实践中。因此，目前，青光眼的临床诊断仅仅基于特征性的视神经杯状并且结合针对患者对应视场缺损的测试。
但是，由于RGC群体的大量丢失可能在可以检测到视场缺损之前发生，并且该结构性丢失可能比可检测到功能丧失早发生长达5年，所以用于青光眼临床诊断的现有方法不足于对疾病进行早期检测。因此，目前，对于能够识别并且量化RGC群体的变化并且在临床设置中也容易管理的检测和预测方法的需求尚未满足。
附图说明
图1示出了根据本发明实施方式的黄斑图7毫米（MM7）扫描图案。该图中示出的示例性扫描在0.58秒内在7毫米的方形区域中扫描了14,944点。
图2示出了黄斑区域中的视网膜横截面的示例性OCT图像。神经节细胞复合体（GCC）也称为内视网膜层（IRL），由神经纤维层（NFL）、神经节细胞层（GCL）和内丛状层（IPL）组成。视网膜在凹窝中比较薄，这充当用于定位凹窝中心的界标。
图3示出了根据本发明实施方式的示例性GCC厚度图。单位是μm。从分析中去除了无法可靠测量GCC的中央凹窝区域（黑圆圈）。
图4示出了根据本发明实施方式的一组示例性GCC厚度图和派生图。左上：被测试的眼睛的GCC厚度图，单位：μm。右上：正常平均GCC厚度图，单位：μm。左下：偏差（D）图，单位：μm。右下：分数偏差（FD）图，单位：%。
图5示出了具有统计上显著的焦点丢失的区域（由阴影线图案标记）的青光眼眼睛的示例性GCC FD图。单位是%。
图6示出了一组示例性GCC厚度图和派生图。左上：被测试的眼睛的GCC厚度图，单位：μm。右上：测试眼睛的图案图，单位%。左下：平均正常眼睛的图案图，单位：%。右下：图案偏差（PD）图，单位%。
图7示出了具有统计上显著的焦点变稀疏的区域（由阴影线图案标记）的青光眼眼睛的示例性GCC FD图。单位是%。
图8示出了对于各种视神经病变的示例性特征PD图组。左上：下方青光眼（inferior glaucoma）（IG），35只眼，右上：均匀青光眼（EG），41只眼。左中：上方青光眼（superior glaucoma）（SG），18只眼，右中：视神经炎（ON），22只眼。左下：上方AION（SAION），11只眼。右下：下方AION（IANON），7只眼。单位是%。
图9示出了从青光眼眼睛获得的黄斑中的视网膜总厚度的示例性偏差图。单位是μm。
图10示出了从青光眼眼睛获得的总和GCC/ORC反射率比的示例性分数偏差图。单位是%。
图11示出了示例性正常乳头周围的NFL厚度图。
图12示出了用于视野计（perimetric）青光眼情况的示例性图组。所有的黄斑GCC厚度（mGCCT）参数都是不正常的（平均值=71μm；p<0.5%；焦点丢失量=12.6%，p<0.5%，总体丢失量=26.5%，p<0.5%；图案方差系数=21%，p<0.5%；上方-下方差=17.0μm，p<0.5%）。（A）示出了示例性mGCCT图。（B）示出了具有由阴影线图案标记的显著焦点丢失的区域的示例性分数偏差图。（C）示出了示例性图案偏差图。（D）示出了表示颞下边缘丢失的示例性神经盘照片。（E）示出了示例性视场（VF）图案偏差（PD）图。VF是不正常的（图案标准偏差=16.5dB，p<0.5%）；青光眼半场测试是在正常限度之外的。椭圆形虚线表示对应于mGCCT图的区域。上VF缺陷对应于下神经节细胞丢失以及神经盘边缘变薄。
图13示出了PG眼的平均mGCCT分数偏差图的示例性图像。这表示青光眼中的神经节细胞丢失的特征图案。
图14示出了来自示例性PPG情况示例的图像。（A）：mGCCT FD图。一些mGCCT参数是不正常的（AVG=82.5μm，p>5%；FLV=4.9%，p<0.5%，GLV=13.2%，p<0.5%；PCV=0.13，P>5%；SID=-12.1μm，p<0.5%）。（B）示出了颞上边缘的早期轻微变薄的神经盘照片。（C）VF PD图。
发明内容
简略地说，本发明提供一种用于检测并诊断受验者的视神经病变的方法，该方法一般包括以下步骤：利用FD-OCT生成黄斑图；基于黄斑图构造黄斑区域的三维厚度图；从所述厚度图计算派生图；通过将图案分析方法应用于图来识别图中的不正常区域；并且基于所述厚度图、所述派生图、在图中识别出的不正常区域、或其组合，来确定诊断参数。
在本发明的方法中，发明人已经发现，黄斑图的质量和厚度图的后续构造很大程度上依靠采用的扫描图案。因此，本发明的发明人已经设计了新的扫描图案以实现快速黄斑扫描，这可以有利于构造高质量三维厚度图。
一旦通过采用根据本发明实施方式的扫描图案获得了黄斑图，然后就可以通过在黄斑图的个体扫描之间进行插值来构造黄斑区域的二维厚度图。根据黄斑区域的该二维图，然后可以构造各种派生图，并且图案分析方法用于提取对检测和诊断视神经病变有用的诊断信息。
在优选实施方式中，通过根据如图1所示的MM7扫描图案来执行一系列横截面扫描，从而首先获取受验者的黄斑区域的FD-OCT图像。接着，对这第一组原始横截面扫描进行处理和插值以构造患者的黄斑区域的三维模型。根据该三维模型，接着可以计算视网膜厚度图和其他派生图。优选的是，根据厚度图计算分数偏差图，并且识别不正常视网膜厚度的区域。利用这些图，接着计算诊断参数。优选的是，如下面的详细描述所限定的，计算参数焦点丢失量（focal loss volume，FLV）和总体丢失量（global loss volume，GLV）并且将其用于辅助诊断的确定。
虽然上述描述的实施方式概述了本发明的一般步骤，但是本领域技术人员将会理解各种修改例都是可能的。根据下面的详细描述和所附权利要求书，本发明的其他方面和优势会变得显而易见。
具体实施方式
如上所述，本发明描述用于获取黄斑的诊断图像以及对图像进行后续分析以产生对做出诊断有用的诊断参数的方法。具体来说，本发明提供利用诸如FD-OCT的非侵入成像技术来产生黄斑图（图像）的方法，以及用于处理并分析所生成的图（图像）的方法。通过以高分辨率测量并且监视黄斑区域周围的解剖结构中的变化，即使在疾病症状表现之前，也可以在早期阶段检测到病理状态。本发明的方法可以应用到影响黄斑区域的所有类型的视神经病变。
通常，根据本发明实施方式的视神经病变检测和诊断方法具有下述步骤：利用傅氏域光学相干断层照相（FD-OCT）来生成以凹窝（黄斑图）或视神经乳头（乳头周围图）为中心的初始图，其中，所述初始图由通过多个横截面扫描所形成的扫描图案组成；基于初始图来构造视网膜特性图；从视网膜特性图来计算派生图；将图案分析方法应用到视网膜特性图或派生图以检测不正常区域；以及基于图、派生图、检测出的不正常区域或其组合来确定诊断参数，其中，所述诊断参数、厚度图和派生图各自可以单独使用或者组合使用以诊断并区分不同类型的视神经病变。
用于生成黄斑图的扫描图案优选是覆盖了黄斑的广大区域的图。所得到的黄斑图应当优选地具有至少10μm的分辨率。扫描应当优选地在大约2秒内完成。在一个优选实施方式中，扫描图案是如图1所示的MM7图案。
从由黄斑的多个横截面图像构成的原始黄斑图来构造厚度图。通过对准各个相邻横截面扫描并且在横截面之间进行插值，可以从黄斑图中的二维轴扫描的集合获得黄斑的三维图像。通过识别图像中的结构边界并且测量解剖结构的距离和大小，可以指定视网膜层的厚度以得到厚度图。取决于所关注的结构，可以构造各种厚度图，包括视网膜总厚度、神经节细胞复合体（或内视网膜层）厚度、神经纤维层厚度、神经节细胞层厚度、内丛状层厚度，但是并不限于此。
优选地通过软件来自动构造厚度图，但是这也可以手动完成。
可以将特性图进一步限制于视网膜的特定空间位置，例如乳头周围的区域（视神经盘附近的区域）。
一旦构造出特性图，就可以通过对特性图应用数据变换算法来获得多种多样的派生图。示例性派生图可以包括偏差图、分数偏差图、图案偏差图、或其组合，但并不限于此。
为了检测图中的不正常区域，可以对特性图或派生图应用各种图案分析算法。图案分析算法优选地是基于统计的算法。本发明的示例性统计图案分析算法可以包括标准偏差比较、总平均、上方平均、下方平均、或其组合，但并不限于此。对图案分析算法的选择将依赖于分析对象。
一旦已经识别了不正常区域，就可以根据图和对不正常区域的了解来定义并且计算诊断参数。示例性诊断参数可以包括焦点丢失量（FLV）、总体丢失量（GLV）、图案方差系数（PCV）或（青光眼）图案关联度（GPCC），但并不限于此。
在做出诊断确定时，单独的诊断参数或者诊断参数与其他参数结合都是有用的。可应用到本发明方法的视神经病变的示例可以包括青光眼、视神经炎和前部缺血性视神经病（AION），但并不限于此。
当将图与某些视神经病变的其他参照图特征相互关联时，可以计算出关联参数。测试受验者图与疾病参照图之间的关联度越高，就表示受验者患有疾病的可能性更高。
为了进一步说明本发明的工作原理和益处，我们将首先考虑青光眼的诊断。
如背景技术中所讨论的，青光眼的一个定义特征是RGC的丢失。本领域中公知的是，RGC群体的很大比例处于黄斑中，因此，黄斑厚度提供了用于检测和预测青光眼的有用的诊断措施。Zeimer等人（参考文献6）首先描述了利用裂隙扫描视网膜厚度分析器（RTA，Tailia技术有限责任公司，Neve-Ilan，以色列）来检测青光眼中的黄斑厚度减小。由于本发明的发明者之一（Huang）及其合作者引入了光学相干断层照相（OCT），因此该技术已经被证实为可用于测量乳头周围的神经纤维层厚度（cpNELT），所述乳头周围的神经纤维层厚度被示出为检测青光眼的有用参数。但是，利用OCT的黄斑视网膜总厚度（mRT）测量还不是像cpNFLT那样精确的诊断参数。早期视网膜OCT系统采用了较慢的时域技术（TD），该技术在几秒内仅可以提供视网膜的几个横截面图像。例如，Stratus OCT系统（Carl Zeiss Meditec，Inc.，Dublin，CA）在2秒内仅仅可以扫描黄斑的6个子午线横截面，这导致视网膜厚度图具有低分辨率。这样的低分辨率的图不足以用于精确的诊断。
目前，利用傅氏域光学相干断层照相（FD-OCT）技术的新一代视网膜OCT系统已经变得可用。FD-OCT比TD-OCT快得多。例如，RTVueFD-OCT系统（Optovue，Inc.，Fremont，CA）比Stratus TD-OCT快65倍。FD-OCT的快得多的扫描速度允许在较短时间段内在较大面积上进行更高密度的视网膜成图。较短的扫描时间减少了运动误差，并且较高的密度和扫描区域允许更详细的图案分析。
但是，FD-OCT的更快扫描速度不能直接转换为更精确的诊断。如果没有适当的方法来对扫描到的图像进行解译，该新技术的潜能就无法充分发挥。鉴于在现有技术中未能满足的需求以及尚未跨过的技术鸿沟，本发明的发明人已经设计了新OCT扫描图案和分析方法以实现FD-OCT在精确地测量视网膜厚度和诊断视神经病变方面的潜能。
通常，根据本发明实施方式的方法具有以下阶段：（1）图像获取；（2）图像处理、变换和分析；（3）诊断参数计算；以及（4）诊断确定。
由于本发明的方法基于非侵入成像技术，因此希望这些技术不限于人类，而且可以应用于除了人类之外的动物。因此，这里使用的术语“受验者”广泛地涉及具有类似于人眼的眼睛的任何有机体。本领域普通技术人员将会理解，具有解剖结构类似于人眼的眼睛的任何有机体都可以认为是本发明方法的可应用受验者。
1.图像获取
在第一阶段中，优选地使用非侵入成像技术来生成受验者眼睛的黄斑区域或乳头周围区域的初始图像。这里，黄斑区域以凹窝为中心，而乳头周围区域以视神经乳头为中心。适当的成像技术必须能够以足够高的分辨率和速度对区域进行成像，使得产生的图像具有足以用于诊断的质量。优选的是，该技术应当能够对不小于黄斑中的大约6mm×6mm或视神经盘的4mm×4mm的区域以不小于大约2秒每图像的速度进行成像。在某些优选的实施方式中，FD-OCT用于生成初始黄斑图。
本领域技术人员会理解，虽然这里描述的示例性实施方式基于FD-OCT，但是也可以使用其他成像技术（包括未来发展的成像技术），只要这些技术能够产生满足这里所述标准的图像即可。例如，有可能可以使用采用并行工作的检测器二维阵列以提高图像获取速度的新型TD-OCT，例如参考文献7所描述的。
但是，当将FD-OCT用作成像技术时，本发明还提供能够促进三维图像的获取和生成的新扫描图案。由于OCT图像是横截面扫描图像，因此三维模型的构造需要对个体横截面扫描进行对准并在扫描之间进行插值，以将最终的三维表示图拼合在一起。由于各扫描在时间上和空间上都是各自独立地进行的，因此对个体扫描进行对准是困难且费时的任务。横截面扫描的未对准将会导致黄斑区域的表示图不精确，这将限制该方法的诊断能力。
为了解决这个问题，发明人已经发明了能够促进快速且准确的对准的新扫描图案。参考图1，示出了如下的示例性扫描图案：其包括多条垂直扫描线，一条水平扫描线穿过所述多条垂直扫描线。所述多条垂直扫描线优选地按相等距离间隔开。水平扫描线优选地以垂直角度与所述多条垂直扫描线相交。该水平扫描线提供用于对准所有垂直扫描线的公共参照点，这大大地促进了对准和插值处理。
为了测量由青光眼造成的黄斑GCC丢失，更宽的扫描图案（尤其在垂直方向上）可以帮助找到在旁中心凹窝（parafovea）和周中心凹窝（perifovea）区域中的神经节细胞的丢失。由于视场功能测试覆盖了比当前OCT黄斑扫描图案宽得多的区域，因此更宽的OCT扫描图案也有助于找出功能丢失和结构损伤之间的关联。与径向扫描相比，网格扫描（gridscan）或光栅扫描（raster scan）在旁中心凹窝和周中心凹窝区域中可以创建具有更高横向分辨率的GCC厚度图。
在图1中所示的特定实施方式中，在这里将扫描图案称为黄斑图7毫米（MM7）扫描图案。利用RTVue FD-OCT系统获得FD-OCT图像。示例性扫描图案包含视网膜的16个垂直横截面以及1个水平横截面。但是，扫描线的数量并不具体受限于此。根据想要的扫描区域和分辨率，也可以适当地使用其他数量的扫描线。
2.图像处理、变换和分析
在第二阶段中，进一步处理在前一阶段获得的原始图像数据以生成黄斑区域的更精确的三维模型。根据用于获取原始数据的技术，可能需要不同的图像处理量来将原始图像数据变换为适合于进一步分析的格式。例如，如果成像技术需要模拟格式的原始数据，那么优选的是，将数据处理并变换为数字格式以有利于进一步分析。
参考图2，示出了视网膜的示例性横截面图像。这个处理后的图像使得能够测量视网膜和内视网膜层的厚度。本领域技术人员已经知道了测量视网膜和内视网膜层厚度的图像处理步骤（参见参考文献1-2，其内容通过引用而并入到本文中）。在与凹窝相交的垂直和水平横截面上识别凹窝的中心。在受验眼睛注视视觉目标的情况下获得扫描。但是，注视点可以稍稍偏离目标。凹窝的中心充当将厚度图的中心定位在其上的解剖学界标。
在其中原始数据由黄斑区域的OCT横截面扫描组成的优选实施方式中，通过在OCT横截面之间的插值，可以构造示例性GCC厚度图（图3）。由于在角落区域的测量是不可靠的，因此优选地剪切该图以仅仅保留中央的7mm圆形区域。优选的是，基于视网膜厚度图，该图以凹窝为中心。如从图3可以看到的，GCC厚度在凹窝中很薄或完全不存在。因此，从GCC的进一步分析中去除凹窝区域。如果在整个视网膜厚度上进行分析，那么可以使用凹窝区域。
一旦以适当格式处理并且设置原始图像数据，就可以设计各种图案分析算法并将其应用到数据以提取各种有用信息。从这些处理后的图像数据，可以派生出用于辅助诊断疾病状态的诊断参数。
本发明的一个示例性实施方式是测量视网膜的不正常地薄的区域内的视网膜组织丢失。这里，该参数称为焦点丢失量（FLV）。在下面的部分中，我们将首先描述基于黄斑GCC的分数丢失图和图案偏差图对FLV的一个示例性计算。接着，我们将进一步描述本发明的另选实施方式。
分数偏差图
在一个优选实施方式中，在这里称为“分数偏差图”的新派生图是根据GCC厚度图而计算出的。为了计算分数偏差图，需要标准参照物。
在示例性实施方式中，用于青光眼成像技术进展研究（AIGS）中的46个正常受验者用作标准参照物。本领域技术人员将会理解，这是一个统计学过程，也可以使用其他适当的数据集。
在AIG研究网站【3】上可以获得的AIGS程序手册（MOP）中限定了受验者招募和测试过程。根据没有患有青光眼、视神经疾病或视网膜疾病的一组正常眼睛来测量GCC厚度图。对正常参照群体的GCC厚度图进行平均以获得正常均值图（图4）。
接着，通过将进行测试的GCC图减去正常均值图，来计算偏差（D）图。
D_map=Map-Normal_mean_map
通过将D图除以正常均值图来计算分数偏差（FD）图。
FD_map=D_map/Normal_mean_map
3.诊断参数计算
在第三阶段中，根据在前两个阶段中获得的图像数据来定义并计算诊断参数。如上所述，可以定义并且计算各种诊断参数以辅助诊断疾病状态。在一个优选实施方式中，计算焦点丢失量（FLV）。
优选的是，利用统计学标准来检测不正常GCC变薄的区域。一个示例性标准是检测变薄到正常采样群体的5个百分点以下。在正常群体的FD图上进行逐点统计学计算以获得标准偏差（SD）图。正态分布的5个百分点以下的值（均值-1.64SD）被认为是很不正常的。根据GCC厚度、D图或FD图，可以计算不正常变薄的区域，并且结果是完全一样的。在一个示例性实施方式中，在FD图（图5）上识别焦点丢失的不正常区域，并且对不正常区域中的FD值进行合计（面积积分）以获得焦点丢失量（FLV）。当完全由FD来定义FLV时，FLV被称为FD-FLV。优选的是，通过划分图区域来对FD-FLV进行归一化，所以FD-FLV可以表示为百分数。因此，9%的FD_FLV（分数偏差焦点丢失量，fractional deviation focalloss volume）表示9%的神经节细胞丢失。
尽管根据GCC厚度FD图来计算上述FD-FLV示例，但是本领域技术人员将会理解，也可以根据D图来计算FD-FLV。诸如反射率图和乳头周围NFL图的其他类型的图也可以用作计算FD-FLV的基础。根据PD图计算FLV不太直接，但是提供了最好的诊断能力。
下面我们将简略描述这些另选实施方式。
图案偏差焦点丢失量（PD_FLV，pattern deviation focal loss volume）
在一个另选实施方式中，可以使用这里被称为图案偏差焦点丢失量（PD_FLV）的另一个示例性诊断参数。其计算描述如下。
1.导出图案偏差图
通过将GCC厚度图除以其平均值来从该图导出图案图（图6）。
Pattern_map=Map/Average
接着，将进行测试的眼睛的图案图减去正常参照群体的均值图案图，以获得图案偏差（PD）图（图6）。
PD_map=Pattern_map-Normal_mean_pattern_map
2.计算图案偏差量焦点丢失量（PD_FLV）
如在先前部分所述，利用GCC厚度、D或FD图来识别GCC变薄的不正常区域。在PD图（图7）上识别不正常区域，并且对不正常区域中的PD的负值进行合计（面积积分），以获得图案偏差焦点丢失量（PD_FLV）。这个过程与根据D或FD图计算FLV略有不同：在合计过程中PD的正值被设置为0（有效地忽略）。优选的是，通过对图区域进行划分来对PD_FLV进行归一化，所以其可以表示为百分数。
尽管根据GCC厚度图来计算上述PD-FLV示例，但是也可以根据诸如反射率图的其他类型图来计算PD-FLV示例。乳头周围NFL图也可以用作计算PD-FLV的基础。
焦点丢失量（FLV）和总体丢失量（GLV）
为了组合FD-FLV图和PD-FLV，将焦点丢失量（FLV）定义为由图案偏差识别出的不正常区域中的负分数偏差的总和。通常，通过对图区域进行划分来对FLV进行归一化，所以FLV可以表示为百分数。
FLV比FD_FLV和PD-FLV更具体，因为其仅仅合计其中GCC在绝对值和相对值方面都变薄的区域。
总体丢失量（GLV）被定义为整个区域中的负分数偏差的总和。通常，通过对图区域进行划分来对GLV进行归一化，所以GLV可以表示为百分数。
由于GLV在计算中需要较少的标准，因此其比FLV具有更好的可重复性。
尽管根据GCC厚度图来计算上述FLV和GLV，但是也可以根据诸如反射率图或神经纤维厚度图的其他类型图来计算FLV和GLV。乳头周围NFL图也可以用作计算FLV和GLV的基础。
用于检测异常的其他类型的图案分析
均值
可以对厚度图、偏差图或FD图计算总平均值和部分平均值。
总平均值是从图的所有有效区域平均的。
上方平均值是从图的上半球中的所有有效区域平均的。
下方平均值是从图的下半球中的所有有效区域平均的。
非对称
大多数情况下青光眼更严重地影响眼睛的下半部。但是，在很少的情况中，青光眼也可能更严重地影响眼睛的上半部分。因此，为了检测青光眼，优选地根据标准针对上方下方差（SID）计算绝对偏差。
SID=Superior_average-Inferior_average
上方下方差的绝对偏差（ADSID）是测试眼和平均正常眼的SID之间的差的绝对值。
AD_SID=Abs(SID-normal_mean_SID)
图案方差系数（均方根）
可以根据偏差图、FD图和PD图来计算图案方差系数（PCV）或均方根（RMS）。图案方差系数（PCV）或均方根（RMS）对于PD图最有用，并且提供了与正常图案的偏差的总和。例如，根据下述公式计算均方根图案偏差（RMS-PD）参数：
RMS_PD=(1/A)×(∫APD2dxdy)1/2，
其中，A是图面积，PD是图上的图案偏差值，x是图的水平尺寸，并且y是图的垂直尺寸。
其他类型的图
不同于上述图的视网膜图也可以用于计算如上所述的诊断参数。
其他厚度图
根据本发明，通过计算图案参数，视网膜总厚度图（图9）和NFL、GCL、IPL、INL的厚度图及其所有可能的组合都可以用于诊断青光眼。已经在青光眼诊断【1-2】中探索了这些组合的平均厚度。
反射率图
青光眼不仅造成内视网膜层变薄，而且减小了来自这些层的反射信号的幅度。因此，内视网膜反射率图也是有价值的。优选的是，通过取得内视网膜层中的平均信号率并将其除以来自参照层的平均信号，消除由于视网膜的外部因素（聚焦较差，介质混浊）而导致的反射率变化。参照层可以是亮度感光器内段、外段和它们的交叉处（IS-OS）、视网膜色素上皮细胞（PRE）、脉络膜血管层或它们的某一组合。IS-OS和RPE的组合称为外视网膜复合体（ORC）。可以对来自这些层的反射信号进行合计或平均以提供诊断信息。它们分别称为合计反射率图（图10）以及平均反射率图。合计反射率的分数偏差图是优选的。
乳头周围图
视神经乳头周围（乳头周围）的NFL或视网膜厚度图也提供用于青光眼诊断的信息。也可以使用乳头周围NFL反射率图。可以应用本发明中描述的所有图案参数。在RTVue OCT系统上，优选地利用视神经乳头4mm（ONH4）扫描图案来测量NFL厚度图（图11）。
4.诊断确定
在第四阶段中，基于计算出的诊断参数来确定诊断。在这个阶段中，可以有利地采用与其他参照数据的比较。在一个示例性实施方式中，进行关联分析以区分不同类型的视神经病变。
图案关联以区分各种类型的视神经病变
不同类型的视神经疾病造成GCC变薄的不同图案。青光眼相对地留出中心盲点区域并且具有可变的上域或下域（通常是下）。前部缺血性视神经病（AION）通常影响黄斑的上半部分（大多数情况）或下半部分。视神经炎经常与多发性硬化（MS）相关联。其产生弥漫性丢失。优选地使用利用关联的图案匹配来区分这些类型的视神经病变。
可以利用D图、FD图或PD图进行分析。优选地使用PD图。已经通过对各疾病类别内的眼睛图进行平均来搜集用于下述类型的视神经病变的特征图（图8）。
1.下方青光眼（IG）：视野青光眼的均值图，SID值>平均值+标准的1SD。
2.均匀青光眼（EG）：视野青光眼的均值图，SID值在标准的SD±1内。
3.上方青光眼（SG）：视野青光眼的均值图，SID值<平均值-标准的1SD。
4.视神经炎（ON）：视神经炎情况的均值图。
5.上方AION（SA）：AION眼睛的均值图，SID<正常值。
6.下方AION（IA）：AION眼睛的均值图，SID>正常值。
还可以类似地导出用于其他类型的视神经病变的特征图。例如，毒性、新陈代谢、和营养性视神经病变应当产生中央或中心盲点GCC丢失。
用所有特征性视神经病变的特征图，通过测试下的图的交叉关联来计算图案关联（PCC，pattern cross-correlation）值。接着，基于最高的PCC值选择诊断。可以根据诸如分数偏差图、图案偏差图和偏差图的多个图中的一个图来定义PCC。例如，PD-PCC可以被定义为：
PCC=∫A（PD×CPD）/[(∫A PD×PDdxdy)1/2(∫A CPD×CPDdxdy)1/2]其中，A是图的面积，PD是进行测试的眼睛的图案偏差，CPD是在考虑的疾病的特征图案偏差，x是图的水平尺寸，并且y是图的垂直尺寸。使用FD和D图，我们可以创建称为FD-PCC和D-PCC的类似参数。对于GCC图，我们优选使用分数偏差图来计算用于青光眼分析的PCC。对于乳头周围的NFL图，我们优选使用偏差图来计算用于青光眼分析的PCC。
诊断能力的比较
为了评估GCC导出的参数的诊断能力，我们使用受试者工作特性曲线下面积（AROC），该区域概括了在可用诊断阈值的整个范围上的诊断敏感性和特异性。数据来自AIGS，利用经受了用RTVue和Stratus测试二者的测试的正常和视野青光眼受验者的子集。表1示出了很多GCC导出的参数优于平均视网膜厚度。较复杂的图案分析优于简单的平均。用PD_FLV获得最好的性能。根据从偏差或分数偏差图计算的FLV图，也获得几乎相同的诊断能力。
提供下述特定示例以进一步说明本发明。
示例：青光眼中黄斑神经节细胞丢失的检测
方法
1.临床研究
包括在2003和2007的时期之间的预期的青光眼成像技术进展研究（AIGS）的参与者。这些参与者被分为四类：正常（N）、视野青光眼（PG）、疑似青光眼（GS）和预视野青光眼（PPG）。仅使用来自基线访问的数据。由于成员的青光眼状态是不确定的，因此该研究中不使用GS组。我们仅使用来自在研究期间采用FD-OCT的AIGS中心的数据。下面简单描述对于所分析的三个组的合格标准。
N组参与者的两只眼具有小于21mm汞柱的眼压（IOP）、正常的Humphrey SITA 24-2视场（VF）【平均偏差（MD）和图案标准偏差（PSD）在正态参照的95%界限之内，并且青光眼半视场测试（GHT）在97%界限之内】、大于等于500μm的中央角膜厚度、正常出现的视神经乳头、正常的神经纤维层、开放前房角、并且没有慢性眼睛或系统性皮质类固醇使用的历史。
PG组参与者的至少一只眼睛满足下述标准：青光眼（不正常）VF丢失【PSD（P<0.05）或者GHT（P<1%）在两个合格VF上的一致图案中的正常界限之外】、以及视神经乳头（ONH）变化，如弥漫性或局部性边缘变薄，盘（裂片）出血，垂直杯/盘比例比另一眼高出>0.2，边缘中的缺口，或者盘的步进性开洞的先前照相文档，视网膜神经沿的步进性变薄、或在裂隙灯生物显微镜上可见的NFL缺陷，或NFL的步进性丢失。
PPG组参与者具有与为PG组限定的ONH变化相同的标准。但是，PPG参与者眼睛的VF不满足用于PG组的合格标准。
用于AIGS中所有组的排除标准是：最佳校正视敏度低于20/40；年龄小于40或者大于79岁；等量球镜屈光不正大于+3.00D或者小于-7.00D；患有可能造成视野丢失或视神经盘畸形的糖尿病视网膜病变或其他疾病；或者先前进行过除了用后房IOL移植的不复杂白内障摘除术之外的眼内手术。
根据赫尔辛基宣言进行研究。在讨论了研究目的和参与结果之后，从所有参与者获得知情同意。该研究涉及的各机构的机构审查委员会批准了研究协议。在AIGS程序手册（该手册可从AIG研究网站下载，其内容以引用的方式结合在本文中）中可以找到AIG研究协议的进一步描述。
2.傅氏域光学相干断层照相
使用RTVue FD-OCT系统（Optovue，Inc.Fremont,CA）对患者进行扫描，该系统每秒钟获取26,000个轴位扫描（a扫描），并且具有5μm的深度分辨率（全幅宽半极大值）。相比之下，标准Stratus TD-OCT系统（Carl Zeiss Meditec，Dublin，CA）每秒钟获取400个a扫描并且具有10μm的分辨率。利用FD-OCT的更高速度，我们发明了黄斑区域的三维扫描，被称为黄斑图7毫米（MM7）扫描，其在7mm方形区域上均匀地采样黄斑（图1）。MM7协议的中心临时偏移0.75mm以改善临时外围的采样。MM7图案由从1条水平线和间隔0.5mm的15条垂直线的14928个a扫描构成。MM7图案的扫描时间是0.6秒。在各AIGS参与者的基线访问上获取三个MM7扫描。导出原始数据以用于进一步的图像处理。
3.图像处理
我们开发了自动的软件以绘制mGCCT图。首先，通过关联将15个垂直OCT横截面图像（参照用于示例性横截面图像的图2）与水平图像对准，以建立对齐的三维（3D）图像组。用中值滤波器和高斯滤波器的组合使图像平滑至较低分辨率，以抑制背景噪声和斑点噪声。接着，以较低的清晰度对它们进行重采样来加快计算速度。后续步骤使用所选择的各种分辨率和清晰度的图像，以最佳化处理的稳健性和速度。感光器色素上皮细胞复合体（PPC）带包括感光器内段-外段（IS/OS）交叉处的亮度带和视网膜色素上皮细胞，被检测为低分辨率图像中的第二（从内侧开始数）最大峰值。接着，IS/OS交叉处被检测为PPC内的第一最大强度峰值。由于来自叠加血管的阴影，PPC的小部分具有低信号；用相邻的a扫描替代这些有阴影的a扫描，以避免边界检测的干扰。在IS/OS交叉处对准图像以促进侧部平滑。内界膜（ILM）被识别为各a扫描的第一正梯度峰值。使用邻居约束和知识模型来将ILM峰值从假性噪声或分离的玻璃面区分出来。接着，识别内丛状层（IPL）的外边界。为了提高边界检测的稳健性，应用步进性的改进过程。该过程开始于对低分辨率（高度低通滤波）3D数据集的边界检测，接着对分辨率逐渐提高的数据逐渐改善边界。从ILM至外IPL边界测量GCC厚度。从ILM至IS/OS交叉处测量视网膜厚度。通过对来自MM7 3D数据集中的16个OCT横截面图像的厚度轮廓进行插值来计算mGCCT和mRT图。在mRT图上识别凹窝的位置，并且将该位置用于重定位图的垂直位置。对图进行剪切以移除分段不太可靠的外围区域。剩余区域是在7mm直径圆内部并且距中央水平线在3mm内的区域。对于mGCCT图，由于GCC过于薄而无法可靠测量，因此也不包括凹窝中心的0.5mm内的区域（1mm直径圆）。
4.诊断参数的导出
图12示出了视野青光眼（PG）情况中的图集。我们基于mGCCT图（图12A)计算多个青光眼诊断参数。最简单的是总平均厚度（mGCCT-AVG）。由于青光眼往往在下方产生的损伤较多，因此我们还计算上半球平均值和下半球平均值之间的差（mGCCT-SID）。
为了从mGCCT图提取甚至更多的诊断信息，我们开发了分析mGCCT丢失图案的方法。为此，我们计算mGCCT丢失图：分数偏差（FD）图和图案偏差（PD）图。首先，对所有正常眼睛的GCC图逐点进行平均，以创建正常参照图。FD图（图12B）是考虑的mGCCT图减去正常参照图，除以正常参照图。图案图是按照GCC厚度图自身的总平均值进行了归一化（除以GCC厚度图自身的总平均值）的GCC厚度图。图案偏差（PD）图（图12C）是考虑的图案图减去正常参照图案。FD图示出了GCC丢失的百分比。PD图示出了mGCCT图案与正常图案的不同之处。
接着，根据这两个派生图计算三个基于图案的诊断参数。焦点丢失量（FLV）是有显著焦点丢失的区域中的总和FD。显著焦点丢失被定义为比正常平均值低1.65标准偏差（SD）以上（在正态分布的第五个百分点之下）的FD。总体丢失量（GLV）是FD为负的区域中的FD的总和。图案方差系数（PCV）是PD图的均方根。
在MATLAB 7.0中对图像处理和诊断参数计算进行编程。
5.时域光学相干断层照相
也通过Stratus OCT（Carl Zeiss Meditec,Inc,Dublin,CA），利用标准的快速视网膜神经纤维层（RNFL）扫描和更快的黄斑厚度图扫描，来扫描所有的参与者。使用标准Stratus 4.0软件来计算cpNFLT和mRT的总平均值。
6.统计学分析
分析各参与者的两只眼。通过使用广义估计方程（GEE）方法或线性混合模型，在统计学测试中给出眼间关联。
使用类内关联、共同SD（pooled SD）和方差系数（CV）来评估诊断参数的可再现性。根据线性混合模型来计算这些指标，在线性混合模型中，分别使用对于受验者和眼睛的方差成分来说明重复的测量和眼间关联。
使用受验者工作特征（AROC）曲线下面积来比较诊断能力。为了给出眼间关联，基于Obuchowski的公式来计算AROC，该公式扩展了应用于聚类数据的Delong等人的非参数方法。在眼科学中的先前研究中已经使用了该方法以处理眼间关联。
为了调整N、PG和PPG组之间的年龄不平衡，使用以年龄和诊断参数为协变量的GEE逻辑回归模型来生成AROC。弥补年龄不平衡的这样方法已经用在眼科学的先前研究中。
为了比较均值，我们对满足正态分布的参数使用t-检验。发现了满足Gamma概率分布（非正态分布）的多个诊断参数。使用对于适当Gamma分布具有广义线性模型的华尔德测试，来比较这些参数的均值。对检验使用眼间关联的GEE调整。由于我们假定了患有疾病的组中的均值比正常组中的低，因此，以单尾方式进行这些检验。
AROC计算写入在MATLAB 7.0软件中，并且用SAS 9.1软件进行其他统计计算。统计学意义的临界阿尔法水平被设定在0.05。
结果
从AIG中心数据库中识别出具有可用RTVue FD-OCT MM7扫描和有效的Stratus TD-OCT扫描的总共180个参与者（328只眼）。由于对于所有三个重复MM7扫描的可见地不精确的分段，因此剔除了14个参与者的15只眼。分析来自179个参与者的剩余313只合格眼。在表1中概括了各组的人口信息和临床信息。预视野青光眼和PG参与者比N参与者（P<0.0001）的年龄大。在后续分析中适当地处理年龄不平衡，如在方法部分中所述的。与PG组相比，在N组中有更多的白种人。但是，从N组中诊断参数的均值的角度来说，在种族组之间没有显著差别。如所期待的眼睛状态和对眼睛状态的分类，从MD和PSD测量的角度来说，在VF测试中N眼表现得比PPG和PG眼要好。N眼具有比PPG和PG眼更低的IOP和更厚的中央角膜厚度（CCT）。在PPG眼的IOP和PG眼的CCT中差别很显著。
为了基于MD将PG眼分类在青光眼的不同阶段，79只眼（70.5%）具有早期青光眼（MD≥-6.0dB），25只眼（22.3%）具有中度青光眼（MD在-6.01至-12.0dB之间），并且8只眼（7.1%）具有高度青光眼（MD<-12dB）。
表2按组概括了各诊断参数的分布统计信息。与N组（P<0.001）相比，所有参数在PPG和PG组中明显更差。因为SID、PCV、FLV和GLV具有非正态分布（Gamma分布），因此使用华尔德测试来比较这些参数，如在方法部分中所述。
青光眼中的mGCCT丢失的特征图案是双拱形（bi-arcuate）的，具有较大的下方缺陷（图13）。
通过三个测度：在相同时期取得的重复测量的IC、共同SD和CV（表3），来评估可再现性。PPG和PG组中的可再现性很重要，因为其提供了参数可以以何种程度在疾病的阶段中跟踪疾病发展的指示。总体来说，FD-OCT mGCCT和mRT均值以及GLV具有优秀的可再现性（在PG和PPG组中ICC=0.99并且CV<1.3%）。尽管TD-OCT mRT和cpNFLT均值也具有很好的再现性，但是它们不如类似的FD-OCT参数好。
AROC提供了对诊断能力的测量（表3）。通过FD-OCT和TD-OCT测量的mRT均值具有相同的AROC。通过分离内视网膜，mGCCT-AVG改善了PG的诊断（AROC=0.90）。这明显好于mRT的情况（P=0.01）。基于FLV和GLV图案的参数在诊断PG时表现得甚至更好。对于GLV，增加是显著的（P=0.01）。黄斑参数mGCCT-AVG、mGCCT-FLV和mGCCT-GLV具有与cpNFLT-AVG类似的诊断能力。对于PPG（与N相对）的诊断，我们发现与mGCCT参数相比mRT没有优势。
对于每10μm的组织丢失，患有青光眼的几率（95%置信区间）对于mGCCT-AVG是7.43（4.13，13.36），对于cpNFL是4.88（2.64，9.03），对于FD-OCT是2.68（1.96，3.65），并且对于TD-OCT mRT是2.48（1.83，3.35）。我们注意到对于每10μm的组织丢失，丢失GCCT-AVG而患有青光眼的几率大约是丢失cpNFL-AVG而患有青光眼的几率的1.5倍。
在PG和PPG的示例情况中示出了用神经盘照相术和VF发现的mGCCT的关联。在这两种情况下，mGCCT FD图示出了丢失的典型双拱形图案。在PG的情况下，显著的下方GCC丢失与下方神经盘边缘丢失和上方VF缺陷很好地关联。在PPG的情况下，GCC丢失是焦点，并且通过基于图案的参数而不是均值来找出不正常。
结论
在这个示例中，我们示出了根据本发明的新诊断参数的应用以寻找黄斑中的青光眼。FD-OCT（65x Stratus TD-OCT）的较快速度使得能够在黄斑的较大区域上进行高密度扫描，而运动伪影较少。RTVue FD-OCT装置的分辨率也比传统时域OCT（例如，Stratus TD-OCT）好两倍。更高的清晰度（更密的采样）和更高分辨率的组合提高了mGCCT测量的精度和稳健性。
我们还发现，在这个示例中，与由FD-OCT或TD-OCT测量的mRT均值相比，由RTVue FD-OCT所测量的mGCCT均值在诊断PG组中的青光眼时明显要好。因此，从外视网膜分离出GCC提高了诊断能力。虽然不旨在限于任何特定理论或说明，但是我们相信这可以由下述事实来解释：受青光眼影响不大的外视网膜占据总视网膜厚度的65%至70%，因此提供了降低判别能力的厚度方差。在PPG和N眼之间进行区分时mGCCT的诊断能力也比mRT的诊断能力高，但是优势在统计上不显著。这可以由下述事实来解释：PPG组的大小较小，并且PPG组中的某些眼有可能实际上并没有患青光眼（PPG眼具有正常的或临界VF）。
而且，与Stratus cpNFL厚度相比，通过FD-OCT或Stratus TD-OCT的mRT对于青光眼检测是较不敏感（具有较低AROC）的参数。与用于青光眼检测的Stratus mRT相比，包括Wollestein等人【4】和Guedes等人【5】的其他研究者已经报告了对于Stratus cpNFL厚度的更高AROC。在当前的研究中，对于对总黄斑厚度的测量，FD-OCT不提供优于TD-OCT的任何显著优势。不同于（由于必须对视网膜层分段而）需要由FD-OCT提供的更高分辨率以及细节的GCC的边界检测，由于ILM和IS/OS交叉处的边界是清楚定义的，因此可以通过较低分辨率的StratusTD-OCT很好地进行对mRT的边界检测。
通过FD-OCT的更高速度，使得对黄斑区域的更宽更清楚的采样是可能的。这有助于对GCC丢失的图案分析。我们设计了着眼于GCC丢失图案的不同方面的几个基于图案的参数，这些参数可以以补充方式来使用。基于青光眼经常具有下域不对称性的观察结果，设计SID参数以检测GCC丢失不对称的情况。GLV和FLV参数以不同的特异性水平概括了黄斑中的GCC丢失量。FLV参数是更具体的，因为其仅合计了其中GCC在绝对（GCCT<正常值）和相对（PD<5个百分点）角度都很薄的区域中的丢失。PCV参数完全基于PD图，并且检测GCCT图案中的任何变化。我们发现，对于PG的诊断，FLV和GLV的诊断能力比简单平值的高。这可以由包括某些情况来解释：其中眼睛可以开始具有上述总体平均GCCT；因此，着眼于总体平均值只能在后期检测到青光眼。在这些情况中，在GCCT图案中寻找不正常可以更早地检测到青光眼。在图14中存在这样的情况，其中mGCCT平均值是正常的，但是基于图案的参数是不正常的。由于沿黄斑丘疹束（maculapapular bundle）的GCCT仍然在平均值之上（图14），所以在发展为青光眼之前，这只眼可能比平均GCC更厚。这只眼具有焦点区域，其中mGCCT比正常的薄30%，这足以被我们的软件识别为不正常的焦点丢失。这个区域对应于PD在-3和-4dB之间的VF上的区域，仍然在正常方差范围之内。这种情况示出了利用mGCCT图案分析作为在患者出现确切的VF缺陷之前诊断青光眼的早期检测方法。其还示出了在VF缺陷是临界值的情况下mGCCT图案可以是有用的相关参数。读者应当意识到，在视网膜上的每一毫米对应于VF上的大约3.5°。因此，MM7mGCCT图（7mm×6mm）向上和向下对着大约11°，向鼻对着10°，在时间上对着15°。它覆盖了标准Humphrey 24-2VF的大约一半区域（图12E和14C），当然，它由于眼睛的光学投影而相对于VF是上/下颠倒的。
综上所述，我们已经示出了根据本发明实施方式的利用FD-OCT的更高速度和更高分辨率的宽黄斑扫描图案。我们还开发了GCC成图软件、以及基于图案的新诊断参数。与从总视网膜厚度导出的参数相比，新mGCCT参数能够按更高灵敏度和特异性将青光眼与非青光眼区分开来。就自身来说，mGCCT参数的诊断能力类似于cpNFL参数，并且可以以补充方式来使用。mGCCT图可以直接与VF图的中央部分关联。某些mGCCT参数是高度可再现的，从而可用在跟踪青光眼发展中。
尽管已经从具体示例性实施方式和示例的角度描述了本发明，但是将会理解，这里公开的实施方式仅用于说明目的，并且本领域技术人员可以做出各种修改和变化，而不脱离所附权利要求所述的本发明的精神和范围。
参考文献
通过引用并入以下各篇参考文献
1.Tan O,Li G，Lu A.等人,Mapping of Macular Substructures withOptical Coherence Tomography for Glaucoma Diagnosis,Ophthalmology2008；l 15(6)：949-56
2.Ishikawa H,Stein D,Wollstein G等人,Macular segmentation withoptical coherence tomography.Invest Ophthalmol Vis Sci 2005；46：2012-7.
3.URL=http://www.aigstudy.net/
4.Wollstein G，Schuman JS,Price LL等人,Optical coherencetomography(OCT)macular and peripapillary retinal nerve fiber layermeasurements and automated visual fields.Am J Ophthalmol 2004；138(2)：218-25.
5.Guedes V,Schuman JS,Hertzmark E等人,Optical coherencetomography measurement of macular and nerve fiber layer thickness innormal and glaucomatous human eyes.Ophthalmology 2003；110(l)：177-89.
6.Zeimer R,Asrani S,Zou S等人，Quantitative detection ofglaucomatous damage at the posterior pole by retinal thickness mapping.Apilot study.Ophthalmology 1998；105(2)：224-31.
7.Watanabe Y,Sato M.Quasi-single shot axial-lateral parallel timedomain optical coherence tomography with Hubert transformation.OptExpress 2008；16(2)：524-34.
表1 研究群体的特征
Figure GDA00001875003200241
缩写：MD=视场平均偏差；PSD=视场图案标准偏差；IOP=眼压；CCT=中央角膜厚度
P*：用于比较N和PPG组的P值
P**：用于比较N和PG组的P值
***某些参与者的一只眼睛被诊断为PPG并且另一只眼睛被诊断为PG
****均值±标准偏差
表2 诊断参数的按组分布
Figure GDA00001875003200251
缩写：SD=标准偏差；mGCCT=黄斑神经节细胞复合体厚度；AVG=均值；SID=上方-下方差；PCV=图案方差系数；FLV=焦点丢失量；GLV=总体丢失量；mRT=黄斑视网膜厚度；cpNFLT=乳头周围神经纤维层厚度
表3 诊断参数的可再现性
Figure GDA00001875003200261
缩写：ICC=类内关联，CV=方差系数；诊断参数的缩写与表2中相同。
表4 参数的诊断能力
Figure GDA00001875003200271
缩写：AROC=受验者工作曲线下面积；SE=标准误差；诊断参数的缩写与表2中相同。
*由于mGCCT-SID在广义估计方程逻辑回归模型（P值=0.11）中不重要，因此不产生AROC。

Claims (11)

1.一种用于分析视网膜组织结构的光学相干断层照相OCT图像以获得特性参数的值的方法，所述特性参数是从焦点丢失量FLV、总体丢失量GLV或者焦点丢失量FLV和总体丢失量GLV二者中选择的，所述方法包括以下步骤：
（1）使用OCT来获得从包括以凹窝为中心的黄斑图和以视神经乳头为中心的乳头周围图的组中选择的初始图，其中，所述初始图由光栅或网格扫描图案构成，并且其中，所述黄斑图覆盖至少6mm×6mm的区域，所述乳头周围图覆盖至少4mm×4mm的区域；
（2）基于所述初始图来构造神经节细胞复合体厚度的图；
（3）计算所述神经节细胞复合体厚度的偏差图；
（4）对所述神经节细胞复合体厚度的图或所述偏差图应用统计学方法；以及
（5）基于对所述偏差图的分析，计算FLV、GLV或者FLV和GLV二者的值。
2.根据权利要求1所述的方法，其中，所述OCT是从傅氏域OCT或时域OCT中选择的一个。
3.根据权利要求1所述的方法，其中，所述扫描图案是黄斑图7毫米扫描图案。
4.根据权利要求1所述的方法，其中，构造神经节细胞复合体厚度的图的步骤包括：在所述初始图的扫描线之间进行插值，以形成所扫描区域的三维模型。
5.根据权利要求1所述的方法，其中，所述统计学方法是从标准偏差比较、总平均、上方平均、下方平均、上方-下方差、均方根或这些统计学方法的组合中选择的一个。
6.根据权利要求5所述的方法，其中，所述标准偏差比较是逐数据点比较，其中识别值高于或低于标准值的预定百分比的数据点。
7.根据权利要求6所述的方法，其中，所述预定百分比是从0.5%至5%。
8.根据权利要求6所述的方法，其中，所述预定百分比是5%。
9.一种用于分析视网膜组织结构的光学相干断层照相OCT图像以获得特性参数的值的方法，所述特性参数是从焦点丢失量FLV、总体丢失量GLV或者焦点丢失量FLV和总体丢失量GLV二者中选择的，该方法包括以下步骤：
获取步骤，获取视网膜组织结构的OCT黄斑图，其中，所述黄斑图的尺寸至少是6mm×6mm；
处理步骤，对所述黄斑图进行处理，以产生黄斑区域的三维模型以及视网膜厚度图；
产生步骤，根据所述视网膜厚度图产生一个或更多个派生图，其中，所述一个或更多个派生图是从包含偏差值的偏差图、包含分数偏差值的分数偏差图、或者包含图案偏差值的图案偏差图中选择的；
识别步骤，分析所述视网膜厚度图和/或所述一个或更多个派生图；以及
计算FLV或GLV的值，其中，通过对所述偏差图中的负分数偏差进行求和来计算FLV，通过对整个所述偏差图中的负分数偏差进行求和来计算GLV。
10.根据权利要求9所述的方法，其中，通过傅氏域OCT来产生所述OCT黄斑图。
11.根据权利要求10所述的方法，其中，所述获取步骤包括使用网格扫描图案或光栅扫描图案来进行傅氏域OCT扫描。
CN2008801033823A 2007-06-15 2008-06-13 用于分析光学相干断层照相图像的方法 Active CN101778593B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US94444907P 2007-06-15 2007-06-15
US60/944,449 2007-06-15
PCT/US2008/066987 WO2008157406A1 (en) 2007-06-15 2008-06-13 Pattern analysis of retinal maps for diagnosis of optic nerve diseases by optical coherence tomography

Publications (2)

Publication Number Publication Date
CN101778593A CN101778593A (zh) 2010-07-14
CN101778593B true CN101778593B (zh) 2013-03-20

Family

ID=40131971

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008801033823A Active CN101778593B (zh) 2007-06-15 2008-06-13 用于分析光学相干断层照相图像的方法

Country Status (7)

Country Link
US (1) US8474978B2 (zh)
EP (1) EP2162054B1 (zh)
JP (1) JP5832744B2 (zh)
CN (1) CN101778593B (zh)
CA (1) CA2690834C (zh)
ES (1) ES2544585T3 (zh)
WO (1) WO2008157406A1 (zh)

Families Citing this family (86)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080312552A1 (en) * 2007-06-18 2008-12-18 Qienyuan Zhou Method to detect change in tissue measurements
WO2009126112A1 (en) * 2008-04-08 2009-10-15 National University Of Singapore Retinal image analysis systems and methods
US8718743B2 (en) * 2008-04-24 2014-05-06 Duke University Methods for single-pass volumetric bidirectional blood flow imaging spectral domain optical coherence tomography using a modified hilbert transform
WO2009139722A1 (en) * 2008-05-14 2009-11-19 Agency For Science, Technology And Research Automatic cup-to-disc ratio measurement system
US8896682B2 (en) 2008-12-19 2014-11-25 The Johns Hopkins University System and method for automated detection of age related macular degeneration and other retinal abnormalities
JP5543126B2 (ja) * 2009-04-16 2014-07-09 キヤノン株式会社 医用画像処理装置及びその制御方法
WO2010129544A1 (en) 2009-05-04 2010-11-11 Duke University Methods and computer program products for quantitative three-dimensional image correction and clinical parameter computation in optical coherence tomography
JP5436076B2 (ja) 2009-07-14 2014-03-05 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム
US20110190657A1 (en) * 2009-08-10 2011-08-04 Carl Zeiss Meditec, Inc. Glaucoma combinatorial analysis
US8419186B2 (en) * 2009-09-30 2013-04-16 Nidek Co., Ltd. Fundus observation apparatus
US8322853B2 (en) * 2009-10-02 2012-12-04 Optos Plc Diagnostic method and apparatus for predicting potential preserved visual acuity
JP5943338B2 (ja) * 2010-01-20 2016-07-05 学校法人慶應義塾 糖尿病網膜症未発症の時期の糖尿病患者の網膜の外顆粒層の厚みに対する内網状層の相対的な厚みの変化、前記外顆粒層の厚みに対する内顆粒層の相対的な厚みの変化、及び神経節細胞層の厚みの変化の測定方法
KR101095302B1 (ko) * 2010-05-14 2011-12-16 가톨릭대학교 산학협력단 부위별 시신경 두께 패턴을 이용한 시신경질환 자동진단 프로그램장치 및 그 프로그램방법
EP2609851B1 (en) * 2010-08-24 2020-02-26 Kowa Company, Ltd. Visual field examination system
JP6180073B2 (ja) * 2010-08-31 2017-08-16 キヤノン株式会社 画像処理装置及びその制御方法、並びに、プログラム
US20120063656A1 (en) * 2010-09-13 2012-03-15 University Of Southern California Efficient mapping of tissue properties from unregistered data with low signal-to-noise ratio
JP5588291B2 (ja) * 2010-09-29 2014-09-10 キヤノン株式会社 情報処理装置、情報処理方法、情報処理システム、及びプログラム
JP2011200632A (ja) * 2010-11-04 2011-10-13 Tomohiro Tsuta 適応光学眼底観察装置
DE102010051281A1 (de) * 2010-11-12 2012-05-16 Carl Zeiss Meditec Ag Verfahren zur modellbasierten Bestimmung der Biometrie von Augen
JP5702991B2 (ja) 2010-11-19 2015-04-15 キヤノン株式会社 画像処理装置及び画像処理方法
EP2646768B1 (en) 2010-12-03 2020-09-09 Optovue, Inc. Method and imaging system of generating a total corneal power map
WO2012078636A1 (en) 2010-12-07 2012-06-14 University Of Iowa Research Foundation Optimal, user-friendly, object background separation
US8801187B1 (en) 2010-12-23 2014-08-12 Carl Zeiss Meditec, Inc. Methods to reduce variance in OCT analysis of the macula
AU2012207076A1 (en) 2011-01-20 2013-08-15 University Of Iowa Research Foundation Automated determination of arteriovenous ratio in images of blood vessels
WO2012100213A2 (en) 2011-01-21 2012-07-26 Duke University Systems and methods for complex conjugate artifact resolved optical coherence tomography
JP5904711B2 (ja) 2011-02-01 2016-04-20 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム
EP2691007A4 (en) * 2011-03-29 2014-09-10 Steven Verdooner DEVICE AND METHOD FOR IDENTIFYING ONE OR MORE AMYLOID BETA PLAQUES IN SEVERAL SEPARATE OCT NETWORK LAYERS
FR2973997B1 (fr) * 2011-04-15 2014-08-22 Univ Paris Curie Dispositif d'aide a la detection de caracteristiques anatomiques d'au moins une portion d'un tissu. procede d'aide a la detection de caracteristiques anatomiques d'au moins une portion d'un tissu
US9357911B2 (en) 2011-05-09 2016-06-07 Carl Zeiss Meditec, Inc. Integration and fusion of data from diagnostic measurements for glaucoma detection and progression analysis
JP6025311B2 (ja) * 2011-08-01 2016-11-16 キヤノン株式会社 眼科診断支援装置および方法
JP2014527434A (ja) * 2011-08-09 2014-10-16 オプトビュー，インコーポレーテッド 光干渉断層法におけるフィーチャの動き補正及び正規化
JP5017491B2 (ja) * 2011-11-22 2012-09-05 キヤノン株式会社 眼科装置
JP6146952B2 (ja) 2012-01-27 2017-06-14 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム。
JP5936254B2 (ja) * 2012-02-29 2016-06-22 国立大学法人京都大学 眼底観察装置及び眼底画像解析装置
JP6025349B2 (ja) * 2012-03-08 2016-11-16 キヤノン株式会社 画像処理装置、光干渉断層撮像装置、画像処理方法および光干渉断層撮像方法
US9545196B2 (en) 2012-05-04 2017-01-17 University Of Iowa Research Foundation Automated assessment of glaucoma loss from optical coherence tomography
US9192294B2 (en) * 2012-05-10 2015-11-24 Carl Zeiss Meditec, Inc. Systems and methods for faster optical coherence tomography acquisition and processing
US9117121B2 (en) * 2012-05-21 2015-08-25 The Chinese University Of Hong Kong Detection of disease-related retinal nerve fiber layer thinning
US9579016B2 (en) 2012-06-15 2017-02-28 Oregon Health & Science University Non-invasive 3D imaging and measuring of anterior chamber angle of the eye
US20140160263A1 (en) * 2012-11-30 2014-06-12 Kabushiki Kaisha Topcon Data visualization method and apparatus utilizing receiver operating characteristic analysis
US9677869B2 (en) 2012-12-05 2017-06-13 Perimeter Medical Imaging, Inc. System and method for generating a wide-field OCT image of a portion of a sample
JP6147001B2 (ja) * 2012-12-28 2017-06-14 キヤノン株式会社 画像処理装置および画像処理方法
WO2014160454A1 (en) * 2013-03-13 2014-10-02 Riverside Research Institute Methods for diagnosing vitreo-retinal disease
US10497124B2 (en) * 2013-03-15 2019-12-03 Kabushiki Kaisha Topcon Optic disc image segmentation method and apparatus
WO2014194317A1 (en) 2013-05-31 2014-12-04 Covidien Lp Surgical device with an end-effector assembly and system for monitoring of tissue during a surgical procedure
WO2015017536A1 (en) * 2013-07-31 2015-02-05 The Board Of Trustees Of The Leland Stanford Junior University Method and system for evaluating progression of age-related macular degeneration
CN103405215B (zh) * 2013-08-16 2015-09-30 江苏省人民医院 一种活体小鼠视网膜多参数测量的方法
CN103699901A (zh) * 2013-12-17 2014-04-02 苏州大学 一种基于支持向量机的3d oct视网膜图像中is/os缺失区域自动检测方法
WO2015171566A1 (en) * 2014-05-06 2015-11-12 Oregon Health & Science University Aqueous cell differentiation in anterior uveitis using optical coherence tomography
WO2016109016A1 (en) * 2014-11-04 2016-07-07 Oregon Health & Science University Spectral fractionation detection of gold nanorod contrast agents using optical coherence tomography
CN104434026B (zh) * 2014-12-17 2016-08-17 深圳市斯尔顿科技有限公司 视网膜注视点偏离黄斑中心凹的检测方法
WO2016109750A1 (en) 2014-12-30 2016-07-07 Optovue, Inc. Methods and apparatus for retina bloodvessel assesment with oct angiography
US9943223B2 (en) 2015-02-13 2018-04-17 University Of Miami Retinal nerve fiber layer volume analysis for detection and progression analysis of glaucoma
KR101667279B1 (ko) 2015-02-25 2016-10-19 연세대학교 산학협력단 빛 간섭 단층촬영을 이용한 망막 내 병변 검출 장치 및 방법
JP6701587B2 (ja) * 2015-03-31 2020-05-27 国立大学法人東北大学 画像表示装置、処理方法、及び、プログラム
KR101683573B1 (ko) 2015-03-31 2016-12-07 연세대학교 산학협력단 안구 단층 촬영 영상에서의 시세포 층 정량화 검출 장치 및 그 방법
US10115194B2 (en) 2015-04-06 2018-10-30 IDx, LLC Systems and methods for feature detection in retinal images
WO2017040705A1 (en) 2015-09-01 2017-03-09 Oregon Health & Science University Systems and methods of gluacoma diagnosis based on frequency analysis of inner retinal surface profile measured by optical coherence tomography
JP6580448B2 (ja) * 2015-10-16 2019-09-25 株式会社トプコン 眼科撮影装置及び眼科情報処理装置
WO2018094381A1 (en) * 2016-11-21 2018-05-24 Tecumseh Vision, Llc System and method for automatic assessment of disease condition using oct scan data
EP3255605B1 (en) * 2015-11-16 2021-06-02 Huawei Technologies Co., Ltd. Method and apparatus for smoothing video
US9706917B1 (en) 2016-03-24 2017-07-18 Don C. Bagwell, Sr. Eye examination template
CN105919550B (zh) * 2016-06-16 2017-09-29 温州医科大学 一种基于光学相干层析成像技术的黄斑中心凹定位测量偏心量的方法
US20180000342A1 (en) * 2016-06-30 2018-01-04 Oregon Health And Science University Diagnostic classification of corneal shape abnormalities
US20180014724A1 (en) * 2016-07-18 2018-01-18 Dariusz Wroblewski Method and System for Analysis of Diagnostic Parameters and Disease Progression
US9968251B2 (en) 2016-09-30 2018-05-15 Carl Zeiss Meditec, Inc. Combined structure-function guided progression analysis
CN109843146B (zh) * 2016-10-14 2022-06-14 爱尔康公司 光学相干断层成像术交叉视图成像
KR20230010843A (ko) * 2016-12-05 2023-01-19 더 잭슨 래보라토리 인간에서의 신경퇴화 및 녹내장에 대한 진단적 마커로서의 망막 및 시신경 내 지방 알갱이
EP3558091A4 (en) 2016-12-21 2020-12-02 Acucela, Inc. MINIATURIZED AFFORDABLE OPTICAL COHERENCE TOMOGRAPHY SYSTEM FOR OPHTHALMIC APPLICATIONS IN THE HOME
WO2019014767A1 (en) 2017-07-18 2019-01-24 Perimeter Medical Imaging, Inc. SAMPLE CONTAINER FOR STABILIZING AND ALIGNING EXCISED ORGANIC TISSUE SAMPLES FOR EX VIVO ANALYSIS
US10694995B2 (en) 2017-12-05 2020-06-30 Renegade Optophysics, Llc Diagnostic eye goggle system
US11080850B2 (en) 2018-01-16 2021-08-03 Electronics And Telecommunications Research Institute Glaucoma diagnosis method using fundus image and apparatus for the same
WO2019178185A1 (en) * 2018-03-13 2019-09-19 The Uab Research Foundation Colocalized detection of retinal perfusion and optic nerve head deformations
EP3809948A4 (en) 2018-06-20 2022-03-16 Acucela Inc. MINIATURIZED MOBILE, LOW COST OPTICAL COHERENCE TOMOGRAPHY SYSTEM FOR HOME OPHTHALMIC APPLICATIONS
US11432719B2 (en) * 2019-03-14 2022-09-06 Oregon Health & Science University Visual field simulation using optical coherence tomography and optical coherence tomographic angiography
RU2713832C1 (ru) * 2019-04-18 2020-02-07 федеральное государственное бюджетное образовательное учреждение высшего образования "Московский государственный медико-стоматологический университет имени А.И. Евдокимова" Министерства здравоохранения Российской Федерации Способ диагностики компрессионной оптической нейропатии при патологии хиазмально-селлярной области
JP6736734B2 (ja) * 2019-07-23 2020-08-05 株式会社トプコン 眼科撮影装置及び眼科情報処理装置
CN111179226B (zh) * 2019-12-14 2024-04-19 中国科学院深圳先进技术研究院 视野图的识别方法、装置及计算机存储介质
EP4081096A4 (en) 2019-12-26 2024-01-10 Acucela Inc OPTICAL COHERENCE TOMOGRAPHY PATIENT ALIGNMENT SYSTEM FOR HOME-BASED OPHTHALMIC APPLICATIONS
JP2021104282A (ja) * 2019-12-27 2021-07-26 株式会社ニコン 眼科装置
CN113781380B (zh) * 2020-06-05 2023-09-29 中山大学中山眼科中心 一种鉴别视神经脊髓炎和原发性开角型青光眼的系统
US10959613B1 (en) 2020-08-04 2021-03-30 Acucela Inc. Scan pattern and signal processing for optical coherence tomography
AU2021324968A1 (en) 2020-08-14 2023-03-02 Acucela Inc. System and method for optical coherence tomography a-scan decurving
US11393094B2 (en) 2020-09-11 2022-07-19 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
CN116322471A (zh) 2020-09-30 2023-06-23 奥克塞拉有限公司 近视预测、诊断、计划和监测设备
WO2022204622A1 (en) 2021-03-24 2022-09-29 Acucela Inc. Axial length measurement monitor

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6104828A (en) * 1994-03-24 2000-08-15 Kabushiki Kaisha Topcon Ophthalmologic image processor
US6055451A (en) * 1997-12-12 2000-04-25 Spectrx, Inc. Apparatus and method for determining tissue characteristics
US20030135122A1 (en) * 1997-12-12 2003-07-17 Spectrx, Inc. Multi-modal optical tissue diagnostic system
CA2343401C (en) * 1998-09-11 2009-01-27 Spectrx, Inc. Multi-modal optical tissue diagnostic system
WO2003105678A2 (en) * 2002-06-12 2003-12-24 Advanced Research And Technology Institute, Inc. Method and apparatus for improving both lateral and axial resolution in ophthalmoscopy
US7347548B2 (en) * 2003-05-01 2008-03-25 The Cleveland Clinic Foundation Method and apparatus for measuring a retinal sublayer characteristic
WO2006022045A1 (ja) * 2004-08-26 2006-03-02 National University Corporation Nagoya University 光干渉断層装置
US7301644B2 (en) * 2004-12-02 2007-11-27 University Of Miami Enhanced optical coherence tomography for anatomical mapping
JP4827745B2 (ja) * 2005-01-19 2011-11-30 株式会社日立メディコ 画像診断支援装置及び画像診断支援プログラム
EP2417903A1 (en) * 2005-01-21 2012-02-15 Massachusetts Institute of Technology Methods and apparatus for optical coherence tomography scanning
JP4850495B2 (ja) * 2005-10-12 2012-01-11 株式会社トプコン 眼底観察装置及び眼底観察プログラム
US7768652B2 (en) * 2006-03-16 2010-08-03 Carl Zeiss Meditec, Inc. Methods for mapping tissue with optical coherence tomography data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特表2003-512085A 2003.04.02

Also Published As

Publication number Publication date
US20080309881A1 (en) 2008-12-18
US8474978B2 (en) 2013-07-02
JP2010529896A (ja) 2010-09-02
WO2008157406A1 (en) 2008-12-24
CA2690834C (en) 2015-08-18
ES2544585T3 (es) 2015-09-01
EP2162054A1 (en) 2010-03-17
EP2162054B1 (en) 2015-05-06
EP2162054A4 (en) 2011-09-07
JP5832744B2 (ja) 2015-12-16
CN101778593A (zh) 2010-07-14
CA2690834A1 (en) 2008-12-24

Similar Documents

Publication Publication Date Title
CN101778593B (zh) 用于分析光学相干断层照相图像的方法
DeLeón-Ortega et al. Discrimination between glaucomatous and nonglaucomatous eyes using quantitative imaging devices and subjective optic nerve head assessment
US7997728B2 (en) Mapping and diagnosis of macular edema by optical coherence tomography
Zangwill et al. A comparison of optical coherence tomography and retinal nerve fiber layer photography for detection of nerve fiber layer damage in glaucoma
Garway-Heath et al. Quantitative evaluation of the optic nerve head in early glaucoma
Mwanza et al. Reproducibility of peripapillary retinal nerve fiber layer thickness and optic nerve head parameters measured with cirrus HD-OCT in glaucomatous eyes
Jeoung et al. Comparison of Cirrus OCT and Stratus OCT on the ability to detect localized retinal nerve fiber layer defects in preperimetric glaucoma
Hu et al. Automated segmentation of neural canal opening and optic cup in 3D spectral optical coherence tomography volumes of the optic nerve head
Bogunović et al. Relationships of retinal structure and humphrey 24-2 visual field thresholds in patients with glaucoma
Lee et al. Ability of Stratus OCT to detect progressive retinal nerve fiber layer atrophy in glaucoma
US9775510B2 (en) Methods to reduce variance in OCT analysis of the macula
Ortega et al. Effect of glaucomatous damage on repeatability of confocal scanning laser ophthalmoscope, scanning laser polarimetry, and optical coherence tomography
JP7159534B2 (ja) 角膜微小層の三次元厚さマッピングおよび角膜診断のための方法およびシステム
Belghith et al. Does the location of Bruch's membrane opening change over time? Longitudinal analysis using San Diego Automated Layer Segmentation Algorithm (SALSA)
Rhodes et al. Variation of laminar depth in normal eyes with age and race
Tun et al. Determinants of optical coherence tomography–derived minimum neuroretinal rim width in a normal chinese population
US11432719B2 (en) Visual field simulation using optical coherence tomography and optical coherence tomographic angiography
Schütze et al. Lesion size detection in geographic atrophy by polarization-sensitive optical coherence tomography and correlation to conventional imaging techniques
Xu et al. Automated assessment of the optic nerve head on stereo disc photographs
Ganekal Ganglion cell complex scan in the early prediction of glaucoma
Weinreb et al. Ophthalmic diagnostic imaging: glaucoma
Duan et al. Evaluation of automated segmentation algorithms for optic nerve head structures in optical coherence tomography images
WO2022175945A1 (en) System and method for monitoring optic disc conditions
Abdellaoui et al. Quantitative assessment of optic disc photographs in normal and open-angle glaucoma patients
Phu et al. Intra‐session repeatability of anterior chamber depth across the chamber width using Pentacam Scheimpflug imaging in healthy subjects

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant