CN113456018B - 眼睛图像数据处理 - Google Patents
眼睛图像数据处理 Download PDFInfo
- Publication number
- CN113456018B CN113456018B CN202110339136.5A CN202110339136A CN113456018B CN 113456018 B CN113456018 B CN 113456018B CN 202110339136 A CN202110339136 A CN 202110339136A CN 113456018 B CN113456018 B CN 113456018B
- Authority
- CN
- China
- Prior art keywords
- eye
- scan
- model
- values
- value
- 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
Links
- 238000012545 processing Methods 0.000 title claims abstract description 25
- 238000003384 imaging method Methods 0.000 claims abstract description 216
- 238000006243 chemical reaction Methods 0.000 claims abstract description 19
- 238000012014 optical coherence tomography Methods 0.000 claims description 147
- 210000001525 retina Anatomy 0.000 claims description 90
- 238000000034 method Methods 0.000 claims description 54
- 230000008859 change Effects 0.000 claims description 38
- 238000013507 mapping Methods 0.000 claims description 24
- 238000005094 computer simulation Methods 0.000 claims description 21
- 238000004088 simulation Methods 0.000 claims description 15
- 238000011156 evaluation Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 13
- 230000001902 propagating effect Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 8
- 238000013500 data storage Methods 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000000875 corresponding effect Effects 0.000 description 120
- 230000006870 function Effects 0.000 description 52
- 230000003287 optical effect Effects 0.000 description 17
- 230000002207 retinal effect Effects 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 7
- 230000015654 memory Effects 0.000 description 4
- 238000002583 angiography Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000003936 working memory Effects 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- MOFVSTNWEDAEEK-UHFFFAOYSA-M indocyanine green Chemical compound [Na+].[O-]S(=O)(=O)CCCCN1C2=CC=C3C=CC=CC3=C2C(C)(C)C1=CC=CC=CC=CC1=[N+](CCCCS([O-])(=O)=O)C2=CC=C(C=CC=C3)C3=C2C1(C)C MOFVSTNWEDAEEK-UHFFFAOYSA-M 0.000 description 2
- 229960004657 indocyanine green Drugs 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004256 retinal image Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000284466 Antarctothoa delta Species 0.000 description 1
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 208000010412 Glaucoma Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- GNBHRKFJIUUOQI-UHFFFAOYSA-N fluorescein Chemical compound O1C(=O)C2=CC=CC=C2C21C1=CC=C(O)C=C1OC1=CC(O)=CC=C21 GNBHRKFJIUUOQI-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 208000002780 macular degeneration Diseases 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000002747 omentum Anatomy 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/0016—Operational features thereof
- A61B3/0025—Operational features thereof characterised by electronic signal processing, e.g. eye models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/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/1005—Objective 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
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
- G06T7/55—Depth or shape recovery from multiple images
- G06T7/586—Depth or shape recovery from multiple images from multiple light sources, e.g. photometric stereo
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Ophthalmology & Optometry (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Quality & Reliability (AREA)
- Eye Examination Apparatus (AREA)
- Image Processing (AREA)
Abstract
本申请涉及眼睛图像数据处理。公开了一种用于通过下列操作来生成使眼睛成像装置的指示在眼睛中的眼睛成像装置获取眼睛的被成像区域的数字图像的扫描位置的扫描参数的值与用于计算在被成像区域中的指定眼睛特征之间的距离的转换因子的相应值相关的映射:模拟(S110)光线传播以使在扫描参数值的序列中的扫描参数的每个值与在模型眼睛中的相对应的位置相关;对于每个扫描参数值,计算(S120)在模型眼睛中的相对应的位置和在模型眼睛中的对应于在序列中的相邻值的位置之间的距离;以及使用(S130)计算出的距离来生成关于每个扫描参数值的转换因子的相应值。
Description
技术领域
本文的示例方面总体涉及数据处理的领域,且更具体地涉及用于处理眼睛成像数据的技术。
背景
在基于扫描仪的眼睛成像系统中,光束是经由扫描元件横越受验者的眼睛的一部分(例如视网膜)进行扫描的,并且从眼睛反射的返回光由成像系统针对扫描元件在扫描中覆盖的多个扫描角中的每一个扫描角进行采样。基于所采样的返回光来生成眼睛的该部分的被成像区域的数字图像。多种眼睛成像系统根据这个原理来操作,例如可用于几种视网膜成像方式(包括例如荧光素眼底血管造影(FFA)、吲哚菁绿(ICG)血管造影术和自发荧光眼底血管造影(FAF))的扫描激光检眼镜(SLO)以及光学相干断层扫描(OCT)成像系统(等等)。所获取的数字图像可以产生关于受验者的健康的有价值的信息。例如,视网膜的数字图像可以允许眼睛的疾病(例如黄斑变性和青光眼)以及全身性疾病(例如糖尿病、高血压和其他可能使受验者痛苦的心血管疾病)的并发症被检测、诊断和管理。
在对所获取的数字图像的分析中能够测量感兴趣的被成像特征的尺寸(例如血管直径或病变的尺寸)使得特征的尺寸随时间的任何变化可以被检测常常是很有价值的。作为示例,US 2010/0150415 A1公开了一种坐标重映射方法,该坐标重映射方法用于计算被成像特征的尺寸,其中对在二维宽视场视网膜图像中的特征的测量是根据在二维图像中的两个或更多个指定坐标进行定义的,这些坐标然后被映射到在视网膜的三维模型中的等效坐标,以及特征的尺寸是使用在三维模型中的等效坐标进行计算的。
概述
上面所述的计算眼睛特征尺寸的已知坐标重映射方法的问题在于重映射和后续距离计算操作对计算资源要求高。这些缺点使已知方法在例如视网膜或眼睛的其他部分的窄视场图像以高帧速率被获取的应用中变得特别成问题,使得许多相对复杂的重映射和后续距离计算需要被执行。发明人已经认识到,在这种类型的应用(以及其他应用)中,给每个数字图像提供用于将在数字图像中以多个图像像素(或用于量化数字图像中的指定距离的其他单位)表示的点之间的距离转换成在眼睛的被成像部分中的相对应的点之间可能以长度单位(例如,米制单位,诸如毫米)表示的物理距离的信息将是合乎需要的。
根据本文的第一示例方面,本发明人设计了一种用于处理由眼睛成像装置的控制器生成的眼睛成像数据以评估用于计算在眼睛的一部分的被成像区域中的指定眼睛特征之间的距离的转换因子的装置。该装置包括扫描参数获得模块,该扫描参数获得模块被配置为从眼睛成像数据获得眼睛成像装置的指示由眼睛成像装置执行以获取眼睛的该部分的被成像区域的数字图像的扫描的在眼睛的部分中的扫描位置的扫描参数的值。该装置还包括转换因子评估模块,该转换因子评估模块被配置为使用扫描参数的所获得的值以及在扫描参数的值和转换因子的相应值之间的映射来确定用于将在数字图像中的像素之间的距离转换为在眼睛的区域中的相对应的位置之间的距离的转换因子的相应值,其中转换因子的每个值指示针对扫描参数的相应值横越眼睛的部分的距离关于扫描参数的函数的模拟变化率,所述像素的相应像素值是在眼睛成像装置使用扫描参数的不同相应值对该区域的成像期间获取的。该装置还包括数据存储模块,该数据存储模块被配置为将数字图像与转换因子的所确定的值相关联地存储。
根据本文的第二示例方面,本发明人设计了一种处理由眼睛成像装置的控制器生成的眼睛成像数据以评估用于计算在眼睛的一部分的被成像区域中的指定眼睛特征之间的距离的转换因子的计算机实现的方法。该方法包括从眼睛成像数据获得眼睛成像装置的指示由眼睛成像装置执行以获取定义眼睛的该部分的被成像区域的图像的图像数据的扫描在眼睛中的扫描位置的扫描参数的值。该方法还包括使用扫描参数的所记录的值以及在扫描参数的值和转换因子的相应值之间的映射来确定用于将在图像中的像素之间的距离转换为在眼睛的区域中的相对应的位置之间的距离的转换因子的相应值,其中,转换因子的每个值指示针对扫描参数的相应值横越眼睛的部分的距离关于扫描参数的函数的模拟变化率,所述像素的相应像素值是在眼睛成像装置使用扫描参数的不同相应值对该区域的成像期间获取的。该方法还包括将数字图像与转换因子的所确定的值相关联地存储。
根据本文的第三示例方面,本发明人设计了一种生成使眼睛成像装置的扫描参数的值与用于计算在眼睛的部分的相应被成像区域中的指定眼睛特征之间的距离的转换因子的相应值相关的映射的计算机实现的方法,所述扫描参数的值指示在眼睛的部分中眼睛成像装置可操作来获取眼睛的该部分的相应被成像区域的数字图像的相应扫描位置。该方法包括运行对穿过眼睛成像装置的模型和眼睛的模型的光线传播的计算机模拟以生成使在扫描参数的值的序列中的扫描参数的每个值与在眼睛的模型的一部分中的相对应的位置相关的模拟结果,当眼睛成像装置的模型根据扫描参数值操作时穿过眼睛成像装置的模型传播的光线射到所述相对应的位置上。该方法还包括对于扫描参数的每个值计算在眼睛的模型的该部分中的相对应的位置和在眼睛的模型的该部分中对应于在值的序列中的相邻值的位置之间的距离。该方法还包括将每个计算出的距离除以在扫描参数的值的函数的相对应的值和扫描参数在序列中的相邻值的函数的值之间的差,以对于扫描参数的每个值生成指示横越眼睛的部分的距离关于扫描参数的函数的变化率的转换因子的相应值。
根据本文的第四示例方面,本发明人还设计了一种包括计算机可读指令的计算机程序,计算机可读指令当由计算机执行时使计算机执行根据上面阐述的第二示例方面和第三示例方面中的至少一个方面的方法。计算机程序可以存储在非暂时性计算机可读存储介质上或者由计算机可读信号携带。
附图简述
现在将仅通过非限制性示例的方式参考下面描述的附图来详细解释示例实施例。在所述附图中的不同附图中出现的相似的参考数字可以表示相同的或在功能上相似的元素,除非另有指示。
图1是根据本文的第一示例实施例的用于处理眼睛成像数据的装置的示意图。
图2是示出用于获取由图1的装置处理的成像数据的组合SLO-OCT成像装置的示例性形式的眼睛成像装置的实现方式的框图。
图3是示出第一示例实施例的装置在可编程信号处理硬件中的示例实现方式的框图。
图4是示出根据本文的第一示例实施例的处理眼睛成像数据的计算机实现方法的流程图。
图5以查找表的形式示出在第一示例实施例中使用的扫描参数的值和转换因子的值之间的映射的示例。
图6是示出根据示例实施例的生成扫描参数的值和转换因子的值之间的映射的计算机实现方法的流程图。
图7A是穿过眼睛成像装置的模型和眼睛的模型的光线传播的计算机模拟的示意图。
图7B示出根据示例实施例的基于图7A中的计算机模拟而生成的查找表的示例,该查找表可用于导出在扫描参数的值和转换因子的值之间的映射。
图8是穿过眼睛成像装置的模型和眼睛的模型的光线传播的示意图。
图9A是根据H-检流计镜倾斜角θ和V-检流计镜倾斜角的所计算的水平角转换因子的三维曲线图。
图9B是根据H-检流计镜倾斜角θ和V-检流计镜倾斜角的所计算的水平角转换因子的灰度级曲线图。
图10A是作为H-检流计镜倾斜角θ和V-检流计镜倾斜角的函数的所计算的垂直角转换因子的三维曲线图。
图10B是作为H-检流计镜倾斜角θ和V-检流计镜倾斜角的函数的所计算的垂直角转换因子的灰度图。
图11是示出根据示例实施例的使用转换因子的值来计算在眼睛的被成像区域中的指定眼睛特征之间的距离的计算机实现的方法的流程图。
图12示出了具有对应于在眼睛的被成像区域中的指定眼睛特征的两个指定像素的OCT B-扫描图像。
示例实施例的详细描述
第一示例实施例
图1是根据第一示例实施例的被配置成处理由眼睛成像装置的控制器生成的眼睛成像数据的装置100的示意图。装置100被配置为基于眼睛成像数据来评估用于计算在眼睛的一部分的由眼睛成像装置已成像(或将成像)的区域中的指定眼睛特征之间的距离的转换因子。如在本示例实施例中的,眼睛的该部分可以是眼睛的视网膜。然而,本文所述的技术可以可选地例如用于评估用于与眼睛的另外的部分的被成像区域(例如在眼的睛前段中的一个或多个结构)相关的距离计算的转换因子。
如图1所示,装置100包括扫描参数获得模块110、转换因子评估模块120和数据存储模块130。如在本示例实施例中的,装置100还可以包括映射生成器模块140和距离计算器模块150。下面更详细地描述装置100的示例实现方式及它的所示部件模块的功能。
图2示出了用于生成将由扫描参数获得模块110处理的眼睛成像数据的眼睛成像装置的示例实现方式。在图2的示例中,眼睛成像装置采取将扫描激光检眼镜(SLO)功能与光学相干断层扫描(OCT)功能组合的组合SLO-OCT眼睛成像装置200的形式。然而,应该注意的是,眼睛成像装置不需要是组合SLO-OCT眼睛成像装置,并且可以可选地采取采用一个或更多个扫描元件来对视网膜或眼睛的其他部分成像的任何其他成像系统(例如SLO成像系统或OCT成像系统)的形式。
组合SLO-OCT成像装置200包括具有OCT引擎210、聚焦机构212和V-检流计镜214的OCT成像部分。此外,SLO-OCT成像装置200包括具有SLO引擎220和多面镜230的SLO成像部分。SLO-OCT成像装置200还具有包括分色镜240、狭缝镜250、H-检流计镜260和椭球镜270的共享光学系统。OCT引擎210和SLO引擎220的操作由SLO-OCT成像装置200的控制器205控制,控制器205可以是被配置为与OCT引擎210和SLO引擎220(如图2所示)都通信的信号处理硬件的单个件或者分别控制OCT引擎210和SLO引擎220的这样的硬件的两个单独件。
在图2中,SLO引擎220包括被配置成(与控制器205结合地)生成受验者的眼睛280的视网膜的数字SLO图像的SLO分束器、SLO光源、SLO光电探测器。SLO光源发射SLO光,如在本示例中,其可以包括红光和绿光以及近红外光。分束器将来自SLO光源的SLO光引导到多面镜230,并将从眼睛280反射的SLO光(即,由眼睛280反射的源自入射SLO光的光)引导回到SLO光电探测器。
多面镜230将来自SLO分束器的SLO光反射到共享光学系统。特别是,如在图2中的示例中所示,当由驱动马达(未示出)驱动时,多面镜230通过绕所示轴C旋转来在Y方向上扫描SLO光。从眼睛280反射并经由共享光学系统被引导到SLO成像部分的光被SLO光电探测器探测到,并且来自SLO光电探测器的信号被控制器205处理以生成SLO图像,该SLO图像可以是例如覆盖视网膜的多达约80%的视网膜的超宽视场(UWF)图像。
OCT引擎210被用于生成眼睛280的视网膜的一部分的数字断层扫描图像,并且如在本示例中可以包括OCT光源、光学耦合器、参考反射镜和OCT光电探测器(未示出)。OCT光源发射被馈送到光学耦合器内并被分成参考光和信号光的低相干光。参考光被引导到参考反射镜,并被参考反射镜反射到光学耦合器。信号光被聚焦机构212引导到V-检流计镜上,V-检流计镜将信号光反射到共享光学系统。如在图2中的示例中所示,当被驱动电机驱动时,V-检流计镜214通过绕轴B旋转来在Y方向(即,垂直方向)上扫描信号光。V-检流计镜214还将被眼睛280反射并源自信号光的OCT光引导回到OCT引擎210。从眼睛280反射的OCT光在光学耦合器处与参考光叠加以生成干涉光。干涉光被OCT光电探测器探测到,并且探测结果被控制器205使用来生成数字OCT图像。
在图2中,分色镜240通过使来自SLO成像部分的多面镜230的SLO光被透射而将SLO光引导到狭缝镜250,并且通过使来自OCT成像部分的V-检流计镜214的OCT信号光被反射而将OCT信号光引导到狭缝镜250。
狭缝镜250朝着H-检流计镜260反射入射光。H-检流计镜260将来自狭缝镜250的光反射到椭球镜270。然后,如在图2的示例中所示,H-检流计镜260在由驱动电机(未示出)驱动时通过绕轴C旋转来在X方向(即,水平方向)上扫描发射光。
椭球镜270将来自H-检流计镜260的光引导到眼睛280的视网膜。由眼睛280反射的来自椭球镜270的光沿着与发射光相同的光路被引导到在共享光学系统中的分色镜240。分色镜240将从眼睛280反射的SLO光引导到SLO成像部分,并将从眼睛280反射的OCT光引导到OCT成像部分。
SLO-OCT成像装置200可在SLO成像模式、OCT模式和组合SLO-OCT成像模式中操作,在SLO成像模式中只有SLO成像系统操作来获取视网膜的数字SLO图像,在OCT模式中只有OCT成像系统操作来获取视网膜的数字OCT图像,以及在组合SLO-OCT成像模式中SLO和OCT成像系统都操作来在同一时间获取视网膜的数字SLO图像和视网膜的数字OCT图像。
在本示例实施例中,在眼睛280中的扫描位置以两个扫描参数值指示,扫描在该扫描位置处由SLO-OCT成像装置200已执行(或将执行)以生成OCT图像或SLO图像。每个扫描参数值提供视网膜的被扫描区域(例如被扫描区域的预定参考点,例如在被扫描区域的中心或在被扫描区域的预定边界处的点)沿着两个不同轴的相对应的轴所位于的位置的指示,视网膜上的任何点可以根据被扫描区域的坐标来被指定。然而应当注意,在其他示例实施例中可以使用单个扫描参数值来提供沿着横越视网膜延伸的线(直的或弯曲的)的扫描位置的指示,该扫描参数值提供视网膜的被扫描区域(例如被扫描区域的预定参考点,例如被扫描区域的中心或在被扫描区域的预定边界处的点)沿着该线所位于的位置的指示。还应当注意,扫描参数不需要是连续可变的,而是可以可选地是离散的,使得扫描参数值可以是离散的标识符,该标识符根据来自视网膜的表面被划分成的多个不同的预定区域当中的视网膜的预定区域来标识扫描的位置。
从前文将认识到,扫描参数可以采取许多不同的形式。作为示例,对于由SLO-OCT成像装置200执行的OCT扫描,扫描参数可以分别提供扫描元件214和260的倾斜角的指示,该倾斜角使OCT光束在旋转时横越眼睛280的视网膜被扫描。扫描元件214和260中的每一个的倾斜角是扫描元件围绕它的旋转轴相对于扫描元件的参考定向的旋转角。这些倾斜角中的每一个倾斜角的指示可以以几种不同的方式被提供,例如作为用于控制设置扫描元件的倾斜角的相对应的马达的控制信号的控制值,该值具有与相对应的扫描元件的倾斜角和因而在视网膜上的扫描位置的预定关系。这个关系可以通过例如模拟在SLO-OCT成像装置200的模型和眼睛280的模型中的光线传播来确定。
更特别地,参考图2,由SLO-OCT成像装置200获取的OCT C-扫描在眼睛280的视网膜上的扫描位置在本示例实施例中是通过H-检流计镜260和V-检流计镜214的倾斜角的值确定的,其中,如上面所提到的,当这些反射镜被定向成使得由此反射的OCT光束在(二维)扫描区域中的几何中心或另一参考点上时,θ表示H-检流计镜260的倾斜角,以及/>表示V-检流计镜214的倾斜角。为了解释的容易,在本示例实施例中,倾斜角θ和/>本身被认为是OCT扫描的扫描参数,但是从前文将理解,扫描参数可以是指示这些角度的其他参数(例如控制信号值)。
注意,由SLO-OCT成像装置200捕获的OCT图像不需要是C-扫描(包括在由SLO-OCT成像装置200执行的对视网膜的二维OCT扫描中生成的A-扫描的二维阵列),而是可以可选地是B-扫描(包括在由SLO-OCT成像装置200执行的对视网膜的一维OCT扫描中生成的A-扫描的一维阵列)。在任一情况下,在视网膜上的扫描位置可以由指示角θ和的值指示,H-检流计镜260和V-检流计镜214以角度θ和/>被设置,以便将OCT光引导到扫描区域的中心(或其他预定参考点)。因此,在视网膜上的任何地方的扫描位置可以通过两个扫描参数θ和/>的值指示。然而,SLO-OCT成像装置200可以由控制器205控制以执行扫描(B-扫描和/或C-扫描),扫描的位置可以仅沿着在视网膜上的预定线变化,并且沿着该线的扫描位置可以由单个扫描参数的值指示。
另一方面,当使用SLO-OCT成像装置200执行SLO扫描时,SLO光的被扫描到眼睛280中的位置取决于H-检流计镜260的倾斜角θ和多面镜230的倾斜角α。因此,对于由SLO-OCT成像装置200捕获的二维SLO图像,扫描参数可以包括指示H-检流计镜260的倾斜角的水平扫描参数和指示多面镜230的倾斜角α的垂直扫描参数。
尽管扫描参数指示扫描元件的倾斜角(该倾斜角由SLO-OCT成像装置200使用来执行扫描以便生成OCT或SLO图像),但是注意,扫描参数不限于这些示例,并且可以是其值可以提供在眼睛280的视网膜或其他部分上的扫描位置的指示的一个或多个其他参数中的任一个。例如,扫描参数可以可选地提供由被扫描光束在传播到眼睛280的视网膜时所覆盖(subtended)的角度的指示、以及眼睛的参考方向(例如,眼睛280的晶状体的光轴)或在SLO-OCT成像装置200中定义的参考方向。下面参考第二示例实施例更详细地描述这个种类的备选方案。
图3是装置100可以如何在可编程信号处理硬件中被实现的示意图。图3所示的信号处理装置300包括用于接收成像数据的接口模块310。信号处理装置300还包括用于控制装置100的处理器(CPU)320、工作存储器330(例如随机存取存储器)和存储包括计算机可读指令的计算机程序的指令存储库340,计算机可读指令在由处理器320执行时使处理器320执行装置100的处理操作。指令存储库340可以包括预加载有计算机可读指令的ROM(例如,电可擦除可编程只读存储器(EEPROM)或闪存的形式)。可选地,指令存储库340可以包括RAM或类似类型的存储器,并且计算机可读指令可以从计算机程序产品(例如计算机可读存储介质350(例如CD-ROM等)或者携带计算机可读指令的计算机可读信号360)被输入到指令存储库340。
在本示例实施例中,在图3中所示的包括处理器320、工作存储器330和指令存储库340的硬件部件的组合被配置为实现装置100的扫描参数获得模块110、转换因子评估模块120、映射生成器模块140和距离计算器模块150的功能,而工作存储器330或指令存储库340提供装置100的数据存储模块130。
图4是示出根据本文的第一示例实施例的计算机实现的方法的流程图,通过该方法,装置100处理由SLO-OCT成像装置200的控制器205生成的眼睛成像数据以评估用于计算在眼睛280的视网膜的被成像区域中的指定眼睛特征之间的距离的转换因子。
在图4的步骤S10中,扫描参数获得模块110从眼睛成像数据获得眼睛成像装置200的指示由眼睛成像装置200执行以获取眼睛280的该部分的被成像区域的数字图像的扫描在眼睛280中的扫描位置的扫描参数的值。在步骤S10中,如在本示例实施例中,扫描参数获得模块110可以通过接收作为由SLO-OCT成像装置200生成的眼睛成像数据的部分的第一扫描参数值θc和第二扫描参数值来获得这些值,第一扫描参数值θc和第二扫描参数值/>指示在眼睛280的视网膜上由SLO-OCT成像装置200执行OCT C-扫描的位置。如在本示例实施例中,眼睛成像数据还可以包括表征SLO-OCT成像装置200通过其来获取眼睛280的视网膜的被成像区域的数字断层扫描图像的图像获取过程的其他方面的数据,例如由对该区域的C-扫描覆盖的扫描参数值的范围Δθ和/>然而,应当注意的是,第一扫描参数值θc和第二扫描参数值/>可以替代地指示将由SLO-OCT成像装置200执行的C-扫描在视网膜上的位置;在这种情况下,由转换因子评估模块120计算的转换因子可以用于将在随后由SLO-OCT成像装置200获取的数字断层扫描图像中的像素之间的距离转换成在眼睛280的视网膜中的相对应的位置之间的距离。
眼睛280的视网膜的被成像部分的数字图像也可以在步骤S10中由装置100接收。眼睛成像数据和也可能在步骤S10中被接收的数字图像可以以本领域中的技术人员已知的任何合适的文件结构和/或文件格式被提供。例如,在一个示例性实现方式中,第一扫描参数值θc、第二扫描参数值和由对视网膜的该区域的C-扫描覆盖的扫描参数值的范围Δθ和/>(可以被包括在眼睛成像数据中的定义/表征C-扫描的信息的任何其他项)可以在包含视网膜的该区域的数字图像的数字图像文件的头部中或者在与数字图像文件分离并存储在数字图像文件的文件夹(或子文件夹)中的文件中被提供。
在图4的步骤S20中,转换因子评估模块120使用扫描参数的所获得的值以及在扫描参数的值和转换因子的相应值之间的映射125来确定用于将在数字断层图像中的像素之间的距离转换为在眼睛280的区域中的相对应的位置之间的距离的转换因子的相应值,其中转换因子的每个值指示针对扫描参数的相应值横越眼睛280的视网膜的距离关于扫描参数的函数的模拟变化率,所述像素的相应像素值是在对于扫描参数的不同相应值由眼睛成像装置200对该区域的成像期间获取的。如在本示例实施例中,扫描参数的函数可以是恒等函数,使得转换因子的每个值指示针对扫描参数的相应值横越眼睛280的视网膜的距离关于扫描参数的模拟变化率。然而,如下面所讨论的,函数不限于恒等函数,并且可以采取各种其他形式。
如在本示例实施例中,转换因子评估模块120可以使用在步骤S10中获得的第一扫描参数值θc和第二扫描参数值的值连同在扫描参数θ和/>的值与转换因子的相应值之间的映射125(其可以使用穿过SLO-OCT成像装置200和眼睛280的光线传播的计算机模拟来获得,如下所述)来针对第一扫描参数值θc和第二扫描参数值/>的接收到的值确定:(i)用于将在接收到的数字C-扫描图像中的像素之间的沿着在数字图像中的水平方向的距离转换为在眼睛280的视网膜的被成像区域中的相对应的位置之间的在视网膜上的相对应方向上的距离的水平转换因子CFh的相对应的值;以及(ii)用于将在接收到的数字C-扫描图像中的像素之间的沿着在数字图像中的垂直方向的距离转换为在眼睛280的视网膜的被成像区域中的相对应的位置之间的在视网膜上的相对应的方向上的距离的垂直转换因子CFv的相对应的值。对于扫描参数θ的给定值,水平转换因子CFh指示针对扫描参数θ的给定值横越眼睛280的视网膜的距离关于扫描参数θ的函数(在这个示例中,恒等函数)的模拟变化率。类似地,对于扫描参数/>的给定值,垂直转换因子CFv指示针对扫描参数/>的给定值横越眼睛280的视网膜的距离关于扫描参数/>的函数(在这个示例中,恒等函数)的模拟变化率。应当注意的是,一般来说,扫描参数θ的函数可以与扫描参数/>的函数相同或可以不与扫描参数的函数相同。
在由SLO-OCT成像装置200执行的扫描是OCT B-扫描的本示例实施例的变形中,转换因子评估模块120可以使用在步骤S10中获得的第一扫描参数值θc和第二扫描参数值的值连同在扫描参数θ和/>的值与转换因子的相应值之间的映射125(其可以使用穿过SLO-OCT成像装置200和眼睛280的光线传播的计算机模拟来获得,如下所述)来针对第一扫描参数值θc和第二扫描参数值/>的接收到的值确定用于将在数字B-扫描图像中的像素之间的在A-扫描排列方向上的距离(在A-扫描排列方向上由SLO-OCT成像装置200获取的A-扫描被排列以形成B-扫描)转换为在眼睛280的视网膜的被成像区域中的相对应的位置之间的在扫描方向上在视网膜上的距离的单个转换因子CF的相对应的值。
如在本示例实施例中,映射125可以采取如图5所示的查找表400的形式,该查找表400使扫描参数θ和的值与水平转换因子CFh和垂直转换因子CFv的相对应的值相关。然而,映射125可以以扫描参数θ和/>的函数的替代形式来被提供。此外,如在本示例实施例中,映射125可以使跨越由SLO-OCT成像装置200在视网膜的宽视场扫描中覆盖的扫描参数值的范围的扫描参数的值与水平转换因子CFh和垂直转换因子CFv的相应值(或者在上面提到的示例实施例的变形中的相应的单个转换因子CF)相关。宽视场扫描可以覆盖眼睛280的视网膜的表面的至少50%,并且更优选地覆盖视网膜的表面的至少60%、70%或80%。
在图4的步骤S30中,数据存储模块130例如使用指针来存储与在图4的步骤S20中确定的水平转换因子CFh和垂直转换因子CFv的值(或者在所描述的变形中确定的单个转换因子CF)相关联的数字图像,所述指针将存储所确定的转换因子的数据存储模块130的存储元件链接到数字图像文件。
上面参考图4描述的数据处理操作在它们的可应用性中并不限于与对OCT C-扫描的获取相关的由控制器205生成的眼睛成像数据,并且另外不仅可应用于(如上面所提到的)与OCT B-扫描的获取相关的由控制器205生成的眼睛成像数据,而且还可应用于SLO图像和表示视网膜或眼睛280的其他部分的一维或二维扫描的结果的其他眼睛图像。
图6示出了本示例实施例的(可选的)映射生成器模块140生成被转换因子评估模块120使用来确定转换因子的映射125的总体过程。例如,在装置100中对映射生成器模块140的包括提供了以下优点:允许映射125被更新以考虑到在SLO-OCT成像装置200中的光学元件或光学元件在SLO-OCT成像装置200中的布置的任何变化的原因,和/或使映射125适应于由SLO-OCT成像装置200成像的受验者的眼睛。然而,相同的过程可以在工厂中对SLO-OCT成像装置200的初始设置期间被使用来生成随后由转换因子评估模块120使用的映射125。
参考图6,在步骤S110中,映射生成器模块140通过运行穿过SLO-OCT成像装置200(或者更一般地,任何其他基于扫描仪的眼睛成像装置)的模型和眼睛280的模型的光线传播的计算机模拟来生成映射125以生成模拟结果。模拟结果使在扫描参数的值的序列中的扫描参数的每个值与在眼睛280的模型中的相对应的位置相关,当眼睛成像装置的模型根据扫描参数的值操作时穿过眼睛成像装置的模型传播的光线射到所述相对应的位置上。
如在本示例实施例中,计算机模拟可以使用将由眼睛成像装置成像的眼睛280的测量结果来对穿过眼睛280的模型和眼睛成像装置的模型的光线传播建模。这些测量结果可以提供关于眼睛的尺寸和/或形状的信息,并且可以使用眼睛成像技术(例如计算机断层扫描、超声波检查法和磁共振成像)来获得。该信息可以提高对眼睛建模的准确度,并因此使模拟结果变得更接近现实。
在图6的步骤S120中,映射生成器模块140针对扫描参数的每个值计算在眼睛280的模型中的相对应的位置(即,当眼睛成像装置的模型根据扫描参数的值操作时穿过眼睛成像装置的模型传播的光线射到的在眼睛280的模型中的位置)和在眼睛280的模型中的对应于在该值的序列中的相邻值的位置(即,当眼睛成像装置的模型根据与扫描参数的值相邻的扫描参数的第二值操作时穿过眼睛成像装置的模型传播的光线射到的在眼睛280的模型中的位置,例如,该相邻值是在序列中的下一个值或前一个值、或者在序列中的下下一个值或上上一个值)之间的距离。
在图6的步骤S130中,映射生成器模块140将每个计算出的距离除以在扫描参数的值的函数的相对应的值和扫描参数在序列中的相邻值的函数的值之间的差,以对于扫描参数的每个值生成指示横越眼睛280的该部分的距离关于扫描参数的函数的模拟变化率的转换因子的相应值。在函数是恒等函数的情况下,如在本示例实施例中,在图6的步骤S130中,映射生成器模块140可以将每个计算出的距离除以扫描参数的相对应的值和扫描参数在序列中的相邻值之间的差,以对于扫描参数的每个值生成指示横越眼睛280的部分的距离关于扫描参数的模拟变化率的转换因子的相应值。
因为图6的步骤S110、S120和S130被执行以生成由图3的步骤S20使用的映射125,步骤S110、S120和S130是在图4中的步骤S20由转换因子评估模块120执行之前由映射生成器模块140执行。
现在将参考图7A和图7B描述图6的过程的更详细的示例实现方式。在图6的步骤S110中,H-检流计镜260的倾斜角θ和V-检流计镜214的倾斜角可以分别被当作是第一扫描参数和第二扫描参数,并且计算机模拟可以被运行以对于θ和/>的多个不同的可能组合中的θ和/>的每对的值追踪光束(光线)604穿过SLO-OCT成像装置200的模型606到在眼睛280的视网膜的模型602上的位置上的传播,如在图7A中所示。因此,当SLO-OCT成像装置200的模型根据θ和/>的值对操作时,模拟结果将θ和/>的每对的值映射到光束入射到的在视网膜模型602上的相对应的位置。如在本示例实施例中,在θ和/>的每对的值与在视网膜模型602上的相对应的位置之间的映射可以存储在表中,如在图7B中所示。
如在本示例实施例中,扫描参数θ和的多对的值可以对应于跨越SLO-OCT成像装置200扫描实质上整个视网膜(例如视网膜的80%或更多)所需的θ和/>的相应(“最大”)范围的θ和/>的值的所有可能组合,其中在这些范围中的每一个范围内的相邻角值之间的差是相应扫描元件(H-检流计镜260或V-检流计镜214)可以被移动的最小增量。
然而在一些示例实施例中,用于光线传播的模拟运行的θ和的多对的值可以替代地是上面提到的所有可能组合的子集。例如,用于光线传播的模拟运行的θ和/或/>的值的组合可以跨越小于上面提到的“最大”范围的θ和/>的相应范围(或者这些角度之一的范围小于上面提到的相应“最大”范围)。
映射生成器模块140可以使上面提到的倾斜角值对(不管它们是否跨越“最大”范围或上面提到的更有限的范围中的至少一个)骤减以创建仅包含这些倾斜角值对的子集的稀疏选择。当存储空间减小时,以这种方式选择θ和值对是有利的。此外,模拟结果可以被内插以对于θ和/>的中间值(没有模拟针对中间值运行)确定光束604在模型视网膜602上的入射的位置。
图7A提供了在步骤S110中执行的光线传播的计算机模拟的示意图,并且示出了由SLO-OCT成像装置200的模型606发射的三个不同的光束,每个光束是在SLO-OCT成像装置200的模型606根据θ和的相应(不同)对的值操作时被发射的。如图7A所示,θ和/>的每对的值导致相应的模拟光束604在视网膜上的不同位置处入射在视网膜的模型602上。在图7A中,这些位置中的每个位置用相应的笛卡尔坐标(x,y,z)表示。
图7B示出了使θ和的每对的值与标识在视网膜模型602上的相对应的位置的坐标的相对应的集合相关的二维表600,当SLO-OCT成像装置200的模型606根据θ和/>的对的值操作时光束604在所述相对应的位置处入射在视网膜模型602上。如图7B所示,θ和/>的值的对可以布置在二维表600中,每行包含从/>到/>的/>的不同值,以及每列包含从θ1到θm的θ的不同值,其中m和n可以是相等的整数或可以是不相等的整数。在表600中的每个单元610存储针对对应于该单元610的θ和/>的值被计算的在视网膜模型602上的位置(x,y,z)。
如图6的步骤S120的实现方式,映射生成器模块140可以首先针对表600的每个单元610确定指示横越视网膜的距离关于的变化的值的模拟变化率的相应垂直角转换因子CFv。为了这个目的,映射生成器模块140可以针对表600中的每个单元610确定在由单元610中的坐标值(x,y,z)指示的在视网膜模型602上的位置和由具有相同的θ值但相邻的/>值的单元中的坐标值(x,y,z)指示的在视网膜模型602上的位置之间的距离。在此上下文中,词“相邻”指如与单元610的/>值相比的在该行中的下一最高或下一最低(或下一最高/最低)的/>值。例如,令/>表示/>的多对的值,其中i=1,……m,j=1,……n,每对的值具有在存储相对应的位置坐标(x,y,z)的表600中的相对应的单元610。因此,图6的步骤S120实现针对在表600中的所有/>值的/>和/>的相应(x,y,z)值之间的距离的计算,i=1,……m,j=1,……n。为了方便起见,在本示例中,对应于/>的(x,y,z)点在本文被表示为/>
作为示例,在图7B中,映射生成器模块140计算分别与和/>的/>值对应的在视网膜模型602上的点(x1,y1,z1)和(x2,y2,z2)之间的距离。如在本示例实施例中,可以使用在三维欧几里得空间中的勾股定理来计算该距离:
在本示例中,上面计算的距离与对应于值的对的单元610相关联地被存储。换句话说,对于表600中的每对值/>在针对角度值的对/>和角度值的对/>计算的相应(x,y,z)值之间的计算出的距离可以在表600中的对应于/>的单元610中被存储为/>
在图6的步骤S130中,映射生成器模块140还针对对应于扫描参数值的每个单元610计算/>的值。映射生成器模块140还将针对每个单元610的计算的距离除以对应的/>的值以生成关于该单元610的垂直角转换因子CFv,该垂直角转换因子CFv指示对于由单元610的θ和/>值指示的在视网膜上的位置横越视网膜的距离关于/>的模拟变化率。
在本示例实施例中,映射生成器模块140还确定关于每对值的水平角转换因子CFh,该水平角转换因子CFh指示横越视网膜的一部分的距离关于θ的模拟变化率。出于这个目的,在图6的步骤S120中,映射生成器模块140可以附加地执行关于在表中的所有值的/>和/>的相应(x,y,z)值之间的距离的计算。在图7B的示例中,计算出的距离对照每个条目/>被存储为/>
此外,在步骤S130中,映射生成器模块140还针对表600中的对应于相应角度值的对的每个单元计算值θi+1-θi。映射生成器模块140还将对应于/>的每个单元的所存储的计算出的距离/> 除以θi+1-θi的相对应的值,以便生成指示横越视网膜的部分的距离关于θ的模拟变化率的水平角转换因子CFh。
在本示例实施例中,转换因子指示横越眼睛的一部分的距离关于扫描参数(在步骤S20中是恒等函数的函数)的模拟变化率。因此,映射生成器模块140针对每对θ和值确定指示横越视网膜的一部分的距离关于θ的模拟变化率的水平角转换因子CFh,并且还针对每对θ和/>值确定指示横越视网膜的一部分的距离关于/>的模拟变化率的垂直角转换因子CFv。然而,在可选的示例实施例中,在步骤S20中的函数可以不是恒等函数,而是可以例如是使扫描参数的变化与另一参数(例如由眼睛成像装置扫描的光束的角度)的变化相关的函数。
在一些示例实施例中,转换因子可以指示横越眼睛的部分的距离关于在模拟中的光束的外部角的模拟变化率,该外部角由在被旋转以在模拟中执行扫描的眼睛成像装置的模型606中的检流计镜的倾斜角确定。该外部角可以是在步骤S110的计算机模拟中进入眼睛的模型602的模拟光束604和模型眼睛602的参考方向(例如模型眼睛602的晶状体的光轴)之间的角度。然而,外部角并不一定相对于模型眼睛602被定义,且例如,可以替代地,相对于眼睛成像装置的模型606的参考方向或平面被定义。
图8示出了在计算机模拟中穿过眼睛成像装置的模型606和眼睛的模型602的光线传播,其中外部角是模拟光束604相对于模型眼睛602中的眼睛晶状体的光轴的方位角α或仰角β。在图8中,晶状体的光轴被视为z轴。入射光束604的方位角α可以被定义为在眼睛602的晶状体的光轴和代表入射光束604的向量到水平(x-z)面上的投影801之间的角度。仰角可以被定义为在入射光束604和代表入射光束604的向量到水平面上的投影801之间的角度。参考图2的眼睛成像系统,方位角α可以是H-检流计镜的倾斜角θ的线性或者非线性函数。类似地,仰角β可以是V-检流计镜的倾斜角的线性或非线性函数。此外,α也可以是/>的函数,而β也可以是θ的函数。因此,入射光束604的方位角α和仰角β一起定义在模型中由眼睛成像系统606扫描的光束604射在眼睛602的位置。
作为示例,在第一示例实施例的变形中,映射生成器模块可以生成在倾斜角对和水平角转换因子CFh’之间的映射,水平角转换因子CFh’指示横越视网膜的部分的距离关于方位角α的模拟变化率。此外,该变形的映射生成器模块还可以生成在倾斜角对和垂直角转换因子CFv’之间的映射,垂直角转换因子CFv’指示横越视网膜的部分的距离关于仰角β(如上面所定义的)的模拟变化率。
该变形的映射生成器模块可以生成由转换因子评估模块120使用来以用类似于第一示例实施例的方式来确定转换因子的映射。更特别地,该变形的映射生成器模块可以执行与在图6的步骤S110中的示例实施例相同的过程,即选择关于的多对的值并运行计算机模拟来针对/>的每对的值计算在视网膜上的相对应的位置(x,y,z)。此外,如上面参考图7B所述,可以以表的形式提供在/>的值和相对应的位置(x,y,z)之间的映射。也可以以与在示例实施例中的相同的方式来执行在图6中的步骤S120的过程。
然而,本变形与示例实施例的不同体现因在步骤S130中执行的过程。在示例实施例(其中转换因子CFh和CFv的值被计算出)中,在视网膜上与相邻对值相对应的点之间的距离除以在相邻对之间的θ值的差(以生成水平角转换因子CFh)或除以在相邻对之间的/>值的差(以生成垂直角转换因子CFv)。然而,在本变形中,因为转换因子表示横越视网膜的距离关于方位角α或仰角β的变化率,所以该变形的映射生成器模块执行一个额外的步骤,即针对多对选择的/>值中的每对/>值导出关于(α,β)相对应的一对值。为了从一对/>值导出一对相对应的(α,β)值,该变形的映射生成器模块可以首先执行计算机模拟以确定对应于该对/>值的单位向量(L,M,N),该单位向量(在笛卡尔坐标中)表示当眼睛成像系统的水平检流计镜和垂直检流计镜的倾斜角根据/>值被设置时由眼睛成像系统扫描到模型眼睛602中的模拟入射光束604。该变形的映射生成器模块还可以将每个单位向量(L,M,N)从笛卡尔坐标转换为球面坐标以获得表示对应于单位向量的方位角θ和仰角/>的/>值。
在获得对应于在多对选择的值中的每对/>值的(α,β)值之后,该变形的映射生成器模块可以针对每对/>值计算/> 对的α值和在在多个选择的/>对具有的相同/>值但下一最高或下一最低θ值的/>对的α值之间的差。
作为示例,回来参考图7B中的表600,其中每个单元610对应于一对值,该变形的映射生成器可以针对表600中的每个/>对计算相对应的一对(αi,βj)值,其中i=1,……m,j=1,……n。在修改的步骤S130中,该变形的映射生成器模块还可以针对表600中的对应于相应角度值对/>的每个单元610计算值αi+1-αi,其中αi+1是对应于/>的方位角α。
该变形的映射生成器模块还可以将在修改的步骤S120中针对每个对(i=1,……m,j=1,……n)计算的距离/>(即,在/>和/>的相应(x,y,z)值之间的距离)除以对应的αi+1-αi的值,以便生成指示横越视网膜的部分的距离关于方位角α的模拟变化率的水平角转换因子CFh’。通过针对在多个选择的/>对中的每对值计算CFh’,可以针对倾斜角的期望范围建立在/>和CFh’之间的映射。
此外,在修改的步骤S130中,该变形的映射生成器模块还可以针对表600中的对应于相应角度值对的每个单元计算值βj+1-βj,其中βi+1是对应于/>的仰角β。该变形的映射生成器模块还可以将针对每个/>对计算的距离/>(即,如在修改的步骤S120中获得的在/>和/>的相应(x,y,z)值之间的距离)除以对应的βj+1-βj的值,以便生成指示横越视网膜的部分的距离关于仰角β的模拟变化率的垂直转换因子CFv’。通过对在选择的多个对/>中的所有/>对执行该计算,可以建立在/> 和相对应的垂直转换因子CFv’之间的映射。
图9A示出了根据H-检流计镜倾斜角θ和V-检流计镜倾斜角计算出的水平角转换因子CFh’的三维曲线图。图9B是图9A所示的计算结果的灰度级表示。
图10A示出了根据H-检流计镜倾斜角θ和V-检流计镜倾斜角计算出的垂直角转换因子CFv’的三维曲线图。图10B是图10A所示的计算结果的灰度级表示。/>
应当理解的是,尽管图7A和图7B的示例采用两个扫描参数并因此使用映射的两个相对应的集合,但是装置100可以类似地仅利用一个扫描参数来工作。例如,对于通过仅改变V-检流计镜倾斜角获得的OCT B-扫描,映射生成器模块140仅需要计算针对/>的垂直角转换因子。然后,该垂直角转换因子可以被用于将在所生成的OCT数字图像中的像素之间的在第一方向(例如垂直方向)上的距离转换成在视网膜的表面上的相对应的位置之间的在相对应的方向(例如垂直方向)上的距离。
图11是示出一个过程的流程图,装置100的(可选的)距离计算器模块150通过该过程使用在图4的步骤S20中确定的转换因子的值将在步骤S10中获取的数字图像中的两个像素之间的距离转换成在数字图像中示出的眼睛280的图像区域中的相对应的眼睛特征之间的物理距离来确定在眼睛280的被成像区域中的指定眼睛特征之间的距离。
在图11的步骤S40中,距离计算器模块150接收对在数字图像中的像素中的第一像素的指定,该第一像素表示在视网膜的被成像区域(或眼睛280的另一成像部分)中的第一眼睛特征。距离计算器模块150还接收对在数字图像中的像素中的第二像素的指定,该第二像素表示在被成像区域中的第二眼睛特征。应当注意的是,第一眼睛特征和第二眼睛特征中的每一个眼睛特征可以与视网膜的在视网膜的被成像区域中的任何部分相关。
如在本示例实施例中,对第一像素和第二像素的指定可以由观看在屏幕上的数字图像的用户进行,例如,该用户可以使用输入设备(例如计算机键盘或鼠标)来使用在屏幕上的光标等来在数字图像中指定感兴趣的特征。然而,该指定可以可选地由图像处理软件(例如运行自动指定在数字图像中的感兴趣点的模式识别算法的软件)执行。然后,在数字图像中的指定像素的坐标被用于距离计算。
更特别地,在图11的步骤S50中,距离计算器模块150使用在图4的步骤S20中确定的转换因子的值来将在数字图像中的第一像素和第二像素之间的距离转换成在眼睛280的被成像区域中的第一眼睛特征和第二眼睛特征之间的物理距离(即,以长度单位(例如,米制单位(例如毫米)计的距离)。
作为示例,图12示出了视网膜的OCT B-扫描图像的形式的数字图像132,其包括如上所述由用户指定的第一像素136-1和第二像素136-2。在指定像素136-1和136-2之间的距离138由距离计算器模块150转换成在视网膜中的相对应的眼睛特征之间的距离,其可以在屏幕上被显示给用户。
例如,在OCT图像是通过在预定范围内改变V-检流计镜214的倾斜角而生成的OCT B-扫描的情况下,B-扫描具有与V-检流计镜214的扫描方向相对应的x维以及与眼睛280的深度方向相对应的y维。在这种情况下,通过使用关于图7B中的示例而解释的垂直角转换因子CFv或CFv’,可以将在OCT B-扫描图像中的任两个像素之间的沿着x维的距离转换成在视网膜上的两个相对应的特征之间的物理距离。
更具体地,距离计算器模块150可以首先确定在第一指定像素和第二指定像素之间的像素距离。距离计算器150还可以确定与所确定的像素距离对应的V-检流计镜214的倾斜角的变化/>更具体地,角/>是V-检流计镜214为了生成关于在第一指定像素和第二指定像素之间的确定的像素距离的像素数据而已旋转(或将旋转)的角度。例如,可以基于OCT光电探测器的采样率和V-检流计镜的角速度来计算这个角变化/>
当垂直转换因子指示距离关于V-检流计镜的倾斜角的变化率时,距离计算器模块150因此可以通过将可应用于正被处理的OCT B-扫描图像的垂直角转换因子CFv的值乘以所确定的角变化/>来计算在对应于两个指定像素的视网膜特征之间的真实(物理)距离。如关于图4所解释的,垂直转换因子CFv可以在图4的步骤S20中由转换因子评估模块120基于在步骤S10中由扫描参数获得模块110接收的/>的值(即/>)来确定。
虽然上面的示例涉及通过旋转V-检流计镜214而生成的OCT B-扫描图像,但是应该明显的是,OCT B-扫描图像可以通过改变H-检流计镜260的倾斜角θ来类似地生成,在这种情况下,OCT B-扫描图像将具有与H-检流计镜260的扫描方向对应的x维。在这样的情形中,关于图7B中的示例解释的水平角转换因子CFh可以被用于将在OCT B-扫描图像中的任意两个像素之间的沿着x维的距离转换成在视网膜上的两个相对应的特征之间的物理距离。
在上面的示例中,V-检流计镜214和H-检流计镜260中的仅仅一者被旋转以获取B-扫描。在SLO-OCT成像系统200通过改变H-检流计镜260和V-检流计镜214两者的倾斜角来捕获OCT C-扫描的情况下,OCT C-扫描图像可以具有与H-检流计镜260的扫描方向对应的第一维(例如x维)、与V-检流计镜214的扫描方向对应的第二维(例如y维)和与眼睛280的深度方向对应的第三方向(例如z维)。对于仅具有第一维和第二维(即,x维和y维)的OCT C-扫描的OCT图像切片,距离计算器模块150可以将在C-扫描图像切片中的两个指定像素之间的距离转换成在眼睛280中的相对应的视网膜特征之间的物理距离。特别是,如果第一指定像素和第二指定像素在x维和y维上都相对于彼此偏移,则距离计算器模块150可以使用水平角转换因子CFh和垂直角转换因子CFv来计算在与两个指定像素对应的相应视网膜特征之间的物理距离。
作为示例,如果第一指定像素和第二指定像素沿着OCT图像切片的x维间隔开a个像素以及沿着OCT图像切片的y维间隔开b个像素,则距离计算器模块150可以确定对应于距离a的H-检流计镜260的角δθ,并且还确定对应于距离b的V-检流计镜214的角角δθ是H-检流计镜260为了生成关于在OCT图像切片的x维上的像素距离a的像素数据而需要旋转的角度。类似地,角/>是V-检流计镜为了生成关于在OCT图像切片的y维上的像素距离b的像素数据而需要旋转的角度。例如,δθ和/>的值可以基于OCT光电探测器的采样率和H-检流计镜和V-检流计镜的已知角速度来被计算,或者它们可以在扫描参数获得模块110接收上面提到的范围Δθ和/>作为成像数据的部分的情况下分别被计算为a/A·Δθ和/>其中数字OCT图像切片包括具有A×B个像素的阵列。
在获得δθ和之后,距离计算器模块150可以将δθ乘以在步骤S20中获得的水平角转换因子CFh以计算在对应于两个指定像素的视网膜特征之间的沿着第一方向(例如水平方向)的第一物理距离。距离计算器模块150还可以将/>乘以也在步骤S20中获得的垂直角转换因子CFv以计算在对应于两个指定像素的视网膜特征之间的沿着第二方向(例如垂直方向)的第二物理距离。最后,距离计算器模块150可以通过使用第一物理距离和第二物理距离来确定在第一视网膜特征(对应于第一指定像素)和第二视网膜特征(对应于第二指定像素)之间的距离。例如,在对应于第一物理距离的第一方向实质上正交于对应于第二物理距离的第二方向的情况下,距离计算器模块150可以通过将勾股定理应用于计算出的第一和第二物理距离来计算在两个视网膜特征之间的距离。这个计算出的距离还可以由装置100输出。例如,装置100可以控制计算机屏幕或其他视觉显示单元以向用户显示计算出的距离。
应当注意的是,当转换因子指示距离关于除恒等函数之外的扫描参数的函数的变化率时,则对应于(在两个指定像素之间)所确定的像素距离的扫描参数的变化首先被转换成扫描参数的函数的值的相对应的变化。转换因子(其指示距离关于扫描参数的函数的变化率)然后可以与扫描参数的函数的值的相对应的变化相乘,以便获得在眼睛280中的对应于在数字图像中的指定像素之间的距离的物理距离。
根据本文的示例实施例,在前文中描述了如下面在E1至E7中阐述的装置以及如下面在E8至E14中阐述的计算机实现的方法以及如下面在E15中阐述的计算机程序。
E1.一种用于处理由眼睛成像装置200的控制器205生成的眼睛成像数据以评估用于计算在眼睛(280)的一部分的被成像区域中的指定眼特征之间的距离的转换因子的装置(100),该装置包括:
扫描参数获得模块110,其被配置为从眼睛成像数据获得眼睛成像装置(200)的指示由眼睛成像装置(200)执行以获取眼睛(280)的部分的被成像区域的数字图像的扫描在眼睛(280)中的扫描位置的扫描参数的值;
转换因子评估模块(120),其被配置为使用扫描参数的所获得的值以及在扫描参数/>的值和转换因子(CFh,CFv)的相应值之间的映射(125)来确定用于将在数字图像中的像素之间的距离转换为在眼睛(280)的区域中的相对应的位置之间的距离的转换因子(CFh,CFv)的相应值,其中转换因子的每个值指示针对扫描参数/>的相应值横越眼睛(280)的部分的距离关于扫描参数/>的函数的模拟变化率,像素的相应像素值是在由眼睛成像装置(200)使用扫描参数/>的不同相应值对该区域的成像期间获取的;以及
数据存储模块(130),其被配置为将数字图像(132)与转换因子(CFh,CFv)的所确定的值相关联地存储。
E2.根据E1所述的装置(100),还包括被配置为通过下列操作来生成映射的映射生成器模块(140):
运行穿过眼睛成像装置(200)的模型和眼睛(280)的模型的光线传播的计算机模拟以生成使在扫描参数的值的序列中的扫描参数的每个值与在眼睛(280)的模型的一部分中的相对应的位置相关的模拟结果,当眼睛成像装置(200)的模型根据扫描参数的值操作时穿过眼睛成像装置(200)的模型传播的光线射到该相对应的位置上;
对于扫描参数的每个值,计算在眼睛(280)的模型的该部分中的相对应的位置和在眼睛(280)的模型的该部分中的对应于在值的序列中的相邻值的位置之间的距离;以及
将每个计算出的距离除以在扫描参数的值的函数的相对应的值和扫描参数/>在序列中的相邻值的函数的值之间的差,以对于扫描参数/>的每个值生成指示横越眼睛(280)的该部分的距离关于扫描参数/>的函数的模拟变化率的转换因子(CFh,CFv)的相应值。/>
E3.根据E2所述的装置(100),其中,映射生成器模块(140)被配置为在计算机模拟中使用眼睛(280)的生物测量结果来对穿过眼睛(280)的模型和眼睛成像装置(200)的模型的光线传播建模。
E4.根据E1至E3中的任一项所述的装置(100),还包括距离计算器模块(150),该距离计算器模块被配置为通过下列操作来确定在眼睛(280)的被成像区域中的指定眼睛特征之间的距离:
接收对在数字图像(132)中像素中表示在眼睛(280)的区域中的第一眼睛特征的第一像素(136-1)的指定以及对在数字图像(132)中像素中表示在眼睛(280)的区域中的第二眼睛特征的第二像素(136-2)的指定;以及
使用转换因子(CFh,CFv)的所确定的值来将在数字图像(132)中的第一像素(136-1)和第二像素(136-2)之间的距离(138)转换成在眼睛(280)的区域中的第一眼睛特征和第二眼睛特征之间的距离。
E5.根据E1至E4中的任一项所述的装置(100),其中,眼睛(280)的部分包括眼睛(280)的视网膜。
E6.根据E5所述的装置(100),其中,所述映射(125)使跨越由眼睛成像装置(200)在视网膜的宽视场扫描中覆盖的扫描参数值的范围的扫描参数的值与转换因子(CFh,CFv)的相应值相关。
E7.根据E1至E6中的任一项所述的装置(100),其中
装置(100)被配置为处理由光学相干断层扫描OCT成像装置的控制器(205)生成的眼科OCT成像数据作为眼睛成像数据,以及
扫描参数获得模块(110)被配置为从OCT成像数据获得OCT成像装置的指示由OCT成像装置执行以获取眼睛(280)的该部分的被成像区域的数字OCT图像的OCT扫描在眼睛(280)中的OCT扫描位置的扫描参数的值。
E8.一种处理由眼睛成像装置(200)的控制器(205)生成的眼睛成像数据以评估用于计算在眼睛(280)的一部分的被成像区域中的指定眼睛特征之间的距离的转换因子(CFh,CFv)的计算机实现的方法,该方法包括:
从眼睛成像数据获得(S10)眼睛成像装置(200)的指示由眼睛成像装置(200)执行以获取眼睛(280)的该部分的被成像区域的数字图像的扫描在眼睛(280)中的扫描位置的扫描参数的值;
使用(S20)扫描参数的所获得的值以及在扫描参数/>的值和转换因子的相应值之间的映射(125)来确定用于将在数字图像(132)中的像素之间的距离(138)转换为在眼睛(280)的区域中的相对应的位置之间的距离的转换因子(CFh,CFv)的相应值,其中转换因子(CFh,CFv)的每个值指示针对扫描参数/>的相应值横越眼睛(280)的该部分的距离关于扫描参数/>的函数的模拟变化率,像素的相应像素值是在由眼睛成像装置(200)使用扫描参数/>的不同相应值对该区域的成像期间获取的;以及
将数字图像与转换因子(CFh,CFv)的所确定的值相关联地存储(S30)。
E9.根据E8的计算机实现的方法,还包括通过下列操作来生成映射:
运行(S110)穿过眼睛成像装置(200)的模型和眼睛(280)的模型的光线传播的计算机模拟以生成使在扫描参数的值的序列中的扫描参数/>的每个值与在眼睛(280)的模型的一部分中的相对应的位置相关的模拟结果,当眼睛成像装置(200)的模型根据扫描参数/> 的值操作时穿过眼睛成像装置(200)的模型传播的光线射到该相对应的位置上;
对于扫描参数的每个值,计算(S120)在眼睛(280)的模型的该部分中的相对应的位置和在眼睛(280)的模型的该部分中的对应于在值的序列中的相邻值的位置之间的距离;以及
将每个计算出的距离除以(S130)在扫描参数的相对应的值的函数的值和扫描参数/>在序列中的相邻值的函数的值之间的差,以对于扫描参数/>的每个值生成指示横越眼睛(280)的该部分的距离关于扫描参数/>的函数的模拟变化率的转换因子(CFh,CFv)的相应值。
E10.根据E9的计算机实现的方法,其中,计算机模拟使用眼睛(280)的生物测量结果来对穿过眼睛(280)的模型和眼睛成像装置(200)的模型的光线传播建模。
E11.根据E8至E10中的任一项所述的计算机实现的方法,还包括通过下列操作来确定在被成像区域中的指定眼睛特征之间的距离:
接收(S40)对在数字图像(132)中的像素中的表示在眼睛(280)的区域中的第一眼睛特征的第一像素(136-1)的指定以及对在数字图像(132)中的像素中的表示在眼睛(280)的区域中的第二眼睛特征的第二像素(136-2)的指定;以及
使用转换因子(CFh,CFv)的所确定的值来将在数字图像(132)中的第一像素(136-1)和第二像素(136-2)之间的距离(138)转换成在眼睛(280)的区域中的第一眼睛特征和第二眼睛特征之间的距离。
E12.根据E8至E11中的任一项所述的计算机实现的方法,其中,眼睛的该部分包括眼睛(280)的视网膜,以及所述映射将跨越由眼睛成像装置(200)在视网膜的宽视场扫描中覆盖的扫描参数值的范围的扫描参数/>的值与转换因子(CFh,CFv)的相应值相关。
E13.根据E8至E12中的任一项所述的计算机实现的方法,其中,由光学相干断层扫描OCT成像装置的控制器(205)生成的眼科OCT成像数据被作为眼睛成像数据处理,以及
所述获得(S10)包括从OCT成像数据获得OCT成像装置的指示由OCT成像装置执行以获取眼睛(280)的该部分的被成像区域的数字OCT图像的OCT扫描在眼睛(280)中的OCT扫描位置的扫描参数的值。
E14.一种生成使眼睛成像装置(200)的指示在眼睛(280)的一部分中的相应扫描位置的扫描参数的值与用于计算在眼睛(280)的部分的相应被成像区域中的指定眼睛特征之间的距离的转换因子(CFh,CFv)的相应值相关的映射的计算机实现的方法,眼睛成像装置(200)能够操作来在所述相应扫描位置处获取眼睛(280)的该部分的相应被成像区域的数字图像,该方法包括:
运行(S110)穿过眼睛成像装置(200)的模型和眼睛(280)的模型的光线传播的计算机模拟,以生成使在扫描参数的值的序列中的扫描参数/>的每个值与在眼睛(280)的模型的一部分中的相对应的位置相关的模拟结果,当眼睛成像装置(200)的模型根据扫描参数/>的值操作时穿过眼睛成像装置(200)的模型传播的光线射到该相对应的位置上;
对于扫描参数的每个值,计算(S120)在眼睛(280)的模型的该部分中的相对应的位置和在眼睛(280)的模型的该部分中的对应于在值的序列中的相邻值的位置之间的距离;以及
使每个计算出的距离除以在扫描参数的相对应的值的函数的值和扫描参数在序列中的相邻值的函数的值之间的差,以对于扫描参数/>的每个值生成指示横越眼睛(280)的部分的距离关于扫描参数/>的函数的变化率的转换因子(CFh,CFv)的相应值。
E15.一种包括计算机可读指令的计算机程序(390),所述计算机可读指令当由计算机处理器(320)执行时使计算机处理器(320)执行根据E8至E14中的至少一项所述的方法。
本文描述的示例方面避免了与用于在眼睛图像处理中计算眼睛特征尺寸的常规技术相关的限制,特别是避免了来源于计算机技术的限制。例如,使用已知的坐标重映射方法来计算眼睛特征尺寸可能涉及缓慢的、资源密集的计算,缓慢的、资源密集的计算在例如视网膜或眼睛的其他部分的窄视场图像以高帧速率被获取的应用中是特别成问题的,使得许多相对复杂的重映射和后续距离计算需要被执行。另一方面,相对于常规技术,通过对所确定的转换因子值的使用,本文所述的技术使在视网膜图像中的指定眼睛特征之间的距离明显更快地且以在计算上更有效的方式使用更少的计算机处理资源和存储器资源被实现。而且,借助于来源于计算机技术的在本文描述的示例方面的前述能力,在本文描述的示例方面改进了计算机和计算机处理/功能,并且还改进了至少图像处理、OCT和数据处理的领域。
在前述描述中,参考几个示例实施例描述了示例方面。因此,说明书应被视为说明性的而不是限制性的。类似地,在附图中示出的突出示例实施例的功能和优点的附图仅仅是为了示例目的而被呈现的。示例实施例的体系结构是足够灵活的和可配置的,使得它可以以除了在附图中所示的方式以外的方式被利用。
虽然上面描述了本发明的各种示例实施例,但是应当理解,它们是作为示例而不是限制被呈现。对在相关领域中的技术人员将明显的是,可以在形式和细节上对其做出各种改变。因此,本发明不应受上述示例实施例中的任一个实施例的限制,而应仅根据随附的权利要求及其等同物来被限定。
此外,摘要的目的通常是使专利局和公众、以及尤其是不熟悉专利或法律术语或措辞的本领域中的科学家、工程师和从业人员能够根据粗略的检查快速确定本申请的技术公开的性质和本质。摘要并不意欲以任何方式关于在本文呈现的示例实施例的范围进行限制。还应该理解的是,在权利要求中叙述的任何过程不需要以所呈现的顺序来被执行。
虽然本说明书包含很多特定实施例细节,但这些不应被解译为对任何发明的范围或可能被要求保护的内容的限制,而是作为特定于本文描述的特定实施例的特征的描述。在单独的实施例的背景下在本说明书中所描述的某些特征也可在单个实施例中结合实现。相反地,也可以在多个实施例中单独地或以任何合适的子组合实现在单一实施例的背景下描述的各种特征。此外,尽管特征可能上文被描述为在特定组合中起作用并甚至起初如此要求保护的,但是来自所要求保护的组合的一个或更多个特征可以在一些情况下从组合中删除,且所要求保护的组合可以针对子组合或子组合的变形。
现在已经描述了一些说明性实施例和实施例,显然,前面的实施例是说明性的而不是限制性的,已经通过示例的方式被给出。
本文描述的装置可以体现在其他特定的形式中而不偏离其特征。前述实施例是说明性的,而不是对所描述的系统和方法的限制。因此,本文所描述的装置的范围由所附权利要求而不是前述描述指示,以及在权利要求的等价物的含义和范围内的变化被包括在其中。
Claims (15)
1.一种用于处理由眼睛成像装置(200)的控制器(205)生成的眼睛成像数据以评估用于计算在眼睛(280)的部分的被成像区域中的指定眼睛特征之间的距离的转换因子的装置(100),所述装置包括:
扫描参数获得模块(110),所述扫描参数获得模块被配置为从所述眼睛成像数据获得所述眼睛成像装置(200)的扫描参数的值,所述扫描参数/>指示由所述眼睛成像装置(200)执行来获取所述眼睛(280)的所述部分的所述被成像区域的数字图像的扫描在所述眼睛(280)中的扫描位置;
转换因子评估模块(120),所述转换因子评估模块被配置为使用所述扫描参数的所获得的值以及在所述扫描参数/>的值和所述转换因子(CFh,CFv)的相应值之间的映射(125)来确定用于将在所述数字图像中的像素之间的距离转换为在所述眼睛(280)的所述被成像区域中的相对应的位置之间的距离的所述转换因子(CFh,CFv)的相应值,其中所述转换因子的每个值指示针对所述扫描参数/>的相应值横越所述眼睛(280)的所述部分的距离关于所述扫描参数/>的函数的模拟变化率,像素的相应像素值是在由所述眼睛成像装置(200)使用所述扫描参数/>的不同相应值对所述成像区域的成像期间获取的;
数据存储模块(130),所述数据存储模块被配置为将所述数字图像(132)与所述转换因子(CFh,CFv)的所确定的值相关联地存储;以及
映射生成器模块(140),所述映射生成器模块被配置为通过以下操作来生成所述映射:
运行穿过所述眼睛成像装置(200)的模型和所述眼睛(280)的模型的光线传播的计算机模拟以生成使模拟结果,所述模拟结果将在所述扫描参数的值的序列中的扫描参数的每个值与在所述眼睛(280)的模型的部分中的相对应的位置相关,当所述眼睛成像装置(200)的模型根据所述扫描参数的值操作时,穿过所述眼睛成像装置(200)的模型传播的光线射到所述相对应的位置上;
对于所述扫描参数的值中的每个值,计算在所述眼睛(280)的模型的所述部分中的所述相对应的位置和在所述眼睛(280)的模型的所述部分中的对应于在值的所述序列中的相邻值的位置之间的距离;以及
将计算出的距离中的每个距离除以在所述扫描参数的值的函数的相对应的值和所述扫描参数/>在所述序列中的所述相邻值的函数的值之间的差,以对于所述扫描参数/>的值中的每个值生成所述转换因子(CFh,CFv)的指示横越所述眼睛(280)的所述部分的距离关于所述扫描参数/>的函数的模拟变化率的相应值。
2.根据权利要求1所述的装置(100),其中,所述映射生成器模块(140)被配置为在所述计算机模拟中使用对所述眼睛(280)的生物测量结果来对穿过所述眼睛(280)的模型和所述眼睛成像装置(200)的模型的光线传播建模。
3.根据权利要求1至2中的任一项所述的装置(100),还包括距离计算器模块(150),所述距离计算器模块被配置为通过下列操作来确定在所述眼睛(280)的所述被成像区域中的所述指定眼睛特征之间的距离:
接收对在所述数字图像(132)中像素中表示在所述眼睛(280)的所述被成像区域中的第一眼睛特征的第一像素(136-1)的指定以及对在所述数字图像(132)中的像素中表示在所述眼睛(280)的所述被成像区域中的第二眼睛特征的第二像素(136-2)的指定;以及
使用所述转换因子(CFh,CFv)的所确定的值来将在所述数字图像(132)中的所述第一像素(136-1)和所述第二像素(136-2)之间的距离(138)转换成在所述眼睛(280)的所述被成像区域中的所述第一眼睛特征和所述第二眼睛特征之间的距离。
4.根据权利要求1至2中的任一项所述的装置(100),其中,所述眼睛(280)的所述部分包括所述眼睛(280)的视网膜。
5.根据权利要求4所述的装置(100),其中,所述映射(125)将所述扫描参数的值与所述转换因子(CFh,CFv)的相应值相关,所述扫描参数的值跨越由所述眼睛成像装置(200)在所述视网膜的宽视场扫描中覆盖的扫描参数值的范围。
6.根据权利要求1至2和5中的任一项所述的装置(100),其中
所述装置(100)被配置为处理由光学相干断层扫描OCT成像装置的控制器(205)生成的眼科OCT成像数据作为眼睛成像数据,以及
所述扫描参数获得模块(110)被配置为从所述OCT成像数据获得所述OCT成像装置的指示由所述OCT成像装置执行来获取所述眼睛(280)的所述部分的所述被成像区域的数字OCT图像的OCT扫描在所述眼睛(280)中的OCT扫描位置的扫描参数的值。
7.一种处理由眼睛成像装置(200)的控制器(205)生成的眼睛成像数据以评估用于计算在眼睛(280)的部分的被成像区域中的指定眼睛特征之间的距离的转换因子(CFh,CFv)的计算机实现的方法,所述计算机实现的方法包括:
从所述眼睛成像数据获得(S10)所述眼睛成像装置(200)的扫描参数的值,所述扫描参数/>指示由所述眼睛成像装置(200)执行来获取所述眼睛(280)的所述部分的所述被成像区域的数字图像的扫描在所述眼睛(280)中的扫描位置;
使用(S20)所述扫描参数的所获得的值以及在所述扫描参数/>的值和所述转换因子的相应值之间的映射(125)来确定用于将在所述数字图像(132)中的像素之间的距离(138)转换为在所述眼睛(280)的所述被成像区域中的相对应的位置之间的距离的所述转换因子(CFh,CFv)的相应值,其中所述转换因子(CFh,CFv)的每个值指示针对所述扫描参数/>的相应值横越所述眼睛(280)的所述部分的距离关于所述扫描参数/>的函数的模拟变化率,所述像素的相应像素值是在由所述眼睛成像装置(200)使用所述扫描参数/>的不同相应值对所述被成像区域的成像期间获取的;
将所述数字图像与所述转换因子(CFh,CFv)的所确定的值相关联地存储(S30);以及
通过以下操作来生成映射:
运行(S110)穿过所述眼睛成像装置(200)的模型和所述眼睛(280)的模型的光线传播的计算机模拟,以生成使在所述扫描参数的值的序列中的所述扫描参数/>的每个值与在所述眼睛(280)的模型的部分中的相对应的位置相关的模拟结果,当所述眼睛成像装置(200)的模型根据所述扫描参数/>的值操作时,穿过所述眼睛成像装置(200)的模型传播的光线射到所述相对应的位置上;
对于所述扫描参数的值中的每个值,计算(S120)在所述眼睛(280)的模型的所述部分中的所述相对应的位置和在所述眼睛(280)的模型的所述部分中的对应于在值的所述序列中的相邻值的位置之间的距离;以及
将所计算出的距离中的每个距离除以(S130)在所述扫描参数的相对应的值的函数的值和所述扫描参数/>在所述序列中的所述相邻值的函数的值之间的差,以对于所述扫描参数/>的值中的每个值生成指示横越所述眼睛(280)的所述部分的距离关于所述扫描参数/>的函数的模拟变化率的所述转换因子(CFh,CFv)的相应值。
8.根据权利要求7所述的计算机实现的方法,其中,所述计算机模拟使用对所述眼睛(280)的生物测量结果来对穿过所述眼睛(280)的模型和所述眼睛成像装置(200)的模型的光线传播建模。
9.根据权利要求7至8中的任一项所述的计算机实现的方法,还包括通过以下操作来确定在所述被成像区域中的所述指定眼睛特征之间的距离:
接收(S40)对在所述数字图像(132)中的像素中表示在所述眼睛(280)的所述被成像区域中的第一眼睛特征的第一像素(136-1)的指定以及对在所述数字图像(132)中的像素中表示在所述眼睛(280)的所述被成像区域中的第二眼睛特征的第二像素(136-2)的指定;以及
使用所述转换因子(CFh,CFv)的所确定的值来将在所述数字图像(132)中的所述第一像素(136-1)和所述第二像素(136-2)之间的距离(138)转换成在所述眼睛(280)的所述被成像区域中的所述第一眼睛特征和所述第二眼睛特征之间的距离。
10.根据权利要求7至8中的任一项所述的计算机实现的方法,其中,所述眼睛的所述部分包括所述眼睛(280)的视网膜,以及所述映射将所述扫描参数的值与所述转换因子(CFh,CFv)的相应值相关,所述扫描参数/>的值跨越由所述眼睛成像装置(200)在所述视网膜的宽视场扫描中覆盖的扫描参数/>值的范围。
11.根据权利要求7至8中的任一项所述的计算机实现的方法,其中,由光学相干断层扫描OCT成像装置的控制器(205)生成的眼科OCT成像数据被作为所述眼睛成像数据处理,以及
所述获得(S10)包括从所述OCT成像数据获得所述OCT成像装置的指示由所述OCT成像装置执行以获取所述眼睛(280)的所述部分的所述被成像区域的数字OCT图像的OCT扫描在所述眼睛(280)中的OCT扫描位置的扫描参数的值。
12.一种非暂时性计算机可读存储介质,其存储包括计算机可读指令的计算机程序(390),所述计算机可读指令当由计算机处理器(320)执行时使计算机处理器(320)执行根据权利要求7至11中的任一项所述的方法。
13.一种映射生成器模块(140),所述映射生成器模块被配置成生成映射,所述映射将眼睛成像装置(200)的扫描参数的值与转换因子(CFh,CFv)的相应值相关,所述扫描参数/>指示在眼睛(280)的部分中的相应扫描位置,所述眼睛成像装置(200)能够操作来在所述相应扫描位置处获取所述眼睛(280)的所述部分的相应被成像区域的数字图像,所述转换因子(CFh,CFv)用于计算在所述眼睛(280)的所述部分的所述相应被成像区域中的指定眼睛特征之间的距离,所述映射生成器模块(140)被配置成通过执行下列操作来生成所述映射:
运行(S110)穿过所述眼睛成像装置(200)的模型和所述眼睛(280)的模型的光线传播的计算机模拟以生成模拟结果,所述模拟结果将在所述扫描参数的值的序列中的所述扫描参数/>的每个值与在所述眼睛(280)的模型的部分中的相对应的位置相关,当所述眼睛成像装置(200)的模型根据所述扫描参数/>的值操作时,穿过所述眼睛成像装置(200)的模型传播的光线射到所述相对应的位置上;
对于所述扫描参数的值中的每个值,计算(S120)在所述眼睛(280)的模型的所述部分中的所述相对应的位置和在所述眼睛(280)的模型的所述部分中的对应于在值的所述序列中的相邻值的位置之间的距离;以及
将计算出的距离中的每个距离除以(S130)在所述扫描参数的相对应的值的函数的值和所述扫描参数在所述序列中的所述相邻值的函数的值之间的差,以对于所述扫描参数/>的值中的每个值生成所述转换因子(CFh,CFv)的指示横越所述眼睛(280)的所述部分的距离关于所述扫描参数/>的函数的变化率的相应值。
14.一种生成将眼睛成像装置(200)的指示在眼睛(280)的部分中的相应扫描位置的扫描参数的值与用于计算在所述眼睛(280)的所述部分的相应被成像区域中的指定眼睛特征之间的距离的转换因子(CFh,CFv)的相应值相关的映射的计算机实现的方法,所述眼睛成像装置(200)能够操作来在所述相应扫描位置处获取所述眼睛(280)的所述部分的所述相应被成像区域的数字图像,所述方法包括:
运行(S110)穿过所述眼睛成像装置(200)的模型和所述眼睛(280)的模型的光线传播的计算机模拟,以生成使在所述扫描参数的值的序列中的所述扫描参数/>的每个值与在所述眼睛(280)的模型的部分中的相对应的位置相关的模拟结果,当所述眼睛成像装置(200)的模型根据所述扫描参数/>的值操作时,穿过所述眼睛成像装置(200)的模型传播的光线射到所述相对应的位置上;
对于所述扫描参数的值中的每个值,计算(S120)在所述眼睛(280)的模型的所述部分中的所述相对应的位置和在所述眼睛(280)的模型的所述部分中的对应于在值的所述序列中的相邻值的位置之间的距离;以及
将所计算出的距离中的每个距离除以在所述扫描参数的相对应值的函数的值和所述扫描参数在所述序列中的所述相邻值的函数的值之间的差,以对于所述扫描参数的值中的每个值生成指示横越所述眼睛(280)的所述部分的距离关于所述扫描参数/>的函数的变化率的所述转换因子(CFh,CFv)的相应值。
15.一种非暂时性计算机可读存储介质,其存储包括计算机可读指令的计算机程序(390),所述计算机可读指令当由计算机处理器(320)执行时使所述计算机处理器(320)执行根据权利要求14所述的方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP20166753.2A EP3889889A1 (en) | 2020-03-30 | 2020-03-30 | Ocular image data processing |
EP20166753.2 | 2020-03-30 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113456018A CN113456018A (zh) | 2021-10-01 |
CN113456018B true CN113456018B (zh) | 2024-04-05 |
Family
ID=70058245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110339136.5A Active CN113456018B (zh) | 2020-03-30 | 2021-03-30 | 眼睛图像数据处理 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20210298592A1 (zh) |
EP (1) | EP3889889A1 (zh) |
JP (1) | JP7267337B2 (zh) |
CN (1) | CN113456018B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113436211B (zh) * | 2021-08-03 | 2022-07-15 | 天津大学 | 一种基于深度学习的医学图像活动轮廓分割方法 |
TWI800991B (zh) * | 2021-11-18 | 2023-05-01 | 宏碁股份有限公司 | 視野檢測方法及視野檢測系統 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7278989B2 (en) * | 2002-04-30 | 2007-10-09 | Nidek Co. Ltd. | Method for analysis of an ablated surface, analysis on inconsistency in a laser beam, and calibration of laser beam irradiation data |
EP1858402B1 (en) * | 2005-01-21 | 2017-11-29 | Massachusetts Institute Of Technology | Methods and apparatus for optical coherence tomography scanning |
PL1753373T3 (pl) * | 2005-04-05 | 2008-11-28 | Alcon Inc | Soczewka wewnątrzgałkowa |
US10362959B2 (en) * | 2005-12-06 | 2019-07-30 | St. Jude Medical, Atrial Fibrillation Division, Inc. | System and method for assessing the proximity of an electrode to tissue in a body |
DE102006002001B4 (de) * | 2006-01-16 | 2009-07-23 | Sensomotoric Instruments Gmbh | Verfahren zur Bestimmung der räumlichen Relation eines Auges einer Person bezüglich einer Kameravorrichtung |
US7758189B2 (en) * | 2006-04-24 | 2010-07-20 | Physical Sciences, Inc. | Stabilized retinal imaging with adaptive optics |
GB0622325D0 (en) | 2006-11-09 | 2006-12-20 | Optos Plc | Improvements in or relating to retinal scanning |
US8550625B2 (en) * | 2007-08-01 | 2013-10-08 | Amo Development, Llc | Systems and methods for fine-tuning refractive surgery |
US8348429B2 (en) * | 2008-03-27 | 2013-01-08 | Doheny Eye Institute | Optical coherence tomography device, method, and system |
US8820931B2 (en) * | 2008-07-18 | 2014-09-02 | Doheny Eye Institute | Optical coherence tomography-based ophthalmic testing methods, devices and systems |
ES2346175B1 (es) * | 2009-04-08 | 2011-09-30 | Consejo Superior De Investigaciones Científicas (Csic) | Instrumento para la simulacion de correcciones oftalmicas multifocales. |
JP5297415B2 (ja) * | 2010-04-30 | 2013-09-25 | キヤノン株式会社 | 眼科装置及び眼科方法 |
DE102010046500A1 (de) * | 2010-09-24 | 2012-03-29 | Carl Zeiss Meditec Ag | Verfahren und Vorrichtung zur Aufnahme und Darstellung eines OCT-Ganzaugenscans |
WO2012178054A1 (en) * | 2011-06-23 | 2012-12-27 | Amo Development, Llc | Ophthalmic range finding |
GB201307990D0 (en) * | 2013-05-02 | 2013-06-12 | Optos Plc | Improvements in and relating to imaging of the eye |
EP3054907B1 (en) * | 2013-10-08 | 2019-04-03 | Optimedica Corporation | Laser eye surgery system calibration |
US10660519B2 (en) * | 2014-01-30 | 2020-05-26 | Duke University | Systems and methods for eye tracking for motion corrected ophthalmic optical coherence tomography |
JP6444666B2 (ja) * | 2014-09-08 | 2018-12-26 | 株式会社トプコン | 眼科撮影装置および眼科情報処理装置 |
JP6624641B2 (ja) * | 2016-03-18 | 2019-12-25 | 株式会社トプコン | 眼科装置 |
JP6740177B2 (ja) * | 2017-06-14 | 2020-08-12 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
EP3566637A1 (en) * | 2018-05-11 | 2019-11-13 | Optos PLC | Oct image processing |
JP7213048B2 (ja) * | 2018-09-25 | 2023-01-26 | 株式会社トプコン | 眼科情報処理装置、眼科装置、及び眼科情報処理方法 |
-
2020
- 2020-03-30 EP EP20166753.2A patent/EP3889889A1/en active Pending
-
2021
- 2021-03-19 US US17/207,077 patent/US20210298592A1/en active Pending
- 2021-03-30 CN CN202110339136.5A patent/CN113456018B/zh active Active
- 2021-03-30 JP JP2021058560A patent/JP7267337B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
EP3889889A1 (en) | 2021-10-06 |
US20210298592A1 (en) | 2021-09-30 |
CN113456018A (zh) | 2021-10-01 |
JP2021168909A (ja) | 2021-10-28 |
JP7267337B2 (ja) | 2023-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8442356B2 (en) | Methods, systems and computer program products for analyzing three dimensional data sets obtained from a sample | |
US20050267366A1 (en) | Ultrasonic diagnostic apparatus and image processing method | |
CN110473265B (zh) | Oct图像处理 | |
CN113456018B (zh) | 眼睛图像数据处理 | |
US11259696B2 (en) | Imaging of the eye | |
JP6436442B2 (ja) | 光音響装置および画像処理方法 | |
CN111629670A (zh) | 用于超声系统的回波窗口伪影分类和视觉指示符 | |
JP2016154858A (ja) | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム | |
US20170251999A1 (en) | Ultrasound diagnostic image generating device and method | |
US11278259B2 (en) | Thrombus detection during scanning | |
JP2018020192A (ja) | 画像処理システム、処理方法及びプログラム | |
JPH11151206A (ja) | 眼底画像の解析方法及びその装置 | |
JP6892234B2 (ja) | 眼科撮影装置 | |
JP2018046959A (ja) | 眼科撮影装置 | |
JP6336271B2 (ja) | 検眼鏡の改良 | |
JP2021166906A (ja) | 眼科撮影装置 | |
EP3781014B1 (en) | Post-processing method to improve lso-based tracking in oct | |
JP2018068778A (ja) | 眼科oct解析装置及び眼科解析プログラム | |
JP6598963B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
US20230320583A1 (en) | Optical coherence tomography angiography data processing for reducing projection artefacts | |
JP6720267B2 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
JP2021049205A (ja) | 画像処理装置、画像処理方法 |
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 |