CN106963424A - Detect the viscoelastic system and method for arteries - Google Patents
Detect the viscoelastic system and method for arteries Download PDFInfo
- Publication number
- CN106963424A CN106963424A CN201710153465.4A CN201710153465A CN106963424A CN 106963424 A CN106963424 A CN 106963424A CN 201710153465 A CN201710153465 A CN 201710153465A CN 106963424 A CN106963424 A CN 106963424A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- skin
- phase
- estimation
- signal
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4411—Device being modular
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- 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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
Abstract
The viscoelastic method of arteries is detected the present invention relates to one kind, including:The outside oscillator of control is producing vibration on the skin of blood vessel;Ultrasonic probe detects pulse to skin transmitting, record the pulse echo signal of vascular tissue in skin under each passage, and the displacement for obtaining vascular tissue in the skin of each passage is calculated, and the pulse wave velocity in displacement is filtered out, while electrocardiogram acquisition machine gathers electrocardiosignal;The instant exposure speed of strain wave is tried to achieve, the instantaneous viscoelastic coefficient of blood vessel is estimated;Under cardiac cycle, the repeatedly instantaneous viscoelastic coefficient of estimation blood vessel, and the viscoelastic coefficient of estimation and the electrocardiosignal of collection is registering in time obtains under cardiac cycle not the modulus of elasticity and viscous modulus of blood vessel in the same time.The invention further relates to one kind detection viscoelastic system of arteries.The present invention replaces ultrasonic probe using outside oscillator, safer, can both form an independent equipment, and can also be used as an additional function modules.
Description
Technical field
Detect that the viscoelastic system and method for arteries, more particularly to one kind are quantitatively detected in body the present invention relates to one kind
The viscoelastic system and method for arteries.
Background technology
Cardiovascular and cerebrovascular disease is the general designation of cardiovascular and cranial vascular disease, is referred to because hyperlipidemia, blood are sticky, dynamic
Ischemic or hemorrhagic disease that heart, brain and body tissue caused by pulse atherosclerosis, hypertension etc. occur.Many institute's weeks
Know, angiocardiopathy is the No.1 killer of the mankind, the annual crowd for dying from cardiovascular and cerebrovascular disease in the whole world is up to 15,000,000 people, China
Cardiovascular and cerebrovascular sufferer is up to 2.9 hundred million, year Hospitalization expenses more than 70,000,000,000 yuan, the Disease Spectrum of cardiovascular disease is increasingly aggravated, into
For the great public health problem of China.
The mechanical characteristic of artery plays the role of a nucleus in cardiovascular function, between arterial stiffness and various angiocardiopathies
There is direct relation, be the important predictive factor of cardiovascular and cerebrovascular morbidity and mortality.Although arterial vascular viscoelastic change
Change usually systemic, and influence whole vascular system, but the arterial vascular hardness of diverse location increased the age, height
The reaction of the factors such as blood pressure, gestation is also different, and often arteries is hardened, and aorta vessel softens, so accurate measurement
Local intra-arterial seems particularly significant.It is this non-because the viscoplasticity of arterial blood tube wall is with the nonlinear change that is changing into of blood pressure
Linear change is most important to the healthy of whole cardiovascular system, and can also detect that unstable in atherosclerosis
Fixed plate.
In order to obtain local, instantaneous, quantitative hardness measurement, it is impossible to the compression produced by cardiac pumping to blood vessel
Ripple, is once fluctuated because heart is about produced for one second, can be met using shearing wave elastogram of the ultrasound based on acoustic radiation force
Multiple transient measurement is carried out to vascular wall under cardiac cycle, but this method has certain limitation:First, acoustic radiation force
The security detected for blood vessel needs to be investigated;Second, the method must possess the THE ADIABATIC SHEAR IN of autonomous property right by the team
Ripple imaging device could realize that system complex involves great expense, and be not easy to promote.
The content of the invention
In view of this, it is necessary to provide a kind of security it is higher body quantitatively detect the viscoelastic system of arteries and
Method.
One kind detection viscoelastic system of arteries, what the system was electrically connected with including main control computer, with the main control computer
Transmit/receive module and electrocardiogram acquisition module, the outside oscillator being electrically connected with the transmit/receive module and ultrasound are visited
Head, the electrocardiogram acquisition machine being electrically connected with the electrocardiogram acquisition module, the main control computer includes processing module, estimation module
And registration module, wherein:The transmit/receive module is used to control outside oscillator under cardiac cycle in the skin close to blood vessel
It is upper to produce vibration;The transmit/receive module is additionally operable to drive ultrasonic probe under cardiac cycle to skin transmitting detection pulse,
Record the pulse echo signal of vascular tissue in the skin under each passage of ultrasonic probe;The processing module is used for according to above-mentioned
The pulse echo signal of vascular tissue in skin under each passage of record, calculating obtains blood vessel group in the skin of each passage
The displacement knitted, and filter out the pulse wave velocity in the displacement;The electrocardiogram acquisition module is used to trigger the heart under cardiac cycle
Electric harvester gathers electrocardiosignal;The estimation module is used to filter out each passage after pulse wave velocity according to obtained above
The displacement of vascular tissue in skin, tries to achieve the instant exposure speed of strain wave, estimates the instantaneous viscoelastic coefficient of blood vessel;It is described to match somebody with somebody
Quasi-mode block obtains week aroused in interest for the instantaneous viscoelastic coefficient of the blood vessel of estimation and the electrocardiosignal gathered is registering in time
The not modulus of elasticity of blood vessel, viscous modulus in the same time under phase.
Wherein, the ultrasonic probe is popped one's head in for channel ultrasonic.
Described processing module is used to calculate the displacement for obtaining vascular tissue in the skin of each passage, is specially:From inspection
Radiofrequency signal that vibration measuring member is received by ultrasonic front-end circuit amplification and the processing of analog-to-digital conversion, then by quadrature demodulation at
Reason, obtains the complex envelope of radiofrequency signal;To a certain test point, the complex envelope signal of echo-signal consecutive frame is carried out at cross-correlation
Reason, extracts the signal of tissue vibration in the different depth;Calculate the phase of test point vibration signal in some depth
Position, the phase of each frame subtracts the phase of consecutive frame, obtains the phase difference of each frame, and then obtains changing over time for vascular tissue
Displacement curve, at a distance of Z two test points, phase can be asked for by the displacement curve of a certain depth of two test points
Speed cL。
Described estimation module is based on Lamb wave mathematical modeling f (c, ω, μ1,μ2)=0, is specifically included:
I. assume<μ1,μ2>One group of initial value;
Ii. the estimation of phase velocity spectrum is obtained based on " criterion of least squares "I.e.:
Iii. the phase velocity spectrum of estimation is calculatedWith actual measurement phase velocity spectrum cphaseError sum of squares between (ω), table
It is shown as:
If iv. α (μ1,μ2) default precision is not up to, updateValue, return the i-th i step;
Otherwise calculate and stop, gainedIt is the estimate of the instantaneous viscoelastic coefficient of blood vessel.
Described registration module specifically for:The instantaneous viscoelastic coefficient of the blood vessel of estimation and the electrocardiosignal of collection are existed
Registration is done on time, to obtain blood vessel elasticity modulus not in the same time and viscous modulus under cardiac cycle.
The present invention also provides a kind of detection viscoelastic method of arteries, and this method comprises the following steps:A. control outer
Portion's oscillator is producing vibration on the skin of blood vessel;B. the ultrasonic probe of multichannel is to skin transmitting detection pulse, and record is every
Under individual passage in skin vascular tissue pulse echo signal, and calculate and obtain the position of vascular tissue in the skin of each passage
Move, and filter out the pulse wave velocity in the displacement, while electrocardiogram acquisition machine starts to gather electrocardiosignal;C. obtained according to above-mentioned
The displacement for filtering out vascular tissue in the skin of each passage after pulse wave velocity, try to achieve the instant exposure speed of strain wave, estimate
Count the instantaneous viscoelastic coefficient of blood vessel;D. by the instantaneous viscoelastic coefficient of the blood vessel of estimation and the electrocardiosignal of collection in time
Registration, obtains under cardiac cycle not the modulus of elasticity and viscous modulus of blood vessel in the same time.
Wherein, described step b includes:From the amplification for detecting first ultrasonic front-end circuit of radiofrequency signal process received that shakes
With the processing of analog-to-digital conversion, then process quadrature demodulation processing, the complex envelope of radiofrequency signal is obtained;To a certain test point, to echo
The complex envelope signal of signal consecutive frame carries out cross correlation process, extracts the signal of tissue vibration in the different depth;Calculate
The phase of test point vibration signal in some depth, the phase of each frame subtracts the phase of consecutive frame, obtains each frame
Phase difference, and then the displacement curve changed over time of vascular tissue is obtained, for two test points at a distance of Z, it can pass through
The displacement curve of a certain depth of two test points asks for phase velocity cL。
This method also includes:Within a cardiac cycle, step a is repeated to step c.
Described step c idiographic flows are as follows:
I. assume<μ1,μ2>One group of initial value;
Ii. the estimation of phase velocity spectrum is obtained based on " criterion of least squares "I.e.:
Iii. the phase velocity spectrum of estimation is calculatedWith actual measurement phase velocity spectrum cphaseError sum of squares between (ω), table
It is shown as:
If iv. α (μ1,μ2) default precision is not up to, updateValue, return the i-th i step;
Otherwise calculate and stop, gainedIt is the estimate of the instantaneous viscoelastic coefficient of blood vessel.
Described step d is specifically included:By the instantaneous viscoelastic coefficient of the blood vessel of estimation with collection electrocardiosignal when
Between on do registration, to obtain blood vessel elasticity modulus not in the same time and viscous modulus in a cardiac cycle.
The present invention detection viscoelastic system and method for arteries, is existed using outside oscillator with the mode that ultrasound is combined
Body quantitatively detects arteries viscoplasticity, and advantage is as follows:(1) present invention replaces ultrasonic probe using outside oscillator, quantitative in body
Measurement has higher security;(2) present invention can both form an independent equipment, simple portable, can also be used as one
Individual additional function modules, are carried in existing color ultrasound system and use, cost-effective, are easy to promote.
Brief description of the drawings
Fig. 1 detects the hardware architecture diagram of the viscoelastic system of arteries for the present invention;
Fig. 2 detects the flow chart of the viscoelastic method of arteries for the present invention.
Embodiment
The physical basis of the present invention are introduced as follows first:
The present invention serves as the motivational techniques of vibration source using outside oscillator (20HZ-1000HZ), i.e., put low frequency vibrator
Integrated in body surface, or by vibrator and ultrasonic probe, with convenient to operate.Shaken by vibration source to internal low frequency of propagating
It is dynamic, velocity information when detecting that the caused shearing wave of vibration is propagated in vascular tissue with ultrasonic probe.Active force puts on people
The body surface of body, to encourage subcutaneous vascular tissue, directly simulates doctor and presses tissue reality with finger with certain frequency
The process of inspection is applied, it is more directly perceived simple and security is higher, because blood vessel is a kind of tissue of thin layer tubular structure, shearing
When ripple is propagated in it, constantly the up-and-down boundary of cornea is had an effect, and produces refraction, reflection and the mould between shear wave and compressional wave
State is changed and forms Lamb wave.
For the Viscoelastic Thin Plates of leaching in a liquid, it is assumed that thin plate and the similar density of liquid and the compressible ripple in them
Wave number is much smaller than Lamb wave wave number, then the asymmetric Lamb wave Dispersion Equation in thin plate is as follows:
Wherein, kL=ω/cL, ω is angular frequency, cLIt is the Lamb wave phase velocity of frequency dependence,It is shearing
The wave number of ripple, μ is modulus of shearing, ρmIt is the density (similar density with water) of sample,H is equal to thickness of sample
Half.
The mechanical behavior of Viscoelastic Thin Plates is described by introducing Voigt models, modulus of shearing μ is then expressed as μ=μ1+i
ωμ2, wherein μ1And μ2It is modulus of elasticity and viscous modulus respectively.Therefore, equation (1) can be turned on Lamb wave phase velocity cL、
Frequencies omega, modulus of elasticity μ1With viscous modulus μ2Mathematical modeling, it is as follows.
This formula shows that Lamb wave speed is the function of elasticity, viscosity and Lamb wave frequency.By measure Lamb wave and
The Lamb wave speed c obtained in each frequency is decomposed on frequency domainL(ω), the modulus of elasticity of sample can be estimated by nonlinear fitting
μ1With viscous modulus μ2, this is the physical basis of the present invention.In the specific implementation, corneal it is viscoelastic estimation using equation (2) as
Mathematical modeling, or carry out certain approximate and simplify to equation (2), the empirical equation for obtaining artery viscoelasticity estimation (must be through
Cross the checking of experiment), it is used as the mathematical modeling of estimation.
Theoretical and experimental study shows that the speed of endovascular Lamb wave is general in 3~10m/s, in medical supersonic
On, such spread speed can be arrived with existing hardware detection completely.
Below in conjunction with the accompanying drawings and specific embodiment the present invention is further detailed explanation.
As shown in fig.1, being the hardware architecture diagram of the present invention detection viscoelastic system of arteries.
The system includes main control computer 101, transmit/receive module 102, electrocardiogram acquisition module 103, outside oscillator 104, ultrasound
Probe 105 and electrocardiogram acquisition machine 106.Wherein:Transmit/receive module 102 respectively with main control computer 101, outside oscillator 104, ultrasound
Probe 105 is electrically connected with;Electrocardiogram acquisition module 103 is electrically connected with main control computer 101, electrocardiogram acquisition machine 106 respectively.
The main control computer 101 controls the running of whole system, is responsible for the signal control and data processing of whole system.It is described
Main control computer 101 includes processing module 107, estimation module 108 and registration module 109.
The transmit/receive module 102 includes radiating circuit and receiving circuit, and radiating circuit includes signal generator and work(
Rate amplifier.The electrocardiogram acquisition module 103 is used to trigger electrocardiogram acquisition machine 106.The outside oscillator 104 is used in skin table
Face produces vibration.The ultrasonic probe 105 includes multiple passages, i.e., multiple detections are shaken member, also referred to as test point, in the present embodiment
The ultrasonic probe 105 includes having small spacing between 2 passages, i.e. passage A and passage B, the passage.
The transmit/receive module 102 is used to control outside oscillator 104 producing certain frequency on the skin of blood vessel
Vibration.Specifically:
Outside oscillator 104 is placed on the skin by blood vessel first in the present embodiment, by main control computer 101 triggering transmitting/
The emission control of receiving module 102, control radiating circuit makes outside oscillator produce within 104 every 0.1 seconds the vibration of 1 secondary frequencies, leads
Blood vessel is caused to have a small vibration, the vibration is propagated in the blood vessel in the form of Lamb wave.
The transmit/receive module 102 is additionally operable to driving ultrasonic probe 105 to skin transmitting detection pulse, and record is each
The pulse echo signal of vascular tissue in skin under passage.
The processing module 107 is used for blood in the skin under each passage according to the transmit/receive module 102 record
The pulse echo signal of tubing, calculates the displacement for obtaining vascular tissue in the skin of each passage, and filter out in the displacement
Pulse wave velocity.Specifically:
The echo-signal r (t, k) of the detection pulse received is expressed as a 2D signal.Launch one-time detection
Pulse, then t echo-signal does not represent the echo at different depth in tissue in the same time.In the present embodiment, passage A, B is with one
Fixed pulse recurrence frequency transmitting detection pulse, detects the displacement of vascular tissue, obtains the curve of the displacement versus time, it is former
Reason and the principle of pulse Doppler are similar.
Further, launch repeated detection pulse by certain pulse recurrence frequency, obtain the echo of a detection pulse
Sequence, k represents the frame number of echo in the sequence.The amplitude and phase of echo-signal are modulated by the tissue vibration, are passed through
Certain algorithm can therefrom extract the phase place change of vascular tissue's vibration.When implementing, shake what member was received from detection
Radio frequency (RF) signal passes through the amplification and the processing of analog-to-digital conversion of ultrasonic front-end circuit, then by quadrature demodulation processing, obtains RF
The complex envelope of signal.To a certain test point, the complex envelope signal of echo-signal consecutive frame carries out cross correlation process, extracts the point
The signal of tissue vibration in different depth.
The phase of test point vibration signal in some depth is calculated, the phase of each frame subtracts the phase of consecutive frame
Position, obtains the phase difference of each frame, and then obtains the displacement curve changed over time of vascular tissue.For two inspections at a distance of Z
Measuring point, can ask for phase velocity c by the displacement curve of a certain depth of two test pointsL。
The electrocardiogram acquisition module 103 is used to trigger the collection electrocardiosignal of electrocardiogram acquisition machine 106.
The estimation module 108 is used to filter out after pulse wave velocity blood vessel in the skin of each passage according to obtained above
The displacement of tissue, tries to achieve the phase velocity of the instant exposure of strain wave, estimates the instantaneous viscoelastic coefficient of blood vessel.
In the present embodiment, according to the curve of the displacement versus time of vascular tissue in the above-mentioned skin obtained under passage A, B
Phase velocity of the vibration in 2 points of propagation can be calculated, the viscoelastic coefficient of blood vessel is estimated by mathematical modeling.
The present embodiment is based on mathematical modeling f (c, ω, μ1,μ2)=0 pair coefficient of elasticity μ1With viscosity μ2Carrying out estimation is
The typical Non-linear parameter estimation problem based on model, is solved using Nonlinear least squares fitting method.Simply retouch below
State specific flow:
I. assume<μ1,μ2>One group of initial value;
Ii. the estimation of phase velocity spectrum is obtained based on " criterion of least squares "I.e.:
Iii. the phase velocity spectrum of estimation is calculatedWith actual measurement phase velocity spectrum cphaseError sum of squares between (ω), table
It is shown as:
If iv. α (μ1,μ2) default precision is not up to, updateValue, return the i-th i step;
Otherwise calculate and stop, gainedIt is the estimate of the instantaneous viscoelastic coefficient of blood vessel.
It can be obtained in implementing using optimized algorithms such as iterative method or search methods<μ1,μ2>Optimal estimation.
The registration module 109 was used for the instantaneous viscoelastic coefficient of the blood vessel of estimation and the electrocardiosignal of collection in the time
Upper registration, obtains under cardiac cycle not the modulus of elasticity of blood vessel, viscous modulus in the same time.Specifically:
In the present embodiment, under the cardiac cycle of about one second, repeat above-mentioned mode and encourage vascular vibration and gather back
Ripple signal repeatedly, does registering, so that the blood vessel not in the same time under obtaining cardiac cycle in time with the electrocardiosignal of collection
Modulus of elasticity and viscous modulus.
As shown in fig.2, being the operation process chart of the present invention detection viscoelastic method preferred embodiment of arteries.
Step S400, the outside oscillator 104 of control is producing the vibration of certain frequency on the skin of blood vessel.It is specific and
Speech, main control computer 101 controls outside oscillator 104 to produce the vibration of certain frequency by transmit/receive module 102, acts on outside
Vascular tissue below oscillator 104 makes it produce vibration, and is propagated to horizontal around.
Step S401, the ultrasonic probe 105 of multichannel records under each passage blood in skin to skin transmitting detection pulse
The pulse echo signal of tubing, calculates the displacement for obtaining vascular tissue in the skin of each passage, and filter out in the displacement
Pulse wave velocity, while electrocardiogram acquisition machine 106 start gather electrocardiosignal.Specifically:
The echo-signal r (t, k) of the detection pulse received is expressed as a 2D signal.Launch one-time detection
Pulse, then t echo-signal does not represent the echo at different depth in tissue in the same time.In the present embodiment, passage A, B is with one
Fixed pulse recurrence frequency transmitting detection pulse, detects the displacement of vascular tissue, obtains the curve of the displacement versus time, it is former
Reason and the principle of pulse Doppler are similar.
Further, launch repeated detection pulse by certain pulse recurrence frequency, obtain the echo of a detection pulse
Sequence, k represents the frame number of echo in the sequence.The amplitude and phase of echo-signal are modulated by the tissue vibration, are passed through
Certain algorithm can therefrom extract the phase place change of vascular tissue's vibration.When implementing, shake what member was received from detection
Radio frequency (RF) signal passes through the amplification and the processing of analog-to-digital conversion of ultrasonic front-end circuit, then by quadrature demodulation processing, obtains RF
The complex envelope of signal.To a certain test point, the complex envelope signal of echo-signal consecutive frame carries out cross correlation process, extracts the point
The signal of tissue vibration in different depth.
The phase of test point vibration signal in some depth is calculated, the phase of each frame subtracts the phase of consecutive frame
Position, obtains the phase difference of each frame, and then obtains the displacement curve changed over time of vascular tissue.For two inspections at a distance of Z
Measuring point, can ask for phase velocity c by the displacement curve of a certain depth of two test pointsL。
Step S402, according to the position obtained above for filtering out vascular tissue in the skin of each passage after pulse wave velocity
Move, try to achieve the phase velocity of the instant exposure of strain wave, estimate the instantaneous viscoelastic coefficient of blood vessel.
In the present embodiment, according to the curve of the displacement versus time of vascular tissue in the above-mentioned skin obtained under passage A, B
Phase velocity of the vibration in 2 points of propagation can be calculated, the viscoelastic coefficient of blood vessel is estimated by mathematical modeling.
The present embodiment is based on mathematical modeling f (c, ω, μ1,μ2)=0 pair coefficient of elasticity μ1With viscosity μ2Carrying out estimation is
The typical Non-linear parameter estimation problem based on model, is solved using Nonlinear least squares fitting algorithm.Simply retouch below
State specific flow:
I. assume<μ1,μ2>One group of initial value;
Ii. the estimation of phase velocity spectrum is obtained based on " criterion of least squares "I.e.:
Iii. the phase velocity spectrum of estimation is calculatedWith actual measurement phase velocity spectrum cphaseError sum of squares between (ω), table
It is shown as:
If iv. α (μ1,μ2) default precision is not up to, updateValue, return the i-th i step;
Otherwise calculate and stop, gainedIt is the estimate of the instantaneous viscoelastic coefficient of blood vessel.
It can be obtained in implementing using optimized algorithms such as iterative method or search methods<μ1,μ2>Optimal estimation.
Step S403, the instantaneous viscoelastic coefficient of the blood vessel of estimation and the electrocardiosignal gathered is registering in time, obtain
The not modulus of elasticity of blood vessel, viscous modulus in the same time under to cardiac cycle.Specifically:
In the present embodiment, under the cardiac cycle of about one second, repeat above-mentioned mode and encourage vascular vibration and gather back
Ripple signal repeatedly, does registering, so that the blood vessel not in the same time under obtaining cardiac cycle in time with the electrocardiosignal of collection
Modulus of elasticity and viscous modulus.
The present invention is repeatedly measured under cardiac cycle.Because in the blood vessel of pulsation, blood pressure is in change, blood pressure difference viscoelastic
Property also can be different, such as rubber band, tension can be different with the consistency and elasticity under relaxation state.Therefore, this implementation
Example carries out repeatedly (>=10) measurement under the cardiac cycle of about one second, can obtain under cardiac cycle not in the same time under it is instantaneous
Viscoelastic Spectral, it is significant to diagnosis human vas lesion.The pulse recurrence frequency of outside a period of time excitation in the present embodiment
For 10HZ, i.e., once, vibration every time can all launch a certain amount of detection pulse detection displacement and change over time for vibration in every 0.1 second,
Obtain the viscoelastic coefficient of blood vessel.
The present invention detection viscoelastic system of arteries, can be independent device, can also be used as additional function
Module, is carried in common color ultrasound system.When as additional functional module, it is only necessary to loaded in common color ultrasound system outer
Portion's oscillator 104 and electrocardiogram acquisition machine 106 can be achieved, cost-effective, convenient and simple.
Although the present invention is described with reference to current better embodiment, those skilled in the art should be able to manage
Solution, above-mentioned better embodiment is only used for illustrating the present invention, any in the present invention not for limiting protection scope of the present invention
Spirit and spirit within, any modification, equivalence replacement, improvements for being done etc., should be included in the present invention right guarantor
Within the scope of shield.
Claims (10)
1. one kind detection viscoelastic system of arteries, it is characterised in that the system includes main control computer and the master control is electromechanical
Property connection transmit/receive module and electrocardiogram acquisition module, the outside oscillator being electrically connected with the transmit/receive module
And ultrasonic probe, the electrocardiogram acquisition machine being electrically connected with the electrocardiogram acquisition module, the main control computer include processing module,
Estimation module and registration module, wherein:
The transmit/receive module is used to control outside oscillator producing vibration on the skin of blood vessel under cardiac cycle;
The transmit/receive module is additionally operable to drive ultrasonic probe to skin transmitting detection pulse under cardiac cycle, and record is super
The pulse echo signal of vascular tissue in skin under each passage of sonic probe;
The processing module is used for the pulse echo signal of vascular tissue in the skin under each passage according to above-mentioned record, meter
The displacement for obtaining vascular tissue in the skin of each passage is calculated, and filters out the pulse wave velocity in the displacement;
The electrocardiogram acquisition module is used to trigger electrocardiogram acquisition machine collection electrocardiosignal under cardiac cycle;
The estimation module is used to filter out after pulse wave velocity vascular tissue in the skin of each passage according to obtained above
Displacement, tries to achieve the phase velocity of the instant exposure of strain wave, estimates the instantaneous viscoelastic coefficient of blood vessel;
The registration module be used for by the instantaneous viscoelastic coefficient of the blood vessel of estimation with collection electrocardiosignal it is registering in time,
Obtain under cardiac cycle the not modulus of elasticity of blood vessel, viscous modulus in the same time.
2. the system as claimed in claim 1, it is characterised in that the ultrasonic probe is popped one's head in for channel ultrasonic.
3. system as claimed in claim 2, it is characterised in that described processing module is used to calculate the skin for obtaining each passage
The displacement of Fu Nei vascular tissues, be specially:
From the amplification and the processing of analog-to-digital conversion for detecting first ultrasonic front-end circuit of radiofrequency signal process received that shakes, then by just
Demodulation process is handed over, the complex envelope of radiofrequency signal is obtained;
To a certain test point, cross correlation process is carried out to the complex envelope signal of echo-signal consecutive frame, the different depths of point are extracted
The signal of tissue vibration in degree;
The phase of test point vibration signal in some depth is calculated, the phase of each frame subtracts the phase of consecutive frame, obtained
The phase difference of each frame is obtained, and then obtains the displacement curve that vascular tissue changes over time, for two test points at a distance of Z, is led to
The displacement curve for crossing a certain depth of two test points asks for phase velocity cL。
4. system as claimed in claim 3, it is characterised in that described estimation module be based on Lamb wave mathematical modeling f (c,
ω,μ1,μ2)=0, is specifically included:
I. assume<μ1,μ2>One group of initial value;
Ii. the estimation of phase velocity spectrum is obtained based on " criterion of least squares "I.e.:
Iii. the phase velocity spectrum of estimation is calculatedWith actual measurement phase velocity spectrum cphaseError sum of squares between (ω), is represented
For:
If iv. α (μ1,μ2) default precision is not up to, updateValue, return the i-th i step;
Otherwise calculate and stop, gainedIt is the estimate of the instantaneous viscoelastic coefficient of blood vessel.
5. system as claimed in claim 4, it is characterised in that described registration module specifically for:
The instantaneous viscoelastic coefficient of the blood vessel of estimation is done registering in time with the electrocardiosignal gathered, to obtain cardiac cycle
Under blood vessel elasticity modulus not in the same time and viscous modulus.
6. one kind detection viscoelastic method of arteries, it is characterised in that this method comprises the following steps:
A. outside oscillator is controlled to produce vibration on the skin of blood vessel;
B. the ultrasonic probe of multichannel is to skin transmitting detection pulse, and the pulse for recording vascular tissue in skin under each passage is returned
Ripple signal, and calculate and obtain the displacement of vascular tissue in the skin of each passage, and filter out the pulse wave velocity in the displacement,
Electrocardiogram acquisition machine starts to gather electrocardiosignal simultaneously;
C. according to the displacement obtained above for filtering out vascular tissue in the skin of each passage after pulse wave velocity, strain wave is tried to achieve
Instant exposure phase velocity, estimate blood vessel instantaneous viscoelastic coefficient;
D. it is the instantaneous viscoelastic coefficient of the blood vessel of estimation and the electrocardiosignal gathered is registering in time, obtain under cardiac cycle
The not modulus of elasticity and viscous modulus of blood vessel in the same time.
7. method as claimed in claim 6, it is characterised in that described step b includes:
From the amplification and the processing of analog-to-digital conversion for detecting first ultrasonic front-end circuit of radiofrequency signal process received that shakes, then by just
Demodulation process is handed over, the complex envelope of radiofrequency signal is obtained;
To a certain test point, cross correlation process is carried out to the complex envelope signal of echo-signal consecutive frame, the different depths of point are extracted
The signal of tissue vibration in degree;
The phase of test point vibration signal in some depth is calculated, the phase of each frame subtracts the phase of consecutive frame, obtained
Each frame phase difference is obtained, and then obtains the displacement curve that vascular tissue changes over time, for two test points at a distance of Z, is passed through
The displacement curve of a certain depth of two test points asks for phase velocity cL。
8. method as claimed in claim 7, it is characterised in that this method also includes:
Within a cardiac cycle, step a is repeated to step c.
9. method as claimed in claim 8, it is characterised in that described step c idiographic flows include:
I. assume<μ1,μ2>One group of initial value;
Ii. the estimation of phase velocity spectrum is obtained based on " criterion of least squares "I.e.:
Iii. the phase velocity spectrum of estimation is calculatedWith actual measurement phase velocity spectrum cphaseError sum of squares between (ω), is represented
For:
If iv. α (μ1,μ2) default precision is not up to, updateValue, return the i-th i step;
Otherwise calculate and stop, gainedIt is the estimate of the instantaneous viscoelastic coefficient of blood vessel.
10. method as claimed in claim 9, it is characterised in that described step d is specifically included:
The instantaneous viscoelastic coefficient of the blood vessel of estimation is done registering in time with the electrocardiosignal gathered, it is aroused in interest to obtain one
Blood vessel elasticity modulus not in the same time and viscous modulus in cycle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710153465.4A CN106963424A (en) | 2017-03-15 | 2017-03-15 | Detect the viscoelastic system and method for arteries |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710153465.4A CN106963424A (en) | 2017-03-15 | 2017-03-15 | Detect the viscoelastic system and method for arteries |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106963424A true CN106963424A (en) | 2017-07-21 |
Family
ID=59329547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710153465.4A Pending CN106963424A (en) | 2017-03-15 | 2017-03-15 | Detect the viscoelastic system and method for arteries |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106963424A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109222924A (en) * | 2018-11-05 | 2019-01-18 | 广州丰谱信息技术有限公司 | Cardiovascular property non-invasive quantitative assessment method and device based on acousto-optic combined spectrum |
CN109394200A (en) * | 2018-12-16 | 2019-03-01 | 冯兴怀 | Microcirculation pulse blood flow for liquid undergoing treatment, capacity management monitors system and method |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004290408A (en) * | 2003-03-27 | 2004-10-21 | Aloka Co Ltd | Ultrasonic diagnosing device |
CN1805713A (en) * | 2003-06-13 | 2006-07-19 | 松下电器产业株式会社 | Ultrasonic diagnostic device |
CN1838914A (en) * | 2004-05-21 | 2006-09-27 | 松下电器产业株式会社 | Ultrasonic diagnostic apparatus and method for controlling the same |
CN101039627A (en) * | 2004-10-19 | 2007-09-19 | 松下电器产业株式会社 | Ultrasonographic device and ultrasonographic device control method |
JP2011234863A (en) * | 2010-05-10 | 2011-11-24 | Toshiba Corp | Ultrasonic diagnostic apparatus and image processor |
CN102805650A (en) * | 2011-06-03 | 2012-12-05 | 富士胶片株式会社 | Ultrasound diagnostic apparatus |
CN103006274A (en) * | 2012-12-21 | 2013-04-03 | 深圳大学 | Method and system for ultrasonic detection of cornea viscoelasticity |
CN104055541A (en) * | 2014-06-26 | 2014-09-24 | 中国科学院苏州生物医学工程技术研究所 | Method for intravascular ultrasound multi-slice shear wave elastography |
CN104135944A (en) * | 2012-02-20 | 2014-11-05 | 富士胶片株式会社 | Ultrasound diagnosis device and ultrasound image generation method |
CN105877717A (en) * | 2016-03-30 | 2016-08-24 | 云南大学 | Subcutaneous superficial layer artery blood vessel local elastic distribution detection probe |
-
2017
- 2017-03-15 CN CN201710153465.4A patent/CN106963424A/en active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004290408A (en) * | 2003-03-27 | 2004-10-21 | Aloka Co Ltd | Ultrasonic diagnosing device |
CN1805713A (en) * | 2003-06-13 | 2006-07-19 | 松下电器产业株式会社 | Ultrasonic diagnostic device |
CN1838914A (en) * | 2004-05-21 | 2006-09-27 | 松下电器产业株式会社 | Ultrasonic diagnostic apparatus and method for controlling the same |
CN101039627A (en) * | 2004-10-19 | 2007-09-19 | 松下电器产业株式会社 | Ultrasonographic device and ultrasonographic device control method |
JP2011234863A (en) * | 2010-05-10 | 2011-11-24 | Toshiba Corp | Ultrasonic diagnostic apparatus and image processor |
CN102805650A (en) * | 2011-06-03 | 2012-12-05 | 富士胶片株式会社 | Ultrasound diagnostic apparatus |
CN104135944A (en) * | 2012-02-20 | 2014-11-05 | 富士胶片株式会社 | Ultrasound diagnosis device and ultrasound image generation method |
CN103006274A (en) * | 2012-12-21 | 2013-04-03 | 深圳大学 | Method and system for ultrasonic detection of cornea viscoelasticity |
CN104055541A (en) * | 2014-06-26 | 2014-09-24 | 中国科学院苏州生物医学工程技术研究所 | Method for intravascular ultrasound multi-slice shear wave elastography |
CN105877717A (en) * | 2016-03-30 | 2016-08-24 | 云南大学 | Subcutaneous superficial layer artery blood vessel local elastic distribution detection probe |
Non-Patent Citations (1)
Title |
---|
MATHIEU COUADE: "《QUANTITATIVE ASSESSMENT OF ARTERIALWALL BIOMECHANICAL》", 《ULTRASOUND IN MED. & BIOL.》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109222924A (en) * | 2018-11-05 | 2019-01-18 | 广州丰谱信息技术有限公司 | Cardiovascular property non-invasive quantitative assessment method and device based on acousto-optic combined spectrum |
CN109394200A (en) * | 2018-12-16 | 2019-03-01 | 冯兴怀 | Microcirculation pulse blood flow for liquid undergoing treatment, capacity management monitors system and method |
CN109394200B (en) * | 2018-12-16 | 2021-10-19 | 冯兴怀 | Microcirculation pulse blood flow monitoring system and method for liquid treatment and volume management |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3652791B2 (en) | Ultrasonic diagnostic equipment | |
Beulen et al. | Toward noninvasive blood pressure assessment in arteries by using ultrasound | |
CN100548224C (en) | Utilize ultrasound detection inner narrow with the equipment of identification by its tissue that causes vibration | |
Rabben et al. | An ultrasound-based method for determining pulse wave velocity in superficial arteries | |
EP2675349B1 (en) | A method of and a system for determining a cardiovascular quantity of a mammal | |
Seo et al. | Noninvasive arterial blood pressure waveform monitoring using two-element ultrasound system | |
CN103479343B (en) | Central aortic pressure detection system and method based on oscillating sphygmomanometer signals | |
Sahani et al. | Automated system for imageless evaluation of arterial compliance | |
Tzschätzsch et al. | In vivo time harmonic elastography of the human heart | |
Bambi et al. | A novel ultrasound instrument for investigation of arterial mechanics | |
CN104042200B (en) | A kind of non-invasive monitoring devices and methods therefor of arteriotony by shooting | |
Lillie et al. | Pulse wave velocity prediction and compliance assessment in elastic arterial segments | |
Pereira et al. | Characterization of a double probe for local pulse wave velocity assessment | |
Huang et al. | Continuous measurement of arterial diameter using wearable and flexible ultrasonic sensor | |
CN107106125A (en) | System and method for measuring arterial parameters | |
CN106963424A (en) | Detect the viscoelastic system and method for arteries | |
US7160250B2 (en) | Method and equipment for analyzing biological signals representing intracranial and blood pressure fluctuations | |
Li et al. | High frame rate and high line density ultrasound imaging for local pulse wave velocity estimation using motion matching: A feasibility study on vessel phantoms | |
Nenadic et al. | In vivo open-and closed-chest measurements of left-ventricular myocardial viscoelasticity using lamb wave dispersion ultrasound vibrometry (LDUV): a feasibility study | |
Olesen et al. | Noninvasive estimation of 2-D pressure gradients in steady flow using ultrasound | |
Thorne et al. | In vivo Doppler ultrasound quantification of turbulence intensity using a high-pass frequency filter method | |
McQuilkin et al. | Digital thermal monitoring (DTM) of vascular reactivity closely correlates with Doppler flow velocity | |
Lombardi et al. | Flow Rate Profiler: an instrument to measure blood velocity profiles | |
Bikia | Non-invasive monitoring of key hemodynamical and cardiac parameters using physics-based modelling and artificial intelligence | |
JP2021514281A (en) | A method and apparatus for simultaneously performing 4D ultrafast Doppler imaging of cardiac blood flow and tissue and acquisition of quantification parameters. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170721 |
|
RJ01 | Rejection of invention patent application after publication |