The method and apparatus for identifying radial artery wave shape dicrotic notch characteristic point
Technical field
This application involves waveforms detection fields, more particularly to a kind of identification radial artery wave shape dicrotic notch characteristic point
Method and apparatus.
Background technique
In today's society, with the continuous development in epoch, people's lives level is continuously improved, and more and more people are more
Focus on the health status of itself.In China, non-contagious chronic diseases have become the No.1 problem for hindering health, center
Vascular diseases are even more to rank the first, and incidence and mortality has become China city and rural crowd often suffers from the disease and lethal original
Because first of.Therefore, how to detect and prevent cardiovascular disease and become for people pay close attention to the most instantly a problem.
People usually know self health status by measuring the method for blood pressure in daily life.Common measurement side
There are two types of formulas, direct-type and indirect type, and wherein direct-type measurement is invasive continuous mode method, is usually used in hospital surgical implementation process
In.Measurement indirectly the parameters such as changes by arterial blood tube wall fluctuation, volumetric blood and measures blood pressure indirectly, generally stethoscopy or shows
Wave method, technology is more mature, is widely used in clinical setting.Usually utilize pertinent instruments measurement human body upper arm or the wrist upper arm
Artery or radial artery pulse wave, and therefrom analytical calculation to obtain pressure value.
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 dicrotic notch characteristic point is provided,
Include:
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: the starting of the radial pulse waveform is determined based on the maximum of the differential signal
Point;
Period determines step: being based on the starting point, determines the period of the radial pulse waveform;
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;
Feature point recognition step: the dicrotic notch characteristic point of calibrated waveform is identified.
This method solve being mentioned the problem of there is larger limitations in existing measurement method, i.e., it is needed to acquired
Human body radial artery pulse wave data carry out that frequency-domain transform, small echo processing, second differnce to be to determine human body radial artery pulse wave wave
It is more to occupy hardware computing resource for the characteristic point of shape, the relatively complicated complexity of data handling procedure.Method benefit provided herein
Convolutional filtering processing is carried out to pulse wave with a specific two-dimensional array;Maximum is then obtained using first-order difference method
Point;After going baseline drift and calibration waveform to radial pulse waveform, finally according to different pulse cycles corresponding
Wave period scans in region, so that it is determined that dicrotic notch F point.The data processing method reduce to a certain extent for
The treatment process of pulse wave data reduces the consumption of computing resource, it is easier to the extraction of wave character point, so that subsequent
It is more convenient using the analysis treatment process of pulse wave.
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.
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, wherein the calculation formula of slope are as follows:
The waveform offset calculation formula are as follows:
Δ vtSlope=vtSlope [i] * (j-m_vtFootPos [i])
Go baseline drift treated waveform are as follows:
Normalization [j]=vtSmoothedData [j]-vtSmoothedData [m_vtFootPos [i]]-Δ
vtSlope
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.
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, the Feature point recognition step includes: that the region of search dicrotic notch characteristic point is determined according to wave period,
In the case where there is differential signal positive zero crossing in this region, using the minimum point of near zero-crossing point as dicrotic notch feature
Point;In the case where zero crossing is not present, using the region incurvature signal maximum point as dicrotic notch characteristic point.
Optionally, the calculation formula of the curvature signal are as follows:
Wherein, vtSmoothedData indicates filtered radial pulse waveform;VtSecDrtSig indicates that the oar is dynamic
The second differnce signal of arteries and veins pulse wave, i indicate i-th bit signal value.
Further aspect of the application provides a kind of device for identifying radial artery wave shape dicrotic notch characteristic point, packet
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 pulse based on the maximum of the differential signal
The starting point of waveform;
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 dicrotic notch characteristic point of calibrated waveform.
The device reduces the treatment process for pulse wave data to a certain extent, reduces disappearing for computing resource
Consumption, it is easier to the extraction of wave character point, so that subsequent more convenient using the analysis treatment process of pulse wave.
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 method for running the identification radial artery wave shape dicrotic notch characteristic point according to the application one embodiment
Computer installation hardware structural diagram;
Fig. 2 is showing according to the method for the identification radial artery wave shape dicrotic notch characteristic point of one embodiment of the application
Meaning property flow chart;
Fig. 3 is the method according to the identification radial artery wave shape dicrotic notch characteristic point of another embodiment of the application
Schematic flow chart;
Fig. 4 is the method according to the identification radial artery wave shape dicrotic notch characteristic point of another embodiment of the application
Schematic flow chart;
Fig. 5 is dicrotic notch F point in the application one embodiment in the upper schematic diagram set of radial pulse waveform;
Fig. 6 is showing according to the device of the identification radial artery wave shape dicrotic notch characteristic point of one embodiment of the application
Meaning property block diagram;
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 reality of method for identifying radial artery wave shape dicrotic notch characteristic point is additionally provided
Apply example, it should be noted that step shown in the flowchart of the accompanying drawings can be in the meter of such as a group of computer-executable instructions
It is executed in calculation machine system, although also, logical order is shown in flow charts, and it in some cases, can be with difference
Shown or described step is executed in sequence herein.
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 is the identification radial artery wave shape dicrotic notch characteristic point run according to the application one embodiment
The computer installation hardware structural diagram of method.As shown in Figure 1, computer installation 10 (or mobile device 10) may include one
A or multiple processors (use 102a, 102b ... ..., 102n is shown, and processor can include but is not limited to micro process in figure
The processing unit of device MCU or programmable logic device FPGA etc.), memory 104 for storing data and for communicating function
The transmission module of energy.It in addition to this, can also include: display, input/output interface (I/O interface), universal serial bus
(USB) port (a port that can be used as in the port of I/O interface is included), network interface, power supply and/or camera.This
Field those of ordinary skill is appreciated that structure shown in FIG. 1 is only to illustrate, and does not cause to the structure of above-mentioned electronic device
It limits.For example, computer installation 10 may also include than shown in Fig. 1 more perhaps less component or have with shown in Fig. 1
Different configurations.
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, it is special that this application provides a kind of identification radial artery wave shape characteristic point, especially gorges
The method for levying point F point.Fig. 2 is the side according to the identification radial artery wave shape dicrotic notch characteristic point of one embodiment of the application
The schematic flow chart of method.This method may comprise steps of:
S100 filter step: radial pulse waveform is filtered;
S200 differential signal generation step: generating the differential signal of the radial pulse waveform, determines the difference letter
Number maximum;
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: being based on the starting point, determines the period of the radial pulse waveform;
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;
S600 Feature point recognition step: the dicrotic notch characteristic point of calibrated waveform is identified.
A kind of method identifying radial artery wave shape characteristic point dicrotic notch F point provided by the present application, this method solve
Mentioned in existing measurement method the problem of there is larger limitations, i.e., it is needed to obtained human body radial pulse wave number
The characteristic point of human body radial artery pulse wave waveform, data processing are determined according to progress frequency-domain transform, small echo processing, second differnce
It is more to occupy hardware computing resource for the relatively complicated complexity of journey.Method provided herein utilizes a specific two-dimensional array
Convolutional filtering processing is carried out to pulse wave;Maximum point is then obtained using first-order difference method;To radial artery pulse wave
After shape removes baseline drift and calibration waveform, finally searched in corresponding wave period region according to different pulse cycles
Rope, so that it is determined that dicrotic notch F point.The data processing method reduces to a certain extent for the processed of pulse wave data
Journey reduces the consumption of computing resource, it is easier to the extraction of wave character point, so that at the subsequent analysis using pulse wave
Reason process is more convenient.
Fig. 3 is the method according to the identification radial artery wave shape dicrotic notch characteristic point of another embodiment of the application
Schematic flow chart.Optionally, the smothing filtering step includes: by the radial pulse waveform and preset two-dimensional array
Carry out convolutional filtering processing.For example, the array can be the two-dimensional array of a 11*14: QG [11] [14].By pulse wave with
Filter coefficient carries out convolutional filtering processing, to obtain filtered pulse wave.
Wherein, each element can be and specifically count as shown in following table one in the specific two-dimensional array QG [11] [14]
Value:
Table one
35 |
17 |
12 |
-3 |
|
|
|
|
|
|
|
|
|
|
21 |
7 |
6 |
3 |
-2 |
|
|
|
|
|
|
|
|
|
231 |
59 |
54 |
39 |
14 |
-21 |
|
|
|
|
|
|
|
|
429 |
89 |
84 |
69 |
44 |
9 |
-21 |
|
|
|
|
|
|
|
143 |
25 |
24 |
21 |
16 |
9 |
0 |
-11 |
|
|
|
|
|
|
1105 |
167 |
162 |
147 |
122 |
87 |
42 |
-13 |
-78 |
|
|
|
|
|
323 |
43 |
42 |
39 |
34 |
27 |
18 |
7 |
-6 |
21 |
|
|
|
|
2261 |
269 |
264 |
249 |
224 |
189 |
144 |
89 |
24 |
-51 |
-136 |
|
|
|
3059 |
329 |
324 |
309 |
284 |
249 |
204 |
149 |
84 |
9 |
-76 |
-171 |
|
|
8059 |
79 |
78 |
75 |
70 |
63 |
54 |
43 |
30 |
15 |
-2 |
-21 |
-42 |
|
5175 |
467 |
462 |
447 |
422 |
387 |
322 |
287 |
222 |
147 |
62 |
-33 |
-138 |
-253 |
Referring to figs. 2 and 3, optionally, in the differential signal generation step, the calculating for generating differential signal is public
Formula 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.
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.Wherein, it is the smallest minimum point calculation method: to find region of search range value
Point is minimum point.The process for determining starting point by differential signal maximum point is: by differential signal maximum point (one
As appear in starting point after), search for 200ms forward and determine region of search, in the region of search, when there is differential signal just
When to zero crossing, as waveform starting point;When positive zero crossing is not present, using the region of search minimum point as starting
Point.
Referring to figs. 2 and 3, optionally, it is determined in step in the period, it is true by the range difference for calculating two starting points
The fixed 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.
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 the period determines step, this method further includes that the peak point determines step, which includes In
Acquisition waveforms maximum of points in the starting point to 2/3 section in the period, using the waveform maximum of points as the peak
Value 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])
Go baseline drift treated waveform are as follows:
Normalization [j]=vtSmoothedData [j]-vtSmoothedData [m_vtFootPos [i]]-Δ
vtSlope
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.
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 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 step includes: that the region of search dicrotic notch characteristic point is determined according to wave period,
In the case where there is differential signal positive zero crossing in this region, using the minimum point of near zero-crossing point as dicrotic notch feature
Point;In the case where zero crossing is not present, using the region incurvature signal maximum point as dicrotic notch characteristic point.
Fig. 4 is the method according to the identification radial artery wave shape dicrotic notch characteristic point of another embodiment of the application
Schematic flow chart.For example, three kinds of situations can be divided into according to the difference of pulse cycle size, and in each case
Region of search is all different;Then according to whether there are differential signal positive zero crossings to determine dicrotic notch F point: then existing if it exists
Search minimum point is as F point in the region 15ms before and after this zero crossing, if it does not exist then according to relevant calculation formula in the region
Interior search curvature signal maximum of points is as F point.The final dicrotic notch F point for obtaining waveform.
(1) when wave period is when 625ms is between 925ms, first from 0.32 wave period to 0.5 wave period
F point is searched in region;When there are differential signal positive zero crossing, minimum point is searched in the region 15ms before and after this zero crossing
As F point;Zero crossing if it does not exist searches for curvature signal maximum of points as F point;
(2) as wave period > 925ms, F is searched in from 0.28 wave period to 0.45 wave period region first
Point;When there are differential signal positive zero crossing, region of search minimum point is as F in this zero crossing, the front and back region 15ms
Point;Zero crossing if it does not exist searches for curvature signal maximum of points as F point;
(3) as wave period < 625ms, F is searched in from 0.36 wave period to 0.52 wave period region first
Point;When there are differential signal positive zero crossing, with this zero crossing, minimum point is searched in the region 15ms of front and back as F point;When
There is no when differential signal positive zero crossing, search for curvature signal maximum of points as F point.
The calculation formula of the curvature signal are as follows:
Wherein, vtSmoothedData indicates filtered radial pulse waveform;VtSecDrtSig indicates that the oar is dynamic
The second differnce signal of arteries and veins pulse wave, i indicate i-th bit signal value.
Fig. 5 is dicrotic notch F point in the application one embodiment in the upper schematic diagram set of radial pulse waveform.The figure
In include be complete signal period radial pulse waveform, be mainly made of ascending branch and descending branch.Ascending branch and descending branch composition master
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.Tightly
Connecing the appearance of dicrotic notch 3 is dicrotic wave 4, also known as dicrotic wave.Above-described wave and gorge are the main components for constituting pulse wave.
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;F point is 3 valley point of dicrotic notch, represents LV Diastolic starting point, is that main 1 decent of wave and dicrotic wave 4 rise
The downward incisura wave of the waveform that branch is constituted.When human body artery tube wall smooth muscle is less, elastomer is more, pulse wave is formed
Main wave C high and it is sharp, 4 trough F of dicrotic wave is big and obvious due to blood backflow impact strength.Conversely, ductus arteriosus wall smooth muscle
Increase it is less with tube wall elasticity fiber, back wave velocity of wave increase, increase the relatively main 1 wave crest location of C of wave of 4 trough F of dicrotic wave gradually
It is high close, different degrees of fusion is finally presented.It can be seen that the height fluctuations of pulse wave characteristic point are able to reflect human body
The variation of vascular resistence and vessel wall elasticity.
Embodiment 2
According to the embodiment of the present application, a kind of device for identifying radial artery wave shape dicrotic notch characteristic point is additionally provided, it should
Device is device corresponding with method described in embodiment 1.Fig. 6 is the identification pressure of the radial artery according to one embodiment of the application
The schematic block diagram of the device of Reeb shape dicrotic notch characteristic point.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;
Feature point recognition module 600 is disposed for identifying the dicrotic notch characteristic point of calibrated waveform.
The device carries out convolutional filtering processing to pulse wave using a specific two-dimensional array;Then use a scale
Method is divided to obtain maximum point;After going baseline drift and calibration waveform to radial pulse waveform, finally according to different
Pulse cycle scans in corresponding wave period region, so that it is determined that dicrotic notch F point.The data processing method is certain
Reduce the treatment process for pulse wave data in degree, reduces the consumption of computing resource, it is easier to wave character point
Extraction so that subsequent more convenient using the analysis treatment process of pulse wave.
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 differential signal waveform
Maximum point.
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 positive mistake based on the maximum point of the differential signal
Zero point region of search is searched in the region of search, 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])
Go baseline drift treated waveform are as follows:
Normalization [j]=vtSmoothedData [j]-vtSmoothedData [m_vtFootPos [i]]-Δ
vtSlope
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 is used for: the region of search dicrotic notch characteristic point is determined according to wave period,
In the case where there is differential signal positive zero crossing in this region, using the minimum point of near zero-crossing point as dicrotic notch feature
Point;In the case where zero crossing is not present, using curvature signal maximum of points as dicrotic notch characteristic point.
The calculation formula of the curvature signal are as follows:
Wherein, vtSmoothedData indicates filtered radial pulse waveform;VtSecDrtSig indicates that the oar is dynamic
The second differnce signal of arteries and veins pulse wave, i indicate 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.
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 (Read-Only Memory, ROM), arbitrary access are deposited
Reservoir (Random Access Memory, RAM), 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.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
SolidStateDisk (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: magnetictape), floppy disk (English: floppydisk), CD (English:
) and any combination thereof opticaldisc.
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.