CN104132649A - 眼睛成像中的改进以及与眼睛成像相关的改进 - Google Patents

眼睛成像中的改进以及与眼睛成像相关的改进 Download PDF

Info

Publication number
CN104132649A
CN104132649A CN201410185236.7A CN201410185236A CN104132649A CN 104132649 A CN104132649 A CN 104132649A CN 201410185236 A CN201410185236 A CN 201410185236A CN 104132649 A CN104132649 A CN 104132649A
Authority
CN
China
Prior art keywords
dimensional representation
retina
dimensional
retina part
coordinate
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
CN201410185236.7A
Other languages
English (en)
Other versions
CN104132649B (zh
Inventor
J·V·赫默特
M·维罗克
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
Photoelectron Corp
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 Photoelectron Corp filed Critical Photoelectron Corp
Priority to CN201711373469.XA priority Critical patent/CN108371543B/zh
Publication of CN104132649A publication Critical patent/CN104132649A/zh
Application granted granted Critical
Publication of CN104132649B publication Critical patent/CN104132649B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0016Operational features thereof
    • A61B3/0025Operational features thereof characterised by electronic signal processing, e.g. eye models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/14Arrangements specially adapted for eye photography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0008Apparatus for testing the eyes; Instruments for examining the eyes provided with illuminating means
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/1005Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for measuring distances inside the eye, e.g. thickness of the cornea
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/1025Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for confocal scanning
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/64Analysis of geometric attributes of convexity or concavity
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30041Eye; Retina; Ophthalmic

Abstract

一种眼睛成像中的改进以及与眼睛成像相关的改进。一种确定眼睛视网膜的几何测量的方法,包括:获取眼睛视网膜的至少一部分的二维表示(34),导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射(36),利用视网膜部分的二维表示的一个或多个坐标来在二维表示上定义视网膜的将被采取的几何测量(38),利用几何重映射将视网膜部分的二维表示的所述坐标或每个坐标转换成视网膜部分的三维表示的等效坐标(40),以及利用视网膜部分的三维表示的所述等效坐标或每个等效坐标来确定眼睛视网膜的几何测量(42)。

Description

眼睛成像中的改进以及与眼睛成像相关的改进
技术领域
本发明涉及眼睛成像中的改进以及与眼睛成像相关的改进,具体涉及对眼睛的几何测量的确定。
背景技术
存在各种类型的可用来获取眼睛图像的成像装置,例如检眼镜和眼底照相机。具体地,这种装置测量近似为球体的三维结构的眼睛视网膜。成像装置的输出是三维视网膜的二维图像,因为二维图像具有多种优势,而且它们可被呈现在计算机屏幕上并被输出为图片等。成像装置检查(例如扫描)三维视网膜并产生视网膜的二维图像。在检查三维视网膜时,装置引入了失真,该失真由装置的光学和机械特性确定。失真包括视网膜的二维图像的一个或多个特性,例如距离、面积以及夹角。如果要求视网膜的几何测量,例如视网膜的两个结构之间的距离,而且该成像装置产生的二维图像被使用,则结果将会失真并且不精确。
发明内容
根据本发明的一个方面,提供了一种确定眼睛视网膜的几何测量的方法,包括:
获取眼睛视网膜的至少一部分的二维表示,
导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射,
利用视网膜部分的二维表示的一个或多个坐标来在二维表示上定义视网膜的将被采取的几何测量,
利用几何重映射将视网膜部分的二维表示的所述坐标或每个坐标转换成视网膜部分的三维表示的等效坐标,以及
利用视网膜部分的三维表示的所述等效坐标或每个等效坐标来确定眼睛视网膜的几何测量。
导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射可包括:确定针对在产生视网膜部分的二维表示时引入视网膜部分的二维表示的失真的校正。失真可通过用来产生二维表示的成像装置的光学和机械特性而引入视网膜部分的二维表示。
确定针对二维表示中的失真的校正可包括对成像装置的光学和机械特性建模。
对成像装置的光学和机械特性建模可包括:
(i)构建包括成像装置和模型眼的成像系统的光学描述,
(ii)通过成像系统将射线传递至模型眼的视网膜表面上,
(iii)计算视网膜表面处的射线的实际测量结果,
(iv)针对射线确定成像系统的水平扫描角度和垂直扫描角度,
(v)利用成像系统的水平扫描角度和垂直扫描角度计算视网膜表面处的射线的期望测量结果,
(vi)针对多个其它射线重复步骤(ii)至(v),以及
(vii)将视网膜表面处的射线的实际测量结果与视网膜表面处的射线的相应期望测量结果进行比较以确定针对成像装置的眼睛表示中的失真的校对。
构建包括成像装置和模型眼的成像系统的光学描述可包括确定成像装置的光路特性,确定模型眼的光路特性并且将光路特性关联起来以给出成像系统的光路特性。确定成像装置的光路特性可包括确认成像装置的具有光路效应的组件,确认成像装置中组件的顺序,建立数学函数以描述每个组件随时间变化的光路特性,按顺序关联组件的光路特性以给出成像装置的光路特性。确定模型眼的光路特性可包括确认模型眼中具有光路效应的组件,确认模型眼中组件的顺序,建立数学函数以描述每个组件的光路特性,按顺序关联组件的光路特性以给出模型眼的光路特性。
通过成像系统将射线传递至模型眼的视网膜表面上可包括使用射线追踪系统。射线追踪系统可以是可商用的射线追踪系统,例如Zemax。使用射线追踪系统可包括将成像系统的光学描述载入射线追踪系统并针对每个射线确定通过成像系统的路径。计算在视网膜表面处的每个射线的实际测量结果可包括计算每个射线与视网膜表面的交叉点的坐标。
确定针对射线的水平扫描角度可包括确认用来产生射线的成像装置的水平扫描元件的角位置,并且利用其来计算针对射线的水平扫描角度。确定针对射线的垂直扫描角度可包括确认用来产生射线的成像装置的垂直扫描元件的角位置,并且利用其来计算针对射线的垂直扫描角度。
算视网膜表面处的每个射线的期望测量结果可包括使用成像系统的水平扫描角度和垂直扫描角度来计算个射线与视网膜表面的交叉点的坐标。
确定针对成像装置的眼睛表示中的失真的校对可包括导出解析变换,该解析变换将模型眼的视网膜表面处的射线的实际测量结果映射至模型眼的视网膜表面处的射线的相应期望测量结果。确定针对成像装置的眼睛表示中的失真的校对可包括构建校正查找表(LUT),该校正查找表包括针对多个射线的每个的模型眼的视网膜表面处的射线的实际位置相对于模型眼的视网膜表面处的射线的实际位置。
用于确定针对二维表示中的失真的校对的方法可进一步包括针对模型眼的多个非零注视角确定针对成像装置的眼睛表示中的失真的校对。对于每个非零注视角,其还可包括模型眼的中央凹(fovial)位置来测量注视角。可针对多个非零注视角中的每一个来构建校正LUT。
确定针对二维表示中的失真的校正可包括:
建立目标,
计算目标的图像,
利用用于产生二维表示的成像装置来扫描目标以产生目标的失真图像,以及
将目标的计算出的图像与目标的失真图像进行比较以确定针对成像装置的图像失真的校正。
计算目标的图像可包括计算目标的图像的多个点的坐标。扫描目标可包括通过成像装置传递多个射线以确定目标的失真图像的多个点的坐标。
将目标的计算出的图像与目标的失真图像进行比较可包括将目标的计算出的图像的点与目标的失真图像相应点进行比较。
确定针对成像装置的图像失真的校正可包括导出解析变换,该解析变换将目标的失真图像的点的坐标映射至目标的计算出的图像的相应点的坐标。确定针对成像装置的图像失真的校正可包括构建校正查找表,该校正查找表包括针对目标的失真图像的多个可能点中的每个的失真图像的可能点的坐标相对于目标的计算出的图像的相应点的坐标。
目标可包括多个对比形状。目标上可具有预定图案。
确定针对二维表示中的失真的校正可包括:接收眼睛的图像,并且变换图像直到利用失真校正已知的成像装置所产生的眼睛图像对其进行登记。
确定针对二维表示中的失真的校正可包括第三方,其确定校正并接收来自第三方的校正。
将视网膜部分的二维表示转换成视网膜部分的三维表示的几何重映射可具有解析方程的形式,该解析方程将二维表示的坐标变换为三维表示的等效坐标。将视网膜部分的二维表示转换成视网膜部分的三维表示的几何重映射可具有查找表的形式,该查找表列出与三维表示的坐标等效的二维表示的坐标。三维表示的坐标可以是Cartesian坐标或者可以是球体坐标。
用来定义几何测量的视网膜部分的二维表示的一个或多个坐标可定义二维表示上的一个或多个点。可利用Cartesian坐标空间来识别二维表示上的所述点或每个点。可利用诸如鼠标之类的点击装置来在屏幕上识别二维表示上的所述点或每个点。
利用几何重映射将视网膜部分的二维表示的一个或多个坐标转换成视网膜部分的三维表示的等效坐标包括:利用几何重映射将二维表示的一个或多个点的坐标映射至二维表示的一个或多个等效点的坐标。
视网膜部分的三维表示的所述等效坐标或每个等效坐标可包括三维表示上的一个或多个点。可利用Cartesian坐标空间来识别三维表示上的所述点或每个点。可利用球体坐标空间来识别三维表示上的所述点或每个点。
视网膜部分的几何测量可包括视网膜部分的结构的距离测量。可利用二维表示的第一和第二点的坐标在视网膜部分的二维表示上定义距离几何测量。几何重映射可被用来将视网膜部分的二维表示的第一和第二点转换成视网膜部分的三维表示的第一和第二等效点。利用三维表示的第一和第二等效点来确定距离几何测量可包括利用三维表示的第一点来定义测量起始位置并利用三维表示的第二点来定义测量结束位置,以及测量三维表示的第一点和第二点之间的距离。视网膜部分的三维表示可被看作是球体以及被测量作为球体的第一和第二点之间的最短距离的第一点和第二点之间的距离。
可利用余弦方法的球形定理以相对于球体的单位(度数或弧度)测量球体上的第一和第二点之间的最短距离,即
其中Δσ是点之间的圆心角,分别是第一和第二点的经纬度,Δλ是经度的绝对差。
以相对于球体的单位测量的球体上的第一和第二点之间的最短距离可利用下式转换成物理单位(mm)
d=rΔσ
其中d是点之间的距离,r是球体的半径,Δσ是点之间的圆心角并且以弧度给出。
可利用vincenty公式以相对于球体的单位(度数或弧度)测量球体上的第一和第二点之间的最短距离,即
Δσ = arctan ( ( cos φ f sin Δλ ) 2 + ( cos φ s sin φ f - sin φ s cos φ f cos Δλ ) 2 sin φ s sin φ f + cos φ s cos φ f cos Δλ )
其中Δσ是点之间的圆心角,分别是第一和第二点的经纬度,Δλ是经度的绝对差。
以相对于球体的单位测量的球体上的第一和第二点之间的最短距离可利用下式转换成物理单位(mm)
d=rΔσ
其中d是点之间的距离,r是球体的半径,Δσ是点之间的圆心角并且以弧度给出。
可向用户输出距离几何测量。可在视网膜部分的二维表示上表达距离几何测量。三维球体上的第一和第二点之间的最短距离(测地线)不是必须对应于视网膜部分的二维表示上的直线。在二维表示上表达距离几何测量可包括使球体上的第一和第二点之间的距离参数化,计算多个中间点并利用该点来在二维表示上表达距离几何测量。
视网膜部分的几何测量可包括视网膜部分的第一和第二结构之间的角度。利用二维表示的第一点和第二点的坐标来标记第一结构而且利用二维表示的第一点和第三点坐标来标记第二结构,角度几何测量可定义在视网膜部分的二维表示上。几何重映射可被用于将视网膜部分的二维表示的第一、第二和第三点转换成视网膜部分的三维表示的第一、第二和第三点等效点。利用三维表示的等效点来确定角度几何测量可包括利用三维表示的第一点和第二点来标记第一结构以及利用三维表示的第一点和第三点来标记第二结构,并且通过三维表示中的第一点来测量第二点和第三点之间的角度。可利用haversines定理来测量角度C,即
haversin(c)=haversin(a-b)+sin(a)sin(b)haversin(C)
其中haversin(x)≡sin2(x/2),x1和x2是标记第一结构的三维表示上的第一和第二点,x1和x3是标记第二结构的三维表示上的点,a是x1和x2之间的测地线的弧长,b是x1和x3之间的测地线的弧长,c是x2和x3之间的测地线的弧长。
可向用户输出角度几何测量。可通过表现二维表示上的角度的边(x1和x2之间的测地线以及x1和x3之间的测地线),在视网膜部分的二维表示上表示角度几何测量。
视网膜部分的几何测量可包括视网膜部分的结构的面积。可利用二维表示的多个点的坐标在视网膜部分的二维表示上定义面积几何测量。几何重映射可被用于将视网膜部分的二维表示的多个点转换成视网膜部分的三维表示的等效多个点。利用三维表示的多个等效点来确定面积几何测量可包括利用三维表示的多个点来定义结构的形状以及测量三维表示中形状的面积。视网膜的一部分的三维表示可被认为是球体,而且结构的形状可被定义为多边形,而且可利用下式测量多边形的面积
A = r 2 ( Σ i = 1 n α i - ( n - 2 ) π )
其中A是多边形的面积,αi(i=1,…,n)是多边形内部的n个角度,r是球体的半径。得到物理单位(例如mm2,如果r的单位是mm)。如果在上述公式中省略r2,则得到相对于球体的单位(球面度(sr),立体角的单位)的结果。
使用三维表示的等效多个点来确定面积几何测量可包括:使用三维表示上的多个像素来定义结构的形状,并且通过求取像素的面积之和来测量形状的面积。可通过假设每个像素是具有四个点的多边形(每个点表示像素的一个角落)来计算三维表示上的像素的面积。
当视网膜的一部分的三维表示被看作是球体时,可通过测量眼睛直径来确定球体的半径。当视网膜的一部分的三维表示被认为是球体时,球体的半径可被设置成人眼的平均半径,这大概是12mm。即使眼睛的半径未知或者不能被近似,也可通过假设眼睛是单位球体来进行相对于视网膜的球体的尺寸的相对测量。可在弧度或度数以及球面度中的面积或平方度中测量球体上的距离。
获取眼睛视网膜的至少一部分的二维表示可包括操作成像装置以产生二维表示。获取眼睛视网膜的至少一部分的二维表示可包括接收成像装置已经产生的二维表示。
获取眼睛视网膜的至少一部分的二维表示可包括接收经过失真校正的视网膜部分的三维表示并且利用投射来从三维表示获取二维表示。投射可以从多个已知投射中选择。投射可以是任意保形投射(保持二维表示中的角度)、等距投影(保持二维表示中的距离)、面积保持投影(保持二维表示中的面积)。导出将视网膜部分的二维表示转换成视网膜部分的三维表示的几何重映射可包括确定获取视网膜部分的二维表示时使用的投影的倒转。
根据本发明的第二方面,提供了一种存储程序指令的计算机可读介质,在执行程序指令时执行根据本发明的第一方面的方法。
根据本发明的第三方面,提供了一种用于确定眼睛视网膜的几何测量的成像系统,包括:
成像装置,其获取眼睛视网膜的至少一部分的二维表示,
导出元件,其导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射,
测量定义元件,其利用视网膜部分的二维表示的一个或多个坐标来在二维表示上定义视网膜的将被采取的几何测量,
重映射元件,其利用几何重映射将视网膜部分的二维表示的所述坐标或每个坐标转换成视网膜部分的三维表示的等效坐标,以及
测量确定元件,其利用视网膜部分的三维表示的所述等效坐标或每个等效坐标来确定眼睛视网膜的几何测量。
测量定义元件可被提供为成像装置中的软件或软件及硬件。导出元件、重映射元件和测量确定元件可被提供为成像装置中的软件。
在不同传统成像装置上成像的相同眼睛将产生不能直接比较或重叠的视网膜表示。不同注视角下相同的传统成像装置产生的相同眼睛表示也是这样。在本发明中,一旦变换成三维球形表示,则可以以物理单位测量在成像装置上等效的视网膜结构,因此允许直接比较和关系的导出。例如,测量不同成像装置产生的同一眼睛的两个图像上的视神经乳头的中心和中央凹之间的距离,应该得到相同的距离。由于不同光学和机械组件以及由于成像时的注视角,每个成像装置将具有不同失真。通过测量相同球形表示中的距离,可以补偿这些差异。
眼睛视网膜的一部分的几何测量可用于:
-视网膜部分中的病变/创伤幻像的归类,
-不同成像装置之间的视网膜部分结构的测量的比较,
-眼睛的视网膜部分中的解剖结构和病理中的几何变化(例如,尺寸、方向、距离中的变化)的径向追踪,
-基于几何测量的治疗计划的制定,
-基于从几何测量导出的几何坐标的诸如激光凝固法之类的治疗的直接传递,
-基于人群解剖结构的几何测量的分布的规范数据库的创建。
附图说明
现在将参考附图仅仅通过实例的方式描述本发明的实施例,其中:
图1是根据本发明的第三方面的成像系统的示意表示,以及
图2是根据本发明的第一方面的方法的流程图。
具体实施方式
参考图1,成像系统包括成像装置10、导出元件26、测量定义元件28、重映射元件30和测量确定元件32。在该实施例中,测量定义元件被提供为软件和硬件,导出元件、几何重映射元件和测量确定元件被提供为成像装置中的软件。
成像装置10包括检眼镜。检眼镜包括发射光束13的光源12、扫描中继元件,扫描中继元件包括第一扫描元件14、第二扫描元件16、扫描补偿元件18和扫描传递元件20。第一扫描元件14包括旋转多边形镜,第二扫描元件16包括振荡平面镜。扫描补偿元件18包括椭球面镜,扫描传递元件20包括非球面镜。
光源12将入射光束13引导至第一扫描元件14。这就产生了第一垂直方向上的光束扫描(如射线A,B和C所示)。入射光束撞击在扫描补偿元件18上并从其反射至第二扫描元件16。这就产生了入射光束在第二水平方向上的扫描。入射光束随后撞击在具有两个焦点的扫描传递元件20上,第二扫描元件16被布置在第一焦点处,对象的眼睛22被布置在第二焦点处。撞击在扫描传递元件20上的来自第二扫描元件16的入射光束将被引导至眼睛22,并且将撞击在眼睛视网膜的一部分上。光源12和检眼镜10的扫描中继元件组合以提供来自外部点光源的入射光束13的二维扫描并将入射光束的二维扫描从外部点光源传递至眼睛视网膜。当在视网膜上扫描入射光束时,其将从中反射以产生反射光束,反射光束通过检眼镜10的元件传递回去并被一个或多个检测器(未示出)接收。为了获取对象的眼睛22的视网膜的一部分的表示,在彼此垂直操作的第一和第二扫描元件14,16提供的光栅扫描模式下在视网膜部分上扫描来自光源12的入射光束,反射光束被一个或多个检测器接收。
参考图1和2,成像装置10获取三维眼睛视网膜22的一部分的二维表示(步骤34)。导出元件26接收视网膜部分的二维表示,并用其来导出几何重映射,几何重映射将视网膜部分的二维表示转换成视网膜部分的三维表示(步骤36)。导出几何重映射包括确定对用于产生二维表示的成像装置的光学和机械特性所引入视网膜部分的二维表示的失真的校正。
在第一实施例中,确定对二维表示中的失真的校正包括对成像装置的光学和机械特性10建模。这包括以下步骤。包括成像装置10和模型眼的成像系统的光学描述被首先构建。这包括确定成像装置10的光路特性,确定模型眼的光路特性并关联光路特性以给出成像系统的光路特性。确定成像装置10的光路特性包括确认成像装置10的具有光路效应的组件,确认成像装置10中组件的顺序,建立数学函数以描述每个组件随时间变化的光路特性,按顺序关联组件的光路特性以给出成像装置10的光路特性。确定模型眼的光路特性可包括确认模型眼中具有光路效应的组件,确认模型眼中组件的顺序,建立数学函数以描述每个组件的光路特性,按顺序关联组件的光路特性以给出模型眼的光路特性。
成像系统的光学描述被加载至射线追踪系统,例如Zemax,射线通过成像系统传递至模型眼的视网膜表面,而且确定射线通过成像系统的路径。视网膜表面处射线的实际测量结果被计算为射线与视网膜表面的交叉点的坐标。
随后针对射线确定成像系统的水平扫描角度和垂直扫描角度。这包括确认用于产生射线的成像装置的水平扫描元件的角位置并用其来计算射线的水平扫描角度,以及确认用于产生射线的成像装置的垂直扫描元件的角位置并用其来计算射线的垂直扫描角度。
利用成像系统的水平扫描角度和垂直扫描角度来计算视网膜表面处射线的期望测量结果。这包括利用成像系统的水平扫描角度和垂直扫描角度来计算射线与视网膜表面的交叉点的坐标。
随后,针对多个其它射线重复上述步骤。视网膜表面处的射线的实际测量结果与视网膜表面处的射线的相应期望测量结果进行比较,以确定对成像装置10的眼睛表示中的失真的校正。失真校正可采取解析变换或校正查找表(LUT)的形式,解析变换将模型眼的视网膜表面处的射线的实际测量结果映射至模型眼的视网膜表面处的射线的相应期望测量结果,校正查找表(LUT)包括针对多个射线的每一个的模型眼的视网膜表面处的射线实际位置相对于模型眼的视网膜表面处的射线期望位置。确定失真校正可进一步包括针对模型眼的多个非零注视角确定失真校正。
在第二实施例中,确定针对二维表示中的失真的校正包括以下步骤。首先,创建包括多个对比形状的目标。目标的图像随后被计算出来。用于产生二维表示的成像装置10被用于扫描目标以创建对象的失真图像。最后,目标的计算出来的图像与目标的失真图像比较以确定对成像装置10的图像失真的校正。
随后,失真校正被用来导出几何重映射,其将视网膜部分的二维表示映射至视网膜部分的三维表示。几何重映射采取解析公式的形式,解析公式将二维表示的坐标映射至三维表示的等效坐标。
视网膜部分的二维表示被测量定义元件28读取,并且距离几何测量被利用二维表示的第一和第二点的坐标而定义在视网膜部分的二维表示上(步骤38)。这包括向用户显示视网膜部分的二维表示,该用户利用诸如鼠标之类的点击装置识别二维表示上的第一和第二点。
随后,重映射元件30使用几何重映射将视网膜部分的二维表示的第一和第二点映射至视网膜部分的三维表示的等效第一和第二点(步骤40)。这包括使用变换来将二维表示的点的坐标映射至三维表示的等效点的坐标。利用球形坐标空间识别三维表示上的点。
随后,测量确定元件32使用三维表示的第一和第二等效点来确定距离几何测量(步骤40)。这包括使用三维表示的第一点来定义测量起始位置以及使用三维表示的第二点来定义测量结束位置,并且测量三维表示的第一点和第二点之间的距离(眼睛视网膜22的结构的长度)。视网膜部分的三维表示被看作是球体,而且第一点和第二点之间的距离被测量作为球体的第一和第二点之间的最短距离。
可利用余弦方法的球形定理以相对于球体的单位(度数或弧度)测量球体上的第一和第二点之间的最短距离,即
其中Δσ是点之间的圆心角,分别是第一和第二点的经纬度,Δλ是经度的绝对差。
以相对于球体的单位测量的球体上的第一和第二点之间的最短距离可利用下式转换成物理单位(mm)
d=rΔσ
其中d是点之间的距离,r是球体的半径,Δσ是点之间的圆心角并且以弧度给出。
通过使球体上的第一和第二点之间的距离参数化,计算多个中间点并投射该点以在二维表示上表达距离几何测量,距离几何测量被表示在视网膜部分的二维表示上。随后,向用户显示视网膜部分的二维表示上展示的距离几何测量。

Claims (18)

1.一种确定眼睛视网膜的几何测量的方法,包括:
获取眼睛视网膜的至少一部分的二维表示,
导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射,
利用视网膜部分的二维表示的一个或多个坐标来在二维表示上定义视网膜的将被采取的几何测量,
利用几何重映射将视网膜部分的二维表示的所述坐标或每个坐标转换成视网膜部分的三维表示的等效坐标,以及
利用视网膜部分的三维表示的所述等效坐标或每个等效坐标来确定眼睛视网膜的几何测量。
2.根据权利要求1所述的方法,其中导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射包括:确定针对在产生视网膜部分的二维表示时引入视网膜部分的二维表示的失真的校正。
3.根据权利要求2所述的方法,其中通过用来产生二维表示的成像装置的光学和机械特性将失真引入视网膜部分的二维表示,而且确定针对二维表示中的失真的校正包括对成像装置的光学和机械特性建模。
4.根据权利要求3所述的方法,其中对成像装置的光学和机械特性建模包括:
(i)构建包括成像装置和模型眼的成像系统的光学描述,
(ii)通过成像系统将射线传递至模型眼的视网膜表面上,
(iii)计算视网膜表面处的射线的实际测量结果,
(iv)针对射线确定成像系统的水平扫描角度和垂直扫描角度,
(v)利用成像系统的水平扫描角度和垂直扫描角度计算视网膜表面处的射线的期望测量结果,
(vi)针对多个其它射线重复步骤(ii)至(v),以及
(vii)将视网膜表面处的射线的实际测量结果与视网膜表面处的射线的相应期望测量结果进行比较以确定针对成像装置的眼睛表示中的失真的校对。
5.根据权利要求4所述的方法,其中确定针对射线的水平扫描角度包括:确认用来产生射线的成像装置的水平扫描元件的角位置,并且利用其来计算针对射线的水平扫描角度。
6.根据权利要求4所述的方法,其中确定针对射线的垂直扫描角度包括:确认用来产生射线的成像装置的垂直扫描元件的角位置,并且利用其来计算针对射线的垂直扫描角度。
7.根据权利要求3至6中任意一项所述的方法,进一步包括针对模型眼的多个非零注视角确定针对成像装置的眼睛表示中的失真的校对。
8.根据权利要求2所述的方法,其中确定针对二维表示中的失真的校正包括:
建立目标,
计算目标的图像,
利用用于产生二维表示的成像装置来扫描目标以产生目标的失真图像,以及
将目标的计算出的图像与目标的失真图像进行比较以确定针对成像装置的图像失真的校正。
9.根据前述任意一项权利要求所述的方法,其中用来定义几何测量的视网膜部分的二维表示的一个或多个坐标可定义二维表示上的一个或多个点。
10.根据权利要求9所述的方法,其中利用几何重映射将视网膜部分的二维表示的一个或多个坐标转换成视网膜部分的三维表示的等效坐标包括:利用几何重映射将二维表示的一个或多个点的坐标映射至二维表示的一个或多个等效点的坐标。
11.根据前述任意一项权利要求所述的方法,其中视网膜部分的几何测量包括对利用二维表示的第一和第二点的坐标定义在视网膜部分的二维表示上的视网膜部分的结构的距离测量。
12.根据权利要求11所述的方法,其中几何重映射被用来将视网膜部分的二维表示的第一和第二点转换成视网膜部分的三维表示的第一和第二等效点,而且三维表示的第一和第二等效点被用来确定距离几何测量。
13.根据权利要求1至10中任意一项所述的方法,其中视网膜部分的几何测量包括:定义在视网膜部分的二维表示上的视网膜部分的第一和第二结构之间的角度,其中利用二维表示的第一点和第二点的坐标来标记第一结构而且利用二维表示的第一点和第三点坐标来标记第二结构。
14.根据权利要求13所述的方法,其中几何重映射被用于将视网膜部分的二维表示的第一、第二和第三点转换成视网膜部分的三维表示的第一、第二和第三点等效点,而且三维表示的第一、第二和第三等效点被用于确定角度几何测量。
15.根据权利要求1至10中任意一项所述的方法,其中视网膜部分的几何测量包括利用二维表示的多个点的坐标定义在视网膜部分的二维表示上的视网膜部分的结构的面积。
16.根据权利要求15所述的方法,其中几何重映射被用于将视网膜部分的二维表示的多个点转换成视网膜部分的三维表示的多个等效点,而且三维表示的多个等效点被用于确定面积几何测量。
17.一种存储程序指令的计算机可读介质,在执行程序指令时执行根据权利要求1至16中任意一项所述的方法。
18.一种用于确定眼睛视网膜的几何测量的成像系统,包括:
成像装置,其获取眼睛视网膜的至少一部分的二维表示,
导出元件,其导出将视网膜部分的二维表示转换为视网膜部分的三维表示的几何重映射,
测量定义元件,其利用视网膜部分的二维表示的一个或多个坐标来在二维表示上定义视网膜的将被采取的几何测量,
重映射元件,其利用几何重映射将视网膜部分的二维表示的所述坐标或每个坐标转换成视网膜部分的三维表示的等效坐标,以及
测量确定元件,其利用视网膜部分的三维表示的所述等效坐标或每个等效坐标来确定眼睛视网膜的几何测量。
CN201410185236.7A 2013-05-02 2014-05-04 眼睛成像中的改进以及与眼睛成像相关的改进 Active CN104132649B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711373469.XA CN108371543B (zh) 2013-05-02 2014-05-04 眼睛成像中的改进以及与眼睛成像相关的改进

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GB1307990.0 2013-05-02
GBGB1307990.0A GB201307990D0 (en) 2013-05-02 2013-05-02 Improvements in and relating to imaging of the eye

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN201711373469.XA Division CN108371543B (zh) 2013-05-02 2014-05-04 眼睛成像中的改进以及与眼睛成像相关的改进

Publications (2)

Publication Number Publication Date
CN104132649A true CN104132649A (zh) 2014-11-05
CN104132649B CN104132649B (zh) 2018-07-10

Family

ID=48627235

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201711373469.XA Active CN108371543B (zh) 2013-05-02 2014-05-04 眼睛成像中的改进以及与眼睛成像相关的改进
CN201410185236.7A Active CN104132649B (zh) 2013-05-02 2014-05-04 眼睛成像中的改进以及与眼睛成像相关的改进

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201711373469.XA Active CN108371543B (zh) 2013-05-02 2014-05-04 眼睛成像中的改进以及与眼睛成像相关的改进

Country Status (8)

Country Link
US (3) US9649031B2 (zh)
EP (1) EP2800058B8 (zh)
JP (3) JP6469360B2 (zh)
CN (2) CN108371543B (zh)
DK (1) DK2800058T3 (zh)
ES (1) ES2713180T3 (zh)
GB (3) GB201307990D0 (zh)
HK (2) HK1201626A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107578449A (zh) * 2017-09-14 2018-01-12 南京方瑞科技有限公司 一种免校准互动白板投影机

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2509131B (en) * 2012-12-21 2017-05-17 Optos Plc Improvements in and relating to ophthalmoscopes
GB201407873D0 (en) * 2014-05-02 2014-06-18 Optos Plc Improvements in and relating to imaging of the eye
GB201501274D0 (en) * 2015-01-26 2015-03-11 Optos Plc Improvements relating to retinal treatment
JP7163921B2 (ja) * 2017-08-01 2022-11-01 株式会社ニコン 画像信号出力装置及び方法、画像データ変換装置及び方法、およびプログラム
WO2019074077A1 (ja) 2017-10-13 2019-04-18 株式会社ニコン 眼科システム、画像信号出力方法、画像信号出力装置、プログラム、及び三次元眼底画像生成方法
JP7014891B2 (ja) 2018-03-30 2022-02-01 株式会社ニコン 眼科装置及び眼科光学系
EP3889889A1 (en) * 2020-03-30 2021-10-06 Optos PLC Ocular image data processing

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5106183A (en) 1987-11-25 1992-04-21 Taunton Technologies, Inc. Topography measuring apparatus
US5159361A (en) 1989-03-09 1992-10-27 Par Technology Corporation Method and apparatus for obtaining the topography of an object
JPH0595903A (ja) * 1991-10-11 1993-04-20 Topcon Corp 眼科装置
JP2000210260A (ja) * 1999-01-22 2000-08-02 Mitsubishi Electric Corp 眼底画像システム
JP3625064B2 (ja) * 2002-06-28 2005-03-02 独立行政法人理化学研究所 眼底の立体表示および座標計測装置
EP1825804A1 (en) * 2004-10-27 2007-08-29 Kowa Company Ltd. Ophtalmological measuring apparatus
WO2008005611A2 (en) 2006-04-21 2008-01-10 Advanced Technology Materials, Inc. Systems and devices for mixing substances and methods of making same
GB0622325D0 (en) * 2006-11-09 2006-12-20 Optos Plc Improvements in or relating to retinal scanning
US7625088B2 (en) * 2007-02-22 2009-12-01 Kowa Company Ltd. Image processing apparatus
EP2232198B1 (en) 2008-01-08 2015-06-24 AMO WaveFront Sciences, LLC Systems and methods for measuring surface shape
US8194936B2 (en) * 2008-04-25 2012-06-05 University Of Iowa Research Foundation Optimal registration of multiple deformed images using a physical model of the imaging distortion
US8483446B2 (en) 2010-06-17 2013-07-09 Mississippi State University Method and system for estimating antler, horn, and pronghorn size of an animal
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
GB2509131B (en) * 2012-12-21 2017-05-17 Optos Plc Improvements in and relating to ophthalmoscopes

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
BARBARA DOBLER ETC.: "Precise modelling of the eye for proton therapy of intra-ocular tumours", 《PHYSICS IN MEDICINE AND BIOLOGY》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107578449A (zh) * 2017-09-14 2018-01-12 南京方瑞科技有限公司 一种免校准互动白板投影机

Also Published As

Publication number Publication date
CN104132649B (zh) 2018-07-10
JP2021003566A (ja) 2021-01-14
EP2800058B8 (en) 2019-03-27
JP6469360B2 (ja) 2019-02-13
CN108371543A (zh) 2018-08-07
GB2515178A (en) 2014-12-17
JP2014217755A (ja) 2014-11-20
CN108371543B (zh) 2020-04-17
DK2800058T3 (en) 2019-03-25
US10357150B2 (en) 2019-07-23
GB201307990D0 (en) 2013-06-12
US20190290119A1 (en) 2019-09-26
EP2800058B1 (en) 2019-02-06
GB201700162D0 (en) 2017-02-22
JP2019088823A (ja) 2019-06-13
GB2543445A (en) 2017-04-19
EP2800058A3 (en) 2014-12-31
GB201407679D0 (en) 2014-06-18
US9649031B2 (en) 2017-05-16
EP2800058A2 (en) 2014-11-05
ES2713180T3 (es) 2019-05-20
US20170215722A1 (en) 2017-08-03
HK1254982A1 (zh) 2019-08-02
US11259696B2 (en) 2022-03-01
JP7050131B2 (ja) 2022-04-07
US20140327877A1 (en) 2014-11-06
JP6763037B2 (ja) 2020-09-30
GB2515178B (en) 2017-02-22
HK1201626A1 (zh) 2015-09-04

Similar Documents

Publication Publication Date Title
CN104132649A (zh) 眼睛成像中的改进以及与眼睛成像相关的改进
US20200166333A1 (en) Hybrid light measurement method for measuring three-dimensional profile
US20130265571A1 (en) System and method for calibrating ultra wide-angle lenses
EP3907702A1 (en) Three-dimensional sensor system and three-dimensional data acquisition method
CN101577002B (zh) 应用于目标检测的鱼眼镜头成像系统标定方法
CN106500596B (zh) 结构光全景测量系统的测量方法
CN104142157A (zh) 一种标定方法、装置及设备
CN104007444A (zh) 一种基于中心投影的地面激光雷达反射强度图像的生成方法
CN109509143A (zh) 一种三维点云转换二维图像的方法
CN103996173A (zh) 一种基于变长轴椭圆拟合的鱼眼图像校正方法
CN103942802A (zh) 基于随机模板的结构光动态场景深度获取方法
JP6625030B2 (ja) 三次元計測システム及び三次元計測方法
US10542884B2 (en) Determining global shape parameters of a cornea using a two-fold fitting from measurement data to an aspheric toric surface model
JP2014134611A (ja) 幾何歪み補正装置、プロジェクタ装置、及び幾何歪み補正方法
CN103876706A (zh) 检眼镜中的改进以及与检眼镜相关的改进
Xing et al. A multi-system weighted fusion method to improve measurement accuracy of structured light 3D profilometry
JP2006194705A (ja) 液体下の表面形状測定方法及びそのシステム
CN108765301A (zh) 一种光学系统及其畸变校正方法及系统
CN115018923A (zh) 一种全景环带镜头参数的检测方法及系统
Lazarenko et al. The algorithm for transformation of images from omnidirectional cameras
KR101877548B1 (ko) 곡면 투사용 영상 생성 시스템 및 그 생성 방법
Muhr et al. A triangulation method for 3D-measurement of specular surfaces
CN205373637U (zh) 一种带电检测电气设备外绝缘配置参数测量装置
Rinaldo et al. Multi-view videogrammetry for the measurement of plate flexural vibration
Arfaoui et al. Calibration of aerial and panoramic camera using cross-diffractive optical elements

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
CP01 Change in the name or title of a patent holder

Address after: British Dan forest

Patentee after: Optos PLC.

Address before: British Dan forest

Patentee before: Photoelectron Corp.

CP01 Change in the name or title of a patent holder