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 PDF

Info

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
Application number
CN201310719832.4A
Other languages
Chinese (zh)
Other versions
CN104720850A (en
Inventor
储霞
桑茂栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Mindray Scientific Co Ltd
Original Assignee
Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Mindray Bio Medical Electronics Co Ltd filed Critical Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority to CN201310719832.4A priority Critical patent/CN104720850B/en
Priority to PCT/CN2014/077326 priority patent/WO2015096352A1/en
Publication of CN104720850A publication Critical patent/CN104720850A/en
Application granted granted Critical
Publication of CN104720850B publication Critical patent/CN104720850B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52077Details 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/481Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • G01S7/52039Details 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

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

Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture
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.
CN201310719832.4A 2013-12-23 2013-12-23 Region detection, the developing method of a kind of ultrasonic contrast imaging method and contrastographic picture Active CN104720850B (en)

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)

* Cited by examiner, † Cited by third party
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
WO2019075621A1 (en) * 2017-10-16 2019-04-25 北京深迈瑞医疗电子技术研究院有限公司 Ultrasonic imaging device and system, and image enhancement method for ultrasonic radiography and 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 深圳迈瑞生物医疗电子股份有限公司 Ultrasound contrast imaging method, ultrasound imaging apparatus, and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
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 杭州弘恩医疗科技有限公司 Device and method for determining target region boundary of medical image
CN102933154A (en) * 2010-06-04 2013-02-13 株式会社日立医疗器械 Ultrasound diagnosis device

Family Cites Families (4)

* Cited by examiner, † Cited by third party
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
WO2006126684A1 (en) * 2005-05-27 2006-11-30 Hitachi Medical Corporation 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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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 杭州弘恩医疗科技有限公司 Device and method for determining target region boundary of medical image

Also Published As

Publication number Publication date
WO2015096352A1 (en) 2015-07-02
CN104720850A (en) 2015-06-24

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
Matrone et al. Depth-of-field enhancement in filtered-delay multiply and sum beamformed images using synthetic aperture focusing
CN104240203A (en) Medical ultrasound image denoising method based on wavelet transform and quick bilateral filtering
Park et al. Combined chirp coded tissue harmonic and fundamental ultrasound imaging for intravascular ultrasound: 20–60 MHz phantom and ex vivo results
Chang et al. Frequency compounded imaging with a high-frequency dual element transducer
CN108135566B (en) Method and system for enhancing ultrasonic contrast image and ultrasonic contrast imaging equipment
CN101496728A (en) Supersonic frequency composite imaging method and device
Zhao et al. Coded excitation for diverging wave cardiac imaging: A feasibility study
Zimbico et al. Eigenspace generalized sidelobe canceller combined with SNR dependent coherence factor for plane wave imaging
US20210077078A1 (en) Ultrasound imaging system for high resolution wideband harmonic imaging
WO2016161574A1 (en) Ultrasound contrast imaging method and apparatus
Chean et al. Mutual conversion between B-mode image and acoustic impedance image
US20200253586A1 (en) Ultrasound imaging device and system, and image enhancement method for contrast enhanced ultrasound imaging
US20040127795A1 (en) Method of ultrasonic speckle reduction using wide-band frequency compounding with tissue-generated harmonics
Sun et al. A wavelet threshold denoising method for ultrasonic signal based on EMD and correlation coefficient analysis
Dei et al. Model-based beamforming with plane wave synthesis in medical ultrasound
Chen et al. Quantitative assessment of effects of phase aberration and noise on high-frame-rate imaging
TW201512686A (en) A method of compensating ultrasound image
Akamatsu et al. Investigation on improving performance of adaptive beamformer by statistical analysis of ultrasonic echoes
CN105675731B (en) Array is the same as hair, the detection signal enhancing method with receipts ultrasonic probe
Kaur et al. Speckle noise reduction using 2-D FFT in ultrasound images
CN111513755A (en) Ultrasonic imaging equipment and ultrasonic imaging method
Chen et al. A method to locate spatial distribution of scattering centers from ultrasonic backscatter signal
Giannoula et al. Estimating the local viscoelastic properties from dispersive shear waves using time–frequency ridge analysis

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