CN110477889A - The method and apparatus for identifying radial artery wave shape reflection point - Google Patents
The method and apparatus for identifying radial artery wave shape reflection point Download PDFInfo
- Publication number
- CN110477889A CN110477889A CN201811168960.3A CN201811168960A CN110477889A CN 110477889 A CN110477889 A CN 110477889A CN 201811168960 A CN201811168960 A CN 201811168960A CN 110477889 A CN110477889 A CN 110477889A
- Authority
- CN
- China
- Prior art keywords
- waveform
- point
- differential signal
- maximum
- starting point
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Abstract
This application discloses a kind of method and apparatus for identifying radial artery wave shape reflection point.This method comprises: being filtered to radial pulse waveform；The differential signal for generating the radial pulse waveform, determines the maximum of the differential signal；The starting point of the radial pulse waveform is determined based on the maximum of the differential signal；Based on the starting point, the period of the radial pulse waveform is determined；It baseline drift is carried out to the radial pulse waveform handles to described based on the starting point and demarcates waveform；With the reflection point for identifying calibrated waveform.It using this method, can be analyzed based entirely on signal waveform, and not have to can be realized using additional hardware, compared with the complex cumbersome analytic process of the prior art, reduced the complexity and difficulty of analysis, also improve precision of analysis.
Description
Technical field
This application involves waveforms detection field, more particularly to a kind of method for identifying radial artery wave shape reflection point and
Device.
Background technique
Pulse is the important biomolecule dynamic signal for characterizing human heart function and blood circulation system, combines heart pump
The various physiological and pathological informations that blood activity and pressure wave carry along systemic arterial tree transmission process.Existing clinical trial results table
It is bright, it is closely connected between the wave character and cardiovascular disease of pulse wave.Pulse wave letter is detected by pulse wave sensor
Number, pulse waveform feature and blood pressure information are combined and analyze every cardiovascular function parameter.Since this method has behaviour
Make the features such as easy, non-invasive, detection parameters are comprehensive, is widely used in practical applications.
Currently, clinically frequently with the main superficial artery of human body, for example, the pulse wave of moon bright artery, radial artery, arteria carotis etc.
Wave character, wherein being more commonly the radial artery measurement of wrist.This method is in clinical setting using relatively simple, survey
Cost is relatively low for amount, and no sequel can be long lasting for detection, and obtained measurement result is accurate.In measured pulse waveform
In can by a series of signal processing method extract individual features parameter, thus obtain assessment cardiovascular function index of correlation
Parameter.Therefore, it is necessary to a kind of practicable sides that waveform feature parameter can be accurately extracted from measured pulse wave
Method.
The measurement blood pressure method based on radial artery pulse wave conduction time that prior art discloses a kind of, this method pass through light
Electric transducer acquires human body radial artery pulse wave signal, is filtered amplification to signal using microprocessor, AD conversion is to reduce
Noise, then frequency-domain transform, small echo processing, second differnce are carried out to determine the characteristic point of human body radial artery pulse wave waveform, it calculates
The exact numerical of human body radial artery pulse wave conduction time establishes the regression equation of measurement blood pressure, finally obtains pressure value.The party
Method used processing method when obtaining wave character is complex cumbersome, and consumed hardware resource is more, therefore needs one kind
The method that can be easier to characteristic point parameter extraction, in order to subsequent analysis treatment process.
Summary of the invention
Aiming to overcome that the above problem or at least being partially solved or extenuate for the application solves the above problems.
According to the one aspect of the application, a kind of method for identifying radial artery wave shape reflection point is provided, comprising:
Filter step: radial pulse waveform is filtered；
Differential signal generation step: the differential signal of the radial pulse waveform is generated, determines the differential signal
Maximum；
Starting point search step: rising for the radial pulse waveform is determined using the maximum and the differential signal
Initial point；
Period determines step: the period of the radial pulse waveform is determined based on the starting point；
Waveform demarcating steps: the radial pulse waveform is carried out at baseline drift to described based on the starting point
It manages and demarcates waveform；With
Feature point recognition step: the reflection point of calibrated waveform is identified.
It using this method, can be analyzed based entirely on signal waveform, and not have to can be realized using additional hardware,
Compared with the complex cumbersome analytic process of the prior art, the complexity and difficulty of analysis are reduced, also improves analysis knot
The accuracy of fruit.Convolutional filtering, first-order difference and the operation of convolution scheduling algorithm are utilized in this method, reduces to a certain extent
The complex process of data processing, the hardware computing resource occupied needed for reducing, so that the characteristic point of radial pulse waveform
Become more easily to extract, is conducive to subsequent analysis treatment process.
Optionally, the smothing filtering step includes: to carry out the radial pulse waveform and preset two-dimensional array
Convolutional filtering processing.
Optionally, in the differential signal generation step, the calculation formula for generating differential signal are as follows:
VtCpDifSig [i]=vtSmoothedData [i+1]-vtSmoothedData [i-1]+2*
(vtSmoothedData[i+2]-vtSmoothedData[i-2])
Wherein, vtCpDifSig indicates differential signal；VtSmoothedData indicates filtered radial pulse waveform；
I is outer circulation variable, indicates i-th bit signal value.
This method calculates differential signal by the way of mathematical computations, and can not have to the hardware using difference channel
It is convenient and efficient to realize operation, it is simple to can be achieved.
Optionally, in the differential signal generation step, the determination method of the waveform maximum includes:
Sectioning search step: the differential signal waveform is divided into eight sections, searches for the maximum point of each section of waveform respectively；
Mean value calculation step: the average value AvgSegmentMax of the maximum of eight sections of waveforms is calculated；With
Maximum determines step: searching for the maximum point in entire differential signal waveform, is greater than in the amplitude of maximum point
In the case where 0.6*AvgSegmentMax, the value of the maximum point is determined as to the maximum point of the differential signal waveform.
Using this method, the search of maximum is carried out by segmentation, and operation speed can be improved by way of parallel processing
Degree, the average value of obtained maximum are used as the maximum of the entire waveform of threshold decision, improve the accuracy rate of judgement.
Optionally, the starting point search step includes: to determine region of search based on the maximum point of the differential signal,
Positive zero crossing is searched in the region of search, there are positive zero crossing, using zero crossing as waveform starting point；
Otherwise, using minimum point as waveform starting point.
Optionally, it is determined in step in the period, the range difference by calculating two starting points determines the waveform week
Phase, the calculation formula of the wave period are as follows:
M_vtPeriod [i]=m_vtFootPos [i]-m_vtFootPos [i-1]
Wherein, m_vtPeriod is wave period；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit
Signal value.
The application due to using the technology described above, compared with prior art, reduces to a certain extent for pulse
The treatment process of Wave data reduces the consumption of computing resource, it is easier to the extraction of wave character point, so that subsequent utilization
The analysis treatment process of pulse wave is more convenient.
Optionally, in the waveform demarcating steps, the slope and waveform offset of point value are originated according to the waveform
It is handled to carry out baseline drift, the calculation formula of slope are as follows:
The waveform offset calculation formula are as follows:
Δ vtSlope=vtSlope [i] * (j-m_vtFootPos [i])
Wherein, vtSlope indicates slope；Δ vtSlope indicates waveform offset；After vtSmoothedData indicates filtering
Radial pulse waveform；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit signal value；J is interior circulation
Variable.
Optionally, the Feature point recognition step includes:
Convolution step: the scale based on setting obtains the calibrated waveform Fourth-Derivative by convolution；
Reflection point determines step: first that Fourth-Derivative is searched in from starting point to 0.67 wave period region is reversed
Zero crossing, the zero crossing are reflection point.
According to further aspect of the application, a kind of device for identifying radial artery wave shape reflection point, packet are additionally provided
It includes:
Filter module is disposed for being filtered radial pulse waveform；
Differential signal generation module is disposed for generating the differential signal of the radial pulse waveform, determines institute
State the maximum of differential signal；
Starting point search module is disposed for determining the radial artery using the maximum and the differential signal
The starting point of pulse wave；
Period determination module is disposed for determining the period of the radial pulse waveform based on the starting point；
Waveform demarcating module is disposed for carrying out the radial pulse waveform to described based on the starting point
It goes baseline drift to handle and demarcates waveform；
Feature point recognition module is disposed for identifying the reflection point of calibrated waveform.
It using the device, can be analyzed based entirely on signal waveform, and not have to can be realized using additional hardware,
Compared with the complex cumbersome analytic process of the prior art, the complexity and difficulty of analysis are reduced, also improves analysis knot
The accuracy of fruit.Convolutional filtering, first-order difference and the operation of convolution scheduling algorithm are utilized in this method, reduces to a certain extent
The complex process of data processing, the hardware computing resource occupied needed for reducing, so that the characteristic point of radial pulse waveform
Become more easily to extract, is conducive to subsequent analysis treatment process.
Optionally, include: in the differential signal generation module
Sectioning search module is disposed for the radial pulse waveform being divided into eight sections, searches for each section respectively
The maximum point of waveform；
Mean value calculation module is disposed for calculating the average value AvgSegmentMax of the maximum of eight sections of waveforms；
With
Maximum determining module is disposed for the maximum point searched in entire radial pulse waveform, very big
In the case that the amplitude of value point is greater than 0.6*AvgSegmentMax, the value of the maximum point is determined as the radial pulse
The maximum point of waveform.
According to the accompanying drawings to the detailed description of the specific embodiment of the application, those skilled in the art will be more
Above-mentioned and other purposes, the advantages and features of the application are illustrated.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present application, constitutes part of this application, this Shen
Illustrative embodiments and their description please are not constituted an undue limitation on the present application for explaining the application.Hereinafter by reference
Some specific embodiments of the application are described in detail by way of example and not limitation in attached drawing.Identical attached drawing mark in attached drawing
Note denotes same or similar part or part.It will be understood by those skilled in the art that these attached drawings be not necessarily by than
What example was drawn.In the accompanying drawings:
Fig. 1 is the computer dress for running the identification radial artery wave shape reflection point methods according to the application one embodiment
Set hardware structural diagram；
Fig. 2 is the schematic flow according to the identification radial artery wave shape reflection point methods of one embodiment of the application
Figure；
Fig. 3 is that the overall procedure for reflecting point methods according to the identification radial artery wave shape of one embodiment of the application shows
It is intended to；
Fig. 4 is the schematic flow chart according to the search reflection point of the method for one embodiment of the application；
Fig. 5 is the position view according to the reflection point of one embodiment of the application on radial pulse waveform；
Fig. 6 is the schematic frame according to the identification radial artery wave shape reflection point device of one embodiment of the application
Figure；
Fig. 7 is the block diagram of one embodiment of the calculating equipment of the application；
Fig. 8 is the block diagram of one embodiment of the computer readable storage medium of the application.
Specific embodiment
In order to make those skilled in the art more fully understand application scheme, below in conjunction in the embodiment of the present application
Attached drawing, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described embodiment is only
The embodiment of the application a part, instead of all the embodiments.Based on the embodiment in the application, ordinary skill people
Member's every other embodiment obtained without making creative work, all should belong to the model of the application protection
It encloses.
It should be noted that the description and claims of this application and term " first " in above-mentioned attached drawing, "
Two " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way
Data be interchangeable under appropriate circumstances, so as to embodiments herein described herein can in addition to illustrating herein or
Sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that cover
Cover it is non-exclusive include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to
Step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, product
Or other step or units that equipment is intrinsic.
Embodiment 1
According to the embodiment of the present application, a kind of embodiment of the method for identifying radial artery wave shape characteristic point is additionally provided, is needed
It is noted that step shown in the flowchart of the accompanying drawings can be in the computer system of such as a group of computer-executable instructions
Middle execution, although also, logical order is shown in flow charts, and it in some cases, can be to be different from herein
Sequence executes shown or described step.
Embodiment of the method provided by the embodiment of the present application one can be in mobile terminal, terminal or similar fortune
It calculates and is executed in device.Fig. 1 shows a kind of hardware configuration frame of computer installation (or mobile device) for the application method
Figure.As shown in Figure 1, computer installation 10 (or mobile device 10) may include one or more processors (in figure using 102a,
102b ... ..., 102n are shown, and processor can include but is not limited to Micro-processor MCV or programmable logic device FPGA etc.
Processing unit), memory 104 for storing data and the transmission module for communication function.In addition to this, may be used also
To include: that display, input/output interface (I/O interface), the port universal serial bus (USB) (can be used as I/O interface
A port in port included), network interface, power supply and/or camera.It will appreciated by the skilled person that Fig. 1
Shown in structure be only illustrate, the structure of above-mentioned electronic device is not caused to limit.For example, computer installation 10 can also wrap
Include than shown in Fig. 1 more perhaps less component or with the configuration different from shown in Fig. 1.
It is to be noted that said one or multiple processors and/or other data processing circuits usually may be used herein
To be referred to as " data processing circuit ".The data processing circuit all or part of can be presented as software, hardware, firmware or its
His any combination.In addition, data processing circuit for single independent processing module or all or part of can be integrated to computer
In any one in other elements in device 10 (or mobile device).As involved in the embodiment of the present application, the number
(such as the selection for the variable resistance end path connecting with interface) is controlled as a kind of processor according to processing circuit.
Memory 104 can be used for storing the software program and module of application software, such as the method in the embodiment of the present application
Corresponding program instruction/data storage device, the software program and mould that processor is stored in memory 104 by operation
Block realizes the method for above-mentioned application program thereby executing various function application and data processing.Memory 104 can wrap
Include high speed random access memory, may also include nonvolatile memory, as one or more magnetic storage device, flash memory or
Other non-volatile solid state memories.In some instances, memory 104 can further comprise remotely located relative to processor
Memory, these remote memories can pass through network connection to computer installation 10.The example of above-mentioned network includes but not
It is limited to internet, intranet, local area network, mobile radio communication and combinations thereof.
Transmitting device is used to that data to be received or sent via a network.Above-mentioned network specific example may include calculating
The wireless network that the communication providers of machine device 10 provide.In an example, transmitting device includes a network adapter
(Network Interface Controller, NIC), can be connected by base station with other network equipments so as to interconnection
Net is communicated.In an example, transmitting device can be radio frequency (Radio Frequency, RF) module, be used to pass through
Wireless mode is communicated with internet.
Display can such as touch-screen type liquid crystal display (LCD), the liquid crystal display aloow user with
The user interface of computer installation 10 (or mobile device) interacts.
Under above-mentioned running environment, this application provides a kind of identification radial artery wave shape characteristic points, especially reflect
The method of point (E point).Fig. 2 is the signal according to the identification radial artery wave shape reflection point methods of one embodiment of the application
Property flow chart.This method may comprise steps of:
S100 filter step: radial pulse waveform is filtered；
S200 differential signal generation step: determining the maximum of the radial pulse waveform, and it is dynamic to generate the oar
The differential signal of arteries and veins pulse wave；
S300 starting point search step: rising for the radial pulse waveform is determined based on the maximum of the differential signal
Initial point；
The S400 period determines step: the period of the radial pulse waveform is determined based on the starting point；
S500 waveform demarcating steps: it baseline is carried out to the radial pulse waveform floats to described based on the starting point
It moves and handles and demarcate waveform；With
S600 Feature point recognition step: the reflection point of calibrated waveform is identified.
It using this method, can be analyzed based entirely on signal waveform, and not have to can be realized using additional hardware,
Compared with the complex cumbersome analytic process of the prior art, the complexity and difficulty of analysis are reduced, also improves analysis knot
The accuracy of fruit.Convolutional filtering, first-order difference and the operation of convolution scheduling algorithm are utilized in this method, reduces to a certain extent
The complex process of data processing, the hardware computing resource occupied needed for reducing, so that the characteristic point of radial pulse waveform
Become more easily to extract, is conducive to subsequent analysis treatment process.
Optionally, the smothing filtering step includes: to carry out the radial pulse waveform and preset two-dimensional array
Convolutional filtering processing.Fig. 3 is the totality according to the identification radial artery wave shape reflection point methods of one embodiment of the application
Flow diagram.For example, the array can be the two-dimensional array of a 11*14: QG [11] [14].By pulse wave and filter
Coefficient carries out convolutional filtering processing, to obtain filtered pulse wave.
Wherein, each element can be the specific value such as following table one in the specific two-dimensional array QG [11] [14]:
Table one
Optionally, in the differential signal generation step, the calculation formula for generating differential signal are as follows:
VtCpDifSig [i]=vtSmoothedData [i+1]-vtSmoothedData [i-1]+2*
(vtSmoothedData[i+2]-vtSmoothedData[i-2])
Wherein, vtCpDifSig indicates differential signal；VtSmoothedData indicates filtered radial pulse waveform；
I is outer circulation variable, indicates i-th bit signal value.
This method calculates differential signal by the way of mathematical computations, and can not have to the hardware using difference channel
It is convenient and efficient to realize operation, it is simple to can be achieved.
Optionally, in the differential signal generation step, the determination method of the waveform maximum includes:
Sectioning search step: the differential signal waveform is divided into eight sections, searches for the maximum point of each section of waveform respectively；
Mean value calculation step: the average value AvgSegmentMax of the maximum of eight sections of waveforms is calculated；
Maximum determines step: searching for the maximum point in entire differential signal waveform, is greater than in the amplitude of maximum point
In the case where 0.6*AvgSegmentMax, the value of the maximum point is determined as to the maximum point of the differential signal waveform.
Using this method, the search of maximum is carried out by segmentation, and operation speed can be improved by way of parallel processing
Degree, the average value of obtained maximum are used as the maximum of the entire waveform of threshold decision, improve the accuracy rate of judgement.
Optionally, the starting point search step includes: to determine region of search based on the maximum point of the differential signal,
Positive zero crossing is searched in the region of search, there are positive zero crossing, using zero crossing as waveform starting point；
Otherwise, using minimum point as waveform starting point.
The process for determining starting point by differential signal maximum point is: (general to occur by differential signal maximum point
After starting point), 200ms is searched for forward and determines region of search, in the region of search, when there is differential signal positive going zeror crossing
When point, as waveform starting point；When positive zero crossing is not present, using the region of search minimum point as starting point.
Optionally, it is determined in step in the period, the range difference by calculating two starting points determines the waveform week
Phase, the calculation formula of the wave period are as follows:
M_vtPeriod [i]=m_vtFootPos [i]-m_vtFootPos [i-1]
Wherein, m_vtPeriod is wave period；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit
Signal value.
The application due to using the technology described above, compared with prior art, reduces to a certain extent for pulse
The treatment process of Wave data reduces the consumption of computing resource, it is easier to the extraction of wave character point, so that subsequent utilization
The analysis treatment process of pulse wave is more convenient.
Referring to figs. 2 and 3, optionally, after the period determines step, this method further includes that the peak point determines step,
The step includes the acquisition waveforms maximum of points in 2/3 section of the starting point to the period, by the waveform maximum value
Point is used as the peak point.
This method uses the acquisition waveforms maximum of points in 2/3 section of the starting point to the period, can subtract
The quantity of few data processing, improves calculating speed.
Optionally, in the waveform demarcating steps, the slope and waveform offset of point value are originated according to the waveform
It is handled to carry out baseline drift, the calculation formula of slope are as follows:
The waveform offset calculation formula are as follows:
Δ vtSlope=vtSlope [i] * (j-m_vtFootPos [i])
Wherein, vtSlope indicates slope；Δ vtSlope indicates waveform offset；After vtSmoothedData indicates filtering
Radial pulse waveform；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit signal value；J is interior circulation
Variable.
Optionally, in one embodiment, the waveform calibration is dynamic to oar according to the systolic pressure and diastolic pressure of arteria brachialis
The amplitude of pulse wave shape is demarcated, and blood pressure unit is made it have.
Referring to figs. 2 and 3, optionally, the Feature point recognition step includes:
Convolution step: the scale based on setting obtains the calibrated waveform Fourth-Derivative by convolution
vtScale2ConvSigNtoP；For example, scale scale1=50, scale2=100 can be set；
Reflection point determines step: first that Fourth-Derivative is searched in from starting point to 0.67 wave period region is reversed
Zero crossing, the zero crossing are reflection point.
Fig. 4 is the schematic flow chart according to the search reflection point of the method for one embodiment of the application.Pass through setting
Good scale is obtained the Fourth-Derivative of waveform using convolution operation, keeps its obtained waveform burr less, and without at subsequent filtering
Reason.The reversed zero crossing for then searching for the Fourth-Derivative in 2/3 region before entire wave period, to obtain the anti-of waveform
Exit point E point.
Wherein, double scale convolution obtain Fourth-Derivative process are as follows:
(1) Wave data after a convolution is extracted；(2) Wave data after secondary convolution is extracted.
Wherein, a convolution is extracted, by the way that Two-way Cycle is arranged, outer circulation variable is i, and interior cyclic variable is j, when i-j's
Absolute value≤0.5*Scale1, a convolution value=(j+ starting point) range value * ((absolute value of -8.0/3.0* (i-j)/
Scale1)+(5.0/6.0))；Work as Scale1 > (i-j) absolute value > 0.5*Scale1, a convolution value=(j+ starting point)
Range value * ((absolute value/Scale1 of 7.0/3.0* (i-j))+(13.0/12.0))；When Scale1≤(i-j) absolute value,
Convolution value=(j+ starting point) range value * ((absolute value/Scale1 of -1.0/6.0* (i-j))+(1.0/4.0)).
Wherein, secondary convolution is extracted, by the way that Two-way Cycle is arranged, outer circulation variable is i, and interior cyclic variable is j, when i-j's
Absolute value≤0.5*Scale2, the range value * of secondary convolution value=mono- time convolution (j+ starting point) be ((- 8.0/3.0*'s (i-j)
Absolute value/Scale2)+(5.0/6.0))；Work as Scale2 > (i-j) absolute value > Scale1/2, a secondary convolution value=secondary volume
The range value * ((absolute value/Scale2 of 7.0/3.0* (i-j))+(13.0/12.0)) of product (j+ starting point)；When Scale2≤
(i-j) absolute value, ((- 1.0/6.0*'s (i-j) is absolute by the range value * of secondary convolution value=mono- time convolution (j+ starting point)
Value/Scale2)+(1.0/4.0))；The secondary convolution value is the Fourth-Derivative vtScale2ConvSigNtoP of waveform.
Fig. 5 is the position view according to the reflection point of one embodiment of the application on radial pulse waveform.It should
Including in figure is complete signal period radial pulse waveform, is mainly made of ascending branch and descending branch.Ascending branch and descending branch composition
Main wave 1 has an incisura to be known as dicrotic notch 3 on descending branch, often occurs winning prewave 2, also known as tidal wave again between main wave 1 and dicrotic notch 3.
The immediately appearance of dicrotic notch 3 is dicrotic wave 4, also known as dicrotic wave.Above-described wave and gorge be constitute pulse wave it is main at
Point.
Specifically, the 1 trough B point of main wave is main artery open point, the pressure turnover that aorta petal starts to open is represented
Point indicates the beginning of heart phase of maximum ejection；C point is main 1 wave crest of wave, is aortic pressure highest point, it is dynamic to represent the systole phase
Pulse pressure highest point；E point is 2 wave crest of dicrotic pulse prewave, i.e. left ventricular ejection halt, the arterial pressure formed by blood back flow
It increases.When human body artery tube wall smooth muscle is less, elastomer is more, the main wave C high of pulse wave is formed and sharp, tidal wave E is not
Obviously.Conversely, ductus arteriosus wall smooth muscle increase it is less with tube wall elasticity fiber, back wave velocity of wave increase, so that tidal wave E is gradually become bright
It is aobvious.It can be seen that the height fluctuations of pulse wave characteristic point can directly reflect the change of human vas resistance and vessel wall elasticity
Change, is the objective important symbol for measuring the variation of human body difference physiological and pathological.
Embodiment 2
According to the embodiment of the present application, additionally provide a kind of identification radial artery wave shape reflection point device, the device be with
The corresponding device of method described in embodiment 1.Fig. 6 is anti-according to the identification radial artery wave shape of one embodiment of the application
The schematic block diagram of exit point device.The apparatus may include:
Filter module 100 is disposed for being filtered radial pulse waveform；
Differential signal generation module 200 is disposed for generating the differential signal of the radial pulse waveform, determines
The maximum of the differential signal；
Starting point search module 300 is disposed for determining the radial artery based on the maximum of the differential signal
The starting point of pulse wave；
Period determination module 400 is disposed for determining the week of the radial pulse waveform based on the starting point
Phase；
Waveform demarcating module 500 is disposed for based on the starting point to described to the radial pulse waveform
Baseline drift is carried out to handle and demarcate waveform；With
Feature point recognition module 600 is disposed for identifying the reflection point of calibrated waveform.
It using the device, can be analyzed based entirely on signal waveform, and not have to can be realized using additional hardware,
Compared with the complex cumbersome analytic process of the prior art, the complexity and difficulty of analysis are reduced, also improves analysis knot
The accuracy of fruit.Convolutional filtering, first-order difference and the operation of convolution scheduling algorithm are utilized in the device, reduces to a certain extent
The complex process of data processing, the hardware computing resource occupied needed for reducing, so that the characteristic point of radial pulse waveform
Become more easily to extract, is conducive to subsequent analysis treatment process.
Optionally, in the differential signal generation module, the calculation formula for generating differential signal are as follows:
VtCpDifSig [i]=vtSmoothedData [i+1]-vtSmoothedData [i-1]+2*
(vtSmoothedData[i+2]-vtSmoothedData[i-2])
Wherein, vtCpDifSig indicates differential signal；VtSmoothedData indicates filtered radial pulse waveform；
I is outer circulation variable, indicates i-th bit signal value.
The device calculates differential signal by the way of mathematical computations, and can not have to the hardware using difference channel
It is convenient and efficient to realize operation, it is simple to can be achieved.
Optionally, the differential signal generation module includes:
Sectioning search module is disposed for the differential signal waveform being divided into eight sections, searches for each section of wave respectively
The maximum point of shape；
Mean value calculation module is disposed for calculating the average value AvgSegmentMax of the maximum of eight sections of waveforms；
Maximum determining module is disposed for the maximum point searched in entire differential signal waveform, in maximum
In the case that the amplitude of point is greater than 0.6*AvgSegmentMax, the value of the maximum point is determined as the radial artery pulse wave
The maximum point of shape.
Using the device, the search of maximum is carried out by segmentation, and operation speed can be improved by way of parallel processing
Degree, the average value of obtained maximum are used as the maximum of the entire waveform of threshold decision, improve the accuracy rate of judgement.
Optionally, the starting point search module is also used to: determining the field of search based on the maximum point of the differential signal
Positive zero crossing is searched in the region of search in domain, there are positive zero crossing, is originated using zero crossing as waveform
Point；Otherwise, using minimum point as waveform starting point.
Optionally, in the period determination module, the range difference by calculating two starting points determines the waveform week
Phase, the calculation formula of the wave period are as follows:
M_vtPeriod [i]=m_vtFootPos [i]-m_vtFootPos [i-1]
Wherein, m_vtPeriod is wave period；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit
Signal value.
The application due to using the technology described above, compared with prior art, reduces to a certain extent for pulse
The treatment process of Wave data reduces the consumption of computing resource, it is easier to the extraction of wave character point, so that subsequent utilization
The analysis treatment process of pulse wave is more convenient.
Optionally, after period determination module, which can also include that the peak point determining module is also used in institute
Acquisition waveforms maximum of points in starting point to 2/3 section in the period is stated, using the waveform maximum of points as the peak value
Point.
The device uses the acquisition waveforms maximum of points in 2/3 section of the starting point to the period, can subtract
The quantity of few data processing, improves calculating speed.
Optionally, in the waveform demarcating module, the slope and waveform offset of point value are originated according to the waveform
It is handled to carry out baseline drift, the calculation formula of slope are as follows:
The waveform offset calculation formula are as follows:
Δ vtSlope=vtSlope [i] * (j-m_vtFootPos [i])
Wherein, vtSlope indicates slope；Δ vtSlope indicates waveform offset；After vtSmoothedData indicates filtering
Radial pulse waveform；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit signal value；J is interior circulation
Variable.
Optionally, in one embodiment, the waveform calibration is dynamic to oar according to the systolic pressure and diastolic pressure of arteria brachialis
The amplitude of pulse wave shape is demarcated, and blood pressure unit is made it have.
Optionally, the Feature point recognition module includes:
Convolution module is disposed for the scale based on setting, obtains the calibrated waveform quadravalence by convolution
Derivative vtScale2ConvSigNtoP；For example, scale scale1=50, scale2=100 can be set；
Reflection point determining module is disposed for searching for Fourth-Derivative in from starting point to 0.67 wave period region
First reversed zero crossing, the zero crossing be reflection point.
Above-mentioned the embodiment of the present application serial number is for illustration only, does not represent the advantages or disadvantages of the embodiments.
In above-described embodiment of the application, all emphasizes particularly on different fields to the description of each embodiment, do not have in some embodiment
The part of detailed description, reference can be made to the related descriptions of other embodiments.
In several embodiments provided herein, it should be understood that disclosed technology contents can pass through others
Mode is realized.Wherein, the apparatus embodiments described above are merely exemplary, such as the division of the unit, only
A kind of logical function partition, there may be another division manner in actual implementation, for example, multiple units or components can combine or
Person is desirably integrated into another system, or some features can be ignored or not executed.Another point, shown or discussed is mutual
Between coupling, direct-coupling or communication connection can be through some interfaces, the INDIRECT COUPLING or communication link of unit or module
It connects, can be electrical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
It, can also be in addition, each functional unit in each embodiment of the application can integrate in one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product
When, it can store in a computer readable storage medium.Based on this understanding, the technical solution of the application is substantially
The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words
It embodies, which is stored in a storage medium, including some instructions are used so that a computer
Equipment (can for personal computer, server or network equipment etc.) execute each embodiment the method for the application whole or
Part steps.And storage medium above-mentioned includes: that USB flash disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), mobile hard disk, magnetic or disk etc. be various to can store program code
Medium.
Embodiment 3
The one aspect of embodiments herein provides a kind of calculating equipment, and referring to Fig. 7, which includes storage
Device 1120, processor 1110 and it is stored in the computer journey that can be run in the memory 1120 and by the processor 1110
Sequence, the computer program are stored in the space 1130 for program code in memory 1120, the computer program by
Reason device 1110 is realized when executing for executing any one according to one of above method step of the application 1131.
The one aspect of embodiments herein additionally provides a kind of computer readable storage medium.Referring to Fig. 8, the calculating
Machine readable storage medium storing program for executing includes the storage unit for program code, which is provided with for executing according to the application's
The program 1131 ' of one of above method step, the program are executed by processor.
The one aspect of the embodiment of the present application additionally provides a kind of computer program product comprising instruction, including computer
Readable code causes the calculating equipment to execute as described above when the computer-readable code is executed by calculating equipment
Method.
In the above-described embodiments, can come wholly or partly by software, hardware, firmware or any combination thereof real
It is existing.When implemented in software, it can entirely or partly realize in the form of a computer program product.The computer program
Product includes one or more computer instructions.When computer loads and executes the computer program instructions, whole or portion
Ground is divided to generate according to process or function described in the embodiment of the present application.The computer can be general purpose computer, dedicated computing
Machine, computer network obtain other programmable devices.The computer instruction can store in computer readable storage medium
In, or from a computer readable storage medium to the transmission of another computer readable storage medium, for example, the computer
Instruction can pass through wired (such as coaxial cable, optical fiber, number from a web-site, computer, server or data center
User's line (DSL)) or wireless (such as infrared, wireless, microwave etc.) mode to another web-site, computer, server or
Data center is transmitted.The computer readable storage medium can be any usable medium that computer can access or
It is comprising data storage devices such as one or more usable mediums integrated server, data centers.The usable medium can be with
It is magnetic medium, (for example, floppy disk, hard disk, tape), optical medium (for example, DVD) or semiconductor medium (such as solid state hard disk
Solid State Disk (SSD)) etc..
Professional should further appreciate that, described in conjunction with the examples disclosed in the embodiments of the present disclosure
Unit and algorithm steps, can be realized with electronic hardware, computer software, or a combination of the two, hard in order to clearly demonstrate
The interchangeability of part and software generally describes each exemplary composition and step according to function in the above description.
These functions are implemented in hardware or software actually, the specific application and design constraint depending on technical solution.
Professional technician can use different methods to achieve the described function each specific application, but this realization
It is not considered that exceeding scope of the present application.
Those of ordinary skill in the art will appreciate that implement the method for the above embodiments be can be with
By program come instruction processing unit completion, the program be can store in computer readable storage medium, and the storage is situated between
Matter is non-transitory (English: non-transitory) medium, such as random access memory, read-only memory, flash
Device, hard disk, solid state hard disk, tape (English: magnetic tape), floppy disk (English: floppy disk), CD (English:
Optical disc) and any combination thereof.
The preferable specific embodiment of the above, only the application, but the protection scope of the application is not limited thereto,
Within the technical scope of the present application, any changes or substitutions that can be easily thought of by anyone skilled in the art,
Should all it cover within the scope of protection of this application.Therefore, the protection scope of the application should be with scope of protection of the claims
Subject to.
Claims (10)
1. a kind of method for identifying radial artery wave shape reflection point, comprising:
Filter step: radial pulse waveform is filtered；
Differential signal generation step: the differential signal of the radial pulse waveform is generated, determines the very big of the differential signal
Value；
Starting point search step: the starting point of the radial pulse waveform is determined based on the maximum of the differential signal；
Period determines step: the period of the radial pulse waveform is determined based on the starting point；
Waveform demarcating steps: baseline drift processing is carried out simultaneously to the radial pulse waveform to described based on the starting point
Demarcate waveform；With
Feature point recognition step: the reflection point of calibrated waveform is identified.
2. the method according to claim 1, wherein the smothing filtering step includes: by the radial artery arteries and veins
Fight waveform and preset two-dimensional array carries out convolutional filtering processing.
3. the method according to claim 1, wherein the generation is poor in the differential signal generation step
The calculation formula of sub-signal are as follows:
VtCpDifSig [i]=vtSmoothedData [i+1]-vtSmoothedData [i-1]+2* (vtSmoothedData
[i+2]-vtSmoothedData[i-2])
Wherein, vtCpDifSig indicates differential signal；VtSmoothedData indicates filtered radial pulse waveform；I is
Outer circulation variable indicates i-th bit signal value.
4. the method according to claim 1, wherein in the differential signal generation step, the waveform pole
The determination method being worth greatly includes:
Sectioning search step: the differential signal waveform is divided into eight sections, searches for the maximum point of each section of waveform respectively；
Mean value calculation step: the average value AvgSegmentMax of the maximum of eight sections of waveforms is calculated；With
Maximum determines step: searching for the maximum point in entire differential signal waveform, is greater than 0.6* in the amplitude of maximum point
In the case where AvgSegmentMax, the value of the maximum point is determined as to the maximum point of the differential signal.
5. the method according to claim 1, which is characterized in that the starting point search step includes: based on the differential signal
Maximum point determine region of search, search for positive zero crossing in the region of search, there are positive zero crossing,
Using zero crossing as waveform starting point；Otherwise, using minimum point as waveform starting point.
6. being risen the method according to claim 1, wherein being determined in step in the period by calculating two
The range difference of initial point determines the wave period, the calculation formula of the wave period are as follows:
M_vtPeriod [i]=m_vtFootPos [i]-m_vtFootPos [i-1]
Wherein, m_vtPeriod is wave period；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit signal
Value.
7. the method according to claim 1, wherein being risen in the waveform demarcating steps according to the waveform
The slope and waveform offset of initial point value are handled to carry out baseline drift, the calculation formula of slope are as follows:
The waveform offset calculation formula are as follows:
Δ vtSlope=vtSlope [i] * (j-m_vtFootPos [i])
Wherein, vtSlope indicates slope；Δ vtSlope indicates waveform offset；VtSmoothedData indicates filtered oar
Arterial pulse waveform；M_vtFootPos is starting point；I is outer circulation variable, indicates i-th bit signal value；J is interior cyclic variable.
8. method according to any one of claim 1 to 7, which is characterized in that the Feature point recognition step includes:
Convolution step: the scale based on setting obtains the calibrated waveform Fourth-Derivative by convolution；With
Reflection point determines step: first reversed zero passage of Fourth-Derivative is searched in from starting point to 0.67 wave period region
Point, the zero crossing are reflection point.
9. a kind of device for identifying radial artery wave shape reflection point, comprising:
Filter module is disposed for being filtered radial pulse waveform；
Differential signal generation module is disposed for generating the differential signal of the radial pulse waveform, determines the difference
The maximum of sub-signal；
Starting point search module is disposed for determining the radial pulse waveform based on the maximum of the differential signal
Starting point；
Period determination module is disposed for determining the period of the radial pulse waveform based on the starting point；
Waveform demarcating module is disposed for carrying out base to the radial pulse waveform to described based on the starting point
Line drift handles and demarcates waveform；With
Feature point recognition module is disposed for identifying the reflection point of calibrated waveform.
10. the apparatus according to claim 1, which is characterized in that include: in the differential signal generation module
Sectioning search module is disposed for the differential signal waveform being divided into eight sections, searches for each section of waveform respectively
Maximum point；
Mean value calculation module is disposed for calculating the average value AvgSegmentMax of the maximum of eight sections of waveforms；With
Maximum determining module is disposed for the maximum point searched in entire differential signal waveform, in maximum point
In the case that amplitude is greater than 0.6*AvgSegmentMax, the value of the maximum point is determined as to the pole of the differential signal waveform
Big value point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811168960.3A CN110477889B (en) | 2018-10-08 | 2018-10-08 | Method and device for identifying reflection point of radial artery pressure waveform |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811168960.3A CN110477889B (en) | 2018-10-08 | 2018-10-08 | Method and device for identifying reflection point of radial artery pressure waveform |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110477889A true CN110477889A (en) | 2019-11-22 |
| CN110477889B CN110477889B (en) | 2021-11-02 |
Family
ID=68545749
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811168960.3A Active CN110477889B (en) | 2018-10-08 | 2018-10-08 | Method and device for identifying reflection point of radial artery pressure waveform |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110477889B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112932475A (en) * | 2021-02-01 | 2021-06-11 | 武汉泰利美信医疗科技有限公司 | Method and device for calculating blood oxygen saturation, electronic equipment and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5027641A (en) * | 1989-02-23 | 1991-07-02 | Costello Jr Leo F | Oscillometric non-invasive blood pressure simulator |
| CN102307520A (en) * | 2009-02-03 | 2012-01-04 | 欧姆龙健康医疗事业株式会社 | Pulse wave analyzing device and pulse wave analyzing method |
| CN103169456A (en) * | 2013-03-29 | 2013-06-26 | 深圳职业技术学院 | Processing method and processing system for pulse wave signals |
| CN106691406A (en) * | 2017-01-05 | 2017-05-24 | 大连理工大学 | Detection method of vascular elasticity and blood pressure based on single probe photoplethysmography pulse wave |
| CN107898443A (en) * | 2017-11-21 | 2018-04-13 | 深圳先进技术研究院 | Dicrotic pulse wave detecting method, device and computer-readable storage medium |
-
2018
- 2018-10-08 CN CN201811168960.3A patent/CN110477889B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5027641A (en) * | 1989-02-23 | 1991-07-02 | Costello Jr Leo F | Oscillometric non-invasive blood pressure simulator |
| CN102307520A (en) * | 2009-02-03 | 2012-01-04 | 欧姆龙健康医疗事业株式会社 | Pulse wave analyzing device and pulse wave analyzing method |
| CN103169456A (en) * | 2013-03-29 | 2013-06-26 | 深圳职业技术学院 | Processing method and processing system for pulse wave signals |
| CN106691406A (en) * | 2017-01-05 | 2017-05-24 | 大连理工大学 | Detection method of vascular elasticity and blood pressure based on single probe photoplethysmography pulse wave |
| CN107898443A (en) * | 2017-11-21 | 2018-04-13 | 深圳先进技术研究院 | Dicrotic pulse wave detecting method, device and computer-readable storage medium |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112932475A (en) * | 2021-02-01 | 2021-06-11 | 武汉泰利美信医疗科技有限公司 | Method and device for calculating blood oxygen saturation, electronic equipment and storage medium |
| CN112932475B (en) * | 2021-02-01 | 2023-02-21 | 武汉泰利美信医疗科技有限公司 | Method and device for calculating blood oxygen saturation, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110477889B (en) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102429649B (en) | Continuous blood pressure measuring device | |
| CN102397064B (en) | Continuous blood pressure measuring device | |
| CN102488503B (en) | Continuous blood pressure measurer | |
| WO2021208490A1 (en) | Blood pressure measuring method and device based on deep neural network | |
| EP2992820B1 (en) | Measuring blood pressure | |
| CN107530005B (en) | Method and device for deriving mean arterial pressure of a subject | |
| CN108498089A (en) | A kind of noninvasive continuous BP measurement method based on deep neural network | |
| CN110251105A (en) | A kind of non-invasive blood pressure measuring method, device, equipment and system | |
| CN106413534A (en) | Blood-pressure continuous-measurement device, measurement model establishment method, and system | |
| CN106037694A (en) | Continuous blood pressure measuring device based on pulse waves | |
| CN109480800B (en) | Apparatus and method for estimating biological information and blood pressure monitoring device | |
| RU2011119486A (en) | SYSTEM AND DEVICE FOR NON-INVASIVE ARTERIAL PRESSURE MEASUREMENT | |
| RU2011122787A (en) | METHODS AND SYSTEMS FOR NON-INVASIVE MEASUREMENT OF GLUCOSE LEVELS | |
| CN108992054B (en) | Pulse signal peak point detection method and device | |
| CN110037668A (en) | The system that pulse signal time-space domain binding model judges age, health status and malignant arrhythmia identification | |
| CN109872820B (en) | Method, device, equipment and storage medium for measuring blood pressure without cuff | |
| CN106798552A (en) | Heart rate with status switch optimization estimates equipment | |
| WO2020125494A1 (en) | Analysis processing method, apparatus, and device for biological signal | |
| KR20200021208A (en) | Apparatus and method for estimating blood pressure | |
| CN106659404A (en) | Continuous blood pressure measurement method, apparatus and device | |
| CN108403101A (en) | A kind of model training method and device, the method and device for measuring blood pressure | |
| CN108124419B (en) | Blood pressure measuring method, intelligent heart rate earphone and system | |
| CN110200642A (en) | A kind of measurement method and terminal of cognitive load and psychological pressure | |
| CN110477889A (en) | The method and apparatus for identifying radial artery wave shape reflection point | |
| CN111248881B (en) | Method and device for predicting blood pressure |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhang Yongliang Inventor after: Ye Jun Inventor after: Li Daoqing Inventor after: Zhang Qilian Inventor before: Li Daoqing Inventor before: Zhang Qilian |
|
| CB03 | Change of inventor or designer information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210414 Address after: 230000 1707, building 1, science and Technology Industrial Park (Zhongwang sub Park), 198 Mingzhu Avenue, high tech Zone, Hefei City, Anhui Province Applicant after: Hefei Yiyang Health Technology Co.,Ltd. Applicant after: BEIJING DONGLIANG HEALTH TECHNOLOGY Co.,Ltd. Address before: 100038 room 2033, 2 / F, 21 YANGFANGDIAN Road, Haidian District, Beijing Applicant before: Heyunli Health Technology (Beijing) Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |