CN106955092A - A kind of measuring method and equipment of pulse distribution - Google Patents

A kind of measuring method and equipment of pulse distribution Download PDF

Info

Publication number
CN106955092A
CN106955092A CN201710211395.3A CN201710211395A CN106955092A CN 106955092 A CN106955092 A CN 106955092A CN 201710211395 A CN201710211395 A CN 201710211395A CN 106955092 A CN106955092 A CN 106955092A
Authority
CN
China
Prior art keywords
pulse
speckle
artery site
distribution information
artery
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.)
Granted
Application number
CN201710211395.3A
Other languages
Chinese (zh)
Other versions
CN106955092B (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

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

This application discloses the measuring method and equipment of a kind of pulse distribution, methods described includes:Speckle interference light path is built, wherein, the artery site of the beat pulse comprising tested organism in the speckle interference light path;According to the speckle interference light path, the pulse distributed intelligence of the artery site is determined, wherein, the pulse distributed intelligence includes:Pulse space distribution information and/or pulse Annual distribution information.According to the method and apparatus of the embodiment of the present application, the measurement being distributed to pulse can be realized, effectively objective pulse data is provided for clinical diagnosis.

Description

A kind of measuring method and equipment of pulse distribution
Technical field
The application is related to whole audience technical field of light detection, more particularly to the measuring method and equipment that a kind of pulse is distributed.
Background technology
Being felt the pulse in traditional Chinese medicine also known as feel the pulse, and refers to traditional Chinese physician by pressing the artery of patient to perceive beat pulse, Pulse condition is recognized, and then is diagnosed the illness according to pulse condition.
Traditional Chinese physician recognizes pulse condition, is often analyzed in terms of arteries and veins position, rapid pulse, arteries and veins shape, arteries and veins gesture four.Arteries and veins position refers to beat pulse The position shown and length, rapid pulse refer to the number of times and the rhythm and pace of moving things of beat pulse, and arteries and veins shape refers to the forms such as the width of beat pulse, and arteries and veins gesture refers to The trend such as power, the smoothness of beat pulse.Pulse condition is the reflection to function of human body, different pulse condition all corresponding illness models Enclose.
Human pulse is a kind of small unstable vibration of amplitude.Electronic sphygmograph of the prior art, for example:Electronic type Electronic sphygmograph or photo-electric electronic sphygmograph, can only be measured, it is impossible to survey to the vibration frequency information (i.e. rapid pulse) of pulse Measure other pulse informations in addition to the vibration frequency information of pulse.
The content of the invention
In view of this, the embodiment of the present application provides a kind of measuring method and equipment of pulse distribution, for solving existing skill The technical problem of other pulse informations in addition to the vibration frequency information of pulse can not be measured in art, visitor is provided for clinical diagnosis The pulse distributed data of sight.
The embodiment of the present application provides a kind of measuring method of pulse distribution, including:
Speckle interference light path is built, wherein, the dynamic of the beat pulse for being tested organism is included in the speckle interference light path Arteries and veins position;
According to the speckle interference light path, the pulse distributed intelligence of the artery site is determined, wherein, the pulse distribution Information includes:Pulse space distribution information and/or pulse Annual distribution information.
Alternatively, according to the speckle interference light path, the pulse space distribution information of the artery site is determined, including:
The first speckle pattern for the artery site is gathered at the first moment, and in the collection of the second moment for described Second speckle pattern of artery site;
The corresponding first phase amount of first speckle pattern is determined, and determines corresponding second phase of second speckle pattern Position amount;
According to the first phase amount and the second phase amount, arteries and veins of the artery site at second moment is determined Fight space distribution information.
Alternatively, according to the speckle interference light path, the pulse Annual distribution information of the artery site is determined, including:
Several speckle patterns of collection for the artery site in preset duration;
It is determined that the corresponding multiple phase masses of described several speckle patterns;
According to the multiple phase mass, pulse Annual distribution letter of the artery site in the preset duration is determined Breath.
Alternatively, the speckle interference light path is spatial carrier digital speckle interference light path, wherein, the spatial carrier number Laser of the word speckle interference light path comprising satisfaction setting coherence length value, beam splitter, beam-expanding element, diaphragm, image-forming component And imaging sensor, the spatial carrier digital speckle interference light path is built, including:
The laser launches laser;
The laser reflects via the beam splitter to be obtained and meets the reference light of set angle and be incident to the figure As sensor;
The laser is transmitted through the beam-expanding element via the beam splitter, and the beam-expanding element is by the laser beam expanding The artery site is illuminated after outgoing, reflecting obtained the first object light through the artery site is incident to the diaphragm, and via The diaphragm is transmitted through the image-forming component, wherein, the diaphragm is used to cause that the spatial frequency satisfaction of first object light to be set Determine scope;
The second object light that the image-forming component imaging is obtained is incident to described image sensor;
Second object light is interfered with the reference light on described image sensor surface.
Alternatively, methods described also includes:
According to the pulse space distribution information and/or the pulse Annual distribution information, the tested organism is analyzed Health status.
The embodiment of the present application also provides a kind of measuring apparatus of pulse distribution, including:Speckle interference optical path unit and determination Unit, wherein:
Speckle interference optical path unit, the artery site of the beat pulse for illuminating tested organism;
Determining unit, the pulse distributed intelligence for determining the artery site, wherein, the pulse distributed intelligence bag Include:Pulse space distribution information and pulse Annual distribution information.
At least one above-mentioned technical scheme that the embodiment of the present application is used can reach following beneficial effect:
By the speckle interference light path for the artery site for building the beat pulse comprising tested organism, and then according to described Speckle interference light path, determines the pulse space distribution information and/or pulse Annual distribution information of the artery site, realizes to arteries and veins Fight the measurement of distribution, effectively objective pulse data is provided for clinical diagnosis.
Brief description of the drawings
Accompanying drawing described herein is used for providing further understanding of the present application, constitutes the part of the application, this Shen Schematic description and description please is used to explain the application, does not constitute the improper restriction to the application.In the accompanying drawings:
A kind of schematic flow sheet of the measuring method for pulse distribution that Fig. 1 provides for the embodiment of the present application;
The light path schematic diagram for the spatial carrier digital speckle interference light path that Fig. 2 provides for the embodiment of the present application;
The schematic diagram of one width speckle pattern of the destination object that Fig. 3 (a) provides for the embodiment of the present application;
The schematic diagram of one width frequency spectrum profile of the destination object that Fig. 3 (b) provides for the embodiment of the present application;
Fig. 3 (c) for the filtered frequency spectrum of destination object that the embodiment of the present application is provided schematic diagram;
The schematic diagram of the deformation phase diagram for the destination object that Fig. 3 (d) provides for the embodiment of the present application;
A kind of structural representation of the measuring apparatus for pulse distribution that Fig. 4 provides for the embodiment of the present application.
Embodiment
In order to realize the purpose of the application, the embodiment of the present application provides the measuring method and equipment that a kind of pulse is distributed, By the speckle interference light path for the artery site for building the beat pulse comprising tested organism, and then according to the speckle interference Light path, determines the pulse space distribution information and/or pulse Annual distribution information of the artery site, realizes what pulse was distributed Measurement, effectively objective pulse data is provided for clinical diagnosis.
Technical scheme is clearly and completely retouched with reference to the application specific embodiment and corresponding accompanying drawing State.Obviously, described embodiment is only some embodiments of the present application, rather than whole embodiments.Based in the application Embodiment, the every other embodiment that those of ordinary skill in the art are obtained under the premise of creative work is not made, Belong to the scope of the application protection.
Below in conjunction with accompanying drawing, the technical scheme that each embodiment of the application is provided is described in detail.
Embodiment 1
A kind of schematic flow sheet of the measuring method for pulse distribution that Fig. 1 provides for the embodiment of the present application.Methods described can With as follows.
Step 11:Build speckle interference light path.
Wherein, the artery site of the beat pulse comprising tested organism in speckle interference light path.
In a step 11, the speckle interference light path of the artery site of the beat pulse comprising tested organism is built, to enter The measurement that one step carries out pulse distribution according to the speckle interference light path is prepared.
Wherein, it is tested organism artery site and represents that the beat pulse of the tested organism is obvious, easy measuring point, example Such as:The radial artery of the nearly wrist portion of human body, is not specifically limited here.
Alternatively, speckle interference light path is spatial carrier digital speckle interference light path.
Beat pulse is a kind of small unstable vibration of amplitude, belongs to dynamic deformation, spatial carrier digital speckle interference Light path can be measured to the dynamic deformation of pulse, and then obtain the corresponding spatial distribution of pulse and Annual distribution.
Spatial carrier digital speckle interference light path is believed by the way that reference light and inclined light shaft set angle are produced into spatial carrier Number.After the speckle pattern obtained to measurement implements Fourier transformation, obtain comprising destination object (the pulse jump of tested organism Dynamic artery site) deformation data HFS and represent the low frequency part of bias light, due to the introducing of spatial carrier signal The HFS in frequency spectrum is occurred frequency displacement, realize the separation with low frequency part, so including destination object deformation by extracting The HFS of information and then interferometric phase is solved, finally according to the relation between interferometric phase and destination object deformation quantity, really Set the goal the deformation quantity of object.
It should be noted that spatial carrier digital speckle interference light path can be the spatial light built in space environment Road or the light path being integrated in measuring instrument, are not specifically limited here.
The light path schematic diagram for the spatial carrier digital speckle interference light path that Fig. 2 provides for the embodiment of the present application.
As shown in Fig. 2 alternatively, spatial carrier digital speckle interference light path includes the laser for meeting setting coherence length value Device, beam splitter, beam-expanding element, diaphragm, image-forming component and imaging sensor, build the spatial carrier digital speckle interference light Road, including:
Laser launches laser;
Laser reflects via beam splitter to be obtained and meets the reference light of set angle and be incident to imaging sensor;
Laser is transmitted through beam-expanding element via beam splitter, and beam-expanding element will illuminate artery portion after the laser beam expanding outgoing Position, reflects the first obtained object light through artery site and is incident to diaphragm, and is transmitted through image-forming component via diaphragm, wherein, the light Door screen is used to make it that the spatial frequency of first object light meets setting range;
The second object light that image-forming component imaging is obtained is incident to imaging sensor;
Second object light is interfered with reference light in image sensor surface.
Wherein, set angle (α) meets formula:
Wherein, λ represents that laser launches the wavelength of laser, fMRepresent the spatial frequency of the first object light, fPRepresent image sensing The time sampling frequency of device.
It should be noted that setting coherence length value represents that the laser of laser emitting ensure that reference light and the second thing Light is interfered, and the size of the setting coherence length value can be determined according to actual conditions, be not specifically limited here.
In the embodiment of the present application, the time sampling frequency of imaging sensor specifically takes on the basis of sample frequency is met Value can be determined according to actual conditions, be not specifically limited here.
It should be noted that in order to according to spatial carrier digital speckle interference light path, determine the pulse of artery site Distribution, the spatial frequency f of the first object lightMMeet setting range:
Because diaphragm is used to make it that the spatial frequency of the first object light meets setting range, therefore, the aperture size of the diaphragm Meet:
Wherein, d represents the diaphragm to the distance of image-forming component.
Alternatively, beam splitter is adjustable light splitting plain film.
Beam splitter is in addition to it can use adjustable light splitting plain film, it would however also be possible to employ adjustable light splitting mirror, can also use Other can carry out the optical element of light splitting, be not specifically limited here.
Alternatively, beam-expanding element is microcobjective.
Beam-expanding element is in addition to it can use the microcobjective, it would however also be possible to employ plano-concave lens, can also use it The optical element that he can be expanded, is not specifically limited here.
Preferably, image-forming component is planoconvex spotlight.
Image-forming component is in addition to it can use planoconvex spotlight, it would however also be possible to employ industrial lens, can also use other energy Enough optical elements being imaged, are not specifically limited here.
In the embodiment of the present application, the speckle interference light path except can using spatial carrier digital speckle interference light path it Outside, 4+1 time-phase displacement digital speckle interference light paths can also be used, the speckle of dynamic deformation can also can be measured using other Optical interference circuit, is not specifically limited here.
Step 12:According to speckle interference light path, the pulse distributed intelligence of the artery site is determined.
Wherein, pulse distributed intelligence includes:Pulse space distribution information and/or pulse Annual distribution information.
In step 12, according to the speckle interference light path built in step 11, it is determined that the beat pulse of tested organism The pulse space distribution information and/or pulse Annual distribution information of artery site.
Alternatively, according to speckle interference light path, the pulse space distribution information of artery site is determined, specifically, the first step, The first speckle pattern for artery site is gathered at the first moment, and is dissipated in the collection of the second moment for the second of artery site Spot figure.
For example, by taking spatial carrier digital speckle interference light path as an example, being because image-forming component is imaged the second obtained object light What is obtained afterwards is imaged to artery site, second object light is interfered with reference light in image sensor surface, therefore, The interference speckle pattern (hereinafter referred to as speckle pattern) for the artery site can be collected by the imaging sensor.
Imaging sensor is directed to the first speckle pattern and second of artery site at the first moment and the collection of the second moment respectively Speckle pattern.
Second step, determines the corresponding first phase amount of the first speckle pattern, and determine corresponding second phase of the second speckle pattern Position amount.
3rd step, according to first phase amount and second phase amount, determines artery site in the pulse space at the second moment point Cloth information.
The schematic diagram of first speckle pattern of the destination object that Fig. 3 (a) provides for the embodiment of the present application.
Still by taking above-mentioned spatial carrier digital speckle interference light path as an example, to destination object (for example:Artery site) it is corresponding First speckle pattern (for example, Fig. 3 (a)) carries out Fourier transformation, obtains the corresponding frequency spectrum profile of the first speckle pattern, Fig. 3 (b) The schematic diagram of the frequency spectrum profile of the destination object provided for the embodiment of the present application.
The corresponding spectrum distribution of first speckle pattern is:
I(fx)=FT [I (x)]=A (fx)+B(fx)+C(fx),
Wherein, I (fx) the corresponding spectrum distribution of the first speckle pattern is represented, I (x) represents the luminous intensity of destination object, A (fx) With B (fx) represent two HFSs being conjugated, C (fx) represent low frequency part, fxFrequency domain coordinates are represented, x represents that the one-dimensional space is sat Mark, symbol FT represents Fourier transform.
In frequency spectrum profile:
A(fx)=FT [O (x) R*(x)];
B(fx)=FT [R (x) O*(x)];
C(fx)=FT [| R (x) |2+|O(x)|2]。
Wherein, R (x) represents the luminous intensity of the object light of destination object, and O (x) represents the luminous intensity of reference light, and symbol * is represented Conjugation.
Filtered by spectral window, obtain one of HFS A (fx), filter off other two portions of the spectrum B (fx) and C(fx), Fig. 3 (c) for the filtered frequency spectrum of destination object that the embodiment of the present application is provided schematic diagram.
Inverse Fourier transform is carried out to frequency spectrum after filtering, according to following formula, the first phase containing destination object is obtained The expression formula of amount:
Wherein, AIFTThe expression formula of the first phase amount containing destination object of expression,Represent the first phase of destination object Amount, A (x) represents the intensity signal of interference light of the destination object when obtaining the first speckle pattern, and symbol IFT represents Fourier inversion Change.
Similarly, for corresponding second speckle pattern of destination object, it again may be by performing above-mentioned Fourier transformation, filtering Processing and the operation of inverse Fourier transform, obtain the expression formula of the second phase amount containing destination object:
Wherein, AIFT' represent the second phase amount containing destination object expression formula,Represent the second phase of destination object Position amount, A'(x) intensity signal of interference light of the destination object when obtaining the second speckle pattern is represented, symbol IFT represents Fourier Inverse transformation.
Obtaining the table of the expression formula of the first phase amount containing destination object and second phase amount containing destination object Up to after formula, according to the expression of the expression formula of the first phase amount containing destination object and the second phase amount containing destination object Formula, by following formula, obtains the deformation data of destination object.
Wherein, Im and Re represent imaginary part and real part respectively, and symbol * represents conjugation.
The schematic diagram of the deformation phase diagram for the destination object that Fig. 3 (d) provides for the embodiment of the present application.
According to the deformation data of destination object, space distribution information of the destination object at the second moment is determined, including:
According to the deformation data of destination object, determine that spatial distribution of the destination object at the second moment is believed by below equation Breath:
Wherein,Represent the space distribution information of destination object, i.e. deformation of each position of destination object at the second moment Distribution,The deformation data of destination object is represented,The sensitive vector of displacement is represented,kiArtery site is illuminated in expression The unit vector in direction, koThe unit vector in the first object light direction is represented, λ represents that laser launches the wavelength of laser.
It should be noted that destination object represents the artery site of the beat pulse of tested organism, i.e., using above-mentioned note The scheme of load, it is determined that the pulse space distribution information of the artery site of tested organism.
Alternatively, according to speckle interference light path, the pulse Annual distribution information of artery site is determined, specifically, the first step, Several speckle patterns of collection for artery site in preset duration.
Imaging sensor gathers several speckle patterns for artery site in preset duration.
It should be noted that the quantity of several speckle patterns is relevant with preset duration, the time sampling frequency of imaging sensor.
For example:The time sampling frequency of imaging sensor is 75 frames/second, and preset duration was 10 seconds, in preset duration 10 seconds The quantity of several speckle patterns inside collected is 750.
Second step, determines the corresponding multiple phase masses of several speckle patterns.
3rd, according to multiple phase masses, determine pulse Annual distribution information of the artery site in the preset duration.
Still by taking above-mentioned spatial carrier digital speckle interference light path as an example, by spatial carrier method, such as above-mentioned Fig. 3 (a), Fig. 3 (b) speckle pattern according to destination object recorded with Fig. 3 (c), determines the operation of the corresponding phase mass of speckle pattern, it is determined that for dynamic The corresponding multiple phase masses of several speckle patterns that arteries and veins position is collected;And, by above-mentioned record according to first phase amount With second phase amount, the operation of the deformation data of destination object is determined, two phases of arbitrary neighborhood in multiple phase masses Amount, determines multiple pulse deformation datas of artery site.
For example, collecting 4 spoke speckle patterns for artery site A in preset duration, 4 width speckle patterns corresponding 4 are determined Individual phase mass a, b, c and d.According to phase mass a and phase mass b, the first pulse deformation data of artery site is determinedAccording to Phase mass b and phase mass c, determines the second pulse deformation data of artery siteAccording to phase mass c and phase mass d, it is determined that 3rd pulse deformation data of artery site
It should be noted that determining the method for the corresponding multiple phase masses of several speckle patterns, spatial carrier method can be used, 4+1 temporal phase-shifting methods can be used, the method that can also can measure dynamic deformation using other is not specifically limited here.
For example, when determining multiple phase masses using 4+1 temporal phase-shifting methods, first, by four-stepped switching policy collection for dynamic The initial speckle pattern of 4 width at arteries and veins position;Secondly, several speckle patterns of collection for artery site in preset duration;Finally, according to Any width and the initial speckle pattern of 4 width in several speckle patterns, passage time phase shift algorithm determine that several speckle patterns are corresponding Multiple phase masses.
The deformation data according to destination object recorded by above-mentioned Fig. 3 (d), determines the side of the deformation distribution of destination object Method, determines multiple pulse deformation distribution of artery site.
Still by taking above-mentioned artery site A as an example, according to the first pulse deformation data, determine that the first pulse deformation is distributed;According to Second pulse deformation data, determines that the second pulse deformation is distributed;According to the 3rd pulse deformation data, the 3rd pulse deformation point is determined Cloth.
First pulse deformation distribution, the second pulse deformation distribution and the 3rd pulse deformation distribution together constitute artery site Pulse Annual distribution information in preset duration.
In the alternative embodiment of the application, methods described also includes:
According to the pulse space distribution information and/or the pulse Annual distribution information, the tested organism is analyzed Health status.
According to the pulse space distribution information of the artery site of tested organism, the arteries and veins position of pulse can be obtained;According to quilt The pulse Annual distribution information of the artery site of organism is surveyed, the rapid pulse of pulse can be obtained;The artery of comprehensive tested organism The pulse space distribution information and pulse Annual distribution information at position, can obtain the arteries and veins shape and arteries and veins gesture of pulse.According to pulse The pulse condition informations such as arteries and veins position, rapid pulse, arteries and veins shape, arteries and veins gesture, can analyze the health status of tested organism.
Technical scheme described in the embodiment of the present application, by the artery portion for building the beat pulse comprising tested organism The speckle interference light path of position, and then according to the speckle interference light path, determine the pulse space distribution information of the artery site And/or pulse Annual distribution information, the measurement being distributed to pulse is realized, effectively objective pulse data is provided for clinical diagnosis.
Embodiment 2
A kind of structural representation of the measuring apparatus for pulse deformation that Fig. 4 provides for the embodiment of the present application, equipment 400 is wrapped Include:Speckle interference optical path unit 401 and determining unit 402, wherein:
Speckle interference optical path unit 401, the artery site of the beat pulse for illuminating tested organism;
Determining unit 402, the pulse distributed intelligence for determining artery site, wherein, pulse distributed intelligence includes:Pulse Space distribution information and/or pulse Annual distribution information.
Optionally it is determined that unit 402 determines the pulse space distribution information of artery site, including:
The first speckle pattern for artery site is gathered at the first moment, and collection is directed to artery site at the second moment The second speckle pattern;
The corresponding first phase amount of the first speckle pattern is determined, and determines the corresponding second phase amount of the second speckle pattern;
According to first phase amount and second phase amount, pulse space distribution information of the artery site at the second moment is determined.
Optionally it is determined that unit 402 determines the pulse Annual distribution information of artery site, including:
Several speckle patterns of collection for artery site in preset duration;
Determine the corresponding multiple phase masses of several speckle patterns;
According to multiple phase masses, pulse Annual distribution information of the artery site in preset duration is determined.
Alternatively, speckle interference optical path unit 401 is spatial carrier digital speckle interference optical path unit, wherein, the space The laser for meeting setting coherence length value, beam splitter, beam-expanding element, light are included in carrier wave digital speckle interference optical path unit Door screen, image-forming component and imaging sensor, wherein:
Laser launches laser;
Laser reflects via beam splitter to be obtained and meets the reference light of set angle and be incident to imaging sensor;
Laser is transmitted through beam-expanding element via beam splitter, and beam-expanding element will illuminate artery site after laser beam expanding outgoing, The first obtained object light is reflected through artery site and is incident to diaphragm, and is transmitted through image-forming component via diaphragm, wherein, diaphragm is used for So that the spatial frequency of the first object light meets setting range;
The second object light that image-forming component imaging is obtained is incident to imaging sensor;
Second object light is interfered with reference light in image sensor surface.
Alternatively, equipment 400 also includes:Analytic unit, wherein:
Analytic unit, for according to pulse space distribution information and/or pulse Annual distribution information, analysis to be tested organism Health status.
The measuring apparatus being distributed according to the pulse, the speckle interference optical path unit is used for the arteries and veins for illuminating tested organism Fight the artery site of bounce;The determining unit is used for the pulse space distribution information and/or pulse for determining the artery site Annual distribution information, realizes the measurement being distributed to pulse, and effectively objective pulse data is provided for clinical diagnosis.
It should be understood by those skilled in the art that, embodiments of the invention can be provided as method, system or computer program Product.Therefore, the present invention can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware Apply the form of example.Moreover, the present invention can be used in one or more computers for wherein including computer usable program code The computer program production that usable storage medium is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.) The form of product.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product Figure and/or block diagram are described.It should be understood that every one stream in flow chart and/or block diagram can be realized by computer program instructions Journey and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided The processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce A raw machine so that produced by the instruction of computer or the computing device of other programmable data processing devices for real The device for the function of being specified in present one flow of flow chart or one square frame of multiple flows and/or block diagram or multiple square frames.
These computer program instructions, which may be alternatively stored in, can guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works so that the instruction being stored in the computer-readable memory, which is produced, to be included referring to Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one square frame of block diagram or The function of being specified in multiple square frames.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that in meter Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented processing, thus in computer or The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in individual square frame or multiple square frames.
In a typical configuration, computing device includes one or more processors (CPU), input/output interface, net Network interface and internal memory.
Internal memory potentially includes the volatile memory in computer-readable medium, random access memory (RAM) and/or The forms such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer-readable medium Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method Or technology realizes information Store.Information can be computer-readable instruction, data structure, the module of program or other data. The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM (SRAM), moved State random access memory (DRAM), other kinds of random access memory (RAM), read-only storage (ROM), electric erasable Programmable read only memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc read-only storage (CD-ROM), Digital versatile disc (DVD) or other optical storages, magnetic cassette tape, the storage of tape magnetic rigid disk or other magnetic storage apparatus Or any other non-transmission medium, the information that can be accessed by a computing device available for storage.Define, calculate according to herein Machine computer-readable recording medium does not include temporary computer readable media (transitory media), such as data-signal and carrier wave of modulation.
It should also be noted that, term " comprising ", "comprising" or its any other variant are intended to nonexcludability Comprising so that process, method, commodity or equipment including a series of key elements are not only including those key elements, but also wrap Include other key elements being not expressly set out, or also include for this process, method, commodity or equipment intrinsic want Element.In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that wanted including described Also there is other identical element in process, method, commodity or the equipment of element.
It will be understood by those skilled in the art that embodiments herein can be provided as method, system or computer program product. Therefore, the application can be using the embodiment in terms of complete hardware embodiment, complete software embodiment or combination software and hardware Form.Deposited moreover, the application can use to can use in one or more computers for wherein including computer usable program code The shape for the computer program product that storage media is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.) Formula.
Embodiments herein is the foregoing is only, the application is not limited to.For those skilled in the art For, the application can have various modifications and variations.It is all any modifications made within spirit herein and principle, equivalent Replace, improve etc., it should be included within the scope of claims hereof.

