[go: up one dir, main page]

CN106955092B - Method and equipment for measuring pulse distribution - Google Patents

Method and equipment for measuring pulse distribution Download PDF

Info

Publication number
CN106955092B
CN106955092B CN201710211395.3A CN201710211395A CN106955092B CN 106955092 B CN106955092 B CN 106955092B CN 201710211395 A CN201710211395 A CN 201710211395A CN 106955092 B CN106955092 B CN 106955092B
Authority
CN
China
Prior art keywords
pulse
speckle
distribution information
arterial
optical path
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
Application number
CN201710211395.3A
Other languages
Chinese (zh)
Other versions
CN106955092A (en
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.)
Beijing Information Science and Technology University
Original Assignee
Beijing Information Science and Technology University
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 Beijing Information Science and Technology University filed Critical Beijing Information Science and Technology University
Priority to CN201710211395.3A priority Critical patent/CN106955092B/en
Publication of CN106955092A publication Critical patent/CN106955092A/en
Application granted granted Critical
Publication of CN106955092B publication Critical patent/CN106955092B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4854Diagnosis based on concepts of alternative medicine, e.g. homeopathy or non-orthodox

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Alternative & Traditional Medicine (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

本申请公开了一种脉搏分布的测量方法和设备,所述方法包括:构建散斑干涉光路,其中,所述散斑干涉光路中包含被测生物体的脉搏跳动的动脉部位;根据所述散斑干涉光路,确定所述动脉部位的脉搏分布信息,其中,所述脉搏分布信息包括:脉搏空间分布信息和/或脉搏时间分布信息。根据本申请实施例的方法和设备,能够实现对脉搏分布的测量,为临床诊断提供有效客观的脉搏数据。

Figure 201710211395

The present application discloses a method and device for measuring pulse distribution. The method includes: constructing a speckle interference optical path, wherein the speckle interference optical path includes an arterial part of the pulse beating of a measured organism; The speckle interferes with the optical path to determine the pulse distribution information of the arterial part, wherein the pulse distribution information includes: pulse spatial distribution information and/or pulse time distribution information. According to the method and device of the embodiments of the present application, the measurement of pulse distribution can be realized, and effective and objective pulse data can be provided for clinical diagnosis.

Figure 201710211395

Description

一种脉搏分布的测量方法和设备Method and device for measuring pulse distribution

技术领域technical field

本申请涉及全场光检测技术领域,尤其涉及一种脉搏分布的测量方法和设备。The present application relates to the technical field of full-field light detection, and in particular, to a method and device for measuring pulse distribution.

背景技术Background technique

中医学中的把脉又称为切脉,是指中医师通过按压患者的动脉来感知脉搏跳动,辨识脉象,进而根据脉象来诊断疾病。In TCM, pulse-taking, also known as pulse-cutting, means that TCM physicians sense the pulse beat by pressing the patient's arteries, identify the pulse, and then diagnose diseases based on the pulse.

中医师辨识脉象,常从脉位、脉数、脉形、脉势四个方面进行分析。脉位指脉搏跳动显现的部位和长度,脉数指脉搏跳动的次数和节律,脉形指脉搏跳动的宽度等形态,脉势指脉搏跳动的强弱、流畅等趋势。脉象是对人体机能的反映,不同的脉象都有所对应的病症范围。Traditional Chinese medicine practitioners often analyze pulse status from four aspects: pulse position, pulse number, pulse shape, and pulse potential. Pulse position refers to the position and length of the pulse, pulse number refers to the frequency and rhythm of the pulse, pulse shape refers to the width and other forms of the pulse, and pulse potential refers to the strength and smoothness of the pulse. Pulse is a reflection of human body functions, and different pulses have corresponding disease ranges.

人体脉搏是一种振幅微小的不稳定振动。现有技术中的脉搏测量仪,例如:电子式脉搏测量仪或光电式脉搏测量仪,只能对脉搏的振动频率信息(即脉数)进行测量,无法测量除脉搏的振动频率信息之外的其他脉搏信息。The human pulse is an unstable vibration with a small amplitude. Pulse measuring instruments in the prior art, such as electronic pulse measuring instruments or photoelectric pulse measuring instruments, can only measure the vibration frequency information of the pulse (that is, the pulse number), and cannot measure other than the vibration frequency information of the pulse. Additional pulse information.

发明内容SUMMARY OF THE INVENTION

有鉴于此,本申请实施例提供一种脉搏分布的测量方法和设备,用于解决现有技术中无法测量除脉搏的振动频率信息之外的其他脉搏信息的技术问题,为临床诊断提供客观的脉搏分布数据。In view of this, embodiments of the present application provide a method and device for measuring pulse distribution, which are used to solve the technical problem in the prior art that other pulse information other than the vibration frequency information of the pulse cannot be measured, and provide an objective for clinical diagnosis. Pulse distribution data.

本申请实施例提供一种脉搏分布的测量方法,包括:Embodiments of the present application provide a method for measuring pulse distribution, including:

构建散斑干涉光路,其中,所述散斑干涉光路中包含被测生物体的脉搏跳动的动脉部位;constructing a speckle interference light path, wherein the speckle interference light path includes an arterial part where the pulse of the tested organism is beating;

根据所述散斑干涉光路,确定所述动脉部位的脉搏分布信息,其中,所述脉搏分布信息包括:脉搏空间分布信息和/或脉搏时间分布信息。According to the speckle interference optical path, the pulse distribution information of the arterial part is determined, wherein the pulse distribution information includes: pulse spatial distribution information and/or pulse time distribution information.

可选地,根据所述散斑干涉光路,确定所述动脉部位的脉搏空间分布信息,包括:Optionally, determining the pulse spatial distribution information of the arterial site according to the speckle interference optical path, including:

在第一时刻采集针对所述动脉部位的第一散斑图,以及在第二时刻采集针对所述动脉部位的第二散斑图;acquiring a first speckle image for the arterial site at a first time, and acquiring a second speckle map for the arterial site at a second time;

确定所述第一散斑图对应的第一相位量,以及确定所述第二散斑图对应的第二相位量;determining a first phase quantity corresponding to the first speckle pattern, and determining a second phase quantity corresponding to the second speckle pattern;

根据所述第一相位量和所述第二相位量,确定所述动脉部位在所述第二时刻的脉搏空间分布信息。According to the first phase quantity and the second phase quantity, the pulse spatial distribution information of the arterial part at the second moment is determined.

可选地,根据所述散斑干涉光路,确定所述动脉部位的脉搏时间分布信息,包括:Optionally, determining the pulse time distribution information of the arterial site according to the speckle interference optical path, including:

在预设时长内采集针对所述动脉部位的多幅散斑图;collecting multiple speckle images for the arterial site within a preset time period;

确定所述多幅散斑图对应的多个相位量;determining multiple phase quantities corresponding to the multiple speckle patterns;

根据所述多个相位量,确定所述动脉部位在所述预设时长内的脉搏时间分布信息。According to the plurality of phase quantities, pulse time distribution information of the arterial part within the preset time period is determined.

可选地,所述散斑干涉光路为空间载波数字散斑干涉光路,其中,所述空间载波数字散斑干涉光路包含满足设定相干长度值的激光器、分光元件、扩束元件、光阑、成像元件和图像传感器,构建所述空间载波数字散斑干涉光路,包括:Optionally, the speckle interference light path is a space carrier digital speckle interference light path, wherein the space carrier digital speckle interference light path includes a laser, a light splitting element, a beam expanding element, a diaphragm, An imaging element and an image sensor for constructing the space carrier digital speckle interference optical path, including:

所述激光器发射激光;the laser emits laser light;

所述激光经由所述分光元件反射得到并且满足设定角度的参考光入射至所述图像传感器;The laser light is reflected by the spectroscopic element and the reference light that meets the set angle is incident on the image sensor;

所述激光经由所述分光元件透射至所述扩束元件,所述扩束元件将所述激光扩束出射后照亮所述动脉部位,经所述动脉部位反射得到的第一物光入射至所述光阑,并经由所述光阑透射至所述成像元件,其中,所述光阑用于使得所述第一物光的空间频率满足设定范围;The laser beam is transmitted to the beam expanding element through the beam splitting element, and the beam expanding element emits the laser beam expanding and then illuminates the arterial part, and the first object light reflected by the arterial part is incident on the arterial part. the diaphragm, which is transmitted to the imaging element via the diaphragm, wherein the diaphragm is used to make the spatial frequency of the first object light meet a set range;

所述成像元件成像得到的第二物光入射至所述图像传感器;The second object light obtained by imaging by the imaging element is incident on the image sensor;

所述第二物光与所述参考光在所述图像传感器表面发生干涉。The second object light and the reference light interfere on the surface of the image sensor.

可选地,所述方法还包括:Optionally, the method further includes:

根据所述脉搏空间分布信息和/或所述脉搏时间分布信息,分析所述被测生物体的健康状况。According to the pulse spatial distribution information and/or the pulse time distribution information, the health status of the measured organism is analyzed.

本申请实施例还提供一种脉搏分布的测量设备,包括:散斑干涉光路单元和确定单元,其中:The embodiment of the present application further provides a pulse distribution measurement device, including: a speckle interference optical path unit and a determination unit, wherein:

散斑干涉光路单元,用于照亮被测生物体的脉搏跳动的动脉部位;The speckle interference optical path unit is used to illuminate the arterial part of the pulse beating of the tested organism;

确定单元,用于确定所述动脉部位的脉搏分布信息,其中,所述脉搏分布信息包括:脉搏空间分布信息和脉搏时间分布信息。A determination unit, configured to determine the pulse distribution information of the arterial part, wherein the pulse distribution information includes: pulse spatial distribution information and pulse time distribution information.

本申请实施例采用的上述至少一个技术方案能够达到以下有益效果:The above-mentioned at least one technical solution adopted in the embodiments of the present application can achieve the following beneficial effects:

通过构建包含被测生物体的脉搏跳动的动脉部位的散斑干涉光路,进而根据所述散斑干涉光路,确定所述动脉部位的脉搏空间分布信息和/或脉搏时间分布信息,实现对脉搏分布的测量,为临床诊断提供有效客观的脉搏数据。By constructing the speckle interference optical path of the arterial part containing the pulse beating of the tested organism, and then determining the pulse spatial distribution information and/or pulse time distribution information of the arterial part according to the speckle interference optical path, so as to realize the pulse distribution information. It can provide effective and objective pulse data for clinical diagnosis.

附图说明Description of drawings

此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:The drawings described herein are used to provide further understanding of the present application and constitute a part of the present application. The schematic embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation of the present application. In the attached image:

图1为本申请实施例提供的一种脉搏分布的测量方法的流程示意图;1 is a schematic flowchart of a method for measuring pulse distribution according to an embodiment of the present application;

图2为本申请实施例提供的空间载波数字散斑干涉光路的光路示意图;2 is a schematic diagram of an optical path of a space carrier digital speckle interference optical path provided by an embodiment of the present application;

图3(a)为本申请实施例提供的目标对象的一幅散斑图的示意图;3(a) is a schematic diagram of a speckle image of a target object provided by an embodiment of the present application;

图3(b)为本申请实施例提供的目标对象的一幅频谱分布图的示意图;3(b) is a schematic diagram of a spectrum distribution diagram of a target object provided by an embodiment of the present application;

图3(c)为本申请实施例提供的目标对象滤波后的频谱的示意图;FIG. 3(c) is a schematic diagram of a filtered spectrum of a target object according to an embodiment of the present application;

图3(d)为本申请实施例提供的目标对象的形变相位图的示意图;FIG. 3(d) is a schematic diagram of a deformation phase diagram of a target object provided by an embodiment of the present application;

图4为本申请实施例提供的一种脉搏分布的测量设备的结构示意图。FIG. 4 is a schematic structural diagram of a pulse distribution measurement device provided by an embodiment of the present application.

具体实施方式Detailed ways

为了实现本申请的目的,本申请实施例提供了一种脉搏分布的测量方法和设备,通过构建包含被测生物体的脉搏跳动的动脉部位的散斑干涉光路,进而根据所述散斑干涉光路,确定所述动脉部位的脉搏空间分布信息和/或脉搏时间分布信息,实现对脉搏分布的测量,为临床诊断提供有效客观的脉搏数据。In order to achieve the purpose of the present application, the embodiments of the present application provide a method and device for measuring pulse distribution. , determine the pulse space distribution information and/or pulse time distribution information of the arterial part, realize the measurement of the pulse distribution, and provide effective and objective pulse data for clinical diagnosis.

下面结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。The technical solutions of the present application will be clearly and completely described below with reference to the specific embodiments of the present application and the corresponding drawings. Obviously, the described embodiments are only a part of the embodiments of the present application, but not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present application.

以下结合附图,详细说明本申请各实施例提供的技术方案。The technical solutions provided by the embodiments of the present application will be described in detail below with reference to the accompanying drawings.

实施例1Example 1

图1为本申请实施例提供的一种脉搏分布的测量方法的流程示意图。所述方法可以如下所示。FIG. 1 is a schematic flowchart of a method for measuring pulse distribution according to an embodiment of the present application. The method can be as follows.

步骤11:构建散斑干涉光路。Step 11: Construct the speckle interference light path.

其中,散斑干涉光路中包含被测生物体的脉搏跳动的动脉部位。Among them, the speckle interference optical path includes the arterial part where the pulse of the living body to be measured is beating.

在步骤11中,构建包含被测生物体的脉搏跳动的动脉部位的散斑干涉光路,为进一步根据该散斑干涉光路进行脉搏分布的测量做准备。In step 11 , a speckle interference optical path including the arterial part of the pulse beating of the living body is constructed, and preparations are made for further measuring the pulse distribution according to the speckle interference optical path.

其中,被测生物体动脉部位表示该被测生物体的脉搏跳动明显、易测量部位,例如:人体近手腕部的桡动脉,这里不做具体限定。Wherein, the arterial part of the measured body refers to the part where the pulse of the measured body is obvious and easy to measure, such as the radial artery of the human body near the wrist, which is not specifically limited here.

可选地,散斑干涉光路为空间载波数字散斑干涉光路。Optionally, the speckle interference light path is a space carrier digital speckle interference light path.

脉搏跳动是一种振幅微小的不稳定振动,属于动态形变,空间载波数字散斑干涉光路能够对脉搏的动态形变进行测量,进而得到脉搏对应的空间分布和时间分布。Pulse beat is a kind of unstable vibration with small amplitude, which belongs to dynamic deformation. The space carrier digital speckle interference optical path can measure the dynamic deformation of pulse, and then obtain the corresponding spatial distribution and time distribution of pulse.

空间载波数字散斑干涉光路通过将参考光与光轴倾斜设定角度产生空间载波信号。在对测量得到的散斑图实施傅里叶变换之后,得到包含目标对象(被测生物体的脉搏跳动的动脉部位)形变信息的高频部分和代表背景光的低频部分,由于空间载波信号的引入使频谱中的高频部分发生频移,实现与低频部分的分离,所以通过提取包含目标对象形变信息的高频部分进而求解干涉相位,最后根据干涉相位与目标对象形变量之间的关系,确定目标对象的形变量。The space carrier digital speckle interference optical path generates a space carrier signal by tilting the reference light with the optical axis at a set angle. After Fourier transform is performed on the measured speckle pattern, the high-frequency part and the low-frequency part representing the background light containing the deformation information of the target object (the pulsed arterial part of the measured organism) are obtained. The introduction makes the high frequency part in the spectrum shift the frequency to realize the separation from the low frequency part, so by extracting the high frequency part containing the deformation information of the target object and then solving the interference phase, and finally according to the relationship between the interference phase and the deformation variable of the target object, Determines the amount of deformation of the target object.

需要说明的是,空间载波数字散斑干涉光路可以是在空间环境中搭建的空间光路,也可以是集成在测量仪器中的光路,这里不做具体限定。It should be noted that the space carrier digital speckle interference optical path may be a space optical path built in a space environment, or may be an optical path integrated in a measuring instrument, which is not specifically limited here.

图2为本申请实施例提供的空间载波数字散斑干涉光路的光路示意图。FIG. 2 is a schematic diagram of an optical path of a space carrier digital speckle interference optical path provided by an embodiment of the present application.

如图2所示,可选地,空间载波数字散斑干涉光路包含满足设定相干长度值的激光器、分光元件、扩束元件、光阑、成像元件和图像传感器,构建该空间载波数字散斑干涉光路,包括:As shown in Figure 2, optionally, the space carrier digital speckle interference optical path includes lasers, beam splitting elements, beam expanding elements, diaphragms, imaging elements and image sensors that satisfy the set coherence length value to construct the space carrier digital speckle Interfering optical paths, including:

激光器发射激光;Laser emits laser light;

激光经由分光元件反射得到并且满足设定角度的参考光入射至图像传感器;The laser is reflected by the spectroscopic element and the reference light that meets the set angle is incident on the image sensor;

激光经由分光元件透射至扩束元件,扩束元件将该激光扩束出射后照亮动脉部位,经动脉部位反射得到的第一物光入射至光阑,并经由光阑透射至成像元件,其中,该光阑用于使得该第一物光的空间频率满足设定范围;The laser is transmitted to the beam expanding element through the beam splitting element, and the beam expanding element emits the laser beam expanding and then illuminates the arterial part, the first object light obtained by the reflection of the arterial part is incident on the diaphragm, and is transmitted to the imaging element through the diaphragm, wherein , the diaphragm is used to make the spatial frequency of the first object light meet the set range;

成像元件成像得到的第二物光入射至图像传感器;The second object light obtained by imaging by the imaging element is incident on the image sensor;

第二物光与参考光在图像传感器表面发生干涉。The second object light interferes with the reference light on the surface of the image sensor.

其中,设定角度(α)满足公式:Among them, the set angle (α) satisfies the formula:

Figure BDA0001260431890000051
Figure BDA0001260431890000051

其中,λ表示激光器发射激光的波长,fM表示第一物光的空间频率,fP表示图像传感器的时间采样频率。Among them, λ represents the wavelength of laser light emitted by the laser, f M represents the spatial frequency of the first object light, and f P represents the time sampling frequency of the image sensor.

需要说明的是,设定相干长度值表示激光器出射的激光能够保证参考光与第二物光发生干涉,该设定相干长度值的大小可以根据实际情况确定,这里不做具体限定。It should be noted that setting the coherence length value means that the laser emitted by the laser can ensure that the reference light and the second object light interfere with each other.

本申请实施例中,图像传感器的时间采样频率在满足采样频率的基础上,具体取值可以根据实际情况确定,这里不做具体限定。In the embodiment of the present application, the time sampling frequency of the image sensor is based on satisfying the sampling frequency, and the specific value may be determined according to the actual situation, which is not specifically limited here.

需要说明的是,为了能够根据空间载波数字散斑干涉光路,确定动脉部位的脉搏分布,第一物光的空间频率fM满足设定范围:It should be noted that, in order to determine the pulse distribution of the arterial site according to the spatial carrier digital speckle interference optical path, the spatial frequency f M of the first object light satisfies the set range:

Figure BDA0001260431890000061
Figure BDA0001260431890000061

由于光阑用于使得第一物光的空间频率满足设定范围,因此,该光阑的孔径尺寸满足:Since the diaphragm is used to make the spatial frequency of the first object light meet the set range, the aperture size of the diaphragm satisfies:

Figure BDA0001260431890000062
Figure BDA0001260431890000062

其中,d表示该光阑到成像元件的距离。Among them, d represents the distance from the diaphragm to the imaging element.

可选地,分光元件为可调分光平片。Optionally, the light-splitting element is a tunable light-splitting plate.

分光元件除了可以采用可调分光平片之外,也可以采用可调分光镜,还可以采用其他能够进行分光的光学元件,这里不做具体限定。In addition to the adjustable beam splitter plate, the beam splitting element can also be an adjustable beam splitter, and other optical elements capable of beam splitting can also be used, which is not specifically limited here.

可选地,扩束元件为显微物镜。Optionally, the beam expander element is a microscope objective.

扩束元件除了可以采用所述显微物镜之外,也可以采用平凹透镜,还可以采用其他能够进行扩束的光学元件,这里不做具体限定。In addition to the microscope objective lens, the beam expander can also use a plano-concave lens, and other optical elements capable of beam expansion can also be used, which is not specifically limited here.

优选地,成像元件为平凸透镜。Preferably, the imaging element is a plano-convex lens.

成像元件除了可以采用平凸透镜之外,也可以采用工业镜头,还可以采用其他能够进行成像的光学元件,这里不做具体限定。In addition to a plano-convex lens, an industrial lens can also be used as the imaging element, and other optical elements capable of imaging can also be used, which is not specifically limited here.

本申请实施例中,所述散斑干涉光路除了可以采用空间载波数字散斑干涉光路之外,还可以采用4+1时间相移数字散斑干涉光路,还可以采用其他能够测量动态变形的散斑干涉光路,这里不做具体限定。In the embodiment of the present application, in addition to the space carrier digital speckle interference light path, the speckle interference light path may also use a 4+1 time phase shift digital speckle interference light path, and other speckle interference light paths capable of measuring dynamic deformation may also be used. The speckle interference optical path is not specifically limited here.

步骤12:根据散斑干涉光路,确定所述动脉部位的脉搏分布信息。Step 12: Determine the pulse distribution information of the arterial site according to the speckle interference optical path.

其中,脉搏分布信息包括:脉搏空间分布信息和/或脉搏时间分布信息。The pulse distribution information includes: pulse spatial distribution information and/or pulse time distribution information.

在步骤12中,根据步骤11中构建的散斑干涉光路,确定被测生物体的脉搏跳动的动脉部位的脉搏空间分布信息和/或脉搏时间分布信息。In step 12, according to the speckle interference optical path constructed in step 11, the pulse spatial distribution information and/or pulse time distribution information of the arterial part of the pulse beating of the measured organism is determined.

可选地,根据散斑干涉光路,确定动脉部位的脉搏空间分布信息,具体地,第一步,在第一时刻采集针对动脉部位的第一散斑图,以及在第二时刻采集针对动脉部位的第二散斑图。Optionally, according to the speckle interference optical path, the pulse spatial distribution information of the arterial site is determined. Specifically, in the first step, the first speckle map for the arterial site is collected at the first moment, and the arterial site is collected at the second moment. of the second speckle pattern.

例如,以空间载波数字散斑干涉光路为例,由于成像元件成像得到的第二物光是对动脉部位进行成像之后得到的,该第二物光与参考光在图像传感器表面发生干涉,因此,能够通过该图像传感器采集到针对该动脉部位的干涉散斑图(以下简称散斑图)。For example, taking the space carrier digital speckle interference optical path as an example, since the second object light obtained by imaging by the imaging element is obtained after imaging the artery, the second object light and the reference light interfere on the surface of the image sensor, therefore, An interference speckle pattern (hereinafter referred to as a speckle pattern) for the arterial part can be collected by the image sensor.

图像传感器分别在第一时刻和第二时刻采集针对动脉部位的第一散斑图和第二散斑图。The image sensor acquires a first speckle pattern and a second speckle pattern for the arterial site at a first time instant and a second time instant, respectively.

第二步,确定第一散斑图对应的第一相位量,以及确定第二散斑图对应的第二相位量。In the second step, the first phase quantity corresponding to the first speckle pattern is determined, and the second phase quantity corresponding to the second speckle pattern is determined.

第三步,根据第一相位量和第二相位量,确定动脉部位在第二时刻的脉搏空间分布信息。In the third step, according to the first phase quantity and the second phase quantity, determine the pulse spatial distribution information of the arterial part at the second moment.

图3(a)为本申请实施例提供的目标对象的第一散斑图的示意图。FIG. 3( a ) is a schematic diagram of a first speckle pattern of a target object according to an embodiment of the present application.

仍以上述空间载波数字散斑干涉光路为例,对目标对象(例如:动脉部位)对应的第一散斑图(例如,图3(a))进行傅里叶变换,得到该第一散斑图对应的频谱分布图,图3(b)为本申请实施例提供的目标对象的频谱分布图的示意图。Still taking the above space carrier digital speckle interference optical path as an example, Fourier transform is performed on the first speckle pattern (for example, Fig. 3(a)) corresponding to the target object (for example, the arterial part) to obtain the first speckle Fig. 3(b) is a schematic diagram of a spectrum distribution diagram of a target object provided by an embodiment of the present application.

第一散斑图对应的频谱分布为:The spectral distribution corresponding to the first speckle pattern is:

I(fx)=FT[I(x)]=A(fx)+B(fx)+C(fx),I(f x )=FT[I(x)]=A(f x )+B(f x )+C(f x ),

其中,I(fx)表示第一散斑图对应的频谱分布,I(x)表示目标对象的光强度,A(fx)和B(fx)表示两个共轭的高频部分,C(fx)表示低频部分,fx表示频域坐标,x表示一维空间坐标,符号FT表示傅立叶变换。Among them, I(f x ) represents the spectral distribution corresponding to the first speckle pattern, I(x) represents the light intensity of the target object, A(f x ) and B(f x ) represent the two conjugated high-frequency parts, C(f x ) represents the low-frequency part, f x represents the frequency domain coordinate, x represents the one-dimensional space coordinate, and the symbol FT represents the Fourier transform.

在频谱分布图中:In the spectrum distribution plot:

A(fx)=FT[O(x)R*(x)];A(f x )=FT[O(x)R * (x)];

B(fx)=FT[R(x)O*(x)];B(f x )=FT[R(x)O * (x)];

C(fx)=FT[|R(x)|2+|O(x)|2]。C(f x )=FT[|R(x)| 2 +|O(x)| 2 ].

其中,R(x)表示目标对象的物光的光强度,O(x)表示参考光的光强度,符号*表示共轭。Among them, R(x) represents the light intensity of the object light of the target object, O(x) represents the light intensity of the reference light, and the symbol * represents the conjugate.

经过频谱窗口滤波,得到其中一个高频部分A(fx),滤去其他两个频谱部分B(fx)和C(fx),图3(c)为本申请实施例提供的目标对象滤波后的频谱的示意图。After spectral window filtering, one of the high-frequency parts A(f x ) is obtained, and the other two spectral parts B(f x ) and C(f x ) are filtered out. FIG. 3( c ) is the target object provided by the embodiment of the present application Schematic of the filtered spectrum.

对滤波后频谱进行傅里叶逆变换,根据下述公式,得到含有目标对象的第一相位量的表达式:Inverse Fourier transform is performed on the filtered spectrum, and the expression containing the first phase quantity of the target object is obtained according to the following formula:

Figure BDA0001260431890000081
Figure BDA0001260431890000081

其中,AIFT表示含有目标对象的第一相位量的表达式,

Figure BDA0001260431890000082
表示目标对象的第一相位量,A(x)表示目标对象在得到第一散斑图时的干涉光的光强信息,符号IFT表示傅立叶逆变换。where A IFT represents an expression containing the first phase quantity of the target object,
Figure BDA0001260431890000082
represents the first phase quantity of the target object, A(x) represents the light intensity information of the interference light of the target object when the first speckle pattern is obtained, and the symbol IFT represents the inverse Fourier transform.

同理,对于目标对象对应的第二散斑图,同样可以通过执行上述傅里叶变换、滤波处理以及傅里叶逆变换的操作,得到含有目标对象的第二相位量的表达式:Similarly, for the second speckle pattern corresponding to the target object, the expression containing the second phase quantity of the target object can also be obtained by performing the operations of the above Fourier transform, filtering processing and inverse Fourier transform:

Figure BDA0001260431890000083
Figure BDA0001260431890000083

其中,AIFT'表示含有目标对象的第二相位量的表达式,

Figure BDA0001260431890000084
表示目标对象的第二相位量,A'(x)表示目标对象在得到第二散斑图时的干涉光的光强信息,符号IFT表示傅立叶逆变换。where A IFT ' represents an expression containing the second phase quantity of the target object,
Figure BDA0001260431890000084
represents the second phase quantity of the target object, A'(x) represents the light intensity information of the interference light of the target object when the second speckle pattern is obtained, and the symbol IFT represents the inverse Fourier transform.

在得到含有目标对象的第一相位量的表达式和含有目标对象的第二相位量的表达式后,根据含有目标对象的第一相位量的表达式和含有目标对象的第二相位量的表达式,通过下述公式,得到目标对象的形变信息。After obtaining the expression containing the first phase quantity of the target object and the expression containing the second phase quantity of the target object, according to the expression containing the first phase quantity of the target object and the expression containing the second phase quantity of the target object The following formula is used to obtain the deformation information of the target object.

Figure BDA0001260431890000091
Figure BDA0001260431890000091

其中,Im和Re分别表示虚部和实部,符号*表示共轭。Among them, Im and Re represent the imaginary part and the real part, respectively, and the symbol * represents the conjugation.

图3(d)为本申请实施例提供的目标对象的形变相位图的示意图。FIG. 3( d ) is a schematic diagram of a deformation phase diagram of a target object according to an embodiment of the present application.

根据目标对象的形变信息,确定目标对象在第二时刻的空间分布信息,包括:According to the deformation information of the target object, determine the spatial distribution information of the target object at the second moment, including:

根据目标对象的形变信息,通过以下公式确定目标对象在第二时刻的空间分布信息:According to the deformation information of the target object, the spatial distribution information of the target object at the second moment is determined by the following formula:

Figure BDA0001260431890000092
Figure BDA0001260431890000092

其中,

Figure BDA0001260431890000093
表示目标对象的空间分布信息,即目标对象的各个位置在第二时刻的形变分布,
Figure BDA0001260431890000094
表示目标对象的形变信息,
Figure BDA0001260431890000095
表示位移敏感矢量,
Figure BDA0001260431890000096
ki表示照亮动脉部位方向的单位矢量,ko表示第一物光方向的单位矢量,λ表示激光器发射激光的波长。in,
Figure BDA0001260431890000093
Represents the spatial distribution information of the target object, that is, the deformation distribution of each position of the target object at the second moment,
Figure BDA0001260431890000094
Represents the deformation information of the target object,
Figure BDA0001260431890000095
represents the displacement-sensitive vector,
Figure BDA0001260431890000096
k i represents the unit vector in the direction of the illuminated arterial site, k o represents the unit vector in the direction of the first object light, and λ represents the wavelength of the laser emitted by the laser.

需要说明的是,目标对象表示被测生物体的脉搏跳动的动脉部位,即采用上述记载的方案,确定了被测生物体的动脉部位的脉搏空间分布信息。It should be noted that the target object represents the arterial part of the measured body's pulse, that is, the above-described scheme is used to determine the pulse spatial distribution information of the arterial part of the measured body.

可选地,根据散斑干涉光路,确定动脉部位的脉搏时间分布信息,具体地,第一步,在预设时长内采集针对动脉部位的多幅散斑图。Optionally, the pulse time distribution information of the arterial part is determined according to the speckle interference optical path. Specifically, in the first step, a plurality of speckle images for the arterial part are collected within a preset time period.

图像传感器在预设时长内采集针对动脉部位的多幅散斑图。The image sensor acquires multiple speckle images for the arterial site within a preset time period.

需要说明的是,多幅散斑图的数量与预设时长、图像传感器的时间采样频率有关。It should be noted that the number of the multiple speckle patterns is related to the preset duration and the time sampling frequency of the image sensor.

例如:图像传感器的时间采样频率为75帧/秒,预设时长为10秒,在预设时长10秒内采集得到的多幅散斑图的数量为750。For example, the time sampling frequency of the image sensor is 75 frames/second, the preset duration is 10 seconds, and the number of multiple speckle images acquired within the preset duration of 10 seconds is 750.

第二步,确定多幅散斑图对应的多个相位量。In the second step, multiple phase quantities corresponding to multiple speckle patterns are determined.

第三,根据多个相位量,确定动脉部位在该预设时长内的脉搏时间分布信息。Third, according to a plurality of phase quantities, the pulse time distribution information of the arterial part within the preset time period is determined.

仍以上述空间载波数字散斑干涉光路为例,通过空间载波法,如上述图3(a)、图3(b)和图3(c)记载的根据目标对象的散斑图,确定散斑图对应的相位量的操作,确定针对动脉部位采集得到的多幅散斑图对应的多个相位量;以及,通过上述记载的根据第一相位量和第二相位量,确定目标对象的形变信息的操作,根据多个相位量中的任意相邻两个相位量,确定动脉部位的多个脉搏形变信息。Still taking the above space carrier digital speckle interference optical path as an example, through the space carrier method, the speckle is determined according to the speckle diagram of the target object as described in the above Figures 3(a), 3(b) and 3(c). By operating the phase quantities corresponding to the image, determine multiple phase quantities corresponding to the multiple speckle patterns collected for the arterial site; and determine the deformation information of the target object according to the first phase quantity and the second phase quantity as described above operation, according to any two adjacent phase quantities in the plurality of phase quantities, to determine a plurality of pulse deformation information of the arterial part.

例如,在预设时长内针对动脉部位A采集得到4辐散斑图,确定4幅散斑图对应的4个相位量a、b、c和d。根据相位量a和相位量b,确定动脉部位的第一脉搏形变信息

Figure BDA0001260431890000101
根据相位量b和相位量c,确定动脉部位的第二脉搏形变信息
Figure BDA0001260431890000102
根据相位量c和相位量d,确定动脉部位的第三脉搏形变信息
Figure BDA0001260431890000103
For example, 4-radial speckle patterns are acquired for the arterial site A within a preset time period, and 4 phase quantities a, b, c, and d corresponding to the 4 speckle patterns are determined. According to the phase quantity a and the phase quantity b, the first pulse deformation information of the arterial part is determined
Figure BDA0001260431890000101
According to the phase quantity b and the phase quantity c, the second pulse deformation information of the arterial part is determined
Figure BDA0001260431890000102
According to the phase quantity c and the phase quantity d, the third pulse deformation information of the arterial part is determined
Figure BDA0001260431890000103

需要说明的是,确定多幅散斑图对应的多个相位量的方法,可以采用空间载波法,可以采用4+1时间相移法,还可以采用其他能够测量动态变形的方法,这里不做具体限定。It should be noted that the method of determining the multiple phase quantities corresponding to the multiple speckle patterns can be the space carrier method, the 4+1 time phase shift method, or other methods that can measure dynamic deformation, which are not described here. Specific restrictions.

例如,采用4+1时间相移法确定多个相位量时,首先,通过四步相移法采集针对动脉部位的4幅初始散斑图;其次,在预设时长内采集针对动脉部位的多幅散斑图;最后,根据该多幅散斑图中的任一幅与该4幅初始散斑图,通过时间相移算法确定多幅散斑图对应的多个相位量。For example, when the 4+1 time phase shift method is used to determine multiple phase quantities, first, four initial speckle images for the arterial site are collected by the four-step phase shift method; and finally, according to any one of the multiple speckle images and the four initial speckle images, a time phase shift algorithm is used to determine multiple phase quantities corresponding to the multiple speckle images.

通过上述图3(d)记载的根据目标对象的形变信息,确定目标对象的形变分布的方法,确定动脉部位的多个脉搏形变分布。Through the method of determining the deformation distribution of the target object based on the deformation information of the target object described in FIG. 3( d ), a plurality of pulse deformation distributions of the arterial site are determined.

仍以上述动脉部位A为例,根据第一脉搏形变信息,确定第一脉搏形变分布;根据第二脉搏形变信息,确定第二脉搏形变分布;根据第三脉搏形变信息,确定第三脉搏形变分布。Still taking the above arterial part A as an example, the first pulse deformation distribution is determined according to the first pulse deformation information; the second pulse deformation distribution is determined according to the second pulse deformation information; and the third pulse deformation distribution is determined according to the third pulse deformation information .

第一脉搏形变分布、第二脉搏形变分布和第三脉搏形变分布共同构成了动脉部位在预设时长内的脉搏时间分布信息。The first pulse deformation distribution, the second pulse deformation distribution and the third pulse deformation distribution together constitute the pulse time distribution information of the arterial part within the preset time period.

在本申请的可选实施例中,所述方法还包括:In an optional embodiment of the present application, the method further includes:

根据所述脉搏空间分布信息和/或所述脉搏时间分布信息,分析所述被测生物体的健康状况。According to the pulse spatial distribution information and/or the pulse time distribution information, the health status of the measured organism is analyzed.

根据被测生物体的动脉部位的脉搏空间分布信息,可以得到脉搏的脉位;根据被测生物体的动脉部位的脉搏时间分布信息,可以得到脉搏的脉数;综合被测生物体的动脉部位的脉搏空间分布信息和脉搏时间分布信息,可以得到脉搏的脉形和脉势。根据脉搏的脉位、脉数、脉形、脉势等脉象信息,能够分析被测生物体的健康状况。The pulse position of the pulse can be obtained according to the pulse space distribution information of the arterial part of the tested organism; the pulse number of the pulse can be obtained according to the pulse time distribution information of the arterial part of the tested organism; the arterial part of the tested organism can be integrated The pulse spatial distribution information and pulse time distribution information can be obtained, and the pulse shape and pulse potential of the pulse can be obtained. According to the pulse information such as pulse position, pulse number, pulse shape, pulse potential, etc., the health status of the tested organism can be analyzed.

本申请实施例所记载的技术方案,通过构建包含被测生物体的脉搏跳动的动脉部位的散斑干涉光路,进而根据所述散斑干涉光路,确定所述动脉部位的脉搏空间分布信息和/或脉搏时间分布信息,实现对脉搏分布的测量,为临床诊断提供有效客观的脉搏数据。According to the technical solutions described in the embodiments of the present application, by constructing a speckle interference optical path of the arterial part containing the pulse beating of the organism under test, and then according to the speckle interference optical path, the pulse spatial distribution information and/or the arterial part can be determined. Or pulse time distribution information, realize the measurement of pulse distribution, and provide effective and objective pulse data for clinical diagnosis.

实施例2Example 2

图4为本申请实施例提供的一种脉搏形变的测量设备的结构示意图,设备400包括:散斑干涉光路单元401和确定单元402,其中:FIG. 4 is a schematic structural diagram of a pulse deformation measurement device provided by an embodiment of the application. The device 400 includes: a speckle interference optical path unit 401 and a determination unit 402, wherein:

散斑干涉光路单元401,用于照亮被测生物体的脉搏跳动的动脉部位;The speckle interference optical path unit 401 is used to illuminate the arterial part of the pulse beating of the tested organism;

确定单元402,用于确定动脉部位的脉搏分布信息,其中,脉搏分布信息包括:脉搏空间分布信息和/或脉搏时间分布信息。The determining unit 402 is configured to determine pulse distribution information of the arterial site, wherein the pulse distribution information includes: pulse spatial distribution information and/or pulse time distribution information.

可选地,确定单元402确定动脉部位的脉搏空间分布信息,包括:Optionally, the determining unit 402 determines the pulse spatial distribution information of the arterial site, including:

在第一时刻采集针对动脉部位的第一散斑图,以及在第二时刻采集针对动脉部位的第二散斑图;acquiring a first speckle image for the arterial site at a first time instant, and acquiring a second speckle map for the arterial site at a second time instant;

确定第一散斑图对应的第一相位量,以及确定第二散斑图对应的第二相位量;determining a first phase quantity corresponding to the first speckle pattern, and determining a second phase quantity corresponding to the second speckle pattern;

根据第一相位量和第二相位量,确定动脉部位在第二时刻的脉搏空间分布信息。According to the first phase quantity and the second phase quantity, the pulse spatial distribution information of the arterial part at the second moment is determined.

可选地,确定单元402确定动脉部位的脉搏时间分布信息,包括:Optionally, the determining unit 402 determines the pulse time distribution information of the arterial site, including:

在预设时长内采集针对动脉部位的多幅散斑图;Acquire multiple speckle images for the arterial site within a preset time period;

确定多幅散斑图对应的多个相位量;Determine multiple phase quantities corresponding to multiple speckle patterns;

根据多个相位量,确定动脉部位在预设时长内的脉搏时间分布信息。According to the multiple phase quantities, the pulse time distribution information of the arterial part within the preset time period is determined.

可选地,散斑干涉光路单元401为空间载波数字散斑干涉光路单元,其中,该空间载波数字散斑干涉光路单元中包含满足设定相干长度值的激光器,分光元件、扩束元件、光阑、成像元件和图像传感器,其中:Optionally, the speckle interference optical path unit 401 is a spatial carrier digital speckle interference optical path unit, wherein the spatial carrier digital speckle interference optical path unit includes a laser that satisfies the set coherence length value, a light splitting element, a beam expanding element, a light diaphragms, imaging elements, and image sensors, including:

激光器发射激光;Laser emits laser light;

激光经由分光元件反射得到并且满足设定角度的参考光入射至图像传感器;The laser is reflected by the spectroscopic element and the reference light that meets the set angle is incident on the image sensor;

激光经由分光元件透射至扩束元件,扩束元件将激光扩束出射后照亮动脉部位,经动脉部位反射得到的第一物光入射至光阑,并经由光阑透射至成像元件,其中,光阑用于使得第一物光的空间频率满足设定范围;The laser is transmitted to the beam expanding element through the beam splitting element, and the beam expanding element illuminates the arterial part after expanding the laser beam, and the first object light reflected by the arterial part is incident on the diaphragm, and is transmitted to the imaging element through the diaphragm, wherein, The diaphragm is used to make the spatial frequency of the first object light meet the set range;

成像元件成像得到的第二物光入射至图像传感器;The second object light obtained by imaging by the imaging element is incident on the image sensor;

第二物光与参考光在图像传感器表面发生干涉。The second object light interferes with the reference light on the surface of the image sensor.

可选地,设备400还包括:分析单元,其中:Optionally, the device 400 further includes: an analysis unit, wherein:

分析单元,用于根据脉搏空间分布信息和/或脉搏时间分布信息,分析被测生物体的健康状况。The analyzing unit is configured to analyze the health status of the measured organism according to the pulse spatial distribution information and/or the pulse time distribution information.

根据所述脉搏分布的测量设备,所述散斑干涉光路单元用于照亮被测生物体的脉搏跳动的动脉部位;所述确定单元用于确定所述动脉部位的脉搏空间分布信息和/或脉搏时间分布信息,实现对脉搏分布的测量,为临床诊断提供有效客观的脉搏数据。According to the pulse distribution measuring device, the speckle interference optical path unit is used for illuminating the arterial part of the pulse beating of the measured body; the determining unit is used for determining the pulse space distribution information and/or Pulse time distribution information, realize the measurement of pulse distribution, and provide effective and objective pulse data for clinical diagnosis.

本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block in the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to the processor of a general purpose computer, special purpose computer, embedded processor or other programmable data processing device to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing device produce Means for implementing the functions specified in a flow or flow of a flowchart and/or a block or blocks of a block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer readable memory capable of directing a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory result in an article of manufacture comprising instruction means, the instructions The apparatus implements the functions specified in the flow or flow of the flowcharts and/or the block or blocks of the block diagrams.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded on a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that The instructions provide steps for implementing the functions specified in the flow or blocks of the flowcharts and/or the block or blocks of the block diagrams.

在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.

内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。Memory may include non-persistent memory in computer readable media, random access memory (RAM) and/or non-volatile memory in the form of, for example, read only memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。Computer-readable media includes both persistent and non-permanent, removable and non-removable media, and storage of information may be implemented by any method or technology. Information may be computer readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read only memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), Flash Memory or other memory technology, Compact Disc Read Only Memory (CD-ROM), Digital Versatile Disc (DVD) or other optical storage, Magnetic tape cartridges, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory computer-readable media, such as modulated data signals and carrier waves.

