CN115886721B - 一种眼球活动度的评估方法、系统和存储介质 - Google Patents
一种眼球活动度的评估方法、系统和存储介质 Download PDFInfo
- Publication number
- CN115886721B CN115886721B CN202210989508.3A CN202210989508A CN115886721B CN 115886721 B CN115886721 B CN 115886721B CN 202210989508 A CN202210989508 A CN 202210989508A CN 115886721 B CN115886721 B CN 115886721B
- Authority
- CN
- China
- Prior art keywords
- eye
- eyeball
- bitmap
- straight line
- angle
- 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
- 210000005252 bulbus oculi Anatomy 0.000 title claims abstract description 117
- 230000000694 effects Effects 0.000 title claims abstract description 33
- 238000011156 evaluation Methods 0.000 title claims abstract description 10
- 210000001508 eye Anatomy 0.000 claims abstract description 178
- 210000001747 pupil Anatomy 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000006073 displacement reaction Methods 0.000 claims description 28
- 230000008859 change Effects 0.000 claims description 15
- 238000003384 imaging method Methods 0.000 claims description 13
- 230000004424 eye movement Effects 0.000 claims description 10
- 208000014733 refractive error Diseases 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 5
- 238000013527 convolutional neural network Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 13
- 210000004087 cornea Anatomy 0.000 description 12
- 210000000554 iris Anatomy 0.000 description 11
- 210000003128 head Anatomy 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 210000005036 nerve Anatomy 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 210000003205 muscle Anatomy 0.000 description 4
- 210000003786 sclera Anatomy 0.000 description 4
- 230000008602 contraction Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- XUMBMVFBXHLACL-UHFFFAOYSA-N Melanin Chemical compound O=C1C(=O)C(C2=CNC3=C(C(C(=O)C4=C32)=O)C)=C2C4=CNC2=C1C XUMBMVFBXHLACL-UHFFFAOYSA-N 0.000 description 2
- 206010033799 Paralysis Diseases 0.000 description 2
- 208000004350 Strabismus Diseases 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000004418 eye rotation Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 208000003164 Diplopia Diseases 0.000 description 1
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 208000019749 Eye movement disease Diseases 0.000 description 1
- 208000029549 Muscle injury Diseases 0.000 description 1
- 208000021642 Muscular disease Diseases 0.000 description 1
- 201000009623 Myopathy Diseases 0.000 description 1
- 208000036110 Neuroinflammatory disease Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 208000029444 double vision Diseases 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003959 neuroinflammation Effects 0.000 description 1
- 208000028780 ocular motility disease Diseases 0.000 description 1
- 239000000049 pigment Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000011895 specific detection Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition 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/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/113—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining or recording eye movement
-
- 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/14—Arrangements specially adapted for eye photography
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Human Computer Interaction (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明涉及一种眼球活动度的评估方法、系统和存储介质。所述评估方法包括:在700‑1200nm的近红外光光场中,在眼睛前方位置拍摄获取用户平视前方的第一眼位图;在眼睛前方位置拍摄获取用户眼球沿待测方向运动到极限位置的第二眼位图和第三眼位图;将第一眼位图分别与第二眼位图和第三眼位图进行对比,计算眼球活动度的活动角。本发明提供的技术方案在近红外光光场中拍摄眼部图像,有效地区分虹膜、瞳孔,从而准确地识别瞳孔边缘,再采用不同眼位图的对比计算眼球活动角,并通过补偿进一步提高活动角的计算精度,从而准确地计算出眼球在各个待测方向上的活动角。
Description
技术领域
本申请属于眼部检测技术领域,具体涉及一种眼球活动度的评估方法、装置和存储介质。
背景技术
动眼、滑车、外展神经具有支配眼球眼外肌运动的功能,称为眼球运动神经。当上述神经或神经核单独或合并受损时,可出现眼球运动不能或复视,完全损害时出现眼外肌全部瘫痪,眼球固定不动。眼外肌损伤、感染或肌病导致的眼外肌麻痹,也可出现眼球运动不能,临床统称为眼球运动障碍。眼球运动障碍还可能与内科疾病相关,包括眼眶疾病、糖尿病、神经炎症等,因此如何评估单眼及双眼的眼球转动能力,对于眼部检测具有重要的意义。
现有技术中,眼球活动度的检测通常采用医生手动测量的方式,具体检测方法如下:1)单眼的眼球运动检查:遮盖对侧眼,自第一眼位开始,嘱患者注视手电筒,沿注视的诊断方向运动;2)双眼的眼球运动检查:自第一眼位开始,嘱患者注视手电筒,沿注视的诊断方向运动。双眼运动检查,可以评估眼球运动的双眼相对位置,获取单眼检查不同的信息。但是,上述检测方法,没有客观的数据,且无法记录患者每次的变化,无法形成连续数据。
此外,现有技术中还公开了采用自动检测方式来检测眼球活动度。如中国发明专利“一种眼球活动度检测仪”(申请号:202011260253.4,公开日:2021年2月19日)所述,通过头部固定架固定患者的头部,指示灯设置在头部固定架前方,所述指示灯朝向远离受试者眼部的方向上发出可见光束,可见光束延伸到受试者静视野之外,成像镜头设置在头部固定架前方,成像镜头用于拍摄受试者眼球沿着可见光束的方向转动时的图像,通过成像镜头拍摄的眼部图像获得受试者的眼球活动情况,但是,上述方案中,未对眼部图像做具体的区分和处理,导致眼球活动情况的检测精度不够;如中国专利“一种基于计算机的眼球运动距离及双眼运动一致性偏差检测装置及方法”(申请号:201710054692.1,公开日:2017年6月13日)所述,通过角膜缘提取单元提取照相机发来的基准照片和各个测试照片中的角膜缘,通过角膜缘上提取的点来计算眼球在各个方向上的运动距离,从而进行双眼运动一致性偏差计算,但是,由于角膜缘为角膜和巩膜之间的过渡带,提取到的角膜缘的边界不够清晰,也会影响到最终结果的精度。
发明内容
本申请的目的是提供一种眼球活动度的评估方法和系统,旨在提高眼球活动度评估的精度,消除测量过程中的误差。通过采集用户在近红外光光场的多个不同眼位图,再对采集的多个不同眼位图像进行分析,计算眼球活动度的活动角,进而实现对眼球活动度的准确评估。
为了实现上述目的,一方面,本申请提供一种眼球活动度的评估方法,所述评估方法包括以下步骤:
步骤S1:在700-1200nm的近红外光光场中,在眼球前方位置拍摄获取用户平视前方的第一眼位图;
步骤S2:在眼球前方位置拍摄获取用户眼球沿待测方向运动到极限位置的第二眼位图和第三眼位图;
步骤S3:将所述第一眼位图分别与所述第二眼位图和第三眼位图进行对比,计算眼球活动度的活动角。
进一步地,所述待测方向包括眼球向上、向下、向内、向外、向内上、向外上、向内下、向外下。
进一步地,所述步骤S3包括:采用卷积神经网络,从所述第一眼位图、第二眼位图、第三眼位图的图像中分割瞳孔;
取出第一眼位图的圆心C1、第二眼位图或第三眼位图的圆心C2,并用直线连接圆心C1和圆心C2;
沿所述直线找到所述第一眼位图与所述第二眼位图或第三眼位图的瞳孔的同一点运动前后的边缘点P1、P2;
将所述边缘点P1和P2连接成直线P1P2;
通过垂直于所述直线P1P2并经过所述圆心C1的直线将所述直线P1P2分割成两段,即直线A和直线B;
通过眼球半径r,获取直线A对应的圆心角α以及直线B对应的圆心角β,将所述圆心角α和所述圆心角β相加得到所述活动角。
进一步地,通过第一眼位图、第二眼位图或第三眼位图中瞳孔变化的幅度、眼球中心位移的幅度、角膜折射误差补偿对所述活动角进行补偿。
进一步地,通过第一眼位图、第二眼位图或第三眼位图中瞳孔变化的幅度对所述活动角进行补偿具体为:通过获取第一眼位图、第二眼位图或第三眼位图中瞳孔变化的幅度补偿所述边缘点P2的位置。
具体地,计算第一眼位图、第二眼位图或第三眼位图中瞳孔直径沿着与待测方向垂直的方向变化的幅度,利用所述幅度等比例地补偿所述边缘点P2的位置。
进一步地,通过第一眼位图、第二眼位图或第三眼位图中眼球中心位移的幅度对所述活动角进行补偿具体为:获取第一眼位图、第二眼位图或第三眼位图中眼球中心不同程度位移形成的重叠成像,将所述重叠成像进行逆向平移以抵消所述位移对活动角的影响。
进一步地,通过第一眼位图、第二眼位图或第三眼位图中角膜折射误差补偿对所述活动角进行补偿具体为:通过角膜折射误差补偿所补偿的角度与所述活动角成线性关系。
一方面,本申请提供了一种眼球活动度的评估装置,包括:
第一获取模块,用于在700-1200nm的近红外光光场中,在眼球前方位置拍摄获取用户平视前方的第一眼位图;
第二获取模块,用于在眼球前方位置拍摄获取用户眼球沿待测方向运动到极限位置的第二眼位图和第三眼位图;
图像处理模块,用于将该第一眼位图分别与第二眼位图和第三眼位图进行对比,计算眼球活动度的活动角。
一方面,本申请提供了一种计算机可读存储介质,该计算机可读存储介质中存储有至少一条程序代码,所述至少一条程序代码由处理器加载并执行以实现该眼球活动度的评估方法所执行的操作。
本申请的技术方案相较于现有技术,至少具有如下有益效果:
在近红外光光场中拍摄眼部图像,可以有效地区分虹膜、瞳孔,从而准确地识别瞳孔边缘,通过采用瞳孔边缘点的位移计算眼球活动角,避免了现有技术中采用虹膜边缘计算眼球活动角带来的误差大的缺点。
进一步地,本申请中采用眼球运动前后瞳孔变化的幅度、眼球中心位移的幅度、角膜折射误差补偿对所述活动角进行补偿,可以近一步提高对于活动角的计算精度,从而准确地计算出眼球在各个待测方向上的活动角,从而为眼部检测提供更为准确的判定依据。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例或现有技术描述中所需要使用的附图作简单的介绍。
图1为本申请中的眼球活动度的评估方法的流程图;
图2为可见光与近红外光下拍摄到的眼部图像;
图3为在各个待测方向上拍摄的第二眼位图或第三眼位图与第一眼位图的重叠图像;
图4为眼球活动角的计算图;
图5为眼球运动前后的瞳孔变化图;
图6为眼球运动前后眼球成像的重叠图像;
图7为眼球运动前后眼球中心的位移图;
图8是角膜晶体对瞳孔成像的折射图;
图9为采用瞳孔边缘和虹膜边缘计算的活动角;
图10为眼球活动度的评估装置的结构示意图。
具体实施方式
下面结合附图和实施例,对本申请的具体实施方式作进一步详细描述。以下实施例用于说明本申请,但不用来限制本申请的范围。
一方面,如图1所示,本实施例1提供了一种眼球活动度的评估方法,该评估方法包括以下步骤:
步骤S1:在700-1200nm的近红外光光场中,用户平视前方,获取在用户的眼睛前方拍摄的第一眼位图;
步骤S2:获取在所述用户的眼睛前方拍摄的眼球沿待测方向运动到极限位置的第二眼位图、第三眼位图;
步骤S3:将所述第一眼位图分别与所述第二眼位图、第三眼位图进行对比,计算眼球活动度的活动角。
在具体步骤中,采用获取装置获取眼部图像时,先固定用户的头部,如通过额托、颌托、眼角固定点实现用户的头部固定。固定好之后,在700-1200的红外红光光场中拍摄用户眼睛正视前方的第一眼位图,以及用户随着各个待测方向的指示灯指示方向将眼球转动到极限位置的第二眼位图、第三眼位图。
需要说明的是,眼位是指在眼科检查中眼球的位置,分为第一眼位、第二眼位、第三眼位。第一眼位是指两眼在水平面上平视前方无限远处的眼位,第二眼位是指眼球向上、下、内、外转动时的眼位,第三眼位是指眼球向内上、内下、外上、外下作斜向转动时的眼位,即向鼻上、鼻下、颞上、颞下转动时的眼位。相应地,第一、第二、第三眼位图分别指眼球在各个眼位时拍摄的图像。
如图2所示,在正常可见光的波段400-700nm中,眼睛不同部分,即瞳孔、虹膜、巩膜的颜色对成像几乎不起作用,由于虹膜与巩膜接触部分的角膜缘的渐变式结构,导致无法准确地识别眼球的球心。人类黑色素色素的吸收峰值发生在335nm左右,而对于波长超过700nm的波段几乎完全不吸收,而虹膜的反射率在波长超过700nm的近红外线波段内是相当稳定的,因此,采用近红外光光场,可以很好地区分巩膜、虹膜、瞳孔边界,从而可以更好地提升算法的准确性与稳定性。
如图3所示,待测方向包括眼球向上、向下、向内、向外、向内上、向外上、向内下、向外下,每只眼睛分别在八个待测方向上测量活动角,每次测量一只眼睛时遮挡另一只眼睛,采用各个待测方向上设置的指示灯分别指示左眼、右眼重复上述各个方向的转动,并获取多张第二眼位图、第三眼位图,以分别与作为基准的第一眼位图进行对比,从而计算出眼球活动度的活动角。
需要说明的是,上述活动角是指眼球活动前后眼球体绕眼球球心转动的最大角度。
实施例2
如图4所示,基于上述方法的实施例1,本实施例2进一步限定了活动角的计算方式。具体为:
采用卷积神经网络,从所述第一眼位图、眼球向下运动到极限位置的第二眼位图的图像中分割瞳孔;
取出第一眼位图的圆心C1和第二眼位图的圆心C2;用直线连接圆心C1和圆心C2;
沿所述直线找到所述第一眼位图与所述第二眼位图的瞳孔的同一点运动前后的边缘点P1、P2;
将所述边缘点P1和P2连接成直线P1P2;
通过垂直于所述直线P1P2并经过所述圆心C1的直线将所述直线P1P2分割成两段,即直线A和直线B;
通过眼球半径r,获取直线A对应的圆心角α以及直线B对应的圆心角β,将所述圆心角α和所述圆心角β相加得到所述活动角θ。
上述计算方式中,通过简单的解析几何关系,通过计算球面上弧长所对应的圆心角,即可求出眼球运动前后的活动角。需要说明的是,由于眼球为规则的球面结构,各个待测方向,例如眼球斜向运动的内上、内下、外上、外下方向,均可在获取第一、第二、第三眼位图后,通过图像和数据处理步骤实现上述活动角的计算。
实施例3
如图5-7所示,基于上述方法的实施例2,本实施例3进一步限定了活动角的补偿方式,从而能够提高活动角的计算精度。在实际测量中发现,眼球运动前后瞳孔的变化、眼球中心的位移、角膜折射的影响,均会对活动角的计算带来误差,下面将具体解释对上述几种误差的补偿方式。
首先,如图5所示,当眼球转动偏离了正视时,瞳孔成像会变成一个椭圆,变扁方向在转动的方向上,而垂直于转动的方向成像不受影响。所以如果垂直于转动的方向上的瞳孔直径有变化,则来自于瞳孔的变化。如果正视比斜视的直径小,则瞳孔转动时有放大,如果正视比斜视的直径大,则瞳孔转动时有缩小。例如,如果垂直于转动方向上的椭圆有缩小,就来自于瞳孔的收缩,将分割出来的两个圆心连出直线,取出两个椭圆上垂直于直线的直径,拿上述两个直径计算瞳孔收缩的比例,等比例地矫正直线上用来计算的椭圆边缘点,即可通过瞳孔变化的幅度补偿所述边缘点P2的位置,借此消除瞳孔收缩对活动角计算精度的影响,从而通过第一眼位图、第二眼位图或第三眼位图中瞳孔变化的幅度对所述活动角进行补偿。
如图6所示,当眼睛在转动时,眼睛并非在空间中悬浮不动,由于周遭都是软组织还有眼肌的发力不均,不同方向的眼球转动会引起不同程度的眼球中心的位移。眼球转动角度越大,眼球中心的位移就越多。眼球中心的位移对算法的影响相当于在重叠对齐正视和侧视的成像时发生了错位,导致计算出来的活动角过大或过小。补偿的方法是知道眼球中心位移多少后,重叠成像时有意反方向错位,抵消掉位移的影响。
如图7所示,图7的左图表示从人体颅顶方向所拍摄的眼球向内、外方向转动前后的变化图,图7的右图表示从人体颅顶方向所拍摄的眼球向上、下方向转动前后的变化图。由于眼球转动前后,前后方向的眼球中心的位移对正前方相机拍出的成像没有影响,因此可以忽略不计,只取眼球投影到横截面上的图像,统计显示眼球运动的方向和眼球在横截面上的中心位移的关系如下:
眼球向内旋转30°的位移:眼球中心向内0.69mm;
眼球向外旋转30°的位移:眼球中心向外0.45mm;
眼球向上旋转20°的位移:眼球中心向下0.43mm;
眼球向下旋转20°的位移:眼球中心向上0.43mm;
并且,眼球中心位移和眼球转动角度为线性规律,如果眼球转动角度为0°时(正视)则眼球中心位移为0mm,可以通过预先统计的上述位移和角度的关系来修正活动角计算中眼球球心C2的位置,从而通过第一眼位图、第二眼位图或第三眼位图中眼球中心位移的幅度对所述活动角进行补偿。
如图8所示,在计算活动角时,瞳孔的边缘点受到角膜晶体的折射影响,画面中瞳孔边缘点并不是现实真实的位置,而虹膜边缘点则因不穿透角膜而不会受其影响。同照片采用瞳孔边缘点与虹膜边缘点两种算法,可得出不同的活动角,其中采用瞳孔边缘点进行计算是带有角膜折射误差的。通过实验测量可得,通过瞳孔和虹膜计算出的活动角的差别成高度线性,遂得对角膜折射误差进行补偿的活动角计算公式如下:
活动角θ的补偿=0.13134×活动角θ+0.52704
修正后的活动角θ=活动角θ+活动角θ的补偿
上述公式中,活动角θ为未经角膜折射误差补偿前的计算活动角。
将上述眼球运动前后瞳孔的变化、眼球中心的位移、角膜折射的影响对活动角的影响均考虑在内并进行相应的补偿后,即可消除活动角计算中的误差,最终得到精确的活动角。
一方面,本申请提供了一种眼球活动度的评估装置,如图10所示,包括:第一获取模块1001、第二获取模块1002和图像处理模块(1003),其中
第一获取模块1001,用于在700-1200nm的近红外光光场中,在眼球前方位置拍摄获取用户平视前方的第一眼位图;
第二获取模块1002,用于在眼球前方位置拍摄获取用户眼球沿待测方向运动到极限位置的第二眼位图和第三眼位图;
图像处理模块1003,用于将该第一眼位图分别与第二眼位图和第三眼位图进行对比,计算眼球活动度的活动角。
一方面,本申请提供了一种计算机可读存储介质,该计算机可读存储介质中存储有至少一条程序代码,所述至少一条程序代码由处理器加载并执行以实现该眼球活动度的评估方法所执行的操作。
关于上述实施例中的眼球活动度的评估装置,其中各个模块执行操作的具体方式已经在方法实施例中进行了详细描述,相关之处参见方法实施例的部分说明。
在示例性实施例中,还提供了一种计算机可读存储介质,包括存储有至少一条程序代码的存储器,所述至少一条程序代码由处理器加载并执行以完成上述实施例中的眼球活动度的评估方法。例如,该计算机可读存储介质可以是只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、只读光盘(Compact Disc Read-OnlyMemory,CDROM)、磁带、软盘和光数据存储设备等。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来至少一条程序代码相关的硬件完成,该程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本申请的较佳实施例而已,并不用以限制本申请,凡在本申请的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本申请的保护范围之内。
Claims (9)
1.一种眼球活动度的评估方法,其特征在于,所述评估方法包括:
步骤S1:在700-1200nm的近红外光光场中,在眼球前方位置拍摄获取用户平视前方的第一眼位图;
步骤S2:在眼球前方位置拍摄获取用户眼球沿待测方向运动到极限位置的第二眼位图和第三眼位图;
步骤S3:将所述第一眼位图分别与所述第二眼位图和第三眼位图进行对比,计算眼球活动度的活动角;
所述步骤S3包括:
采用卷积神经网络,从所述第一眼位图、第二眼位图、第三眼位图的图像中分割瞳孔;
取出第一眼位图的圆心C1 、第二眼位图或第三眼位图的圆心C2,并用直线连接圆心C1和圆心C2;
沿所述直线找到所述第一眼位图与所述第二眼位图或第三眼位图的瞳孔的同一点运动前后的边缘点P1、P2;
将所述边缘点P1和P2连接成直线P1P2;
通过垂直于所述直线P1P2并经过所述圆心C1的直线将所述直线P1P2分割成两段,即直线A和直线B;
通过眼球半径r,获取直线A对应的圆心角α以及直线B对应的圆心角β,将所述圆心角α和所述圆心角β相加得到所述活动角。
2.根据权利要求1所述的眼球活动度的评估方法,其特征在于,
所述待测方向包括眼球向上、向下、向内、向外、向内上、向外上、向内下、向外下。
3.根据权利要求1或2所述的眼球活动度的评估方法,其特征在于,通过第一眼位图、第二眼位图或第三眼位图中瞳孔变化的幅度、眼球中心位移的幅度、角膜折射误差补偿对所述活动角进行补偿。
4.根据权利要求3所述的眼球活动度的评估方法,其特征在于,通过第一眼位图、第二眼位图或第三眼位图中瞳孔变化的幅度对所述活动角进行补偿具体为:通过获取所述第一眼位图、所述第二眼位图或第三眼位图中瞳孔变化的幅度补偿所述边缘点P2的位置。
5.根据权利要求4所述的眼球活动度的评估方法,其特征在于,计算第一眼位图、第二眼位图或第三眼位图中瞳孔直径沿着与待测方向垂直的方向变化的幅度,利用所述幅度等比例地补偿所述边缘点P2的位置。
6.根据权利要求3所述的眼球活动度的评估方法,其特征在于,通过第一眼位图、第二眼位图或第三眼位图中眼球中心位移的幅度对所述活动角进行补偿具体为:获取所述第一眼位图、所述第二眼位图或第三眼位图中眼球中心不同程度位移形成的重叠成像,将所述重叠成像进行逆向平移以抵消所述位移对活动角的影响。
7.根据权利要求3所述的眼球活动度的评估方法,其特征在于,通过第一眼位图、第二眼位图或第三眼位图中角膜折射误差补偿对所述活动角进行补偿具体为:通过角膜折射误差补偿所补偿的角度与所述活动角成线性关系。
8.一种眼球活动度的评估系统,其特征在于,包括:
第一获取模块,用于在700-1200nm的近红外光光场中,在眼球前方位置拍摄获取用户平视前方的第一眼位图;
第二获取模块,用于在眼球前方位置拍摄获取用户眼球沿待测方向运动到极限位置的第二眼位图和第三眼位图;
图像处理模块,用于将所述第一眼位图分别与所述第二眼位图和第三眼位图进行对比,计算眼球活动度的活动角;
所述计算包括:
采用卷积神经网络,从所述第一眼位图、第二眼位图、第三眼位图的图像中分割瞳孔;
取出第一眼位图的圆心C1 、第二眼位图或第三眼位图的圆心C2,并用直线连接圆心C1和圆心C2;
沿所述直线找到所述第一眼位图与所述第二眼位图或第三眼位图的瞳孔的同一点运动前后的边缘点P1、P2;
将所述边缘点P1和P2连接成直线P1P2;
通过垂直于所述直线P1P2并经过所述圆心C1的直线将所述直线P1P2分割成两段,即直线A和直线B;通过眼球半径r,获取直线A对应的圆心角α以及直线B对应的圆心角β,将所述圆心角α和所述圆心角β相加得到所述活动角。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有至少一条程序代码,所述至少一条程序代码由处理器加载并执行以实现如权利要求1至7任一项所述的眼球活动的评估方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210989508.3A CN115886721B (zh) | 2022-08-18 | 2022-08-18 | 一种眼球活动度的评估方法、系统和存储介质 |
PCT/CN2023/107803 WO2024037267A1 (zh) | 2022-08-18 | 2023-07-18 | 一种眼球活动度的评估方法、系统和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210989508.3A CN115886721B (zh) | 2022-08-18 | 2022-08-18 | 一种眼球活动度的评估方法、系统和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115886721A CN115886721A (zh) | 2023-04-04 |
CN115886721B true CN115886721B (zh) | 2023-08-01 |
Family
ID=86486135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210989508.3A Active CN115886721B (zh) | 2022-08-18 | 2022-08-18 | 一种眼球活动度的评估方法、系统和存储介质 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN115886721B (zh) |
WO (1) | WO2024037267A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115886721B (zh) * | 2022-08-18 | 2023-08-01 | 上海佰翊医疗科技有限公司 | 一种眼球活动度的评估方法、系统和存储介质 |
CN117373103B (zh) * | 2023-10-18 | 2024-05-07 | 北京极溯光学科技有限公司 | 一种图像特征提取方法、装置、设备及存储介质 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5158842B2 (ja) * | 2007-04-24 | 2013-03-06 | 学校法人中部大学 | 眼球運動計測方法および眼球運動計測装置 |
CN102125422A (zh) * | 2010-01-12 | 2011-07-20 | 北京科技大学 | 视线追踪系统中基于瞳孔-角膜反射的视线估计方法 |
CN104814717B (zh) * | 2015-04-14 | 2016-09-07 | 赵桂萍 | 一种补偿式消除变体位误差的眼震全图的检测方法和装置 |
ES2887976T3 (es) * | 2017-09-08 | 2021-12-29 | Tobii Ab | Compensación del radio de la pupila |
CN108354584B (zh) * | 2018-03-06 | 2020-12-29 | 京东方科技集团股份有限公司 | 眼球追踪模组及其追踪方法、虚拟现实设备 |
JP7198023B2 (ja) * | 2018-09-27 | 2022-12-28 | 株式会社アイシン | 眼球情報推定装置、眼球情報推定方法および眼球情報推定プログラム |
CN114529976A (zh) * | 2020-10-30 | 2022-05-24 | 奥佳华瑞(厦门)医疗科技有限公司 | 一种基于图像处理的眼球活动度检测方法 |
CN112426123A (zh) * | 2020-11-23 | 2021-03-02 | 上海交通大学医学院附属第九人民医院 | 基于眼球运动度的眼外肌检查方法、装置、设备、及仪器 |
CN114391805A (zh) * | 2022-02-09 | 2022-04-26 | 复旦大学附属中山医院 | 一种实时眼球生物数据测量装置及测量方法 |
CN115886721B (zh) * | 2022-08-18 | 2023-08-01 | 上海佰翊医疗科技有限公司 | 一种眼球活动度的评估方法、系统和存储介质 |
-
2022
- 2022-08-18 CN CN202210989508.3A patent/CN115886721B/zh active Active
-
2023
- 2023-07-18 WO PCT/CN2023/107803 patent/WO2024037267A1/zh unknown
Also Published As
Publication number | Publication date |
---|---|
CN115886721A (zh) | 2023-04-04 |
WO2024037267A1 (zh) | 2024-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dubbelman et al. | The shape of the aging human lens: curvature, equivalent refractive index and the lens paradox | |
CN115886721B (zh) | 一种眼球活动度的评估方法、系统和存储介质 | |
Eggert | Eye movement recordings: methods | |
ES2861259T3 (es) | Método y sistema de prescripción de gafas | |
Otero-Millan et al. | Knowing what the brain is seeing in three dimensions: A novel, noninvasive, sensitive, accurate, and low-noise technique for measuring ocular torsion | |
US20130235346A1 (en) | System and methods for documenting and recording of the pupillary red reflex test and corneal light reflex screening of the eye in infants and young children | |
Ramasubramanian et al. | Objective measurement of accommodative biometric changes using ultrasound biomicroscopy | |
JP7166473B2 (ja) | 眼検査 | |
Sousa de Almeida et al. | Computer-aided methodology for syndromic strabismus diagnosis | |
Demir et al. | Comparison of normal and keratoconic corneas by Galilei Dual-Scheimpflug Analyzer | |
US7110582B1 (en) | Method for determining binocular balance and disorders of binocularity of an individual or clinical groups of individuals | |
CN110367925A (zh) | 主客观一体式诊断性验光装置及验光方法 | |
CN115590462A (zh) | 一种基于摄像头的视力检测方法和装置 | |
Reddy et al. | Corneal topography and corneal thickness in children | |
Hirota et al. | Automatic recording of the target location during smooth pursuit eye movement testing using video-oculography and deep learning-based object detection | |
KR102209573B1 (ko) | 사시 중등도 평가장치 및 방법 | |
Yazıcı et al. | Measurements of anterior segment parameters using three different non-contact optical devices in keratoconus patients | |
Upadhyaya et al. | Pentacam–A corneal tomography system | |
CN115670370A (zh) | 一种去除眼底图像玻璃体混浊斑的视网膜成像方法和装置 | |
CN115969308A (zh) | 基于超广角眼底成像的远程眼检系统 | |
CN113384236A (zh) | 一种眼部光学相干断层扫描的方法及其装置 | |
Sahu | Corneal topography: sirius | |
Takegami et al. | An algorithm for an eye tracking system with self‐calibration | |
CN112493983B (zh) | 一种间接实现分析人眼内外及全眼波前像差方法 | |
CN112043235B (zh) | 一种便携式眼球静态旋转测量仪及利用其测量眼球旋转角度的方法 |
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 |