CN104720850B - Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture - Google Patents
Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture Download PDFInfo
- Publication number
- CN104720850B CN104720850B CN201310719832.4A CN201310719832A CN104720850B CN 104720850 B CN104720850 B CN 104720850B CN 201310719832 A CN201310719832 A CN 201310719832A CN 104720850 B CN104720850 B CN 104720850B
- Authority
- CN
- China
- Prior art keywords
- signal
- contrast
- tissue
- contrastographic picture
- region
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 238000001514 detection method Methods 0.000 title claims abstract description 17
- 239000002872 contrast media Substances 0.000 claims abstract description 67
- 231100000732 tissue residue Toxicity 0.000 claims abstract description 63
- 230000009467 reduction Effects 0.000 claims abstract description 35
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000009471 action Effects 0.000 claims abstract description 13
- 238000002601 radiography Methods 0.000 claims abstract description 9
- 230000011218 segmentation Effects 0.000 claims abstract description 8
- 238000009826 distribution Methods 0.000 claims description 11
- 238000012805 post-processing Methods 0.000 claims description 10
- 230000008520 organization Effects 0.000 claims description 6
- 239000000523 sample Substances 0.000 claims description 5
- 230000008569 process Effects 0.000 description 14
- 230000000694 effects Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 239000003550 marker Substances 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000025366 tissue development Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agents, e.g. microbubbles introduced into the bloodstream
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
- G01S7/52039—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target exploiting the non-linear response of a contrast enhancer, e.g. a contrast agent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Nonlinear Science (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention provides a kind of ultrasonic contrast imaging method and region detection, the developing method of contrastographic picture, the method for detecting area of the contrastographic picture includes：Tissue signal and reduction of contrast signal during acquisition ultrasonic contrast imaging；Tissue signal and reduction of contrast signal are pre-processed respectively；Pretreated tissue signal and reduction of contrast signal are compared processing；The image that is formed after processing will be compared by way of Threshold segmentation and be partitioned into the contrast agent zone of action, the tissue residue zone of action and noise contributions region；Corresponding position on contrastographic picture is respectively labeled as contrast medium areas, tissue residue region and noise region.One aspect of the present invention substantially increases the CTR of ultrasonic contrast image, is conducive to clinical observation；On the other hand the signal of contrast agent is enhanced, the radiography dosage injected during ultrasonic contrast is substantially reduced, effectively reduces check cost.
Description
Technical field
The present invention relates to ultrasound contrast imaging technology, and in particular to the area of a kind of ultrasonic contrast imaging method and contrastographic picture
Domain detection, developing method.
Background technology
Ultrasonic imaging is the Non-invasive detection mode of current comparative maturity, compared to CT and MRI (Magnetic Resonance
Imaging, magnetic resonance imaging), the radiation of Ultrasonic Diagnosis is few, and expense is cheap, and the advantage such as easy to operate causes Ultrasonic Diagnosis to turn into
Clinical primary detection methods.
Ultrasonic contrast is to strengthen scattered echoes using contrast agent, can effectively improve resolving power, the sensitivity of Ultrasonic Diagnosis
Property and specificity.The core of ultrasound contrast imaging technology is how to extract the contrast agent echo component and suppression group in echo-signal
The echo component knitted.At present, the imaging for carrying out contrast agent microbubble in non-linear detection method is the main side of current ultrasonic contrast
Method.
However, using ultrasonic contrast method of the prior art, unavoidable be mixed into is organized into the contrast agent signals detected
Point.
One is due to the cancellation performance of circuit system in itself, and two amplitudes, frequencies are the same in theory, phase difference 180 degree
Waveform adder be 0, but the direct impulse launched in actual circuit and the reverse impulse that 180 degree is differed with direct impulse phase
It is not full symmetric, it behind 0, that is, forward and reverse pulse cancellation is not 0 to cause after forward and reverse impulse summation and be not.Phase
Disappearing property influences the CTR (contrast to tissue ratio, radiography tissue ratio) of contrastographic picture, CTR to be to weigh contrastographic picture
One major criterion of quality, system cancellation performance is better, and CTR is higher, and the gap of contrast agent and tissue is bigger, image table
It is exactly now that tissue development is weaker, the level of contrast agent images is apparent.On the contrary, the poorer CTR of the destructive of system is poorer, contrast agent
Gap with tissue is smaller, i.e., the structural constituent remained in contrastographic picture is more, and the contrast resolution of contrastographic picture is poorer.
Two be due to that tissue can also produce non-linear component, and no matter how low MI (Mechanical Index, mechanical index) is,
All the non-linear component of tissue can be contained during non-linear detection with the presence of a small amount of harmonic components in tissue.
Three be due to that signal is too big, and each system has signal after the cancellation performance of itself, small-signal cancellation small, strong signal
Tissue residue after cancellation is also very strong, and strong tissue residue causes reduction of contrast signal and tissue residue cannot be distinguished by.
Four be due to that amplitude modulation(PAM) is used in saturation, non-linear detection method, launches the different pulse of several amplitudes, receives
When echo-signal will be readjusted by a narrow margin according to the amplitude scale between high and low pulse, and subtract each other with high-amplitude echo-signal, so
The fundamental signal of tissue can be eliminated；Saturation makes high and low width signal be unsatisfactory for proportionate relationship, causes the linear components of tissue can not phase
Disappear, remain in contrastographic picture.
Five be due to motion, and motion causes the phase of transmitted waveform several times, amplitude to be unsatisfactory for linear cancellation condition so that line
Property structural constituent eliminate stay in contrastographic picture.So, in current ultrasonic contrast imaging, it can not all avoid tissue residue
In the presence of.
The content of the invention
It is an object of the present invention to provide region detection, the imaging side of a kind of ultrasonic contrast imaging method and contrastographic picture
Method, based on nonlinear fundamental contrastographic picture, with reference to nonlinear fundamental and linear fundamental information, non-thread is distinguished in accurately identification
Property fundamental wave contrastographic picture in contrast medium areas, tissue residue and noise region, and further suppress tissue residue, strengthen radiography
Agent information, and then improve the CTR of contrastographic picture.
To achieve the above object, the present invention uses following technical scheme：
A kind of method for detecting area of contrastographic picture, including：
Tissue signal and reduction of contrast signal during acquisition ultrasonic contrast imaging；
The tissue signal and reduction of contrast signal are compared processing, the signal value for comparing the image formed after processing is obtained
Histogram；
The threshold value of contrast agent signals and tissue residue signal is determined by histogram, and by way of Threshold segmentation
The image that is formed after processing will be compared and be partitioned into the contrast agent zone of action, the tissue residue zone of action and noise contributions region.
A kind of developing method of contrastographic picture, including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and noise region using the process described above；
Contrast medium areas and tissue residue region are multiplied by after different coefficients respectively, new contrastographic picture is combined into and carries out
Display.
A kind of developing method of contrastographic picture, including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and noise region using the process described above；
Contrast medium areas, tissue residue region and noise region are multiplied by after different coefficients respectively, new make is combined into
Shadow image is shown.
A kind of developing method of contrastographic picture, including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and noise region using the process described above；
Contrast medium areas, tissue residue region, noise region and original reduction of contrast signal are multiplied by after different coefficients respectively,
New contrastographic picture is combined into be shown.
A kind of developing method of contrastographic picture, including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and noise region using the process described above；
Contrast medium areas, tissue residue region and original reduction of contrast signal are multiplied by after different coefficients respectively, are combined into new
Contrastographic picture shown.
A kind of ultrasonic contrast imaging method, including：
Ultrasonic wave is launched to tested individual by ultrasonic probe, and obtains echo-signal；
Echo-signal is handled, the linear components and non-linear component of echo-signal are extracted；
Linear components and non-linear component respectively to echo-signal carry out the post processing of multiple links, and tissue is generated respectively
Image and contrastographic picture, and show；
Region detection is carried out to contrastographic picture using the method for detecting area of above-described contrastographic picture.
A kind of ultrasonic contrast imaging method and region detection, the developing method of contrastographic picture that the present invention is provided, on the one hand
The CTR of ultrasonic contrast image is substantially increased, is conducive to clinical observation；On the other hand the signal of contrast agent is enhanced, makes ultrasound
The radiography dosage injected during radiography can be substantially reduced, and effectively reduce check cost.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the handling process schematic diagram of the first implementation of embodiment one.
Fig. 2 is the handling process schematic diagram of second of implementation of embodiment one.
Fig. 3 is the handling process schematic diagram of the third implementation of embodiment one.
Fig. 4 is the handling process schematic diagram of embodiment two.
Fig. 5 is the Nogata distribution of the signal value for the image that reduction of contrast signal and tissue signal are formed after being divided by embodiment two
Figure.
Fig. 6 is a kind of handling process schematic diagram of implementation of embodiment three.
Fig. 7 is the handling process schematic diagram of another implementation of embodiment three.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
Embodiment one
As shown in Figure 1 to Figure 3, a kind of ultrasonic contrast imaging method provided in an embodiment of the present invention, specifically includes following step
Suddenly：
S101, ultrasonic wave launched to tested individual by ultrasonic probe, and obtain echo-signal；
S102, echo-signal is handled, extract the linear components and non-linear component of echo-signal；
S103, linear components respectively to echo-signal and non-linear component carry out the post processing of multiple links, give birth to respectively
Into organization chart picture and contrastographic picture, and show；
S104, to contrastographic picture carry out region detection.
Specifically, in S101, tested individual is usually human body, and ultrasonic wave is in human body communication process by reflection, folding
Penetrate and scatter（Mainly reflect）, the echo with tissue anatomical features is propagated back to ultrasonic probe again, and ultrasonic probe is connecing
The echo received is converted into electric signal, i.e. echo-signal.
Echo-signal includes many information, it is necessary to handle echo-signal in S102, extracts useful information,
It is main including Beam synthesis with pass through the different signal component of a variety of wave filters extraction.In ultrasonic contrast imaging process
In, generally require to extract two kinds of different signals, a kind of linear components of echo-signal reflect anatomic tissue feature, also known as
For tissue signal；Another is the non-linear component of echo-signal, reflects the information of contrast agent microbubble, also known as reduction of contrast signal.
In embodiments of the present invention, the linear components are linear fundamental signal, and the non-linear component is nonlinear fundamental signal.
Linear components also need to just obtain final organization chart picture by post processing, in S103, after linear components
Processing is main to be included demodulation, seeks several links such as envelope, dynamic range transform；Likewise, non-linear component also will be through locating later
Reason can just obtain final contrastographic picture, and in S103, the post processing link of non-linear component can be similar with linear components,
The selection of parameter is different when simply handling.
In S104, linear components are compared first（Linear fundamental signal）And non-linear component（Nonlinear fundamental signal）,
The image after comparing is divided into contrast agent, three regions of tissue residue and noise followed by the method for Threshold segmentation, finally
In the corresponding position mark contrast agent of former contrastographic picture, three regions of tissue residue and noise, realize nonlinear fundamental
Contrastographic picture is divided into tissue residue, three regions of contrast agent and noise, has accurately detected contrast agent microbubble information.Wherein, institute
It both can be organization chart picture and the contrastographic picture generated after post processing to state linear components and non-linear component, as shown in Figure 1；
Can also be the M signal of a certain identical link in last handling process, as shown in Figure 2；It can also be and obtained according to echo-signal
Initial linear fundamental signal and nonlinear fundamental signal, as shown in Figure 3.
Specifically, the method that contrastographic picture carries out region detection will be described in detail in embodiment two.
As an improvement, as shown in figure 3, the embodiment of the present invention also includes：S105, S104 area detection result is shown to
On contrastographic picture, i.e., contrast agent, tissue residue and the trizonal segmentation situation of noise are shown on contrastographic picture.It is described
S105 is not the necessary step of ultrasonic contrast imaging, can be decided whether to carry out according to actual application demand.
Embodiment two
The embodiments of the invention provide a kind of method for detecting area of contrastographic picture, it is applied to nonlinear fundamental radiography figure
Picture, second harmonic contrastographic picture, and the contrastographic picture that other any nonlinear detection techniques are formed.
As shown in figure 4, a kind of method for detecting area of contrastographic picture provided in an embodiment of the present invention comprises the following steps：
Tissue signal and reduction of contrast signal during S201, acquisition ultrasonic contrast imaging；
S202, tissue signal and reduction of contrast signal are pre-processed respectively；
S203, pretreated tissue signal and reduction of contrast signal be compared to processing, obtain and compare what is formed after processing
The histogram of the signal value of image；
S204, the threshold value for determining by histogram contrast agent signals and tissue residue signal, and pass through Threshold segmentation
Mode will compare the image that is formed after processing in S203 and be partitioned into the contrast agent zone of action, the tissue residue zone of action and make an uproar
The sound zone of action；
S205, will the correspondence contrast agent zone of action on contrastographic picture, the tissue residue zone of action and noise contributions region
Position is respectively labeled as contrast medium areas, tissue residue region and noise region.
In S201, the tissue signal and reduction of contrast signal are respectively the line of the echo-signal during ultrasonic contrast imaging
Property composition and non-linear component.Specifically, the tissue signal both can be the linear of the echo during ultrasonic contrast imaging
Fundamental signal or carried out post-processing obtained organization chart picture according to linear fundamental signal, can also be last handling process
In any one link M signal；The reduction of contrast signal both can be the non-linear of the echo during ultrasonic contrast imaging
Fundamental signal or carried out post-processing obtained contrastographic picture according to nonlinear fundamental signal, can also be and post-processed
The M signal of the link of any one in journey.It should be noted that according to embodiment one, post processing is main include demodulation,
Seek several links such as envelope, dynamic range transform；Therefore, the middle letter in tissue signal and reduction of contrast signal are last handling process
Number when, both should come from post processing in identical link.
In S202, pretreatment has main steps that denoising, it is therefore an objective to recovering signal itself size, it is to avoid TGC（Time
Gain Compensate, time gain compensation）And noise brings influence on subsequent treatment.It should be noted that in correspondence S201
The tissue signal and reduction of contrast signal in different link sources, the mode of denoising is different.
In S203, the implementation for comparing processing has two kinds, and one kind is by phase in reduction of contrast signal and tissue signal
Signal with position is divided by, and can be reduction of contrast signal divided by tissue signal or tissue signal divided by reduction of contrast signal；It is another
It is that, by the signal subtraction of same position in reduction of contrast signal and tissue signal, can seem reduction of contrast signal and subtract tissue signal to plant,
Can be that tissue signal subtracts reduction of contrast signal.In the embodiment of the present invention, at the comparison being divided by using reduction of contrast signal and tissue signal
Reason mode.As shown in figure 5, the histogram of the signal value of the image formed after being divided by for reduction of contrast signal and tissue signal.By
In noise signal after the Denoising in Pretreatment in S2 be mostly zero, it is impossible to be divided by；Therefore when obtaining histogram, it will need to believe
The noise signal that number value is zero carries out quantification treatment, and finally, the signal value of noise is the signal of -60dB in Fig. 5.
From fig. 5, it can be seen that the signal distributions for comparing the image formed after processing are made up of two crests and a trough,
One of peak is larger, and another peak is smaller.Fig. 5 histogram is approximately to be made up of two normal distributions,
One is normal distribution of the less distribution center of peak value in -40dB tissue residue, during another is the larger distribution of peak value
Normal distribution of the heart in -15dB contrast agent signals.For different picture materials, distributed area can be variant, but based on super
The characteristic of sound contrast imaging, this approximate two normal distributions are characterized in inevitable.
In S204, the key of Threshold segmentation is to find suitable threshold range, and the embodiment of the present invention provides two kinds of determinations
The method of threshold range.
Wherein, first method is：Find the relatively small leak and larger peak value of histogram；Will be compared with small leak or so 2
Individual or 3 standard deviations scope as tissue residue distributed area, by the scope of 2 or 3 standard deviations of larger peak value or so
It is used as the distributed area of contrast agent；Take the upper limit signal value corresponding with the lower limit that contrast agent distribution is interval of tissue residue distributed area
Intermediate value be used as contrast agent signals and the threshold value of tissue residue signal.
Second method is：Find the relatively small leak and larger peak value of histogram；Will be compared with small leak and larger peak value
Between valley value be used as contrast agent signals and the threshold value of tissue residue signal.
After above two method threshold value, it will compare after processing that signal value is more than institute in the image that is formed in S203
The zone marker for stating threshold value is the contrast agent zone of action, and signal value is less than the threshold value and the zone marker more than noise signal is
The tissue residue zone of action, remaining zone marker is noise contributions region.Due to carrying out explanation, the signal of noise in S401
Value is the signal of -60dB in Fig. 5, therefore signal value is less than or equal to -60dB region for noise contributions region.
The signal magnitude distribution of contrast agent and tissue residue is overlapped in general contrastographic picture, it is impossible to is distinguished and is made from signal
In shadow agent and tissue residue, the embodiment of the present invention, nonlinear fundamental reduction of contrast signal and linear fundamental tissue signal are compared place
In the image of formation after reason, contrast agent signals are different with the distributed areas of tissue residue signal, can be split with suitable threshold value
It is which position of contrast agent signals is tissue residue signal to go out which position, then same stroke of progress on original contrastographic picture
Point, contrastographic picture can be just divided into three kinds of different regions of contrast agent, tissue residue and noise.Compare the mode of processing not
Together, the method for Threshold segmentation is different, and the separation of contrast agent and tissue residue can be realized by no matter being divided by or subtracting each other.
Embodiment three
The embodiment of the present invention, there is provided a kind of developing method of contrastographic picture, passes through control on the basis of embodiment two
The display weight of contrast agent, tissue residue and noise in contrastographic picture, realizes the contrast for improving contrast agent and tissue residue, improves
Contrastographic picture CTR effects, can also weaken the display of noise, improve contrastographic picture SNR (Signal to Noise Ratio, noise
Than) effect.
As shown in fig. 6, a kind of developing method of contrastographic picture provided in an embodiment of the present invention, including the use of the institute of embodiment two
Contrastographic picture is divided into contrast medium areas, tissue residue region and noise region by the method for detecting area for the contrastographic picture stated,
Also include：Contrast medium areas, tissue residue region and noise region are multiplied by different FACTOR Ps 1, P2, P3 respectively, then group again
The new contrastographic picture of synthesis is shown.
Wherein, the different values of each coefficient will bring different display effects.Coefficient is more than 1, shows that the composition is emphasised
It has been shown that, coefficient is less than 1, shows that the composition is weakened display, coefficient is equal to 0, shows that the composition is not displayed.If P1=P2=P3
=1, the effect finally shown is as original reduction of contrast signal；If P1>The contrast of P2, contrast agent and tissue residue is widened, so
The effect for improving image CTR can be reached.P1, P2 and P3 selection have many kinds, in order to reach different image effects, Ke Yiling
P1, P2 and P3 value are controlled livingly.
For the contrastographic picture presented, image CTR effect is improved, the contrast of contrast agent and tissue residue is widened, this
Inventive embodiments choose P1>1, improve the intensity of contrast agent；Choose P2<1, weaken the development of tissue residue；P3=1 is chosen, makes to make an uproar
Sound is more natural, without the entire gain for adjusting image.By the adjustment of P1 and P2 parameters, flexibly control contrastographic picture can be reached
CTR effect, allows doctor to more easily discriminate contrast agent and tissue residue, is preferably diagnosed according to contrastographic picture.
As shown in fig. 7, being a kind of improved implementation of the embodiment of the present invention.It is in the contrastographic picture shown in Fig. 6
On the basis of developing method, the display weight of original reduction of contrast signal is added, it is therefore an objective to make the contrastographic picture finally shown more
It is natural.It is also possible to neatly control P1, P2, P3 and P4 value, different image display effects are reached.When P4=0, then image
Display effect is as Fig. 6；As long as P4 ≠ 0, then tissue residue will necessarily be shown.Proved according to experiment repeatedly, in order to carry
Hi-vision CTR effect, makes contrastographic picture more naturally, preferably a set of weighted value is：The ＜ P4 of P2=P3=0, P1+P4 ＞ 1,0
＜ 1.
Finally it should be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
The present invention is described in detail with reference to the foregoing embodiments, it will be understood by those within the art that：It still may be used
To be modified to the technical scheme described in previous embodiment, or equivalent substitution is carried out to which part technical characteristic；And
These modifications are replaced, and the essence of appropriate technical solution is departed from the spirit and model of various embodiments of the present invention technical scheme
Enclose.
Claims (11)
1. a kind of method for detecting area of contrastographic picture, it is characterised in that including：
Tissue signal and reduction of contrast signal during acquisition ultrasonic contrast imaging；
The tissue signal and reduction of contrast signal are compared processing, the straight of the signal value of the image formed after processing is compared in acquisition
Square distribution map；
The threshold value of contrast agent signals and tissue residue signal is determined by histogram, and will be compared by way of Threshold segmentation
The image formed after relatively handling is partitioned into the contrast agent zone of action, the tissue residue zone of action and noise contributions region.
2. the method for detecting area of contrastographic picture according to claim 1, it is characterised in that the tissue signal and radiography
Signal is respectively the linear components and non-linear component of the echo-signal during ultrasonic contrast imaging.
3. the method for detecting area of contrastographic picture according to claim 1, it is characterised in that described by the tissue signal
The step of being compared processing with reduction of contrast signal be：By the signal subtraction of same position on radiography figure and organization chart or it is divided by.
4. the method for detecting area of contrastographic picture according to claim 1, it is characterised in that determined by histogram
The step of contrast agent signals and the threshold value of tissue residue signal is：
Find the relatively small leak and larger peak value of histogram；
Using compared with the scope of 2 or 3 standard deviations of small leak or so as tissue residue distributed area, by larger peak value or so 2
Individual or 3 standard deviations scope as contrast agent distributed area；
The intermediate value of upper limit signal value corresponding with the lower limit that contrast agent distribution is interval of tissue residue distributed area is taken as contrast agent
The threshold value of signal and tissue residue signal.
5. the method for detecting area of contrastographic picture according to claim 1, it is characterised in that determined by histogram
The step of contrast agent signals and the threshold value of tissue residue signal is：
Find the relatively small leak and larger peak value of histogram；
Contrast agent signals and the threshold value of tissue residue signal will be used as compared with the valley value between small leak and larger peak value.
6. the method for detecting area of contrastographic picture according to claim 1, it is characterised in that by the tissue signal and
Also include before the step of reduction of contrast signal is compared processing：The tissue signal and reduction of contrast signal are pre-processed respectively.
7. a kind of developing method of contrastographic picture, it is characterised in that including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and made an uproar by any described method of usage right requirement 1 to 6
Sound area domain；
Contrast medium areas and tissue residue region are multiplied by after different coefficients respectively, new contrastographic picture is combined into and is shown
Show.
8. a kind of developing method of contrastographic picture, it is characterised in that including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and made an uproar by any described method of usage right requirement 1 to 6
Sound area domain；
Contrast medium areas, tissue residue region and noise region are multiplied by after different coefficients respectively, new radiography figure is combined into
As being shown.
9. a kind of developing method of contrastographic picture, it is characterised in that including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and made an uproar by any described method of usage right requirement 1 to 6
Sound area domain；
Contrast medium areas, tissue residue region, noise region and original reduction of contrast signal are multiplied by after different coefficients respectively, combined
Cheng Xin contrastographic picture is shown.
10. a kind of developing method of contrastographic picture, it is characterised in that including：
Contrastographic picture is divided into contrast medium areas, tissue residue region and made an uproar by any described method of usage right requirement 1 to 6
Sound area domain；
Contrast medium areas, tissue residue region and original reduction of contrast signal are multiplied by after different coefficients respectively, new make is combined into
Shadow image is shown.
11. a kind of ultrasonic contrast imaging method, it is characterised in that including：
Ultrasonic wave is launched to tested individual by ultrasonic probe, and obtains echo-signal；
Echo-signal is handled, the linear components and non-linear component of echo-signal are extracted；
Linear components and non-linear component respectively to echo-signal carry out the post processing of multiple links, and organization chart picture is generated respectively
And contrastographic picture, and show；
The method for detecting area of any described contrastographic picture of usage right requirement 1 to 6 carries out region detection to contrastographic picture.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310719832.4A CN104720850B (en) | 2013-12-23 | 2013-12-23 | Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture |
| PCT/CN2014/077326 WO2015096352A1 (en) | 2013-12-23 | 2014-05-13 | Contrast-enhanced ultrasound imaging method and method for regional detection and development of imaged picture |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310719832.4A CN104720850B (en) | 2013-12-23 | 2013-12-23 | Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104720850A CN104720850A (en) | 2015-06-24 |
| CN104720850B true CN104720850B (en) | 2017-10-03 |
Family
ID=53445681
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310719832.4A Active CN104720850B (en) | 2013-12-23 | 2013-12-23 | Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN104720850B (en) |
| WO (1) | WO2015096352A1 (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018000359A1 (en) * | 2016-06-30 | 2018-01-04 | 北京深迈瑞医疗电子技术研究院有限公司 | Method and system for enhancing ultrasound contrast images and ultrasound contrast imaging device |
| CN111278363B (en) * | 2017-10-16 | 2022-07-22 | 北京深迈瑞医疗电子技术研究院有限公司 | Ultrasound imaging equipment, system and image enhancement method for contrast-enhanced ultrasound imaging |
| CN109124687B (en) * | 2018-08-30 | 2021-05-14 | 飞依诺科技(苏州)有限公司 | Ultrasonic device and medical equipment for simultaneously carrying out ultrasonic diagnosis and treatment |
| CN109498057B (en) * | 2018-12-29 | 2021-09-28 | 深圳开立生物医疗科技股份有限公司 | Ultrasonic contrast imaging method, system, control equipment and storage medium |
| CN115811958A (en) * | 2020-07-07 | 2023-03-17 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic contrast imaging method, ultrasonic imaging device and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5577505A (en) * | 1996-02-06 | 1996-11-26 | Hewlett-Packard Company | Means for increasing sensitivity in non-linear ultrasound imaging systems |
| CN1476310A (en) * | 2001-07-31 | 2004-02-18 | �ʼҷ����ֵ�������˾ | Apparatus and method of frequency compounding to perform contrast imaging |
| EP1953563A1 (en) * | 1997-10-02 | 2008-08-06 | Sunnybrook Health Science Centre | Pulse inversion doppler ultrasonic diagnostic imaging |
| CN102920477A (en) * | 2012-03-05 | 2013-02-13 | 杭州弘恩医疗科技有限公司 | Apparatus and method for determining boundary of target area in medical imaging |
| CN102933154A (en) * | 2010-06-04 | 2013-02-13 | 株式会社日立医疗器械 | Ultrasound diagnosis device |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4662380A (en) * | 1985-10-31 | 1987-05-05 | General Electric Company | Adaptive time gain compensation system for ultrasound imaging |
| EP1884195A4 (en) * | 2005-05-27 | 2011-03-30 | Hitachi Medical Corp | Ultrasonograph and ultrasonic image display method |
| CN101791231A (en) * | 2010-01-28 | 2010-08-04 | 深圳市妇幼保健院 | Image processing method for automatically judging fetal hydrocephalus from ultrasonic images |
| US9095273B2 (en) * | 2011-09-26 | 2015-08-04 | Sunnybrook Research Institute | Systems and methods for automated dynamic contrast enhancement imaging |
-
2013
- 2013-12-23 CN CN201310719832.4A patent/CN104720850B/en active Active
-
2014
- 2014-05-13 WO PCT/CN2014/077326 patent/WO2015096352A1/en active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5577505A (en) * | 1996-02-06 | 1996-11-26 | Hewlett-Packard Company | Means for increasing sensitivity in non-linear ultrasound imaging systems |
| EP1953563A1 (en) * | 1997-10-02 | 2008-08-06 | Sunnybrook Health Science Centre | Pulse inversion doppler ultrasonic diagnostic imaging |
| CN1476310A (en) * | 2001-07-31 | 2004-02-18 | �ʼҷ����ֵ�������˾ | Apparatus and method of frequency compounding to perform contrast imaging |
| CN102933154A (en) * | 2010-06-04 | 2013-02-13 | 株式会社日立医疗器械 | Ultrasound diagnosis device |
| CN102920477A (en) * | 2012-03-05 | 2013-02-13 | 杭州弘恩医疗科技有限公司 | Apparatus and method for determining boundary of target area in medical imaging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104720850A (en) | 2015-06-24 |
| WO2015096352A1 (en) | 2015-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104720850B (en) | Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture | |
| DE10224234B4 (en) | System and method for phase reversal ultrasound imaging | |
| US9002080B2 (en) | Singular value filter for imaging or detection | |
| Adam et al. | The combined effect of spatial compounding and nonlinear filtering on the speckle reduction in ultrasound images | |
| Zhao et al. | Subarray coherence based postfilter for eigenspace based minimum variance beamformer in ultrasound plane-wave imaging | |
| Park et al. | Combined chirp coded tissue harmonic and fundamental ultrasound imaging for intravascular ultrasound: 20–60 MHz phantom and ex vivo results | |
| WO2016161574A1 (en) | Ultrasound contrast imaging method and apparatus | |
| CN108135566B (en) | Method and system for enhancing contrast-enhanced ultrasound images, and contrast-enhancing ultrasound imaging equipment | |
| US12004908B2 (en) | Ultrasound imaging system for high resolution wideband harmonic imaging | |
| CN101496728A (en) | Supersonic frequency composite imaging method and device | |
| Zimbico et al. | Eigenspace generalized sidelobe canceller combined with SNR dependent coherence factor for plane wave imaging | |
| Zhao et al. | Coded excitation for diverging wave cardiac imaging: A feasibility study | |
| US11737734B2 (en) | Ultrasound imaging device and system, and image enhancement method for contrast enhanced ultrasound imaging | |
| Chang et al. | Frequency compounded imaging with a high-frequency dual element transducer | |
| Lee et al. | Characterization of limb lymphedema using the statistical analysis of ultrasound backscattering | |
| Zhang et al. | Partial Hadamard encoded synthetic transmit aperture for high frame rate imaging with minimal l 2-norm least squares method | |
| CN105078512B (en) | A method and device for generating composite ultrasonic blood flow images | |
| Akamatsu et al. | Investigation on improving performance of adaptive beamformer by statistical analysis of ultrasonic echoes | |
| Zheng et al. | Signal eigenvalue factor for synthetic transmit aperture ultrasound imaging | |
| CN110009585A (en) | A kind of multi-scale image Enhancement Method | |
| CN105675731B (en) | Array is the same as hair, the detection signal enhancing method with receipts ultrasonic probe | |
| Sun et al. | A wavelet threshold denoising method for ultrasonic signal based on EMD and correlation coefficient analysis | |
| Giannoula et al. | Estimating the local viscoelastic properties from dispersive shear waves using time–frequency ridge analysis | |
| Chen et al. | Quantitative assessment of effects of phase aberration and noise on high-frame-rate imaging | |
| TW201512686A (en) | Ultrasonic image compensation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180518 Address after: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong Co-patentee after: SHENZHEN MINDRAY SCIENTIFIC Co.,Ltd. Patentee after: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. Address before: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong Patentee before: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20150624 Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd. Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. Contract record no.: X2022440020009 Denomination of invention: Contrast-enhanced ultrasound imaging method and area detection and imaging method of contrast-enhanced image Granted publication date: 20171003 License type: Common License Record date: 20220804 |
|
| EE01 | Entry into force of recordation of patent licensing contract |