CN110473265A - Oct图像处理 - Google Patents

Oct图像处理 Download PDF

Info

Publication number
CN110473265A
CN110473265A CN201910383653.5A CN201910383653A CN110473265A CN 110473265 A CN110473265 A CN 110473265A CN 201910383653 A CN201910383653 A CN 201910383653A CN 110473265 A CN110473265 A CN 110473265A
Authority
CN
China
Prior art keywords
coordinate
pixel
oct
coordinate system
image data
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.)
Granted
Application number
CN201910383653.5A
Other languages
English (en)
Other versions
CN110473265B (zh
Inventor
艾伦·邓肯·弗莱明
G·穆约
迈克尔·维尔霍克
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.)
Optos PLC
Original Assignee
Optos PLC
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 Optos PLC filed Critical Optos PLC
Publication of CN110473265A publication Critical patent/CN110473265A/zh
Application granted granted Critical
Publication of CN110473265B publication Critical patent/CN110473265B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/10—Geometric effects
    • G06T15/20—Perspective computation
    • G06T15/205—Image-based rendering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/80—Geometric correction
    • 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
    • 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/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
    • 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
    • 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
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/04—Context-preserving transformations, e.g. by using an importance map
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/18—Image warping, e.g. rearranging pixels individually
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • 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/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Ophthalmology & Optometry (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Quality & Reliability (AREA)
  • Computer Graphics (AREA)
  • Computing Systems (AREA)
  • Pathology (AREA)
  • Epidemiology (AREA)
  • Signal Processing (AREA)
  • Primary Health Care (AREA)
  • General Business, Economics & Management (AREA)
  • Business, Economics & Management (AREA)
  • Eye Examination Apparatus (AREA)
  • Image Processing (AREA)

Abstract

公开了OCT图像处理。提供了渲染由在一定范围的扫描位置扫描眼睛的视网膜的OCT扫描仪获取的光学相干断层扫描OCT视网膜图像数据的装置(100‑1)，其中每个扫描位置与在第一坐标系中的相应坐标关联，且渲染的OCT视网膜图像数据的各像素与在不同于第一坐标系的第二坐标系中的相应坐标关联。该装置包括：布置为接收由OCT扫描仪获取的OCT视网膜图像数据的通信模块(110)；及布置为用从在第一坐标系中的坐标到第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值的坐标确定模块(120‑1)。该装置还包括布置为在第二坐标系中的坐标的确定的值处的像素的值之间插值以计算渲染的OCT视网膜图像数据的像素的值的插值模块(130‑1)。

Description

OCT图像处理
技术领域
本发明总体上涉及图像处理的领域，且更具体地，涉及对视网膜的光学相干断层扫描图像的渲染以便减少几何失真。
背景
对视网膜的光学相干断层扫描(OCT)能够生成视网膜的(在成像平面内的)二维断层扫描或三维图像。这种图像在临床上广泛用于观察视网膜的内部结构、视神经乳头和眼睛的其他结构。OCT扫描仪可用于进行视网膜层厚度的测量并观察异常，例如水肿、视网膜脱离、黄斑裂孔和肿瘤。
OCT扫描仪照惯例通过跨越一定范围的偏转角扫描激光束来生成在眼睛的区域上的图像。例如，可以有被配置成在名义上水平和垂直的方向上偏转光束的两个光束偏转设备。
当激光束穿透视网膜时，在激光束的瞬时位置处获得的数据被称为A扫描。照惯例，来自这种设备的图像是通过对准A扫描数据以产生二维或三维图像而生成的。照惯例，在每个A扫描的定时和所捕获的数据被渲染所在的像素坐标之间存在线性的关系。
概述
上面概述的传统图像渲染将不生成实际视网膜的在几何上真实(即，角度保持、等距缩放)的图像，因为(1)激光束在视网膜上的入射角在扫描过程期间变化、(2)在成像设备的光学器件中存在非线性度、以及(3)在眼睛的光学特性中存在空间和角度变化。以在所渲染的位置和捕获的时间之间的线性关系渲染的图像将包含一些失真。当在眼睛的后极处对狭窄的视网膜视场成像时，这种失真可能小到足以被忽略。然而，当对较宽的视网膜视场或远离后极移动的视场成像时，由这种稚拙派的图像渲染产生的几何失真可能是显著的。
本发明的发明人设计了一种手段来产生视网膜的OCT图像的几何上精确的渲染，而不管视场的宽度或者扫描的位置是否是中心的或在视网膜的周边中。
更具体地，发明人设计了一种渲染OCT视网膜图像数据的方法，该OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个扫描位置与在第一坐标系(“捕获坐标系”)中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于第一坐标系的第二坐标系(“渲染坐标系”)中的相应坐标相关联。该方法包括接收由OCT扫描仪获取的OCT视网膜图像数据，使用从在第一坐标系中的坐标到在第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值，并且在第二坐标系中的确定的坐标的值处的像素的值之间插值以计算所渲染的OCT视网膜图像数据的像素的值。
发明人还设计了一种渲染OCT视网膜图像数据的方法，该OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于第一坐标系的第二坐标系中的相应坐标相关联。该方法包括接收由OCT扫描仪获取的OCT视网膜图像数据，并通过下列操作来使用接收到的OCT视网膜图像数据处理像素阵列以生成所渲染的OCT视网膜图像数据：使用从在第一坐标系中的坐标到在第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值；在第二坐标系中确定的坐标的值之间插值以计算像素阵列中的像素在第一坐标系中的坐标的值；以及使用所计算的在像素阵列中的像素的在第一坐标系中的坐标的值来在接收到的OCT视网膜图像数据的像素的值之间插值，以计算在像素阵列中的像素的值。
发明人还设计了一种渲染OCT视网膜图像数据的方法，该OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于第一坐标系的第二坐标系中的相应坐标相关联。该方法包括接收(S10)由OCT扫描仪获取的OCT视网膜图像数据，并通过下列操作来使用接收到的OCT视网膜图像数据处理像素阵列以生成所渲染的OCT视网膜图像数据：对于像素阵列的每个像素，使用从在第二坐标系中的坐标到在第一坐标系中的坐标的变换来确定在第一坐标系中的相应坐标；以及通过使用在第一坐标系中的所确定的坐标以在接收到的OCT视网膜图像数据的像素的值之间插值来计算像素阵列的像素的值。
发明人还设计了一种非暂时性计算机可读存储介质或携带计算机程序指令的信号，当计算机程序指令由计算机执行时使计算机执行一种或更多种上面阐述的方法。
发明人还设计了一种用于渲染OCT视网膜图像数据的装置，该OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于第一坐标系的第二坐标系中的相应坐标相关联。该装置包括：通信模块，其被布置为接收由OCT扫描仪获取的OCT视网膜图像数据；以及坐标确定模块，其被布置为使用从在第一坐标系中的坐标到在第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值。该装置还包括插值模块，该插值模块被布置为在第二坐标系中确定的坐标的值处的像素的值之间插值以计算所渲染的OCT视网膜图像数据的像素的值。
发明人还设计了一种用于渲染OCT视网膜图像数据的装置，该OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于第一坐标系的第二坐标系中的相应坐标相关联。该装置包括：通信模块，其被布置为接收由OCT扫描仪获取的OCT视网膜图像数据；以及坐标确定模块，其被布置为使用从在第一坐标系中的坐标到在第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值。该装置还包括插值模块，该插值模块被布置为通过下列操作来使用接收到的OCT视网膜图像数据处理像素阵列以生成所渲染的OCT视网膜图像数据：在第二坐标系中的确定的坐标的值之间插值以计算像素阵列中的像素在第一坐标系中的坐标的值；以及使用所计算的像素阵列中的像素的在第一坐标系中的坐标的值来在接收到的OCT视网膜图像数据的像素的值之间插值，以计算在像素阵列中的像素的值。
发明人还设计了一种用于渲染OCT视网膜图像数据的装置，该OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于第一坐标系的第二坐标系中的相应坐标相关联。该装置包括：通信模块，其被布置为接收由OCT扫描仪获取的OCT视网膜图像数据；以及坐标确定模块，其被布置为：对于被处理以生成所渲染的OCT视网膜图像数据的像素阵列的每个像素，使用从在第二坐标系中的坐标到在第一坐标系中的坐标的变换来确定在第一坐标系中的相应坐标。该装置还包括插值模块，该插值模块被布置为通过使用在第一坐标系中的所确定的坐标以在像素阵列的像素的值之间插值来计算像素阵列的像素的值。
附图说明
现在通过非限制性示例参考附图来详细解释本发明的实施例，其中：
图1是根据本发明的第一实施例的用于渲染OCT视网膜图像数据的装置的示意图；
图2是示出根据实施例的用于渲染OCT视网膜图像数据的装置的硬件配置的示例的框图；
图3是示出根据第一实施例的渲染OCT视网膜图像数据的方法的流程图；
图4是示出图3中的过程S20的细节的流程图；
图5是根据第一实施例的变形的用于渲染OCT视网膜图像数据的装置的示意图；
图6是示出根据变形的渲染OCT视网膜图像数据的方法的流程图；
图7是根据本发明的第二实施例的用于渲染OCT视网膜图像数据的装置的示意图；
图8是示出根据第二实施例的渲染OCT视网膜图像数据的方法的流程图；
图9是根据本发明的第三实施例的用于渲染OCT视网膜图像数据的装置的示意图；以及
图10是示出根据第三实施例的渲染OCT视网膜图像数据的方法的流程图。
具体实施方式
现在将参考附图更详细地描述本发明的实施例。
下面的符号用在下面提供的实施例的描述中：
本文描述的OCT扫描仪的配置是通过扫描参数H、V、和D定义的，其中，H和V是水平和垂直扫描位置或控制信号，而D是光路距离。OCT扫描仪的配置利用这些参数被表示为(H,V,D)。
由本文描述的用于渲染OCT视网膜图像数据的装置接收的所捕获的图像数据被表示为Icap，并且在所捕获的图像的像素坐标系中的像素坐标被表示为(x,y,z)。
所渲染的图像数据被表示为Iren，以及在所渲染的图像的像素坐标系中的像素坐标被表示为(X,Y,Z)。
在一些实施例中，在系统建模中使用的真实世界坐标系被表示为(A,B,C)，并且用于图像渲染的真实世界坐标系被表示为(A’,B’,C’)，以笛卡尔或球面坐标系为示例。
对于上述坐标系中的每一个，下标“cap”和“ren”在下面提供的实施例的描述中用于表示所捕获的(cap)或所渲染的(ren)像素的在相应坐标系中的坐标的值。例如，使用在所捕获的图像中的像素的坐标系(A’,B’,C’)表达的坐标值被表示为(A’cap,B’cap,C’cap)。作为另一个示例，在所渲染的图像的像素坐标系中的像素坐标(x,y,z)的值被表示为(xren,yren,zren)。
实施例1
图1是根据本发明的第一实施例的用于渲染OCT视网膜图像数据的装置100-1的示意图。由装置100-1处理的OCT视网膜图像数据可以由本领域中的技术人员已知的任何合适的OCT扫描仪获得，OCT扫描仪在其执行的OCT扫描的过程期间获取OCT视网膜图像数据，其中该扫描覆盖在患者的眼睛中的一定范围的扫描位置。
例如，OCT扫描仪可以以组合的扫描激光检眼镜(SLO)和OCT成像系统的形式被提供，如在国际专利申请号PCT/BG2013/052556(被公布为WO2014/053824 A1)中所述，该专利申请的内容通过引用被全部并入本文。所描述的宽视场SLO-OCT成像系统包括两个扫描元件，该两个扫描元件可以围绕相应的旋转轴独立地旋转，以沿着不同的扫描轴偏转来自光源(例如超辐射发光二极管(SLD)或诸如此类)的准直光的采样光束，以执行对视网膜的表面及其内部的“切片”的二维扫描。扫描轴通常是正交的，并且包括水平轴(H轴)和垂直轴(V轴)。对于从光源到视网膜的光路距离D的多个值中的每一个，OCT扫描仪用参考光束干涉样本光束从视网膜的反射，并且对于由扫描覆盖的多个二维扫描位置中的每一个二维扫描位置测量被干涉光束的强度。
这些扫描位置中的每一个扫描位置都与指示相应扫描位置的相对应OCT扫描参数值相关联。更具体地，在沿H轴或V轴执行一维扫描的情况下，例如，沿H轴(或V轴，视情况而定)的每个扫描位置与指示该扫描位置的扫描参数的值H(或V轴，视情况而定)(例如沿H轴(或V轴，视情况而定)扫描采样光束的扫描元件的相对应的偏转角、用于使扫描元件实现这种偏转角的控制信号、或类似的操作参数)相关联。例如，可以根据校准或光学系统建模来测量或确定偏转角。引申开来，在沿着H轴和V轴执行二维扫描的情况下，每个扫描位置与指示该扫描位置的相应的一对扫描参数值((再次，如根据校准或光学系统建模测量或确定的)例如两个扫描元件的相对应的偏转角、用于使扫描元件实现这种偏转角的控制信号、或诸如此类)相关联。
OCT扫描仪因此在视网膜的深度方向上(即，沿着光束的方向，在D的不同值处)对于由H和V的相应值表示的每个二维扫描位置获得多个OCT测量，(对于H的任何给定值和V的任何给定值获取的)OCT测量的每个这样的组被广泛地称为“A扫描”，如上面所提到的。通过在沿着二维扫描轴之一(例如，沿着指定的H轴的方向)的点处获取多个A扫描，OCT扫描仪可以获取被广泛称为“B扫描”的测量，并且通过在沿着剩余的二维扫描轴(在这个示例中，指定的V轴)的点处获取多个B扫描，OCT扫描仪可以获取被称为“C扫描”的测量。
下面描述的图像处理操作利用在由OCT扫描仪获取的每个所测量的强度值(即像素值)和指示测量在扫描期间被执行所在的位置的OCT扫描参数值的相对应集合之间的已知关系。因此，在OCT视网膜图像数据中的每个像素值都与在OCT扫描仪的机器坐标系中的相应坐标相关联，该坐标由D的值与H的值和/或V的值相结合来定义，当像素值被测量时，OCT扫描仪将根据这些值来进行操作。如下提供这个关系可如何被确定的示例。
在下文中描述的视网膜图像处理操作涉及不仅以在C扫描的采集期间获得的四维数据集(其中每个数据元素包括测量值和在三维空间中的相关联的点)的形式提供的OCT视网膜图像数据的处理而且还有与B扫描相关联的三维数据集(其中每个数据元素包括测量值和在二维空间中的相关联的点)的处理。B扫描的三维数据集可以被认为代表定义通过视网膜的切片的二维图像的像素的二维阵列，而C扫描数据的四维数据集可以被认为代表定义视网膜的三维图像的体素的三维阵列。在下文中描述的图像处理方法可以将这个原始二维/三维图像变换成在几何上真实的渲染，根据该渲染可以进行更准确的测量。
装置100-1包括可操作来接收由OCT扫描仪获取的OCT视网膜图像数据的通信模块110。如在本实施例中的，OCT视网膜图像数据可以在由OCT扫描仪获取之后被存储在任何合适的存储介质(例如CD或DVD或计算机硬盘)上，并且随后在计算机从存储介质读取OCT视网膜图像数据时或之后(例如，直接地或通过网络例如局域网(LAN)或因特网)从计算机或诸如此类进行接收。可选地，OCT视网膜图像数据可以在它在视网膜扫描的执行期间被获取时直接从OCT扫描仪接收。
如在下文中将进一步解释的，每个扫描位置与由在第一坐标系(在这里也被称为“捕获坐标系”)中的坐标的值定义的相应点相关联，如在本实施例中，第一坐标系可以是由(A,B,C)表示的第一真实世界坐标系，并且在对光在OCT扫描仪中的传播建模时被使用。此外，如在本实施例中，所渲染的OCT视网膜图像数据的每个像素可以与由在第二坐标系(在这里也被称为“渲染坐标系”)中的坐标的值定义的相应点相关联，第二坐标系是由(A’,B’,C’)表示的另一个真实世界坐标系。真实世界坐标系的示例是笛卡尔坐标系和球面坐标系，笛卡尔坐标系的坐标照惯例被表示为(x,y,z)，以及球面坐标系的坐标照惯例被表示为
装置100-1还包括坐标确定模块120-1，其被布置成使用在第一坐标系中的坐标和在第二坐标系中的坐标之间的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值。如在本实施例中，坐标确定模块120-1可以被布置成不仅使用所述变换而且使用(i)在OCT视网膜图像数据的采集期间使用的至少一组OCT扫描参数值(每一所述一组OCT扫描参数值指示在扫描内的相应位置)和(ii)在对应于至少一组OCT扫描参数值的第一坐标系中的坐标之间的映射来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值。下面描述坐标确定模块120-1可如何使用这种变换和映射来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值的示例。
装置100-1还包括插值模块130-1，该插值模块130-1被布置为在由坐标确定模块120-1确定的第二坐标系中的坐标(其通常在第二坐标系中不规则地间隔开)的值处的像素的值之间插值，以计算所渲染的OCT视网膜图像数据的像素(其通常在第二坐标系中有规则地间隔开，被分布以定义点的网格)的值。换句话说，插值模块130-1被布置成使用接收到的OCT视网膜图像数据来处理像素阵列，以通以下操作来生成所渲染的OCT视网膜图像数据：在由坐标确定模块120-1确定的第二坐标系中的坐标处的像素的值之间插值以计算像素阵列的像素的值，从而生成所渲染的OCT视网膜图像数据。下面还提供了该插值过程的另外的细节。
图2示出了装置100-1可如何在可编程信号处理硬件中被实现的示例。图2中所示的信号处理装置200包括通信接口(I/F)210，其用于接收由OCT扫描仪获取的OCT视网膜图像数据并输出所渲染的OCT视网膜图像数据，使得它可以由信号处理装置或外部设备显示或以另外方式处理。信号处理装置200还包括处理器(CPU)220、工作存储器230(例如随机存取存储器)和存储计算机可读指令的指令储存装置240，计算机可读指令当由处理器220执行时使处理器220执行下文描述的处理操作，以渲染接收到的OCT视网膜图像数据(即，处理接收到的OCT数据以生成表示已经被扫描以生成OCT数据的视网膜的部分的二维或三维图像的图像数据)。指令储存装置240可以包括预加载有计算机可读指令的ROM(例如，以电可擦除可编程只读存储器(EEPROM)或闪存的形式)。可选地，指令储存装置240可以包括RAM或类似类型的存储器，并且计算机可读指令可以从计算机程序产品(例如计算机可读存储介质250(例如CD-ROM等)或者携带计算机可读指令的计算机可读信号260向其进行输入。然而应当注意，本文描述的实施例的装置可以可选地在不可编程硬件(例如专用集成电路(ASIC))中被实现。
在本实施例中，通信I/F 210被配置为执行图1所示的通信模块110的功能，而图2中的包括处理器220、工作存储器230和指令储存装置240的所示的硬件部件的组合270被配置为执行坐标确定模块120-1和插值模块130-1的功能，其现在将参考图3和图4来进行详细描述。
图3是示出一个过程的流程图，装置100-1通过该过程来渲染OCT视网膜图像数据以生成视网膜的被扫描部分的在几何上真实(换句话说，角度保持、等距缩放)的图像。
在过程S10中，通信模块110接收由OCT扫描仪获取的OCT视网膜图像数据。在本实施例中，接收到的数据包括指示通过用来自OCT扫描仪的参考臂的光干涉从视网膜反射的采样光而产生的干涉信号的强度的像素值。
通过示例的方式，在本实施例中，由通信模块110接收由OCT扫描仪捕获的表示二维(B扫描)OCT图像Icap的OCT视网膜图像数据。然而，如上面所提到的，可选地可以由通信模块110接收表示由OCT扫描仪获取的三维(C扫描)OCT图像的OCT视网膜图像数据。接收到的OCT视网膜图像数据的每一列对应于A扫描，每一列具有H和V的已知和恒定值。H和V的这些值连同所捕获的图像Icap的相关联的像素值一起可以由通信模块110接收，使得装置100-1可以将接收到的像素值与H和V的相应的接收到的值相关联地进行存储。可选地，通信模块110可以仅接收像素值本身，下面描述的图像渲染算法基于在接收到的图像数据的像素坐标系中的像素的像素位置(这些像素位置由像素坐标系中的坐标x、y和z表示)以及与获取OCT视网膜图像数据的OCT扫描仪相关的预先存储的校准数据来确定每个像素的H和V值。在这种情况下，校准数据使在图像数据中的每个像素位置与H和V的相应值相关联。坐标zc表示沿着对应于激光束方向的图像轴的坐标，而x和y表示沿着与该图像轴正交的轴的坐标。
在过程S20中，坐标确定模块120-1使用从在第一坐标系中的坐标到在第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在第二坐标系中的坐标的值。如在本实施例中的，坐标确定模块120-1还可以在该过程中使用在OCT视网膜图像数据的获取期间使用的至少一组OCT扫描参数值和在对应于该至少一组OCT扫描参数值的在第一坐标系(A,B,C)中的坐标之间的映射。
可以通过使用光学系统建模连同OCT扫描仪和假设的眼睛模型的组合模型来导出这个映射，以针对为激光束扫描元件设置的偏转角的范围来找到在OCT扫描仪内的光线到视网膜的路径。可以以多种方式(例如固定几何形状和尺寸的特定的理想化眼睛或者针对特定的眼睛从眼睛生物统计数据导出的参数化模型)来定义假设的眼睛模型。光学系统建模的终点是在OCT扫描参数H、V和D的值与在第一(真实世界)坐标系(A,B,C)中的坐标的相对应的值之间的映射。光学系统建模可以基于OCT扫描仪的光学结构的知识和/或基于人造目标的结构的知识以及当该人造目标的结构被OCT扫描仪成像时生成的图像。
对于参数H和V的可能值(即，对于可能的光线)，可以使用光学系统建模来确定下面的数据：
1.真实世界坐标(其中光线射到视网膜)的值。
2.光线在眼睛的玻璃体内的方向，作为在真实世界坐标中的矢量。
3.从光发射器到光线与视网膜的交叉点的距离，其被称为光路距离(D)。
然而，应该注意，并非所有上述数据都需要在光学系统建模中被获得，并且光学系统建模的结果可以以许多不同的形式进行表示。
关于在OCT扫描参数D的值和在像素坐标系中的相对应的坐标z之间的映射，这可以如下被确定。
D的值对于具有z的给定值的所有像素是恒定的，并且在zcap方向上的像素之间的间距以真实距离单位计是已知的。因此，D＝(z–z1)*p，其中，p是等于在z方向上的像素间间距的已知值，以及z1是需要被确定的偏移。可以用多种不同的方式之一来确定z1的值。
例如，可以分析(例如，通过计算机图像分析)在接收到的OCT视网膜图像数据中的像素值以检测视网膜的层，并且可以找到在视网膜的该层处的图像列的子集的z的值。在沿着这个检测到的层的点处，假设从光学系统建模的输出中已知D的值。这允许确定针对z1的估计。
作为替代方案，可以通过使用来自OCT扫描仪的控制系统的信息来确定z1。在参考臂中的光路长度由该设备精确地控制，且因此在位置控制参数和z1之间存在已知的关系。这个关系可以从OCT扫描仪的设计的知识中被确定或从校准中被获得。
然而，应当注意，z1不需要通过上述(或其他)方法中的任何一种来明确地被确定；优选获得的是针对zcap的每个值的D值的确定，以及因此在OCT扫描参数H、V和D的值和相对应坐标(x,y,z)之间的映射的确定。
为了将OCT视网膜图像数据像素的坐标从第一坐标系重新投影到第二坐标系中，坐标确定模块120-1使用在OCT扫描参数H、V和D的值与对应的坐标(x,y,z)之间的所确定的第一映射以及在OCT扫描参数H、V和D的值与在第一真实世界坐标系(A,B,C)中的对应的坐标之间的所确定的第二映射，现在将参考图4进行描述。
图4是示出坐标确定模块120-1可如何确定在接收到的OCT视网膜图像数据的像素的在第二坐标系中的坐标的值的示例的流程图。
在过程S21中，选择接收到的OCT视网膜图像数据的第一像素。然后，在过程S22中，坐标确定模块120-1使用上面所述的第一映射来为在当前迭代中考虑的第一像素确定OCT扫描参数H、V和D的对应值。在过程S23中，使用上面所述的在OCT扫描参数的值和在第一真实世界坐标系(A,B,C)中的对应坐标之间的(第二)映射来将这些扫描参数值转换成在第一真实世界坐标系(A,B,C)中的坐标的对应值。因此，在图像渲染过程的这个阶段确定在接收到的OCT视网膜图像数据中的像素的真实世界坐标的值。
根据在光学系统建模的输出中可用的数据，可以简化像素的真实世界坐标的计算。例如，在眼睛的玻璃体内的光线的方向被称为在真实世界坐标系中的矢量的实施例的变形中，可以对在每个图像列中的仅一个像素执行上述迭代(A扫描)。然后，可以通过在迭代列中的像素上迭代并加上通过像素间距离(其是已知的)按比例缩放的光线的方向矢量来为列中的剩余像素计算在第一真实世界坐标系(Acap,Bcap,Ccap)中的像素的坐标的值。
选择第二真实世界坐标系(A’,B’,C’)，其中渲染OCT视网膜图像数据。这个选择可以以各种方式被通知或完成。例如，可以选择在多个图像捕获事件中被标准化的或者与在另一个捕获事件中使用的坐标系相同的笛卡尔或球面坐标系。这可能适合于例如(a)便于由用户将当前图像与另一图像进行比较，或者(b)便于由自动系统使用当前图像和另一图像进行疾病变化检测或图像配准的比较。可选地，可以使用以某种最佳方式与所捕获的数据对准的笛卡尔或球面坐标系。第二坐标系(A’,B’,C’)的一个或更多个轴可以与所捕获的数据的一个或更多个轴紧密对准。例如，如果使用笛卡尔坐标系，一个轴可以以某种最佳方式与扫描方向H对准，而另一个轴可以以某种最佳方式与扫描方向V对准。可选地，如果使用球形坐标系，则径向轴可以与光传播方向对准，D是沿着该光传播方向被测量的。
无论什么方法用于决定第二真实世界坐标系(A’,B’,C’)，在第一真实世界坐标系(A,B,C)和第二真实世界坐标系(A’,B’,C’)之间都必须有已知的关系，使得存在用于将坐标值从一个坐标系变换到另一个坐标系的变换例如，该步骤可能需要均匀化矩阵P(homogenized matrix，P)的计算，以便可以通过矩阵乘法来执行转换：
(A’,B’,C’,1)T＝P(A,B,C,1)T
对于图像渲染，选择在将被处理以产生所渲染的图像的像素阵列中的像素的像素坐标值(Xren,Yren,Zren)的范围以及在第二真实世界坐标系(A’,B’,C’)和所渲染的像素网格坐标(X,Y,Z)之间的线性映射。每组像素坐标可以简单地是标识在像素阵列/网格中的像素的位置的索引的一组(整数)值或者通过任何适当的一对一映射从其导出的一组值。选择范围的一种方法是，在所捕获的图像数据的每个轴上取在第二真实世界坐标系(A’,B’,C’)中的它的坐标的最小值和最大值min(A’cap)、max(A’cap)、min(B’cap)、max(B’cap)、min(C’cap)、max(C’cap)。如果上述线性映射被选择成使得在对应于相邻的所捕获的像素的Xcap、Ycap和Zcap之间的间距的平均值在每个维上近似为1，则所捕获的信息将在所渲染的图像中被很好地渲染。
在过程S24中，使用上面提到的变换来将在第一真实世界坐标系(A,B,C)中的像素的所确定的坐标值(Acap,Bcap,Ccap)转换成在第二真实世界坐标系(A’,B’,C’)中的相对应的坐标值(A’cap,B’cap,C’cap)。
在过程S25中，坐标确定模块120-1确定在接收到的视网膜图像数据中的所有像素(或指定的感兴趣子集)是否已经被处理。如果否，则在过程S26中选择下一个像素，并且处理循环回到过程S22。另一方面，当在过程S25中确定不再有像素要处理时，S20的处理完成，并且处理继续进行到图3中的S30。以这种方式，接收到的OCT视网膜图像数据的像素被迭代地处理，以确定它们在第一真实世界坐标系中的相应坐标，并且然后将这些坐标重新投影到第二真实世界坐标系中。在以这种方式确定的第二真实世界坐标系((A’,B’,C’)中的坐标(A’cap,B’cap,C’cap)的值通常在该坐标系中不规则地间隔开。
再次参考图3，在过程S30中，插值模块130-1在第二坐标系(A’,B’,C’)中的坐标(A’cap,B’cap,C’cap)的所确定的值处的像素的值之间插值，以计算将形成所渲染的OCT视网膜图像数据的像素(在坐标值(A’ren,B’ren,C’ren)处)的值。
因此，对于将被处理以生成所渲染的OCT图像数据的像素阵列的像素(其为了方便起见在本文被简单地称为所渲染的OCT视网膜图像数据的像素)，通过在被投影到第二真实世界坐标系中的接收到的OCT视网膜图像数据Icap的像素的不规则地间隔开的坐标(A’cap,B’cap,C’cap)处的已知值之间进行插值来确定在第二真实坐标系中的有规则地间隔开的坐标(A’ren,B’ren,C’ren)处的所渲染的OCT视网膜图像数据的像素的强度值Iren。要做到这一点，插值模块130-1执行散点插值(scatter interpolation)；与规则插值(其中值在有规则地间隔开的位置处是已知的以及找到在介于中间的位置处的值是简单的)不同，对于散点插值，值仅在不规则地间隔开的位置处是已知的，且从这些找到在有规则地间隔开的位置处的值更具挑战性。然而，对这样的数据执行插值的方法是众所周知的，并包括例如基于Delauney三角测量和薄板样条的三角网。
在上面提到的实施例的变形中，可通过下列操作来在过程S20中确定在A扫描的接收到的OCT视网膜图像数据中的像素的在第二真实坐标系中的坐标：使用变换和第一映射以确定在A扫描的接收到的OCT视网膜图像数据中的第一像素的在第二真实坐标系中的坐标，该第一像素与第一组OCT扫描参数值相关联；以及使用变换 和在第一真实世界坐标系中的指示OCT扫描仪根据第一组OCT扫描参数值操作时入射在眼睛的视网膜上或在眼睛的视网膜中传播的光束的方向的矢量以确定在A扫描的接收到的OCT视网膜图像数据中的剩余像素的在第二真实世界坐标系中的坐标。
然后，由通信模块110输出所渲染的OCT视网膜图像数据，使得它可以显示在任何合适的显示设备(例如计算机屏幕或诸如此类)上和/或自动地或在用户指导下被处理以确定在所渲染的OCT视网膜图像数据中的一个或更多个感兴趣特征的尺寸测量。
现在将参考图5和图6描述根据第一实施例的变形的用于渲染OCT视网膜图像数据的装置100-2。
在图5中示意性示出的变形的装置100-2与第一实施例的装置100-1的不同之处仅在于插值模块130-2的功能。装置100-2在信号处理硬件中的实现及其剩余部件的功能与在第一实施例中的相同，且因此不再次被描述。现在将参考图6描述插值模块130-3的功能方面的差异。
图6是示出一个过程的流程图，装置100-2可通过该过程来渲染OCT视网膜图像数据以生成视网膜的被扫描部分的在几何上真实的图像。过程S10和S20与第一实施例的过程相同。然而，该过程(该变形的插值模块130-2通过该过程来使用接收到的OCT视网膜图像数据处理像素阵列以生成所渲染的OCT视网膜图像数据)不同于由第一实施例的插值模块130-1执行的处理，如将在下面解释的。
在过程S40中，插值模块130-2在接收到的OCT视网膜图像数据的像素的在第二真实世界坐标系(A’,B’,C)中的坐标之间插值，以计算在第一真实世界坐标系(A,B,C)中的像素阵列的像素的坐标。在第二真实世界坐标系(A’,B’,C’)中的接收到的OCT视网膜图像数据的像素的坐标通常是不规则地间隔开的。该变形的插值模块130-2例如使用例如基于Delauney三角测量或薄板样条的三角网来对该数据执行射点插值，以确定在第一真实坐标系(A,B,C)中的像素阵列的像素的坐标。
在过程S50中，插值模块130-2使用在第一真实世界坐标系(A,B,C)中的像素阵列的像素的计算出的坐标来在接收到的OCT视网膜图像数据的像素的值之间插值，以计算像素阵列的像素的值，从而生成所渲染的OCT视网膜图像数据。
实施例2
现在将参考图7和图8描述根据本发明的第二实施例的用于渲染OCT视网膜图像数据的装置100-3。
在图7中示意性示出的本实施例的装置100-3与第一实施例的装置100-1的不同之处仅在于坐标确定模块120-3的功能。装置100-3在信号处理硬件中的实现、其剩余部件的功能以及可以对装置做出的可能变化与在第一实施例中的相同，且因此不再进行描述。
现在将参考图8描述该过程，本实施例的坐标确定模块120-3通过该过程来在作为图3所示的过程S20的替代方案的过程S20’中确定在接收到的OCT视网膜图像数据中的像素的在第二坐标系中的坐标。
在过程S21’中，坐标确定模块120-3选择在接收到的OCT视网膜图像数据中基本上均匀地分布(在像素位置当中)的在接收到的OCT视网膜图像数据中的像素的子集。
在过程S22’中，坐标确定模块120-3使用在接收到的OCT视网膜图像数据的在(第一)坐标系中的坐标(x,y,z)和将被处理以给出所渲染的图像的像素阵列的在(第二)坐标系中的坐标(X,Y,Z)之间的变换来为像素的子集的每个像素确定在第二坐标系中的相应坐标。因此，使用像素的子集来获得“理想”映射。
如在本实施例中的，坐标确定模块120-3可以使用下列项来确定像素的子集的每个像素的在第二坐标系中的相应坐标：第一映射，其用于将接收到的OCT图像数据的像素的坐标(x,y,z)从第一坐标系转换到第一真实世界坐标系(A,B,C)；第二映射，其用于将在第一真实世界坐标系(A,B,C)中的像素的坐标转换成在第二真实世界坐标系(A’,B’,C’)中的坐标；以及第三映射，其用于将在像素阵列中的像素的坐标(X,Y,Z)转换成在第二真实世界坐标系(A’,B’,C’)中的对应的坐标。
在过程S23’中，坐标确定模块120-3使用在第二坐标系中的像素的子集的像素的所确定的坐标来确定在第一坐标系中的坐标(x,y,z)和在第二坐标系中的坐标(X,Y,Z)之间的参数化映射的参数。关于该映射的参数可以通过精确的拟合例如通过拟合样条(分段多项式)来生成。在这种情况下，在应用理想映射后的点的子集与在应用参数化映射后的点的子集基本上是相同的。可选地，回归技术(regression technique)可以用于基于在应用理想映射后的点的子集和在应用参数化映射后的点的子集之间的差异来最小化损失函数(loss function)。
最后，在过程S24’中，坐标确定模块120-3使用具有所确定的参数的参数化映射来确定在接收到的OCT视网膜图像数据中的像素的在第二坐标系中的坐标(X,Y,Z)。
实施例3
现在将参考图9和图10描述根据本发明的第三实施例的用于渲染OCT视网膜图像数据的装置100-4。
在图9中示意性示出的本实施例的装置100-4与第一实施例的装置100-1的不同之处仅在于坐标确定模块120-4和插值模块130-4的功能。装置100-4在信号处理硬件中的实现和通信模块110的功能与在第一实施例中的相同，且因此不再进行描述。
图10是示出一个过程的流程图，装置100-4通过该过程来渲染OCT视网膜图像数据以生成视网膜的被扫描部分的在几何上真实的图像。过程S10与在第一实施例中的相同。然而，由坐标确定模块120-4和插值模块130-4执行的操作不同于由第一实施例的坐标确定模块120-1和插值模块130-1执行的操作，且这些现在将被描述。
在过程S60中，坐标确定模块120-4使用从作为渲染OCT视网膜图像数据的坐标的在第二坐标系中的坐标到在第一坐标系中的坐标的变换来为将被处理以生成所渲染的OCT视网膜图像数据的渲染像素阵列的每个像素确定作为接收到的OCT视网膜图像数据的坐标的在第一坐标系中的相应坐标(x,y,z)。因此，在第二坐标系中的坐标(Xren,Yren,Zren)被转换成在第一坐标系中的相对应的坐标(xren,yren,zren)。在该实施例中的第一坐标系和第二坐标系可以以任何合适的方式被定义。例如，接收到的OCT视网膜图像数据的每个像素可以具有在第一坐标系中的像素坐标，其是表示像素在构成接收到的OCT视网膜图像数据的像素的二维阵列或三维阵列或网格中的位置的整数值，每组整数值与在三维阵列的情况下的OCT扫描参数H、V和D(或者在二维阵列的情况下的H和D(可选地V和D))的相应值相关联。
如在本实施例中的，坐标确定模块120-4可以使用下列项来确定像素阵列数据的每个像素的在第一坐标系中的相应坐标：(i)第一映射，其用于将在第二坐标系中的像素阵列的像素的坐标值(Xren,Yren,Zren)转换成在上面讨论的真实世界坐标系(A’,B’,C’)中的坐标值(A’ren,B’ren,C’ren)；(ii)第二映射，其用于将在真实世界坐标系(A’,B’,C’)中的像素阵列的像素的坐标值(A’ren,B’ren,C’ren)转换成在上面讨论的另一个真实世界坐标系(A,B,C)中的坐标值(Aren,Bren,Cren)；(iii)第三映射，其用于将在真实世界坐标系(A,B,C)中的像素阵列的像素的坐标值(Aren,Bren,Cren)转换成指示在由OCT扫描仪执行的扫描内的位置的OCT扫描参数值(Hren,Vren,Dren)的相对应的一组值；以及(iv)第四映射，其用于将该组OCT扫描参数值转换成在第一坐标系中的像素阵列的像素的对应坐标值(xren,yren,zren)。如在本实施例中的，可以选择要由所渲染的图像覆盖的真实世界坐标的范围，使得在所捕获的数据中的感兴趣区域(或者也许全部)出现在所渲染的图像中。假设所捕获的数据或其感兴趣区域在概念上是数据点的长方体，一种获得要由所渲染的图像覆盖的在(A’,B’,C’)中的坐标的范围的方法是估计在所捕获的数据的这个概念长方体的八个角处的数据点的(A’ren,B’ren,C’ren)。(A’ren,B’ren,C’ren)的这些值的极值是要由所渲染的图像覆盖的(A’,B’,C’)范围的极值的合适值。
在过程S70中，插值模块130-4通过使用在过程S60中由坐标确定模块120-4确定的在第一真实坐标系(A,B,C)中的坐标或在第一坐标系中的像素阵列的像素的对应的坐标(xren,yren,zren)来计算像素阵列的像素的值，以在接收到的OCT视网膜图像数据的像素的值之间插值，从而生成所渲染的OCT视网膜图像数据。
在该实施例中执行的插值基于源数据的规则网格。因此，该实施例可以提供比上面所述的其他实施例更有效的实现，因为这些实施例每个包括其源数据在不规则网格上的插值。
尽管上面参考附图给出了本发明的示例性实施例的描述，但是示例性实施例的具体配置不限于此，并且包括在不偏离本发明的精神和范围的范围内进行设计等。
在本说明书中提及的所有出版物、专利申请和技术标准在好像每一个单独的出版物、专利申请或技术标准被特别和单独地指示以通过引用被并入的相同程度上通过引用被并入本说明书中。

Claims (20)

1.一种渲染光学相干断层扫描OCT视网膜图像数据的方法，所述OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个所述扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于所述第一坐标系的第二坐标系中的相应坐标相关联，所述方法包括：
接收(S10)由所述OCT扫描仪获取的所述OCT视网膜图像数据；
使用从在所述第一坐标系中的坐标到在所述第二坐标系中的坐标的变换来确定(S20；S20’)在所接收到的OCT视网膜图像中的像素的在所述第二坐标系中的坐标的值；以及
在所述第二坐标系中的确定的坐标的值处的像素的值之间插值(S30)以计算所述渲染的OCT视网膜图像数据的像素的值。
2.一种渲染光学相干断层扫描OCT视网膜图像数据的方法，所述OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个所述扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于所述第一坐标系的第二坐标系中的相应坐标相关联，所述方法包括：
接收(S10)由所述OCT扫描仪获取的所述OCT视网膜图像数据；以及
使用接收到的OCT视网膜图像数据处理像素阵列以通过下列操作来生成所渲染的OCT视网膜图像数据：
使用从在所述第一坐标系中的坐标到在所述第二坐标系中的坐标的变换来确定(S20)在所述接收到的OCT视网膜图像中的像素的在所述第二坐标系中的坐标的值；
在所述第二坐标系中的确定的坐标的值之间插值(S40)以计算所述像素阵列中的像素的在所述第一坐标系中的坐标的值；以及
使用(S50)所述像素阵列中的像素的在所述第一坐标系中的坐标的所计算的值来在所述接收到的OCT视网膜图像数据的像素的值之间插值，以计算在所述像素阵列中的像素的值。
3.根据权利要求1或权利要求2所述的方法，其中，使用所述变换和(i)在所述OCT视网膜图像数据的获取期间使用的至少一组OCT扫描参数值和(ii)对应于所述至少一组OCT扫描参数值的在所述第一坐标系中的坐标之间的映射来确定(S20)在所述接收到的OCT视网膜图像中的像素的在所述第二坐标系中的坐标的值，每组OCT扫描参数值指示在所述扫描内的相应位置。
4.根据权利要求3所述的方法，其中：
所述接收到的OCT视网膜图像数据包括由所述OCT扫描仪执行的对所述视网膜的A扫描的OCT视网膜图像数据；以及
通过下列操作来确定(S20)在所述A扫描的所述接收到的OCT视网膜图像数据中的像素的在所述第二坐标系中的坐标的值：
使用所述变换和所述映射来确定在所述A扫描的所述接收到的OCT视网膜图像数据中的第一像素的在所述第二坐标系中的坐标的值，所述第一像素与第一组OCT扫描参数值相关联；以及
使用所述变换和在所述第一坐标系中的指示所述OCT扫描仪根据所述第一组OCT扫描参数值操作时入射在所述眼睛的视网膜上或在所述眼睛的视网膜中传播的光束的方向的矢量来确定在所述A扫描的所述接收到的OCT视网膜图像数据中的剩余像素的在所述第二坐标系中的坐标的值。
5.根据权利要求1所述的方法，其中，确定(S20’)在所述接收到的OCT视网膜图像数据中的像素的在所述第二坐标系中的坐标的值包括：
选择(S21’)在所述接收到的OCT视网膜图像数据中基本上均匀地分布在所述接收到的OCT视网膜图像数据中的像素的子集；
使用从在所述第一坐标系中的坐标到在所述第二坐标系中的坐标的变换来为所述像素的子集的每个像素确定(S22’)在所述第二坐标系中的坐标的相应值；
使用(S23’)在所述第二坐标系中的像素的子集的像素的所确定的坐标的值来确定在所述第一坐标系中的坐标和在所述第二坐标系中的坐标之间的参数化映射的参数；以及
使用(S24’)具有所确定的参数的所述参数化映射来确定在所述接收到的OCT视网膜图像数据中的像素的在所述第二坐标系中的坐标的值。
6.根据权利要求5所述的方法，其中，使用下列项来确定(S22’)所述像素的子集的每个像素的在所述第二坐标系中的相应坐标：
第一映射，其用于将所述接收到的OCT图像数据的像素的坐标从所述第一坐标系转换成第一真实世界坐标系；
第二映射，其用于将在所述第一真实世界坐标系中的像素的坐标转换成在第二真实世界坐标系中的坐标；以及
第三映射，其用于将在所述像素阵列中的像素的坐标转换成在所述第二真实世界坐标系中的对应的坐标。
7.一种渲染光学相干断层扫描OCT视网膜图像数据的方法，所述OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个所述扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于所述第一坐标系的第二坐标系中的相应坐标相关联，所述方法包括：
接收(S10)由所述OCT扫描仪获取的所述OCT视网膜图像数据；以及
使用接收到的OCT视网膜图像数据处理像素阵列以通过下列操作来生成所述渲染的OCT视网膜图像数据：
对于所述像素阵列的每个像素，使用从在所述第二坐标系中的坐标到在所述第一坐标系中的坐标的变换来确定(S60)在所述第一坐标系中的相应坐标；以及
通过使用在所述第一坐标系中的所确定的坐标以在所述接收到的OCT视网膜图像数据的像素的值之间插值来计算(S70)所述像素阵列的像素的值。
8.根据权利要求7所述的方法，其中，使用下列项来确定(S60)所述像素阵列的每个像素的在所述第一坐标系中的相应坐标：
第一映射，其用于将所述像素阵列的像素的在所述第二坐标系中的坐标转换成第一真实世界坐标系中的坐标；
第二映射，其用于将所述像素阵列的像素的在所述第一真实世界坐标系中的坐标转换成在第二真实世界坐标系中的坐标；
第三映射，其用于将所述像素阵列的像素的在所述第二真实坐标系中的坐标转换成指示在由所述OCT扫描仪执行的所述扫描内的位置的一组OCT扫描参数值；以及
第四映射，其用于将所述一组OCT扫描参数值转换成所述像素阵列的像素的在所述第一坐标系中的坐标。
9.根据任一前述权利要求所述的方法，其中，使用公共的第二坐标系来渲染由所述OCT扫描仪获取的多个OCT视网膜图像的OCT视网膜图像数据。
10.一种存储计算机程序指令的非暂时性计算机可读存储介质(250)，所述计算机程序指令当由计算机(270)执行时使所述计算机(270)执行根据任一前述权利要求的方法。
11.一种携带计算机程序指令的信号(260)，所述计算机程序指令在由计算机(270)执行时使所述计算机(270)执行根据权利要求1到9中的任一项的方法。
12.一种用于渲染光学相干断层扫描OCT视网膜图像数据的装置(100-1；100-2)，所述OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个所述扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于所述第一坐标系的第二坐标系中的相应坐标相关联，所述装置包括：
通信模块(110)，其被布置为接收由所述OCT扫描仪获取的所述OCT视网膜图像数据；以及
坐标确定模块(120-1；120-2)，其被布置为使用从在所述第一坐标系中的坐标到在所述第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在所述第二坐标系中的坐标的值；以及
插值模块(130-1)，其被布置为在所述第二坐标系中的所确定的坐标的值处的像素的值之间插值以计算所述渲染的OCT视网膜图像数据的像素的值。
13.一种用于渲染光学相干断层扫描OCT视网膜图像数据的装置(100-3)，所述OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个所述扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于所述第一坐标系的第二坐标系中的相应坐标相关联，所述装置包括：
通信模块(110)，其被布置为接收由所述OCT扫描仪获取的所述OCT视网膜图像数据；
坐标确定模块(120-1)，其被布置为使用从在所述第一坐标系中的坐标到在所述第二坐标系中的坐标的变换来确定在接收到的OCT视网膜图像中的像素的在所述第二坐标系中的坐标的值；以及
插值模块(130-3)，其被布置为使用所述接收到的OCT视网膜图像数据处理像素阵列以通过下列操作来生成所述渲染的OCT视网膜图像数据：
在所述第二坐标系中的确定的坐标的值之间插值以计算所述像素阵列中的像素的在所述第一坐标系中的坐标的值；以及
使用所述像素阵列中的像素的在所述第一坐标系中的坐标的所计算的值来在所述接收到的OCT视网膜图像数据的像素的值之间插值，以计算在所述像素阵列中的像素的值。
14.根据权利要求12或权利要求13所述的装置(100-1)，其中，所述坐标确定模块(120-1)被布置成使用所述变换和(i)在所述OCT视网膜图像数据的获取期间使用的至少一组OCT扫描参数值和(ii)对应于所述至少一组OCT扫描参数值的所述第一坐标系中的坐标之间的映射来确定在所述接收到的OCT视网膜图像中的像素的在所述第二坐标系中的坐标的值，每组OCT扫描参数值指示在所述扫描内的相应位置。
15.根据权利要求14所述的装置(100-1)，其中：
所述通信模块(110)被布置成接收由所述OCT扫描仪执行的所述视网膜的A扫描的OCT视网膜图像数据；以及
所述坐标确定模块(120-1)被布置成通过下列操作来确定在所述A扫描的所述接收到的OCT视网膜图像数据中的像素的在所述第二坐标系中的坐标的值：
使用所述变换和所述映射来确定在所述A扫描的所述接收到的OCT视网膜图像数据中的第一像素的在所述第二坐标系中的坐标的值，所述第一像素与第一组OCT扫描参数值相关联；以及
使用所述变换和在所述第一坐标系中的指示所述OCT扫描仪根据所述第一组OCT扫描参数值操作时入射在所述眼睛的所述视网膜上的或在所述眼睛的所述视网膜中传播的光束的方向的矢量来确定在所述A扫描的所述接收到的OCT视网膜图像数据中的剩余像素的在所述第二坐标系中的坐标的值。
16.根据权利要求12所述的装置(100-2)，其中，所述坐标确定模块(120-2)被布置成通过下列操作来确定在所述接收到的OCT视网膜图像数据中的像素的在所述第二坐标系中的坐标的值：
选择(S21’)在所述接收到的OCT视网膜图像数据中基本上均匀地分布在所述接收到的OCT视网膜图像数据中的像素的子集；
使用从在所述第一坐标系中的坐标到在所述第二坐标系中的坐标的变换来为所述像素的子集的每个像素确定(S22’)在所述第二坐标系中的坐标的相应值；
使用(S23’)在所述第二坐标系中的像素的子集的像素的所确定的坐标的值来确定在所述第一坐标系中的坐标和在所述第二坐标系中的坐标之间的参数化映射的参数；以及
使用(S24’)具有所确定的参数的所述参数化映射来确定在所述接收到的OCT视网膜图像数据中的像素的在所述第二坐标系中的坐标的值。
17.根据权利要求16所述的装置(100-2)，其中，所述坐标确定模块(120-2)使用下列项来确定所述像素的子集的每个像素的在所述第二坐标系中的坐标的相应值：
第一映射，其用于将所述接收到的OCT图像数据的像素的坐标从所述第一坐标系转换成第一真实世界坐标系；
第二映射，其用于将在所述第一真实世界坐标系中的像素的坐标转换成在第二真实世界坐标系中的坐标；以及
第三映射，其用于将在所述像素阵列中的像素的坐标转换成在所述第二真实世界坐标系中的对应的坐标。
18.一种用于渲染光学相干断层扫描OCT视网膜图像数据的装置(100-4)，所述OCT视网膜图像数据由在一定范围的扫描位置上扫描眼睛的视网膜的OCT扫描仪获取，其中，每个所述扫描位置与在第一坐标系中的相应坐标相关联，并且所渲染的OCT视网膜图像数据的每个像素与在不同于所述第一坐标系的第二坐标系中的相应坐标相关联，所述装置包括：
通信模块(110)，其被布置为接收由所述OCT扫描仪获取的所述OCT视网膜图像数据；以及
坐标确定模块(120-4)，其被布置为：对于要被处理以生成所渲染的OCT视网膜图像数据的像素阵列的每个像素，使用从在所述第二坐标系中的坐标到在所述第一坐标系中的坐标的变换来确定在所述第一坐标系中的坐标的相应值；以及
插值模块(130-4)，其被布置为通过使用在所述第一坐标系中的所确定的坐标的值以在所述像素阵列的像素的值之间插值来计算所述像素阵列的像素的值。
19.根据权利要求18所述的装置(100-4)，其中，所述坐标确定模块(120-4)被布置成使用下列项来确定所述像素阵列的每个像素的在所述第一坐标系中的坐标的相应值：
第一映射，其用于将所述像素阵列的像素的在所述第二坐标系中的坐标转换成在第一真实世界坐标系中的坐标；
第二映射，其用于将所述像素阵列的像素的在所述第一真实世界坐标系中的坐标转换成在第二真实世界坐标系中的坐标；
第三映射，其用于将所述像素阵列的像素的在所述第二真实世界坐标系中的坐标转换成指示在由所述OCT扫描仪执行的所述扫描内的位置的一组OCT参数值；以及
第四映射，其用于将所述一组OCT扫描参数值转换成所述像素阵列的像素的在所述第一坐标系中的坐标。
20.根据权利要求12至19中的任一项所述的装置(100-1；100-2；100-3；100-4)，其中，所述装置被布置成使用公共的第二坐标系来渲染由所述OCT扫描仪获取的多个OCT视网膜图像的OCT视网膜图像数据。
CN201910383653.5A 2018-05-11 2019-05-09 Oct图像处理 Active CN110473265B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP18171806.5A EP3566637A1 (en) 2018-05-11 2018-05-11 Oct image processing
EPEP18171806.5 2018-05-11

Publications (2)

Publication Number Publication Date
CN110473265A true CN110473265A (zh) 2019-11-19
CN110473265B CN110473265B (zh) 2023-06-16

Family

ID=62152436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910383653.5A Active CN110473265B (zh) 2018-05-11 2019-05-09 Oct图像处理

Country Status (10)

Country Link
US (1) US10929963B2 (zh)
EP (2) EP4085819B1 (zh)
JP (3) JP6841863B2 (zh)
KR (1) KR102199134B1 (zh)
CN (1) CN110473265B (zh)
AU (1) AU2019202426B2 (zh)
CA (2) CA3091570C (zh)
DK (1) DK4085819T3 (zh)
ES (1) ES2969279T3 (zh)
FI (1) FI4085819T3 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257282A (zh) * 2020-01-22 2020-06-09 杭州捷诺飞生物科技股份有限公司 Oct成像系统、电子设备和机器可读存储介质
CN114549575A (zh) * 2020-11-12 2022-05-27 奥普托斯股份有限公司 体积oct图像数据处理

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4052022A1 (en) 2019-10-28 2022-09-07 Ricoh Company, Ltd. Measuring apparatus and biological information measuring apparatus
EP3889889A1 (en) * 2020-03-30 2021-10-06 Optos PLC Ocular image data processing
CN112785640B (zh) * 2020-12-28 2022-08-09 宁波江丰生物信息技术有限公司 一种用于扫描仪内部切片位置检测的方法和系统
KR102464056B1 (ko) * 2021-09-10 2022-11-09 가톨릭대학교 산학협력단 Mri 영상과 oct 영상의 병합을 이용한 안구 및 시신경의 3차원 모델링 방법
CN114184581B (zh) * 2021-10-25 2024-05-07 北京理工大学 基于oct系统的图像优化方法、装置、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150182111A1 (en) * 2013-11-29 2015-07-02 Nidek Co., Ltd. Optical tomographic image photographing apparatus and optical tomographic image photographing program
CN104958061A (zh) * 2015-07-28 2015-10-07 北京信息科技大学 双目立体视觉三维成像的眼底oct成像方法及其系统
JP2017113146A (ja) * 2015-12-22 2017-06-29 株式会社トプコン 眼科撮影装置
JP2018051177A (ja) * 2016-09-30 2018-04-05 キヤノン株式会社 眼科撮影装置、及び眼科撮影装置の作動方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003011764A2 (en) 2001-08-03 2003-02-13 Volker Westphal Real-time imaging system and method
JP5355074B2 (ja) * 2008-12-26 2013-11-27 キヤノン株式会社 ３次元形状データ処理装置、３次元形状データ処理方法及びプログラム
US8693745B2 (en) 2009-05-04 2014-04-08 Duke University Methods and computer program products for quantitative three-dimensional image correction and clinical parameter computation in optical coherence tomography
PL2563206T3 (pl) * 2010-04-29 2018-12-31 Massachusetts Institute Of Technology Sposób i urządzenie do korekcji ruchu i poprawy jakości obrazu w optycznej tomografii koherencyjnej
JP5794664B2 (ja) 2011-01-20 2015-10-14 キヤノン株式会社 断層画像生成装置及び断層画像生成方法
ES2391510B1 (es) 2011-04-29 2013-11-11 Consejo Superior De Investigaciones Científicas (Csic) Procedimiento de calibracion y correccion de la distorsion de barrido de un sistema de tomografia de coherencia optica
JP6115007B2 (ja) 2012-01-31 2017-04-19 株式会社ニデック 眼科画像処理装置及びプログラム
GB201217538D0 (en) 2012-10-01 2012-11-14 Optos Plc Improvements in or relating to scanning laser ophthalmoscopes
US10568586B2 (en) * 2012-10-05 2020-02-25 Volcano Corporation Systems for indicating parameters in an imaging data set and methods of use
JP6808381B2 (ja) 2016-07-07 2021-01-06 キヤノン株式会社 保持装置、投影光学系、露光装置、および物品製造方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150182111A1 (en) * 2013-11-29 2015-07-02 Nidek Co., Ltd. Optical tomographic image photographing apparatus and optical tomographic image photographing program
CN104958061A (zh) * 2015-07-28 2015-10-07 北京信息科技大学 双目立体视觉三维成像的眼底oct成像方法及其系统
JP2017113146A (ja) * 2015-12-22 2017-06-29 株式会社トプコン 眼科撮影装置
JP2018051177A (ja) * 2016-09-30 2018-04-05 キヤノン株式会社 眼科撮影装置、及び眼科撮影装置の作動方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257282A (zh) * 2020-01-22 2020-06-09 杭州捷诺飞生物科技股份有限公司 Oct成像系统、电子设备和机器可读存储介质
CN111257282B (zh) * 2020-01-22 2022-11-15 杭州捷诺飞生物科技股份有限公司 Oct成像系统、电子设备和机器可读存储介质
CN114549575A (zh) * 2020-11-12 2022-05-27 奥普托斯股份有限公司 体积oct图像数据处理

Also Published As

Publication number Publication date
JP2021079144A (ja) 2021-05-27
CA3091570A1 (en) 2019-11-11
KR20190129760A (ko) 2019-11-20
US20190347774A1 (en) 2019-11-14
JP7251028B2 (ja) 2023-04-04
JP6841863B2 (ja) 2021-03-10
DK4085819T3 (da) 2024-01-22
CA3040403C (en) 2022-03-22
EP4085819B1 (en) 2023-12-20
US10929963B2 (en) 2021-02-23
EP4085819A1 (en) 2022-11-09
CN110473265B (zh) 2023-06-16
KR102199134B1 (ko) 2021-01-06
JP2019195636A (ja) 2019-11-14
AU2019202426A1 (en) 2019-11-28
AU2019202426B2 (en) 2020-02-13
CA3040403A1 (en) 2019-11-11
ES2969279T3 (es) 2024-05-17
EP3566637A1 (en) 2019-11-13
CA3091570C (en) 2023-07-18
FI4085819T3 (fi) 2024-01-17
JP2022033344A (ja) 2022-02-28

Similar Documents

Publication Publication Date Title
CN110473265A (zh) Oct图像处理
CN108926355A (zh) 用于站立对象的x射线系统和方法
JPH05509015A (ja) 物体のトポグラフィーを得るための方法及び装置
JP2013504752A (ja) 非接触物体検査
JP2019088823A (ja) 眼の画像処理の改良
US8526570B2 (en) Computed tomography method, computer software, computing device and computed tomography system for determining a volumetric representation of a sample
CN113456018B (zh) 眼睛图像数据处理
JP6336271B2 (ja) 検眼鏡の改良
WO2020037582A1 (en) Graph-based key frame selection for 3-d scanning
WO2014096835A1 (en) Improvements in and relating to ophthalmoscopes
WO2020090439A1 (ja) 画像処理装置、画像処理方法およびプログラム
CN113641318B (zh) 显示数据校准方法及系统
WO2023181357A1 (ja) 光干渉断層画像生成装置、光干渉断層画像生成方法、及び、記録媒体
JP6701277B2 (ja) 情報処理装置、情報処理方法、プログラム
Sadhu et al. Automatic calibration of time of flight based non-line-of-sight reconstruction
CN104997482B (zh) 一种直线高清扫描随访实现方法及装置
WO2020152338A1 (en) Method for selecting an illumination pattern for a fluorescence molecular tomography measurement

Legal Events

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