Claims (10)

1. a kind of measuring method of pulse distribution, it is characterised in that including:
Speckle interference light path is built, wherein, the artery portion of the beat pulse comprising tested organism in the speckle interference light path Position;
According to the speckle interference light path, the pulse distributed intelligence of the artery site is determined, wherein, the pulse distributed intelligence Including:Pulse space distribution information and/or pulse Annual distribution information.
2. the method as described in claim 1, it is characterised in that according to the speckle interference light path, determine the artery site Pulse space distribution information, including:
The first speckle pattern for the artery site is gathered at the first moment, and collection is directed to the artery at the second moment Second speckle pattern at position;
The corresponding first phase amount of first speckle pattern is determined, and determines the corresponding second phase of second speckle pattern Amount;
According to the first phase amount and the second phase amount, determine that the artery site is empty in the pulse at second moment Between distributed intelligence.
3. the method as described in claim 1, it is characterised in that according to the speckle interference light path, determine the artery site Pulse Annual distribution information, including:
Several speckle patterns of collection for the artery site in preset duration;
It is determined that the corresponding multiple phase masses of described several speckle patterns;
According to the multiple phase mass, pulse Annual distribution information of the artery site in the preset duration is determined.
4. the method as described in claim 1, it is characterised in that the speckle interference light path is spatial carrier digital speckle interference Light path, wherein, the spatial carrier digital speckle interference light path includes laser, the light splitting member for meeting setting coherence length value Part, beam-expanding element, diaphragm, image-forming component and imaging sensor, build the spatial carrier digital speckle interference light path, including:
The laser launches laser;
The laser via the beam splitter reflect obtain and meet set angle reference light be incident to described image pass Sensor;
The laser is transmitted through the beam-expanding element via the beam splitter, and the beam-expanding element is by the laser beam expanding outgoing After illuminate the artery site, reflecting obtained the first object light through the artery site is incident to the diaphragm, and via described Diaphragm is transmitted through the image-forming component, wherein, the diaphragm is used to cause that the spatial frequency of first object light to meet setting model Enclose;
The second object light that the image-forming component imaging is obtained is incident to described image sensor;
Second object light is interfered with the reference light on described image sensor surface.
5. the method as described in claim any one of 1-4, it is characterised in that methods described also includes:
According to the pulse space distribution information and/or the pulse Annual distribution information, the strong of the tested organism is analyzed Health situation.
6. a kind of measuring apparatus of pulse distribution, it is characterised in that the equipment includes:Speckle interference optical path unit and determination are single Member, wherein:
Speckle interference optical path unit, the artery site of the beat pulse for illuminating tested organism;
Determining unit, the pulse distributed intelligence for determining the artery site, wherein, the pulse distributed intelligence includes:Arteries and veins Fight space distribution information and/or pulse Annual distribution information.
7. equipment as claimed in claim 6, it is characterised in that the determining unit determines the pulse space of the artery site Distributed intelligence, including:
The first speckle pattern for the artery site is gathered at the first moment, and collection is directed to the artery at the second moment Second speckle pattern at position;
The corresponding first phase amount of first speckle pattern is determined, and determines the corresponding second phase of second speckle pattern Amount;
According to the first phase amount and the second phase amount, determine that the artery site is empty in the pulse at second moment Between distributed intelligence.
8. equipment as claimed in claim 6, it is characterised in that the determining unit determines the pulse space of the artery site Distributed intelligence, including:
Several speckle patterns of collection for the artery site in preset duration;
It is determined that the corresponding multiple phase masses of described several speckle patterns;
According to the multiple phase mass, Annual distribution information of the artery site in the preset duration is determined.
9. equipment as claimed in claim 6, it is characterised in that the speckle interference optical path unit is spatial carrier digital speckle Optical interference circuit unit, wherein, the spatial carrier digital speckle interference optical path unit swashs comprising satisfaction setting coherence length value Light device, beam splitter, beam-expanding element, diaphragm, image-forming component and imaging sensor, build the spatial carrier digital speckle interference Light path, including:
The laser launches laser;
The laser via the beam splitter reflect obtain and meet set angle reference light be incident to described image pass Sensor;
The laser is transmitted through the beam-expanding element via the beam splitter, and the beam-expanding element is by the laser beam expanding outgoing After illuminate the artery site, reflecting obtained the first object light through the artery site is incident to the diaphragm, and via described Diaphragm is transmitted through the image-forming component, wherein, the diaphragm is used to cause that the spatial frequency of first object light to meet setting model Enclose;
The second object light that the image-forming component imaging is obtained is incident to described image sensor;
Second object light is interfered with the reference light on described image sensor surface.
10. the equipment as described in claim any one of 6-9, it is characterised in that the equipment also includes:Analytic unit, wherein:
Analytic unit, for according to the pulse space distribution information and/or the pulse Annual distribution information, analyzing the quilt Survey the health status of organism.
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 true CN106955092A (en) 2017-07-18
CN106955092B 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)

