CN107095692A - Ultrasonic imaging system, method for ultrasonic imaging and one dimension displacement scan method - Google Patents
Ultrasonic imaging system, method for ultrasonic imaging and one dimension displacement scan method Download PDFInfo
- Publication number
- CN107095692A CN107095692A CN201610094686.4A CN201610094686A CN107095692A CN 107095692 A CN107095692 A CN 107095692A CN 201610094686 A CN201610094686 A CN 201610094686A CN 107095692 A CN107095692 A CN 107095692A
- Authority
- CN
- China
- Prior art keywords
- signal
- data
- ultrasonic
- mrow
- ultrasonic imaging
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 47
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000010008 shearing Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 31
- 230000008859 change Effects 0.000 claims abstract description 26
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 238000002604 ultrasonography Methods 0.000 claims abstract description 16
- 230000000704 physical effect Effects 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 230000011664 signaling Effects 0.000 claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims description 32
- 238000001914 filtration Methods 0.000 claims description 13
- 238000004891 communication Methods 0.000 claims description 12
- 238000005314 correlation function Methods 0.000 claims description 10
- 239000000523 sample Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008520 organization Effects 0.000 claims description 4
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 210000001519 tissue Anatomy 0.000 description 21
- 230000005540 biological transmission Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 210000000577 adipose tissue Anatomy 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
Abstract
The invention discloses a kind of ultrasonic imaging system, method for ultrasonic imaging and one dimension displacement scan method, belong to ultrasonic imaging system field, designed to solve the problems such as existing apparatus Detection results are poor.Ultrasonic imaging system of the present invention includes front end displacement transducer, signal-processing board and host computer, and two-way signaling is connected between front end displacement transducer and signal-processing board, and two-way signaling is connected between signal-processing board and host computer.Method for ultrasonic imaging of the present invention uses ultrasonic echo signal as the carrier of shearing wave echo, velocity variations and Young's modulus of the shearing wave in object to be detected is obtained by handling ultrasound echo signal calculating, to know the change in physical properties of object to be detected.Ultrasonic imaging system of the present invention and method for ultrasonic imaging, one dimension displacement scan method and fat content detection method based on the imaging system, the stability of a system are good, and testing result is more accurate.
Description
Technical field
The present invention relates to ultrasonic imaging system field, more particularly to a kind of ultrasonic imaging system and based on the imaging
Method for ultrasonic imaging, one dimension displacement scan method and the fat content detection method of system.
Background technology
The principle of ultrasonic elastograph imaging is the excitation by applying an outside to object to be detected (detected tissue),
Detected tissue will produce response (such as displacement, strain, a speed under the effect of the physics law such as Elasticity and biomethanics
Distribution can produce difference), calculate the displacement for obtaining being detected tissue by Digital Signal Processing and digital image processing techniques
Change and strain variation.
Existing one-dimensional elastic imaging system does not account for the influence of transversal wave movement and the difference of fatty attenuation degree
Difference, causes the displacement and the strain value degree of accuracy of elastic calculation and stability not high.
The content of the invention
It is an object of the present invention to propose a kind of ultrasonic imaging system for detecting that signal stabilization is good.
It is another object of the present invention to propose a kind of more accurate method for ultrasonic imaging of testing result.
The present invention's further an object is that propose a kind of one dimension displacement scanning that can be fully corrected to testing result
Method.
Another object of the present invention is to propose a kind of testing result more accurately fat content detection method.
For up to this purpose, on the one hand, the present invention uses following technical scheme:
A kind of ultrasonic imaging system, including front end displacement transducer, signal-processing board and host computer, the preceding end movement
Two-way signaling is connected between sensor and signal-processing board, and two-way signaling is connected between the signal-processing board and host computer;Its
In, front end displacement transducer is used to send shearing wave signal, and the shearing wave is used for according to its speed in object to be detected
Change with Young's modulus to reflect the change in physical properties of object to be detected;Signal-processing board is used to produce ultrasonic wave and receives institute
State the control port information of host computer transmission;Host computer is used to send control port information to signal-processing board, and processing ultrasound is returned
Ripple signal simultaneously calculates velocity variations and Young's modulus of the shearing wave in object to be detected, and then knows described be detected
The change in physical properties of object.
Particularly, the signal-processing board includes being used to produce the fpga chip of ultrasonic wave and for receiving the host computer
The high speed communication protocol chip of the control port information of transmission, fpga chip and high speed communication the protocol chip two-way signaling connects
It is logical.
Further, the fpga chip includes Link Initialization module, ultrasonic wave transmitter module, sampled data receiving module
With shearing wave sending module;Wherein, Link Initialization module is used to send instruction of shaking hands to the host computer by serial ports;Ultrasound
Ripple transmitter module is used to send ultrasonic square-wave signal;Sampled data receiving module is used to pass through after the ultrasonic square-wave signal is sent
Setting time is crossed to start to gather RF data;Shearing wave sending module is used to loading shearing wave signal, simultaneously in ultrasonic square-wave signal
Ultrasonic square-wave signal after loading is sent to the front end displacement transducer.
Particularly, the high speed communication protocol chip includes high speed communication agreement serial ports.
Particularly, the host computer includes initialization module, thread loops module, data read module and data processing mould
Block, wherein, initialization module is used to, by controlling end points to send read command control word, read back Serial Port Information by controlling end points,
To verify whether device link starts success of shaking hands;Thread loops module, which is used to send by the control end points, writes startup life
Order, circulation thread is read for opening block port;Data read module is used for evidence of being fetched from described piece of port;Data processing module
For producing M-mode image information using initial data, for knowing the shearing wave in object to be detected using initial data
In velocity variations and Young's modulus information.
On the other hand, the present invention uses following technical scheme:
A kind of method for ultrasonic imaging based on above-mentioned ultrasonic imaging system, shearing is used as using ultrasonic echo signal
The carrier of ripple echo, velocity variations and Young of the shearing wave in object to be detected are obtained by handling ultrasound echo signal calculating
Modulus, to know the change in physical properties of the object to be detected.
Particularly, RF data are read and data matrix is arranged in, matched filtering is carried out to the RF data;By hibert
The envelope of signal is extracted after filtering, drafting obtains M-mode ultrasonoscopy.
Particularly, read raw RF data and be arranged in data matrix form, each column represents a scan line;To the original
Beginning RF data matrix carries out bandpass filtering, and then the respective data blocks to adjacent scanning lines carry out computing cross-correlation, and calculating is obtained
Object to be detected displacement of tissue data;The displacement of tissue data are carried out with smothing filtering and matched filtering, singular value is replaced,
Obtain revised displacement of tissue data;Calculated according to the revised displacement of tissue data and obtain tissue strain, according to institute
State tissue strain calculation and obtain velocity variations and Young's modulus of the shearing wave in object to be detected.
Another further aspect, the present invention uses following technical scheme:
A kind of one dimension displacement scan method based on above-mentioned ultrasonic imaging system, by the echographic line number evidence collected
It is divided into the data block of particular size;The cross-correlation function of the data block of correspondence position in the adjacent scan line is calculated, and is counted
Calculate the cross-correlation function of the data block of correspondence position in time adjacent scan line.
Especially, adjacent the i-th -1 scan line, i-th scan line and i+1 bar scan line cross-correlation function is calculated,
Calculated results are weighted to the average displacement bias value as corresponding data block in i-th scan line, the displacement is inclined
Shifting is worth velocity variations and Young's modulus for calculating shearing wave, and then reflects the change in physical properties of object to be detected.
Further aspect, the present invention uses following technical scheme:
A kind of fat content detection method based on above-mentioned ultrasonic imaging system, the detection method is returned by ultrasound
The attenuation coefficient of ripple calculates and obtains fat content.
It is preferred that, the detection method comprises the steps:
Step 1, set up attenuation model of the ultrasonic wave in heterogeneous tissue
Wherein, ω is the number of signal
Word frequency, d is the propagation depth of signal, and (ω is d) echo spectrum of signal to S;C (d) is the DC component of signal, and representative is spread out
Penetrate coefficient;G (ω) is that, with the Gaussian Profile frequency spectrum of center probe frequency, β is the attenuation coefficient that frequency is relied on,To decline
Subtract the phase feedback of propagation equation, (ω is d) the diffusing reflection equation of tissue to R;
Step 2, each frame RF signal adding windows are divided into several overlapping parts, take the frequency spectrum S of each lap
(ω, d), by the frequency spectrum S (ω, d) carry out logarithmic transformation,
Wherein, logc (d)
For DC component, calculated by the zero-frequency component of partitioning portion frequency spectrum;LogG (ω) is the corresponding Gauss point of probe intermediate frequency
Cloth frequency spectrum,ω0It is the digital center frequency of probe, the bandwidth of σ representation signals, the bandwidth σ of signal
Obtained by RF signal of change;
The residual spectrum of the Gaussian spectrum of step 3, subtraction signal DC component and ultrasonic pulse, is represented such as in logarithm threshold
Under,
Wherein,To calculate obtained ultrasonic attenuation system
Number, the numerical value changes with the change of frequency and emission depth;R (ω, d) be organization internal diffusing reflection equation, the equation
It is the parameter changed with frequency accidental;
Step 4, to frequency spectrum logSl(ω, d) carrying out linear fit obtains slope K, is calculated by K/ (- 2d) and obtains ultrasound
Attenuation coefficient
Ultrasonic imaging system of the present invention includes front end displacement transducer, signal-processing board and host computer, based on the imaging
The method for ultrasonic imaging of system calculates velocity variations of the shearing wave in object to be detected by handling ultrasound echo signal
With Young's modulus to reflect the change in physical properties of object to be detected, the stability of a system is good, and testing result is more accurate.
One dimension displacement scan method of the present invention is based on above-mentioned ultrasonic imaging system, changes simultaneously in the longitudinal direction of detection tissue
It will laterally change as auxiliary reference amount, testing result be modified, makes testing result more accurate.
Fat content detection method of the present invention is based on above-mentioned ultrasonic imaging system, and detection process is more simplified, testing result
It is more accurate.
Brief description of the drawings
Fig. 1 is the schematic diagram for the ultrasonic imaging system that the preferred embodiment of the present invention one is provided;
Fig. 2 is the principle schematic for the one dimension displacement scan method that the preferred embodiment of the present invention two is provided.
Embodiment
Further illustrate technical scheme below in conjunction with the accompanying drawings and by embodiment.
Preferred embodiment one:
This preferred embodiment discloses a kind of ultrasonic imaging system.As shown in figure 1, the ultrasonic imaging system is main by preceding
End movement sensor, signal-processing board and the part of host computer three composition.Using the characteristic of ultrasonic waves no-wound, pass through displacement transducer
Shearing wave signal is sent, and allows ultrasonic echo signal as the carrier of shearing wave echo, by handling ultrasound echo signal meter
Velocity variations and Young's modulus of the shearing wave in object to be detected are calculated, to reflect the change in physical properties of object to be detected.
Wherein, the single-carrier signal launched by signal-processing board is converted to vector signal by front end displacement transducer,
Motor sends shearing wave signal, while sending and receiving ultrasonic signal;Signal-processing board is main by fpga chip and high speed
Communications protocol chip is constituted, and wherein fpga chip mainly completes to interact with related peripherals, and high speed communication agreement mainly completes association
Discuss the transmission of data and fpga chip is sent instructions to by asynchronous serial port port;Host computer mainly completes elastogram algorithm
Processing and interface display and control.
FPGA module signal processing flow is as follows in signal-processing board:
1st, link initialization, after host computer completes the download of firmware and reads the information state description of firmware, by upper
Reset or hand-reset, instruction of shaking hands is sent by serial ports to host computer;
2nd, ultrasonic wave transmitting is received with sampled data, is opened and is set after the enabled instruction that host computer is sent by serial ports is received
It is standby, ultrasonic square-wave signal is sent, starts the collection of RF data after waiting for a period of time;
3rd, shearing wave is sent, and starts the transmission of shearing wave after certain scan line is collected, and sends the single carrier of low frequency
To displacement transducer, the sampled data hereafter just information containing shearing wave;
4th, the information of control access, the control end message of host computer transmission is transmitted by the serial ports of high speed communication agreement
Breath, including the reading of data path and write command (only need to read the data of data path here by block, it is not necessary to by block to
Data path writes data);The status information that slave computer is transmitted by the serial ports of high speed communication agreement is instructed to host computer;
5th, the information of data path, is stored data into storage medium by the read write command of host computer, wherein, write from memory
Recognize transmission is that reading instruction reads the data collected always.
Epigynous computer section handling process:
1st, software initialization, first initialization firmware, then host computer by control end points send write order control word, open
Open monitoring worker thread circulation;To the effect that by controlling end points to send read command control word, being read back string by controlling end points
Message breath is used for verifying whether device link starts success of shaking hands;
2nd, start triggering thread loops, startup order is write by controlling end points to send, open block port and read circulation thread;
3rd, digital independent, reads in circulation thread in block port using the block read method reading data responded to host computer;
4th, data processing, including two independent sectors, a part are to produce M-mode image information using initial data, separately
A part is to produce shear wave velocity change information and Young's modulus information using initial data.
Method for ultrasonic imaging based on the ultrasonic imaging system is to use ultrasonic echo signal to be returned as shearing wave
The carrier of ripple, velocity variations and Young mould of the shearing wave in object to be detected are obtained by handling ultrasound echo signal calculating
Amount, to know the change in physical properties of object to be detected.
Specifically, by that can draw M-mode ultrasonoscopy respectively to RF data progress processing and calculate instantaneous elasticity parameter.
In M-mode ultrasonoscopy is drawn, by reading the RF data that slave computer is gathered, data matrix is arranged in, then to RF data
Matched filtering is carried out to improve the signal to noise ratio of signal, the envelope of signal is extracted after hibert is filtered, M-mode is finally drawn and surpasses
Acoustic image.In instantaneous elasticity algorithm, raw RF data is read, data matrix form is arranged in, each column represents a scan line;
The signal to noise ratio that bandpass filtering improves data is carried out to initial data, then the data block to adjacent scanning lines carries out cross-correlation fortune
Calculate, calculate the displacement of tissue caused by being propagated by shearing wave, smothing filtering and matched filtering are carried out to displacement of tissue, and replace
Singular value, then according to revised displacement of tissue computation organization strain, finally according to tissue strain calculation shear wave velocity and
Young's modulus.
Preferred embodiment two:
This preferred embodiment discloses a kind of one dimension displacement scan method.As shown in Fig. 2 by the echographic line number collected
According to the data block for being divided into particular size, every scan line is made up of certain sampled data, and correspondence depth is object to be detected (quilt
Detect tissue) detection zone.Generally only consider that the longitudinal direction of tissue changes, and calculates adjacent in conventional one dimension displacement is calculated
The cross-correlation function of the data block of correspondence position in scan line, without considering that the horizontal of tissue changes, causes displacement to calculate knot
Fruit stability is not high.Except the cross-correlation function of the data block of correspondence position in the adjacent scan line of calculating in the present embodiment,
The data block cross-correlation function of the correspondence position in time adjacent scan line is calculated, influence of the cross directional variations to result is taken into full account.Examine
Consider and the i-th -1 adjacent scan line, i-th scan line and i+1 bar scan line are calculated in the real-time of imaging, the present embodiment
Cross-correlation function, and result of calculation is weighted to the average displacement bias value as corresponding data block in i-th scan line.
Preferred embodiment three:
This preferred embodiment discloses a kind of fat content detection method, and application includes but is not limited at liver region
Fat content is detected.
Fat content is estimated by the attenuation coefficient of ultrasonic echo, and described detection method is according to ultrasound different
The attenuation model of matter tissue:Estimation ultrasound
Attenuation coefficient in different adipose tissues.ω is the numerical frequency of signal, and d is the propagation depth of signal, and (ω d) is represented S
The echo spectrum of signal.C (d) is that the DC component of signal represents diffraction coefficient, and G (ω) is represented with the height of center probe frequency
This spread spectrum, β is the attenuation coefficient that frequency is relied on,It is the phase feedback of fading propagation equation, R (ω, d) then generation
The diffusing reflection equation of table tissue.
In order to improve the sampling density of estimated spectral, the RF signal adding windows of each frame are divided into several overlapping portions
Point, take each lap frequency spectrum S (ω, d), by frequency spectrum S (ω, d) carry out logarithmic transformation obtain it is as follows:
Frequency spectrum per part subtracts DC component logc (d) Gaussian Profile frequency spectrums corresponding with probe intermediate frequency in logarithm threshold
LogG (ω), wherein logc (d) are estimated that Gaussian Profile frequency spectrum logG (ω) can by the zero-frequency component of partitioning portion frequency spectrum
To write:ω0It is the digital center frequency of probe, the bandwidth of σ representation signals, the bandwidth σ of signal
It is also to be estimated by RF signals.Think that ultrasonic pulsative signal meets Gaussian Profile to a certain extent, pulsed ultrasonic wave
It is able to can be changed with the change of investigation depth in the centre frequency of gaussian frequency spectrum, the size of bandwidth is fixed.Pass through
The average frequency spectrum of multiple adjacent battle arrays is calculated, the spectrum width of 0.6 times of peak value of amplitude regards the estimation of signal bandwidth as.
Subtraction signal DC component and the residual spectrum of the Gaussian spectrum of ultrasonic pulse are expressed as in logarithm threshold: Represent the estimation of ultrasonic attenuation coefficient, be with frequency and
Emission depth change and change, R (ω, d) be organization internal diffusing reflection equation, be the parameter changed with frequency accidental, it is right
Frequency spectrum logSl(ω, d) carrying out linear fit obtains slope K, and the attenuation coefficient of ultrasound can be calculated by K/ (- 2d)
Note, the technical principle that above are only presently preferred embodiments of the present invention and used.Those skilled in the art can manage
Solution, the invention is not restricted to specific embodiment described here, can carry out various obvious changes for a person skilled in the art
Change, readjust and substitute without departing from protection scope of the present invention.Therefore, although the present invention is entered by above example
Go and be described in further detail, but the present invention is not limited only to above example, without departing from the inventive concept,
Other more equivalent embodiments can also be included, and the scope of the present invention is determined by scope of the appended claims.
Claims (10)
1. a kind of ultrasonic imaging system, it is characterised in that including front end displacement transducer, signal-processing board and host computer, institute
State two-way signaling between front end displacement transducer and signal-processing board to connect, two-way letter between the signal-processing board and host computer
Number connection;Wherein,
Front end displacement transducer, for sending shearing wave signal, the shearing wave is used for according to its speed in object to be detected
Degree change reflects the change in physical properties of object to be detected with Young's modulus;
Signal-processing board, for producing ultrasonic wave and receiving the control port information that the host computer is sent;
Host computer, for sending control port information to signal-processing board, handles ultrasound echo signal and calculates the shearing
Velocity variations and Young's modulus of the ripple in object to be detected, and then know the change in physical properties of the object to be detected.
2. ultrasonic imaging system according to claim 1, it is characterised in that the signal-processing board includes being used to produce
The fpga chip of ultrasonic wave and the high speed communication protocol chip for receiving the control port information that the host computer is sent, it is described
Fpga chip is connected with high speed communication protocol chip two-way signaling.
3. ultrasonic imaging system according to claim 2, it is characterised in that the fpga chip includes link initialization
Module, ultrasonic wave transmitter module, sampled data receiving module and shearing wave sending module;Wherein,
Link Initialization module, for sending instruction of shaking hands to the host computer by serial ports;
Ultrasonic wave transmitter module, for sending ultrasonic square-wave signal;
Sampled data receiving module, for starting to gather RF data by setting time after the ultrasonic square-wave signal is sent;
Shearing wave sending module, for loading shearing wave signal in ultrasonic square-wave signal and believing the ultrasonic square wave after loading
Number send to the front end displacement transducer.
4. ultrasonic imaging system according to claim 2, it is characterised in that the high speed communication protocol chip includes height
Fast communications protocol serial ports.
5. according to any described ultrasonic imaging system of Claims 1-4, it is characterised in that the host computer includes initial
Change module, thread loops module, data read module and data processing module, wherein,
Initialization module, for by controlling end points to send read command control word, being read back Serial Port Information by controlling end points, to test
Whether card device link starts success of shaking hands;
Thread loops module, startup order is write for being sent by the control end points, and circulation thread is read for opening block port;
Data read module, for evidence of being fetched from described piece of port;
Data processing module, for producing M-mode image information using initial data, for knowing described cut using initial data
Cut velocity variations and Young's modulus information of the ripple in object to be detected.
6. a kind of method for ultrasonic imaging based on the ultrasonic imaging system as described in claim 1 to 5 is any, its feature exists
In, using ultrasonic echo signal as the carrier of shearing wave echo, by handle ultrasound echo signal calculate obtain shearing wave
Velocity variations and Young's modulus in object to be detected, to know the change in physical properties of the object to be detected.
7. method for ultrasonic imaging according to claim 6, it is characterised in that read RF data and be arranged in data square
The RF data are carried out matched filtering by battle array;The envelope of signal is extracted after hibert is filtered, drafting obtains M-mode ultrasound
Image.
8. the method for ultrasonic imaging according to claim 6 or 7, it is characterised in that read raw RF data and be arranged in
Data matrix form, each column represents a scan line;Bandpass filtering is carried out to the raw RF data matrix, then swept to adjacent
The respective data blocks for retouching line carry out computing cross-correlation, and calculating obtains object to be detected displacement of tissue data;To the displacement of tissue
Data carry out smothing filtering and matched filtering, replace singular value, obtain revised displacement of tissue data;According to the amendment
Displacement of tissue data afterwards, which are calculated, obtains tissue strain, and shearing wave is obtained in object to be detected according to the tissue strain calculation
Velocity variations and Young's modulus.
9. a kind of one dimension displacement scan method based on the ultrasonic imaging system as described in claim 1 to 5 is any, its feature exists
In by the echographic line number collected according to the data block for being divided into particular size;Calculate correspondence position in the adjacent scan line
The cross-correlation function for the data block put, and calculate the cross-correlation function of the data block of correspondence position in time adjacent scan line;
Especially, adjacent the i-th -1 scan line, i-th scan line and i+1 bar scan line cross-correlation function is calculated, by institute
Obtain result of calculation and be weighted the average displacement bias value as corresponding data block in i-th scan line, the displacement bias value
Velocity variations and Young's modulus for calculating shearing wave, and then reflect the change in physical properties of object to be detected.
10. a kind of fat content detection method based on the ultrasonic imaging system as described in claim 1 to 5 is any, its feature
It is, the detection method is to calculate to obtain fat content by the attenuation coefficient of ultrasonic echo;
It is preferred that, the detection method comprises the steps:
Step 1, set up attenuation model of the ultrasonic wave in heterogeneous tissue
Wherein, ω is the numeral frequency of signal
Rate, d is the propagation depth of signal, and (ω is d) echo spectrum of signal to S;C (d) is the DC component of signal, represents diffraction system
Number;G (ω) is that, with the Gaussian Profile frequency spectrum of center probe frequency, β is the attenuation coefficient that frequency is relied on,Passed for decay
The phase feedback of equation is broadcast, (ω is d) the diffusing reflection equation of tissue to R;
Step 2, each frame RF signal adding windows are divided into several overlapping parts, take each lap frequency spectrum S (ω,
D), by the frequency spectrum S (ω, d) carry out logarithmic transformation,
<mrow>
<mi>log</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>,</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>log</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>log</mi>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<mover>
<mi>&beta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>d</mi>
<mo>|</mo>
<mfrac>
<mi>&omega;</mi>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mo>+</mo>
<mi>log</mi>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>,</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein, logc (d) is straight
Flow component, is calculated by the zero-frequency component of partitioning portion frequency spectrum;LogG (ω) is the corresponding Gaussian Profile frequency of probe intermediate frequency
Spectrum, ω 0 is the digital center frequency of probe, and the bandwidth of σ representation signals, the bandwidth σ of signal passes through
RF signal of change is obtained;
The residual spectrum of the Gaussian spectrum of step 3, subtraction signal DC component and ultrasonic pulse, is expressed as follows in logarithm threshold,
<mrow>
<msub>
<mi>logS</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>,</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mover>
<mi>&beta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>d</mi>
<mo>|</mo>
<mfrac>
<mi>&omega;</mi>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mo>+</mo>
<mi>log</mi>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>,</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein, should to calculate obtained ultrasonic attenuation coefficient
Numerical value changes with the change of frequency and emission depth;R (ω, d) be organization internal diffusing reflection equation, the equation be with frequency
The parameter that rate changes at random;
Step 4, to frequency spectrum logSl(ω, d) carrying out linear fit obtains slope K, and the decay for obtaining ultrasound is calculated by K/ (- 2d)
Coefficient
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610094686.4A CN107095692B (en) | 2016-02-19 | 2016-02-19 | Ultrasonic imaging system, ultrasonic imaging method and one-dimensional displacement scanning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610094686.4A CN107095692B (en) | 2016-02-19 | 2016-02-19 | Ultrasonic imaging system, ultrasonic imaging method and one-dimensional displacement scanning method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107095692A true CN107095692A (en) | 2017-08-29 |
CN107095692B CN107095692B (en) | 2024-02-09 |
Family
ID=59658977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610094686.4A Active CN107095692B (en) | 2016-02-19 | 2016-02-19 | Ultrasonic imaging system, ultrasonic imaging method and one-dimensional displacement scanning method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107095692B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107684439A (en) * | 2017-09-12 | 2018-02-13 | 乐普(北京)医疗器械股份有限公司 | A kind of strain figure recognition methods based on ultrasonic Transient elastography technology |
CN109044402A (en) * | 2017-09-12 | 2018-12-21 | 乐普(北京)医疗器械股份有限公司 | A kind of Diagnosis of Fatty liver system based on Ultrasonic Elasticity Imaging |
CN110988764A (en) * | 2019-12-11 | 2020-04-10 | 深圳先进技术研究院 | Tissue parameter monitoring method, device, imaging system and medium |
CN111789632A (en) * | 2019-04-04 | 2020-10-20 | 株式会社日立制作所 | Ultrasonic diagnostic apparatus, signal processing apparatus, and storage medium |
CN114376603A (en) * | 2022-01-07 | 2022-04-22 | 乐普(北京)医疗器械股份有限公司 | Two-dimensional spectrum Doppler ultrasonic cardiogram image processing method and device |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1586411A (en) * | 2004-08-06 | 2005-03-02 | 清华大学 | Two dimension complex interrelative biological tissue displacement evaluating method |
CN102283679A (en) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | Ultrasonic imaging system for elasticity measurement and method for measuring elasticity of biological tissue |
US20120108968A1 (en) * | 2010-10-27 | 2012-05-03 | Siemens Medical Solutions Usa, Inc | Tissue Density Quantification Using Shear Wave Information in Medical Ultrasound Scanning |
CN104605891A (en) * | 2014-12-31 | 2015-05-13 | 中国科学院苏州生物医学工程技术研究所 | Method for detecting transmission speed of shear wave in biological tissue, method for detecting elasticity of biological tissue and method for biological tissue elasticity imaging |
CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
CN205683094U (en) * | 2016-02-19 | 2016-11-16 | 乐普(北京)医疗器械股份有限公司 | Ultrasonic imaging system |
-
2016
- 2016-02-19 CN CN201610094686.4A patent/CN107095692B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1586411A (en) * | 2004-08-06 | 2005-03-02 | 清华大学 | Two dimension complex interrelative biological tissue displacement evaluating method |
US20120108968A1 (en) * | 2010-10-27 | 2012-05-03 | Siemens Medical Solutions Usa, Inc | Tissue Density Quantification Using Shear Wave Information in Medical Ultrasound Scanning |
CN102283679A (en) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | Ultrasonic imaging system for elasticity measurement and method for measuring elasticity of biological tissue |
CN104605891A (en) * | 2014-12-31 | 2015-05-13 | 中国科学院苏州生物医学工程技术研究所 | Method for detecting transmission speed of shear wave in biological tissue, method for detecting elasticity of biological tissue and method for biological tissue elasticity imaging |
CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
CN205683094U (en) * | 2016-02-19 | 2016-11-16 | 乐普(北京)医疗器械股份有限公司 | Ultrasonic imaging system |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107684439A (en) * | 2017-09-12 | 2018-02-13 | 乐普(北京)医疗器械股份有限公司 | A kind of strain figure recognition methods based on ultrasonic Transient elastography technology |
CN109044402A (en) * | 2017-09-12 | 2018-12-21 | 乐普(北京)医疗器械股份有限公司 | A kind of Diagnosis of Fatty liver system based on Ultrasonic Elasticity Imaging |
CN111789632A (en) * | 2019-04-04 | 2020-10-20 | 株式会社日立制作所 | Ultrasonic diagnostic apparatus, signal processing apparatus, and storage medium |
CN111789632B (en) * | 2019-04-04 | 2023-06-27 | 富士胶片医疗健康株式会社 | Ultrasonic diagnostic apparatus, signal processing apparatus, and storage medium |
CN110988764A (en) * | 2019-12-11 | 2020-04-10 | 深圳先进技术研究院 | Tissue parameter monitoring method, device, imaging system and medium |
CN110988764B (en) * | 2019-12-11 | 2021-08-31 | 深圳先进技术研究院 | Tissue parameter monitoring method, device, imaging system and medium |
CN114376603A (en) * | 2022-01-07 | 2022-04-22 | 乐普(北京)医疗器械股份有限公司 | Two-dimensional spectrum Doppler ultrasonic cardiogram image processing method and device |
CN114376603B (en) * | 2022-01-07 | 2023-11-28 | 乐普(北京)医疗器械股份有限公司 | Processing method and device for two-dimensional spectrum Doppler ultrasound cardiac image |
Also Published As
Publication number | Publication date |
---|---|
CN107095692B (en) | 2024-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107095692A (en) | Ultrasonic imaging system, method for ultrasonic imaging and one dimension displacement scan method | |
Whittaker et al. | Principles of ultrasound and measurement of intramuscular fat | |
EP0123427B1 (en) | Ultrasonic medium characterization | |
JP3432204B2 (en) | Ultrasound diagnostic equipment | |
JPS63122923A (en) | Ultrasonic thermometric apparatus | |
EP1967866A1 (en) | Ultrasonic apparatus and ultrasonic diagnostic method | |
CN101176675A (en) | Ultrasonic diagnostic apparatus, ultrasonic diagnostic method, and image processing program for ultrasonic diagnostic apparatus | |
CN105300856B (en) | Based on ultrasonic impedances spectrum to granule density and the measuring method of size | |
CN102540192B (en) | Thing mark sniffer and thing mark detection method | |
US5589209A (en) | Method for a non-destructive determination of quality parameters in fresh produce | |
CN106382898A (en) | Thickness measuring method and system for icing layer outside tunnel lining of cold region | |
CN107278264A (en) | Determine device, assay method and program | |
CN205683094U (en) | Ultrasonic imaging system | |
CN104995530A (en) | Subject information acquisition apparatus, subject information acquisition method, and program | |
CN105116371B (en) | A kind of object localization method and device based on Continuous Transmission Frequency Modulation signal | |
CN109541026A (en) | A kind of the non-linear ultrasonic detection system and detection method of reflective contact metal croop property | |
CN201035147Y (en) | Wireless ultrasound wave shoal detector | |
CN105723214A (en) | Method for the non-destructive ultrasonic testing of a part by echo analysis | |
CN107064304A (en) | A kind of fruit structure the cannot-harm-detection device and method | |
CN113109825B (en) | Long bone phased ultrasonic signal characterization and bone evaluation system based on Radon transformation | |
CN109044402A (en) | A kind of Diagnosis of Fatty liver system based on Ultrasonic Elasticity Imaging | |
EP2366997A1 (en) | Method and device for determining the structural organization of an object with ultrasounds | |
CN207318423U (en) | Electromagnetic ultrasonic transducer control device and device | |
CN114072066A (en) | Elasticity imaging method, system and computer readable storage medium | |
JPH0854379A (en) | Method and apparatus for processing image signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |