CN1302751C - Method for detecting myocardial blood supply state based on myocardial edge tracking - Google Patents
Method for detecting myocardial blood supply state based on myocardial edge tracking Download PDFInfo
- Publication number
- CN1302751C CN1302751C CNB2004100627654A CN200410062765A CN1302751C CN 1302751 C CN1302751 C CN 1302751C CN B2004100627654 A CNB2004100627654 A CN B2004100627654A CN 200410062765 A CN200410062765 A CN 200410062765A CN 1302751 C CN1302751 C CN 1302751C
- Authority
- CN
- China
- Prior art keywords
- myocardium
- frame
- particle
- edge
- interest
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 230000036770 blood supply Effects 0.000 title claims abstract description 24
- 230000002107 myocardial effect Effects 0.000 title claims abstract description 21
- 210000004165 myocardium Anatomy 0.000 claims abstract description 131
- 239000002245 particle Substances 0.000 claims abstract description 59
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000010586 diagram Methods 0.000 claims abstract description 4
- 210000005240 left ventricle Anatomy 0.000 claims description 15
- 238000007789 sealing Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 210000001174 endocardium Anatomy 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 4
- 238000000205 computational method Methods 0.000 claims 1
- 230000000302 ischemic effect Effects 0.000 abstract description 4
- 239000003086 colorant Substances 0.000 abstract 1
- 238000003384 imaging method Methods 0.000 description 10
- 208000031225 myocardial ischemia Diseases 0.000 description 8
- 241000270295 Serpentes Species 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 208000029078 coronary artery disease Diseases 0.000 description 3
- 239000012467 final product Substances 0.000 description 3
- 208000028867 ischemia Diseases 0.000 description 3
- 208000010125 myocardial infarction Diseases 0.000 description 3
- 210000003516 pericardium Anatomy 0.000 description 3
- 238000005267 amalgamation Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 208000024891 symptom Diseases 0.000 description 2
- 230000002861 ventricular Effects 0.000 description 2
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 206010025482 malaise Diseases 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention relates to a method for detecting myocardial blood supply states based on myocardial edge tracks, which belongs to the technical field of medical ultrasonics. The method comprises the steps that a mapping table for integral backscattering fluctuating degrees and colors is firstly established; the detected heart is scanned and ultrasonic radio frequency signals are continuously collected to reconstruct a B-mode image sequence reflecting the structure of a scanned area; a focused rectangular area is selected from the first frame and the frame positions of the myocardial edge are tracked; the frame positions of all myocardial particle points in the focused area are calculated; the integral backscattering curve and the integral backscatter fluctuating degree of each myocardial particle point is calculated and the blood supply states diagram of each myocardial particle point in the focused area is obtained according to the mapping table. The method provided by the present invention can effectively distinguish a normal myocardium from an ischemic myocardium. In addition, the ischemic myocardium is imaged for CVIB by means of the method for detecting myocardial blood supply states provided by the present invention and the method of the existing patent and erroneous judgements can be effectively avoided.
Description
Technical field
The present invention relates to a kind of method of measuring myocardial blood supply state, relate in particular to the method for myocardial back scattering integral fluctuation degree (hereinafter to be referred as CVIB) being carried out measuring after the imaging myocardial blood supply state, belong to the medical ultrasonic technical field, the ultrasonic tissue that is used in particular for cardiac muscle is levied the i.e. detection of myocardial ischemia surely.
Background technology
Coronary heart disease and the myocardial infarction sickness rate in the crowd is very high, in mid-aged population, is the main cause of disease that causes death of current society especially.Find that through clinical research most of coronary heart disease or myocardial infarction patient at the ill initial stage, all have the myocardial ischemia symptom.Therefore,, find out the position of myocardial ischemia, just may take the treatment measure more in time, provide technical help the patient of prevention and treatment coronary heart disease or myocardial infarction if can in time find their myocardial ischemia symptom.
Lot of experiments shows that the ultrasonic backscattering integration (hereinafter to be referred as IB) of cardiac muscle can become the new argument of a quantitative Diagnosis myocardial ischemia.Can observe in experiment, IB periodic the fluctuating occur with the cycle aroused in interest, and the amplitude that generally claims this periodic undulations is CVIB.In general, the CVIB value of normal myocardium is bigger, and the CVIB value of ischemic myocardium is less.If can access cardiac muscle at the IB of whole cardiac cycle curve, and calculate the CVIB value, just can a basic understanding be arranged to myocardial blood supply state.
Chinese patent ZL98203234.X discloses a kind of two-dimentional CVIB formation method, and its principle is: with color to different CVIB values encode (be gradient to red corresponding CVIB value from blueness and change to 12dB) from 3dB such as using; On B ultrasonic figure, choose region-of-interest (ROI), and calculate the CVIB value of each myocardium particle in a cardiac cycle in the region-of-interest; Choosing the frame B ultrasonic image that the corresponding heart relaxes latter stage then, is that the imaging factor is undertaken painted by predefined coding to the myocardium particle in the ROI with the CVIB value corresponding with it.Color showing myocardium blood supply situation, since it be with the B ultrasonic image co-registration represented with gray scale together, what B ultrasonic figure comprised is the structural information of heart, the blood supply information of heart and structural information organically merge in a figure like this, can show relative position and the blood supply function situation thereof of region-of-interest in heart very intuitively.
Said method should be followed the trail of respectively each myocardium particle in the region-of-interest theoretically, calculates separately IB curve and CVIB, last imaging.Yet, owing to follow the trail of relatively difficulty of myocardium particle, so existing patented technology is supposed the motion of cardiac muscle and can be ignored when imaging, when calculating the IB curve of certain myocardium particle A in the region-of-interest, it is the pairing radiofrequency signal of point of finding out each frame B ultrasonic picture position and the coincidence of A point, form the IB curve, calculate CVIB, such account form can produce erroneous judgement in some cases.
Summary of the invention
The objective of the invention is to propose a kind of method of measuring myocardial blood supply state, by the myocardium edge in the region-of-interest is followed the trail of automatically, on this basis the position of all myocardium particles in the region-of-interest on each frame B ultrasonic image followed the trail of, to realize two-dimentional CVIB imaging, improve the accuracy of CVIB imaging, the situation of avoiding traditional C VIB formation method to occur judging by accident.
The method based on the myocardium marginal test myocardial blood supply state of tracking that the present invention proposes may further comprise the steps:
(1) sets up the mapping table of Integral Backscatter fluctuating quantity and color;
(2) scanning measurand heart, the ultrasonic radiofrequency signal of continuous acquisition;
(3), rebuild the B ultrasonic image sequence of reflection scanning area structure according to above-mentioned ultrasonic radiofrequency signal;
(4) in first frame of above-mentioned B ultrasonic image sequence, choose the rectangle region-of-interest;
(5) the myocardium edge in the above-mentioned region-of-interest is followed the trail of in the position of each frame;
(6) according to the position of the above-mentioned myocardium edge that tracks in each frame B ultrasonic image, all myocardium particles in the region-of-interest are calculated in each frame position;
(7) according to each myocardium particle position in each two field picture in the above-mentioned region-of-interest, find out each myocardium particle in the corresponding constantly radiofrequency signal section of difference, calculate the Integral Backscatter curve of each myocardium particle according to the radiofrequency signal section;
(8) according to the Integral Backscatter curve of each myocardium particle in the region-of-interest, calculate the Integral Backscatter fluctuating quantity of each myocardium particle, according to above-mentioned mapping table, obtain with region-of-interest in the corresponding color image of blood supply situation of each myocardium particle;
(9) above-mentioned color image is merged mutually with first frame of above-mentioned B ultrasonic image sequence, obtain the blood supply condition diagram of each myocardium particle in the region-of-interest.
Follow the trail of the method for the position of myocardium edge on all the other each frame B ultrasonic images in the said method, its process may further comprise the steps:
(1) the rectangle region-of-interest of choosing on the B ultrasonic image sequence is separated from original image, obtain one group of subimage;
(2) delineate out the initial profile line at myocardium edge from first frame of above-mentioned subimage;
(3) subimage and initial profile line are carried out coordinate transform, image upset and merge the myocardium edge and the initial profile line of formation sealing;
(4) use the active contour model method, the myocardium edge of corresponding sealing is carried out rim detection according to the initial profile line of above-mentioned sealing;
(5) to above-mentioned detected edge, the inverse transformation of carrying out coordinate transform, image upset and merging obtains the myocardium edge in the original subimage;
(6) with the initial profile line of this detected myocardium edge as the subimage of next frame, repeating step (3) obtains the myocardium edge in the region-of-interest on each frame B ultrasonic image to (5).
In the said method, according to the position of myocardium edge in each frame B ultrasonic image, all myocardium particles in the region-of-interest are carried out Calculation Method in each frame position, when the scanning view was left-ventricular short-axis, its process may further comprise the steps:
(1) according to the myocardium edge that tracks, match obtains the center of left ventricle at each frame B ultrasonic image;
(2) calculate in the first frame B ultrasonic image myocardium particle to be followed the trail of and the line at left ventricle center and the angle theta of horizontal direction;
(3) calculate in each frame with the left ventricle center be end points with the horizontal direction angle be ray and endocardium and the epicardial intersection point of θ;
(4) establish cardiac muscle radially for evenly shrinking, calculate and wait to follow the trail of the position of myocardium particle on each frame B ultrasonic image.
When the scanning view was long axis of left ventricle, its process may further comprise the steps:
(1) according to the myocardium edge that tracks, match obtain left ventricle at the center of each frame B ultrasonic image, the fan-shaped angle and the position, segmental two border of myocardium edge and left ventricle center formation;
(2) calculate in the first frame B ultrasonic image the myocardium particle to be followed the trail of and the angle α of line between the left ventricle center and fan-shaped left margin
1
(3) according to above-mentioned angle α
1, calculate wait to follow the trail of myocardium particle and left ventricle center line in all the other each frames with the angle α of fan-shaped left margin
i, i represents the frame ordinal number of B ultrasonic image sequence;
(4) calculating in each frame, is that angle end points and fan-shaped left margin is α with the center of left ventricle
iRay and endocardium and epicardial intersection point;
(5) establish cardiac muscle radially for evenly shrinking, calculate and wait to follow the trail of the position of myocardium particle on each frame B ultrasonic image.
Utilize the method for the mensuration myocardial blood supply state that the present invention proposes, respectively to the normal person, the myocardial ischemia patient carries out the CVIB imaging according to the said CVIB formation method of the present invention.Imaging results shows, for the normal person, the cardiac muscle of choosing any position is as region-of-interest, and the color of cardiac muscle mark mainly all is red and green in the region-of-interest, shows that cardiac muscle is normal, and to the patient of myocardial ischemia, imaging results shows that the myocardium color in the ischemia position has been marked blueness, shows ischemia, other position shows normal based on red and green.This illustrates that method proposed by the invention can distinguish normal myocardium and ischemic myocardium effectively.In addition, put forward the method for measuring myocardial blood supply state and the method for existing patent is carried out the CVIB imaging to the cardiac muscle of ischemia object respectively with the present invention, the result shows have patent report method may cause situation about judging by accident to oneself, can effectively avoid erroneous judgement with institute of the present invention extracting method.
Description of drawings:
Fig. 1 is a reconstructed B-mode image in the inventive method, the region-of-interest that the rectangle frame representative is selected, and wherein Fig. 1 (a) is a long axis view, Fig. 1 (b) short axis view.
The sketch map of Fig. 2 for non-closed boundary being detected with the active contour model method.Wherein (a) is the part of original image, and rectangle frame is a region-of-interest, and curve is an initial profile; (b) comprise the honest rectangular image that subimage forms for region-of-interest after rotating; (c) for being transformed into polar result from figure (b) rectangular coordinate; (d) be to scheme (c) upset, merging sketch map; (e) be amalgamation result; (f) with the Snakes method figure (e) being followed the trail of the border (g)-(i) that obtains is the inverse transformation of (c) to (e); (i) visceral pericardium of curved section for tracking at last among the figure.
Fig. 3 is in the left-ventricular short-axis view, utilize myocardium edge in each frame B ultrasonic picture position to the region-of-interest heart myoplasm point sketch map that follow the trail of the position in each frame B ultrasonic image.
Fig. 4 is in the long axis of left ventricle view, utilize myocardium edge in each frame B ultrasonic picture position to the region-of-interest heart myoplasm point sketch map that follow the trail of the position in each frame B ultrasonic image.
The specific embodiment
The method based on the myocardium marginal test myocardial blood supply state of tracking that the present invention proposes may further comprise the steps:
(1) sets up the mapping table of Integral Backscatter fluctuating quantity and color.Be gradient to the green red next corresponding CVIB value variation from low to high that is gradient to again from blueness.The promptly blue myocardial ischemia of representing, redness represents myocardial blood supply state normal.
(2) scanning measurand heart, the ultrasonic radiofrequency signal of continuous acquisition.In one embodiment of the invention, probe with B ultrasonic equipment commonly used scans measurand, radiofrequency signal from B ultrasonic device interior circuit extraction wave beam after synthetic, through amplifying, gather with high-speed data acquisition card, generally require acquisition time to surpass 1 second, sample rate 10M, what one embodiment of the invention was used is the PCI-9812 data collecting card of AD-LINK company.
(3), rebuild the B ultrasonic image sequence of reflection scanning area structure according to above-mentioned ultrasonic radiofrequency signal.The method of rebuilding the B ultrasonic image sequence is: after radiofrequency signal is removed direct current, do the HILBERT conversion, according to the DSC algorithm in the B ultrasonic system, can form last image then.
(4) choose region-of-interest in first frame of above-mentioned B ultrasonic image sequence, generally to short axis view, choose region-of-interest and should comprise whole left ventricular wall, to long axis view, region-of-interest is chosen left ventricular posterior wall.See the rectangle frame of Fig. 1 (a) and Fig. 1 (b) about the region-of-interest chosen position.
(5) the myocardium edge in the above-mentioned region-of-interest is followed the trail of in the position of each frame.Myocardium edge in the region-of-interest used based on the method for Snakes follow the trail of, myocardium edge to the short axis view sealing, carry out rim detection, follow the trail of and to get final product with traditional Snakes method, to long axis view, the cardiac muscle edge is non-sealing, can not directly use the Snakes method, so follow the trail of with the Snakes method at non-closed boundary of the present invention's proposition.
(6), all myocardium particles in the institute in the region-of-interest are calculated in each frame position according to the position of the above-mentioned myocardium edge that tracks in each frame B ultrasonic image.
(7) according to each myocardium particle position in each two field picture in the above-mentioned region-of-interest, find out each myocardium particle in the corresponding constantly radiofrequency signal section of difference, calculate the Integral Backscatter curve of each myocardium particle according to the radiofrequency signal section.The formula that calculates IB is:
In the formula (1), x (t) is meant radiofrequency signal.
(8) according to the Integral Backscatter curve of each myocardium particle in the region-of-interest, calculate the Integral Backscatter fluctuating quantity of each myocardium particle, according to above-mentioned mapping table, obtain with region-of-interest in the corresponding color image of blood supply situation of each myocardium particle.Calculate any method that CVIB can take open source literature, present embodiment uses the single order Fourier series method, and its process is:
CVIB＝4*max(|X(k)|)/N (2)
Wherein, N is the length of IB sequence, and X (k) is that length is the CZT conversion that the IB sequence of N is ordered at frequency domain 10*N.
(9) above-mentioned color image is merged mutually with first frame of above-mentioned B ultrasonic image sequence, obtain the blood supply condition diagram of each myocardium particle in the region-of-interest.Amalgamation mode is: the point in that first frame of above-mentioned B ultrasonic image sequence is marked, carry out painted getting final product with the color that its CVIB value is corresponding.
Follow the trail of the method for the position of nonocclusive myocardium edge on all the other each frame B ultrasonic images in the said method, its process may further comprise the steps:
(1) the rectangle region-of-interest of choosing on the B ultrasonic image sequence is separated from original image, obtain one group of subimage.Soon the part in the red rectangle frame is separated separately among Fig. 2 (a), if the non-honest rectangular area of selecting can form honest rectangle through rotation, if shown in Fig. 2 (b).
(2) delineate out the initial profile line at myocardium edge from first frame of above-mentioned subimage.Generally can utilize conventional computer to paint software and carry out manual cropping.
(3) subimage and initial profile line are carried out coordinate transform, image upset and merge the myocardium edge and the initial profile line of formation sealing.At first the rectangle region-of-interest is implemented from rectangular coordinate to polar conversion, the angular range of image is 90 degree (being similar to the DSC process in the fan sweeping B ultrasonic) in polar coordinate.In order to change conveniently, at first region-of-interest is rotated, make the rectangle of inclination become vertical rectangle (Fig. 2 (b)), rectangular coordinate is seen Fig. 2 (c) to polar transformation result.About the figure of 2 (c) carried out respectively, about, 180 degree upsets, obtain four width of cloth figure (seeing Fig. 2 (d)), four width of cloth figure are merged, obtain figure conversion end product (Fig. 2 (e)).For the initialization profile, we have also carried out same operation, the results are shown in Figure the dotted line among 2 (e).Through this step, the border of sealing and the initial profile line of sealing have just been arranged on the image that we obtain at last.
(4) use the active contour model method, the myocardium edge of corresponding sealing is carried out rim detection, the results are shown in Figure 2 (f) according to the initial profile line of above-mentioned sealing.
(5) to above-mentioned detected edge, the inverse transformation of carrying out coordinate transform, image upset and merging obtains the myocardium edge in the original subimage.This process and just in time opposite to the process of Fig. 2 (e) by Fig. 2 (a), each goes on foot the result of conversion referring to Fig. 2 (g), Fig. 2 (h), epicardial position is shown in Fig. 2 (i) in Fig. 2 (i), the region-of-interest that finally tracks.
(6) with the initial profile line of this detected myocardium edge as the subimage of next frame, repeating step (3) obtains the myocardium edge in the region-of-interest on each frame B ultrasonic image to (5).
Can use above method to nonocclusive endocardium of left heart major axis and visceral pericardium, carry out endocardium and visceral pericardium are followed the trail of in the position of different frame.
In the said method, according to the position of myocardium edge in each frame B ultrasonic image, all myocardium particles in the region-of-interest are carried out Calculation Method in each frame position, when the scanning view was left-ventricular short-axis, referring to Fig. 3, its process may further comprise the steps:
(1) according to the myocardium edge that tracks, match obtains left ventricle at the center of each frame B ultrasonic image (Ox[i], Oy[i]).The match mode is: the point of finding out with the distance variance minimum of myocardium edge each point promptly is a central point.(Ox[1], Oy[1]) among Fig. 3 and (Ox[i], Oy[i]) are respectively the left ventricle center that the 1st frame and the match of i frame obtain.
(2) calculate in the first frame B ultrasonic image myocardium particle to be followed the trail of and the line at left ventricle center and the angle theta of horizontal direction., as shown in Figure 3, computing formula is:
θ＝atan2(y[1]-Oy[1]，x[1]-Ox[1]) (1)
In the formula (1), x[1], y[1] be meant and wait to follow the trail of the coordinate of myocardium particle at first frame, Ox[1], Oy[1] be the centre coordinate of left ventricle at first frame, atan2 is the MATLAB function, what return is the angle of 0 ~ 2 π.
(3) calculate in each frame with the left ventricle center be end points with the horizontal direction angle be the ray of θ and endocardium and epicardial intersection point (xi[i], yi[i]) and (xo[i], yo[i]., i is the frame ordinal number of B ultrasonic image sequence herein.We suppose that myocardium integral body is less than rotation, so myocardium particle to be followed the trail of in all the other each frames and the line of left ventricle and the angle of horizontal direction also be θ, therefore, wait to follow the trail of myocardium particle should (xi[i], yi[i]) and (xo[i], yo[i]) line on.Referring to Fig. 3 (b).
(4) establish cardiac muscle radially for evenly shrinking, calculate and wait to follow the trail of the position of myocardium particle on each frame B ultrasonic image.Computing formula is:
x[i]＝xi[i]+(xo[i]-xi[i])·(x[1]-xi[1])/(xo[1]-xi[1]) (2)
y[i]＝yi[i]+(yo[i]-y[i])·(y[1]-yi[1])/(yo[1]-yi[1]) (3)
More than be example with myocardium particle to be followed the trail of at first frame coordinate (x[1], y[1]), described and utilized its process of myocardium edge calculations in the short axis view in all the other frame position.All the other myocardium particles in the region-of-interest are implemented same process, just finished the tracking of all myocardium particles in the region-of-interest.
When the scanning view was long axis of left ventricle, its process may further comprise the steps:
(1) according to the myocardium edge that tracks, match obtain left ventricle at the center of each frame B ultrasonic image, the fan-shaped angle θ [i] (i represents the frame ordinal number of B ultrasonic image sequence) and the position, segmental two border of myocardium edge and left ventricle center formation.Fig. 4 (a) and Fig. 4 (b) have shown the 1st frame and i frame left ventricle center, and it is fan-shaped that myocardium edge and left ventricle center constitute.
(2) calculate in the first frame B ultrasonic image the myocardium particle to be followed the trail of and the angle α of line between the left ventricle center and fan-shaped left margin
1
(3) according to above-mentioned angle α
1, calculate wait to follow the trail of myocardium particle and left ventricle center line in all the other each frames with the angle α of fan-shaped left margin
i, i represents the frame ordinal number of B ultrasonic image sequence.Suppose that it is uniformly that segmental radian changes, then:
α[i]＝α[1]*θ[1]/θ[i] (4)
(4) calculating in each frame, is that angle end points and fan-shaped left margin is α with the center of left ventricle
iRay and endocardium and epicardial intersection point (xi[i], yi[i]), (xo[i], yo[i]).
(5) establish cardiac muscle radially for evenly shrinking, calculate and wait to follow the trail of the position of myocardium particle on each frame B ultrasonic image.Computing formula gets final product according to formula (2) and formula (3).
More than be example with myocardium particle to be followed the trail of at first frame coordinate (x[1], y[1]), described and utilized its process of myocardium edge calculations in the long axis view in all the other frame position.All the other myocardium particles in the region-of-interest are implemented same process, just finished the tracking of all myocardium particles in the region-of-interest.
Claims (1)
1, a kind of based on the method for following the trail of myocardium marginal test myocardial blood supply state, it is characterized in that this method may further comprise the steps:
(1) sets up the mapping table of Integral Backscatter fluctuating quantity and color;
(2) scanning measurand heart, the ultrasonic radiofrequency signal of continuous acquisition;
(3), rebuild the B ultrasonic image sequence of reflection scanning area structure according to above-mentioned ultrasonic radiofrequency signal;
(4) in first frame of above-mentioned B ultrasonic image sequence, choose the rectangle region-of-interest;
(5) the myocardium edge in the above-mentioned region-of-interest is followed the trail of in the position of each frame, the method for tracking may further comprise the steps:
(a) the rectangle region-of-interest of choosing on the B ultrasonic image sequence is separated from original image, obtain one group of subimage;
(b) delineate out the initial profile line at myocardium edge from first frame of above-mentioned subimage;
(c) subimage and initial profile line are carried out coordinate transform, image upset and merge the myocardium edge and the initial profile line of formation sealing;
(d) use the active contour model method, the myocardium edge of corresponding sealing is carried out rim detection according to the initial profile line of above-mentioned sealing;
(e) to above-mentioned detected edge, the inverse transformation of carrying out coordinate transform, image upset and merging obtains the myocardium edge in the original subimage;
(f) with the initial profile line of this detected myocardium edge as the subimage of next frame, repeating step (c) obtains the myocardium edge in the region-of-interest on each frame B ultrasonic image to (e);
(6) according to the position of the above-mentioned myocardium edge that tracks in each frame B ultrasonic image, all myocardium particles in the region-of-interest are calculated in each frame position, its computational methods are;
When the scanning view was left-ventricular short-axis, its computational process may further comprise the steps:
(a) according to the myocardium edge that tracks, match obtains the center of left ventricle at each frame B ultrasonic image;
(b) calculate in the first frame B ultrasonic image myocardium particle to be followed the trail of and the line at left ventricle center and the angle theta of horizontal direction;
(c) calculate in each frame with the left ventricle center be end points with the horizontal direction angle be ray and endocardium and the epicardial intersection point of θ;
(d) establish cardiac muscle radially for evenly shrinking, calculate and wait to follow the trail of the position of myocardium particle on each frame B ultrasonic image;
When the scanning view was long axis of left ventricle, its computational process may further comprise the steps:
(a) according to the myocardium edge that tracks, match obtain left ventricle at the center of each frame B ultrasonic image, the fan-shaped angle and the position, segmental two border of myocardium edge and left ventricle center formation;
(b) calculate in the first frame B ultrasonic image the myocardium particle to be followed the trail of and the angle α of line between the left ventricle center and fan-shaped left margin
1
(c) according to above-mentioned angle α
1, calculate wait to follow the trail of myocardium particle and left ventricle center line in all the other each frames with the angle α of fan-shaped left margin
i, i represents the frame ordinal number of B ultrasonic image sequence;
(d) calculating in each frame, is that angle end points and fan-shaped left margin is α with the center of left ventricle
iRay and endocardium and epicardial intersection point;
(e) establish cardiac muscle radially for evenly shrinking, calculate and wait to follow the trail of the position of myocardium particle on each frame B ultrasonic image;
(7) according to each myocardium particle position in each two field picture in the above-mentioned region-of-interest, find out each myocardium particle in the corresponding constantly radiofrequency signal section of difference, calculate the Integral Backscatter curve of each myocardium particle according to the radiofrequency signal section;
(8) according to the Integral Backscatter curve of each myocardium particle in the region-of-interest, calculate the Integral Backscatter fluctuating quantity of each myocardium particle, according to above-mentioned mapping table, obtain with region-of-interest in the corresponding color image of blood supply situation of each myocardium particle;
(9) above-mentioned color image is merged mutually with first frame of above-mentioned B ultrasonic image sequence, obtain the blood supply condition diagram of each myocardium particle in the region-of-interest.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100627654A CN1302751C (en) | 2004-07-09 | 2004-07-09 | Method for detecting myocardial blood supply state based on myocardial edge tracking |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100627654A CN1302751C (en) | 2004-07-09 | 2004-07-09 | Method for detecting myocardial blood supply state based on myocardial edge tracking |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1586404A CN1586404A (en) | 2005-03-02 |
| CN1302751C true CN1302751C (en) | 2007-03-07 |
Family
ID=34603776
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2004100627654A Expired - Fee Related CN1302751C (en) | 2004-07-09 | 2004-07-09 | Method for detecting myocardial blood supply state based on myocardial edge tracking |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1302751C (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009512486A (en) * | 2005-10-20 | 2009-03-26 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ultrasound imaging system and method |
| CN100418478C (en) * | 2006-06-08 | 2008-09-17 | 上海交通大学 | Virtual endoscope surface color mapping method based on blood flow imaging |
| JP5443719B2 (en) * | 2007-09-07 | 2014-03-19 | 株式会社東芝 | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program |
| JP5259267B2 (en) * | 2008-06-19 | 2013-08-07 | 株式会社東芝 | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program |
| JP5438936B2 (en) * | 2008-08-29 | 2014-03-12 | 株式会社東芝 | Ultrasonic diagnostic apparatus, image processing apparatus, and image processing program |
| US10729406B2 (en) * | 2015-03-10 | 2020-08-04 | Koninklijke Philips N.V. | Ultrasonic diagnosis of cardiac performance using heart model chamber segmentation with user control |
| WO2017008202A1 (en) * | 2015-07-10 | 2017-01-19 | 深圳迈瑞生物医疗电子股份有限公司 | Monitoring system, method and device |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5097836A (en) * | 1989-02-16 | 1992-03-24 | Fujitsu Limited | Untrasound diagnostic equipment for calculating and displaying integrated backscatter or scattering coefficients by using scattering power or scattering power spectrum of blood |
| JPH08249450A (en) * | 1995-03-14 | 1996-09-27 | Fujitsu Ltd | Image tracking device |
| US6005609A (en) * | 1996-11-15 | 1999-12-21 | Daewoo Electronics Co., Ltd. | Method and apparatus for digital correlation object tracker using a shape extraction focalization technique |
| CN1270794A (en) * | 2000-06-02 | 2000-10-25 | 清华大学 | Ultrasonic fusion imaging method integrating cardiac muscle's backward scattering |
| US6445832B1 (en) * | 2000-10-10 | 2002-09-03 | Lockheed Martin Corporation | Balanced template tracker for tracking an object image sequence |
| US6592522B2 (en) * | 2001-06-12 | 2003-07-15 | Ge Medical Systems Global Technology Company, Llc | Ultrasound display of displacement |
-
2004
- 2004-07-09 CN CNB2004100627654A patent/CN1302751C/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5097836A (en) * | 1989-02-16 | 1992-03-24 | Fujitsu Limited | Untrasound diagnostic equipment for calculating and displaying integrated backscatter or scattering coefficients by using scattering power or scattering power spectrum of blood |
| JPH08249450A (en) * | 1995-03-14 | 1996-09-27 | Fujitsu Ltd | Image tracking device |
| US6005609A (en) * | 1996-11-15 | 1999-12-21 | Daewoo Electronics Co., Ltd. | Method and apparatus for digital correlation object tracker using a shape extraction focalization technique |
| CN1270794A (en) * | 2000-06-02 | 2000-10-25 | 清华大学 | Ultrasonic fusion imaging method integrating cardiac muscle's backward scattering |
| US6445832B1 (en) * | 2000-10-10 | 2002-09-03 | Lockheed Martin Corporation | Balanced template tracker for tracking an object image sequence |
| US6592522B2 (en) * | 2001-06-12 | 2003-07-15 | Ge Medical Systems Global Technology Company, Llc | Ultrasound display of displacement |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1586404A (en) | 2005-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103889337B (en) | Diagnostic ultrasound equipment and ultrasonic diagnosis apparatus control method | |
| CN104797199B (en) | The standard flat assessed for real-time heart of fetus is automatically positioned | |
| Rohling et al. | Automatic registration of 3-D ultrasound images | |
| Leung et al. | Automated border detection in three-dimensional echocardiography: principles and promises | |
| JP5670324B2 (en) | Medical diagnostic imaging equipment | |
| US7542622B1 (en) | Spatio-temporal treatment of noisy images using brushlets | |
| US9717474B2 (en) | Image processing apparatus, ultrasound diagnosis apparatus, and image processing method | |
| US8107703B2 (en) | Quantitative real-time 4D stress test analysis | |
| US20100215238A1 (en) | Method for Automatic Segmentation of Images | |
| US20110257545A1 (en) | Imaging based symptomatic classification and cardiovascular stroke risk score estimation | |
| US20110257505A1 (en) | Atheromatic?: imaging based symptomatic classification and cardiovascular stroke index estimation | |
| CN104584074A (en) | Coupled segmentation in 3D conventional ultrasound and contrast-enhanced ultrasound images | |
| CN1859872A (en) | Ultrasonic cardiac volume quantification | |
| JP2009502354A (en) | Heart imaging system and method | |
| CN1862596A (en) | System and method for fused PET-CT visualization for heart unfolding | |
| US11191520B2 (en) | Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method | |
| US20180192987A1 (en) | Ultrasound systems and methods for automatic determination of heart chamber characteristics | |
| Ma et al. | Model driven quantification of left ventricular function from sparse single-beat 3D echocardiography | |
| CN1302751C (en) | Method for detecting myocardial blood supply state based on myocardial edge tracking | |
| US7376254B2 (en) | Method for surface-contouring of a three-dimensional image | |
| Hao et al. | Segmenting high-frequency intracardiac ultrasound images of myocardium into infarcted, ischemic, and normal regions | |
| McPherson et al. | Ultrasound characterization of acute myocardial ischemia by quantitative texture analysis | |
| Hafizah et al. | 3D ultrasound image reconstruction based on VTK | |
| Deng et al. | Automated detection of fetal cardiac structure from first-trimester ultrasound sequences | |
| Gastounioti et al. | Computerized analysis of ultrasound images: potential associations between texture and motion properties of the diseased arterial wall |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070307 Termination date: 20100709 |