Cited By (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

Citations (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
CN101116087A (en) * 2005-02-08 2008-01-30 皇家飞利浦电子股份有限公司 Optical speckle pattern investigation
CN101784227A (en) * 2007-07-06 2010-07-21 工业研究有限公司 laser speckle imaging systems and methods
JP2010210399A (en) * 2009-03-10 2010-09-24 Saitama Univ Minute displacement measuring method and apparatus
CN102144917A (en) * 2006-05-24 2011-08-10 泰瑞连激光技术有限公司 Optical vital sign detection method and measurement device
CN102221342A (en) * 2011-04-02 2011-10-19 北京交通大学 Method for measuring object deformation by time-domain multi-wavelength heterodyne speckle interference
WO2016029166A1 (en) * 2014-08-22 2016-02-25 Thomas Adrian Furness, Iii Automated diagnosis based at least in part on pulse waveforms

Patent Citations (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
CN101116087A (en) * 2005-02-08 2008-01-30 皇家飞利浦电子股份有限公司 Optical speckle pattern investigation
CN102144917A (en) * 2006-05-24 2011-08-10 泰瑞连激光技术有限公司 Optical vital sign detection method and measurement device
CN101784227A (en) * 2007-07-06 2010-07-21 工业研究有限公司 laser speckle imaging systems and methods
JP2010210399A (en) * 2009-03-10 2010-09-24 Saitama Univ Minute displacement measuring method and apparatus
CN102221342A (en) * 2011-04-02 2011-10-19 北京交通大学 Method for measuring object deformation by time-domain multi-wavelength heterodyne speckle interference
WO2016029166A1 (en) * 2014-08-22 2016-02-25 Thomas Adrian Furness, Iii Automated diagnosis based at least in part on pulse waveforms

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘凯等: "空间载波数字散斑干涉光路关键参数研究", 《自动化仪表》 *

Cited By (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

Also Published As

Publication number Publication date
CN106955092B (en) 2021-02-23

Similar Documents

Publication Publication Date Title
Xu et al. Pattern recognition based on time-frequency analysis and convolutional neural networks for vibrational events in φ-OTDR
Arizaga et al. Speckle time evolution characterization by the co-occurrence matrix analysis
US7123363B2 (en) Speckle pattern analysis method and system
CN109310337A (en) Skin diagnosis device, skin condition output method, program and recording medium
CN104634699B (en) Laguerre-gaussian beam-based speckle contrast imaging measurement device and laguerre-gaussian beam-based speckle contrast imaging measurement method
Smith et al. CoLux: Multi-object 3d micro-motion analysis using speckle imaging
CN110088711A (en) Magnetic disturbance detection and correction
CN110461215A (en) Health mark is determined using portable device
CN106955092A (en) A kind of measuring method and equipment of pulse distribution
CN109645954A (en) The elasticity measurement system and method for multiple beam optical coherence based on microlens array
Moukadem et al. Time-Frequency domain for segmentation and classification of Non-stationary signals: The Stockwell transform applied on Bio-signals and electric signals
Zhou et al. Vibration measurement with video processing based on alternating optimization of frequency and phase shifts
JP2017524138A (en) Instantaneous time-domain optical coherence tomography
Harvey et al. Visualization of acoustic waves in air and subsequent audio recovery with a high-speed schlieren imaging system: Experimental and computational development of a schlieren microphone
Morton et al. Transverse motions in sunspot super-penumbral fibrils
US9442459B2 (en) Making holographic data of complex waveforms
MacRitchie et al. Efficient tracking of pianists’ finger movements
US11625003B2 (en) Coherent optical holographic imaging
JP2017146696A (en) Image processing device, image processing method, and image processing system
Chen et al. The glottaltopogram: A method of analyzing high-speed images of the vocal folds
Silva et al. Spatio-temporal decomposition of 2D travelling waves from video measurements
Moro et al. Defining parametric dependencies for the correct interpretation of speckle dynamics in photon Doppler velocimetry
Cajas et al. Modeling of motion artifacts on PPG signals for heart-monitoring using wearable devices
CN109674441A (en) Multiple beam elasticity measurement system and method based on photoswitch and microlens array
Rosenberg et al. Hybrid method combining orthogonal projection Fourier transform profilometry and laser speckle imaging for 3D visualization of flow profile

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