还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。It should also be noted that the terms "comprising", "comprising" or any other variation thereof are intended to encompass a non-exclusive inclusion such that a process, method, article or device comprising a series of elements includes not only those elements, but also Other elements not expressly listed, or which are inherent to such a process, method, article of manufacture, or apparatus are also included. Without further limitation, an element qualified by the phrase "comprising a..." does not preclude the presence of additional identical elements in the process, method, article of manufacture or device that includes the element.

本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。It will be appreciated by those skilled in the art that the embodiments of the present application may be provided as a method, a system or a computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。The above descriptions are merely examples of the present application, and are not intended to limit the present application. Various modifications and variations of this application are possible for those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of this application shall be included within the scope of the claims of this application.

Claims (5)

1.一种脉搏分布的测量方法,其特征在于,包括:1. a measuring method of pulse distribution, is characterized in that, comprises: 构建散斑干涉光路,其中,所述散斑干涉光路中包含被测生物体的脉搏跳动的动脉部位、激光器、分光元件、扩束元件、成像元件和图像传感器;构建所述散斑干涉光路包括:Constructing a speckle interference optical path, wherein the speckle interference optical path includes the pulsed arterial part, laser, beam splitting element, beam expanding element, imaging element and image sensor of the tested organism; constructing the speckle interference optical path includes: : 所述激光器发射激光;the laser emits laser light; 所述激光经由所述分光元件反射得到参考光,所述参考光入射至所述图像传感器;The laser light is reflected by the spectroscopic element to obtain reference light, and the reference light is incident on the image sensor; 所述激光经由所述分光元件透射至所述扩束元件,所述扩束元件将所述激光扩束出射后照亮所述动脉部位,经所述动脉部位反射得到第一物光,所述第一物光入射至所述成像元件,所述成像元件成像得到的第二物光入射至所述图像传感器;The laser beam is transmitted to the beam expanding element through the beam splitting element, the beam expanding element emits the laser beam expanding and then illuminates the arterial part, and is reflected by the arterial part to obtain the first object light, and the The first object light is incident on the imaging element, and the second object light obtained by imaging by the imaging element is incident on the image sensor; 所述第二物光与所述参考光在所述图像传感器表面发生干涉;the second object light and the reference light interfere on the surface of the image sensor; 根据所述散斑干涉光路,确定所述动脉部位的脉搏分布信息,其中,所述脉搏分布信息包括:脉搏空间分布信息和脉搏时间分布信息;Determine the pulse distribution information of the arterial part according to the speckle interference optical path, wherein the pulse distribution information includes: pulse spatial distribution information and pulse time distribution information; 其中,根据所述散斑干涉光路,确定所述动脉部位的脉搏空间分布信息,包括:Wherein, determining the pulse spatial distribution information of the arterial site according to the speckle interference optical path, including: 在第一时刻采集针对所述动脉部位的第一散斑图,以及在第二时刻采集针对所述动脉部位的第二散斑图;acquiring a first speckle image for the arterial site at a first time, and acquiring a second speckle map for the arterial site at a second time; 确定所述第一散斑图对应的第一相位量,以及确定所述第二散斑图对应的第二相位量;determining a first phase quantity corresponding to the first speckle pattern, and determining a second phase quantity corresponding to the second speckle pattern; 根据所述第一相位量和所述第二相位量,确定所述动脉部位在所述第二时刻的脉搏空间分布信息;According to the first phase amount and the second phase amount, determine the pulse space distribution information of the arterial part at the second moment; 其中,根据所述散斑干涉光路,确定所述动脉部位的脉搏时间分布信息,包括:Wherein, according to the speckle interference optical path, the pulse time distribution information of the arterial part is determined, including: 在预设时长内采集针对所述动脉部位的多幅散斑图;collecting multiple speckle images for the arterial site within a preset time period; 确定所述多幅散斑图对应的多个相位量;determining multiple phase quantities corresponding to the multiple speckle patterns; 根据所述多个相位量,确定所述动脉部位在所述预设时长内的脉搏时间分布信息。According to the plurality of phase quantities, pulse time distribution information of the arterial part within the preset time period is determined. 2.如权利要求1所述的方法,其特征在于,所述散斑干涉光路为空间载波数字散斑干涉光路,其中,所述激光器为满足设定相干长度值的激光器,入射至所述图像传感器的参考光为满足设定角度的参考光;所述空间载波数字散斑干涉光路还包含光阑;2 . The method according to claim 1 , wherein the speckle interference optical path is a spatial carrier digital speckle interference optical path, wherein the laser is a laser satisfying a set coherence length value and is incident on the image. 3 . The reference light of the sensor is the reference light that satisfies the set angle; the space carrier digital speckle interference light path also includes a diaphragm; 所述第一物光入射至所述成像元件包括:The first object light incident on the imaging element includes: 经所述动脉部位反射得到的第一物光入射至所述光阑,并经由所述光阑透射至所述成像元件,其中,所述光阑用于使得所述第一物光的空间频率满足设定范围。The first object light reflected by the arterial part is incident on the diaphragm, and is transmitted to the imaging element through the diaphragm, wherein the diaphragm is used to make the spatial frequency of the first object light meet the setting range. 3.一种脉搏分布的测量设备,其特征在于,所述设备包括:散斑干涉光路单元和确定单元,其中:3. A measurement device for pulse distribution, characterized in that the device comprises: a speckle interference optical path unit and a determination unit, wherein: 散斑干涉光路单元,用于照亮被测生物体的脉搏跳动的动脉部位;The speckle interference optical path unit is used to illuminate the arterial part of the pulse beating of the tested organism; 其中,所述散斑干涉光路单元包含激光器、分光元件、扩束元件、成像元件和图像传感器;构建所述散斑干涉光路包括:Wherein, the speckle interference optical path unit includes a laser, a light splitting element, a beam expanding element, an imaging element and an image sensor; constructing the speckle interference optical path includes: 所述激光经由所述分光元件反射得到参考光,所述参考光入射至所述图像传感器;The laser light is reflected by the spectroscopic element to obtain reference light, and the reference light is incident on the image sensor; 所述激光经由所述分光元件透射至所述扩束元件,所述扩束元件将所述激光扩束出射后照亮所述动脉部位,经所述动脉部位反射得到第一物光,所述第一物光入射至所述成像元件,所述成像元件成像得到的第二物光入射至所述图像传感器;The laser beam is transmitted to the beam expanding element through the beam splitting element, the beam expanding element emits the laser beam expanding and then illuminates the arterial part, and is reflected by the arterial part to obtain the first object light, and the The first object light is incident on the imaging element, and the second object light obtained by imaging by the imaging element is incident on the image sensor; 所述第二物光与所述参考光在所述图像传感器表面发生干涉;the second object light and the reference light interfere on the surface of the image sensor; 确定单元,用于确定所述动脉部位的脉搏分布信息,其中,所述脉搏分布信息包括:脉搏空间分布信息和脉搏时间分布信息;a determining unit, configured to determine the pulse distribution information of the arterial part, wherein the pulse distribution information includes: pulse spatial distribution information and pulse time distribution information; 其中,所述确定单元确定所述动脉部位的脉搏空间分布信息,包括:Wherein, the determining unit determines the pulse spatial distribution information of the arterial part, including: 在第一时刻采集针对所述动脉部位的第一散斑图,以及在第二时刻采集针对所述动脉部位的第二散斑图;acquiring a first speckle image for the arterial site at a first time, and acquiring a second speckle map for the arterial site at a second time; 确定所述第一散斑图对应的第一相位量,以及确定所述第二散斑图对应的第二相位量;determining a first phase quantity corresponding to the first speckle pattern, and determining a second phase quantity corresponding to the second speckle pattern; 根据所述第一相位量和所述第二相位量,确定所述动脉部位在所述第二时刻的脉搏空间分布信息;According to the first phase amount and the second phase amount, determine the pulse space distribution information of the arterial part at the second moment; 其中,所述确定单元确定所述动脉部位的脉搏时间分布信息,包括:Wherein, the determining unit determines the pulse time distribution information of the arterial part, including: 在预设时长内采集针对所述动脉部位的多幅散斑图;collecting multiple speckle images for the arterial site within a preset time period; 确定所述多幅散斑图对应的多个相位量;determining multiple phase quantities corresponding to the multiple speckle patterns; 根据所述多个相位量,确定所述动脉部位在所述预设时长内的时间分布信息。According to the plurality of phase quantities, time distribution information of the arterial part within the preset time period is determined. 4.如权利要求3所述的设备,其特征在于,所述散斑干涉光路单元为空间载波数字散斑干涉光路单元,其中,所述激光器为满足设定相干长度值的激光器,入射至所述图像传感器的参考光为满足设定角度的参考光;所述空间载波数字散斑干涉光路单元还包含光阑;4. The device according to claim 3, wherein the speckle interference optical path unit is a spatial carrier digital speckle interference optical path unit, wherein the laser is a laser that satisfies the set coherence length value, and is incident on the The reference light of the image sensor is the reference light that satisfies the set angle; the space carrier digital speckle interference optical path unit also includes a diaphragm; 所述第一物光入射至所述成像元件包括:The first object light incident on the imaging element includes: 经所述动脉部位反射得到的第一物光入射至所述光阑,并经由所述光阑透射至所述成像元件,其中,所述光阑用于使得所述第一物光的空间频率满足设定范围。The first object light reflected by the arterial part is incident on the diaphragm, and is transmitted to the imaging element through the diaphragm, wherein the diaphragm is used to make the spatial frequency of the first object light meet the setting range. 5.如权利要求3-4任一项所述的设备,其特征在于,所述设备还包括:分析单元,其中:5. The device according to any one of claims 3-4, wherein the device further comprises: an analysis unit, wherein: 分析单元,用于根据所述脉搏空间分布信息和/或所述脉搏时间分布信息,分析所述被测生物体的健康状况。An analysis unit, configured to analyze the health status of the measured organism according to the pulse spatial distribution information and/or the pulse time distribution information.
CN201710211395.3A 2017-03-31 2017-03-31 Method and equipment for measuring pulse distribution Active CN106955092B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710211395.3A CN106955092B (en) 2017-03-31 2017-03-31 Method and equipment for measuring pulse distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710211395.3A CN106955092B (en) 2017-03-31 2017-03-31 Method and equipment for measuring pulse distribution

