CN104127207A - Method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging - Google Patents
Method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging Download PDFInfo
- Publication number
- CN104127207A CN104127207A CN201410413077.1A CN201410413077A CN104127207A CN 104127207 A CN104127207 A CN 104127207A CN 201410413077 A CN201410413077 A CN 201410413077A CN 104127207 A CN104127207 A CN 104127207A
- Authority
- CN
- China
- Prior art keywords
- blood flow
- region
- image
- flow rate
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention provides a method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging. The method comprises the steps that a blood flow speed image is obtained; a blood flow area of the blood flow speed image is obtained according to the blood flow speed image; a blood flow energy image is obtained; a blood flow area of the blood flow energy image is obtained according to the blood flow energy image; the blood flow area of the blood flow speed image and the blood flow area of the blood flow energy image are fused to form a new blood flow area image. Compared with the prior art, the method and system for detecting the sensitivity of the small blood flow signals based on ultrasonic imaging are a post-processing method for spatial features of the blood flow energy image and the blood flow speed image based on an ultrasonic Doppler, the blood flow signals and tissue signals can be distinguished effectively, and therefore the reflection echo is weak, the small blood flow signals with low energy are reserved, the sensitivity of color blood flow of an ultrasonic system is improved, and the sensitivity for detecting the small blood flow signals is further improved.
Description
Technical field
The present invention relates to a kind of method and system that detect little blood flow signal sensitivity based on ultra sonic imaging, be mainly used in ultrasonic color Doppler mapping technology field.
Background technology
Ultrasound investigation is to utilize human body to observe hyperacoustic reflection; , scan by certain direction to one group of ultrasound wave of human-body emitting by ultrasonic probe; According to the time delay of its echo of monitoring, strong and weak distance and the character that just can judge internal organs, pass through afterwards the processing of electronic circuit and computer, just form the B ultrasonic image at this position; Ultrasonographic cost ratio CT and MRI are cheap, have again without wound and obtain in real time people's in-vivo tissue feature of image, have become the main supplementary means of diagnosis, have a very wide range of applications clinically.
Color doppler ultrasonography (being color ultrasound) is generally to carry out Doppler signal processing with autocorrelation technique, and the blood flow signal that autocorrelation technique is obtained is superimposed upon in real time on two dimensional image after coloud coding, forms Color Doppler Ultrasonic picture.Its major advantage is: energy quicklook shows the two dimensional surface distribution of blood flow; Can show the traffic direction of blood flow; Be conducive to distinguish tremulous pulse and vein; Be conducive to identify vascular lesion and non-vascular lesion; Be conducive to understand the character of blood flow; Can conveniently understand time phase and the speed of blood flow; Can find reliably shunting and backflow; Can carry out quantitative analysis to the origin of blood flow bundle, width, length, area.As can be seen here, color ultrasound had both had advantages of two-dimensional ultrasound structural images, provided again hemodynamic abundant information, practical application to be subject to paying attention to widely and welcoming simultaneously, was described as clinically " atraumatic angiography ".
Because the blood flow in minute blood vessel is less, erythrocyte is also corresponding less, and the echo-signal of its backscattering is very faint; The impact of adding decay and the signal processing precision etc. of signal in tissue, often makes the tiny blood flow in colorful ultrasonic finally can not show, and be there is no blood flow thereby take for; This problem is commonly called blood flow poor sensitivity.Blood flow sensitivity problem shows particularly outstandingly on renal blood flow; The structure of renal blood flow is very complicated, the blood flow (being medically called fourth stage blood flow) of its end often cannot show or show incomplete, make user likely think that tip does not have blood flow, but in fact, this is for no other reason than that ultrasonic image-forming system blood flow poor sensitivity causes the form that cannot reproduce blood vessel itself.
Accordingly, blood flow sensitivity has become one of reference standard of weighing machine class, general, can reach and put forward highly sensitive object by changing and improve the synthetic and signal processing link of ultrasound emission and reception, wave beam, but these methods make the improved while of sensitivity, to the requirement of hardware also corresponding increase, thereby increase the hardware cost of machine.
Accordingly, the judgement of blood flow and tissue in ultrasonic image-forming system, mainly realizes by setting echo amplitude threshold value; Its basic thought is to set noise and two threshold values of tissue, when Doppler's energy is greater than noise threshold, be less than while organizing threshold value, and velocity magnitude is judged to be blood echo while meeting certain condition; This is that a kind of simple and coarse differentiating method all takes the differentiation of clean cut formula to be easy to occur wrong result to all positions, is not quite reasonable.
Therefore, realize that efficient to detect the method and system of little blood flow signal sensitivity based on ultra sonic imaging significantly still very challenging.
Summary of the invention
In order to address the above problem, the object of the present invention is to provide a kind of method and system that detect little blood flow signal sensitivity based on ultra sonic imaging, the method and system are the post-processing approachs of the space characteristics of blood flow energy image based on ultrasonic doppler and blood flow rate image, can effectively distinguish blood flow signal and tissue signal, thereby improve the sensitivity that little blood flow signal detects.
Accordingly, an embodiment of the present invention detect the method for little blood flow signal sensitivity based on ultra sonic imaging, said method comprising the steps of: obtain blood flow rate image; Obtain the blood flow region of described blood flow rate image according to described blood flow rate image;
Obtain blood flow energy image; Obtain the blood flow region of described blood flow energy image according to described blood flow energy image;
The blood flow region of the blood flow region of described blood flow rate image and described blood flow energy image is merged, form new blood flow area image.
As a further improvement on the present invention, described method is further comprising the steps of: described new blood flow area image is carried out to morphology processing, to remove after the noise signal of current described new blood flow area image, form new blood-stream image.
As a further improvement on the present invention, described " obtaining blood flow rate image " specifically comprises the following steps:
M1, obtain the rf echo signal of each sampled point in ultrasonoscopy;
M2, described rf echo signal is carried out to quadrature demodulation, and after low-pass filtering, obtain quadrature component I and in-phase component Q;
Noise signal in M3, filtering quadrature component I and in-phase component Q, to obtain the blood flow doppler signal of each sampled point in described ultrasonoscopy;
M4, according to described blood flow doppler signal, obtain the blood flow rate estimated value of each sampled point, and described blood flow rate estimated value spliced, form described blood flow rate image.
As a further improvement on the present invention, described " obtaining the blood flow region of described blood flow rate image according to described blood flow rate image " specifically comprises the following steps:
Judge in described blood flow rate image, whether the blood flow rate estimated value of each sampled point is null value;
If so, determine that current sampling point region is the background area of described blood flow rate image;
If not, determine that current sampling point region is blood flow region and the noise region of described blood flow rate image.
As a further improvement on the present invention, after described " if not, determining that current sampling point region is blood flow region and the noise region of described blood flow rate image ", described method is further comprising the steps of:
N1, described blood flow rate image is divided into K ROI region;
N2, calculate the variance of all sampled points in each ROI region
;
N3, calculate the variance of asking in described K ROI region
mean variance
;
N4, judge described variance
whether be greater than the first system threshold value,
If so, determine that current described ROI region is the noise region of described blood flow rate image;
If not, determine that current described ROI region is the blood flow region of described blood flow rate image;
Described the first system threshold value is mean variance
's
doubly, described in
for being greater than 0 natural number.
As a further improvement on the present invention, described " obtaining supersonic blood energy diagram picture " specifically comprises the following steps:
Obtain the backward energy value of each sampled point in ultrasonoscopy;
The backward energy value of the above-mentioned each sampled point obtaining is spliced, form described blood flow energy image.
As a further improvement on the present invention, described " obtaining the blood flow region of described blood flow energy image according to described blood flow energy image " specifically comprises the following steps:
Judge, in described blood flow energy image, whether the backward energy value of each sampled point is greater than second system predetermined threshold value, if so, determine that current sampling point region is the blood flow region of described energy diagram picture.
Correspondingly, an embodiment of the present invention detect the system of little blood flow signal sensitivity based on ultra sonic imaging, described system comprises: blood flow rate image generation module, for obtaining blood flow rate image;
The first blood flow region generation module, for obtaining the blood flow region of described blood flow rate image according to described blood flow rate image;
Blood flow energy image generation module, for obtaining blood flow energy image;
The second blood flow region generation module, for obtaining the blood flow region of described blood flow energy image according to described blood flow energy image;
Blood flow regional ensemble module, for the blood flow region of the blood flow region of described blood flow rate image and described blood flow energy image is merged, forms new blood flow area image.
As a further improvement on the present invention, described system also comprises: morphological process module, for described new blood flow area image is carried out to morphology processing, to remove after the noise signal of current described new blood flow area image, forms new blood-stream image.
As a further improvement on the present invention, described blood flow rate image generation module specifically for, obtain the rf echo signal of each sampled point in ultrasonoscopy;
Described rf echo signal is carried out to quadrature demodulation, and after low-pass filtering, obtain quadrature component I and in-phase component Q;
Noise signal in filtering quadrature component I and in-phase component Q, to obtain the blood flow doppler signal of each sampled point in described ultrasonoscopy;
According to described blood flow doppler signal, obtain the blood flow rate estimated value of each sampled point, and described blood flow rate estimated value is spliced, form described blood flow rate image.
As a further improvement on the present invention, described the first blood flow region generation module specifically for, judge in described blood flow rate image, whether the blood flow rate estimated value of each sampled point is null value;
If so, determine that current sampling point region is the background area of described blood flow rate image;
If not, determine that current sampling point region is blood flow region and the noise region of described blood flow rate image.
As a further improvement on the present invention, described the first blood flow region generation module also for, described blood flow rate image is divided into K ROI region;
Calculate the variance of all sampled points in each ROI region
;
Calculate the variance of asking in K ROI region of institute
mean variance
;
Described in judgement
whether be greater than the first system threshold value,
If so, determine that current described ROI region is the noise region of described blood flow rate image;
If not, determine that current described ROI region is the blood flow region of described blood flow rate image;
Described the first system threshold value is mean variance
's
doubly, described in
for being greater than 0 natural number.
As a further improvement on the present invention, described blood flow energy image generation module specifically for, obtain the backward energy value of each sampled point in ultrasonoscopy; The backward energy value of the above-mentioned each sampled point obtaining is spliced, form described blood flow energy image.
As a further improvement on the present invention, described the second blood flow region generation module specifically for, judge, in described blood flow energy image, whether the backward energy value of each sampled point is greater than second system predetermined threshold value, if so, determine that current sampling point region is the blood flow region of described energy diagram picture.
Compared with prior art, the method and system that detect little blood flow signal sensitivity based on ultra sonic imaging of the present invention, the method and system are the post-processing approachs of the space characteristics of blood flow energy image based on ultrasonic doppler and blood flow rate image, can effectively distinguish blood flow signal and tissue signal, thereby a little less than making reflection echo, the little blood flow signal that energy is very low is retained, and has improved the sensitivity of ultrasonic system colorful blood, further improves the sensitivity that little blood flow signal detects.
Brief description of the drawings
Fig. 1 is the flow chart that detects the method for little blood flow signal sensitivity in an embodiment of the present invention based on ultra sonic imaging;
Fig. 2 is the flow chart that obtains blood flow rate image in an embodiment of the present invention;
Fig. 3 is the flow chart that obtains the blood flow region of blood flow rate image in an embodiment of the present invention;
Fig. 4 is the structural representation that detects little blood flow signal sensitivity system in an embodiment of the present invention based on ultra sonic imaging.
Detailed description of the invention
Describe the present invention below with reference to each embodiment shown in the drawings.But these embodiments do not limit the present invention, the conversion in structure, method or function that those of ordinary skill in the art makes according to these embodiments is all included in protection scope of the present invention.
As shown in Figure 1, in an embodiment of the present invention, detect the method for little blood flow signal sensitivity based on ultra sonic imaging, said method comprising the steps of:
S1, obtain blood flow rate image.
Concrete, shown in Fig. 2, obtain blood flow rate image and specifically comprise the following steps:
M1, obtain the rf echo signal of each sampled point in ultrasonoscopy.
Concrete, in present embodiment, described rf echo signal is RF signal.
M2, described rf echo signal is carried out to quadrature demodulation, and after low-pass filtering, obtain quadrature component I and in-phase component Q.
General, the rf echo signal that corresponding each sampled point repeatedly obtains carries out quadrature demodulation, more accurate for the result that makes finally to obtain, generally the described rf echo signal obtaining for 8 to 16 times is carried out to quadrature demodulation, after low-pass filtering, obtain afterwards quadrature component I and in-phase component Q.
Accordingly, the rf echo signal each sampled point repeatedly being obtained carries out quadrature demodulation.
It should be noted that, the quadrature component I here and in-phase component Q are not single locational blood flow signal, but along the blood flow signal of degree of depth scanning direction diverse location, are referred to as direction when fast; In addition, the signal by identical sampling depth, forming is in the same time called signal when slow; Therefore, quadrature component I and in-phase component Q can discrete representation be
with
, the sampling number that n is depth direction, the sampling number that m is same position, the length of signal when slow.
In M3, filtering quadrature component I and in-phase component Q noise signal, to obtain the blood flow doppler signal of each sampled point in described ultrasonoscopy.
Accordingly, more accurate for the result that makes finally to obtain, noise signal in signal while also needing to remove each degree of depth slow.Concrete, during to each degree of depth slow, signal is by wall filtering, English by name: Clutter filter, the noise signal that the static or slow motor tissue of filtering produces, has just obtained blood flow doppler signal.
M4, according to described blood flow doppler signal, obtain the blood flow rate estimated value of each sampled point, and described blood flow rate estimated value spliced, form described blood flow rate image.
Concrete, blood flow doppler signal is carried out to the movement velocity of blood flow and estimate.
Accordingly, blood flow rate formula is:
Wherein, n represents the position at signal place,
represent the speed that ultrasound wave is propagated in tissue,
the mid frequency transmitting,
the angle between probe scanning direction and blood flow direction,
it is the Doppler frequency shift that echo-signal is produced due to blood flow.
Accordingly, obtained the value at n place, signal position by formula (2).
Wherein, T is the repetition period of ultrasonic emitting pulse, and arg represents to get the phase place of complex signal, and span is [π, π].
Wherein, j is imaginary symbols.
Wherein, * represents to get plural conjugation, and average on M express time launched M subpulse continuously at the same space position n, and the M subpulse echo-signal at n place, position is averaged.
In conjunction with above-mentioned 4 formula, can obtain the blood flow rate estimated value of each sampled point in ultrasonoscopy
, afterwards, described blood flow rate estimated value is spliced to form to a width two dimensional image, that is: form described blood flow rate image.
S2, the blood flow region that obtains described blood flow rate image according to described blood flow rate image.
Accordingly, described blood flow rate image comprises: the background area of blood flow rate image, the blood flow region of the noise region of blood flow rate image and blood flow rate image.
The actual reflection of blood flow rate image be the phase diagram of signal, according to the spatial distribution characteristic of described blood flow rate image, phase place is zero region as a setting, region, and adjacent signals phase place changes faster region as noise region, and adjacent signals phase place changes slower region as blood flow region.
Accordingly, specifically comprise the following steps in conjunction with described step S2 shown in Fig. 3:
P1, judge in described blood flow rate image, whether the blood flow rate estimated value of each sampled point is null value;
If so, determine that current sampling point region is the background area of described blood flow rate image;
If not, determine that current sampling point region is blood flow region and the noise region of described blood flow rate image;
Further, in the blood flow region and noise region of isolated described blood flow rate image, the blood flow region that obtains described blood flow rate image.
P2, differentiate by the variance size in ROI region in blood flow rate image blood flow region or the noise region that current region is under the jurisdiction of blood flow rate image.
Described ROI region is the abbreviation of English Region Of Interest, translates into Chinese to be: area-of-interest.
Concrete:
N1, described blood flow rate image is divided into K ROI region;
The sampled point size in described single ROI region is consistent with the size of the space smoothing of system default.
N2, calculate the variance of all sampled points in each ROI region
(k=1,2 ... K);
Accordingly, in whole blood flow rate image, total K ROI region, described K is positive integer.
N3, calculate the variance of asking in described K ROI region
mean variance
.
N4, judge described variance
whether be greater than the first system threshold value,
If so, determine that current described ROI region is the noise region of described blood flow rate image;
If not, determine that current described ROI region is the blood flow region of described blood flow rate image.
Described the first system threshold value is mean variance
's
doubly, described in
for being greater than 0 natural number.
S3, obtain blood flow energy image.
Concrete, obtain the backward energy value of each sampled point in ultrasonoscopy;
In the specific embodiment of the present invention, obtain the backward energy value of each sampled point in ultrasonoscopy by following formula.
Accordingly, described backward energy value is:
Wherein, the sampling number that n is depth direction, the sampling number that m is same position, average on M express time, has launched M subpulse continuously at the same space position n, and the M subpulse echo-signal at n place, position is averaged.
Further, the backward energy value of the above-mentioned each sampled point obtaining is spliced to form to a width two dimensional image, that is: forms described blood flow energy image.
S4, the blood flow region that obtains described blood flow energy image according to described blood flow energy image.
Judge, in described blood flow energy image, whether the backward energy value of each sampled point is greater than second system predetermined threshold value
, if so, determine that current sampling point region is the blood flow region of described energy diagram picture.
Described
be a flexible adjustable parameter, its size is set according to actual needs voluntarily.
S5, the blood flow region of the blood flow region of described blood flow rate image and described blood flow energy image is merged, form new blood flow area image.
Described new blood flow area image can be expressed as the region union in the blood flow region of described blood flow rate image and the blood flow region of described blood flow energy image in one embodiment of the present invention;
It should be noted that, the blood flow region in described blood flow rate image comprises: large blood flow region and little blood flow region; In blood flow rate image, ultrasonic signal is through a series of conversion, finally by cross energy after wall filtering still higher signal be large blood flow signal; Meanwhile, because blood flow signal has distribution seriality, little blood flow signal is distributed in the tip that large blood flow extends conventionally; Therefore, in the described blood flow rate image obtaining, can ensure signal higher energy to remain, the continuous signal that distributes on hodograph can be remained again, this just can ensure to make the weak little blood flow of part energy to be unlikely to lose; Therefore, described new blood flow area image, can effectively distinguish blood flow signal and tissue signal, thereby a little less than making reflection echo, the little blood flow signal that energy is very low is retained, and has improved the sensitivity of ultrasonic system colorful blood, further improves the sensitivity that little blood flow signal detects.
S6, described new blood flow area image is carried out to morphology processing, to remove the noise signal of current described new blood flow area image.
Be understandable that, extract the blood flow region of described blood flow rate image and extract in the blood flow region process of described blood flow energy image at step S4 at step S2, have unavoidably some noise signals residual, therefore, in the final new blood flow area image obtaining, also have unavoidably some noise signals residual.
Accordingly, new blood flow area image is carried out to denoising, to form new blood-stream image.
Concrete, described new blood flow area image is carried out to morphology processing, remove the segment that area is less than color ultrasound imaging resolution, to remove after the noise signal in current described new blood flow area image, form new blood-stream image.
It should be noted that, the above-mentioned step of asking for new blood-stream image is not unique, for example: also can first carry out step S3, S4, carry out at step S1, S2 again, or first carry out S1, S3, carry out again S2, S4 etc., certainly can also adopt other order to carry out, will not enumerate at this, the change of its step, can't affect the final result that will obtain, is not described in detail at this.
Compared with prior art, the method that detects little blood flow signal sensitivity based on ultra sonic imaging of the present invention, the method is the post-processing approach of the space characteristics of blood flow energy image based on ultrasonic doppler and blood flow rate image, can effectively distinguish blood flow signal and tissue signal, thereby a little less than making reflection echo, the little blood flow signal that energy is very low is retained, and has improved the sensitivity of ultrasonic system colorful blood, further improves the sensitivity that little blood flow signal detects.
Accordingly, shown in Fig. 4, Fig. 4 is the structural representation that detects little blood flow signal sensitivity system in an embodiment of the present invention based on ultra sonic imaging.
Accordingly, described system comprises: blood flow rate image generation module 100, the first blood flow region generation modules 200, blood flow energy image generation module 300, the second blood flow region generation modules 400, blood flow regional ensemble module 500, morphological process module 600.
Corresponding blood flow rate image generation module 100 is for obtaining blood flow rate image.
Blood flow rate image generation module 100 is for obtaining the rf echo signal of the each sampled point of ultrasonoscopy.
Concrete, in present embodiment, described rf echo signal is RF signal.
Blood flow rate image generation module 100 is for described rf echo signal is carried out to quadrature demodulation, and after low-pass filtering, obtains quadrature component I and in-phase component Q.
General, corresponding each sampled point, the rf echo signal that blood flow rate image generation module obtains for more than 100 time carries out quadrature demodulation, for the result that makes finally to obtain more accurate, generally the described rf echo signal obtaining for 8 to 16 times is carried out to quadrature demodulation, after low-pass filtering, obtain afterwards quadrature component I and in-phase component Q.
Accordingly, the rf echo signal each sampled point repeatedly being obtained carries out quadrature demodulation.
It should be noted that, the quadrature component I here and in-phase component Q are not single locational blood flow signal, but along the blood flow signal of degree of depth scanning direction diverse location, are referred to as direction when fast; In addition, the signal by identical sampling depth, forming is in the same time called signal when slow; Therefore, quadrature component I and in-phase component Q can discrete representation be
with
, the sampling number that n is depth direction, the sampling number that m is same position, the length of signal when slow.
Blood flow rate image generation module 100 is for the noise signal of filtering quadrature component I and in-phase component Q, to obtain the blood flow doppler signal of each sampled point in described ultrasonoscopy.
Accordingly, for the result that makes finally to obtain more accurate, noise signal in signal when blood flow rate image generation module 100 also needs to remove each degree of depth slow.The noise signal concrete, during to each degree of depth slow, signal produces by the static or slow motor tissue of wall filtering filtering, has just obtained blood flow doppler signal.
Blood flow rate image generation module 100, for according to described blood flow doppler signal, obtains the blood flow rate estimated value of each sampled point, and described blood flow rate estimated value is spliced, and forms described blood flow rate image.
Concrete, the movement velocity that blood flow rate image generation module 100 carries out blood flow to blood flow doppler signal is estimated.
Accordingly, blood flow rate formula is:
Wherein, n represents the position at signal place,
represent the speed that ultrasound wave is propagated in tissue,
the mid frequency transmitting,
the angle between probe scanning direction and blood flow direction,
it is the Doppler frequency shift that echo-signal is produced due to blood flow.
Accordingly, obtained the value at n place, signal position by formula (2).
Wherein, T is the repetition period of ultrasonic emitting pulse, and arg represents to get the phase place of complex signal, and span is [π, π].
Wherein, j is imaginary symbols.
Wherein, * represents to get plural conjugation, and average on M express time launched M subpulse continuously at the same space position n, and the M subpulse echo-signal at n place, position is averaged.
In conjunction with above-mentioned 4 formula, blood flow rate image generation module 100 can obtain the blood flow rate estimated value of each sampled point in ultrasonoscopy
, afterwards, described blood flow rate estimated value is spliced to form to a width two dimensional image, that is: form described blood flow rate image.
Accordingly, the first blood flow region generation module 200 is for obtaining the blood flow region of described blood flow rate image according to described blood flow rate image.
Accordingly, described blood flow rate image comprises: the background area of blood flow rate image, the blood flow region of the noise region of blood flow rate image and blood flow rate image.
The actual reflection of blood flow rate image be the phase diagram of signal, according to the spatial distribution characteristic of described blood flow rate image, phase place is zero region as a setting, region, and adjacent signals phase place changes faster region as noise region, and adjacent signals phase place changes slower region as blood flow region.
Concrete, the first blood flow region generation module 200 judges in described blood flow rate image, and whether the blood flow rate estimated value of each sampled point is nonzero value;
If so, the first blood flow region generation module 200 determines that current sampling point region is the background area of described blood flow rate image;
If not, the first blood flow region generation module 200 determines that current sampling point region is blood flow region or the noise region of described blood flow rate image;
Further, the first blood flow region generation module 200 in the blood flow region and noise region of isolated described blood flow rate image, the blood flow region that obtains described blood flow rate image.
Accordingly, in the present embodiment, the first blood flow region generation module 200 differentiates by the variance size in ROI region in blood flow rate image blood flow region or the noise region that current region is under the jurisdiction of blood flow rate image.
Described ROI region is the abbreviation of English Region Of Interest, translates into Chinese to be: local region of interest.
Concrete: the first blood flow region generation module 200 is for being divided into K ROI region by described blood flow rate image;
The sampled point size in described single ROI region is consistent with the size of the space smoothing of system default.
The first blood flow region generation module 200 is for calculating the variance of all sampled points in each ROI region
(k=1,2 ... K);
Accordingly, in whole blood flow rate image, total K ROI region, described K is positive integer.
The first blood flow region generation module 200 is for calculating the variance of asking in described K ROI region
mean variance
.
The first blood flow region generation module 200 is for judging described variance
whether be greater than the first system threshold value,
If so, the first blood flow region generation module 200 determines that current described ROI region is the noise region of described blood flow rate image;
If not, the first blood flow region generation module 200 determines that current described ROI region is the blood flow region of described blood flow rate image.
Described the first system threshold value is mean variance
's
doubly, described in
for being greater than 0 natural number.Blood flow energy image generation module 300 is for obtaining blood flow energy image.
Concrete, blood flow energy image generation module 300 is for obtaining the backward energy value of the each sampled point of ultrasonoscopy;
In a detailed description of the invention of the present invention, blood flow energy image generation module 300 obtains the backward energy value of each sampled point in ultrasonoscopy by following formula.
Accordingly, described backward energy value is:
Wherein, the sampling number that n is depth direction, the sampling number that m is same position, average on M express time, has launched M subpulse continuously at the same space position n, and the M subpulse echo-signal at n place, position is averaged.
Further, the backward energy value of the above-mentioned each sampled point obtaining is spliced to form a width two dimensional image by blood flow energy image generation module 300, forms described blood flow energy image that is:.
The second blood flow region generation module 400 is for obtaining the blood flow region of described blood flow energy image according to described blood flow energy image.
The second blood flow region generation module 400 is for judging described blood flow energy image, and whether the backward energy value of each sampled point is greater than second system predetermined threshold value
, if so, determine that current sampling point region is the blood flow region of described energy diagram picture.Described
be a flexible adjustable parameter, its size is set according to actual needs voluntarily.
Blood flow regional ensemble module 500, for the blood flow region of the blood flow region of described blood flow rate image and described blood flow energy image is merged, forms new blood flow area image.
Described new blood flow area image can be expressed as the region union in the blood flow region of described blood flow rate image and the blood flow region of described blood flow energy image in one embodiment of the present invention;
It should be noted that, the blood flow region in described blood flow rate image comprises: large blood flow region and little blood flow region; In blood flow rate image, ultrasonic signal is through a series of conversion, finally by cross energy after wall filtering still higher signal be large blood flow signal; Meanwhile, because blood flow signal has distribution seriality, little blood flow signal is distributed in the tip that large blood flow extends conventionally; Therefore, in the described blood flow rate image obtaining, can ensure signal higher energy to remain, the continuous signal that distributes on hodograph can be remained again, this just can ensure to make the weak little blood flow of part energy to be unlikely to lose; Therefore, described new blood flow area image, can effectively distinguish blood flow signal and tissue signal, thereby a little less than making reflection echo, the little blood flow signal that energy is very low is retained, and has improved the sensitivity of ultrasonic system colorful blood, further improves the sensitivity that little blood flow signal detects.
Morphological process module 600, for described new blood flow area image is carried out to morphology processing, to remove after the noise signal in current described new blood flow region, forms new blood-stream image.
Be understandable that, blood flow region and the second blood flow region generation module 400 that the first blood flow region generation module 200 extracts described blood flow rate image extract in the blood flow region process of described blood flow energy image, have unavoidably some noise signals residual, therefore, in the final new blood flow area image obtaining of blood flow regional ensemble module 500, also have unavoidably some noise signals residual.
Accordingly, morphological process module 600 is carried out denoising to new blood flow area image, to form new blood-stream image.
Concrete, morphological process module 600 is carried out morphology processing to described new blood flow area image, removes the segment that area is less than color ultrasound imaging resolution, to remove after the noise signal of current described new blood flow area image, forms new blood-stream image.
Compared with prior art, the method and system that detect little blood flow signal sensitivity based on ultra sonic imaging of the present invention, the method is the post-processing approach of the space characteristics of blood flow energy image based on ultrasonic doppler and blood flow rate image, can effectively distinguish blood flow signal and tissue signal, thereby a little less than making reflection echo, the little blood flow signal that energy is very low is retained, and has improved the sensitivity of ultrasonic system colorful blood, further improves the sensitivity that little blood flow signal detects.
For convenience of description, while describing above device, being divided into various modules with function describes respectively.Certainly, in the time implementing the application, the function of each module can be realized in same or multiple software and/or hardware.
As seen through the above description of the embodiments, those skilled in the art can be well understood to the mode that the application can add essential general hardware platform by software and realizes.Based on such understanding, the part that the application's technical scheme contributes to prior art in essence in other words can embody with the form of software product, this computer software product can be kept in Protector, as ROM/RAM, magnetic disc, CD etc., comprise that some instructions (can be personal computers in order to make a computer equipment, Information Push Server, or the network equipment etc.) carry out the method described in some part of each embodiment of the application or embodiment.
Device embodiments described above is only schematic, the wherein said module as separating component explanation can or can not be also physically to separate, the parts that show as module can be or can not be also physical modules, can be positioned at a place, or also can be distributed on multiple mixed-media network modules mixed-medias.Can select according to the actual needs some or all of module wherein to realize the object of present embodiment scheme.Those of ordinary skill in the art, in the situation that not paying creative work, are appreciated that and implement.
The application can be used in numerous general or special purpose computingasystem environment or configuration.For example: personal computer, Information Push Server computer, handheld device or portable set, laptop device, multi-processing module system, system based on micro treatment module, set top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer, the distributed computing environment that comprises above any system or equipment etc.
The application can describe in the general context of computer executable instructions, for example program module.Usually, program module comprises and carries out particular task or realize routine, program, object, assembly, data structure of particular abstract data type etc.Also can in distributed computing environment, put into practice the application, in these distributed computing environment, be executed the task by the teleprocessing equipment being connected by communication network.In distributed computing environment, program module can be arranged in the local and remote computer Protector including preservation equipment.
Be to be understood that, although this description is described according to embodiment, but be not that each embodiment only comprises an independently technical scheme, this narrating mode of description is only for clarity sake, those skilled in the art should make description as a whole, technical scheme in each embodiment also can, through appropriately combined, form other embodiments that it will be appreciated by those skilled in the art that.
Listed a series of detailed description is above only illustrating for feasibility embodiment of the present invention; they are not in order to limit the scope of the invention, all do not depart from the equivalent embodiment that skill spirit of the present invention does or change and all should be included in protection scope of the present invention within.
Claims (14)
1. a method that detects little blood flow signal sensitivity based on ultra sonic imaging, is characterized in that, said method comprising the steps of:
Obtain blood flow rate image; Obtain the blood flow region of described blood flow rate image according to described blood flow rate image;
Obtain blood flow energy image; Obtain the blood flow region of described blood flow energy image according to described blood flow energy image;
The blood flow region of the blood flow region of described blood flow rate image and described blood flow energy image is merged, form new blood flow area image.
2. the method that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 1, is characterized in that, described method is further comprising the steps of:
Described new blood flow area image is carried out to morphology processing, to remove after the noise signal of current described new blood flow area image, form new blood-stream image.
3. the method that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 1, is characterized in that, described " obtaining blood flow rate image " specifically comprises the following steps:
M1, obtain the rf echo signal of each sampled point in ultrasonoscopy;
M2, described rf echo signal is carried out to quadrature demodulation, and after low-pass filtering, obtain quadrature component I and in-phase component Q;
Noise signal in M3, filtering quadrature component I and in-phase component Q, to obtain the blood flow doppler signal of each sampled point in described ultrasonoscopy;
M4, according to described blood flow doppler signal, obtain the blood flow rate estimated value of each sampled point, and described blood flow rate estimated value spliced, form described blood flow rate image.
4. the method that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 3, is characterized in that, described " obtaining the blood flow region of described blood flow rate image according to described blood flow rate image " specifically comprises the following steps:
Judge in described blood flow rate image, whether the blood flow rate estimated value of each sampled point is null value,
If so, determine that current sampling point region is the background area of described blood flow rate image;
If not, determine that current sampling point region is blood flow region and the noise region of described blood flow rate image.
5. the method that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 4, it is characterized in that, after described " if not, determining that current sampling point region is blood flow region and the noise region of described blood flow rate image ", described method is further comprising the steps of:
N1, described blood flow rate image is divided into K ROI region;
N2, calculate the variance of all sampled points in each ROI region
;
N3, calculate the variance of asking in described K ROI region
mean variance
;
N4, judge described variance
whether be greater than the first system threshold value,
If so, determine that current described ROI region is the noise region of described blood flow rate image;
If not, determine that current described ROI region is the blood flow region of described blood flow rate image;
Described the first system threshold value is mean variance
's
doubly, described in
for being greater than 0 natural number.
6. the method that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 1, is characterized in that, described " obtaining supersonic blood energy diagram picture " specifically comprises the following steps:
Obtain the backward energy value of each sampled point in ultrasonoscopy;
The backward energy value of the above-mentioned each sampled point obtaining is spliced, form described blood flow energy image.
7. the method that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 6, is characterized in that, described " obtaining the blood flow region of described blood flow energy image according to described blood flow energy image " specifically comprises the following steps:
Judge, in described blood flow energy image, whether the backward energy value of each sampled point is greater than second system predetermined threshold value, if so, determine that current sampling point region is the blood flow region of described energy diagram picture.
8. a system that detects little blood flow signal sensitivity based on ultra sonic imaging, is characterized in that,
Described system comprises: blood flow rate image generation module, for obtaining blood flow rate image;
The first blood flow region generation module, for obtaining the blood flow region of described blood flow rate image according to described blood flow rate image;
Blood flow energy image generation module, for obtaining blood flow energy image;
The second blood flow region generation module, for obtaining the blood flow region of described blood flow energy image according to described blood flow energy image;
Blood flow regional ensemble module, for the blood flow region of the blood flow region of described blood flow rate image and described blood flow energy image is merged, forms new blood flow area image.
9. the system that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 8, is characterized in that,
Described system also comprises: morphological process module, for described new blood flow area image is carried out to morphology processing, to remove after the noise signal of current described new blood flow area image, forms new blood-stream image.
10. the system that detects little blood flow signal sensitivity based on ultra sonic imaging according to claim 8, is characterized in that,
Described blood flow rate image generation module specifically for, obtain the rf echo signal of each sampled point in ultrasonoscopy;
Described rf echo signal is carried out to quadrature demodulation, and after low-pass filtering, obtain quadrature component I and in-phase component Q;
Noise signal in filtering quadrature component I and in-phase component Q, to obtain the blood flow doppler signal of each sampled point in described ultrasonoscopy;
According to described blood flow doppler signal, obtain the blood flow rate estimated value of each sampled point, and described blood flow rate estimated value is spliced, form described blood flow rate image.
11. systems that detect little blood flow signal sensitivity based on ultra sonic imaging according to claim 10, is characterized in that,
Described the first blood flow region generation module specifically for, judge in described blood flow rate image, whether the blood flow rate estimated value of each sampled point is null value;
If so, determine that current sampling point region is the background area of described blood flow rate image;
If not, determine that current sampling point region is blood flow region and the noise region of described blood flow rate image.
12. systems that detect little blood flow signal sensitivity based on ultra sonic imaging according to claim 11, is characterized in that,
Described the first blood flow region generation module also for, described blood flow rate image is divided into K ROI region;
Calculate the variance of all sampled points in each ROI region
;
Calculate the variance of asking in K ROI region of institute
mean variance
;
Described in judgement
whether be greater than the first system threshold value,
If so, determine that current described ROI region is the noise region of described blood flow rate image;
If not, determine that current described ROI region is the blood flow region of described blood flow rate image;
Described the first system threshold value is mean variance
's
doubly, described in
for being greater than 0 natural number.
13. systems that detect little blood flow signal sensitivity based on ultra sonic imaging according to claim 8, is characterized in that,
Described blood flow energy image generation module specifically for, obtain the backward energy value of each sampled point in ultrasonoscopy;
The backward energy value of the above-mentioned each sampled point obtaining is spliced, form described blood flow energy image.
14. systems that detect little blood flow signal sensitivity based on ultra sonic imaging according to claim 13, is characterized in that,
Described the second blood flow region generation module specifically for, judge, in described blood flow energy image, whether the backward energy value of each sampled point is greater than second system predetermined threshold value, if so, determine that current sampling point region is the blood flow region of described energy diagram picture.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410413077.1A CN104127207A (en) | 2014-08-21 | 2014-08-21 | Method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410413077.1A CN104127207A (en) | 2014-08-21 | 2014-08-21 | Method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104127207A true CN104127207A (en) | 2014-11-05 |
Family
ID=51800247
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410413077.1A Pending CN104127207A (en) | 2014-08-21 | 2014-08-21 | Method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104127207A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107913082A (en) * | 2015-11-10 | 2018-04-17 | 三星麦迪森株式会社 | Supersonic imaging apparatus and the method for operating the supersonic imaging apparatus |
| CN109360173A (en) * | 2018-11-23 | 2019-02-19 | 珠海医凯电子科技有限公司 | Color Doppler Flow Imaging based on improved variance is as noise-reduction method |
| CN109512466A (en) * | 2018-12-08 | 2019-03-26 | 余姚市华耀工具科技有限公司 | Intelligent gynecological B ultrasound instrument |
| CN110169787A (en) * | 2019-06-28 | 2019-08-27 | 深圳开立生物医疗科技股份有限公司 | A kind of blood flow mixed display method and a kind of ultrasonic device |
| CN111388010A (en) * | 2020-03-26 | 2020-07-10 | 深圳开立生物医疗科技股份有限公司 | Ultrasonic Doppler blood flow imaging method, device, equipment and readable storage medium |
| US10898163B2 (en) | 2015-11-10 | 2021-01-26 | Samsung Medison Co., Ltd. | Ultrasound imaging apparatus and method of operating the same |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102429684A (en) * | 2010-09-28 | 2012-05-02 | 深圳迈瑞生物医疗电子股份有限公司 | Doppler color blood flow imaging method and device |
| CN102805649A (en) * | 2011-06-03 | 2012-12-05 | 深圳迈瑞生物医疗电子股份有限公司 | Color ultrasonic imaging method and device |
| US20130194267A1 (en) * | 2010-10-28 | 2013-08-01 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus and ultrasonic image display method |
| US8649586B2 (en) * | 2010-12-30 | 2014-02-11 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Color flow gain adjustment method and device |
-
2014
- 2014-08-21 CN CN201410413077.1A patent/CN104127207A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102429684A (en) * | 2010-09-28 | 2012-05-02 | 深圳迈瑞生物医疗电子股份有限公司 | Doppler color blood flow imaging method and device |
| US20130194267A1 (en) * | 2010-10-28 | 2013-08-01 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus and ultrasonic image display method |
| US8649586B2 (en) * | 2010-12-30 | 2014-02-11 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Color flow gain adjustment method and device |
| CN102805649A (en) * | 2011-06-03 | 2012-12-05 | 深圳迈瑞生物医疗电子股份有限公司 | Color ultrasonic imaging method and device |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107913082A (en) * | 2015-11-10 | 2018-04-17 | 三星麦迪森株式会社 | Supersonic imaging apparatus and the method for operating the supersonic imaging apparatus |
| US10898163B2 (en) | 2015-11-10 | 2021-01-26 | Samsung Medison Co., Ltd. | Ultrasound imaging apparatus and method of operating the same |
| CN109360173A (en) * | 2018-11-23 | 2019-02-19 | 珠海医凯电子科技有限公司 | Color Doppler Flow Imaging based on improved variance is as noise-reduction method |
| CN109360173B (en) * | 2018-11-23 | 2021-09-03 | 珠海医凯电子科技有限公司 | Color Doppler blood flow image noise reduction method based on improved variance |
| CN109512466A (en) * | 2018-12-08 | 2019-03-26 | 余姚市华耀工具科技有限公司 | Intelligent gynecological B ultrasound instrument |
| CN110169787A (en) * | 2019-06-28 | 2019-08-27 | 深圳开立生物医疗科技股份有限公司 | A kind of blood flow mixed display method and a kind of ultrasonic device |
| CN110169787B (en) * | 2019-06-28 | 2022-04-01 | 深圳开立生物医疗科技股份有限公司 | Blood flow mixing display method and ultrasonic equipment |
| CN111388010A (en) * | 2020-03-26 | 2020-07-10 | 深圳开立生物医疗科技股份有限公司 | Ultrasonic Doppler blood flow imaging method, device, equipment and readable storage medium |
| CN111388010B (en) * | 2020-03-26 | 2022-06-24 | 深圳开立生物医疗科技股份有限公司 | Ultrasonic Doppler blood flow imaging method, device, equipment and readable storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104127207A (en) | Method and system for detecting sensitivity of small blood flow signals based on ultrasonic imaging | |
| US8864671B2 (en) | Methods and systems for color flow imaging | |
| US9351707B2 (en) | Methods and apparatus to determine shear wave propagation property | |
| US8050476B2 (en) | Heart rate demodulation of periodic movement in ultrasound data for object identification | |
| US8684934B2 (en) | Adaptively performing clutter filtering in an ultrasound system | |
| US20170164923A1 (en) | Image Processor, Ultrasound Diagnostic Device Including Same, And Image Processing Method | |
| CN103908297B (en) | For identifying ultrasonic image-forming system from the data of shadow region and method | |
| EP2016905B1 (en) | Ultrasound diagnostic apparatus | |
| EP2830508B1 (en) | Methods and apparatus for ultrasound imaging | |
| EP3621524B1 (en) | Fetal development monitoring | |
| US20140276003A1 (en) | Systems and Methods to Detect and Present Interventional Devices via Ultrasound Imaging | |
| EP2233945A1 (en) | Surface rendering for volume data in an ultrasound system | |
| EP2609868A1 (en) | Providing user interface in ultrasound system | |
| EP2016906A2 (en) | Ultrasound diagnostic apparatus | |
| JP2011152416A (en) | Ultrasonic system and method for processing three-dimensional ultrasonic screen image and measuring size of concerned object | |
| Wells | Ultrasonic colour flow imaging | |
| CN105212964B (en) | Based on RF data ultrasonic imaging processing method and system | |
| US20100280386A1 (en) | Ultrasound System And Method For Adaptively Performing Clutter Filtering | |
| US20170105699A1 (en) | Ultrasonic diagnostic device | |
| WO2021018107A1 (en) | Ultrasonic signal processing method and ultrasonic signal processing apparatus, and device and storage medium | |
| JP5999997B2 (en) | Ultrasonic diagnostic equipment | |
| CN108652666A (en) | A kind of generation method and device of doppler blood flow imaging | |
| JP7440328B2 (en) | Ultrasound diagnostic equipment and programs | |
| EP2380498A1 (en) | Adaptive clutter filtering in an ultrasound system | |
| Wells | Current Doppler technology and techniques |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20141105 |