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
light
distribution information
speckle
determining
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 pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • 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 classified in A61B5/00, 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 traditional oriental medicine

Abstract

The application discloses a method and equipment for measuring pulse distribution, wherein the method comprises the following steps: constructing a speckle interference light path, wherein the speckle interference light path comprises an artery part of pulse pulsation of a detected organism; determining pulse distribution information of the artery part according to the speckle interference light path, wherein the pulse distribution information comprises: pulse spatial distribution information and/or pulse temporal distribution information. According to the method and the device, the pulse distribution can be measured, and effective and objective pulse data are provided for clinical diagnosis.

Description

Method and equipment for measuring pulse distribution
Technical Field
The application relates to the technical field of full-field light detection, in particular to a pulse distribution measuring method and device.
Background
Pulse taking in traditional Chinese medicine is also called pulse feeling, which means that traditional Chinese medical doctors can sense pulse by pressing arteries of patients, identify pulse conditions and diagnose diseases according to the pulse conditions.
TCM can identify the pulse condition and analyze it from four aspects, i.e., pulse position, pulse rate, pulse shape and pulse condition. The pulse position refers to the position and length of the pulse, the frequency and rhythm of the pulse, the width of the pulse, and the like. The pulse is the reflection of the human body's function, and different pulse conditions all have corresponding ranges of symptoms.
The human pulse is an unstable vibration with a minute amplitude. Pulse measuring instruments in the prior art, for example: the electronic pulse measurement instrument or the photoelectric pulse measurement instrument can only measure the oscillation frequency information (i.e., the pulse rate) of the pulse, and cannot measure pulse information other than the oscillation frequency information of the pulse.
Disclosure of Invention
In view of this, embodiments of the present application provide a method and an apparatus for measuring pulse distribution, so as to solve the technical problem that pulse information other than the vibration frequency information of a pulse cannot be measured in the prior art, and provide objective pulse distribution data for clinical diagnosis.
The embodiment of the application provides a method for measuring pulse distribution, which comprises the following steps:
constructing a speckle interference light path, wherein the speckle interference light path comprises an artery part of pulse pulsation of a detected organism;
determining pulse distribution information of the artery part according to the speckle interference light path, wherein the pulse distribution information comprises: pulse spatial distribution information and/or pulse temporal distribution information.
Optionally, determining pulse space distribution information of the artery part according to the speckle interference optical path includes:
acquiring a first speckle pattern for the arterial site at a first time instant and acquiring a second speckle pattern 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;
and determining pulse space distribution information of the artery part at the second moment according to the first phase quantity and the second phase quantity.
Optionally, determining pulse time distribution information of the artery part according to the speckle interference optical path includes:
acquiring a plurality of speckle patterns aiming at the artery part within a preset time length;
determining a plurality of phase quantities corresponding to the plurality of speckle patterns;
and determining pulse time distribution information of the artery part in the preset time length according to the phase quantities.
Optionally, the speckle interference optical path is a spatial carrier digital speckle interference optical path, where the spatial carrier digital speckle interference optical path includes a laser, a beam splitting element, a beam expanding element, a diaphragm, an imaging element and an image sensor that satisfy a set coherence length value, and the spatial carrier digital speckle interference optical path is constructed, including:
the laser emits laser light;
the laser light is reflected by the light splitting element, and reference light which meets a set angle is incident to the image sensor;
the laser is transmitted to the beam expanding element through the light splitting element, the beam expanding element expands the laser beam and emits the laser beam to illuminate the artery part, a first object light obtained by reflection of the artery part enters the diaphragm and is transmitted to the imaging element through the diaphragm, and the diaphragm is used for enabling the spatial frequency of the first object light to meet a set range;
the second object light imaged by the imaging element is incident to the image sensor;
the second light and the reference light interfere at the surface of the image sensor.
Optionally, the method further comprises:
and analyzing the health condition of the tested organism according to the pulse space distribution information and/or the pulse time distribution information.
The embodiment of the present application further provides a pulse distribution measuring apparatus, including: speckle interference optical path unit and confirm the unit, wherein:
the speckle interference optical path unit is used for illuminating the pulse beating artery part of the tested organism;
a determining unit, configured to determine pulse distribution information of the artery portion, wherein the pulse distribution information includes: pulse spatial distribution information and pulse temporal distribution information.
The embodiment of the application adopts at least one technical scheme which can achieve the following beneficial effects:
the pulse space distribution information and/or the pulse time distribution information of the artery part are determined according to the speckle interference light path by constructing the speckle interference light path of the artery part containing the pulse pulsation of the tested organism, so that the pulse distribution is measured, and effective and objective pulse data are provided for clinical diagnosis.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application. In the drawings:
fig. 1 is a schematic flowchart illustrating a method for measuring pulse distribution according to an embodiment of the present disclosure;
fig. 2 is a schematic optical path diagram of a spatial carrier digital speckle interference optical path provided in the embodiment of the present application;
FIG. 3(a) is a schematic diagram of a speckle pattern of a target object provided by an embodiment of the present application;
fig. 3(b) is a schematic diagram of a frequency spectrum distribution diagram of a target object according to an embodiment of the present application;
fig. 3(c) is a schematic diagram of a filtered spectrum of a target object according to an embodiment of the present application;
fig. 3(d) is a schematic diagram of a deformed phase diagram of a target object provided in an embodiment of the present application;
fig. 4 is a schematic structural diagram of a pulse distribution measuring apparatus according to an embodiment of the present disclosure.
Detailed Description
In order to achieve the purpose of the present application, the embodiments of the present application provide a method and an apparatus for measuring pulse distribution, which implement measurement of pulse distribution by constructing a speckle interference optical path of an artery portion containing pulse pulsation of a measured organism, and then determining pulse spatial distribution information and/or pulse temporal distribution information of the artery portion according to the speckle interference optical path, so as to provide effective and objective pulse data for clinical diagnosis.
The technical solutions of the present application will be described clearly and completely below with reference to the specific embodiments of the present application and the accompanying drawings. It should be apparent that the described embodiments are only some of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The technical solutions provided by the embodiments of the present application are described in detail below with reference to the accompanying drawings.
Example 1
Fig. 1 is a flowchart illustrating a method for measuring pulse distribution according to an embodiment of the present disclosure. The method may be as follows.
Step 11: and constructing a speckle interference light path.
Wherein, the speckle interference light path comprises the artery part of the pulse beat of the tested organism.
In step 11, a speckle interference optical path of an artery portion including a pulse beat of the measured living body is constructed to prepare for further measurement of pulse distribution according to the speckle interference optical path.
The pulse position of the measured organism represents a position where the pulse of the measured organism is obvious and easy to measure, such as: the radial artery of the human body near the wrist is not particularly limited.
Optionally, the speckle interference optical path is a spatial carrier digital speckle interference optical path.
The pulse beat is unstable vibration with small amplitude, belongs to dynamic deformation, and the dynamic deformation of the pulse can be measured by the space carrier digital speckle interference light path, so that the space distribution and the time distribution corresponding to the pulse are obtained.
The space carrier digital speckle interference optical path generates a space carrier signal by inclining the reference light and the optical axis by a set angle. After the measured speckle pattern is subjected to Fourier transform, a high-frequency part containing deformation information of a target object (an artery part of pulse beat of a measured organism) and a low-frequency part representing background light are obtained, the high-frequency part in a frequency spectrum is subjected to frequency shift due to introduction of a space carrier signal, and separation from the low-frequency part is realized, so that an interference phase is solved by extracting the high-frequency part containing the deformation information of the target object, and finally the deformation amount of the target object is determined according to the relation between the interference phase and the deformation amount of the target object.
It should be noted that the spatial carrier digital speckle interference optical path may be a spatial optical path established in a spatial environment, or may be an optical path integrated in a measurement instrument, and is not specifically limited herein.
Fig. 2 is an optical path schematic diagram of a spatial carrier digital speckle interference optical path provided in the embodiment of the present application.
As shown in fig. 2, optionally, the spatial carrier digital speckle interference optical path includes a laser, a light splitting element, a beam expanding element, a diaphragm, an imaging element and an image sensor that satisfy a set coherence length value, and the spatial carrier digital speckle interference optical path is constructed by:
the laser emits laser;
the laser light is reflected by the light splitting element and enters the image sensor, wherein the reference light meets the set angle;
the laser transmits to the beam expanding element through the light splitting element, the beam expanding element expands the beam of the laser and emits the laser to illuminate an artery part, a first object light obtained by reflection of the artery part is incident to the diaphragm and transmitted to the imaging element through the diaphragm, and the diaphragm is used for enabling the spatial frequency of the first object light to meet a set range;
the second object light imaged by the imaging element is incident to the image sensor;
the second object light interferes with the reference light at the image sensor surface.
Wherein the set angle (α) satisfies the formula:
Figure BDA0001260431890000051
where λ denotes the wavelength of the laser light emitted by the laser, fMRepresenting the spatial frequency, f, of the first object lightPRepresenting the temporal sampling frequency of the image sensor.
It should be noted that the set coherence length value indicates that the laser emitted by the laser can ensure that the reference light interferes with the second light, and the size of the set coherence length value may be determined according to an actual situation, which is not specifically limited herein.
In the embodiment of the application, the time sampling frequency of the image sensor is determined according to the actual situation on the basis of meeting the sampling frequency, and the time sampling frequency is not particularly limited.
In order to determine the pulse distribution at the artery site from the spatial carrier digital speckle interference optical path, the spatial frequency f of the first object light isMThe set range is satisfied:
Figure BDA0001260431890000061
since the diaphragm is used to make the spatial frequency of the first object light satisfy the set range, the aperture size of the diaphragm satisfies:
Figure BDA0001260431890000062
where d denotes the distance of the diaphragm to the imaging element.
Optionally, the light splitting element is an adjustable light splitting flat sheet.
The light splitting element may be an adjustable light splitting flat sheet, an adjustable light splitting mirror, or other optical elements capable of splitting light, and is not limited specifically herein.
Optionally, the beam expanding element is a microscope objective.
Besides the micro objective, the beam expanding element may also be a plano-concave lens, and may also be other optical elements capable of expanding beams, which are not limited herein.
Preferably, the imaging element is a plano-convex lens.
The imaging element may be a plano-convex lens, an industrial lens, or other optical elements capable of imaging, and is not limited in this respect.
In the embodiment of the application, besides the space carrier digital speckle interference light path, the speckle interference light path can also adopt a 4+1 time phase shift digital speckle interference light path, and other speckle interference light paths capable of measuring dynamic deformation can also be adopted, and no specific limitation is made here.
Step 12: and determining pulse distribution information of the artery part according to the speckle interference light path.
Wherein the pulse distribution information includes: pulse spatial distribution information and/or pulse temporal distribution information.
In step 12, pulse space distribution information and/or pulse time distribution information of an artery portion of the pulse beat of the measured living body is determined based on the speckle interference optical path constructed in step 11.
Optionally, pulse spatial distribution information of the arterial site is determined from the speckle interference optical path, in particular, a first step of acquiring a first speckle pattern for the arterial site at a first time instant and a second speckle pattern for the arterial site at a second time instant.
For example, taking a spatial carrier digital speckle interference optical path as an example, since the second object light imaged by the imaging element is obtained after imaging an artery portion, and the second object light and the reference light interfere with each other on the surface of the image sensor, an interference speckle pattern (hereinafter, referred to as speckle pattern) for the artery portion can be acquired by the image sensor.
An image sensor acquires a first speckle pattern and a second speckle pattern for an arterial site at a first time and a second time, respectively.
And secondly, determining a first phase quantity corresponding to the first speckle pattern and a second phase quantity corresponding to the second speckle pattern.
And thirdly, determining pulse space distribution information of the artery part at a second moment according to the first phase quantity and the second phase quantity.
Fig. 3(a) is a schematic diagram of a first speckle pattern of a target object provided in an embodiment of the present application.
Still taking the above-mentioned spatial carrier digital speckle interference optical path as an example, fourier transform is performed on a first speckle pattern (for example, fig. 3(a)) corresponding to a target object (for example, an artery portion) to obtain a spectrum distribution diagram corresponding to the first speckle pattern, and fig. 3(b) is a schematic diagram of the spectrum distribution diagram of the target object provided in the embodiment of the present application.
The first speckle pattern corresponds to a spectral distribution:
I(fx)=FT[I(x)]=A(fx)+B(fx)+C(fx),
wherein, I (f)x) Representing the spectral distribution corresponding to the first speckle pattern, I (x) representing the light intensity of the target object, A (f)x) And B (f)x) Representing two conjugated high-frequency parts, C (f)x) Representing the low frequency part, fxRepresenting the frequency domain coordinates, x the one-dimensional spatial coordinates, and the notation FT the fourier transform.
In the spectral profile:
A(fx)=FT[O(x)R*(x)];
B(fx)=FT[R(x)O*(x)];
C(fx)=FT[|R(x)|2+|O(x)|2]。
where 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 symbol x represents conjugation.
Filtering the obtained signal by a frequency spectrum window to obtain one high-frequency part A (f)x) The other two spectral portions B (f) are filtered outx) And C (f)x) Fig. 3(c) is a schematic diagram of a filtered spectrum of a target object according to an embodiment of the present application.
And performing inverse Fourier transform on the filtered frequency spectrum, and obtaining an expression containing a first phase quantity of the target object according to the following formula:
Figure BDA0001260431890000081
wherein A isIFTAn expression representing a first phase quantity containing a target object,
Figure BDA0001260431890000082
a first phase amount representing the target object, a (x) light intensity information of the interference light of the target object when the first speckle pattern is obtained, and a symbol IFT representing an inverse fourier transform.
Similarly, for the second speckle pattern corresponding to the target object, the expression containing the second phase amount of the target object can be obtained by performing the operations of fourier transform, filtering processing, and inverse fourier transform as described above:
Figure BDA0001260431890000083
wherein A isIFT' represents an expression containing a second phase amount of the target object,
Figure BDA0001260431890000084
a second phase amount representing the target object, a' (x) representing light intensity information of the interference light of the target object when the second speckle pattern is obtained, and the symbol IFT representing an inverse fourier transform.
After obtaining an expression containing a first phase amount of the target object and an expression containing a second phase amount of the target object, deformation information of the target object is obtained through the following formula according to the expression containing the first phase amount of the target object and the expression containing the second phase amount of the target object.
Figure BDA0001260431890000091
Where Im and Re represent the imaginary and real parts, respectively, and the symbol denotes the conjugate.
Fig. 3(d) is a schematic diagram of a deformed phase diagram of a target object provided in an embodiment of the present application.
According to the deformation information of the target object, determining the spatial distribution information of the target object at the second moment, including:
according to the deformation information of the target object, determining the spatial distribution information of the target object at the second moment by the following formula:
Figure BDA0001260431890000092
wherein the content of the first and second substances,
Figure BDA0001260431890000093
spatial distribution information representing the target object, i.e. a deformation distribution of the respective position of the target object at the second moment in time,
Figure BDA0001260431890000094
deformation information representing the target object is displayed on the screen,
Figure BDA0001260431890000095
a displacement-sensitive vector is represented which is,
Figure BDA0001260431890000096
kiunit vector, k, representing the direction of the illuminated artery regionoAnd a unit vector representing the direction of the first object light, and λ represents the wavelength of the laser light emitted by the laser.
In addition, the target object indicates an artery portion where the pulse of the test subject is pulsating, that is, pulse space distribution information of the artery portion of the test subject is determined by the above-described method.
Optionally, pulse time distribution information of the artery portion is determined according to a speckle interference light path, and specifically, in the first step, a plurality of speckle patterns for the artery portion are acquired within a preset time duration.
The image sensor acquires a plurality of speckle patterns for the artery portion within a preset time period.
It should be noted that the number of speckle patterns is related to the preset time duration and the time sampling frequency of the image sensor.
For example: the time sampling frequency of the image sensor is 75 frames/second, the preset time duration is 10 seconds, and the number of the speckle patterns acquired within 10 seconds of the preset time duration is 750.
And secondly, determining a plurality of phase quantities corresponding to the speckle patterns.
Thirdly, determining pulse time distribution information of the artery part in the preset time length according to the plurality of phase quantities.
Still taking the above-mentioned spatial carrier digital speckle interference optical path as an example, the operation of determining the phase amount corresponding to the speckle pattern according to the speckle pattern of the target object as described in the above-mentioned fig. 3(a), fig. 3(b) and fig. 3(c) by the spatial carrier method determines a plurality of phase amounts corresponding to a plurality of speckle patterns acquired for the arterial site; and determining a plurality of pulse shape information of the artery region from any two adjacent phase amounts of the plurality of phase amounts by the operation of determining the shape change information of the target object from the first phase amount and the second phase amount described above.
For example, a 4-radiation speckle pattern is acquired for the artery part A within a preset time period4 phase quantities a, b, c and d are determined for 4 speckle patterns. According to the phase quantity a and the phase quantity b, determining first pulse deformation information of the artery part
Figure BDA0001260431890000101
Determining second pulse deformation information of the artery part according to the phase quantity b and the phase quantity c
Figure BDA0001260431890000102
Determining the third pulse deformation information of the artery part according to the phase quantity c and the phase quantity d
Figure BDA0001260431890000103
It should be noted that the method for determining the plurality of phase quantities corresponding to the plurality of speckle patterns may adopt a spatial carrier method, a 4+1 time phase shift method, or other methods capable of measuring dynamic deformation, which is not specifically limited herein.
For example, when a plurality of phase quantities are determined by adopting a 4+1 time phase shift method, firstly, 4 initial speckle patterns aiming at an artery part are acquired by a four-step phase shift method; secondly, collecting a plurality of speckle patterns aiming at the artery part within a preset time length; and finally, determining a plurality of phase quantities corresponding to the plurality of speckle patterns through a time phase shift algorithm according to any one of the plurality of speckle patterns and the 4 initial speckle patterns.
The plurality of pulse wave deformation distributions at the artery region are determined by the method for determining the deformation distribution of the target object based on the deformation information of the target object described in fig. 3 (d).
Taking the artery part a as an example, determining a first pulse deformation distribution according to the first pulse deformation information; determining second pulse deformation distribution according to the second pulse deformation information; and determining a third pulse deformation distribution according to the third pulse deformation information.
The first pulse deformation distribution, the second pulse deformation distribution and the third pulse deformation distribution jointly form pulse time distribution information of the artery part within preset time.
In an alternative embodiment of the present application, the method further comprises:
and analyzing the health condition of the tested organism according to the pulse space distribution information and/or the pulse time distribution information.
Obtaining the pulse position of the pulse according to the pulse space distribution information of the artery part of the tested organism; obtaining the pulse number of the pulse according to the pulse time distribution information of the artery part of the tested organism; the pulse shape and pulse condition of the pulse can be obtained by integrating the pulse space distribution information and the pulse time distribution information of the artery part of the tested organism. The health condition of the tested organism can be analyzed according to the pulse condition information of pulse position, pulse number, pulse shape, pulse condition and the like.
According to the technical scheme, the speckle interference light path of the artery part containing the pulse pulsation of the detected organism is constructed, and then the pulse space distribution information and/or the pulse time distribution information of the artery part are determined according to the speckle interference light path, so that the pulse distribution is measured, and effective and objective pulse data are provided for clinical diagnosis.
Example 2
Fig. 4 is a schematic structural diagram of a pulse deformation measurement apparatus according to an embodiment of the present application, where the apparatus 400 includes: a speckle interference optical path unit 401 and a determination unit 402, wherein:
a speckle interference optical path unit 401 for illuminating the pulse beating artery part of the tested organism;
a determining unit 402, configured to determine pulse distribution information of the artery portion, where the pulse distribution information includes: pulse spatial distribution information and/or pulse temporal distribution information.
Optionally, the determining unit 402 determines pulse space distribution information of the artery region, including:
acquiring a first speckle pattern for the arterial site at a first time instant and acquiring a second speckle pattern 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;
and determining pulse space distribution information of the artery part at a second moment according to the first phase quantity and the second phase quantity.
Optionally, the determining unit 402 determines pulse time distribution information of the artery region, including:
acquiring a plurality of speckle patterns aiming at an artery part within a preset time length;
determining a plurality of phase quantities corresponding to the plurality of speckle patterns;
and determining pulse time distribution information of the artery part within a preset time length according to the plurality of phase quantities.
Optionally, the speckle interference optical path unit 401 is a spatial carrier digital speckle interference optical path unit, where the spatial carrier digital speckle interference optical path unit includes a laser, a beam splitting element, a beam expanding element, a diaphragm, an imaging element and an image sensor, which satisfy a set coherence length value, where:
the laser emits laser;
the laser light is reflected by the light splitting element and enters the image sensor, wherein the reference light meets the set angle;
the laser transmits to the beam expanding element through the light splitting element, the beam expanding element expands the beam of the laser and emits the laser to illuminate an artery part, a first object light obtained by reflection of the artery part is incident to the diaphragm and transmitted to the imaging element through the diaphragm, and the diaphragm is used for enabling the spatial frequency of the first object light to meet a set range;
the second object light imaged by the imaging element is incident to the image sensor;
the second object light interferes with the reference light at the image sensor surface.
Optionally, the apparatus 400 further comprises: an analysis unit, wherein:
and the analysis unit is used for analyzing the health condition of the tested organism according to the pulse space distribution information and/or the pulse time distribution information.
The speckle interference optical path unit is used for illuminating an artery part of pulse pulsation of a tested organism according to the measuring equipment of the pulse distribution; the determining unit is used for determining pulse space distribution information and/or pulse time distribution information of the artery part, realizing measurement of pulse distribution and providing effective and objective pulse data for clinical diagnosis.
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, and the like) 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 of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto 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 which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, 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 Discs (DVD) or other optical storage, magnetic cassettes, 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, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or 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, and the like) having computer-usable program code embodied therein.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims (5)

1. A method for measuring pulse distribution, comprising:
constructing a speckle interference optical path, wherein the speckle interference optical path comprises an artery part of pulse pulsation of a detected organism, a laser, a light splitting element, a beam expanding element, an imaging element and an image sensor; constructing the speckle interference light path comprises the following steps:
the laser emits laser light;
the laser light is reflected by the light splitting element to obtain reference light, and the reference light is incident to the image sensor;
the laser is transmitted to the beam expanding element through the light splitting element, the beam expanding element expands the laser beam and emits the laser beam to illuminate the artery part, the laser beam is reflected by the artery part to obtain first object light, the first object light enters the imaging element, and second object light obtained by imaging of the imaging element enters the image sensor;
the second light and the reference light interfere on the surface of the image sensor;
determining pulse distribution information of the artery part according to the speckle interference light path, wherein the pulse distribution information comprises: pulse spatial distribution information and pulse temporal distribution information;
determining pulse space distribution information of the artery part according to the speckle interference light path, wherein the determining comprises the following steps:
acquiring a first speckle pattern for the arterial site at a first time instant and acquiring a second speckle pattern 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;
determining pulse space distribution information of the artery part at the second moment according to the first phase quantity and the second phase quantity;
determining pulse time distribution information of the artery part according to the speckle interference light path, wherein the pulse time distribution information comprises the following steps:
acquiring a plurality of speckle patterns aiming at the artery part within a preset time length;
determining a plurality of phase quantities corresponding to the plurality of speckle patterns;
and determining pulse time distribution information of the artery part in the preset time length according to the phase quantities.
2. The method of 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 the reference light incident to the image sensor is a reference light satisfying a set angle; the space carrier digital speckle interference optical path also comprises a diaphragm;
the incident of the first object light to the imaging element includes:
the first object light obtained by reflection of the artery part enters the diaphragm and is transmitted to the imaging element through the diaphragm, wherein the diaphragm is used for enabling the spatial frequency of the first object light to meet a set range.
3. A device for measuring the distribution of a pulse, characterized in that it comprises: speckle interference optical path unit and confirm the unit, wherein:
the speckle interference optical path unit is used for illuminating the pulse beating artery part of the tested organism;
the speckle interference optical path unit comprises a laser, a light splitting element, a beam expanding element, an imaging element and an image sensor; constructing the speckle interference light path comprises the following steps:
the laser light is reflected by the light splitting element to obtain reference light, and the reference light is incident to the image sensor;
the laser is transmitted to the beam expanding element through the light splitting element, the beam expanding element expands the laser beam and emits the laser beam to illuminate the artery part, the laser beam is reflected by the artery part to obtain first object light, the first object light enters the imaging element, and second object light obtained by imaging of the imaging element enters the image sensor;
the second light and the reference light interfere on the surface of the image sensor;
a determining unit, configured to determine pulse distribution information of the artery portion, wherein the pulse distribution information includes: pulse spatial distribution information and pulse temporal distribution information;
wherein the determining unit determines pulse space distribution information of the artery part, including:
acquiring a first speckle pattern for the arterial site at a first time instant and acquiring a second speckle pattern 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;
determining pulse space distribution information of the artery part at the second moment according to the first phase quantity and the second phase quantity;
wherein the determining unit determines pulse time distribution information of the artery part, including:
acquiring a plurality of speckle patterns aiming at the artery part within a preset time length;
determining a plurality of phase quantities corresponding to the plurality of speckle patterns;
and determining the time distribution information of the artery part in the preset time length according to the plurality of phase quantities.
4. The apparatus 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 satisfying a set coherence length value, and the reference light incident to the image sensor is a reference light satisfying a set angle; the space carrier digital speckle interference light path unit also comprises a diaphragm;
the incident of the first object light to the imaging element includes:
the first object light obtained by reflection of the artery part enters the diaphragm and is transmitted to the imaging element through the diaphragm, wherein the diaphragm is used for enabling the spatial frequency of the first object light to meet a set range.
5. The apparatus of any of claims 3-4, wherein the apparatus further comprises: an analysis unit, wherein:
and the analysis unit is used for analyzing the health condition of the tested organism according to the pulse space 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
WO2006085252A2 (en) * 2005-02-08 2006-08-17 Koninklijke Philips Electronics N.V. Optical speckle pattern investigation
US8343063B2 (en) * 2006-05-24 2013-01-01 Tarilian Laser Technologies, Limited Optical vital sign detection method and measurement device
EP2166934A4 (en) * 2007-07-06 2012-10-17 Ind Res Ltd 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
KR20170087855A (en) * 2014-08-22 2017-07-31 펄스 테크토닉스 엘엘씨 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
KR102225808B1 (en) Skin diagnosis device, skin condition output method, program, and recording medium
JP5166889B2 (en) Quantitative measurement device for fundus blood flow
Thompson et al. Tissue perfusion measurements: multiple-exposure laser speckle analysis generates laser Doppler–like spectra
JP6948311B2 (en) Full-field interference microscope imaging method and system
JP6550745B2 (en) Blood flow measuring device
Gurov et al. High-speed video capillaroscopy method for imaging and evaluation of moving red blood cells
JP6230023B2 (en) Apparatus and method for generating a two-dimensional image from three-dimensional optical coherence tomography interferogram data
US10524663B2 (en) Phase measurement, analysis, and correction methods for coherent imaging systems
JP7154542B2 (en) Apparatus and method for tomographic visualization of tissue viscoelasticity
CN106955092B (en) Method and equipment for measuring pulse distribution
Campo et al. Non-invasive technique for assessment of vascular wall stiffness using laser Doppler vibrometry
Lam Po Tang et al. Non-contact quantification of jugular venous pulse waveforms from skin displacements
Savchenko et al. The use of speckle correlation analysis to determine blood flow velocity
Silva et al. Signal-carrying speckle in optical coherence tomography: a methodological review on biomedical applications
US11625003B2 (en) Coherent optical holographic imaging
Sharath et al. Radial arterial compliance measurement by fiber Bragg grating pulse recorder
Rao et al. Imaging pulsatile retinal blood flow in human eye
JP6648891B2 (en) Apparatus and method for tomographically visualizing substance content
US11957441B2 (en) System and method for depth flow inspection
RU2626310C2 (en) Method of visualizing object areas containing micromotions
Di Lascio et al. Non-invasive assessment of carotid PWV via accelerometric sensors: validation of a new device and comparison with established techniques
US20180018785A1 (en) Processing optical coherence tomography scans
US20190099083A1 (en) Information processing apparatus and information processing method
Horstmann et al. Optical full-field holographic detection system for non-contact photoacoustic tomography
KR102489416B1 (en) Method and apparatus for estimating young's modulus of red blood cell based on speckle pattern of fingernail fold

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