基于运动视觉诱发电位的视敏度检测方法Visual acuity detection method based on motor visual evoked potential
技术领域Technical field
本发明涉及脑-机接口与眼视光学检查技术领域,具体涉及基于运动视觉诱发电位的视敏度检测方法。The invention relates to the technical field of brain-computer interface and optometry inspection, in particular to a visual acuity detection method based on motion visual evoked potential.
背景技术Background technique
世界卫生组织的数据显示,全球有2.53亿人患有视力障碍,并且由于人口增长和老龄化,这一数字可能会增加两倍,而视敏度检查是眼科学最基本的研究任务之一。常见的视敏度测试比较依赖于被试的主观感受和配合,如字母视标视敏度表和翻滚E视敏度表等。对于有沟通问题的人,如语前儿童或幼儿,不适合利用这类方法进行视敏度评估;成年人也可能受到环境因素的影响产生误判;在法医学检验中不能排除被鉴定者的谎报与欺骗。According to data from the World Health Organization, 253 million people worldwide suffer from visual impairment, and due to population growth and aging, this number may triple. Visual acuity testing is one of the most basic research tasks in ophthalmology. Common visual acuity tests rely on subjects’ subjective feelings and cooperation, such as letter visual acuity scale and rolling E visual acuity scale. For people with communication problems, such as pre-lingual children or toddlers, this type of method is not suitable for visual acuity assessment; adults may also be affected by environmental factors and cause misjudgments; false reports by the examinee cannot be ruled out in forensic examinations And deception.
脑-机接口(BCI)和脑电图(EEG),尤其是头皮脑电图,是一种更客观、更直接地评估视觉功能的新方法。近年来,不少学者利用扫描视觉诱发电位(sVEP)、稳态视觉诱发电位(SSVEP)、模式视觉诱发电位(PVEP)等方法完成了对视敏度以及对比敏感度的检测。但目前的研究主要着眼于通过改变一些相关刺激参数以确定视觉功能阈值以及观察脑电图反应。这些研究没有将客观的脑电图测试与主观的心理物理测试相结合,刺激参数与视觉功能状态之间也没有确切的相关性。此外,传统的棋盘格和光栅的SSVEP或PVEP范式,更可能导致视觉疲劳,从而导致信号质量下降。Brain-computer interface (BCI) and electroencephalogram (EEG), especially scalp electroencephalogram, are new methods for more objective and direct assessment of visual function. In recent years, many scholars have used scanning visual evoked potential (sVEP), steady-state visual evoked potential (SSVEP), pattern visual evoked potential (PVEP) and other methods to complete the detection of visual acuity and contrast sensitivity. However, the current research mainly focuses on determining the threshold of visual function and observing the EEG response by changing some relevant stimulus parameters. These studies did not combine objective EEG tests with subjective psychophysical tests, and there was no exact correlation between stimulation parameters and visual function status. In addition, the traditional checkerboard and grating SSVEP or PVEP paradigm is more likely to cause visual fatigue, which leads to reduced signal quality.
发明内容Summary of the invention
为了克服上述现有技术的缺点,本发明的目的在于提供基于运动视觉诱发电位的视敏度检测方法,通过设计一种用于视敏度检查的同心圆环稳态运动视觉诱发电位(SSMVEP)范式,并提出相应的阈值判定标准,得出脑电客观视敏度值,为眼科检查中视敏度的检测提供一种客观且定量的测量方法。In order to overcome the above-mentioned shortcomings of the prior art, the purpose of the present invention is to provide a visual acuity detection method based on motion visual evoked potential, by designing a concentric ring steady-state motion visual evoked potential (SSMVEP) for visual acuity inspection Paradigm, and put forward the corresponding threshold criteria to obtain the EEG objective visual acuity value, which provides an objective and quantitative measurement method for the detection of visual acuity in eye examinations.
为了达到上述目的,本发明采取的技术方案为:In order to achieve the above objectives, the technical solutions adopted by the present invention are:
基于运动视觉诱发电位的视敏度检测方法,包括以下步骤:The visual acuity detection method based on motion visual evoked potential includes the following steps:
1)同心圆环刺激范式:采用脑-机接口中的SSMVEP方法,通过 MATLAB使用Psychophysics Toolbox编程绘制范式图案纹理的周期性收缩扩张运动,能够稳定刺激诱发SSMVEP;运动刺激目标表现为明暗相间同心圆环,明亮区域和暗区域间距相等;在刺激呈现过程中,范式整体亮度与整体尺寸保持不变;圆环相位按正弦方式进行调制,在收缩和扩张两个方向上形成周期往复振荡运动,收缩和扩张运动的翻转频率定义为刺激频率,运动方向在一个周期内改变两次;1) Concentric circle stimulation paradigm: Using the SSMVEP method in the brain-computer interface, the periodic contraction and expansion movement of the paradigm pattern texture is drawn by MATLAB using Psychophysics Toolbox programming, which can stably induce SSMVEP; the movement stimulus target appears as light and dark concentric circles In the ring, the distance between the bright area and the dark area is equal; during the stimulus presentation process, the overall brightness and overall size of the paradigm remain unchanged; the phase of the ring is modulated in a sinusoidal manner, forming periodic reciprocating oscillations in both directions of contraction and expansion, and contraction The flip frequency of the expansion movement is defined as the stimulation frequency, and the movement direction changes twice in one cycle;
2)刺激范式空间频率扫描梯度:对应于标准对数视敏度表,设计一系列同心圆环范式,作为SSMVEP视敏度检测的视标;对应于标准对数视敏度表的1.1 logMAR与-0.1 logMAR视标,通过调整范式的空间频率,使同心圆环的黑白间隔对人眼视角与视敏度表中字母“E”的黑白间隔对人眼视角相等;空间频率从1.1 logMAR以每秒0.1 logMAR向-0.1 logMAR线性递减,总扫描时长为12秒;2) Spatial frequency sweep gradient of stimulus paradigm: Corresponding to the standard logarithmic visual acuity table, a series of concentric circle paradigms are designed as the visual standard for SSMVEP visual acuity detection; corresponding to 1.1 logMAR and the standard logarithmic visual acuity table -0.1 logMAR optotype, by adjusting the spatial frequency of the paradigm, the black-and-white interval of the concentric circles is equal to the human eye angle of view with the black-and-white interval of the letter "E" in the visual acuity table; the spatial frequency is from 1.1 logMAR to each The second 0.1 logMAR decreases linearly to -0.1 logMAR, and the total scan time is 12 seconds;
3)脑-机接口平台搭建:实验前将参考电极放置于受试者左耳垂,地电极放置于受试者前额,测量电极布置在头部枕区;电极与脑电采集模块连接,经放大、滤波、数模转换后输出脑电信号到计算机,进一步数据处理;3) Construction of brain-computer interface platform: Before the experiment, the reference electrode was placed on the subject’s left earlobe, the ground electrode was placed on the subject’s forehead, and the measurement electrode was placed on the occipital area of the head; the electrode was connected to the EEG acquisition module, and enlarged , After filtering, digital-to-analog conversion, output EEG signals to the computer for further data processing;
4)刺激范式呈现:计算机通过显示器扩展使高刷新率显示器上呈现刺激范式;在实验过程中,使用者需要根据屏幕提示注视刺激范式;实验过程中,被试者需要坐在实验室中,每次进行单眼测试;通过脑电采集设备采集使用者注视范式时产生的脑电信号,经放大、滤波与A/D转换后,将处理后的脑电信号输入计算机利用典型相关分析(CCA)进行特征提取;4) Presentation of the stimulus paradigm: The computer expands the display to present the stimulus paradigm on the high refresh rate display; during the experiment, the user needs to follow the screen prompts to watch the stimulus paradigm; during the experiment, the subject needs to sit in the laboratory every time Perform a monocular test once; collect the EEG signals generated by the user's gaze paradigm through EEG acquisition equipment, and after amplifying, filtering and A/D conversion, the processed EEG signals are input to the computer and performed by canonical correlation analysis (CCA) Feature extraction
5)视敏度阈值判定:通过窗长2秒、滑移时长1秒的矩形窗对12秒的SSMVEP信号进行分析,得出典型相关系数随时间变化曲线及对应的信噪比SNR,从而进行范式响应值的判定;采用SNR值对CCA谱进行响应值判定,当SNR大于2.0时,说明视觉诱发电位显著,响应值为1;信噪比不高于2.0时,视觉诱发电位不显著,响应值为0;SSMVEP视敏度值定义为最后脑电响应值为1时对应的刺激范式空间频率值logMAR;5) Visual acuity threshold judgment: Analyze the 12-second SSMVEP signal through a rectangular window with a window length of 2 seconds and a slip duration of 1 second, and obtain the typical correlation coefficient over time curve and the corresponding signal-to-noise ratio SNR, and then Paradigm response value determination; use SNR value to determine the response value of the CCA spectrum. When the SNR is greater than 2.0, the visual evoked potential is significant, and the response value is 1; when the signal-to-noise ratio is not higher than 2.0, the visual evoked potential is not significant, and the response Value is 0; SSMVEP visual acuity value is defined as the spatial frequency value logMAR corresponding to the stimulation paradigm when the final EEG response value is 1;
6)视敏度检测结果反馈:完成所有刺激范式刺激后,使用者的客观脑电视敏度值通过屏幕输出。6) Visual acuity detection result feedback: After completing all stimulation paradigms, the user's objective brain television acuity value is output on the screen.
所述的步骤5)中矩形窗取窗长最中心时间所对应的空间频率,并得出其所对应的CCA值与SNR值,再根据阈值判定算法进行判定。In the step 5), the rectangular window takes the spatial frequency corresponding to the most central time of the window length, and obtains the corresponding CCA value and SNR value, and then determines according to the threshold determination algorithm.
所述的步骤5)中实验得到CCA值与SNR值随范式空间频率logMAR值1.0-0.0的变化趋势,其中CCA值与SNR值都随logMAR值而变小;采用SNR值对CCA频谱进行响应值判定,信噪比定义为目标刺激频率8Hz相关系数的平方值与CCA谱上10个相邻点相关系数平方的平均值之比;当SNR大于2.0时,说明视觉诱发电位显著,响应值为1;信噪比不高于2.0时,视觉诱发电位不显著,响应值为0;SSMVEP视敏度值定义为最后脑电响应值为1时对应的刺激范式空间频率值logMAR。In the step 5), the experiment obtained the trend of CCA value and SNR value with the logMAR value 1.0-0.0 of the paradigm space frequency, where the CCA value and SNR value both decrease with the logMAR value; the SNR value is used to respond to the CCA spectrum It is judged that the signal-to-noise ratio is defined as the ratio of the square value of the correlation coefficient of the target stimulation frequency 8Hz to the average value of the square value of the correlation coefficient of 10 adjacent points on the CCA spectrum; when the SNR is greater than 2.0, it means that the visual evoked potential is significant, and the response value is 1 ; When the signal-to-noise ratio is not higher than 2.0, the visual evoked potential is not significant, and the response value is 0; the SSMVEP visual acuity value is defined as the stimulus paradigm spatial frequency value logMAR corresponding to the final EEG response value of 1.
本发明还提供了一种基于运动视觉诱发电位的视敏度检测系统,包括:The present invention also provides a visual acuity detection system based on motion visual evoked potential, including:
脑-机接口平台,包括参考电极、地电极和测量电极,各电极与脑电采集模块连接,经放大、滤波、数模转换后输出脑电信号到计算机;The brain-computer interface platform includes a reference electrode, a ground electrode and a measurement electrode. Each electrode is connected to the EEG acquisition module, and the EEG signal is output to the computer after amplification, filtering, and digital-to-analog conversion;
计算机,对所述的脑电信号做进一步数据处理,包括:基于所述的同心圆环刺激范式,刺激范式空间频率扫描梯度、刺激范式呈现,并进行视敏度阈值(阈值可预设)判定;The computer performs further data processing on the EEG signal, including: based on the concentric circle stimulation paradigm, the stimulation paradigm spatial frequency sweep gradient, the stimulation paradigm presentation, and the visual acuity threshold (the threshold can be preset) is determined ;
显示器,连接所述计算机,呈现前述的范式;并可根据需要,在完成所有刺激范式刺激后,使用者的客观脑电视敏度值通过屏幕输出,该屏幕可与呈现范式的显示器为同一显示器,也可为另一单独连接的显示器。The display is connected to the computer to present the aforementioned paradigm; and according to needs, after all stimulation paradigms are completed, the user’s objective brain television acuity value is output through the screen, which can be the same display as the display that presents the paradigm, It can also be another separately connected display.
本发明的有益效果为:The beneficial effects of the present invention are:
本发明提出了基于运动视觉诱发电位的视敏度检测方法,操作简单快捷,解决了传统主观视敏度测试不够客观,不适用于幼儿、语前儿童、沟通困难患者以及法医鉴定的问题,显示了如下优越性:The invention proposes a visual acuity detection method based on motor visual evoked potentials, which is simple and quick to operate, and solves the problem that the traditional subjective visual acuity test is not objective enough and is not suitable for infants, pre-lingual children, patients with communication difficulties and forensic identification. The following advantages:
(1)本发明设计基于同心圆环SSMVEP范式,图案整个运动过程中保持亮度恒定,降低了使用者的视觉疲劳,提升了脑-机接口的交互性能。(1) The design of the present invention is based on the concentric ring SSMVEP paradigm, the pattern keeps the brightness constant throughout the movement process, reduces the user's visual fatigue, and improves the interaction performance of the brain-computer interface.
(2)本发明提出视敏度阈值判定方法,分析方法准确率高、鲁棒性好,测试结果客观、定量,灵敏快速。(2) The present invention proposes a visual acuity threshold determination method. The analysis method has high accuracy, good robustness, objective, quantitative test results, and rapid sensitivity.
附图说明Description of the drawings
图1为本发明的同心圆环范式收缩-扩张运动过程。Figure 1 shows the contraction-expansion movement process of the concentric ring paradigm of the present invention.
图2为刺激范式空间频率扫描梯度示意图。Figure 2 is a schematic diagram of the spatial frequency sweep gradient of the stimulation paradigm.
图3为典型相关系数值与信噪比值随范式扫描梯度变化示意图。Fig. 3 is a schematic diagram of the variation of the typical correlation coefficient value and the signal-to-noise ratio with the gradient of the normal scan.
图4为某被试(主观性视敏度0.4 logMAR)的CCA谱图。Figure 4 shows the CCA spectrum of a subject (subjective visual acuity 0.4 logMAR).
图5为客观脑电视敏度检测值与主观视敏度值的比较散点图。Figure 5 is a scatter diagram of the comparison between the objective brain television acuity detection value and the subjective visual acuity value.
具体实施方式Detailed ways
下面结合附图对本发明作详细说明。The present invention will be described in detail below with reference to the drawings.
基于运动视觉诱发电位的视敏度检测方法,包括以下步骤:The visual acuity detection method based on motion visual evoked potential includes the following steps:
1)同心圆环刺激范式:采用脑-机接口中最有实用意义的SSMVEP方法,通过MATLAB使用Psychophysics Toolbox编程绘制范式图案纹理的周期性收缩扩张运动,能够稳定刺激诱发SSMVEP;参照图1,运动刺激目标表现为明暗相间同心圆环,明亮区域和暗区域间距相等;在刺激呈现过程中,范式整体亮度与整体尺寸保持不变;圆环相位按正弦方式进行收缩-扩张运动,其中:1) Concentric ring stimulation paradigm: Using the most practical SSMVEP method in the brain-computer interface, the periodic contraction and expansion movement of the paradigm pattern texture is drawn through MATLAB using Psychophysics Toolbox programming, which can stably stimulate SSMVEP; refer to Figure 1, Movement The stimulus target is expressed as a light and dark concentric ring with the same distance between the bright area and the dark area; during the stimulus presentation process, the overall brightness and overall size of the paradigm remain unchanged; the phase of the ring undergoes a sinusoidal contraction-expansion movement, where:
式中:Where:
f
c——运动频率,即同心圆环收缩-扩张一次所需时间的倒数;
f c ——Motion frequency, that is, the reciprocal of the time required for concentric rings to contract-expand once;
φ(t)为棋盘格收缩-扩张的相位值函数;φ(t) is the phase value function of checkerboard contraction-expansion;
通过改变相位函数φ(t)由0到π时,同心圆环收缩,相位函数φ(t)由π到0时,同心圆环扩张,参照图1,在一个周期中,发生两次运动方向的改变,运动方向改变的频率为运动翻转频率f,为运动频率f
c的2倍;由于SSMVEP主要来源于方向改变激发的大脑活动,能量主要集中在运动反转频率上,因此,采用运动反转频率作为视觉刺激的基频;
By changing the phase function φ(t) from 0 to π, the concentric rings contract, and when the phase function φ(t) changes from π to 0, the concentric rings expand. Referring to Figure 1, in one cycle, two movement directions occur The frequency of the movement direction change is the movement reversal frequency f, which is twice the movement frequency f c ; since SSMVEP mainly comes from the brain activity stimulated by the change of direction, the energy is mainly concentrated on the movement reversal frequency. The frequency is used as the fundamental frequency of visual stimulation;
2)刺激范式空间频率扫描梯度:对应于标准对数视敏度表,设计一系列同心圆环范式,作为SSMVEP视敏度检测的视标;通过调整屏幕距离与范式尺寸,使范式整体对人眼视角为2°;对应于标准对数视敏度表的1.1 logMAR与-0.1 logMAR视标,通过调整范式的空间频率,使同心圆环的黑白间隔对人眼视角与视敏度表中字母“E”的黑白间隔对人眼视角相等;参照图2,空间频率从1.1 logMAR以每秒0.1 logMAR向-0.1 logMAR线性递减,总扫描时长为12秒;2) Spatial frequency sweep gradient of stimulus paradigm: Corresponding to the standard logarithmic visual acuity scale, a series of concentric circle paradigms are designed to serve as the visual standard for SSMVEP visual acuity detection; the paradigm is overall human-oriented by adjusting the screen distance and paradigm size The eye angle of view is 2°; corresponding to 1.1 logMAR and -0.1 logMAR optotypes of the standard logarithmic visual acuity table, by adjusting the spatial frequency of the paradigm, the black-and-white interval of the concentric circles is relative to the human eye’s angle of view and the letters in the visual acuity table. The black-and-white interval of "E" is equal to the viewing angle of the human eye; referring to Figure 2, the spatial frequency decreases linearly from 1.1 logMAR to -0.1 logMAR per second, and the total scanning time is 12 seconds;
3)脑-机接口平台搭建:实验前按照10/20系统法布置电极,将参考电极放置于受试者左耳垂A1,地电极放置于受试者前额Fpz,六块测量电极布置在头部枕区(PO3,PO4,POz,O1,O2,Oz)。给各个测量电极注入导电膏,保证电极与头皮的良好接触;电极与脑电采集模块连接,经放大、滤波、数模转换后输出脑电信号到计算机,进一步数据处理;3) Brain-computer interface platform construction: Before the experiment, the electrodes were arranged according to the 10/20 system method. The reference electrode was placed on the subject's left earlobe A1, the ground electrode was placed on the subject's forehead Fpz, and six measuring electrodes were placed on the head Occipital area (PO3, PO4, POz, O1, O2, Oz). Inject conductive paste into each measuring electrode to ensure good contact between the electrode and the scalp; the electrode is connected to the EEG acquisition module, and the EEG signal is output to the computer after amplification, filtering, and digital-to-analog conversion for further data processing;
4)刺激范式呈现:计算机通过显示器扩展使高刷新率显示器上呈现刺激范式;在实验过程中,使用者需要根据屏幕提示注视刺激范式;实验过程中,被试者需要坐在安静的不被干扰的实验室中,每次进行单眼测试;通过脑电采集设备采集使用者注视范式时产生的脑电信号,经放大、滤波与A/D转换后,将处理后的脑电信号输入计算机利用典型相关分析(CCA)进行特征提取;4) Presentation of stimulus paradigm: The computer expands the display to present the stimulus paradigm on the high refresh rate display; during the experiment, the user needs to follow the screen prompts to watch the stimulus paradigm; during the experiment, the subject needs to sit quietly and not be disturbed In the laboratory, the monocular test is performed every time; the EEG signal generated by the user’s gaze paradigm is collected through EEG acquisition equipment, and after amplification, filtering and A/D conversion, the processed EEG signal is input into the computer using typical Correlation analysis (CCA) for feature extraction;
5)视敏度阈值判定:通过矩形窗(窗长2秒,滑移时长1秒)对12秒的SSMVEP信号进行截断并进行CCA分析,参照图3,取窗长2秒中间的刺激范式空间频率为横坐标,得出典型相关系数随空间频率变化曲线及对应的信噪比SNR,从而进行范式响应值的判定;实验得到CCA值与SNR值随范式空间频率logMAR值(1.0-0.0)的变化趋势,其中CCA值与SNR值都随logMAR值而变小;采用SNR值对CCA频谱进行响应值判定,信噪比定义为目标刺激频率8Hz相关系数的平方值与CCA谱上10个相邻点相关系数平方的平均值之比;当SNR大于2.0时,说明视觉诱发电位显著,响应值为1;信噪比不高于2.0时,视觉诱发电位不显著,响应值为0;SSMVEP视敏度值定义为最后脑电响应值为1时对应的刺激范式空间频率值logMAR;5) Visual acuity threshold determination: The 12-second SSMVEP signal is truncated and analyzed by CCA through a rectangular window (window length of 2 seconds, sliding duration of 1 second), refer to Figure 3, take the stimulus paradigm space in the middle of the window length of 2 seconds The frequency is the abscissa, and the curve of the typical correlation coefficient with the spatial frequency and the corresponding SNR are obtained, so as to determine the paradigm response value; the experiment obtained the CCA value and the SNR value with the paradigm spatial frequency logMAR value (1.0-0.0) Change trend, in which the CCA value and SNR value decrease with the logMAR value; the SNR value is used to determine the response value of the CCA spectrum, and the signal-to-noise ratio is defined as the square value of the target stimulation frequency 8Hz correlation coefficient and 10 adjacent ones on the CCA spectrum The ratio of the average value of the square of the point correlation coefficient; when the SNR is greater than 2.0, the visual evoked potential is significant, and the response value is 1; when the signal-to-noise ratio is not higher than 2.0, the visual evoked potential is not significant, and the response value is 0; SSMVEP visual sensitivity The degree value is defined as the spatial frequency value logMAR of the stimulation paradigm when the final EEG response value is 1;
6)视敏度检测结果反馈:完成所有刺激范式刺激后,使用者的客观脑电视敏度值通过屏幕输出。6) Visual acuity test result feedback: After completing all stimulation paradigms, the user's objective brain television acuity value is output on the screen.
下面结合实施例对本发明进行说明。The present invention will be described below in conjunction with embodiments.
对一名被试者(其主观视敏度为0.4 logMAR)进行了上述实验,按照上述步骤3)对被试安放电极并搭建脑机接口平台,按照上述步骤4)进行范式呈现与脑电信号采集;按照上述步骤5)进行脑电信号分析与客观脑电视敏度阈值判定,参照图4,得到范式响应值:1,1,1,1,1,1,1,0,0,0,0,得到SSMVEP视敏度检测结果: 0.4 logMAR。The above experiment was carried out on a subject (with a subjective visual acuity of 0.4 logMAR). The electrodes were placed on the subject and the brain-computer interface platform was built according to the above step 3), and the paradigm presentation and EEG signals were performed according to the above step 4) Acquisition; according to the above step 5) for EEG signal analysis and objective brain TV acuity threshold determination, refer to Figure 4, get the paradigm response value: 1,1,1,1,1,1,0,0,0, 0, get the SSMVEP visual acuity test result: 0.4 logMAR.
对样本中11名使用者的脑电数据进行分析判断,并将其客观检测结果与主观测试结果比较,参照图5,绘制散点图并计算线性相关性。可以证明客观视敏度检测结果与主观测试结果显著相关(P<0.001)。The EEG data of 11 users in the sample was analyzed and judged, and the objective test results were compared with the subjective test results. Refer to Figure 5, draw a scatter plot and calculate the linear correlation. It can be proved that the objective visual acuity test results are significantly related to the subjective test results (P<0.001).
本发明能够客观、定量地检测使用者的视敏度,与主观心理物理学检查之间建立很好的相关关系,实现了对幼儿、语前儿童、法医被鉴定人等特殊使用者视敏度的客观检测手段。基于运动视觉诱发电位的视敏度检测方法方便快捷,客观定量,具有良好的实用化前景。The present invention can objectively and quantitatively detect the visual acuity of users, establishes a good correlation with subjective psychophysical examination, and realizes the visual acuity of special users such as infants, pre-linguistic children, and forensic examinees. Objective testing methods. The visual acuity detection method based on motion visual evoked potential is convenient, fast, objective and quantitative, and has a good practical prospect.