CN104546012B - Cardiac function assessment method and apparatus - Google Patents
Cardiac function assessment method and apparatus Download PDFInfo
- Publication number
- CN104546012B CN104546012B CN201410853915.7A CN201410853915A CN104546012B CN 104546012 B CN104546012 B CN 104546012B CN 201410853915 A CN201410853915 A CN 201410853915A CN 104546012 B CN104546012 B CN 104546012B
- Authority
- CN
- China
- Prior art keywords
- left ventricle
- vortex
- width
- heart left
- heart
- 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.)
- Active
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/06—Measuring blood flow
-
- 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/48—Diagnostic techniques
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Hematology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention provides a kind of Cardiac function assessment method and apparatus, wherein this method includes：Choose the continuous heart left ventricle's ultrasonoscopy of N width；Texture Matching analysis is carried out to the continuous heart left ventricle's ultrasonoscopy of N width, to obtain the displacement vector figure of blood flow in different moments heart left ventricle；According to the displacement vector figure of blood flow in different moments heart left ventricle and the time interval of ultrasound image acquisition, the blood flow velocity vector distribution map in different moments heart left ventricle is calculated；According to the blood flow velocity vector distribution map in different moments heart left ventricle, the vortex parameter in heart left ventricle is calculated, and assess cardiac function according to vortex parameter.The present invention solves of high cost, the technical problem that time-consuming, big to human injury, precision is relatively low of assessment cardiac function in the prior art, has reached effective reduction assessed cost, has reduced the assessment time, reduces the injury to human body and improve the technique effect of assessment result accuracy.
Description
Technical field
The present invention relates to ultrasound in medicine and biology technical field, more particularly to a kind of Cardiac function assessment method and apparatus.
Background technology
Cardiovascular and cerebrovascular disease seriously endangers human health, it has also become the great social common problem being increasingly urgent to.Heart
Whether function normally also becomes people's issues that need special attention.
Currently used for cardiac function measure method include mainly：Left ventriculography, Radionuclide Ventriculography and super
Sound cardiogram.The generally acknowledged left ventriculography of clinic is " goldstandard " of cardiac function evaluation, and echocardiogram is simple and easy to do, accurate because of its
Really reliable and common method as clinical evaluation cardiac function.
However, although left ventriculography in terms of evaluating cardiac function there are good effect, left ventriculography need by
Conduit gos deep into inside ventricle, belongs to invasive interventional therapy, and conduit or seal wire attempt to travel through transposed valve and can cause to wind
Or the displacement of annulus, and for after myocardial infarction or the chambers of the heart expands and suspects have left ventricle thrombosis or the aneurysm to be formed
Patient's operation still has certain danger.
Radionuclide Ventriculography, because of the method that magnetic resonance imaging is utilized, the speed of imaging is slower, and expense is high
It is expensive, heart can not be imaged in real time, radionuclide secondly is injected to human body, certain harm is had to human body.
It is as a result more accurate although the method for current Echocardiography on Left Ventricular Systolic Function in Mitral Regurgitation is more, clinically still
So there are some actual conditions, it is difficult to Left ventricular systolic function is accurately evaluated, such as：Cardiac dysrhythmia, the especially heart
The dynamic patient of atrial fibrillation, especially heart failure merge the patient of auricular fibrillation, clinically require accurate evaluation heart function, heart failure
Exhaust the Evaluation on Cardiac Function and hypertrophic cardiomyopathy Left Ventricular Ejection Fraction for merging bicuspid valve and the serious incompetence of aorta petal
False normal etc., these are all the difficult points of echocardiogram cardiac function evaluation.
For in the prior art, assess cardiac function it is of high cost, time-consuming, precision lower skill big to human injury
Art problem, currently no effective solution has been proposed.
Invention content
An embodiment of the present invention provides a kind of Cardiac function assessment methods, to solve to assess cardiac function in the prior art
Technical problem of high cost, time-consuming, big to human injury, precision is relatively low, this method include：
Choose the continuous heart left ventricle's ultrasonoscopy of N width, wherein N is more than or equal to 3；
Texture Matching analysis is carried out to the continuous heart left ventricle's ultrasonoscopy of the N width, to obtain different moments heart
The displacement vector figure of blood flow in left ventricle；
According to the displacement vector figure of blood flow and the time of ultrasound image acquisition in the different moments heart left ventricle
Interval, is calculated the blood flow velocity vector distribution map in different moments heart left ventricle；
According to the blood flow velocity vector distribution map in the different moments heart left ventricle, it is calculated in heart left ventricle
Vortex parameter；
According to the vortex parameter in the heart left ventricle being calculated, cardiac function is assessed.
In one embodiment, Texture Matching analysis is carried out to the continuous heart left ventricle's ultrasonoscopy of the N width, to obtain
The displacement vector figure of blood flow in different moments heart left ventricle is taken, including：
Each image in the continuous heart left ventricle's ultrasound figure of the N width is divided into multiple analyses according to same way
Window；
Since the first width ultrasonoscopy, following operation is repeated, until obtaining the two-dimension displacement of N-1 width ultrasonoscopys
Vector distribution map：By each analysis window of current ultrasonoscopy analysis window corresponding with the lower piece image of the ultrasonoscopy
Mouthful, two-dimensional cross correlation operation is carried out, the displacement vector of the blood flow representated by all analysis windows of current ultrasonoscopy is obtained
Then figure carries out successive ignition to obtained displacement vector figure：Displacement gradient is first passed through to accurately calculate displacement, recycles two dimension
Standard cross correlation algorithm obtains the spatial resolution higher than preceding an iteration using the analysis window smaller than preceding iteration, and
After displacement each time calculates, algorithm is all rejected by error vector to obtain the Displacement Estimation of higher precision；
By the two-dimension displacement vector distribution map of the N-1 width ultrasonoscopys of acquisition, as blood in different moments heart left ventricle
The dynamic displacement vector figure of liquid stream.
In one embodiment, the vortex parameter includes：The vorticity of vortex, the depth of vortex, vortex lateral position,
The vortex relative intensity of the length of vortex, the width of vortex, the vortex relative intensity of entire ventricular area and vortex region.
In one embodiment, it according to the blood flow velocity vector distribution map in the different moments heart left ventricle, calculates
The vortex parameter in heart left ventricle is obtained, including：
According to the blood flow velocity vector distribution map in the different moments heart left ventricle, it is calculated according to following formula
The vorticity of every bit in N-1 width ultrasonoscopys：
Wherein, ω indicates vorticity, vxIndicate the velocity component in the certain point directions x, vyIndicate the speed point in the certain point directions y
Amount；
According to following formula, Fu is carried out according to the sequence of the vorticity value composition to the N-1 width ultrasonoscopy same position
In leaf change, obtain the zeroth order harmonic wave ω of each position of ultrasonoscopy0With an order harmonics ω1：
According to zeroth order harmonic wave ω0, calculate the location and shape information of vortex；
According to an order harmonics ω1With zeroth order harmonic wave ω0, calculate vortex relative intensity and the vortex region of entire ventricular area
Vortex relative intensity.
In one embodiment, according to zeroth order harmonic wave ω0The location and shape information being vortexed is calculated, including：
According to following formula, to zeroth order harmonic wave ω0It is integrated：
Ω0=∫ω0(x,y)dxdy
According to integral result Ω0, according to following formula, calculate the depth VD of vortex：
According to integral result Ω0, according to following formula, calculate the lateral position VT that is vortexed：
According to integral result Ω0With the depth VD of vortex the length VL of vortex is calculated according to following formula：
According to integral result Ω0Lateral position VT calculates the width VW of vortex according to following formula with being vortexed：
Wherein, HlvIndicate the distance between long axis of left ventricle direction top and bottom, x, y indicate cartesian coordinate, h, s table
Show biology coordinate, wherein h indicates that the direction between the bottom and top of left ventricle, s indicate the bottom and top with left ventricle
Between the vertical direction in direction.
In one embodiment, according to an order harmonics ω1With zeroth order harmonic wave ω0, calculate the vortex phase in entire left ventricle region
To the vortex relative intensity of intensity and vortex region, including：
According to following formula, to zeroth order harmonic wave ω0It is integrated：
Ω0=∫ ω0(x,y)dxdy：
According to integral result and an order harmonics ω1, according to following formula, the vortex for calculating entire left ventricle region is relatively strong
Spend RS：
Wherein, S1 indicates entire heart left ventricle region；
According to integral result and an order harmonics ω1, according to following formula, calculate the vortex relative intensity VRS in vortex region：
Wherein, S2 indicates there is the region being vortexed in heart left ventricle.
In one embodiment, the continuous cardiac ultrasound images of N width are that N width continuously squeezes into encapsulated microbubble
Heart left ventricle's long axis view image.
In one embodiment, the continuous heart left ventricle's ultrasonoscopy of the N width was covered at least three cardiac cycles
Ultrasonic system acquired image.
The embodiment of the present invention additionally provides a kind of Cardiac function assessment device, to solve to assess cardiac function in the prior art
Of high cost, the technical problem that time-consuming, big to human injury, precision is relatively low, which includes：
Image selection unit, for choosing the continuous heart left ventricle's ultrasonoscopy of N width, wherein N is more than or equal to 3；
Displacement vector acquiring unit, for carrying out Texture Matching point to the continuous heart left ventricle's ultrasonoscopy of the N width
Analysis, to obtain the displacement vector figure of blood flow in different moments heart left ventricle；
Velocity vector computing unit, for the displacement vector according to blood flow in the different moments heart left ventricle
The time interval of figure and ultrasound image acquisition, is calculated the blood flow velocity vector distribution map in different moments heart left ventricle；
Vortex parameter calculation unit, for being distributed according to the blood flow velocity vector in the different moments heart left ventricle
Figure, is calculated the vortex parameter in heart left ventricle；
Assessment unit, for according to the vortex parameter in the heart left ventricle being calculated, assessing cardiac function.
In one embodiment, institute's displacement vector acquiring unit, including：
Window divides subelement, is used for each image in the continuous heart left ventricle's ultrasound figure of the N width according to phase
It is divided into multiple analysis windows with mode；
Correlation computations subelement, for since the first width ultrasonoscopy, repeating following operation, until obtaining N-1
The two-dimension displacement vector distribution map of width ultrasonoscopy：It will be under each analysis window of current ultrasonoscopy and the ultrasonoscopy
The corresponding analysis window of piece image carries out two-dimensional cross correlation operation, obtains all analysis window institutes of current ultrasonoscopy
Then the displacement vector figure of the blood flow of representative carries out successive ignition to obtained displacement vector figure：It is next smart to first pass through displacement gradient
It is true to calculate displacement, recycle two-dimentional standard cross correlation algorithm once to be changed than preceding using the analysis window smaller than preceding iteration to obtain
For high spatial resolution, and after the calculating of displacement each time, algorithm is all rejected by error vector to obtain higher precision
Displacement Estimation；
Displacement vector determination subelement, the two-dimension displacement vector distribution map of the N-1 width ultrasonoscopys for that will obtain, as
The displacement vector figure of different moments intracardiac blood flowing.
In embodiments of the present invention, it is proposed that a kind of to calculate the variation of heart left ventricle's inside vortex parameter using ultrasonic
Method, so as to realize the assessment to cardiac function, calculating the vortex parameter in heart left ventricle using ultrasonoscopy has
It is noninvasive, without it is invasive, take short and at low cost advantage, when the function of heart gets muddled, blood flow in heart left ventricle
Vortex parameter will change, and can be realized to myocardial ischemia, cardiac muscle using the variation of heart left ventricle's inside vortex parameter
The heart diseases such as infraction, heart failure are effectively assessed, to solve the cost for assessing cardiac function in the prior art
The technical problem high, time-consuming, big to human injury, precision is relatively low, when having reached effective reduction assessed cost, having reduced assessment
Between, reduce to the injury of human body and improve the technique effect of assessment result accuracy.
Description of the drawings
Attached drawing described herein is used to provide further understanding of the present invention, and is constituted part of this application, not
Constitute limitation of the invention.In the accompanying drawings：
Fig. 1 is the method flow diagram of Cardiac function assessment method according to the ... of the embodiment of the present invention；
Fig. 2 is the another method flow chart of Cardiac function assessment method according to the ... of the embodiment of the present invention；
Fig. 3 is the schematic diagram of the imitative body of left ventricle according to the ... of the embodiment of the present invention；
Fig. 4 is left ventricle integral domain schematic diagram according to the ... of the embodiment of the present invention；
Fig. 5 is another schematic diagram of left ventricle integral domain according to the ... of the embodiment of the present invention；
Fig. 6 is the connection diagram of experimental provision according to the ... of the embodiment of the present invention；
Fig. 7 is the structure diagram of Cardiac function assessment device according to the ... of the embodiment of the present invention.
Specific implementation mode
To make the objectives, technical solutions, and advantages of the present invention clearer, right with reference to embodiment and attached drawing
The present invention is described in further details.Here, the exemplary embodiment and its explanation of the present invention be for explaining the present invention, but simultaneously
It is not as a limitation of the invention.
Inventor considers：Vortex in heart is one and possesses round and moving in rotation fluidal texture, in left ventricle
During filling, the formation of left ventricle inside vortex can effectively preserve the energy of blood flow, and during left ventricular ejection, vortex can
More easily to change the direction of blood flow so that it is more efficient to penetrate blood.The formation of vortex is related with the function of heart, works as heart
Function when getting muddled, the vortex in heart left ventricle can change, so inventor is in view of that can pass through ultrasonic grain
Subgraph velocity measuring technique quantifies whirlpool in heart left ventricle to calculate blood flow velocity field in heart by blood flow velocity field computation
The parameter of rotation has the function that assist diagnosis and treatment to the function of quantitative evaluation heart.
In order to calculate heart left ventricle's inside vortex parameter, may be used：The heart nmr imaging technique of phase-contrast
With the blood flow vector figure based on color Doppler.Although however, the imaging method of nuclear magnetic resonance three dimensions have it is very high
Spatial resolution, but this method is greatly limited face to face in temporal resolution, and very elapsed time, it is colored more
Although general Le method have the advantages that it is cheap, take it is short, with very high temporal resolution, this method is external not yet
Model is verified, and inaccurate to being measured perpendicular to the vortex of sound beam direction.
Therefore, in embodiments of the present invention, ultrasonic Particle Image Velocimetry is selected to measure intracardiac blood dynamics
The available more accurate and higher blood flow flow image of resolving power of parameter, and the advantage that the calculating time is short, at low cost, specifically
, as shown in Figure 1, the Cardiac function assessment method, includes the following steps：
Step 101：Choose the continuous heart left ventricle's ultrasonoscopy of N width, wherein N is more than or equal to 3；
Step 102：Texture Matching analysis is carried out to the continuous heart left ventricle's ultrasonoscopy of the N width, to obtain difference
The displacement vector figure of blood flow in moment heart left ventricle；
Step 103：It is obtained according to the displacement vector figure of blood flow in the different moments heart left ventricle and ultrasonoscopy
The blood flow velocity vector distribution map in different moments heart left ventricle is calculated in the time interval taken；
Step 104：According to the blood flow velocity vector distribution map in the different moments heart left ventricle, heart is calculated
Vortex parameter in left ventricle；
Step 105：According to the vortex parameter in the heart left ventricle being calculated, cardiac function is assessed.
In the above-described embodiments, it is proposed that a kind of side calculating the variation of heart left ventricle's inside vortex parameter using ultrasound
Method, so as to realize the assessment to cardiac function, calculating the vortex parameter in heart left ventricle using ultrasonoscopy has nothing
It creates, without invasive, time-consuming short and at low cost advantage, when the function of heart gets muddled, the whirlpool of blood flow in heart left ventricle
Rotation parameter will change, and can be realized to myocardial ischemia, cardiac muscle stalk using the variation of heart left ventricle's inside vortex parameter
The heart diseases such as plug, heart failure are effectively assessed, to solve in the prior art assessment cardiac function it is of high cost,
The technical problem that time-consuming, big to human injury, precision is relatively low has reached effective reduction assessed cost, has reduced assessment time, drop
The technique effect of low injury and raising assessment result accuracy to human body.
Specifically, in this example, it is proposed that a kind of to measure whirlpool in heart left ventricle using ultrasonic Particle Image Velocimetry
Rotation parameter is come the method for assessing cardiac function, a series of continuous heart left ventricle's ultrasound figures acquired by ultrasonic image-forming system
Picture, using the Texture Matching analysis method in Particle Image Velocimetry, to obtain blood stream in different moments heart left ventricle
Dynamic displacement vector figure, by the time interval of displacement vector figure and sequential image acquisition, (time interval can be obtained by frequency
To), the blood flow velocity distribution map in different moments heart left ventricle is calculated, is calculated by the velocity contour in heart
To the vorticity of heart left ventricle's inside vortex, the size of vortex, the position of vortex, relative intensity, the relative intensity of vortex etc. these
Vortex parameter can realize the assessment to cardiac function by these parameters.
It is illustrated with reference to a specific embodiment, however it is worth noting that, the specific embodiment is merely to more
Illustrate the present invention well, does not constitute improper limitations of the present invention.As shown in Fig. 2, this method includes mainly following step
Suddenly：
Step 1：First, input N width continuously squeezes into heart left ventricle's long axis view image of encapsulated microbubble,
In, long axis of left ventricle tangent plane picture includes：Blood vessel between left ventricle and atrium sinistrum and aorta vessel, wherein N width is continuous
Tangent plane picture should cover ultrasonic system acquired image at least three cardiac cycles.
For example, the Image Acquisition frame frequency FR of ultrasonic system is 100 frames/second, human heartbeat's frequency f is 60Hz/ minutes, aroused in interest
Cycle T c=60/f=1 seconds, then N=m × FR × Tc=100m frames (m=1,2,3 ...), therefore, N can be 100 integer
Times, by N width ultrasonoscopy number consecutivelies, number can be indicated with n.
Step 2：A ROI is selected on first frame image (i.e. piece image), and (Region of Interest, sense are emerging
Interesting region), image sequence n is set as 1.
Step 3：Texture Matching calculating is carried out to the ROI of n-th frame and the (n+1)th frame image, obtains every bit in n-th frame image
Displacement vector, ROI is divided into multiple analysis windows, two-dimensional cross correlation is carried out to corresponding two analysis windows in two width figures
Operation obtains the displacement vector of the blood flow representated by the analysis window, successively to carrying out two dimension mutually per a pair of analysis window in ROI
Related operation obtains the two-dimension displacement vector distribution map in n-th frame image ROI.
Step 4：N=n+1 is enabled, judges whether n is greater than or equal to N, if "No", return to step 3, if "Yes",
Execute step 5.By the loop computation to step 3, the two-dimension displacement vector of the ultrasonoscopy ROI of N-1 width input can be obtained
Distribution map is also equivalent to obtain the displacement vector of each heartbeat moment intracardiac blood flowing in several cardiac cycles
Figure.
Specifically, in step 3 and step 4, correlation analysis is carried out to the continuous multiple frames image collected, it first will figure
As being divided into multiple analysis windows, then to the analysis window of two field pictures using two-dimentional standard cross correlation algorithm, sub-pix
Method and filter coefficients method obtain the two-dimension displacement at this, and then multiple iterative algorithm accurately calculates position by displacement gradient
It moves, last two dimension standard cross correlation algorithm obtains higher spatial resolution using smaller analysis window, and error vector picks
The Displacement Estimation of degree of precision is obtained except algorithm.Then, above-mentioned operation is carried out to all analysis windows in this two field pictures, with
The two-dimension displacement vector field at moment representated by the frame is obtained, is repeated the above process, is obtained different in one or several cardiac cycles
Blood flow displacement vector field in moment left ventricle.
Step 5：The ultrasound of ROI (i.e. left ventricle) interior the two-dimension displacement polar plot and acquisition of the N-1 width obtained by step 4
The blood flow velocity vector distribution map in N-1 width heart left ventricle is calculated in the time interval of image.
Such as：The displacement at a certain moment in left ventricle divided by time interval (are calculated into blood flow in a certain moment left ventricle
The acquisition time interval of two images required for displacement), VELOCITY DISTRIBUTION is obtained (that is, blood flow velocity point in this moment left ventricle
Cloth), the VELOCITY DISTRIBUTION which can be obtained by different moments in left ventricle is carried out to blood flow displacement in different moments left ventricle
.
Step 6：By the blood flow velocity vector distribution map in heart left ventricle, vorticity point in heart left ventricle is calculated
Butut, specifically, vorticity ω can be calculated according to following formula：
In the formula, ω indicates vorticity, vxIndicate the velocity component in the certain point directions x, vyIndicate the certain point directions y
Velocity component.
Step 7：In view of vorticity is the variable of a cycle, the vorticity to N-1 width image same positions can be passed through
The sequence of value composition carries out Fourier's variation, obtains the zeroth order harmonic wave and an order harmonics of each position of image, i.e., in following formula
ω0And ω1Component.
The vorticity of each position is fluctuation in left ventricle, because using discrete Fourier transformation by same position
Vorticity has resolved into the form of different frequency fluctuation adduction, which is, the signal formula of discrete Fourier transform, discrete Fu
In leaf variation waveform can be decomposed into the i.e. different frequency of each order harmonic component trigonometric function component sum.Wherein, ω is indicated
The different amplitude of each order harmonics, T indicate that each order harmonics different period, φ indicate the different phase of each order harmonics.To the N-1
The sequence of the vorticity value composition of width ultrasonoscopy same position carries out discrete Fourier transformation, and it is each to obtain each position vorticity
The component of order harmonics, however, being only concerned about zeroth order harmonic component amplitude ω here0, an order harmonic component amplitude ω1, to basis
Discrete fourier variation represented by above-mentioned formula can obtain the zeroth order harmonic wave ω of each position vorticity left ventricle Nei0With
One order harmonics ω1。
Pass through the zeroth order harmonic wave ω of vorticity0To calculate the location and shape information of heart inside vortex：
First, to zeroth order harmonic wave ω0It is integrated：
Ω0=∫ ω0(x,y)dxdy：
Then, according to integral result Ω0, according to following formula, calculate the depth VD of vortex：
According to integral result Ω0, according to following formula, calculate the lateral position VT that is vortexed：
According to integral result Ω0With the depth VD of vortex the length VL of vortex is calculated according to following formula：
According to integral result Ω0Lateral position VT calculates the width VW of vortex according to following formula with being vortexed：
In above-mentioned formula, it is illustrated in figure 3 the schematic diagram of the imitative body of a left ventricle, x and y indicate cartesian coordinate,
But calculate shape and during position used in be biology coordinate (h and s), the two coordinate representation left ventricles
The both direction of biology coordinate, one is direction (h) between the bottom and top of left ventricle, and one is perpendicular to the direction
Direction (s), HlvIndicate the distance between long axis of left ventricle direction top and bottom, i.e. at the top of the left ventricle of the directions h and bottom in Fig. 3
The distance between portion.
Step 8：Pass through an order harmonics ω1With zeroth order harmonic wave ω0, the vortex for calculating entire heart left ventricle region is relatively strong
The vortex relative intensity of degree and vortex region：
According to integral result and an order harmonics ω1, according to following formula, calculate the vortex phase in entire heart left ventricle region
To intensity RS：
Wherein, S1 indicates entire heart left ventricle region；
According to integral result and an order harmonics ω1, according to following formula, calculate the vortex relative intensity VRS in vortex region：
Wherein, S2 indicates there is the region being vortexed in heart left ventricle.
Although RS with VRS computational methods are identical, the integral domain of the two is different, and RS (relative intensity) is to entire
Left ventricle region is integrated, and as shown in the white stripes region in Fig. 4, the relative intensity of vortex is to there is the region being vortexed
The integral of progress, as shown in the white stripes region in Fig. 5.
After obtaining vortex parameter, so that it may to be assessed cardiac function by vortex parameter, this assessment institute foundation
Principle be exactly：Such as when occurring some patches in left ventricle, the blood flow of the position can be affected in left ventricle, from
And shape and the position of ventricle inside vortex are changed, or when certain diseases occur for heart, the function reduction of left ventricle,
The intensity of the left indoor vortex of the heart can change.
When specific implementation, Texture Matching method may include：(one or several is covered to the continuous multiple frames collected
A cardiac cycle) image progress correlation analysis：First, ultrasonoscopy is divided into multiple analysis windows, then to analysis window
Two-dimension displacement at this is obtained using two-dimentional standard cross correlation algorithm, subpixel method and filter coefficients method, subsequently, is passed through
Multiple iterative algorithm moves closer to correct, accurate displacement, and finally, two-dimentional standard cross correlation algorithm uses smaller analysis window
Mouthful higher spatial resolution is obtained, algorithm is rejected using error vector to obtain the Displacement Estimation of degree of precision, then to whole
The above-mentioned operation of all analysis windows progress is repeated the above process with obtaining the two-dimension displacement vector field of the frame on a frame image,
Obtain the velocity vector that different moments intracardiac blood flows in one or several cardiac cycles.It is obtained using Texture Matching method
In to heart after the distribution of blood flow velocity field, so that it may which, to calculate the vorticity of every bit in image, by vorticity, this parameter can
To calculate the shape and location information of heart inside vortex, and calculate the relative intensity of heart inside vortex.
Experimental provision is connected to carry out experiment in vitro by Experimental equipment shown in fig. 6, and in figure 6, circle is exactly to connect
Enter the position of left ventricle, left ventricle is placed into the sink for filling water, and imitation is connected at the both ends of left ventricle
In the jerk pump of cardiac pumping, surge flask, flowmeter and device for pressure measurement are connected into left ventricle and jerk pump.
Specifically, left ventricle can make in the following way：(the quality hundred of 4%wt is internally injected in mold
Divide than being agarose solution 4%), forms the filler of ventricular model, then place filler into outer mold, pass through mould
The aperture of tool top is injected into the solution of PVA, is put into -20 DEG C of refrigerator that (freeze thawing is at -20 DEG C by freeze thawing twice
12 hours of refrigerator melt 12 hours in room temperature), finally, by squeezing the agarose pulverized in left ventricle, use water
Agarose is gone out, PVA left ventricles are obtained.
The manufacturing process of PVA may include：Use a concentration of 8wt% of PVA (mass percent 8%), Sigmacell
Cellulose imitates body for the hydrogel of 1wt% to make tissue, accurately weighs PVA polymer first and is added in beaker, and it is fixed to be added
The pure water of amount stirs 1~2 hour at 30 DEG C, PVA polymer is made fully to be swollen, and PVA polymerize when in favor of heating below
The dissolving of object, adjustment temperature are observed the particulate matter for not having to suspend in solution, are become substantially equal to 95 DEG C, 1 hour of heating stirring
One it is transparent when, the loss of moisture in order to prevent during heating can cover one layer of tinfoil paper paper membrane on container.
After PVA polymer is completely dissolved, the Sigmacell Cellulose weighed up are added, continue heating stirring 15 minutes, make
Sigmacell Cellulose are uniformly distributed in PVA aqueous solutions, and at this moment solution is rendered as milky, after stopping heating stirring,
Since viscosity is larger, bubble might have in the solution, container is put into ultrasonic cleaners, shake 15 minutes, remove solution
Internal bubble, is then allowed to stand and is cooled to 30 DEG C.
In the above-described embodiments, it is proposed that a kind of to utilize ultrasound quantification heart left ventricle inside vortex parameter evaluation cardiac function
Technology, which can be as the supplementary means of interventional treatment doctor, physician and diagnosing image technician, in conjunction with heart
The variation of the parameters such as position, shape, the intensity of inside vortex further diagnoses the cardiac health of patient, specifically
, mainly there is the application of the following aspects：
1) this method can provide the variation of real-time, non-intruding heart left ventricle's inside vortex parameter, and obtain heart
The distribution map of left ventricle inside vortex, when cardiac function gets muddled, the blood flow vortex in heart left ventricle can change, and lead to
The variation for crossing heart left ventricle's inside vortex parameter can be with the health status of effective evaluation heart；
2) the ultrasound quantification heart left ventricle inside vortex parameter that this example is proposed can be applied in myocardial ischemia, mental and physical efforts stalk
In the assessment of the heart diseases such as plug, heart failure, there is important application value in terms of the assessment of these heart diseases；3) this example
The ultrasound quantification heart left ventricle inside vortex parameter proposed can be also used for judging curative effect of medication and clinical operation treatment
Recruitment evaluation.
The method for evaluating cardiac function using ultrasound that this example is proposed belongs to non-invasive methods, overcomes in the prior art
Left ventriculography needs the problem of injecting conduit to ventricle, certain damage can be caused to human body, using the method for ultrasonic imaging
It overcomes Radionuclide Ventriculography in the prior art to need to inject radionuclide into human body, some radiation can be brought
The problem of effect, while Radionuclide Ventriculography can not detect heart dynamical state in real time, pass through this example
Cardiac function evaluation method can obtain the higher blood flow flow image of accuracy higher, resolution ratio, and the calculating time is short, cost
It is low, cardiac magnetic resonance imaging somewhat expensive is overcome, the problem of time length is calculated.
Based on same inventive concept, a kind of Cardiac function assessment device is additionally provided in the embodiment of the present invention, it is such as following
Described in embodiment.Since the principle that Cardiac function assessment device solves the problems, such as is similar to Cardiac function assessment method, heart
The implementation of functional assessment device may refer to the implementation of Cardiac function assessment method, and overlaps will not be repeated.It is following to be used
, the combination of the software and/or hardware of predetermined function may be implemented in term " unit " or " module ".Although following embodiment institute
The device of description preferably realized with software, but the realization of the combination of hardware or software and hardware be also may and quilt
Conception.Fig. 7 is a kind of structure diagram of the Cardiac function assessment device of the embodiment of the present invention, as shown in fig. 7, comprises：Image
Selection unit 701, displacement vector acquiring unit 702, velocity vector computing unit 703, vortex parameter calculation unit 704 and assessment
Unit 705 below illustrates the structure.
Image selection unit 701, for choosing the continuous heart left ventricle's ultrasonoscopy of N width, wherein N is more than or equal to 3；
Displacement vector acquiring unit 702, for carrying out Texture Matching to the continuous heart left ventricle's ultrasonoscopy of the N width
Analysis, to obtain the displacement vector figure of blood flow in different moments heart left ventricle；
Velocity vector computing unit 703, for being sweared according to the displacement of blood flow in the different moments heart left ventricle
The blood flow velocity vector distribution in different moments heart left ventricle is calculated in the time interval of spirogram and ultrasound image acquisition
Figure；
Vortex parameter calculation unit 704, for according to the blood flow velocity vector in the different moments heart left ventricle point
The vortex parameter in heart left ventricle is calculated in Butut；
Assessment unit 705, for according to the vortex parameter in the heart left ventricle being calculated, commenting cardiac function
Estimate.
In one embodiment, displacement vector acquiring unit 702 includes：Window divides subelement, is used for the N width
Each image in continuous heart left ventricle ultrasound figure is divided into multiple analysis windows according to same way；Correlometer operator list
Member, for since the first width ultrasonoscopy, repeating following operation, until obtaining the two-dimension displacement of N-1 width ultrasonoscopys
Vector distribution map：By each analysis window of current ultrasonoscopy analysis window corresponding with the lower piece image of the ultrasonoscopy
Mouthful, two-dimensional cross correlation operation is carried out, the displacement vector of the blood flow representated by all analysis windows of current ultrasonoscopy is obtained
Then figure carries out successive ignition to obtained displacement vector figure：Displacement gradient is first passed through to accurately calculate displacement, recycles two dimension
Standard cross correlation algorithm obtains the spatial resolution higher than preceding an iteration using the analysis window smaller than preceding iteration, and
After displacement each time calculates, algorithm is all rejected by error vector to obtain the Displacement Estimation of higher precision；Displacement vector determines
Subelement, the two-dimension displacement vector distribution map of the N-1 width ultrasonoscopys for that will obtain, as different moments intracardiac blood stream
Dynamic displacement vector figure.
In another embodiment, a kind of software is additionally provided, the software is for executing above-described embodiment and preferred reality
Apply the technical solution described in mode.
In another embodiment, a kind of storage medium is additionally provided, above-mentioned software is stored in the storage medium, it should
Storage medium includes but not limited to：CD, floppy disk, hard disk, scratch pad memory etc..
It can be seen from the above description that the embodiment of the present invention realizes following technique effect：Propose a kind of utilization
Ultrasound uses come the method for calculating the variation of heart left ventricle's inside vortex parameter so as to realize the assessment to cardiac function
Ultrasonoscopy calculate the vortex parameter in heart left ventricle have the advantages that it is noninvasive, without it is invasive, take it is short and at low cost, it is careful
When dirty function gets muddled, the vortex parameter of blood flow will change in heart left ventricle, utilize whirlpool in heart left ventricle
The variation of rotation parameter, which can be realized, effectively assesses the heart diseases such as myocardial ischemia, myocardial infarction, heart failure, from
And solve in the prior art assessment cardiac function it is of high cost, time-consuming, precision lower technology big to human injury is asked
Topic has reached effective reduction assessed cost, has reduced the assessment time, reduces the injury to human body and improve assessment result accuracy
Technique effect.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with
It is realized with general computing device, they can be concentrated on a single computing device, or be distributed in multiple computing devices
On the network formed, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it
Store and be performed by computing device in the storage device, and in some cases, can be to be held different from sequence herein
The shown or described step of row, either they are fabricated to each integrated circuit modules or will be multiple in them
Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not limited to it is any specific hard
Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made by
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (6)
1. a kind of Cardiac function assessment device, which is characterized in that including：
Image selection unit, for choosing the continuous heart left ventricle's ultrasonoscopy of N width, wherein N is more than or equal to 3；
Displacement vector acquiring unit, for carrying out Texture Matching analysis to the continuous heart left ventricle's ultrasonoscopy of the N width, with
Obtain the displacement vector figure of blood flow in different moments heart left ventricle；
Velocity vector computing unit, for the displacement vector figure according to blood flow in the different moments heart left ventricle, and
The blood flow velocity vector distribution map in different moments heart left ventricle is calculated in the time interval of ultrasound image acquisition；
Vortex parameter calculation unit, for according to the blood flow velocity vector distribution map in the different moments heart left ventricle, meter
Calculation obtains the vortex parameter in heart left ventricle；
Assessment unit, for according to the vortex parameter in the heart left ventricle being calculated, assessing cardiac function；
Wherein, vortex parameter calculation unit is specifically used for according to the blood flow velocity vector in the different moments heart left ventricle point
The vorticity of every bit in N-1 width ultrasonoscopys is calculated according to following formula for Butut：
Wherein, ω indicates vorticity, vxIndicate the velocity component in the certain point directions x, vyIndicate the velocity component in the certain point directions y；
According to following formula, Fourier is carried out according to the sequence of the vorticity value composition to the N-1 width ultrasonoscopy same position
Variation, obtains the zeroth order harmonic wave ω of each position of ultrasonoscopy0With an order harmonics ω1：
According to zeroth order harmonic wave ω0, calculate the location and shape information of vortex；
According to an order harmonics ω1With zeroth order harmonic wave ω0, calculate the whirlpool of the vortex relative intensity and vortex region of entire ventricular area
Revolve relative intensity；
Wherein, according to zeroth order harmonic wave ω0The location and shape information being vortexed is calculated, including：
According to following formula, to zeroth order harmonic wave ω0It is integrated：
Ω0=∫ ω0(x,y)dxdy
According to integral result Ω0, according to following formula, calculate the depth VD of vortex：
According to integral result Ω0, according to following formula, calculate the lateral position VT that is vortexed：
According to integral result Ω0With the depth VD of vortex the length VL of vortex is calculated according to following formula：
According to integral result Ω0Lateral position VT calculates the width VW of vortex according to following formula with being vortexed：
Wherein, HlvIndicate that the distance between long axis of left ventricle direction top and bottom, x, y indicate that cartesian coordinate, h, s indicate life
Object coordinate, wherein h indicates that the direction between the bottom and top of left ventricle, s indicate between the bottom and top of left ventricle
The vertical direction in direction；φ1Indicate the phase of an order harmonics, φkIndicate that the phase of k order harmonics, T indicate that each order harmonics are different
Period, t indicate time-parameters.
2. device as described in claim 1, which is characterized in that institute's displacement vector acquiring unit, including：
Window divides subelement, is used for each image in the continuous heart left ventricle's ultrasonoscopy of the N width according to identical
Model split is multiple analysis windows；
Correlation computations subelement, for since the first width ultrasonoscopy, repeating following operation, until it is super to obtain N-1 width
The two-dimension displacement vector distribution map of acoustic image：By next width of each analysis window of current ultrasonoscopy and the ultrasonoscopy
The corresponding analysis window of image carries out two-dimensional cross correlation operation, obtains representated by all analysis windows of current ultrasonoscopy
Blood flow displacement vector figure, successive ignition then is carried out to obtained displacement vector figure：Displacement gradient is first passed through accurately to count
Displacement is calculated, recycles two-dimentional standard cross correlation algorithm to be obtained using the analysis window smaller than preceding iteration higher than preceding an iteration
Spatial resolution, and after the calculating of displacement each time, algorithm is all rejected by error vector and obtains the displacement of higher precision
Estimation；
Displacement vector determination subelement, the two-dimension displacement vector distribution map of the N-1 width ultrasonoscopys for that will obtain, as difference
The displacement vector figure of moment intracardiac blood flowing.
3. device as described in claim 1, which is characterized in that the vortex parameter includes：The vorticity of vortex, the depth of vortex
Degree, lateral position, the length of vortex, the width of vortex, the vortex relative intensity of entire ventricular area and the vortex region being vortexed
Vortex relative intensity.
4. device as described in claim 1, which is characterized in that according to an order harmonics ω1With zeroth order harmonic wave ω0, calculate entire left
The vortex relative intensity of ventricular area and the vortex relative intensity in vortex region, including：
According to following formula, to zeroth order harmonic wave ω0It is integrated：
Ω0=∫ ω0(x,y)dxdy：
According to integral result and an order harmonics ω1, according to following formula, calculate the vortex relative intensity RS in entire left ventricle region：
Wherein, S1 indicates entire heart left ventricle region；
According to integral result and an order harmonics ω1, according to following formula, calculate the vortex relative intensity VRS in vortex region：
Wherein, S2 indicates there is the region being vortexed in heart left ventricle.
5. device according to any one of claims 1 to 4, which is characterized in that the continuous heart left ventricle's ultrasound of N width
Image is heart left ventricle's long axis view image that N width continuously squeezes into encapsulated microbubble.
6. device according to any one of claims 1 to 4, which is characterized in that the continuous heart left ventricle's ultrasound of N width
Image covers ultrasonic system acquired image at least three cardiac cycles.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410853915.7A CN104546012B (en) | 2014-12-31 | 2014-12-31 | Cardiac function assessment method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410853915.7A CN104546012B (en) | 2014-12-31 | 2014-12-31 | Cardiac function assessment method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104546012A CN104546012A (en) | 2015-04-29 |
| CN104546012B true CN104546012B (en) | 2018-07-24 |
Family
ID=53063842
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410853915.7A Active CN104546012B (en) | 2014-12-31 | 2014-12-31 | Cardiac function assessment method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104546012B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6808824B2 (en) * | 2016-09-23 | 2021-01-06 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | A system for evaluating the outflow tract stenosis of the subject's heart, how the system operates, and a computer program for evaluating the outflow tract stenosis of the subject's heart. |
| AU2019218655B2 (en) * | 2018-02-07 | 2024-05-02 | Cimon Medical AS - Org.Nr.923156445 | Ultrasound blood-flow monitoring |
| CN117152125B (en) * | 2023-10-25 | 2024-04-12 | 广州爱孕记信息科技有限公司 | Heart abnormality analysis method and device based on ultrasonic image |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080269611A1 (en) * | 2007-04-24 | 2008-10-30 | Gianni Pedrizzetti | Flow characteristic imaging in medical diagnostic ultrasound |
| CN103006272B (en) * | 2013-01-09 | 2014-12-03 | 东南大学 | Velocity distribution measuring method based on ultrasonic interleave programming |
| CN103340620B (en) * | 2013-05-31 | 2016-03-30 | 中国科学院深圳先进技术研究院 | A kind of measuring method of tube wall stress phase angle and system |
| CN104207803B (en) * | 2014-08-13 | 2016-03-02 | 四川省人民医院 | Based on the cardiac flow whirlpool Motion Adaptive visable positioning method of color doppler image information |
-
2014
- 2014-12-31 CN CN201410853915.7A patent/CN104546012B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104546012A (en) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sengupta et al. | Emerging trends in CV flow visualization | |
| Pedrizzetti et al. | On estimating intraventricular hemodynamic forces from endocardial dynamics: a comparative study with 4D flow MRI | |
| Hirsch et al. | In vivo measurement of volumetric strain in the human brain induced by arterial pulsation and harmonic waves | |
| Ebbers et al. | Estimation of relative cardiovascular pressures using time‐resolved three‐dimensional phase contrast MRI | |
| Lotz et al. | Cardiovascular flow measurement with phase-contrast MR imaging: basic facts and implementation | |
| Kefayati et al. | Turbulence intensity measurements using particle image velocimetry in diseased carotid artery models: Effect of stenosis severity, plaque eccentricity, and ulceration | |
| Yang et al. | Flow and myocardial interaction: an imaging perspective | |
| US10383600B2 (en) | Method for analysis and display of blood flow information | |
| US10134127B2 (en) | Method for post-processing flow-sensitive phase contrast magnetic resonance images | |
| CN104546012B (en) | Cardiac function assessment method and apparatus | |
| Soudah et al. | Estimation of wall shear stress using 4D flow cardiovascular MRI and computational fluid dynamics | |
| Kang et al. | High-frequency dual mode pulsed wave Doppler imaging for monitoring the functional regeneration of adult zebrafish hearts | |
| US10776998B1 (en) | Method and system for analysis of 3D deformations and regional function of a heart with 3D SinMod | |
| Tang et al. | The effect of respiration-driven flow waveforms on hemodynamic metrics used in Fontan surgical planning | |
| Riemer et al. | Contrast agent-free assessment of blood flow and wall shear stress in the rabbit aorta using ultrasound image velocimetry | |
| Eriksson et al. | Spatial heterogeneity of four‐dimensional relative pressure fields in the human left ventricle | |
| Yang et al. | Lagrangian-averaged vorticity deviation of spiraling blood flow in the heart during isovolumic contraction and ejection phases | |
| Cygan et al. | Left ventricle phantom and experimental setup for MRI and echocardiography–Preliminary results of data acquisitions | |
| Bourquin et al. | Quantitative pulsatility measurements using 3D dynamic ultrasound localization microscopy | |
| CN103340620B (en) | A kind of measuring method of tube wall stress phase angle and system | |
| Wick et al. | Detection of cardiac quiescence from B-mode echocardiography using a correlation-based frame-to-frame deviation measure | |
| Okrasinski et al. | Assessment of myocardial elastography performance in phantoms under combined physiologic motion configurations with preliminary in vivo feasibility | |
| Kefayati et al. | 3-D flow characterization and shear stress in a stenosed carotid artery bifurcation model using stereoscopic PIV technique | |
| EP3273863A1 (en) | Non-invasive estimation of intravascular pressure changes using vector velocity ultrasound (us) | |
| Yamada et al. | Numerical analysis of the effect of trabeculae carneae models on blood flow in a left ventricle model constructed from magnetic resonance images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |