[go: up one dir, main page]

CN101513342A - Full-view pupil analysis measurement method - Google Patents

Full-view pupil analysis measurement method Download PDF

Info

Publication number
CN101513342A
CN101513342A CNA2009100217046A CN200910021704A CN101513342A CN 101513342 A CN101513342 A CN 101513342A CN A2009100217046 A CNA2009100217046 A CN A2009100217046A CN 200910021704 A CN200910021704 A CN 200910021704A CN 101513342 A CN101513342 A CN 101513342A
Authority
CN
China
Prior art keywords
pupil
image
full
center
circle
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.)
Pending
Application number
CNA2009100217046A
Other languages
Chinese (zh)
Inventor
张作明
耿佳
郭群
李莉
徐世伟
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.)
Air Force Medical University
Original Assignee
Fourth Military Medical University FMMU
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 Fourth Military Medical University FMMU filed Critical Fourth Military Medical University FMMU
Priority to CNA2009100217046A priority Critical patent/CN101513342A/en
Publication of CN101513342A publication Critical patent/CN101513342A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Eye Examination Apparatus (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明是一种全视野瞳孔分析测量方法,包括如下步骤:步骤300,首先通过视频采集卡将眼部区域的图像数据实时采集,得到眼部区域图像;步骤301,通过均值法或最值法或公式法对300步采集到的图像进行灰度化;步骤302,通过直方图法对步骤301灰度化的图像进行二值化处理,当背景与目标差异大,直方图出现双峰,这时的最低谷点作为阈值点;步骤303,对进行二值化处理的图像瞳孔区域反光进行消除;步骤304确定瞳孔中心,计算其直径,提取瞳孔的圆心,以及确定瞳孔的直径,实现瞳孔中心的跟踪测量。该方法为一种实时性好、使用方便的一种全视野瞳孔测量分析方法。

The present invention is a full-view pupil analysis and measurement method, comprising the following steps: step 300, first collect the image data of the eye region in real time through a video capture card, and obtain an image of the eye region; step 301, through the mean value method or the maximum value method Or the formula method carries out gray scale to the image that 300 steps collect; Step 302, carries out binarization processing to the image of step 301 gray scale by the histogram method, when background and target difference are big, the histogram appears bimodal, this The lowest valley point during the time is used as the threshold point; step 303, eliminate the reflection of the image pupil area of the binarized image; step 304 determines the center of the pupil, calculates its diameter, extracts the center of the pupil, and determines the diameter of the pupil to realize the center of the pupil tracking measurement. This method is a full-field pupillometric analysis method with good real-time performance and convenient use.

Description

全视野瞳孔分析测量方法 Full field pupil analysis measurement method

技术领域 technical field

本发明属于图像检测及处理技术,特别是一种全视野瞳孔分析测量方法。The invention belongs to image detection and processing technology, in particular to a full-view pupil analysis and measurement method.

背景技术 Background technique

目前,随着科学技术和医疗水平的不断进步,经研究表明,客观定量瞳孔检测被越来越多的应用于各个不同方向。At present, with the continuous advancement of science and technology and medical level, studies have shown that objective and quantitative pupil detection is more and more applied in various directions.

在医学研究中,准分子激光屈光矫正手术(LASIK),神经眼科(例如球后视神经炎),糖尿病,神经生物学和自主神经系统功能研究,中枢神经系统疾病、脑功能、神经精神疾病(例如颅内压增高监测、脑中风病情观察、脑死亡辅助判定),重症监护和急救医学,麻醉(例如麻醉深度监测),临床药理学,眼动检测、视频眼震检测/红外眼震检测(用于眩晕、前庭器官疾病的检查)等;吸毒研究和检测,可用于吸毒界定、生理和心理成瘾(渴求)程度判定、戒毒效果评定等科研和检测领域,供禁毒执法部门、戒毒所、吸毒鉴定机构及相关科研机构使用,可大大降低检定费用,并大大节省毒品和其他药物的初筛时间;疲劳和饮酒的判定,在工业、交通、运输等部门使用瞳孔检测仪可有效评定工作者或司机的疲劳程度(疲劳驾驶)或者是否有酒后驾驶,有助于降低相关事故和酒后撞车的发生,并可提高工作场所的安全性以及增加工作效率。In medical research, excimer laser refractive surgery (LASIK), neuro-ophthalmology (such as retrobulbar optic neuritis), diabetes, neurobiology and autonomic nervous system function research, central nervous system diseases, brain function, neuropsychiatric diseases ( Such as monitoring of increased intracranial pressure, observation of cerebral apoplexy, auxiliary judgment of brain death), intensive care and emergency medicine, anesthesia (such as anesthesia depth monitoring), clinical pharmacology, eye movement detection, video nystagmus detection/infrared nystagmus detection (using Vertigo, vestibular disease inspection), etc.; drug abuse research and testing, can be used in scientific research and testing fields such as the definition of drug abuse, the determination of the degree of physical and psychological addiction (craving), and the evaluation of drug rehabilitation effects. Appraisal institutions and related scientific research institutions can greatly reduce the cost of verification, and greatly save the initial screening time of drugs and other drugs; the judgment of fatigue and drinking, the use of pupil detectors in industry, transportation, transportation and other departments can effectively evaluate workers or A driver's level of fatigue (drowsy driving) or the presence or absence of drunk driving can help reduce the occurrence of related accidents and drunk crashes, and can improve workplace safety and increase work efficiency.

国内现有的简易瞳孔测量板不能够达到实时测量,双眼瞳孔检测仪又不适宜携带,其他手持瞳孔测量器方便携带,但不能得出实时测量结果。The existing simple pupil measuring boards in China cannot achieve real-time measurement, and the binocular pupil detector is not suitable for carrying. Other hand-held pupil measuring devices are easy to carry, but cannot obtain real-time measurement results.

发明内容 Contents of the invention

本发明的目的是提供一种实时性好、使用方便的一种全视野瞳孔分析测量方法。The object of the present invention is to provide a full-field pupil analysis measurement method with good real-time performance and convenient use.

本发明的目的是由以下技术方案实现的:一种全视野瞳孔分析测量方法,包括如下步骤:The purpose of the present invention is achieved by the following technical solutions: a full-field pupil analysis measurement method, comprising the steps:

从步骤300开始,首先通过视频采集卡将眼部区域的图像数据实时采集,得到眼部区域图像;Starting from step 300, firstly, the image data of the eye area is collected in real time through a video capture card to obtain an image of the eye area;

步骤301通过均值法或最值法或公式法对300步采集到的图像进行灰度化;Step 301 grayscales the image collected in step 300 by the mean value method or the maximum value method or the formula method;

步骤302,通过直方图法对步骤301灰度化的图像进行二值化处理,当背景与目标差异大,直方图出现双峰,这时的最低谷点作为阈值点;Step 302, binarize the grayscaled image in step 301 by the histogram method, when the difference between the background and the target is large, the histogram has double peaks, and the lowest valley point at this time is used as the threshold point;

步骤303,对进行二值化处理的图像瞳孔区域反光进行消除;Step 303, eliminating the reflection in the pupil region of the image undergoing binarization processing;

步骤304确定瞳孔中心,计算其直径,提取瞳孔的圆心,以及确定瞳孔的直径,实现瞳孔中心的跟踪测量。Step 304 determines the center of the pupil, calculates its diameter, extracts the center of the pupil, and determines the diameter of the pupil, so as to realize the tracking measurement of the pupil center.

所述的确定圆心的方法可以通过三条弦的垂直平分线的交点确定。The method for determining the center of the circle can be determined by the intersection of the perpendicular bisectors of the three chords.

所述的确定圆心的方法可以通过两条弦的中点确定或通过圆上的三点,进一步利用如下公式确定:设(a,b)为圆心,R为半径,(x1,y1),(x2,y2),(x3,y3)为圆上的三点。根据下面的公式,最终可以解出a,b,R的值,进而可以确定圆心,计算出面积。The method for determining the center of the circle can be determined by the midpoints of two chords or by three points on the circle, and further determined by the following formula: Let (a, b) be the center of the circle, R be the radius, (x1, y1), ( x2, y2), (x3, y3) are three points on the circle. According to the following formula, the values of a, b, and R can be solved finally, and then the center of the circle can be determined and the area can be calculated.

(x1-a)2+(y1-b)2=R2 (x1-a) 2 +(y1-b) 2 =R 2

(x2-a)2+(y2-b)2=R2 (x2-a) 2 +(y2-b) 2 =R 2

(x3-a)2+(y3-b)2=R2(x3-a) 2 +(y3-b) 2 =R 2 .

所述的步骤301对同一幅彩色图像进行灰度化是利用公式:Y=0.3R+0.59G+0.11B进行处理。In the step 301, the gray scale of the same color image is processed by using the formula: Y=0.3R+0.59G+0.11B.

该方法采用的装置至少包括一个图像检测单元、一图像处理单元,其图像检测单元和图像处理单元之间通过电连接。The device adopted in the method at least includes an image detection unit and an image processing unit, and the image detection unit and the image processing unit are electrically connected.

所述的电连接可以是有线电连接。The electrical connection may be a wired electrical connection.

所述的电连接可以是无线电连接。Said electrical connection may be a radio connection.

所述的图像检测单元电连接有无线发射装置,无线发射装置是一视频无线发射装置。The image detection unit is electrically connected with a wireless transmitting device, and the wireless transmitting device is a video wireless transmitting device.

所述的图像处理单元电连接有无线接收装置,无线接收装置是视频接收卡。The image processing unit is electrically connected with a wireless receiving device, and the wireless receiving device is a video receiving card.

所述的图像处理单元是笔记本电脑。The image processing unit is a notebook computer.

所述的图像检测单元是CMOS摄像机或CCD摄像机。The image detection unit is a CMOS camera or a CCD camera.

本发明中采用了三点确定圆心以及扫描像素得到直径的方法。In the present invention, three points are used to determine the center of the circle and the method of scanning pixels to obtain the diameter is adopted.

本发明的工作过程及优点是:Working process of the present invention and advantage are:

采用有线视频或无线发射装置和视频接收卡进行视频信息传递,可以使摄像机构成的手持式瞳孔测量器以采集频率为30帧每秒图像信号实时传到笔记本电脑内,经笔记本电脑的图像处理软件对瞳孔图像进行快速分析,迅速得到分析结果。Using wired video or wireless transmitting device and video receiving card for video information transmission, the hand-held pupil measuring instrument composed of camera can be transmitted to the notebook computer in real time with the acquisition frequency of 30 frames per second, and the image processing software of the notebook computer Quickly analyze the pupil image and get the analysis result quickly.

由于现有的笔记本电脑具有携带方便的特点,而笔记本电脑对图像处理速度快,因此可实现实时测量和图像分析的任务。Since the existing notebook computer is easy to carry, and the notebook computer has a fast image processing speed, it can realize the tasks of real-time measurement and image analysis.

附图说明 Description of drawings

下面结合实施例附图对本发明作进一步说明。The present invention will be further described below in conjunction with the accompanying drawings of the embodiments.

图1是本发明实施1结构示意图;Fig. 1 is a structural representation of the present invention's implementation 1;

图2是本发明实施2结构示意图;Fig. 2 is the structure schematic diagram of implementing 2 of the present invention;

图3是瞳孔图像提取和快速分析和流程图;Fig. 3 is pupil image extraction and rapid analysis and flowchart;

图4是得到眼部区域图像;Fig. 4 is to obtain eye region image;

图5是二值化后的图像示意图;Fig. 5 is a schematic diagram of an image after binarization;

图6是二值化后的消除瞳孔反光图像示意图;Fig. 6 is a schematic diagram of a pupil reflection image after binarization;

图7是图像处理中膨胀的原理过程图;Fig. 7 is a principle process diagram of expansion in image processing;

图8是将瞳孔的区域转换成一个圆域(二值图像);Fig. 8 is to convert the region of the pupil into a circular domain (binary image);

图9是将二值化后得到的圆域进行边界提取(提取圆域的轮廓)。FIG. 9 is to extract the boundary of the circle domain obtained after binarization (extract the contour of the circle domain).

图中:1、图像处理单元;2、图像检测单元;3、无线电;4、无线接收装置;5、无线发射装置;6、手柄;7、开关;8、摄像头;9、照明灯;10、驱动电路;11、连接线;12、视频接口端。In the figure: 1. Image processing unit; 2. Image detection unit; 3. Radio; 4. Wireless receiving device; 5. Wireless transmitting device; 6. Handle; 7. Switch; 8. Camera; 9. Lighting lamp; 10. Driving circuit; 11. Connecting wire; 12. Video interface terminal.

具体实施方式 Detailed ways

如图1所示,这是一款利用现有带摄像头的笔记本电脑改进的实例方式,与现有的带摄像头的笔记本电脑不同之处,笔记本电脑(图像处理单元1)原装的摄像头8通过连接线11与笔记本电脑电连接;笔记本电脑包括视频接口端12。这样图像检测单元2就应包括摄像头8、驱动电路10、固定摄像头8和驱动电路10的壳体,壳体下端有手柄6,手柄6上有开关7。实际上图像检测单元2构成了一种手持式瞳孔测量仪。As shown in Figure 1, this is an example method that utilizes an existing notebook computer with a camera to improve, and is different from an existing notebook computer with a camera, the original camera 8 of the notebook computer (image processing unit 1) is connected The line 11 is electrically connected with the notebook computer; the notebook computer includes a video interface terminal 12 . In this way, the image detection unit 2 should include a camera 8, a drive circuit 10, a housing for fixing the camera 8 and the drive circuit 10, a handle 6 is arranged at the lower end of the housing, and a switch 7 is arranged on the handle 6. The image detection unit 2 actually constitutes a hand-held pupillometer.

图像检测单元2的电源直接取自笔记本电脑电源;将笔记本电脑原有的摄像头8驱动电路10移自图像检测单元2内,使图像检测单元2直接输出视频信号到笔记本电脑,笔记本电脑内再对视频信号进行处理。The power supply of the image detection unit 2 is directly taken from the notebook computer power supply; the original camera 8 drive circuit 10 of the notebook computer is moved from the image detection unit 2, so that the image detection unit 2 directly outputs the video signal to the notebook computer, and then the notebook computer video signal processing.

这种实施方式由于不需要创造性劳动,实施起来特别方便。为了实现本发明对瞳孔测量和处理的需求,笔记本电脑内有图像处理软件,通过对摄像头8提取的瞳孔图像灰度化处理,边缘提取,扫描像素对视频单帧图像进行分析,可排除或者补偿因眨眼所带来的误差导致瞳孔直径的不确定性,保证了测量结果的可信性和可重复性。Since this implementation mode does not require creative labor, it is particularly convenient to implement. In order to realize the needs of the present invention for pupil measurement and processing, there is image processing software in the notebook computer, through the gray scale processing of the pupil image extracted by the camera 8, edge extraction, and scanning pixels to analyze the single frame image of the video, it can be eliminated or compensated. The uncertainty of the pupil diameter due to the error caused by blinking ensures the reliability and repeatability of the measurement results.

为了不影响暗环境下检查,摄像头8一侧配有光源,该光源为照明灯9,光源可以是外置伸缩光源可控,包括照射单个瞳孔时确保不漏光照到另外一个瞳孔,影响测量结果;保证在暗环境下检查时不受室内外照条件的影响,确保各种环境下测试的一致性。摄像头8可以物距调整和定标校准,在不同患者的眼睛因凸凹程度不同时可按需分别调节摄像头的物距,保证始终得到个体化的最佳图像质量,并保证测量结果不受物距变动的影响。In order not to affect the inspection in a dark environment, the side of the camera 8 is equipped with a light source, which is a lighting lamp 9, and the light source can be an external retractable light source that can be controlled, including ensuring that no light leaks to another pupil when irradiating a single pupil, which will affect the measurement results. ;Ensure that it is not affected by indoor and outdoor lighting conditions when checking in a dark environment, and ensure the consistency of testing in various environments. The camera 8 can adjust and calibrate the object distance. When the eyes of different patients are different due to the degree of convexity and concaveness, the object distance of the camera can be adjusted separately according to the needs, so as to ensure that the best individual image quality is always obtained, and the measurement results are not affected by the object distance. impact of changes.

使用时首先将手持式瞳孔测量仪对准被测人眼(左右眼都可),打开开关,通过调焦按钮调节远近,在屏幕上得到清晰的被测人眼视频图像,按下开始键,录像开始,当录像结束时,按下停止键即可。When using it, first align the hand-held pupil measuring instrument with the eye of the person under test (both left and right eyes are acceptable), turn on the switch, adjust the distance through the focus button, and get a clear video image of the eye under test on the screen, press the start button, The recording starts, and when the recording ends, press the stop button.

图2是一种无线式实施方案,在这种方案中,通过在图像检测单元2(手持式瞳孔测量仪)的视频输出端连接一无线发射装置5,无线发射装置5是一视频无线发射装置,通过无线电3不断向外发送无线视频信号。而在笔记本电脑(图像处理单元1)的USB接口上电连接一无线接收装置6,无线接收装置6是视频接收卡。就可实现像检测单元2(手持式瞳孔测量仪)与笔记本电脑(图像处理单元1)远距离互动。Fig. 2 is a kind of wireless implementation scheme, in this scheme, by connecting a wireless transmitting device 5 at the video output end of image detection unit 2 (hand-held pupil measuring instrument), wireless transmitting device 5 is a video wireless transmitting device , continuously sending out wireless video signals through radio 3. On the USB interface of the notebook computer (image processing unit 1), a wireless receiving device 6 is electrically connected, and the wireless receiving device 6 is a video receiving card. The long-distance interaction between the image detection unit 2 (handheld pupil measuring instrument) and the notebook computer (image processing unit 1) can be realized.

无论是实施例1和实施例2都是利用现有的技术做了结构性的变化,使全视野瞳孔分析变得方便、准确。而(手持式瞳孔测量仪)与笔记本电脑(图像处理单元1)通过有线或无线连接则是实现这一技术的关键。Both Embodiment 1 and Embodiment 2 utilize the existing technology to make structural changes, so that the full-field pupil analysis becomes convenient and accurate. And (handheld pupil measuring instrument) and notebook computer (image processing unit 1) are the key to realize this technology through wired or wireless connection.

图3是利用图1或图2进行瞳孔图像提取和快速分析和流程图:Fig. 3 is to utilize Fig. 1 or Fig. 2 to carry out pupil image extraction and fast analysis and flowchart:

步骤300为瞳孔的提取过程开始,首先通过视频采集卡将眼部区域的图像数据实时采集,得到眼部区域图像,如图4所示。Step 300 is the start of the pupil extraction process. Firstly, the image data of the eye area is collected in real time through a video capture card to obtain an image of the eye area, as shown in FIG. 4 .

步骤301是图象的灰度化处理,通过均值法或最值法或公式法对300步采集到的图像进行灰度化,上述的三种方法对同一幅彩色图像进行灰度化,通过视觉观察对比发现利用公式法处理后的图像的效果比较理想,图像的颜色相对更均匀,视觉效果更好,因此,本软件采用了第三种方法进行灰度化即利用公式:Y=0.3R+0.59G+0.11B。Step 301 is the grayscale processing of the image. The image collected in step 300 is grayscaled by the mean value method or the maximum value method or the formula method. The above three methods grayscale the same color image. Observation and comparison found that the effect of the image processed by the formula method is ideal, the color of the image is relatively more uniform, and the visual effect is better. Therefore, this software adopts the third method for grayscale conversion, that is, using the formula: Y=0.3R+ 0.59G+0.11B.

通过直方图法对步骤301灰度化的图像进行二值化处理。直方图法是当背景与目标差异大并且面积相当时直方图出现双峰,这时的最低谷点为T,这样在进行二值化时就将最低谷点T作为阈值点,进行步骤302眼部区域的二值化处理,如图5所示。The grayscaled image in step 301 is binarized by the histogram method. The histogram method is that when the difference between the background and the target is large and the area is equal, the histogram appears double peaks. At this time, the lowest valley point is T, so when the binarization is performed, the lowest valley point T is used as the threshold point, and step 302 is performed. Binarization processing of the inner region, as shown in Figure 5.

步骤303瞳孔区域反光的消除,在实时采集眼部图像时会在瞳孔部位产生反光区域,二值化后可以清楚的看到反光区域,通过图5的白色区域所示,这些反光区域将会影响下一步的瞳孔的中心定位,因此必须将反光区域进行消除。Step 303 Elimination of reflections in the pupil area. When the eye image is collected in real time, a reflection area will be generated on the pupil. After binarization, the reflection area can be clearly seen. As shown in the white area in Figure 5, these reflection areas will affect The next step is to locate the center of the pupil, so the reflective area must be eliminated.

如图7所示,反光区域进行消除以图像处理中膨胀的原理为基础进行处理。首先膨胀的简单定义:设A为图像区域集合,B为结构元素。

Figure A20091002170400081
为B的关于原点对称的集合,
Figure A20091002170400082
表示平移z. A ⊕ B = { z | [ ( B ^ ) Z ∩ A ] ⊆ A } . B膨胀A实际上就是
Figure A20091002170400085
的位移与A至少有一个非零元素相交时B的原点为位置的集合。As shown in Figure 7, the elimination of the reflective area is based on the principle of expansion in image processing. First a simple definition of inflation: Let A be the set of image regions and B be the structural element.
Figure A20091002170400081
is a set of B that is symmetric about the origin,
Figure A20091002170400082
express translation z. A ⊕ B = { z | [ ( B ^ ) Z ∩ A ] ⊆ A } . B inflate A is actually
Figure A20091002170400085
The set of positions where the origin of B intersects with at least one nonzero element of A.

根据图7可以看出,如果想进行图像内部的空白自动的黑色填充,关键必须区分出背景的白色,和对象(瞳孔内部区域)的白色。对于这个问题解决方案如下:首先逐行扫描,找出左边是黑色,右边是白色的像素点,这样就可以将左部分的背景的白色区域排除掉,留下的区域如图6红色,蓝色标记。其次再判断每一个红色区域的像素点所在的行的右边是否还有黑色的像素点,如果有黑色的像素则说明此点是在内部,否则可以判断出此点就是蓝色区域所在的背景的白色。最后将判断出的内部区域进行自动填充。结果如下图8所示。According to Figure 7, it can be seen that if you want to automatically fill the blank inside the image with black, the key must be to distinguish the white of the background from the white of the object (the inner area of the pupil). The solution to this problem is as follows: First, scan line by line to find out the pixels that are black on the left and white on the right, so that the white area of the background on the left can be excluded, and the remaining area is shown in Figure 6. Red, blue mark. Secondly, judge whether there are black pixels on the right side of the row where each pixel in the red area is located. If there are black pixels, it means that this point is inside, otherwise it can be judged that this point is the background of the blue area. White. Finally, the determined internal area is automatically filled. The result is shown in Figure 8 below.

步骤304确定瞳孔中心,计算其直径,本发明实施例的最终的目的就是提取瞳孔的圆心,以及确定瞳孔的直径(瞳孔暂可近似看成一个圆),实现瞳孔中心的跟踪测量。根据几何的方法,确定圆心的方法比较多,可以通过三条弦的垂直平分线的交点确定,还可以通过两条弦的中点确定,也可以通过圆上的三点,进一步利用公式确定。本发明中采用了三点确定圆心以及扫描像素得到直径的方法。Step 304 determines the center of the pupil and calculates its diameter. The ultimate purpose of the embodiment of the present invention is to extract the center of the pupil and determine the diameter of the pupil (the pupil can be approximated as a circle temporarily), so as to realize the tracking measurement of the pupil center. According to geometric methods, there are many ways to determine the center of a circle. It can be determined by the intersection of the perpendicular bisectors of the three chords, the midpoint of the two chords, or the three points on the circle, and further use formulas to determine. In the present invention, three points are used to determine the center of the circle and the method of scanning pixels to obtain the diameter is adopted.

通过圆心计算直径:Calculate the diameter from the center of a circle:

设(a,b)为圆心,R为半径,(x1,y1),(x2,y2),(x3,y3)为圆上的三点。根据下面的公式,最终可以解出a,b,R的值,进而可以确定圆心,计算出面积。Let (a, b) be the center of the circle, R be the radius, and (x1, y1), (x2, y2), (x3, y3) be three points on the circle. According to the following formula, the values of a, b, and R can be solved finally, and then the center of the circle can be determined and the area can be calculated.

(x1-a)2+(y1-b)2=R2 (x1-a) 2 +(y1-b) 2 =R 2

(x2-a)2+(y2-b)2=R2 (x2-a) 2 +(y2-b) 2 =R 2

(x3-a)2+(y3-b)2=R2 (x3-a) 2 +(y3-b) 2 =R 2

根据上面的公式,如何得到圆上的三点是关键。而原始图像在经过上述的过程处理之后,已经是将瞳孔的区域转换成一个圆域(二值图像)如图8所示。因此如何得到圆域的边界上的点是关键。利用下述的方法得到圆上三点。首先,将二值化后得到的圆域进行边界提取(提取圆域的轮廓)操作如图9所示,得到圆后在圆上取三个不同的点,然后按照上述的公式计算出圆心,半径,进而确定了圆心,计算可得直径。According to the above formula, how to get the three points on the circle is the key. After the original image has been processed through the above process, the area of the pupil has been converted into a circular domain (binary image) as shown in FIG. 8 . So how to get the point on the boundary of the circle domain is the key. Use the following method to get three points on the circle. First, the boundary extraction (extraction of the contour of the circle) of the circle domain obtained after binarization is performed as shown in Figure 9. After the circle is obtained, three different points are taken on the circle, and then the center of the circle is calculated according to the above formula. Radius, and then determine the center of the circle, the diameter can be calculated.

轮廓提取确定圆上3点:Contour extraction determines 3 points on the circle:

轮廓提取的主要实现原理是,逐行扫描图像找出黑色像素的点,然后进一步判断这个点的上,下,左,右都是黑色的像素点,如果是,则将此点黑色像素附予白色,否则不改变此像素点。最后可以通过边缘上3点确定瞳孔中心点,通过在边缘上提取点的坐标可以计算出瞳孔的直径。The main implementation principle of contour extraction is to scan the image line by line to find out the point of black pixels, and then further judge that the top, bottom, left and right of this point are all black pixels, if so, attach this black pixel to White, otherwise do not change this pixel. Finally, the center point of the pupil can be determined by 3 points on the edge, and the diameter of the pupil can be calculated by extracting the coordinates of the points on the edge.

扫描像素得到直径:Scan pixels to get diameter:

在边缘提取过的图像上从图像下方开始向上逐行扫描像素,扫描到的第一个黑色的像素点记录下纵坐标以及横坐标,沿着该点记录的横坐标继续垂直向上扫描直到扫描到第二个黑色的像素点,记录该点的纵坐标,2点的纵坐标之差即为瞳孔直径,单位为像素,由此可以得到瞳孔直径。On the edge-extracted image, scan the pixels from the bottom of the image line by line, and record the ordinate and abscissa of the first black pixel scanned, and continue to scan vertically upwards along the abscissa recorded at this point until it is scanned For the second black pixel, record the ordinate of this point, the difference between the ordinates of the two points is the pupil diameter, and the unit is pixel, so the pupil diameter can be obtained.

每测量出一个值以后,由软件通过判别,如果该测量值与前5个后5个测量值相差超过10个平均值的20%,那么则判断该值为眨眼波,则排除,由此对测量结果进行补偿。After each value is measured, it is judged by the software. If the difference between the measured value and the first 5 and the last 5 measured values exceeds 20% of the average value of 10, then it is judged that the value is a blink wave, and it is excluded. The measurement results are compensated.

Claims (10)

1, full-view pupil Measurement and analysis method comprises the steps:
From step 300, at first the view data of ocular is gathered in real time by video frequency collection card, obtain the ocular image;
Step 301 is by averaging method or be worth method most or equation carries out gray processing to 300 images that collect of step;
Step 302 is carried out binary conversion treatment by histogram method to the image of step 301 gray processing, and when background and target difference are big, rectangular histogram occurs bimodal, and minimum valley point at this moment is as threshold point;
Step 303 is to carrying out reflective elimination of image pupil region of binary conversion treatment;
Step 304 is determined pupil center, calculates its diameter, extracts the center of circle of pupil, and the diameter of definite pupil, realizes the tracking measurement of pupil center.
2, full-view pupil Measurement and analysis method according to claim 1 is characterized in that: the method in described definite center of circle can be passed through the intersection point of the perpendicular bisector of three strings and determine.
3, full-view pupil Measurement and analysis method according to claim 1, it is characterized in that: the method in described definite center of circle can the mid point by two strings be determined or by 3 points on the circle, further utilize following formula to determine: to establish that (a b) is the center of circle, and R is a radius, (x1, y1), (x2, y2), (x3 y3) is 3 points on the circle; According to following formula, finally can solve a, b, the value of R, and then can determine the center of circle, calculate area:
(x1-a) 2+(y1-b) 2=R 2
(x2-a) 2+(y2-b) 2=R 2
(x3-a) 2+(y3-b) 2=R 2
4, full-view pupil Measurement and analysis method according to claim 1 is characterized in that: it is to utilize formula that described step 301 pair same width of cloth coloured image carries out gray processing: Y=0.3R+0.59G+0.11B handles.
5, full-view pupil measurement analysis device which comprises at least an image detecting element (2), a graphics processing unit (1), it is characterized in that: pass through between image detecting element (2) and the graphics processing unit (1) to be electrically connected.
6, full-view pupil measurement analysis device according to claim 5 is characterized in that: described electrical connection can be wired electrical connection or dedicated radio link.
7, full-view pupil measurement analysis device according to claim 5 is characterized in that: described image detecting element (2) is electrically connected with wireless launcher (5), and wireless launcher (5) is a video wireless launcher.
8, full-view pupil measurement analysis device according to claim 5 is characterized in that: described graphics processing unit (1) is electrically connected with radio receiver (4), and radio receiver (4) is the video reception card.
9, full-view pupil measurement analysis device according to claim 5 is characterized in that: described graphics processing unit (2) is a notebook computer.
10, full-view pupil measurement analysis device according to claim 5 is characterized in that: described image detecting element (2) is cmos camera or ccd video camera.
CNA2009100217046A 2009-03-26 2009-03-26 Full-view pupil analysis measurement method Pending CN101513342A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNA2009100217046A CN101513342A (en) 2009-03-26 2009-03-26 Full-view pupil analysis measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNA2009100217046A CN101513342A (en) 2009-03-26 2009-03-26 Full-view pupil analysis measurement method

Publications (1)

Publication Number Publication Date
CN101513342A true CN101513342A (en) 2009-08-26

Family

ID=41038077

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2009100217046A Pending CN101513342A (en) 2009-03-26 2009-03-26 Full-view pupil analysis measurement method

Country Status (1)

Country Link
CN (1) CN101513342A (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102824160A (en) * 2011-06-17 2012-12-19 北京爱生科贸有限公司 Eyeball movement monitoring method and equipment thereof
CN104739367A (en) * 2015-03-14 2015-07-01 中国科学院苏州生物医学工程技术研究所 Binocular pupil light synthetic detection system
CN105125169A (en) * 2014-05-30 2015-12-09 中国人民解放军第三军医大学野战外科研究所 Lung blast injury detecting instrument and method for detecting lung blast injury image
CN105760848A (en) * 2016-03-04 2016-07-13 蒋志平 Pupil positioning method based on annular template convolution
CN108562548A (en) * 2018-01-26 2018-09-21 北京秒康智能科技有限公司 The color identification method and its system of intelligent urine examination closestool
CN110495895A (en) * 2019-08-26 2019-11-26 重庆大学 A fatigue detection method and system based on eye movement tracking
CN110751093A (en) * 2019-10-20 2020-02-04 赵俊杰 Method for automatically identifying pupil and calculating pupil diameter
CN110786823A (en) * 2019-12-12 2020-02-14 新疆医科大学第三附属医院 A kind of infrared pupil identification measurement method, device, storage medium and terminal equipment
CN111633337A (en) * 2020-05-25 2020-09-08 西咸新区大熊星座智能科技有限公司 Reflection eliminating method and device for laser welding seam measurement
WO2021093744A1 (en) * 2019-11-11 2021-05-20 中国科学院深圳先进技术研究院 Method and apparatus for measuring diameter of pupil, and computer-readable storage medium
CN114652480A (en) * 2022-03-10 2022-06-24 深圳市眼科医院 Method for constructing optic nerve injury model of rana grahami, model and application method
CN117056872A (en) * 2023-07-25 2023-11-14 中国人民解放军空军军医大学 Advanced brain fatigue prediction method

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102824160A (en) * 2011-06-17 2012-12-19 北京爱生科贸有限公司 Eyeball movement monitoring method and equipment thereof
CN102824160B (en) * 2011-06-17 2015-01-21 北京爱生科贸有限公司 Eyeball movement monitoring method and equipment thereof
CN105125169A (en) * 2014-05-30 2015-12-09 中国人民解放军第三军医大学野战外科研究所 Lung blast injury detecting instrument and method for detecting lung blast injury image
CN104739367A (en) * 2015-03-14 2015-07-01 中国科学院苏州生物医学工程技术研究所 Binocular pupil light synthetic detection system
CN105760848A (en) * 2016-03-04 2016-07-13 蒋志平 Pupil positioning method based on annular template convolution
CN105760848B (en) * 2016-03-04 2019-05-14 蒋志平 A kind of pupil positioning method based on annular mask convolution
CN108562548A (en) * 2018-01-26 2018-09-21 北京秒康智能科技有限公司 The color identification method and its system of intelligent urine examination closestool
CN110495895A (en) * 2019-08-26 2019-11-26 重庆大学 A fatigue detection method and system based on eye movement tracking
CN110751093A (en) * 2019-10-20 2020-02-04 赵俊杰 Method for automatically identifying pupil and calculating pupil diameter
CN110751093B (en) * 2019-10-20 2023-10-27 赵俊杰 Method for automatically identifying pupil and calculating pupil diameter
WO2021093744A1 (en) * 2019-11-11 2021-05-20 中国科学院深圳先进技术研究院 Method and apparatus for measuring diameter of pupil, and computer-readable storage medium
CN110786823A (en) * 2019-12-12 2020-02-14 新疆医科大学第三附属医院 A kind of infrared pupil identification measurement method, device, storage medium and terminal equipment
CN111633337A (en) * 2020-05-25 2020-09-08 西咸新区大熊星座智能科技有限公司 Reflection eliminating method and device for laser welding seam measurement
CN114652480A (en) * 2022-03-10 2022-06-24 深圳市眼科医院 Method for constructing optic nerve injury model of rana grahami, model and application method
CN117056872A (en) * 2023-07-25 2023-11-14 中国人民解放军空军军医大学 Advanced brain fatigue prediction method
CN117056872B (en) * 2023-07-25 2025-01-24 中国人民解放军空军军医大学 A method for predicting brain fatigue in advance

Similar Documents

Publication Publication Date Title
CN101513342A (en) Full-view pupil analysis measurement method
JP6072798B2 (en) System and method for documenting and recording pupil red reflex examination and corneal light reflex screening of eyes in infants and children
EP1587414B1 (en) A pupilometer
CN103919557B (en) Nystagmus parameter feature acquisition method and device for diagnosing benign paroxysmal positional vertigo
CN105433899A (en) Pupil detection device
EP3127471B1 (en) Automatic vision testing device without human monitoring
CN109662685B (en) Eyeball saliency meter and measuring method thereof
CN110710978A (en) A Multimodal Immersive Synchronous Acquisition System Based on Eye Tracking-Brain Function Activity Detection
KR20150036395A (en) Strabismus Detection
CN101190120A (en) Binocular pupil detecting system
CN110507281A (en) An eyepiece type pupil light reflex automatic detection device
CN106725295A (en) A kind of miniature medical examination equipment, device and using method thereof
CN110786823A (en) A kind of infrared pupil identification measurement method, device, storage medium and terminal equipment
CN111297347A (en) A method and apparatus for generating photoplethysmography signals
CN117617889A (en) A strabismus degree detection system based on covering test
CN211094132U (en) Multi-mode immersive synchronous acquisition system based on eye movement tracking-brain function activities
CN101108120A (en) Eye movement test analysis method
CN201375508Y (en) Full-field pupil analysis measuring device
CN113011286B (en) Strabismus discrimination method and system based on video-based deep neural network regression model
EP3925518A1 (en) Detecting and tracking macular degeneration
CN105380590A (en) Device with eye position detection function and achieving method thereof
TWI511095B (en) A fatigue detecting method with a characteristic of combining an eyeball image and brain waves
CN211432840U (en) Eyepiece formula pupil light reflex automated inspection equipment
CN117530654A (en) A real-time binocular pupil inspection system and detection method
CN211723119U (en) A kind of infrared pupil recognition measuring device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Open date: 20090826