CN105939652B - 具有相位敏感的b扫描配准的光学相干断层扫描(oct)*** - Google Patents

具有相位敏感的b扫描配准的光学相干断层扫描(oct)*** Download PDF

Info

Publication number
CN105939652B
CN105939652B CN201580006926.4A CN201580006926A CN105939652B CN 105939652 B CN105939652 B CN 105939652B CN 201580006926 A CN201580006926 A CN 201580006926A CN 105939652 B CN105939652 B CN 105939652B
Authority
CN
China
Prior art keywords
scan
oct
cluster
offset
phase
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
CN201580006926.4A
Other languages
English (en)
Other versions
CN105939652A (zh
Inventor
杰弗里·P·芬格勒
斯科特·E·弗雷泽
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.)
California Institute of Technology CalTech
University of Southern California USC
Original Assignee
California Institute of Technology CalTech
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 California Institute of Technology CalTech, University of Southern California USC filed Critical California Institute of Technology CalTech
Priority claimed from PCT/US2015/014410 external-priority patent/WO2015120017A1/en
Publication of CN105939652A publication Critical patent/CN105939652A/zh
Application granted granted Critical
Publication of CN105939652B publication Critical patent/CN105939652B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0016—Operational features thereof
    • A61B3/0025—Operational features thereof characterised by electronic signal processing, e.g. eye models
    • 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
    • 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/12—Objective 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/1225—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation
    • A61B3/1233—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation for measuring blood flow, e.g. at the retina
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/0209—Low-coherence interferometers
    • G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B6/00—Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings
    • G02B6/0001—Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings specially adapted for lighting devices or systems
    • G02B6/0005—Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings specially adapted for lighting devices or systems the light guides being of the fibre type
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Ophthalmology & Optometry (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Hematology (AREA)
  • Signal Processing (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

本公开涉及光学相干断层扫描(OCT)领域。本公开尤其涉及OCT***，所述OCT***具有使用相位敏感的B扫描配准方法的配置。在本公开中，OCT***可具有如下配置：扫描物理对象，获取OCT信号来形成B扫描，使用这些B扫描来通过使用B扫描之间的总相位误差确定在轴向方向上的最佳偏移，并且对准B扫描，从而最小化可能在扫描所述物理对象期间发生的运动的影响。

Description

具有相位敏感的B扫描配准的光学相干断层扫描(OCT)***
相关申请的交叉引用
本申请基于2014年2月4日提交的题为“OCT Phase-Sensitive B-scanRegistration Algorithm,”、代理人档案号为CIT-6801-P的美国临时专利申请61/935,431，并且要求所述临时专利申请的优先权。本申请还基于2015年2月4日提交的题为“Optical Coherence Tomography(OCT)with Improved Motion Contrast,”、代理人档案号为064693-0312的专利合作条约(PCT)申请号PCT/US15/14410，并且要求所述PCT申请的优先权。所述临时专利申请和PCT申请的全部内容以引用的方式并入本文。
关于联邦政府资助研究的声明
本发明是在政府支持下进行的(国立卫生研究所(NIH)授予的批准号是NIH STTR1 R41 EY021054)。政府对本发明拥有某些权利。
背景
技术领域
本公开涉及光学相干断层扫描(OCT)领域。本公开尤其涉及具有相位敏感的B扫描配准的OCT***。本公开还尤其涉及用于所述OCT***的相位敏感的B扫描配准方法。
相关技术描述
光学相干断层扫描(OCT)自从在1991年引入以来已经成为不可缺少的临床成像工具。关于OCT技术的背景，参见例如Drexler和Fujimoto等人的“Optical CoherenceTechnology:Technology and Applications”Springer,Heidelberg,Germany,2008。这本书以引用的方式整体并入本文。OCT是基于被称为低相干干涉测量法的光学测量技术。OCT通过将光束引导至物理对象，然后测量并分析反向散射光的量值和时间延迟来执行对物理对象的内部微观结构的高分辨率、横截面成像。
通过执行时间延迟的多次轴向测量(轴向扫描或A扫描)和横向地扫描入射光束来产生横截面成像。这产生A扫描的二维数据集(即，B扫描)，其表示在穿过物理对象的横截平面中的光学反向散射。通过按光栅图案扫描入射光束而获取按顺序的横截面图像，可以产生三维体积数据集(三维OCT或3D-OCT)。这种技术得出物理对象的内部微观结构图像，其具有非常高的清晰度。例如，可在原位并且实时有效地对组织的病理学进行成像，其中分辨率小于15微米。
已经开发了若干类型的OCT***和方法，例如，时域OCT(TD-OCT)和傅立叶域OCT(FD-OCT)。使用FD-OCT允许对视网膜形态进行高分辨率成像，这几乎可以与组织学分析相比较。FD-OCT技术的实例包括频谱域OCT(SD-OCT)和扫频源OCT(SS-OCT)。
OCT可用来识别常见的视网膜血管疾病，诸如年龄相关的黄斑变性(AMD)、糖尿病视网膜病变(DR)以及视网膜血管阻塞。然而，尽管OCT成像快速发展，当前的OCT技术可能无法提供视网膜和脉络膜微脉管***的充分可视化。因此，临床医生经常被迫安排患有视网膜血管疾病的患者进行OCT和荧光素血管造影术(FA)两者。使用在FD-OCT成像期间产生的数据来产生眼底的血管造影图像已经越来越受到关注。这些血管造影片无创地实现而无需注射荧光染料。
近来，已经引入相位方差OCT(PV-OCT)来对视网膜微脉管***进行成像。例如参见：Fingler等人的“Dynamic Motion Contrast and Transverse Flow Estimation UsingOptical Coherence Tomography”的美国专利号7,995,814；Fingler等人的“DynamicMotion Contrast and Transverse Flow Estimation Using Optical CoherenceTomography”的美国专利号8,369,594；Fingler等人的“Mobility and transverse flowvisualization using phase variance contrast with spectral domain opticalcoherence tomography”Opt.Express 2007；15:12636-53；Fingler等人的“Phase-contrast OCT imaging of transverse flows in the mouse retina andchoroid.Invest Ophthalmol.Vis.Sci.2008；49:5055-9；Fingler等人的“Volumetricmicrovascular imaging of human retina using optical coherence tomography witha novel motion contrast technique”Opt.Express[serial online]2009；17:22190-200；Kim等人的“In vivo volumetric imaging of human retinal circulation withphase-variance optical coherence tomography”Biomed Opt Express[serial online]2011；2:1504-13；Kim等人的“Noninvasive imaging of the foveal avascular zonewith high-speed,phase-variance optical coherence tomography”Invest.Ophthalmol.Vis.Sci.2012；53:85-92；以及Kim等人的“Optical imaging of thechorioretinal vasculature in the living human eye”PNAS，2013年8月27日，第110卷，编号35,14354-14359。所有这些公布和专利公开内容以引用的方式整体并入本文。
PV-OCT使用在FD-OCT成像期间正常获取但未使用的数据的软件处理。利用不同于存在于商业器械中的扫描协议，PV-OCT识别在连续的B扫描之间的运动区域，将其与更少的移动区域进行对比。在视网膜和脉络膜中，具有运动的区域对应于脉管***；这些脉管易于与相对静止的其他视网膜组织区分开。
一种用来获取视网膜脉管***的图像的替代性方法是多普勒OCT，其测量在连续深度扫描之间的散射***置变化并且使用这个信息来计算平行于成像方向的流动分量(称作轴向流动)。已经使用多普勒OCT来对视网膜中较大的轴向流动进行成像，但在没有专用扫描协议的情况下，当缓慢流动或取向成横向于成像方向流动时，这种技术受到限制。因为这种技术取决于当成像速度改进继续用于FD-OCT***时测量在连续深度扫描之间的运动变化，所以散射体具有更少时间来在测量之间移动并且最慢的运动因噪声而变得模糊。这进一步减少典型多普勒OCT技术的可视化能力。
相比之下，PV-OCT将能够在增加的FD-OCT成像速度下实现相位测量之间的相同时间分隔，从而保持所展示的能力来独立于脉管取向可视化快速血管流动和缓慢微血管流动。
近年来，若干团体已经开发出OCT成像方法来超越常规的多普勒OCT成像限制。一些方法包括：诸如在2波束扫描中通过FD-OCT机器的硬件修改来增加流动对比，或产生用于提取流动分量的外差频率。其他研究人员已经使用非常规的扫描图案或(诸如PV-OCT中所使用的)重复的B扫描获取来增加相位测量之间的时间分隔并且增强微血管流动的多普勒流动对比。除用来可视化脉管***的基于相位的对比技术之外，基于强度的微脉管***可视化已经被开发用于OCT，其使用分割、基于散斑的时间变化、基于去相关的技术以及基于相位变化和强度变化两者的对比。这些方法中的每一种关于微血管可视化、噪声级以及失真具有变化的能力，同时对在获取期间正在进行典型运动的视网膜组织进行成像。利用体积数据的选择性分割或通过更长成像时间的增加的统计值，可以克服噪声和失真限制中的一些，但需要进一步分析才能够比较来自所有这些不同***的所有可视化能力。
关于OCT方法和***以及其应用的进一步描述例如参见：Schwartz等人的“Phase-Variance Optical Coherence Tomography:A Technique for NoninvasiveAngiography”American Academy of Ophthalmology，第121卷,2014年1月第1期，第180-187页；Sharma等人的“Data Acquisition Methods for Reduced Motion Artifacts andApplications in OCT Angiography”美国专利号8,857,988；Narasimha-Iyer等人的“Systems and Methods for Improved Acquisition of Ophthalmic Optical CoherenceTomography Data”美国专利申请公布号2014/0268046；Everett的“Methods for MappingTissue With Optical Coherence Tomography Data”美国专利号7,768,652。所有这些公布和专利公开内容以引用的方式整体并入本文。
概述
本公开涉及光学相干断层扫描(OCT)领域。本公开尤其涉及具有相位敏感的B扫描配准的OCT***。本公开尤其涉及用于所述OCT***的相位敏感的B扫描配准方法。
例如，所述OCT***可具有如下配置：用具有一定束宽和方向的光束来扫描具有一定表面和深度的物理对象；从所述扫描获取OCT信号；从所获取OCT信号生成离散化A扫描数据和B扫描数据；将所述离散化A扫描数据分配到像素中；使用所获取OCT信号来形成至少一个B扫描集群集，所述至少一个B扫描集群集包括至少一个B扫描集群，所述至少一个B扫描集群包括至少两个B扫描，所述至少两个B扫描相互平行并且形成平行于光束方向的平面；以及由所述至少两个B扫描形成B扫描对。
所述OCT***还可具有如下配置：(a)计算每个B扫描对的总相位误差并且在平行于光束方向的方向上将每个B扫描对中的一个B扫描偏移至少两次，或(b)在平行于光束方向的方向上将每个B扫描对中的一个B扫描偏移至少三次；其中每次偏移是一个像素并且平行于光束。
所述OCT***还可具有如下配置：在每次偏移之后计算每个B扫描对的总相位误差；识别每个B扫描对的所有计算出的总相位误差中最小计算出的总相位误差；识别每个B扫描对的偏移量，所述偏移量导致所识别的最小计算出的总相位误差；计算所述至少一个B扫描集群内的每个B扫描的累积轴向偏移；以及对准所述至少一个B扫描集群内的所有B扫描。
所述OCT***还可具有如下配置：在形成B扫描对之后，计算每个B扫描对的相位差；并且其中可使用计算出的相位差来计算每个总相位误差，其中每个计算出的相位差可构成相位差数据点。所述相位差数据点可用来计算整体运动(bulk motion)相位差。
可针对整体运动的影响对计算出的整体运动相位差进行校正以便提供已校正的相位差。
可将基于强度的阈值方法应用于所述已校正的相位差以便降低噪声的影响。
所述OCT***可具有使用已对准的B扫描来形成所述物理对象的图像的配置。所述物理对象可以是人类组织。
所述OCT***可包括：至少一个光源，其提供至少一个光束；至少一个回射反射器；至少一个光纤耦合器或至少一个自由空间耦合器，其将所述至少一个光束导引至所述物理对象和至少一个回射反射器。导引至所述物理对象的所述至少一个光束可形成至少一个反向散射光束。导引至所述至少一个回射反射器的所述至少一个光束可形成至少一个反射的参考光束。
所述OCT***还可包括：至少一个扫描光学器件，其在所述物理对象上扫描所述至少一个光束；至少一个检测器。所述检测器可组合所述至少一个反向散射光束与所述至少一个反射光束来形成光干涉，检测所述至少一个反向散射光束的量值和时间延迟，并且形成OCT信号。所述至少一个光纤耦合器或所述至少一个自由空间耦合器可将所述至少一个反向散射光束和所述至少一个反射光束导引至所述至少一个检测器。
所述OCT***还可包括：至少一个处理器，其获得并分析由所述至少一个检测器形成的OCT信号，并且形成所述物理对象的至少一个图像；以及至少一个显示器，其显示所述物理对象的所述至少一个图像。
所述OCT***可具有基于B扫描之间的强度或相位变化来识别运动区域的配置。所述OCT***可具有如下配置：例如通过使用相位方差OCT(PV-OCT)方法、相位对比OCT(PC-OCT)方法、强度/散斑方差OCT(IV-OCT)方法、多普勒OCT(D-OCT)方法、多普勒功率频移OCT(PDS-OCT)方法、频谱分离幅度去相关分析(SSADA)方法、光学微血管造影(OMAG)方法、相关映射OCT(cmOCT)方法或其组合来识别所述运动区域。所述OCT***可具有可使用相位方差OCT(PV-OCT)方法的配置。
一种包含指令程序的非暂态、有形计算机可读存储介质(“存储介质”)也在本公开的范围内，所述指令程序可致使运行所述指令程序的计算机***充当光学相干断层扫描(OCT)***。
所述存储介质可具有如下配置：用具有一定束宽和方向的光束来扫描具有一定表面和深度的物理对象；从所述扫描获取OCT信号；从所获取OCT信号生成离散化A扫描数据和B扫描数据；将所述离散化A扫描数据分配到像素中；使用所获取OCT信号来形成至少一个B扫描集群集，所述至少一个B扫描集群集各自包括至少一个B扫描集群，所述至少一个B扫描集群包括至少两个B扫描，所述至少两个B扫描相互平行并且形成平行于光束方向的平面；以及由每个B扫描集群中的所述至少两个B扫描形成B扫描对。所述存储介质还可具有进行以下操作中的任一个的配置：(a)在任何B扫描的任何偏移之前计算每个B扫描对的总相位误差，然后将每个B扫描对中的一个B扫描偏移至少两次，或(b)将每个B扫描对中的一个B扫描偏移至少三次。每次偏移可以是一个像素，并且是在平行于所述光束的方向上。所述存储介质还可具有如下配置：在每次偏移之后计算每个B扫描对的总相位误差；识别每个B扫描对的所有计算出的总相位误差中最小的计算出的总相位误差；识别每个B扫描对的偏移量，所述偏移量导致所识别的最小计算出的总相位误差；计算所述至少一个B扫描集群内的每个B扫描的累积轴向偏移；并且对准所述至少一个B扫描集群内的所有B扫描。
所述指令程序可致使运行所述指令程序的计算机***进行以下操作：获得并分析由至少一个检测器形成的OCT信号；并且形成所述物理对象的至少一个图像；并且显示所述物理对象的所述至少一个图像。
上述***、存储介质和方法的任意组合在本公开的范围内。
这些以及其他部件、步骤、特征、对象、益处以及优点现在将从阅读以下说明性实施方案的详述、附图以及权利要求变得显而易见。
附图简述
附图是说明性实施方案的图。附图并不示出所有实施方案。此外或作为替代，可使用其他实施方案。可以省略可能是显而易见的或不必要的细节，以便节省空间或用于更有效的说明。可在具有另外的部件或步骤的情况下且/或在没有所示出的所有部件或步骤的情况下实践一些实施方案。当相同的数字出现在不同附图中时，它是指相同或相似的部件或步骤。
图1是一般化的OCT***的示意图。
图2示意性地示出用于OCT***的扫描配置。
图3示意性地示出人类左眼的矢状图。
图4示意性地示出视网膜的横截面层。
图5示出视网膜中心凹区域的横截面OCT图像。
图6示出横向于传播方向的光束的强度分布的实例。
图7通过举例的方式示意性地示出可用于计算OCT相位敏感的B扫描配准的四个B扫描、两个B扫描集群以及一个B扫描集群集。
图8示出可如何通过使用相位敏感的B扫描配准方法来校正OCT图像。
图9示意性地示出可如何计算总相位误差的示例性方法。
图10示出由于在z轴方向上的运动而引起的第二B扫描的偏移的实例。
图11示出确定相位差最小值和最佳偏移的的实例。
示例性实施方案的详述
现在描述说明性实施方案。此外或作为替代，可使用其他实施方案。可省略可能是显而易见的或不必要的细节，以便节省空间或用于更有效的说明。可在具有另外的部件或步骤的情况下且/或在没有所描述的所有部件或步骤的情况下实践一些实施方案。
已经论述的部件、步骤、特征、目的、益处和优点仅是说明性的。其中的任何一个或与其相关的论述均不意图以任何方式限制保护范围。也涵盖许多其他实施方案。这些实施方案包括具有更少的、另外的和/或不同的部件、步骤、特征、目的、益处和优点的实施方案。这些实施方案还包括将部件和/或步骤不同地布置和/或排序的实施方案。
本公开涉及光学相干断层扫描(OCT)领域。本公开尤其涉及OCT***，所述OCT***具有使用相位敏感的B扫描配准方法的配置。本公开尤其涉及用于所述OCT***的相位敏感的B扫描配准方法。
本公开涉及OCT***。所述OCT***可包括具有任何光学设计的任何干涉仪，包括诸如迈克尔森干涉仪、马赫-增德尔干涉仪、吉莱-图努瓦干涉仪、基于共同路径的设计或其他干涉仪体系结构。干涉仪中的样品臂和参考臂可由任何类型的光学器件(例如，块状光学器件、光纤、混合块状光学***等)组成。
所述OCT***还可包括任何基础OCT***。基础OCT***的实例可包括时域OCT(TD-OCT)和傅立叶域或频域OCT(FD-OCT)。FD-OCT的实例可包括频谱域OCT(SD-OCT)、扫频源OCT(SS-OCT)和光频域成像(OFDI)。
所述OCT***可具有使用识别且/或可视化运动区域的任何OCT方法(“OCT运动对比方法”)的配置。所述OCT运动对比方法可使用发生在物理对象内的运动来识别且/或可视化具有改进的对比的区域。例如，OCT可以使用由血管中的血流引起的OCT信号变化来通过重要的图像对比识别且/或可视化脉络膜脉管***。因此，获得更好的图像，且/或以前通过典型OCT***不能被识别的结构变得可见。例如，通过使用OCT运动对比方法，脉络膜血管层可变得可见。所述OCT运动对比方法的实例可包括相位方差OCT(PV-OCT)、相位对比OCT(PC-OCT)、强度/散斑方差OCT(IV-OCT)、多普勒OCT(D-OCT)、多普勒功率频移OCT(PDS-OCT)、频谱分离幅度去相关分析(SSADA)、光学微血管造影(OMAG)、相关映射OCT(cmOCT)等。所述PV-OCT方法的实例由Fingler等人的“Dynamic Motion Contrast and Transverse FlowEstimation Using Optical Coherence Tomography”的美国专利号7,995,814；Fingler等人的“Dynamic Motion Contrast and Transverse Flow Estimation Using OpticalCoherence Tomography”的美国专利号8,369,594；Fingler等人的“Mobility andtransverse flow visualization using phase variance contrast with spectraldomain optical coherence tomography”Opt.Express[serial online]2007；15:12636-53公开；所述散斑方差OCT方法的实例由Mariampillai等人的“Speckle variancedetection of microvasculature using swept-source optical coherencetomography,”Opt.Lett.33(13)，1530-1532(2008)公开；所述相关映射OCT方法的实例由Enfield等人的"In vivo imaging of the microcirculation of the volar forearmusing correlation mapping optical coherence tomography(cmOCT)"Biomed.Opt.Express 2,1184-1193(2011)公开；所述OMAG方法的实例由An等人的"In vivovolumetric imaging of vascular perfusion within human retina and choroidswith optical micro-angiography"Opt.Express 16,11438-11452(2008)公开；所述功率多普勒OCT方法的实例由Makita等人的"Optical coherence angiography"Opt.Express14,7821-7840(2006)公开；所述SSADA方法的实例由Jia等人的“Split-spectrumamplitude-decorrelation angiography with optical coherence tomography,”Opt.Express20(4)，4710-4725(2012)公开。这些公开的全部内容以引用的方式并入本文。
OCT***100可包括至少一个光源110、至少一个扫描光学器件200、至少一个回射反射器180、至少一个光纤耦合器220或至少一个自由空间耦合器、至少一个检测器130、至少一个处理单元140以及至少一个显示单元150。
在图1中示意性地示出的一般化的OCT***的实例由Fingler等人的“DynamicMotion Contrast and Transverse Flow Estimation Using Optical CoherenceTomography”的美国专利号7,995,814；Fingler等人的“Dynamic Motion Contrast andTransverse Flow Estimation Using Optical Coherence Tomography”的美国专利号8,369,594；以及Sharma等人的题为“Data Acquisition Methods for Reduced MotionArtifacts and Applications in OCT Angiography”的美国专利号8,857,988公开。这些公开的内容以引用的方式整体并入本文。OCT***100可包括这种一般化的OCT***。
至少一个光源110可包括任何光源，例如，低相干光源。可通常通过使用至少一个光纤220来导引来自光源110的光来照射物理对象210。物理对象210的实例可以是人眼中的任何组织。例如，所述组织可以是视网膜。光源110可以是在SD-OCT情况下具有较短时间相干性长度的宽带低相干光源或在SS-OCT情况下的波长可调谐激光源。可通常用处于光纤220的输出与物理对象210之间的扫描光学器件200来扫描光，使得被导引用于物理对象210的光束(虚线)在将要成像的面积或体积上被横向地(在x轴和/或y轴上)扫描。扫描光学器件200可包括适合于扫描的任何光学元件。扫描光学器件200可包括至少一个部件。扫描光学器件200的至少一个部件可以是光学部件。可通常将从物理对象210散射的光收集到用来导引光来照射物理对象210的相同光纤220中。(在图1中示出物理对象210，仅仅是为了示意性地展示物理对象210与OCT***100相关。物理对象210不是OCT***100的部件。)
OCT***100还可包括光束分离器120，其用来分离由光源110提供的光并且将光导引至参考臂230和物理对象臂240。OCT***还可包括放在光束分离器120与回射反射器180之间的透镜160。OCT***还可包括放在光束分离器120与扫描光学器件200之间的透镜170。
源自同一光源110的参考光250可沿着单独路径行进，所述路径在这种情况下包括光纤220和具有可调光学延迟的回射反射器180。回射反射器180可包括至少一个部件。回射反射器180的至少一个部件可以是光学部件，例如参考镜。还可使用透射性参考路径，并且可调延迟可放在干涉仪100的物理对象臂240或参考臂230中。
可通常在光纤耦合器中组合从物理对象210散射的所收集光260与参考光250以便在检测器130中形成光干涉，从而形成OCT信号。尽管示出通往检测器130的单个光纤端口，但是干涉仪的各种设计可用于对用于SS-OCT的干涉信号或用于SD-OCT的分光计检测器进行平衡或不平衡检测。
可将来自检测器130的输出供应至处理器140。可将结果存储在处理器140中或显示在显示器150上。处理功能和存储功能可以位于OCT***内，或功能可在外部处理单元上执行，所收集的数据被传送到所述外部处理单元。这种外部单元可专用于数据处理或执行相当通用且不专用于OCT***的其他任务。
如本文所使用的光束应被解释为任何仔细地引导的光路径。在时域***中，参考臂230可能需要具有可调谐的光学延迟来产生干涉。在TD-OCT和SS-OCT***中可通常使用平衡检测***，而在用于SD-OCT***的检测端口处可使用分光计。
干涉可致使干涉光的强度在光谱上变化。干涉光的傅里叶变换可揭示散射强度在不同路径长度下的轮廓，并且因此随物理对象中的深度(z轴方向)而变的散射。例如参见Leitgeb等人的“Ultrahigh resolution Fourier domain optical coherencetomography,”Optics Express 12(10):2156，2004。这个公布的全部内容以引用的方式并入本文。
随深度而变的散射的轮廓被称为轴向扫描(A扫描)，如在图2中示意性地示出。在物理对象中的相邻位置处测量的A扫描的集合产生物理对象的横截面图像(断层扫描或B扫描)。在样品上的不同横向位置处收集的单独B扫描的集合构成数据体或立方体。通过组合多个B扫描，可以形成三维C扫描。对于特定的数据体，术语快轴是指沿单个B扫描的扫描方向，而慢轴是指收集多个B扫描所沿的轴线。
B扫描可由在x轴和y轴所表示的平面中的任何横向扫描形成。B扫描可例如沿水平或x轴方向、沿垂直或y轴方向、沿x轴方向和y轴方向的对角方向、按圆形或螺旋形图案及其组合形成。本文所论述的大多数实例可指在x-z轴方向上的B扫描，但本公开可同样地适用于任何横截面图像。
物理对象210可以是任何物理对象。物理对象210可以是如在图3中以简化方式示出的人眼500。人眼包括眼角膜510、瞳孔520、视网膜300、脉络膜540、视网膜中心凹区域550、视神经盘560、视神经570、玻璃体腔580和视网膜血管590。
物理对象210可以是组织。所述组织的实例是视网膜。在图4中示意性地示出视网膜300的各层的简化横截面图像。视网膜层包括神经纤维层(NFL)310、外界膜(ELM)320、光感受器内/外段330、光感受器外段340、视网膜色素上皮细胞(RPE)350、视网膜色素上皮细胞(RPE)/布鲁赫膜复合物360。图4还示意性地示出视网膜中心凹370。图5示出视网膜中心凹区域的横截面OCT图像。
所述OCT***可包括提供至少一个OCT***的方法，所述至少一个OCT***用来通过用至少一个光束扫描物理对象来获取数据，以形成至少一个A扫描和至少一个B扫描。所述至少一个OCT***可以是如上文所公开的任何OCT***。
物理对象可包括如上文所公开的任何物理对象。物理对象具有表面和深度。例如，眼睛的眼底具有通过瞳孔从外侧环境接收光的外表面。眼睛的眼底还具有在其外表面处开始并且从其外表面延伸的深度。
在本公开中，z轴是平行于延伸到物理对象的深度中的光束的轴线(“轴向轴线”)，x轴和y轴(“横向轴线”)是横向的，从而是垂直于z轴的轴线。在图1-5和图7中示出这三个轴线的取向。
由OCT***提供的至少一个光束在物理对象的位置处具有宽度和强度。在一个实例中，在物理对象的这个位置处，光束被聚焦(“聚焦的光束”)。例如，在这个位置处，光束的宽度处于其最小值。光束的横截面积可具有任何形状。例如，所述横截面积可具有圆形形状或椭圆形形状。聚焦的光束的强度沿其横向轴线变化，所述横向轴线垂直于其传播轴线。这个横向光束轴线可以是径向轴线。在光束的中心处的光束强度处于其峰值(即，光束强度处于其最大值)，并且沿其横向轴线降低，从而形成强度分布。可通过高斯函数来大致估计这种分布，如图6中所示。光束的宽度(“束宽”)被定义为在两个相对点处与强度分布相交的线的长度，在所述两个相对点处，强度是其峰值的1/e2倍。光束可包括多于一个的峰值。使用具有最高光束强度的峰值来计算束宽。在图6中示意性地示出束宽的实例。典型OCT***的典型束宽在物理对象位置处可在10微米到30微米范围中变化。
OCT***可使用还可包括获取数据来形成至少一个B扫描集群集的方法。至少一个B扫描集群集的数目P等于或大于1，其中P是整数。例如，P可以是1、2、3、4、5、10、100、1,000、10,000或100,000。
每个B扫描集群集可包括任意数量的B扫描集群，N等于或大于2，其中N是整数。例如，N可以是2、3、4、5、10、100、1,000、10,000或100,000。
每个B扫描集群可包括任意数量的B扫描，M等于或大于2，其中M是整数。例如，M可以是2、3、4、5、10、20、100、1,000、10,000或100,000。
每个B扫描可沿相同的横向轴线或另一个横向轴线(x轴或y轴)来定位，所述横向轴线可平行于B扫描集群集内的其他B扫描的横向轴线。每个B扫描可形成垂直于所述横向轴线中的一个的平面，并且每个B扫描平面可因此平行于另一个B扫描的平面。也就是说，每个B扫描可平行于z轴。
每个B扫描包括例如位于(x-z)平面上的多个数据点。在一段时间内获取每个B扫描、每个B扫描集群以及每个B扫描集群集。也就是说，每个B扫描、每个B扫描集群以及每个B扫描集群集是在分别与所有其他B扫描、所有其他B扫描集群以及所有其他B扫描集群集不同的时间处形成。在本公开中，“首先形成”意思是在时间上首先形成；“第二个形成”意思是在时间上第二个形成；“接下来形成”或“邻近的”意思是在时间上接下来形成；并且“最后形成”意思是在时间上最后形成。
可针对在一段时间内获取的每个(x,z)数据点来计算OCT运动对比。可根据复杂的OCT信号、OCT强度信息、相位信息或其组合来计算这个运动对比。
在每个B扫描集群内的每个B扫描平面之间的空间距离(集群内距离)可在0到一半束宽(以微米为单位)的范围中变化。例如，所述集群内距离可在0微米到15微米的范围中变化。
在每个B扫描集群的最后形成的B扫描与邻近于所述B扫描的另一个接下来的B扫描集群的首先形成的B扫描之间的空间距离(“集群间距离”)可至少等于或大于1微米。例如，所述集群内距离可在1微米到10微米、1微米到100微米或1微米到1,000微米的范围中变化。
在图7中示出的实例中，首先形成的B扫描集群包括两个B扫描并且第二个形成的B扫描集群包括两个B扫描。在这个实例中，所述B扫描集群集包括两个B扫描集群。
所述方法还可包括获取数据来形成至少两个B扫描集群集。每个B扫描集群集可沿相同的横向轴线或另一个横向轴线(x轴或y轴)来定位，所述横向轴线可平行于其他B扫描集群集的横向轴线。在每个B扫描集群集的最后形成的B扫描与邻近于所述B扫描集群集的另一个B扫描集群集的首先形成的B扫描之间的空间距离(“集群集间距离”)可至少等于或大于20微米。
虽然在OCT图像获取期间的样本运动发生在三维中，但是仅可容易地识别两个按顺序的B扫描之间的轴向运动。在获取两个按顺序的B扫描期间发生的横向运动可能会在所测量的OCT信号中引起不可被直接识别的散斑噪声。此外，尽管可通过经由基于OCT强度的对准方法对B扫描进行对准和再配准来最小化所述轴向运动的影响，但是因为B扫描未被最佳地对准，所以这类方法可通过B扫描向OCT图像对比引入散斑和/或相位噪声。
本公开涉及用于最小化样本运动对OCT图像对比的影响的方法(“相位敏感的B扫描配准”)。相位敏感的B扫描配准方法可在z轴方向上最小化样本运动对OCT图像对比的影响。本公开还涉及用于最小化在两个按顺序的B扫描之间的样本运动的影响的方法。OCT信号可包括OCT信号强度和/或相位。
相位敏感的B扫描配准方法可包括对准运动的轴向(z轴方向)分量以最小化运动对OCT信号的影响(“轴向对准方法”)。然后，剩余的相位噪声可与可能已经发生的横向运动直接相关。
本公开涉及对准所有B扫描来针对轴向运动进行校正的光学相干断层扫描(OCT)***。已对准的所有B扫描可用来提供已针对轴向运动进行校正的OCT图像。
用于轴向对准的方法可包括在两个按顺序B扫描之间的基于相位的轴向对准方法。这种方法还可包括：
(a)最小化经由在两个B扫描之间的轴向运动引入的相位噪声，和/或
(b)确定可用来评估可在B扫描期间发生的横向运动的量的定量度量。这种方法可允许改进的处理选项和用于将数据与具有最大运动噪声的B扫描隔离的统计异常值分析。
相位敏感的B扫描配准方法可包括通过使用OCT***来获取数据以形成至少两个B扫描。两个邻近的B扫描可形成一对(“首先形成的B扫描和第二个形成的B扫描”或“B扫描对”)。B扫描对的每个B扫描可包括一定数目的A扫描Q，其中Q是等于或大于1的整数。例如，Q可以是1、10、100、1,000或10,000。B扫描对的两个B扫描可具有相同数目的A扫描。
相位敏感的B扫描配准方法还可包括：通过将B扫描对中的一个相对于另一个B扫描沿z轴方向偏移预定距离来计算在首先形成的B扫描与第二个形成的B扫描之间的至少三个相位差，使用这三个相位差来计算针对每个差的总相位误差，直到得到总相位误差的最小值(“总相位误差最小值”)为止。预定距离可以是一个像素。也就是说，偏移可以是(大约)一个像素。
例如，OCT***可具有如下配置：用具有一定束宽和方向的光束来扫描具有一定表面和深度的物理对象；从所述扫描获取OCT信号；从所获取OCT信号形成像素；使用所获取OCT信号来形成至少一个A扫描；使用所获取OCT信号来形成至少一个B扫描集群集，所述至少一个B扫描集群集包括至少一个B扫描集群，所述至少一个B扫描集群包括至少两个B扫描，所述至少两个B扫描相互平行并且形成平行于光束方向的平面；并且由所述至少两个B扫描形成B扫描对。
OCT***还可具有如下配置：(a)计算每个B扫描对的总相位误差并且在光束的方向上将每个B扫描对中的一个B扫描偏移至少两次，或(b)在光束的方向上将每个B扫描对中的一个B扫描偏移至少三次；其中每次偏移是一个像素并且平行于光束。
OCT***还可具有如下配置：在每次偏移之后计算每个B扫描对的总相位误差；确定每个B扫描对的最小总相位误差；确定每个B扫描对的最佳偏移；计算每个B扫描在所述偏移结束时已经偏移的总距离；并且对准所述至少一个B扫描集群内的所有B扫描。
在本公开中，“相位差最小值”是“计算出的最小总相位误差”；最佳偏移是在所述相位差最小值处识别的偏移的量。
在图8中以简化方式示意性地示出这个***可如何用来校正OCT图像的实例。
实例1.通过使用计算出的相位差来计算总相位误差。
可通过使用任何方法来计算总相位误差。可通过使用计算出的相位差来计算总相位误差。例如，可通过使用以下方程来计算总相位误差。
其中在方程1中，是在属于首先形成的B扫描的(x，z)位置处的数据点的相位，是在属于第二个形成的B扫描的(x，z+Δz)位置处的数据点的相位，所述第二个形成的B扫描轴向偏移了距离Δz距离，mod是在0与2π之间的相位差的模，并且是在首先形成的B扫描的数据点的相位与第二个形成的B扫描的数据点的相位之间的相位差。
每个数据点的计算出的值提供相位差数据点。这些相位差数据点用来通过使用现有技术中可用的任何方法计算整体运动相位差例如，可使用由Makita等人的“Optical coherence angiography”Opt.Express 14,7821-7840(2006)所公开的整体运动相位差计算方法。这个公开的内容以引用的方式整体并入本文。
可通过使用以下方程针对整体运动的影响进行校正：
其中在方程2中，是针对对应于给定x坐标x0的每个单独A扫描的已校正的相位差。
然后，通过使用现有技术中可用的任何方法来将基于强度的阈值方法应用于已校正的相位差。例如，由Fingler等人公开的阈值方法。本申请还基于2015年2月4日提交的题为“Optical Coherence Tomography(OCT)with Improved Motion Contrast,”、代理人档案号为064693-0312的专利合作条约(PCT)申请号PCT/US15/14410，并且主张所述PCT申请的优先权。这个公开的全部内容以引用的方式并入本文。
在这里，是用阈值方法校正过的相位差。
最后，可通过使用以下方程来计算总相位误差。
方程式方程3。
图9中示意性地示出这种方法。
实例2.得到总相位误差最小值。
可通过任何方法得到这些相位差。例如，可通过使用由OCT***获取的首先形成的B扫描和第二个形成的B扫描来计算第一相位差，而无需将B扫描对相对于彼此进行偏移。也就是说，可针对无偏移来计算第一相位差。
可通过首先将第二个形成的B扫描相对于首先形成的B扫描沿z轴偏移来计算其他相位差。偏移可以是一个像素。例如，计算相位差可包括选择对B扫描对中的一个(例如，首先形成的B扫描)进行的偏移数T，其中T是等于或大于3的整数。例如，T可以是3次偏移、4次偏移、5次偏移、10次偏移或100次偏移。
可沿z轴在任何方向上进行偏移。一个这种方向在下文被称为“正方向”，并且与正方向相反的另一个方向在下文被称为“负方向”。例如，可通过以下操作来计算其他相位差：首先将第二个形成的B扫描相对于首先形成的B扫描在正方向上沿z轴一个像素偏移接一个像素偏移进行偏移，并且然后计算针对每次偏移在B扫描对之间的相位差。
这些计算出的相位差中的每一个用来计算针对每次偏移的总相位误差。可例如通过使用实例1中公开的方法来计算针对每次偏移的总相位误差。如果在每次偏移之后总相位误差降低，那么可将第二B扫描沿正方向一个像素接一个像素进行进一步偏移，直到总相位误差增加。从而得到总相位误差的最小值。
如果在正方向上的一次偏移之后总相位误差增加，那么可通过以下操作来计算其他相位误差：将第二B扫描相对于首先形成的B扫描在负方向上沿z轴一个像素接一个像素进行偏移，并且然后计算针对每次偏移在B扫描对之间的总相位差。如果在每次偏移之后总相位误差降低，那么可将第二B扫描沿负方向一个像素接一个像素进行进一步偏移，直到总相位误差增加。从而得到总相位误差的最小值。
通过举例的方式在图10-11中示出这种方法。在这种示范性实例中，由于在获取首先形成的B扫描和第二个形成的B扫描时发生的运动影响，已经将(以简化方式示出的)三个视网膜层的位置在负方向(例如，向后的方向)上偏移了2个像素。参见图10。为了针对这种运动进行校正，针对在正方向(例如，向前的方向)上无偏移、1像素偏移、2像素偏移、3像素偏移以及4像素偏移首先计算在首先形成的B扫描与第二个形成的B扫描之间的总相位误差。然后，因此根据这些计算将总相位误差的位置确定为正方向上的2个像素，如图11中所示。用于这个实例的偏移量是2个像素。
实例3.计算B扫描集群内的每个B扫描的累积轴向偏移。
可通过使用任何方法来计算B扫描集群内的每个B扫描的累积轴向偏移。可选择参考B扫描并且可将集群内的其他B扫描与这个参考B扫描对准。这个参考B扫描可以是B扫描集群内的任何B扫描。例如，参考B扫描可以是属于B扫描集群的首先形成的B扫描。参考这个首先形成的B扫描对以下实例进行解释。然而，可选择任何其他B扫描来实现相同的最终结果。
在这种方法种，例如，相互邻近(即，相邻)的每两个B扫描可形成一个B扫描对。参考B扫描可以是属于B扫描集群的首先形成的B扫描。可以实例2中所公开的方式确定每个B扫描对的总相位误差的最小值的位置。这是在形成B扫描对的B扫描之间的最佳偏移Δz最佳，j，如以下方程所定义：
Δz最佳，j＝在B扫描j-1与B扫描j之间的最佳偏移 方程4。
其中Δz最佳，j是对应于最小总相位误差的Δz。
B扫描集群内的每个B扫描的累积轴向偏移可通过计算每个B扫描对的Δz最佳，j并且通过使用以下方程将它们相加来确定：
Δz累积，j＝B扫描j的累积轴向偏移 方程5。
实例4.
针对7个相对轴向位置对两个B扫描进行比较。可使用布尔强度阈值，所述布尔强度阈值基于针对这些对准位置中的每一个的平均强度图像，从而选择仅比单独OCT强度图像的噪声级高一点的阈值。针对这些轴向再对准位置中的每一个(例如，针对每个z位置)，通过计算两个A扫描之间的相位变化来对所述两个A扫描进行比较，并且可针对这个A扫描对计算整体轴向运动，并且基于相位卷绕使用2π的模来校正并调整所述相位变化。一旦针对每个A扫描对计算了整体轴向运动，就可在整个B扫描上对阈值限定的、已校正的相位变化的绝对值求和。可针对每个轴向变化取向来计算这些总和，并且按这组的最高值进行归一化，从而产生归一化的相位总和。可导致归一化的相位总和的最小值的轴向对准可以是最佳轴向对准以最小化所述相位噪声，并且在所述对准处的最小相位总和值可与在获取期间已经发生的横向运动相关联。
实例5.
如果轴向再对准是显著的，那么有可能的是，针对B扫描的一部分的最佳轴向对准可能与其他部分不同。这种处理算法可单独在整个B扫描的各部分上使用，以总体上验证用于所有部分的相同配准，并且识别需要另外的处理和注意来解决样本轴向运动的这种可变性的情况。
实例6.
在这个实例中，提高处理速度。在上述轴向再对准方法中，针对R轴向取向执行B扫描到B扫描的相位变化计算。不是默认地计算所有这些取向，可通过仅初始地计算较少的情况(诸如针对轴向变化Δz＝R、0、1和-1，其中R为整数)来实现速度提高。一般来讲，Δz＝1、0、-1可给予足够的数据来识别最佳轴向再对准为0的情况，并且在其他情况下可识别需要继续进行哪些另外的计算来实现最佳对准。如果所述集合的最小度量是在Δz＝1处，那么可以计算Δz＝2来识别哪一个实际上可以是最小值或是否可能需要另外的计算。在Δz＝R处的情况应该提供足够的去相关来对用来评估横向运动噪声的度量进行归一化。
实例7.
代替对两个B扫描内的每个A扫描执行轴向再对准和相位校正，可选取对A扫描的选择，来作为在整个B扫描上展开的代表性集合。对于这种变化，可以多种方式实现最佳轴向再对准，包括但不限于：(a)单独地比较A扫描与A扫描轴向再对准，或(b)将这个所选取的选择内的所有A扫描一起进行再对准，以识别最佳轴向对准，并且从而减少由引发相位误差的特征(诸如主要的脉管***)所引发的误差，这将影响单独的A扫描比较中的一些。
上文公开的OCT方法可用于任何OCT相关的应用。例如，这种方法可用于形成具有物理对象的更大视野OCT图像。可将这种方法并入与基于OCT的血管造影术相关的方法和***中。例如，可通过使用OCT运动对比方法来更详细地识别脉络膜脉管***。包括OCT运动对比方法的OCT方法还可用于诊断和/或治疗健康状况(诸如疾病)。例如，包括OCT运动对比方法的OCT方法可用于表征视网膜健康。
上文公开的OCT***可提供述物理对象相关的任何信息。例如，可使用运动对比方法的这个***可提供2D(即，横截面)图像、横断面图像、3D图像、与健康状况相关的度量等等。这种***可与任何其他***一起使用。例如，OCT***可与用于诊断或治疗目的的超声波装置或外科***一起使用。OCT***可用来分析任何物理对象。例如，OCT***可用于分析(例如)任何类型的生命形式和无生命对象的(例如)图像的形成。生命形式的实例可以是动物、植物、细胞等。
除非另外指出，否则本文中已论述的处理单元140可用计算机***来实现，所述计算机***被配置来执行本文中针对这个单元已论述的功能。所述计算机***包括一个或多个处理器、有形存储器(例如，随机存取存储器(RAM)、只读存储器(ROM)和/或可编程只读存储器(PROM))、有形存储装置(例如，硬盘驱动器、CD/DVD驱动器和/或闪存存储器)、***总线、视频处理部件、网络通信部件、输入/输出端口和/或用户接口装置(例如，键盘、指向装置、显示器、麦克风、声音再现***和/或触摸屏)。
用于所述处理单元140的计算机***可包括在相同或不同位置处的一个或多个计算机。当在不同位置处时，计算机可被配置来通过有线和/或无线网络通信***相互通信。
计算机***可包括软件(例如，一个或多个操作***、装置驱动程序、应用程序和/或通信程序)。当包括软件时，软件包括编程指令并且可包括相关联的数据和库。当包括编程指令时，编程指令被配置来实现一种或多种算法，所述一种或多种算法实现所述计算机***的功能中的一个或多个，如本文所述。由每个计算机***执行的每个功能的描述也构成执行所述功能的算法的描述。
软件可存储在一个或多个非暂态、有形存储装置(诸如一个或多个硬盘驱动器、CD、DVD和/或闪存存储器)上或中。软件可以是源代码和/或目标代码格式。相关联的数据可存储在任何类型的易失性和/或非易失性存储器中。软件可被加载到非暂态存储器中并且可由一个或多个处理器执行。
上文公开的方法、装置、***和特征的任意组合在本公开的范围内。
除非另外说明，否则在本说明书中(包括在以下权利要求书中)所阐述的所有测量、值、等级、位置、量值、大小及其他规范均为近似的，而不是准确的。它们意图具有合理的范围，这个范围与它们的有关功能一致并且与它们所属领域中的习惯一致。
在本公开中提到的所有物品、专利、专利申请和其他公布以引用的方式并入本文。
在本公开中，不定冠词“一个”和短语“一个或多个”以及“至少一个”是同义的并且意思是“至少一个”。
短语“用于...的装置”当在权利要求书中使用时意图并且应该被解释为涵盖已经描述的对应结构和材料和它们的等同物。类似地，短语“用于...的步骤”当在权利要求书中使用时意图并且应该被解释为涵盖已经描述的对应动作和它们的等同物。权利要求中不存在这些短语意思是所述权利要求并不意图并且不应该被解释为受限于这些对应结构、材料或动作或它们的等同物。
保护范围仅受到以下权利要求书的限制。当根据本说明书和以下申请过程来解释时，所述范围意图并且应该被解释为与权利要求书中所用语言的普通意思的宽度相一致，除非特定意思已经被阐述并且用来涵盖所有结构和功能等同物。
诸如“第一”和“第二”等关系性术语仅仅可以用来将一个实体或动作与另一个实体或动作区分开，而没有必要要求或暗示它们之间的任何实际关系或顺序。术语“包括”和其任何其他变型在结合说明书或权利要求书中的元素列表使用时意图指示所述列表不是排他的并且可包括其他元素。类似地，在没有进一步约束的情况下，前面有“一个”的元素并不排除存在相同类型的另外的元素。
权利要求中没有一项意图涵盖不能满足专利行为的101部分、102部分或103部分的要求的主题，也没有一项应该以这种方式被解释。因此放弃对这种主题的任何非有意的涵盖。除了如这一段中刚刚所述，已陈述或说明的任何内容都不意图或不应该被解释为致使任何部件、步骤、特征、目标、益处、优点或等同物献给公众，无论它是否在权利要求书中有叙述。
提供摘要来帮助读者快速确定技术性公开的本质。应理解，所提交的摘要不是用来解释或限制权利要求书的范围或含义。此外，在上文详细描述中的各种特征在各种实施方案中被分组在一起以精简本公开。这种公开的方法不应被解释为需要所要求保护的实施方案需要比每个权利要求中所明确表述的特征更多的特征。相反，如以下权利要求书所反映的，本发明的主题存在于单个所公开实施方案的少于所有特征中。因此，以下权利要求书因此并入详述中，其中每项权利要求自身可作为单独要求保护的主题。

Claims (16)

1.一种光学相干断层扫描(OCT)***，其具有如下配置：
用具有一定束宽和方向的光束来扫描具有一定表面和深度的物理对象；
从所述扫描获取OCT信号；
从所获取OCT信号生成离散化A扫描数据和B扫描数据；
将所述离散化A扫描数据分配到像素中；
使用所获取OCT信号来形成至少一个B扫描集群集，所述至少一个B扫描集群集各自包括至少一个B扫描集群，所述至少一个B扫描集群包括至少两个B扫描，所述至少两个B扫描相互平行并且形成平行于所述光束方向的平面；
由每个B扫描集群中的所述至少两个B扫描形成B扫描对；
进行以下操作中的任一个：
在任何B扫描的任何偏移之前计算每个B扫描对的总相位误差，并且然后将每个B扫描对中的一个B扫描偏移至少两次；
其中每次偏移是一个像素，并且是在平行于所述光束的方向上；
在每次偏移之后计算每个B扫描对的总相位误差；
识别每个B扫描对的所有计算出的总相位误差中最小的计算出的总相位误差；
识别每个B扫描对的偏移量，所述偏移量导致所述识别的最小计算出的总相位误差；
计算所述至少一个B扫描集群内的每个B扫描的累积轴向偏移；以及
对准所述至少一个B扫描集群内的所有B扫描。
2.如权利要求1所述的OCT***，其中：
所述OCT***具有如下配置：在形成所述B扫描对之后，计算每个B扫描对的相位差；
使用所述计算出的相位差来计算每个总相位误差；并且
每个计算出的相位差构成相位差数据点。
3.如权利要求2所述的OCT***，其中所述相位差数据点用来计算整体运动相位差。
4.如权利要求3所述的OCT***，其中针对整体运动的影响对所述计算出的整体运动相位差进行校正以便提供已校正的相位差。
5.如权利要求4所述的OCT***，其中将基于强度的阈值方法应用于所述已校正的相位差。
6.如权利要求1所述的OCT***，其中所述OCT***具有如下配置：形成包括A扫描对的所述B扫描对，并且使用所述B扫描对内的所述A扫描对中的全部或一部分来计算每个B扫描对的总相位误差；其中所述A扫描对中的一个A扫描是在所述B扫描对的所述B扫描中的一个内，并且所述A扫描对中的另一个A扫描是在所述B扫描对中的另一个B扫描内。
7.如权利要求1所述的OCT***，其中所述OCT***具有使用所述已对准的B扫描来形成所述物理对象的图像的配置。
8.如权利要求1所述的OCT***，其中所述物理对象是人类组织。
9.如权利要求1所述的OCT***，其中所述OCT***包括：
至少一个光源，其提供所述光束；
至少一个回射反射器；
至少一个光纤耦合器或至少一个自由空间耦合器，其将所述光束导引至所述物理对象并且导引至所述至少一个回射反射器；其中被导引至所述物理对象的所述光束形成至少一个反向散射光束；并且其中被导引至所述至少一个回射反射器的所述光束形成至少一个反射的参考光束；
至少一个扫描光学器件，其在所述物理对象上扫描所述至少一个光束；
至少一个检测器；其：
组合所述至少一个反向散射光束与所述至少一个反射光束来形成光干涉，
检测所述至少一个反向散射光束的量值和时间延迟，以及
形成OCT信号；
其中所述至少一个光纤耦合器或所述至少一个自由空间耦合器将所述至少一个反向散射光束和所述至少一个反射光束导引至所述至少一个检测器；
至少一个处理器，其获得并分析由所述至少一个检测器形成的所述OCT信号，并且形成所述物理对象的至少一个图像；以及
至少一个显示器，其显示所述物理对象的所述至少一个图像。
10.如权利要求9所述的OCT***，其中所述OCT***具有基于所述B扫描之间的强度或相位变化来识别运动区域的配置。
11.如权利要求10所述的OCT***，其中所述OCT***具有使用相位方差OCT(PV-OCT)方法、相位对比OCT(PC-OCT)方法、强度/散斑方差OCT(IV-OCT)方法、多普勒OCT(D-OCT)方法、多普勒功率频移OCT(PDS-OCT)方法、频谱分离幅度去相关分析(SSADA)方法、光学微血管造影(OMAG)方法、相关映射OCT(cmOCT)方法或其组合来识别所述运动区域的配置。
12.如权利要求10所述的OCT***，其中所述OCT***具有使用相位方差OCT(PV-OCT)方法的配置。
13.如权利要求10所述的OCT***，其中所述OCT***具有如下配置：
将每个B扫描对中的一个B扫描偏移至少三次。
14.一种包含指令程序的非暂态、有形计算机可读存储介质，所述指令程序致使运行所述指令程序的计算机***充当光学相干断层扫描(OCT)***，包括进行以下操作：
用具有一定束宽和方向的光束来扫描具有一定表面和深度的物理对象；
从所述扫描获取OCT信号；
从所获取OCT信号生成离散化A扫描数据和B扫描数据；
将所述离散化A扫描数据分配到像素中；
使用所获取OCT信号来形成至少一个B扫描集群集，所述至少一个B扫描集群集各自包括至少一个B扫描集群，所述至少一个B扫描集群包括至少两个B扫描，所述至少两个B扫描相互平行并且形成平行于所述光束方向的平面；
由每个B扫描集群中的所述至少两个B扫描形成B扫描对；
进行以下操作中的任一个：
在任何B扫描的任何偏移之前计算每个B扫描对的总相位误差，并且然后将每个B扫描对中的一个B扫描偏移至少两次；
其中每次偏移是一个像素，并且是在平行于所述光束的方向上；
在每次偏移之后计算每个B扫描对的总相位误差；
识别每个B扫描对的所有计算出的总相位误差中最小的计算出的总相位误差；
识别每个B扫描对的偏移量，所述偏移量导致所述识别的最小计算出的总相位误差；
计算所述至少一个B扫描集群内的每个B扫描的累积轴向偏移；以及
对准所述至少一个B扫描集群内的所有B扫描。
15.如权利要求14所述的存储介质，其中所述指令程序致使运行所述指令程序的所述计算机***进行以下操作：
获得并分析由至少一个检测器形成的所述OCT信号，并且形成所述物理对象的至少一个图像；并且
显示所述物理对象的所述至少一个图像。
16.如权利要求14所述的存储介质，其中所述指令程序致使运行所述指令程序的所述计算机***进行以下操作：
将每个B扫描对中的一个B扫描偏移至少三次。
CN201580006926.4A 2014-02-04 2015-02-04 具有相位敏感的b扫描配准的光学相干断层扫描(oct)*** Active CN105939652B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201461935431P 2014-02-04 2014-02-04
US61/935,431 2014-02-04
PCT/US2015/014410 WO2015120017A1 (en) 2014-02-04 2015-02-04 Optical coherence tomography (oct) system with improved motion contrast
PCT/US2015/014468 WO2015120055A1 (en) 2014-02-04 2015-02-04 Optical coherence tomography (oct) system with phase-sensitive b-scan registration
USPCT/US2015/014410 2015-02-04

Publications (2)

Publication Number Publication Date
CN105939652A CN105939652A (zh) 2016-09-14
CN105939652B true CN105939652B (zh) 2018-07-24

Family

ID=53778404

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580006926.4A Active CN105939652B (zh) 2014-02-04 2015-02-04 具有相位敏感的b扫描配准的光学相干断层扫描(oct)***

Country Status (4)

Country Link
US (1) US9763570B2 (zh)
EP (1) EP3102090B1 (zh)
CN (1) CN105939652B (zh)
WO (1) WO2015120055A1 (zh)

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9700206B2 (en) * 2015-02-05 2017-07-11 Carl Zeiss Meditec, Inc. Acquistion and analysis techniques for improved outcomes in optical coherence tomography angiography
US10231619B2 (en) 2015-12-09 2019-03-19 Oregon Health & Science University Systems and methods to remove shadowgraphic flow projections in OCT angiography
EP3216387B1 (en) 2016-03-10 2018-10-10 Canon Kabushiki Kaisha Method and system for motion artefacts removal in optical coherence tomograpy
US20180012359A1 (en) * 2016-07-06 2018-01-11 Marinko Venci Sarunic Systems and Methods for Automated Image Classification and Segmentation
AU2017382218B2 (en) 2016-12-21 2023-05-11 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11241155B2 (en) * 2017-01-27 2022-02-08 The Regents Of The University Of California Optical coherence tomography device for characterization of atherosclerosis with a 1.7 micron swept laser source
JP2020513983A (ja) 2017-03-23 2020-05-21 ドヘニー アイ インスティテュート 光干渉断層撮影の多重ｅｎ−ｆａｃｅ血管造影平均化システム、方法および装置
JP2019047841A (ja) 2017-09-07 2019-03-28 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
JP2019047839A (ja) * 2017-09-07 2019-03-28 キヤノン株式会社 画像処理装置、位置合わせ方法及びプログラム
JP7027076B2 (ja) * 2017-09-07 2022-03-01 キヤノン株式会社 画像処理装置、位置合わせ方法及びプログラム
WO2019075376A1 (en) 2017-10-13 2019-04-18 The Research Foundation For The State University Of New York OPTICAL DOPPLER TOMOGRAPHY WITH WAVELENGTH DIVISION MULTIPLEXING SCANNED SOURCE
CN108309245B (zh) * 2018-02-23 2020-12-01 东北大学秦皇岛分校 一种谱域光学相干层析的相位补偿修正方法及***
CN112638233B (zh) 2018-06-20 2024-06-14 奥克塞拉有限公司 基于家庭的眼科应用的微型移动低成本光学相干断层扫描***
US10581210B2 (en) 2018-07-30 2020-03-03 Te Connectivity Corporation Receptacle assembly having cabled receptacle connectors
EP3608625B1 (de) 2018-08-07 2023-10-25 Hexagon Technology Center GmbH Oct vermessung
CN110361869B (zh) * 2019-07-15 2021-02-23 北京理工大学 一种基于光学斩波器的oct散斑噪声抑制装置
WO2021134087A1 (en) 2019-12-26 2021-07-01 Acucela Inc. Optical coherence tomography patient alignment system for home based ophthalmic applications
JP7512643B2 (ja) 2020-03-30 2024-07-09 株式会社ニデック 眼科画像処理装置および眼科画像処理プログラム
US10959613B1 (en) 2020-08-04 2021-03-30 Acucela Inc. Scan pattern and signal processing for optical coherence tomography
CN116390683A (zh) 2020-08-14 2023-07-04 奥克塞拉有限公司 用于光学相干断层扫描a扫描去弯曲的***和方法
US11393094B2 (en) 2020-09-11 2022-07-19 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
CN112155511B (zh) * 2020-09-30 2023-06-30 广东唯仁医疗科技有限公司 一种基于深度学习补偿oct采集过程中人眼抖动的方法
AU2021352417A1 (en) 2020-09-30 2023-04-06 Acucela Inc. Myopia prediction, diagnosis, planning, and monitoring device
EP4000500A1 (en) * 2020-11-12 2022-05-25 Optos PLC Volumetric oct image data processing
DK4000502T3 (da) * 2020-11-12 2024-02-05 Optos Plc Volumetrisk oct-billeddatabearbejdning
WO2022165526A1 (en) * 2021-01-28 2022-08-04 Cornell University Resolution-enhanced optical coherent tomography
JP2024511085A (ja) 2021-03-24 2024-03-12 アキュセラ インコーポレイテッド 軸長測定モニタ
CN114909987B (zh) * 2022-05-09 2024-02-02 东南大学 一种b扫描分割法抑制样品整体抖动的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103002794A (zh) * 2010-02-08 2013-03-27 奥勒冈保健科学大学 用于超高灵敏光学微血管显影术之设备及方法
CN103063129A (zh) * 2011-10-24 2013-04-24 佳能株式会社 多波长干涉仪、测量设备以及测量方法
WO2013103080A1 (ja) * 2012-01-05 2013-07-11 国立大学法人筑波大学 波長走査型光干渉断層計及びその位相安定化方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1872084A1 (en) * 2005-03-21 2008-01-02 EPFL Ecole Polytechnique Fédérale de Lausanne Phase sensitive fourier domain optical coherence tomography
US7995814B2 (en) * 2006-06-26 2011-08-09 California Institute Of Technology Dynamic motion contrast and transverse flow estimation using optical coherence tomography
US8244334B2 (en) 2007-04-10 2012-08-14 University Of Southern California Methods and systems for blood flow measurement using doppler optical coherence tomography
JP5448353B2 (ja) * 2007-05-02 2014-03-19 キヤノン株式会社 光干渉断層計を用いた画像形成方法、及び光干渉断層装置
US8625104B2 (en) * 2009-10-23 2014-01-07 Bioptigen, Inc. Systems for comprehensive fourier domain optical coherence tomography (FDOCT) and related methods
US8363225B2 (en) * 2010-01-14 2013-01-29 University Of Rochester Optical coherence tomography (OCT) apparatus, methods, and applications
US8433393B2 (en) 2011-07-07 2013-04-30 Carl Zeiss Meditec, Inc. Inter-frame complex OCT data analysis techniques
US10215551B2 (en) * 2012-07-27 2019-02-26 Praevium Research, Inc. Agile imaging system
CN104769481B (zh) * 2012-10-12 2018-12-18 统雷有限公司 紧凑、低色散以及低像差自适应光学扫描***
US9247874B2 (en) * 2013-02-01 2016-02-02 Carl Zeiss Meditec, Inc. Systems and methods for sub-aperture based aberration measurement and correction in interferometric imaging

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103002794A (zh) * 2010-02-08 2013-03-27 奥勒冈保健科学大学 用于超高灵敏光学微血管显影术之设备及方法
CN103063129A (zh) * 2011-10-24 2013-04-24 佳能株式会社 多波长干涉仪、测量设备以及测量方法
WO2013103080A1 (ja) * 2012-01-05 2013-07-11 国立大学法人筑波大学 波長走査型光干渉断層計及びその位相安定化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
common approach for compensation of axial motion artifacts in swept source OCT and dispersion in Fourier domain OCT;Hillmann et.al.;《Optics Express》;20120803;第20卷(第6期);全文 *
In vivo volumetric imaging of vascular perfusion within human retina and choroids with optical micro-angiography;Lin An et.al.;《OPTICS EXPRESS》;20080721;第16卷(第15期);全文 *
Mobility and transverse flow visualization using phase variance contrast with spectral domain optical coherence tomography;Fingler et.al.;《Optics Express》;20070918;第15卷(第20期);全文 *

Also Published As

Publication number Publication date
US9763570B2 (en) 2017-09-19
US20170020387A1 (en) 2017-01-26
EP3102090B1 (en) 2020-03-25
EP3102090A1 (en) 2016-12-14
CN105939652A (zh) 2016-09-14
EP3102090A4 (en) 2017-09-13
WO2015120055A1 (en) 2015-08-13

Similar Documents

Publication Publication Date Title
CN105939652B (zh) 具有相位敏感的b扫描配准的光学相干断层扫描(oct)***
CN106028912B (zh) 具有改进的运动对比的光学相干断层扫描(oct)***
Kashani et al. Optical coherence tomography angiography: a comprehensive review of current methods and clinical applications
US10398302B2 (en) Enhanced vessel characterization in optical coherence tomograogphy angiography
EP2852317B1 (en) Analysis and visualization of oct angiography data
Schwartz et al. Phase-variance optical coherence tomography: a technique for noninvasive angiography
DeBuc A review of algorithms for segmentation of retinal image data using optical coherence tomography
US20160317026A1 (en) Optical coherence tomography system for health characterization of an eye
JP6843125B2 (ja) 高感度フロー視覚化方法
US20160287071A1 (en) Methods of measuring total retinal blood flow using en face doppler oct
US20160331227A1 (en) Extended duration optical coherence tomography (oct) system
JP2018038689A (ja) 眼科撮影装置及び眼科画像処理装置
JP7111874B2 (ja) 眼科撮影装置
Smid Optical coherence tomography in chorioretinal vascular diseases
Lindberg Enhanced Depth Imaging OCT and OCT Angiography for Quantifying Optic Disc Drusen and Related Complications
Li Computational Methods for Enhancements of Optical Coherence Tomography
Huang et al. In vivo microvascular network imaging of the human retina
Chu Ophthalmic Application of Optical Coherence Tomography Based Angiography for Choroidal Imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant