A kind of pulse wave transfer rate detection method and system based on image ghost image
Technical field
The present invention relates to medical signals process field, more particularly, to a kind of pulse wave transfer rate based on image ghost image
Detection method and system.
Background technology
Pulse wave transmission speed (Pulse Wave Velocity, PWV) refers to that pulse wave is passed by an ad-hoc location of tremulous pulse
Cast to the speed of another ad-hoc location.The index arterial stiffness of arterial elasticity and expansivity increases rises can PWV numerical value, and
Arterial elasticity weakens declines can PWV numerical value.So by measuring PWV numerical value, we can estimate the elasticity of tremulous pulse and can expand
Extensional.The very useful Noninvasive index of PWV numerical value or one so that arteriosclerosiss, hypertension, hyperlipemia, glycosuria
The early discovery and diagnosis of various blood vessel kind diseases such as disease and nephropathy is possibly realized.
The detection of arteriosclerosis is increasingly paid attention to by medical circle, and its testing result can remind patient oneself that tremulous pulse is hard
The situation of change, in time treatment or the corresponding exercise of increase, change dietary structure.The computing formula of PWV is very simple, the meter of PWV
Calculation method such as formula:PTT represents pulse wave transmission time, and L represents distance, does not have anything to be worth research, can
It is simpler formula, the requirement to calculating parameter is higher, and the measurement of wherein time parameter requires that detecting instrument is provided
Good detection waveform, this point are extremely difficult, and the corresponding signal of human body carries very big noise, how to carry critical noisy
The problem that actual signal is medical detecting Instrument facing is extracted in signal.The multiplex general filter of traditional PWV detecting instruments
Wave technology, or probability statistics technology is eliminating noise as far as possible.But filtering technique can bring very big time-lag, PWV
Detection is mainly the time in counting period, and retardation can bring very big calculation error, and probability statistics need substantial amounts of sample data,
Time of measuring can be not only greatly increased, and will effectively remove the sample size of signal noise and be practically impossible to what is realized.
The content of the invention
The purpose of the present invention is exactly to provide to overcome above-mentioned prior art to there is not accurate enough the defect of testing result
A kind of pulse wave transfer rate detection method and system based on image ghost image, carries out PWV detections using image ghost image technology.It is first
Pretreatment is carried out to input signal (ECG signal and pulse wave signal) first;Then the peak position feature for carrying out ECG signal is carried
Take, so that image ghost image technology below is used;Continuous signal is divided into by corresponding independence using image ghost image technology finally
Module, carries out time delay estimation with cross-power spectrum method to standalone module, not only easy to operate, and testing result is accurately reliable.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of pulse wave transfer rate detection method based on image ghost image includes:
Step S1:Collection electrocardiosignal and pulse wave signal;
Step S2:After respectively step S1 center telecommunications number and pulse wave signal are carried out with pretreatment and period position mark,
Extract multiple periodic sequences of electrocardiosignal and pulse wave signal;
Step S3:The peak position feature of multiple periodic sequences of difference obtaining step S2 center telecommunications number and pulse wave signal,
And in obtaining electrocardiosignal and pulse wave signal according to peak position feature each periodic sequence position;
Step S4:Based on image ghost image technology according to peak position feature respectively by step S3 center telecommunications number and pulse wave signal
In each periodic sequence starting point align, the value superposition of each the periodic sequence relevant position after alignment obtains electrocardiosignal and arteries and veins
The respective monocycle sequence waveform of ripple signal of fighting, obtains the power spectrum peak of two monocycle sequence waveforms, respectively by electrocardio
Time interval corresponding to signal and pulse wave signal power spectrum peak propagates transmission time PTT as pulse wave, by formula
(1) obtain pulse wave transfer rate PWV:
Wherein, L is the distance between tremulous pulse had both been pinpointed at two, obtains the value by measuring body surface.
Pretreatment in step S2 is included at the amplitude of the amplitude to electrocardiosignal, Phase Processing and pulse wave signal
Reason;
The amplitude of electrocardiosignal, Phase Processing include:Take one group of electrocardiosignal, obtain this group of electrocardiosignal meansigma methodss,
Maximum and minima, judge meansigma methodss to the width of maximum whether more than the width of meansigma methodss to minima, the if so, then heart
Signal of telecommunication input is forward signal, otherwise, carries out value to electrocardiosignal and negates, and the electrocardiosignal after judgement is normalized place
Reason;
The amplitude of pulse wave signal is processed to be included:Take one group of pulse wave signal to be normalized.
In step S2, period position mark includes:According to the Changing Pattern of signal, at the end of each signal period
It is identified, erasure signal cycle incomplete part obtains complete periodic sequence.
Corresponding 10~20 periodic sequences are extracted in step S2.
In step S3, peak position feature includes number and the corresponding position of crest and trough.
In step S4, power spectrum peak is highest crest characteristic point.
It is a kind of to realize including based on the pulse wave transfer rate detecting system of image ghost image for said method:
Data acquisition module, for gathering electrocardiosignal and pulse wave signal;
Pretreatment module, for the output of receiving data acquisition module, is carried out to electrocardiosignal and pulse wave signal respectively
After pretreatment and period position mark, multiple periodic sequences of electrocardiosignal and pulse wave signal are extracted;
Crest detection module, for receiving the output of pretreatment module, obtains electrocardiosignal and pulse wave signal respectively
The peak position feature of multiple periodic sequences, and each periodic sequence in obtaining electrocardiosignal and pulse wave signal according to peak position feature
Position;
Image ghost image module, for the output of received wave blob detection module, based on image ghost image technology according to peak position feature
Respectively each periodic sequence starting point in electrocardiosignal and pulse wave signal is alignd, each the periodic sequence corresponding positions after alignment
The value superposition put obtains electrocardiosignal and the respective monocycle sequence waveform of pulse wave signal, obtains two monocycle sequence waveforms
Power spectrum peak, using the time corresponding to electrocardiosignal and pulse wave signal power spectrum peak as pulse wave propagate pass
Defeated time PTT, and then obtain pulse wave transfer rate PWV.
The data acquisition module includes:
Sensor unit, for gathering electrocardiosignal and pulse wave signal;
Signal condition unit, for receiving the output of sensor unit, is carried out successively to electrocardiosignal and pulse wave signal
One-level differential amplification, low-pass filtering, high-pass filtering and two grades of amplifications;
AD conversion unit, manages the output of unit for receiving signal, AD is carried out to electrocardiosignal and pulse wave signal and is turned
Pretreatment module is sent to after changing.
The pretreatment module, crest detection module and image ghost image module are by the host computer reality containing virtual instrument
It is existing.
Compared with prior art, the present invention has advantages below:
1) the inventive method once only needs the input signal sample of limited periodic sequence, and sample is few, so as to fast
Speed realizes the collection of signal, while rejecting incomplete periodic sequence, Jing after image ghost image technical finesse, output waveform change is bright
It is aobvious, and with comparative, phase is carried out using the monocyclic waveform of crosspower spectrum function pair ECG signal and pulse wave signal
Answer the time interval of characteristic point to estimate, from the successful waveform of ghost image can be accurately obtained two signal highest crest characteristic points when
Between be spaced, the time interval is exactly required pulse propagation time, and the method compares traditional PWV detection methods, Neng Gouyou
Effect overcomes effect of noise, will not cause temporal retardation again, improve PWV accuracy of detection.
2) the inventive method carries out phase judgment and periodic sequence station location marker to electrocardiosignal and pulse wave signal, prevents
There is Electrocardiographic reverse signal and affect phase place positioning and reject incomplete periodic sequence, reduce calculation error.
3) peak position feature extraction of the inventive method to electrocardiogram and pulse wave signal, obtain input signal crest ripple
Paddy number and corresponding position, for dividing the position of each periodic sequence, it is ensured that monocycle sequence divides accurate.
4) the inventive method extracts electrocardiosignal respectively and the multiple complete periodic sequence of pulse wave signal is overlapped
Process, obtain amplitude and expand many times of electrocardiosignal and the monocycle waveform of pulse wave signal, the monocycle letter obtained after addition
The change of number wave-shape amplitude is obvious, realizes the reinforcement to useful signal, weakens impact of the noise to signal.
5) the inventive method devises host computer structure, practicality using image ghost image technology and Virtual instrument LabVIEW
By force, it is portable good.
6) present system adopts electrocardioelectrode and pulse wave sensor collection electrocardiosignal and pulse wave signal, realizes people
The collection of body millivolt level signal;Devising signal conditioning circuit carries out denoising and processing and amplifying to signal, is met AD samplings
The analogue signal of condition;Using the USB-6008 data collecting cards of NI companies, gather while realizing multiple signals and be sent to
Host computer.The host computer based on virtual instrument is devised simultaneously, makes full use of virtual instrument strong on signal testing and analysis
Big advantage so that the detection means construction cycle is short, flexible and efficient, realizes image ghost image technology using virtual instrument, effectively overcomes
Effect of noise, will not cause temporal retardation again, improve PWV accuracy of detection.
Description of the drawings
Fig. 1 is data acquisition module structural representation in present system;
Fig. 2 is the host computer structural representation containing virtual instrument in present system;
Fig. 3 is signal condition unit circuit diagram in present system;
Fig. 4 is the signal waveform schematic diagram that collects in the inventive method;
Fig. 5 is the signal waveform schematic diagram after the inventive method acceptance of the bid " | " is processed;
Fig. 6 is the electrocardiosignal and pulse wave signal waveform diagram after amplitude phase process in the inventive method;
Fig. 7 is the inventive method medium wave peak detection function schematic diagram;
Fig. 8 is waveform alignment function schematic diagram in the inventive method;
Fig. 9 is crosspower spectrum function schematic diagram in the inventive method;
Figure 10 is the signal waveform schematic diagram after ghost image in the inventive method.
In figure:1st, data acquisition module, 2, host computer, 11, sensor unit, 12, signal condition unit, 13, AD conversion
Unit, 21, pretreatment module, 22, crest detection module, 23, image ghost image module, 121, one-level differential amplifier circuit, 122,
Low-pass filter circuit, 123, high-pass filtering circuit, 124, second amplifying circuit.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.The present embodiment is with technical solution of the present invention
Premised on implemented, give detailed embodiment and specific operating process, but protection scope of the present invention be not limited to
Following embodiments.
As shown in Figure 1 and Figure 2, a kind of pulse wave transfer rate detecting system based on image ghost image includes data acquisition module
2 two large divisions of block 1 and host computer.Data acquisition module 1 is used to gather electrocardiosignal and pulse wave signal, comprising three units,
Respectively sensor unit 11, signal condition unit 12 and AD conversion unit 13;Host computer 2 containing virtual instrument includes three
Individual processing module, respectively pretreatment module 21, crest detection module 22 and image ghost image module 23.
Operation principle:First, the collection electrocardiosignal of sensor unit 11 and pulse wave signal, output signal is with noise
With the analogue signal of DC component, need to obtain pure electrocardial analog signal and pulse wave analogue through signal condition unit 12
Signal, the signal after conditioning are AD converted by AD conversion unit 13 again, and AD conversion unit 13 is by the electrocardio for collecting and pulse
The digital signal of ripple uploads to host computer 2, and host computer 2 is carried out to electrocardiosignal, pulse wave signal using integrated virtual instrument
Process and analyze, and the result after Treatment Analysis is exported.Below each group is specifically described into part:
Sensor unit 11 includes ECG electrode and fluctuation sensor, and fluctuation sensor adopts PVDF piezoelectric membranes, surveys
Amount 0~50g of scope, trueness error≤5%F.S, sensitivity >=20mV/F.S, the pulse wave signal of extraction are truer, are surveyed
Human pulse ripple signal is millivolt level signal.
The outfan of the sensor input connection unit 11 of signal condition unit 12.12 points of signal condition unit is electrocardio
Signal conditioning circuit and pulse wave signal modulate circuit, two signal conditioning circuit structures are similar.Pulse wave signal modulate circuit
As shown in figure 3, including the one-level differential amplifier circuit 121, low-pass filter circuit 122,123 and of high-pass filtering circuit that are sequentially connected
Second amplifying circuit 124.One-level differential amplifier circuit 121 includes AD620 difference amplifiers and its peripheral circuits (including resistance R1
With electric capacity C1, C2).Low-pass filter circuit 122 include OPA2277 chips and its peripheral circuits (including resistance R2, R3, R4, R5 and
Electric capacity C3), high-pass filtering circuit 123 include OPA2277 chips and its peripheral circuits (including resistance R6, R7, R8 and electric capacity C4,
C5), second amplifying circuit 124 includes OPA2277 chips and its peripheral circuits (including resistance R9, R10, R11, R12 and electric capacity
C6、C7).Wherein, OPA2277 chips are double operational circuits, and U2A and U2B is the different piece in same chip.Two signals
The collection signal of ECG electrode, fluctuation sensor by differential amplifier pre-amp output signal, is improved by modulate circuit
Signal to noise ratio, then the high-frequency noise in signal is removed by low-pass filter circuit 122, then letter is removed by high-pass filtering circuit 123
DC component in number, then the analogue signal of AD sampling conditions is met by second amplifying circuit 124.
The input of AD conversion unit 13 connects the outfan of signal condition unit 12, and outfan connects the defeated of host computer 2
Enter end.USB-6008 data collecting card of the AD conversion unit 13 using NI companies.The card is many work(based on USB2.0 serial ports
Energy data collecting card, has 8 simulation input terminals to can be used to receive electrocardiosignal and pulse wave signal, per road input range be ±
10V, 10Ks/s, 12 bit resolution.It converts analog signals into digital signal and sends into host computer 2 by USB data line
USB interface, is further processed.
Host computer 2 is the computer for being integrated with the LabVIEW software platforms based on virtual instrument Programming in Digital Image Processing, is realized to letter
The work(such as number process, analysis, display, storage, wave shape playback, signal smoothing filtering, Periodic identification, baseline adjustment, Feature point recognition
Can, make full use of powerful advantages of the host computer 2 containing virtual instrument on signal testing and analysis so that detection means is developed
It is cycle is short, flexible and efficient.As shown in Fig. 2
Pretreatment module 21, for the output of receiving data acquisition module 1, enters to electrocardiosignal and pulse wave signal respectively
After row pretreatment and period position mark, electrocardiosignal and respective 10~20 periodic sequences of pulse wave signal are extracted, so
Accuracy of detection is may insure, reduces time of measuring, in the present embodiment, choose 15 periodic sequences;
Crest detection module 22, for receiving the output of pretreatment module 21, to many of electrocardiosignal and pulse wave signal
The peak position feature of individual periodic sequence, and the position of each periodic sequence is obtained according to peak position feature;
Image ghost image module 23, for the output of received wave blob detection module 22, based on image ghost image technology according to peak position
Each periodic sequence starting point of electrocardiosignal and pulse wave signal is alignd by feature respectively, each the periodic sequence phase after alignment
Answer the value superposition of position to obtain electrocardiosignal and the respective monocycle sequence waveform of pulse wave signal, obtain two monocycle sequences
The power spectrum peak of waveform, passed the time corresponding to electrocardiosignal and pulse wave signal power spectrum peak as pulse wave
Broadcast transmission time PTT and (time of individual features point is obtained from the oscillogram of monocyclic electrocardiosignal and pulse wave signal
Interval), and then obtain pulse wave transfer rate PWV.
As shown in Figure 1 and Figure 2, the pulse wave transfer rate that a kind of utilization said detecting system is carried out based on image ghost image is examined
Survey method includes:
Step S1:Data acquisition module 1 gathers electrocardiosignal and pulse wave signal.Electrocardiosignal is produced when referring to cardiomotility
Electro-physiological signals, pulse wave signal refer to radial artery blood flow produce physiological signal.
The present invention acquires the signal of the upper and lower limb of human body and aorta, and the signal for collecting is as shown in figure 4, from top to bottom
It is followed successively by carotid artery waveform, radial artery waveform, femoral artery waveform, ankle red artery waveform, phonocardiogram, electrocardiogram.Can be with from waveform
Find out that collectable signal has very big noise.The method that the present invention uses electrocardiograph pulse, it is only necessary to electrocardiosignal and
Tremulous pulse end pulse wave signal at one is processed.The present invention extracts electrocardiosignal and noise is relatively small and easy measurement
The signal of the pulse wave of radial artery, that is, extract the second row and the 6th row signal.
Step S2:Pretreatment module 21 is carried out to the electrocardiosignal and pulse wave signal of the output of data acquisition module 1 respectively
After pretreatment and period position mark, electrocardiosignal and respective 15 periodic sequences of pulse wave signal are extracted.Wherein, pretreatment
Process including the amplitude of the amplitude to electrocardiosignal, Phase Processing and pulse wave signal, specially:
Take one group of ECG signal and one group of pulse wave signal, obtain the meansigma methodss of this group of ECG signal, maximum and
Minima, which width is big to maximum and meansigma methodss to minima to compare meansigma methodss, if meansigma methodss are larger to maximum width,
Then for forward signal, output to next step otherwise carries out value to input ECG signal and negates for input, finally to electrocardiogram and
The amplitude that pulse wave signal is normalized is processed.
Period position mark includes:Electrocardiosignal and pulse wave signal are compared to numerical value in respective signal respectively, are obtained
To the Changing Pattern of signal, it is identified at the end of each signal period, erasure signal cycle incomplete part, final difference
Obtain the complete periodic sequence of electrocardiosignal and pulse wave signal.As shown in figure 5, periodic sequence is designated to periodic sequence mark
" | " is processed, and pulse wave train that will be original is divided into several subsequences, removes two subsequences of head and the tail, cannot reject completely
Whole part, each middle subsequence is a complete pulse cycle, can individually extract a complete cycle
Sequence, realizes the Periodic identification to original pulse wave.Corresponding 15 complete periodic sequences are extracted for image ghost image technical office
Reason is prepared, and so may insure accuracy of detection, reduces time of measuring.
To electrocardiosignal and radial artery signal, the result of the two signal ghost images after pretreatment is as shown in Figure 6.
Step S3:Crest detection module 22 obtains the electrocardiosignal and pulse wave signal of the output of pretreatment module 21 respectively
The peak position feature of periodic sequence, and the position of each periodic sequence is obtained according to peak position feature.Peak position feature includes crest and ripple
The number of paddy and corresponding position.
Crest detection function of the crest detection module 22 using standard, is the signal processing instrument of labview virtual softwares
A waveforms detection control in bag, can be provided the crest number of waveform, the amplitude vector of corresponding crest by the control
With the position vector of crest, crest detection function is as shown in Figure 7.
Step S4:Image ghost image module 23 is based on image ghost image Technology design, according to peak position characteristic use waveform alignment letter
The electrocardiosignal and each periodic sequence of pulse wave signal that crest detection module 22 is exported by number respectively carries out starting point alignment,
Time is identical with component value, is easy to two signals to compare, and the value superposition of each the periodic sequence relevant position after alignment obtains electrocardio
Signal and the respective monocycle sequence waveform of pulse wave signal, then two monocycle sequence waveforms are obtained by crosspower spectrum function
Power spectrum peak (power spectrum peak be highest crest characteristic point), by electrocardiosignal and pulse wave signal power spectrum highest
The point corresponding time propagates transmission time PTT as pulse wave, obtains pulse wave transfer rate PWV by formula (1):
Wherein, L is the distance between tremulous pulse had both been pinpointed at two, and distance of such as radial artery to two measurement point of femoral artery is logical
Cross measurement body surface and obtain the value.
Image ghost image technology refers on one or more addition of waveforms to another waveform, forms the effect of ghost image, fold
Each waveform of shadow must be alignment in abscissa time point, before ghost image must also carry out denoising to waveform.
Waveform alignment function is as shown in figure 8, control for the waveform conditioning in the signal processing tool kit of labview virtual softwares
Part, as a result middle waveform A outputs are identical with component value with the time of waveform B output.Waveform alignment function is by electrocardiosignal and pulse
The starting point alignment of the monocycle waveform of ripple signal, in order to avoid image offset occur and cause resultant error.
Crosspower spectrum function is as shown in figure 9, for the analysis of spectrum control in the signal processing tool kit of labview virtual softwares
Part, crosspower spectrum react mutual relation of the input signal in different frequency point, and its phase place is a kind of conventional delay time estimation method,
The method is in the environment of small noise and moderate reverberation once, it is possible to obtain more accurate time delay estimated value.Cross-power spectrum method
Obtain the time interval of electrocardiosignal and pulse wave signal individual features point.
15 complete periodic sequences of electrocardiogram and pulse wave signal are overlapped process, obtain amplitude and expand 15 times
The monocycle waveform of ECG signal and pulse wave signal, the monocycle signal wave-shape amplitude change obtained after addition are obvious, real
Show the reinforcement to useful signal, weaken impact of the noise to signal.It is superimposed successful waveform as shown in Figure 10, solid line is pulse
Ripple signal, dotted line are electrocardiosignal, it can be seen that last oscillogram signal is very clear, also no retardation, are so calculated
The phase place retardation of signal is just very accurate, so as to substantially increase certainty of measurement and repeatability.Obtained by crosspower spectrum function
The time interval of electrocardiosignal and pulse wave signal individual features point can just greatly improve computational accuracy calculating PTT.
To sum up, the method that adopts of the present invention is for electrocardiograph pulse method, by collection electrocardiosignal and tremulous pulse end pulse signal,
The time interval of electrocardio and pulse signal individual features point is calculated, transmission time is represented.The treatment technology for adopting is folded for image
Shadow technology, to gather electrocardiosignal and pulse wave signal carry out a series of process, finally obtain monocyclic electrocardiosignal and
Pulse wave signal, obtains the time interval of individual features point, i.e. pulse wave using crosspower spectrum function after alignment image ghost image
Transmission time.