Publications (2)

Publication Number Publication Date
CN106955092A CN106955092A (en) 2017-07-18
CN106955092B true CN106955092B (en) 2021-02-23

Family

ID=59483936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710211395.3A Active CN106955092B (en) 2017-03-31 2017-03-31 Method and equipment for measuring pulse distribution

Country Status (1)

Country Link
CN (1) CN106955092B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107468222A (en) * 2017-09-15 2017-12-15 王庆亚 High-precision fast-response dynamic 3 D pulse detector
CN110101373A (en) * 2019-06-14 2019-08-09 杭州电子科技大学 A kind of pulses measure method and system based on laser speckle

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003139515A (en) * 2001-11-02 2003-05-14 Fukuoka Prefecture Method for measuring absolute value of deformation quantity using speckle
CN101116087B (en) * 2005-02-08 2011-08-31 皇家飞利浦电子股份有限公司 Optical speckle pattern investigation
AU2007267633B2 (en) * 2006-05-24 2014-07-31 Tarilian Laser Technologies, Limited Optical vital sign detection method and measurement device
WO2009008745A2 (en) * 2007-07-06 2009-01-15 Industrial Research Limited Laser speckle imaging systems and methods
JP5036005B2 (en) * 2009-03-10 2012-09-26 国立大学法人埼玉大学 Minute displacement measuring method and apparatus
CN102221342B (en) * 2011-04-02 2013-04-24 北京交通大学 Method for measuring object deformation by time-domain multi-wavelength heterodyne speckle interference
US10299735B2 (en) * 2014-08-22 2019-05-28 Pulse Tectonics Llc Automated diagnosis based at least in part on pulse waveforms

Also Published As

Publication number Publication date
CN106955092A (en) 2017-07-18

Similar Documents

Publication Publication Date Title
JP5787255B2 (en) Program for correcting measurement data of PS-OCT and PS-OCT system equipped with the program
JP5626687B2 (en) 2-beam optical coherence tomography system
JP5623028B2 (en) Imaging method and apparatus for taking optical coherence tomographic image
JP6948311B2 (en) Full-field interference microscope imaging method and system
US20080025570A1 (en) Dynamic motion contrast and transverse flow estimation using optical coherence tomography
JP6256879B2 (en) Polarization-sensitive optical image measurement system and program installed in the system
JP6230023B2 (en) Apparatus and method for generating a two-dimensional image from three-dimensional optical coherence tomography interferogram data
JP2009165710A (en) Quantitative measurement device for fundus blood flow
US20150216415A1 (en) Method and apparatus for performing multidimensional velocity measurements using amplitude and phase in optical interferometry
CN105595968A (en) Subject information obtaining device and subject information obtaining method
CN109645954B (en) Multi-beam optical coherent elastic measurement system and method based on microlens array
CN105865613A (en) Underwater optical detection and imaging sensing method and system used for ocean stereo monitoring
WO2018212115A1 (en) Device and method for tomographic visualization of tissue viscoelasticity
Qureshi et al. Advances in laser speckle imaging: From qualitative to quantitative hemodynamic assessment
CN106955092B (en) Method and equipment for measuring pulse distribution
Savchenko et al. The use of speckle correlation analysis to determine blood flow velocity
JP2017524138A (en) Instantaneous time-domain optical coherence tomography
JP2010088499A (en) Device and method of processing biological information
Zhang et al. Precise detection of wrist pulse using digital speckle pattern interferometry
Bazan et al. A New High‐Resolution Spectral Approach to Noninvasively Evaluate Wall Deformations in Arteries
Horstmann et al. Optical full-field holographic detection system for non-contact photoacoustic tomography
CN111568385B (en) Diagnostic detection method, device and system for optical coherence tomography equipment
KR102489416B1 (en) Method and apparatus for estimating young's modulus of red blood cell based on speckle pattern of fingernail fold
JP2023059459A (en) Image processing device, image processing method and program
JP2009031230A (en) Displaying measurement data

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
OL01 Intention to license declared
OL01 Intention to license declared