CN106691499A - Method of using scatterer distribution statistics to analyze ultrasonic echo signal - Google Patents
Method of using scatterer distribution statistics to analyze ultrasonic echo signal Download PDFInfo
- Publication number
- CN106691499A CN106691499A CN201510420848.4A CN201510420848A CN106691499A CN 106691499 A CN106691499 A CN 106691499A CN 201510420848 A CN201510420848 A CN 201510420848A CN 106691499 A CN106691499 A CN 106691499A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic wave
- ultrasonograph
- value
- echo signal
- those
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000009826 distribution Methods 0.000 title description 20
- 239000006185 dispersion Substances 0.000 claims abstract description 30
- 241000208340 Araliaceae Species 0.000 claims description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 4
- 235000008434 ginseng Nutrition 0.000 claims description 4
- 238000002592 echocardiography Methods 0.000 claims 1
- 238000002604 ultrasonography Methods 0.000 abstract description 15
- 230000014509 gene expression Effects 0.000 abstract description 11
- 210000000056 organ Anatomy 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 3
- 210000004204 blood vessel Anatomy 0.000 abstract description 2
- 230000017531 blood circulation Effects 0.000 abstract 1
- 230000011218 segmentation Effects 0.000 abstract 1
- 208000019425 cirrhosis of liver Diseases 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 230000008520 organization Effects 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 206010016654 Fibrosis Diseases 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004761 fibrosis Effects 0.000 description 3
- 230000002440 hepatic effect Effects 0.000 description 3
- 210000004185 liver Anatomy 0.000 description 3
- 230000011664 signaling Effects 0.000 description 3
- KRQUFUKTQHISJB-YYADALCUSA-N 2-[(E)-N-[2-(4-chlorophenoxy)propoxy]-C-propylcarbonimidoyl]-3-hydroxy-5-(thian-3-yl)cyclohex-2-en-1-one Chemical compound CCC\C(=N/OCC(C)OC1=CC=C(Cl)C=C1)C1=C(O)CC(CC1=O)C1CCCSC1 KRQUFUKTQHISJB-YYADALCUSA-N 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000000790 scattering method Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention provides a method and system for obtaining ultrasonic images of sound wave scatterers. The method comprises the following steps that an ultrasonic echo signal point of ultrasonic image data is regarded as a center, dispersion measures, location measures or any combination calculation of the two measures from signal image values of the ultrasonic echo signal points in a window range are counted and then a parameter value is obtained; segmentation of an distance is repeated, the above steps are repeated until all the ultrasonic image data calculations are completed; the parameter values are regarded as the sound wave scatterers values of the ultrasonic echo signal points; an ultrasound image is output in a common expression way of color gradation so that blood vessel wall information and internal blood flow information can be distinguished, doctors help to confirm the outline boundary of an observed organ or the relative position with an adjacent organ and the structural characteristics of tissue variation of different degrees are identified.
Description
Technical field
Present invention system relevant a kind of ultrasonograph generation system and its method, particularly a kind of acquirement sound
The ultrasonograph system and method for scattering of wave.
Background technology
Refer to Fig. 1, flow Figure 100 of Fig. 1 systems known techniques ultrasonic wave grey-tone image generation of the present invention.
Because the echo signal that ultrasonic probe is received is very faint and complicated, it is necessary to by many signals
Treatment could be shown with image.In step 110, image rf data is obtained through ultrasonic probe
(Radio Frequency Data), but echo signal is still belonging to the signal that shakes up and down, will be negative
The signal of phase change on the occasion of, and take the maximum of its waveform, this step 120 is referred to as signal demodulation, image
Rf data is demodulated to produce package (Envelope Data) ultrasound data.Due to the energy of ultrasonic wave
Amount can decay with tissue depth, so needing penetratingdepth gain compensation (Depth Gain again
Compensation), amplification is compensated according to depth；Again because the data of image depth are that reaction is being returned
In sound time length, the time gain compensation for manually adjusting also is provided in addition in general ultrasound machine
(Time Gain Compensation) is making up the deficiency of image contrast.Because electronic system device pair
Differ in effective dynamic frequency response (Dynamic Frequency Response) of input signal, mesh
The responding range (Dynamic Range) of preceding clinical ultrasound wave system system is normal in 100~more than 140dB,
Represent that the signal of maximum in echo differs 10~10,000,000 times with most weak cognizable signal, be to expand frequency
Response range, package ultrasound data forms GTG (Gray Scale) through overcompression (Compression)
Ultrasound data, this step 130 is referred to as data compression, for example logarithmic compression of common method
(Logarithmic Compression), adjust each position echo amplification proportion (that is, adjustment image
Color range is contrasted) to highlight fainter signal, in step 140, signal processor is a series of through this
Treatment, depth is converted to by echo time point according to the velocity of sound, and the echo signal is then aobvious with corresponding grey decision-making
Show ultrasonograph.
But, although it has been the skill popularized very much through ultrasonic system pick-up image offer medical judgment
Art, but same aura may be different through the quality of image that different ultrasonic systems is captured,
For example, ultrasound data is made up of spot, the presence of mottled effect, some low contrasts
Pathological characters will be unable to be read in ultrasonograph, to avoid its influence quality of image, it is many with
The method for the purpose of spot occurs is reduced also to propose in succession.For example, it is original apply to radar return it
Nakagami distributions also start to be applied to the statistical analysis of ultrasonic wave signal.
Ultrasonic wave Nakagami images are a kind of features for assessing tissue scattering-in arrangement situation
Ultrasonograph.When diameter of the ultrasonic wave incident wavelength much larger than organization internal scattering, it may occur that
The scattering phenomenon of ultrasonic wave, and random light and shade spot can occur in grey-tone image in ultrasonic wave echo signal
The bright spot of shape, here it is so-called spot.The presence of spot has obscured the image of minute structure, reduces
The contrast and resolution ratio of ultrasonograph.Therefore, general clinical ultrasound wave system system, can allow to use
Family adjusts the display condition of image, and provides the selection of image processing techniques to reduce the presence of speckle.
Although spot is filtered or smoothed by this measure favorable clinical using observation, and that is extended is another
One problem is that physiology and the pathological information of micro-structural are possible to therefore lose.This namely why
Traditional grey-tone image is only capable of qualitative display institutional framework, and the reason for more difficult reflection tissue characteristics.
Further, since ultrasonic wave is easily difficult to confirm aura according to shadow because of radio reception angle and noise
Feature and the relative position organized.Some problems that the noise of ultrasonic wave mainly occurs are distortion, one
The loss of a little physiologic informations, incorrect or unsuitable positioning, unsuitable brightness, unsuitable shape
Shape and unsuitable size, so that allowing doctor in actually diagnosis, it is impossible to provide effectively and correct
Relevant information.However, the user of ultrasonic system can adjust various parameters, its corresponding echo is strong
Because user's adjustment default or different ultrasonic waves scan instrument, itself image procossing is different for weak meeting,
So that relevant information is then often subject to omit in grey-tone image.
United States Patent (USP) US8092387 discloses a kind of for reflecting that it is super that organization internal scatters sub- distribution scenario
Acoustic imaging technology, it is, using the m parameter of Nakagami distributions, to further include correction program and utilize ash
Rank pattern or virtual color pattern are presented the display program of m parameter array.
Republic of China patent TWI320705 discloses a kind of for reflecting that organization internal scatters sub- distribution scenario
Ultrasound imaging techniques, it is using correction program to obtain optimal size windows and virtual color shows
Nakagami is distributed m parameter array to reflect the organization internal distribution scenario of scattering.
China's Mainland patent CN102379721 discloses a kind of side of the dynamic analysis sub- changes in distribution of scattering
Method, arranges in pairs or groups moving window skill to analyze the two-dimentional or three-dimensional of ultrasound data using probability density function
Son distribution and change in concentration situation are scattered in space, it is to calculate each seat using a probability density function
The statistical parameter of punctuate, statistical parameter represents speckle signal size of each coordinate points in moving window
Statistical distribution, and the system for comparing each corresponding coordinate points in the first and second ultrasound datas
Meter parameter, situation is changed with the distribution of scattering and concentration of dynamic analysis sample.
But US8092387, TWI320705 and CN102379721, applied entirely without the present invention is disclosed
Various sub- distribution statisticses amount modes of scattering are divided in ultrasonic wave echo signal with assessing tissue scattering-in
The situation of cloth and structural arrangement, the present invention can assist to provide different clinical information, especially be directed to
Vessel position, is distinguish between the information of vascular wall and internal flow, understands vessel position, can aid in doctor
Teacher is confirmed the profile border of its observation organ or the relative position with adjacent organs, contributes to doctor
The facility of teacher's distinguishing image.
The content of the invention
The invention provides the ultrasonograph method that one kind obtains sound scattering, in its step, first
A ultrasonograph data are obtained, the ultrasonograph data have a plurality of ultrasonic wave echo signal points,
It is respectively provided with a signal image values.
Then, centered on one first ultrasonic wave echo signal point, in counting the ultrasonograph data
Ultrasonic wave echo signal point all in one window ranges, and with a dispersion amount number (Measure of
Dispersion, MD), a location measures (Measure of Location, ML) or aforementioned quantities number it is any
Combination calculates those signal image values, obtains one first parameter value；It is spaced the first ultrasonic wave echo
Signaling point counts the ultrasound to a little less distance, then centered on one second ultrasonic wave echo signal point
Those signal image numbers of ultrasonic wave echo signal point all in the window ranges in ripple image data
Value, obtains another second parameter value.
Repeat with point distance to separate, in the window centered on calculated complex ultrasonic wave echo signal point
Those signal image values of all ultrasonic wave echo signal points in the range of mouthful, obtain plurality of parameters
Value, is comprising the ultrasonograph data until calculating obtains first parameter value to one n-th parameter value
In each ultrasonic wave echo signal point signal image values.
Finally, using first parameter value to one n-th parameter value as the first ultrasonic wave echo signal
Point exports a ultrasonograph, with normal to the sound scattering subnumber value of the n-th ultrasonic wave echo signal point
With color range representation show those sound scattering subnumber values.
It is with foregoing side invention also provides the ultrasonograph system that one kind obtains sound scattering
Based on method, by various sub- distribution statisticses amount modes of scattering in the signal image of ultrasonograph data
Numerical statistic, and then the situation of tissue scattering-in distribution and structural arrangement is assessed, assist difference to face
The interpretation of bed information.
A kind of ultrasonograph system for obtaining sound scattering of the present invention, comprising a receiving device,
One resolution unit and a display unit.
The receiving device, for obtaining a ultrasonograph data, the ultrasonograph data have multiple
Several ultrasonic wave echo signal points, each ultrasonic wave echo signal point has a signal image values.
The resolution unit, connects the receiving device, is used to analyze the ultrasonograph data, is with one
Centered on first ultrasonic wave echo signal point, institute in a window ranges is counted in the ultrasonograph data
The ultrasonic wave echo signal point having, and with a dispersion amount number (Measure of Dispersion, MD), one
Any combination of location measures (Measure of Location, ML) or aforementioned quantities number calculates those signal shadows
As numerical value, one first parameter value is obtained；Furthermore, it is spaced the first ultrasonic wave echo signal point at least
Point distance, centered on one second ultrasonic wave echo signal point, counting should in the ultrasonograph data
Those signal image values of ultrasonic wave echo signal point all in window ranges, obtain one second and join
Numerical value；Duplicate block is centered on the point distance, a plurality of ultrasonic wave echo signal points of statistics in the window
In the range of all ultrasonic wave echo signal point those signal image values, obtain plurality of parameters value,
Until calculating obtains first parameter value to one n-th parameter value, be comprising in the ultrasonograph data it
The signal image values of those ultrasonic wave echo signal points.
The display unit, connects the resolution unit, and first parameter value to one n-th parameter value is made respectively
It is the first ultrasonic wave echo signal point to the sound scattering subnumber value of the n-th ultrasonic wave echo signal point,
A ultrasonograph is exported, those sound scattering subnumber values are shown with the color range representation commonly used.
Brief description of the drawings
Fig. 1 systems illustrate the flow chart of known techniques gray scale ultrasound ripple image generation；
Fig. 2 systems explanation present invention obtains the system reality designed by the ultrasonograph method of sound scattering
Shi Tu；
Fig. 3 systems illustrate the operation chart of window movement in ultrasonograph of the invention.
Fig. 4 systems explanation present invention is according to three kinds of forms of different percentiles of proposition with formula four；
Fig. 5 A to Fig. 5 E systems are according to produced sound scattering with 4 first groups of parameters of formula of the present invention
Ultrasonograph and its original ultrasonograph.
Fig. 6 A to Fig. 6 E systems are according to produced sound scattering with 4 second groups of parameters of formula of the present invention
Ultrasonograph and its original ultrasonograph.
Fig. 7 A to Fig. 7 E systems are according to produced sound scattering with 3rd group of parameter of formula four of the present invention
Ultrasonograph and its original ultrasonograph.
The explanation of Fig. 8 systems is according to four kinds of forms of different percentiles of proposition with formula of the present invention six；
Fig. 9 A to Fig. 9 D systems are according to the sub ultrasonic wave shadow of produced sound scattering with formula of the present invention six
Picture, and Fig. 9 E are original ultrasonograph；
Figure 10 systems are according to the sub ultrasonograph of produced sound scattering with formula of the present invention nine；
Figure 11 systems are according to the sub ultrasonograph of produced sound scattering with formula of the present invention ten；
Figure 12 systems are according to the sub ultrasonograph of produced sound scattering with formula of the present invention 11；
And
Figure 13 systems present invention obtains an embodiment flow chart of the ultrasonograph method of sound scattering.
Specific embodiment
The invention provides one kind obtain sound scattering ultrasonograph method and its system, be by
Counted in the signal image values of ultrasonograph data by various sub- distribution statisticses amount modes of scattering, entered
And the sub situation about being distributed with structural arrangement of tissue scattering-in is assessed, assist the interpretation of different clinical information.
Refer to shown in Fig. 2 and Fig. 3, Fig. 2 systems illustrate the present invention to obtain the ultrasonograph of sound scattering
The designed system of method implements figure, and the ultrasonograph system 200 of the present invention includes a receiving device
210th, a resolution unit 220 and a display unit 230；Fig. 3 systems illustrate a window (window of the invention
Block) 330 in ultrasonograph movement operation chart.
Ultrasonograph system 200 passes through receiving device 210, and such as ultrasonic probe is carried out to a tissue
Shooting obtains a ultrasonograph data 300, wherein, ultrasonograph 300 includes a plurality of ultrasounds
Ripple echo signal point, each ultrasonic wave echo signal point has a signal image values.
Then, resolution unit 220 connects the receiving device 210, is used to analyze the ultrasonograph data
300, user obtains and a window 330 is set in ultrasonograph 300, the setting of its Aspect Ratio of window 330
Launch the N multiples of incident pulse length for ultrasonic scanner, N is the natural number of non-zero, then with super
The upper left corner of sound wave image 300 is original position, ultrasonic wave echo signal point all in statistical window 330
Those signal image values, and with a dispersion amount number (Measure of Dispersion, MD), one
Any combination of the amount of putting number (Measure of Location, ML) or aforementioned quantities number calculates those signal images
Numerical value, obtains a parameter value.
Operation of the resolution unit of the present invention 220 to ultrasonograph 300, first with one first ultrasonic wave echo
Centered on signaling point, the ultrasonic wave for counting all in the range of a window 330 in the ultrasonograph data is returned
Those signal image values of sound signaling point, obtain one first parameter value.Then, the profit of resolution unit 220
Moved with fixing point distance or on-fixed point range range mode with window 330, wherein, with on-fixed
As a example by point distance, vertical or horizontal mobile point distance is all different every time, and it is moved with fixed lateral
Point distance 310, such as 6 point distances of ultrasonic wave echo signal point, and the point that fixation is vertically moved
Distance 320, such as 36 point distances of ultrasonic wave echo signal point.Lateral or longitudinal movement one each time
The point of each self-retaining apart from when, by resolution unit 220 recalculate acquirement one new parameter value.
For example, resolution unit 220 obtain first parameter value after, window 330 be spaced this first
6 points of ultrasonic wave echo signal point are moved simultaneously apart from horizontal direction and 36 points apart from longitudinal direction
And one second ultrasonic wave echo signal point is chosen, resolution unit 220 is with the second ultrasonic wave echo signal
Centered on point, circle selects the window ranges 330, and counts the window ranges 330 in the ultrasonograph data
Those signal image values of interior all ultrasonic wave echo signal points, obtain one second parameter value.Window
Mouth 330 is repeated with point distance to separate, and obtains plurality of parameters value through resolution unit 220, directly meter
Calculation obtains first parameter value to one n-th parameter value, be comprising in the ultrasonograph data those surpass
The signal image values of sound wave echo signal point, when the point distance at ultrasonic wave echo signal point interval is more than
During one ultrasonic wave echo signal point, those acquired parameter values can be carried out an interpolating function calculating
To obtain those complete parameter values.If by taking fixing point distance as an example, each longitudinal direction 320 or horizontal 310
The point distance of movement is all.
Wherein, The location measures (Measure of Location, ML) Fan Zhi Capital material value set position it
Amount number, comprising a mode value, mono- System Meter percentiles or an average value；In a Real applies example,
Parsing Single 220 Meter of unit calculate the ultrasonic wave of a little ultrasonic wave echo Xin No. Dot As centers what The window ranges 330 of The
The a little News No. image values of The of echo Xin No. Dot, obtain the Ginseng numbers of one first mode value to one n-th mode value
Value, make As sound scattering subnumbers value (expression of S) , Department equation below one：
S=Mode (window block)
In another embodiment, resolution unit 220 calculate those ultrasonic wave echo signal points centered in this
Those signal image values of the ultrasonic wave echo signal point of window ranges 330, obtain one and count hundredths
Numerical value, statistics percentile can be made by oneself by user, and X=25,50,75 ... etc. are equation below two
Represent it：
S=Percentile X (window block)
As X=50, statistics percentile is I d median, is with the first I d median to one n-th
The parameter value of I d median, as sound scattering subnumber value (S), is further represented by below equation：
S=Percentile 50 (window block)=Median (window block)
In another embodiment, resolution unit 220 calculate those ultrasonic wave echo signal points centered in the window
Those signal image values of the ultrasonic wave echo signal point of mouth scope 330, obtain one first average value extremely
The parameter value of one n-th average value, is the expression of equation below three as sound scattering subnumber value (S)：
S=Mean (window block)
Furthermore, the dispersion amount number (Measure of Dispersion, MD) refers to data values and is dissipated from center
The amount number of cloth situation, comprising a standard deviation value or a statistics range numerical value；In one embodiment,
Resolution unit 220 calculates all super in the window ranges 330 centered on those ultrasonic wave echo signal points
Those signal image values of sound wave echo signal point, obtain the parameter value of the statistics range numerical value,
Used as sound scattering subnumber value (S), foregoing statistics range numerical value is by the first statistics percentiles
Subtract each other acquired with the second statistics percentiles, be the expression of equation below four, count hundredths
Number can be made by oneself by user：
S=(Percentile X1 (window block)-Percentile X2 (window block)
Refer to Fig. 4, Fig. 4 systems with formula of the present invention four be according to propose three kinds of different percentiles it
Form, wherein, the form can propose two kinds of different data for variable X 1 and X2.Refer to figure
5A to Fig. 5 E its be with formula of the present invention four according to first group of variable (X1=95 and X2=5) be statistics
Amount mode can understand the distribution of hepatic tissue scattering-in and structural arrangement feelings for assessing variation in various degree
Condition is in the liver fibrosis case of different times.
Refer again to Fig. 6 A to Fig. 6 E its be with formula of the present invention four be according to second group of variable (X1=95
And X2=50) it is that statistic mode can understand hepatic tissue scattering-in for assessing variation in various degree
Distribution and structural arrangement situation are in the liver fibrosis case of different times.Refer again to Fig. 7 A to Fig. 7 E its
It according to the 3rd group of variable (X1=50 and X2=5) is that statistic mode can be with that system is with formula of the present invention four
The hepatic tissue scattering-in distribution of the clear different degrees of variation of assessment and structural arrangement situation are when different
The liver fibrosis case of phase.
In another embodiment, resolution unit 220 calculate those ultrasonic wave echo signal points centered in the window
Those signal image values of all ultrasonic wave echo signal points of mouth scope 330, obtain one first standard
Deviation to one n-th standard deviation value parameter value, be equation below as sound scattering subnumber value (S)
Five expressions：
S=std (window block)
Further any combinations plural number difference of the invention dispersion amount number (Measure of Dispersion,
MD) and location measures (Measure of Location, ML), the window is calculated by the combined type of the tittle number
Those signal image values of all ultrasonic wave echo signal points of mouth scope 330, carry out a parameter value
Obtain.
System such as obtains the one first dispersion amount number and one second dispersion amount number of those signal image values, and
With being divided by after value for above-mentioned dispersion amount number, then a weighting constant (C) and first power constant (m) are given as ginseng
Numerical value, as sound scattering subnumber value (S), with S=C* (MD1/MD2)mIt is represented, wherein, C and m
It can be arbitrary constant.In one embodiment, as C=1 and m=1, resolution unit 220 is by those signal shadows
After range numerical value being counted as the first statistics range numerical value of numerical value divided by the second of those signal image values
Parameter value is obtained, is the expression of equation below six as sound scattering subnumber value (S), wherein, those
Percentiles, can be made by oneself by user：
S=(Percentile X1 (window block)-Percentile X2 (window block))/
(Percentile X3(window block)–Percentile X4(window block))
Fig. 8 is refer to, Fig. 8 systems are according to four kinds of tables of different percentiles of proposition with formula of the present invention six
Lattice, wherein, the form is directed to variable X 1, X2, and X3 and X4 includes four kinds of different data.Please join again
Examine Fig. 9 A to Fig. 9 D its be with formula of the present invention six be according to Fig. 8 forms disclosed in percentile,
Fig. 9 E are original ultrasonograph, therefore the ultrasonograph of generation sound scattering can be clearly
Pick out the Structural Characteristics of tissue.Fig. 9 A its be with formula of the present invention six be according to first group of variable
(X1=75, X2=5, X3=95 and X4=5) is the ultrasonic wave shadow that statistic mode produces sound scattering
Picture.Fig. 9 B its be with formula of the present invention six be according to second group of variable (X1=50, X2=5, X3=95 and
X4=5 the ultrasonograph of sound scattering) is produced for statistic mode.Fig. 9 C its be with formula of the present invention
Six is that the 3rd group of variable (X1=25, X2=5, X3=95 and X4=5) of foundation is statistic mode generation sound
The ultrasonograph of scattering of wave.Fig. 9 D its be with formula of the present invention six be according to the 4th group of variable
(X1=95, X2=5, X3=50 and X4=0) is the ultrasonic wave shadow that statistic mode produces sound scattering
Picture, the scattering-in distribution of assessment tissue and structural arrangement situation.
In another embodiment, as C=1 and m=1, resolution unit 220 be by those signal image values it
One statistics range numerical value divided by obtaining parameter value after the standard deviation value of those signal image values, as
Sound scattering subnumber value (S), wherein, those percentiles can be made by oneself by user, be following public affairs
The expression of formula seven：
S=(Percentile X1 (window block)-Percentile X2 (window block))/
std(window block)
Measure number, a second place amount number and in the first position for and for example obtaining those signal image values
Dispersion amount number, and after stating the difference of location measures divided by the dispersion amount number above, then to give the weighting normal
(C) and power constant (m) are counted as parameter value (S), with S=C* ((ML1-ML2)/MD)mRepresent it,
Wherein, C and m can be arbitrary constant.In one embodiment, as C=1 and m=1, resolution unit 220 is
After the mode value of those signal image values is subtracted into the average value of those signal image values, then divided by
Parameter value is obtained after one of those signal image values statistics range numerical value, as sound scattering subnumber value
(S), wherein, those percentiles can be made by oneself by user, be the expression of equation below eight：
S=(Mode (window block)-Mean (window block))/(Percentile X1 (window
block)–Percentile X2(window block))
In another embodiment, as C=1 and m=1, resolution unit 220 be by those signal image values it
After mode value subtracts the average value of those signal image values, then divided by the mark of those signal image values
Parameter value is obtained after quasi- deviation, is the expression of equation below nine as sound scattering subnumber value (S)：
S=(Mode (window block)-Mean (window block))/std (window block)
Figure 10 is refer to, Figure 10 systems are according to the super of produced sound scattering with formula of the present invention nine
Sound wave image, display unit 230 by acquirement parameter value as ultrasonic wave echo signal point sound scattering
Subnumber value, is that can consider the original ultrasonographs of Fig. 9 E in light of actual conditions, and assessment tissue scattering-in is distributed and ties
Structure arranges situation.
In another embodiment, as C=1 and m=1, resolution unit 220 be by those signal image values it
After I d median subtracts the average value of those signal image values, then divided by those signal image values it
Parameter value is obtained after standard deviation value, is the expression of equation below ten as sound scattering subnumber value (S)：
S=(Median (window block)-Mean (window block))/std (window block)
Figure 11 is refer to, Figure 11 systems are according to the super of produced sound scattering with formula of the present invention ten
Sound wave image, display unit 230 by acquirement parameter value as ultrasonic wave echo signal point sound scattering
Subnumber value, is that can consider the original ultrasonographs of Fig. 9 E in light of actual conditions, and assessment tissue scattering-in is distributed and ties
Structure arranges situation.
After and for example obtaining a dispersion amount number of those signal image values and a location measures phase divider value,
The weighting constant (C) and power constant (m) are given again as parameter value (S), with S=C* (MD/ML)m
It is represented, wherein, C and m can be arbitrary constant.In one embodiment, as C=1 and m=1, parsing
Unit 220 is that the standard deviation value of those signal image values is flat divided by those signal image values
Average, obtains parameter value, is the expression of equation below 11 as sound scattering subnumber value (S)：
S=Std (window block)/Mean (window block)
Figure 12 is refer to, Figure 12 systems are according to produced sound scattering with formula of the present invention 11
Ultrasonograph, display unit 230 is dissipated by the parameter value of acquirement as the sound wave of ultrasonic wave echo signal point
Penetrate subnumber value, be that can consider the original ultrasonographs of Fig. 9 E in light of actual conditions, assessment tissue scattering-in distribution and
Structural arrangement situation.
The present invention contrasts original ultrasonograph by the generation of the ultrasonograph of sound scattering,
The Structural Characteristics of linked groups can be clearly picked out, also for vessel position, blood vessel is distinguish between
The information of wall and internal flow, and can be confirmed its profile border for observing organ or with adjacent device
The relative position of official, there is provided different clinical information (Fig. 8 to Figure 12).In addition liver also can effectively be distinguished
The order of severity of fibrosis, it can be found that fibrosis from sound scattering tissue statistics image
The image brilliance of its corresponding sound scattering subnumber value higher is relatively low, wherein, Fig. 5 A, Fig. 6 A and
Fig. 7 A represent normal liver, and Fig. 5 B to Fig. 5 E, Fig. 6 B to Fig. 6 E and Fig. 7 B to Fig. 7 E generations
The Liver fibrosis tissue of table variation in various degree, sound scattering subnumber value is smaller represents fibrosis for this
It is higher.
Refer to Figure 13, Figure 13 its be with the present invention obtain sound scattering ultrasonograph method it
The flow chart 400 of one embodiment.
First, in step 410, it is irradiated through 210 pairs one tissues of receiving device, obtains one surpassing
Sound wave image data 300, the ultrasonograph 300 has a plurality of ultrasonic wave echo signals, each to be somebody's turn to do
Ultrasonic wave echo signal point has a signal image values.
Then, at step 420, resolution unit 220 is centered on one first ultrasonic wave echo signal point,
Those numbers of all ultrasonic wave echo signal points of a window ranges 330 in statistics ultrasonograph data
Value, is used to obtain one first parameter value.
In step 430, resolution unit 220 is counted centered on one second ultrasonic wave echo signal point
Those news of ultrasonic wave echo signal point all in the window ranges 330 in the ultrasonograph data
Number image values, are used to obtain one second parameter value, wherein, the second ultrasonic wave echo signal point is
It is spaced the first ultrasonic wave echo signal point to a little less distance；For example, obtain using window 330 with
Direction laterally or longitudinally moves, wherein, the point distance 310 of transverse shifting, such as 6 ultrasounds
The distance of ripple echo signal point, and the point distance 320 for vertically moving, such as 36 ultrasonic wave echo letters
Number point distance.
In step 440, resolution unit 220 is repeated with point distance to separate, and counts a plurality of ultrasounds
Ultrasonic wave echo signal point all in the range of the window 330 centered on ripple echo signal point those
Signal image values, obtain plurality of parameters value, until calculating obtains first parameter value to one n-th ginseng
Numerical value, is the signal image number comprising those ultrasonic wave echo signal points in the ultrasonograph data
Value, when the point distance at ultrasonic wave echo signal point interval is more than a ultrasonic wave echo signal point, meeting
Those acquired parameter values are carried out into an interpolating function calculating to obtain those complete parameter values, its
In, those parameter values are a dispersion amount number, a location measures or the aforementioned quantities of those signal image values
Any combination of number.
Finally, in step 450, display unit 230 distinguishes first parameter value to one n-th parameter value
As the sound scattering subnumber of the first ultrasonic wave echo signal point to the n-th ultrasonic wave echo signal point
Value, exports a ultrasonograph, and show with the color range representation commonly used.
Disclosed herein the ultrasonograph system and method for obtaining sound scattering, by user elder generation
A window 330 is set in ultrasonograph, with the image upper left corner as original position, further according to using it
Different formulas calculate its parameter value, and its parameter Data-Statistics can be one or more or its combination, with rear window
Mouthfuls 330 movements are fixed or point distance of on-fixed, often move one after, just calculating parameter value again,
Progressively compute repeatedly untill window 330 completes whole ultrasonograph scanning, eventually the colour to commonly use
Color range is presented sound scattering subnumber value array.
Therefore through the ultrasonograph system and method disclosed herein acquirement sound scattering, can
Apply Different Organs ultrasonograph diagnose, no matter with traditional grey-tone image or ultrasonic wave
Nakagami images or compared to or using ultrasonograph produced by logarithmic compression, can all provide more
Information clinically, especially for vessel position, is distinguish between the information of vascular wall and internal flow,
Understand vessel position, doctor can be aided in be confirmed the profile of its observation organ (not only limitation is only in liver)
Border or the relative position with adjacent organs, contribute to the facility of doctor's distinguishing image.
Above-described embodiment principle and its effect only to illustrate the invention, not limits the present invention.Therefore
The personage practised in this technology above-described embodiment is modified and change still not take off the present invention spirit.This
Listed by the claim that the interest field of invention should be as be described hereinafter.
Claims (12)
1. one kind obtains the ultrasonograph method of sound scattering, and its step is included：
A ultrasonograph data are obtained, there are the ultrasonograph data a plurality of ultrasonic wave echoes to believe
Number point, each ultrasonic wave echo signal point has a signal image values；
Centered on one first ultrasonic wave echo signal point, a window in the ultrasonograph data is counted
In the range of all ultrasonic wave echo signal point those signal image values, obtain one first parameter
Value；
Centered on one second ultrasonic wave echo signal point, the window in the ultrasonograph data is counted
In the range of all ultrasonic wave echo signal point those signal image values, obtain one second parameter
Value, wherein, the second ultrasonic wave echo signal point is to be spaced the first ultrasonic wave echo signal point extremely
A little less distance；
Repeat with point distance to separate, in the window centered on calculated complex ultrasonic wave echo signal point
Those signal image values of all ultrasonic wave echo signal points in the range of mouthful, obtain a plurality of ginsengs
Numerical value, is comprising the ultrasonograph number until obtaining first parameter value to one n-th parameter value
The signal image values of those all ultrasonic wave echo signal points in, wherein, those parameter values
It is any combination of a dispersion amount number, a location measures or aforementioned quantities number of those signal image values；
And
Using first parameter value to one n-th parameter value as the first ultrasonic wave echo signal point extremely
The sound scattering subnumber value of the n-th ultrasonic wave echo signal point, exports a ultrasonograph.
2. the ultrasonograph method of sound scattering as described in claims 1, wherein, those parameter values are
One first dispersion amount number of those signal image values and the value of being divided by of one second dispersion amount number, and take
Its first power constant is multiplied by a weighting constant.
3. the ultrasonograph method of sound scattering as described in claims 1, wherein, those parameter values are
The one first position amount number and the difference of a second place amount number of those signal image values, then divided by
The dispersion amount number, and take its power constant and be multiplied by the weighting constant.
4. the ultrasonograph method of sound scattering as described in claims 1, wherein, those parameter values are
The dispersion amount number of those signal image values and the value of being divided by of the location measures, and take its power
Constant is multiplied by the weighting constant.
5. the ultrasonograph method of sound scattering as described in claims 1, wherein, the position value is
One mode value, statistics percentiles or an average value.
6. the ultrasonograph method of sound scattering as described in claims 1, wherein, the dispersion value is
One standard deviation value or a statistics range numerical value.
7. one kind obtains the ultrasonograph system of sound scattering, and it includes:
One receiving device, is used to obtain a ultrasonograph data, and the ultrasonograph data have plural number
Individual ultrasonic wave echo signal, each ultrasonic wave echo signal point has a signal image values；
One resolution unit, connects the receiving device, is used to analyze the ultrasonograph data, is with one
Centered on 1 ultrasonic wave echo signal point, count in the ultrasonograph data in a window ranges
Those signal image values of all ultrasonic wave echo signal points, obtain one first parameter value；
Furthermore, centered on one second ultrasonic wave echo signal, counting should in the ultrasonograph data
Those signal image values of ultrasonic wave echo signal point all in window ranges, obtain one the
Two parameter values, the second ultrasonic wave echo signal point is interval the first ultrasonic wave echo signal point
To a little less distance；Repeat with point distance to separate, count a plurality of ultrasonic wave echo signals
In those signal image numbers of ultrasonic wave echo signal point all in the window ranges centered on point
Value, obtains plurality of parameters value, until calculating obtains first parameter value to one n-th parameter value,
Signal image of the system comprising those ultrasonic wave echo signal points all in the ultrasonograph data
Numerical value, wherein, those parameter value systems be a dispersion amount number of those numerical value, a location measures or
Any combination of aforementioned quantities number；And
One display unit, connects the computing unit, and first parameter value to one n-th parameter value is made respectively
It is sound scattering of the first ultrasonic wave echo signal point to the n-th ultrasonic wave echo signal point
Numerical value, exports a ultrasonograph.
8. the ultrasonograph system of sound scattering as described in claims 7, wherein, those parameter values are
One first dispersion amount number of those signal image values and the value of being divided by of one second dispersion amount number, and take
Its first power constant is multiplied by a weighting constant.
9. the ultrasonograph system of sound scattering as described in claims 7, wherein, those parameter values are
The one first position amount number and the difference of a second place amount number of those signal image values, then divided by
The dispersion amount number, and take its power constant and be multiplied by the weighting constant.
10. the ultrasonograph system of sound scattering as described in claims 7, wherein, those parameter values are
The dispersion amount number of those signal image values and the value of being divided by of the location measures, and take its power
Constant is multiplied by the weighting constant.
The ultrasonograph system of 11. sound scattering as described in claims 7, wherein, the position value is
One mode value, statistics percentiles or an average value.
The ultrasonograph system of 12. sound scattering as described in claims 7, wherein, the dispersion value is
One standard deviation value or a statistics range numerical value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510420848.4A CN106691499B (en) | 2015-07-17 | 2015-07-17 | Method for analyzing an ultrasonic echo signal using scatterer distribution statistics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510420848.4A CN106691499B (en) | 2015-07-17 | 2015-07-17 | Method for analyzing an ultrasonic echo signal using scatterer distribution statistics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106691499A true CN106691499A (en) | 2017-05-24 |
| CN106691499B CN106691499B (en) | 2020-09-18 |
Family
ID=58895047
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510420848.4A Active CN106691499B (en) | 2015-07-17 | 2015-07-17 | Method for analyzing an ultrasonic echo signal using scatterer distribution statistics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106691499B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110368023A (en) * | 2019-08-01 | 2019-10-25 | 无锡海斯凯尔医学技术有限公司 | Signal processing method, device, equipment and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TW200820946A (en) * | 2006-11-10 | 2008-05-16 | Univ Nat Taiwan | Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue |
| CN102379721A (en) * | 2010-08-30 | 2012-03-21 | 国立清华大学 | Method for dynamically analyzing scatterer distribution change and application thereof |
| US9050022B2 (en) * | 2010-06-15 | 2015-06-09 | National Tsing Hua University | Method for dynamically analyzing distribution variation of scatterers and application using the same |
| CN104739448A (en) * | 2015-04-03 | 2015-07-01 | 深圳先进技术研究院 | Ultrasonic imaging method and device |
-
2015
- 2015-07-17 CN CN201510420848.4A patent/CN106691499B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TW200820946A (en) * | 2006-11-10 | 2008-05-16 | Univ Nat Taiwan | Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue |
| US9050022B2 (en) * | 2010-06-15 | 2015-06-09 | National Tsing Hua University | Method for dynamically analyzing distribution variation of scatterers and application using the same |
| CN102379721A (en) * | 2010-08-30 | 2012-03-21 | 国立清华大学 | Method for dynamically analyzing scatterer distribution change and application thereof |
| CN104739448A (en) * | 2015-04-03 | 2015-07-01 | 深圳先进技术研究院 | Ultrasonic imaging method and device |
Non-Patent Citations (1)
| Title |
|---|
| TSUI 等: "imaging local sctterer concentrations by the Nakagami statistical model", 《ULTRASOUND IN MEDICINE AND BIOLOGY》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110368023A (en) * | 2019-08-01 | 2019-10-25 | 无锡海斯凯尔医学技术有限公司 | Signal processing method, device, equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106691499B (en) | 2020-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11633175B2 (en) | Method and ultrasound system for shear wave elasticity imaging | |
| US8092387B2 (en) | Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue | |
| US20110125016A1 (en) | Fetal rendering in medical diagnostic ultrasound | |
| US7744533B2 (en) | Ultrasonic diagnostic apparatus, image processing apparatus and image processing method | |
| Khairalseed et al. | Real-time H-scan ultrasound imaging using a Verasonics research scanner | |
| JP2020531074A (en) | Ultrasound system with deep learning network for image artifact identification and removal | |
| DE102007020862A1 (en) | User interface and method for displaying information in an ultrasound system | |
| US12471885B2 (en) | Additional diagnostic data in parametric ultrasound medical imaging | |
| JP3802462B2 (en) | Ultrasonic diagnostic equipment | |
| Bahbah et al. | Linear and nonlinear characterization of microbubbles and tissue using the Nakagami statistical model | |
| CN106691499A (en) | Method of using scatterer distribution statistics to analyze ultrasonic echo signal | |
| Loizou et al. | Introduction to ultrasound imaging and speckle noise | |
| EP3179269B1 (en) | Analysis methods of ultrasound echo signals based on statistics of scatterer distributions | |
| CN113491535B (en) | Ultrasonic diagnostic apparatus, signal processing apparatus, and recording medium | |
| JP4347860B2 (en) | Ultrasonic diagnostic equipment | |
| JP7336766B2 (en) | Ultrasonic diagnostic device, ultrasonic diagnostic method and ultrasonic diagnostic program | |
| WO2007101346A1 (en) | Ultrasound simulator and method of simulating an ultrasound examination | |
| US10105123B2 (en) | Analysis methods of ultrasound echo signals based on statistics of scatterer distributions | |
| TWI539178B (en) | Analysis methods of ultrasound echo signals based on statistics of scatterer distributions | |
| CN119423813B (en) | Ultrasonic elastography method and ultrasonic imaging equipment | |
| CN113939236B (en) | Ultrasonic imaging device and method for processing ultrasonic echo signal thereof | |
| Loizou et al. | Introduction to Speckle Noise in Ultrasound Imaging and Video | |
| CN118252528A (en) | Intravascular ultrasound image blood noise suppression method and device and electronic equipment | |
| Gomez et al. | Computerized simulation of breast ultrasound images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |