CN105212961B - A kind of acoustic radiation shear-wave velocity detection method and system - Google Patents
A kind of acoustic radiation shear-wave velocity detection method and system Download PDFInfo
- Publication number
- CN105212961B CN105212961B CN201510514501.6A CN201510514501A CN105212961B CN 105212961 B CN105212961 B CN 105212961B CN 201510514501 A CN201510514501 A CN 201510514501A CN 105212961 B CN105212961 B CN 105212961B
- Authority
- CN
- China
- Prior art keywords
- check bit
- wave
- depth
- shearing wave
- transposed matrix
- 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.)
- Expired - Fee Related
Links
Abstract
This application discloses a kind of acoustic radiation shear-wave velocity detection method and systems, including select impact position, the first check bit and the second check bit respectively, and acoustic radiation force impact is carried out to impact position；It adds up respectively to two groups of time depth shearing wave transposed matrixs along depth direction, obtains two groups of cumulative displacement arrays by adding up and forming；Maximum value search is carried out to two groups of cumulative displacement arrays respectively, using maximum value as the time at shearing wave crest location；According to the shearing wave crest location time of the distance between first check bit and second check bit and first check bit and second check bit, shear wave velocity is calculated.The application is not necessarily to find the focal position in transposed matrix, but first transposed matrix is added up along depth direction, comprehensively utilize the shearing wave wave crest information of different depth, find out with most shearing wave wave crests by time location, find out the lateral position where maximum value, you can calculate shear wave velocity.
Description
Technical field
This application involves ultrasonic imaging technique more particularly to a kind of acoustic radiation shear-wave velocity detection methods and system.
Background technology
The core of acoustic radiation force imaging is to carry out quantitative analysis to tissue hardness by capturing shearing wave.Shear-wave velocity
Estimation is the key that Young's modulus (hardness) calculating, and the basis of imaging.Seminar of Duke University is in paper《Acoustic
Radiation Force Impulse Imaging:in Vivo Demonstration of Clinical
Feasibility》In propose TOF (Time of Flight, flight time) method carry out shear wave velocity estimation.Such as Fig. 1
It is shown, after being impacted to 30 positions of impact using acoustic radiation force, two location points are chosen on tissue 10：First detection
Position and the second check bit, Δ R is denoted as by distance among them, and calculate shearing wave 20 by cross-correlation method draws 2 positions
The displacement risen.20 wave crest of shearing wave is denoted as t respectively by 2 points of time1And t2, 20 speed of shearing wave isTOF methods
Realize very simple, computational efficiency height.
But TOF methods also have following defect：TOF methods require the shearing wave depth location for first finding out most strong displacement (i.e. burnt
Point position), then flight curve is drawn, shear wave velocity is then calculated by the crest location found in flight curve.But it looks for
There is following requirement in quasi- focal position：1, the precision of computing cross-correlation must be increased, i.e., must select small window and small step progress
Row computing cross-correlation.Operation time can be greatly increased in this way, and operation Duplication and redundancy are very big.2, to radiating force signal
Signal-to-noise ratio requires height, and noise spot is otherwise easy to be mistakenly considered focus.TOF also has some shortcomings, such as focus is only utilized
Shearing wave, is assumed to be single wave by the shearing wave information of position, and the shearing wave data useless for comprehensively utilizing other depth in this way
Perfect condition.It is this to assume to receive inside body uniformly imitating with medium, it is inhomogeneous but if uneven having
It is propagated in object, shearing wave has the characteristic of more complexity, such as divides, refraction etc..The single shearing wave for calculating focal position
Error-prone.The information of focal position in transposed matrix is only utilized in TOF methods, not to the shearing wave displacement in other depth
There is good utilization.In the case where radiant force effect is good, TOF methods can accurately estimate shear wave velocity.But in reality
In the application of human body or animal, radiant force effect can be interfered very much, and the focal position in transposed matrix is not easy to find, and
When acoustic radiation force effect is weaker, displacement of tissue unobvious caused by shearing wave are difficult to find crest location accurately and carry out using TOF methods
Wave velocity estimation.
Invention content
A kind of acoustic radiation shear-wave velocity detection method of the application offer and system.
According to the application's in a first aspect, the application provides a kind of acoustic radiation shear-wave velocity detection method, feature exists
In, including：
Impact position, the first check bit and the second check bit are selected respectively, and acoustic radiation force punching is carried out to the impact position
It hits；
The data for acquiring first check bit and second check bit are respectively processed, and obtain m- depth at two groups
Degree-shearing wave transposed matrix；
Along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix add up, obtain two groups by tiring out
Sum it up the cumulative displacement array of composition；
Respectively to described in two groups add up displacement array carry out maximum value search, using maximum value as shearing wave crest location at
Time；
According to the distance between first check bit and second check bit and first check bit and described
The shearing wave crest location time of two check bits, shear wave velocity is calculated.
The above method, it is described along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix tire out
Add, including：Local data is chosen to add up.
The above method, it is described along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix add up
Before, further include：To it is described when m- depth-shearing wave transposed matrix carry out lateral interpolation.
The above method, it is described along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix add up
Before, further include：To it is described when m- depth-shearing wave transposed matrix carry out lateral filtering.
The data of the above method, the acquisition first check bit and second check bit are respectively processed, and are had
Body includes：Using cross-correlation analysis respectively at first check bit and the data of second check bit acquisition
Reason.
According to the second aspect of the application, the application provides a kind of acoustic radiation shear-wave velocity detecting system, including：
Selecting module, for selecting impact position, the first check bit and the second check bit respectively, to the impact position into
Row acoustic radiation force impacts；
Acquisition module, the data for acquiring first check bit and second check bit are respectively processed, obtain
Go out m- depth-shearing wave transposed matrix at two groups；
Accumulator module, for along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix tire out
Add, obtains two groups of cumulative displacement arrays by adding up and forming；
Searching module, for carrying out maximum value search to the displacement array that adds up described in two groups respectively, using maximum value as cutting
Cut the time at wave crest location；
Processing module, for according to the distance between first check bit and second check bit and first inspection
The shearing wave crest location time of location and second check bit, shear wave velocity is calculated.
Above system, the accumulator module are additionally operable to selection local data and add up.
Above system, the accumulator module, be additionally operable to it is described when m- depth-shearing wave transposed matrix laterally insert
Value.
Above system, the accumulator module, be additionally operable to it is described when m- depth-shearing wave transposed matrix carry out lateral filter
Wave.
Above system, the acquisition module are additionally operable to using cross-correlation analysis respectively in first check bit and institute
The data for stating the acquisition of the second check bit are handled.
As a result of above technical scheme, the advantageous effect that the application has is made to be：
In the specific implementation mode of the application, along depth direction respectively to two groups when m- depth-shearing wave transposed matrix
It adds up, obtains two groups of cumulative displacement arrays by adding up and forming, maximum value is carried out to two groups of cumulative displacement arrays respectively
Search, using maximum value as the time at shearing wave crest location, according to the distance between the first check bit and the second check bit
And first check bit and the second check bit the shearing wave crest location time, shear wave velocity is calculated.The application is without seeking
The focal position in transposed matrix is looked for, but first transposed matrix is added up along depth direction, comprehensively utilizes the shearing of different depth
Wave wave crest information, find out with most shearing wave wave crests by time location, find out the lateral position where maximum value, i.e.,
Shear wave velocity can be calculated.
Description of the drawings
Fig. 1 is shear-wave velocity estimation principle figure；
Fig. 2 is the flow chart of the acoustic radiation shear-wave velocity detection method of the application in one embodiment；
Fig. 3 is the structural schematic diagram of the acoustic radiation shear-wave velocity detecting system of the application in one embodiment.
Specific implementation mode
The application is described in further detail below by specific implementation mode combination attached drawing.
Embodiment one：
As shown in Figure 1, the acoustic radiation shear-wave velocity detection method of the application, a kind of embodiment, including following step
Suddenly：
Step 102：Impact position, the first check bit and the second check bit are selected respectively, and acoustic radiation is carried out to impact position
Power is impacted.
Step 104：The data for acquiring the first check bit and the second check bit are respectively processed, and obtain m- depth at two groups
Degree-shearing wave transposed matrix.
In one embodiment, it acquires the first check bit and the data of the second check bit is respectively processed, specifically may be used
To include：Using cross-correlation analysis respectively to being handled in the data of the first check bit and the acquisition of the second check bit.Cross-correlation
Window that need not be too small in calculating is long and step-length, reduction operation time can save system resource.
Step 106：Along depth direction respectively to two groups when m- depth-shearing wave transposed matrix add up, obtain two groups
By the cumulative displacement array for adding up and forming.
In one embodiment, along depth direction respectively to two groups when m- depth-shearing wave transposed matrix tire out
Add, may include choosing local data to add up.The selection of wherein local data position includes：If when m- depth-shearing wave
Transposed matrix size is M rows N row, centered on longitudinal center point, respectively takes the data length of M/4 to add up up and down along longitudinal direction.
Because of the position of focus in transposed matrix) it is likely near longitudinal center, the data near longitudinal center have most strong
Shearing wave displacement.
In one embodiment, along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix carry out
Can also include the steps that the m- depth of clock synchronization-shearing wave transposed matrix carries out lateral interpolation before cumulative.In another embodiment party
In formula, along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix add up before, can also include pair
M- depth-shearing wave transposed matrix is laterally filtered when described.
Step 108：Maximum value search is carried out to two groups of cumulative displacement arrays respectively, using maximum value as shearing wave wave crest position
Set the time at place.
Comprehensively utilize different depth shearing wave wave crest information, find out with most shearing wave wave crests by when meta position
It sets.This way has more statistical significance, makes final shearing wave in the case that radiant force poor signal or medium are non-uniform
Velocity estimation result is more stablized.
Step 110：According to the distance between the first check bit and the second check bit and the first check bit and the second check bit
The shearing wave crest location time, shear wave velocity is calculated.It can specifically be calculated by following formula：Its
Displacements of the middle △ R between the first check bit and the second check bit, △ t are that the wave crest of shearing wave passes through the first check bit and second
The time difference of check bit.
Since the echo data processing mode of the first check bit and the second check bit is the same, so here with the first check bit
Echo data for be described in detail：
1. the first check bit gathered data size of hypothesis is M*N, it will be appreciated that acquire longitudinal length in the first check bit
For the data of M, n times are acquired.
2. pair data matrix adjacent column carries out cross-correlation Displacement Estimation calculating, if step-length is Step, window is a length of
WinLength finally show that cross-correlation transposed matrix size is K* (N-1), wherein K=[(M-WinLength)/Step]+1.
This matrix m- depth-shearing wave transposed matrix when being.
3. the m- depth of clock synchronization-shearing wave transposed matrix carries out data accumulation along depth direction, i.e. column direction, length is obtained
For the cumulative and one-dimension array of N-1.
4. traversing maximizing in cumulative and array, and records and be designated as T1 under maximum value.
5. above method pair the second check bit data are run, T2 is obtained.According to formulaAcquire shear wave velocity
Degree, whereinΔR distances between the first check bit and the second check bit, by being determined when chosen position.
In the ideal situation, the position for carrying out radiant force impact is the focus in transposed matrix；But reality is impossible
, due to the uncertainty of the shearing wave wave source position of excitation, focus (the i.e. strongest propagation position of wave source energy of transposed matrix
Set) it is unknown, it needs to be calculated according to the transposed matrix calculated, estimates focal position.Common method is direct
The position that transposed matrix maximum value is found out along depth direction is judged as focal position, and the application is first transposed matrix edge
Depth direction is cumulative, find out maximum value lateral position, can directly calculate shear wave velocity, settle at one go.The application's
Acoustic radiation shear-wave velocity detection method comprehensively utilizes the shearing wave wave crest information of different depth, finds out with most shearing waves
Wave crest by time location.This way has more statistical significance, in radiant force poor signal or the non-uniform feelings of medium
Final shear wave velocity estimated result is set more to stablize under condition.When weaker present application addresses acoustic radiation force effect, shearing wave draws
The displacement of tissue unobvious risen are difficult to find crest location accurately and the problem of carry out wave velocity estimation using TOF methods.
Embodiment two：
As shown in figure 3, the acoustic radiation shear-wave velocity detecting system of the application, a kind of embodiment, including selection mould
Block, acquisition module and accumulator module.Selecting module is right for selecting impact position, the first check bit and the second check bit respectively
Impact position carries out acoustic radiation force impact.Acquisition module, for acquire the data of the first check bit and the second check bit respectively into
Row processing, obtains m- depth-shearing wave transposed matrix at two groups.Accumulator module, for along depth direction respectively to two groups when
M- depth-shearing wave transposed matrix adds up, and obtains two groups of cumulative displacement arrays by adding up and forming.Searching module is used
In respectively to described in two groups add up displacement array carry out maximum value search, using maximum value as at shearing wave crest location when
Between.Processing module, for according to the distance between the first check bit and the second check bit and the first check bit and the second check bit
The shearing wave crest location time, shear wave velocity is calculated.
The acoustic radiation shear-wave velocity detecting system of the application, accumulator module can be also used for choosing local data's progress
It is cumulative.The selection of wherein local data position includes：If when m- depth-shearing wave transposed matrix size be M rows N row, with longitudinal direction
Centered on central point, the data length of M/4 is respectively taken to add up up and down along longitudinal direction.Because of the position of focus in transposed matrix)
It is likely near longitudinal center, so the data near longitudinal center have strongest shearing wave displacement.
In one embodiment, accumulator module can be also used for the m- depth of clock synchronization-shearing wave transposed matrix and carry out cross
To interpolation.In another embodiment, accumulator module, can be also used for it is described when m- depth-shearing wave transposed matrix into
Row laterally filtering.
The acoustic radiation shear-wave velocity detecting system of the application, acquisition module can be also used for utilizing cross-correlation analysis point
It is other that the data acquired in the first check bit and the second check bit are handled.Need not be too small in cross-correlation calculation window length and
Step-length reduces operation time, can save system resource.
Since the echo data processing mode of the first check bit and the second check bit is the same, so here with the first check bit
Echo data for be described in detail：
1. the first check bit gathered data size of hypothesis is M*N, it will be appreciated that acquire longitudinal length in the first check bit
For the data of M, n times are acquired.
2. pair data matrix adjacent column carries out cross-correlation Displacement Estimation calculating, if step-length is Step, window is a length of
WinLength finally show that cross-correlation transposed matrix size is K* (N-1), wherein K=[(M-WinLength)/Step]+1.
This matrix m- depth-shearing wave transposed matrix when being.
3. the m- depth of clock synchronization-shearing wave transposed matrix carries out data accumulation along depth direction, i.e. column direction, length is obtained
For the cumulative and one-dimension array of N-1.
4. traversing maximizing in cumulative and array, and records and be designated as T1 under maximum value.
5. above method pair the second check bit data are run, T2 is obtained.According to formulaAcquire shear wave velocity
Degree, whereinΔR distances between the first check bit and the second check bit, by being determined when chosen position.
In the ideal situation, the position for carrying out radiant force impact is the focus in transposed matrix；But reality is impossible
, due to the uncertainty of the shearing wave wave source position of excitation, focus (the i.e. strongest propagation position of wave source energy of transposed matrix
Set) it is unknown, it needs to be calculated according to the transposed matrix calculated, estimates focal position.Common method is direct
The position that transposed matrix maximum value is found out along depth direction is judged as focal position, and the application is first transposed matrix edge
Depth direction is cumulative, find out maximum value lateral position, can directly calculate shear wave velocity, settle at one go.The application's
Acoustic radiation shear-wave velocity detection method comprehensively utilizes the shearing wave wave crest information of different depth, finds out with most shearing waves
Wave crest by time location.This way has more statistical significance, in radiant force poor signal or the non-uniform feelings of medium
Final shear wave velocity estimated result is set more to stablize under condition.When weaker present application addresses acoustic radiation force effect, shearing wave draws
The displacement of tissue unobvious risen are difficult to find crest location accurately and the problem of carry out wave velocity estimation using TOF methods.
The foregoing is a further detailed description of the present application in conjunction with specific implementation manners, and it cannot be said that this Shen
Specific implementation please is confined to these explanations.For those of ordinary skill in the art to which this application belongs, it is not taking off
Under the premise of conceiving from the application, a number of simple deductions or replacements can also be made.
Claims (8)
1. a kind of acoustic radiation shear-wave velocity detection method, which is characterized in that including：
Impact position, the first check bit and the second check bit are selected respectively, and acoustic radiation force impact is carried out to the position；
The data for acquiring first check bit and second check bit are respectively processed, and m- depth-is cut when obtaining two groups
Cut wave transposed matrix；
Along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix add up, obtain two groups by cumulative and
The cumulative displacement array of composition；
Respectively to described in two groups add up displacement array carry out maximum value search, using maximum value as at shearing wave crest location when
Between；
According to the distance between first check bit and second check bit and first check bit and second inspection
The shearing wave crest location time of location, shear wave velocity is calculated；
The data of the acquisition first check bit and second check bit are respectively processed, and are specifically included：
Using cross-correlation analysis respectively to being handled in the data of first check bit and second check bit acquisition,
In, using cross-correlation analysis respectively to being handled in the data of first check bit and second check bit acquisition, wrap
It includes and cross-correlation Displacement Estimation calculating is carried out to data matrix adjacent column, if step-length is Step, a length of WinLength of window, obtain mutually
Correlation displacement matrix size is K* (N-1), and wherein K=(M-WinLength)/Step+1, M are in the first check bit gathered data
Longitudinal length, N is times of collection, this matrix m- depth-shearing wave transposed matrix when being.
2. acoustic radiation shear-wave velocity detection method as described in claim 1, which is characterized in that described to distinguish along depth direction
M- depth-shearing wave transposed matrix adds up when to described in two groups, including：
Local data is chosen to add up.
3. acoustic radiation shear-wave velocity detection method as described in claim 1, which is characterized in that described to distinguish along depth direction
Before m- depth-shearing wave transposed matrix adds up when to described in two groups, further include：
To it is described when m- depth-shearing wave transposed matrix carry out lateral interpolation.
4. acoustic radiation shear-wave velocity detection method as described in claim 1, which is characterized in that described to distinguish along depth direction
Before m- depth-shearing wave transposed matrix adds up when to described in two groups, further include：
To it is described when m- depth-shearing wave transposed matrix carry out lateral filtering.
5. a kind of acoustic radiation shear-wave velocity detecting system, which is characterized in that including：
Selecting module, for selecting impact position, the first check bit and the second check bit respectively, to the impact position carry out sound
Radiant force impacts；
Acquisition module, the data for acquiring first check bit and second check bit are respectively processed, and obtain two
M- depth-shearing wave transposed matrix when group；
Accumulator module, for along depth direction respectively to described in two groups when m- depth-shearing wave transposed matrix add up, obtain
To two groups of cumulative displacement arrays by adding up and forming；
Searching module, for carrying out maximum value search to the displacement array that adds up described in two groups respectively, using maximum value as shearing wave
Time at crest location；
Processing module, for according to the distance between first check bit and second check bit and first check bit
With the shearing wave crest location time of second check bit, shear wave velocity is calculated；
The acquisition module is additionally operable to using cross-correlation analysis respectively to being adopted in first check bit and second check bit
The data of collection are handled, wherein using cross-correlation analysis respectively to being adopted in first check bit and second check bit
The data of collection are handled, including carry out cross-correlation Displacement Estimation calculating to data matrix adjacent column, if step-length is Step, window is long
For WinLength, show that cross-correlation transposed matrix size is K* (N-1), wherein K=(M-WinLength)/Step+1, this square
Battle array m- depth-shearing wave transposed matrix when being.
6. acoustic radiation shear-wave velocity detecting system as claimed in claim 5, which is characterized in that the accumulator module is also used
It adds up in choosing local data.
7. acoustic radiation shear-wave velocity detecting system as claimed in claim 5, which is characterized in that the accumulator module is also used
M- depth-shearing wave transposed matrix carries out lateral interpolation when to described.
8. acoustic radiation shear-wave velocity detecting system as claimed in claim 5, which is characterized in that the accumulator module is also used
M- depth-shearing wave transposed matrix carries out lateral filtering when to described.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510514501.6A CN105212961B (en) | 2015-08-20 | 2015-08-20 | A kind of acoustic radiation shear-wave velocity detection method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510514501.6A CN105212961B (en) | 2015-08-20 | 2015-08-20 | A kind of acoustic radiation shear-wave velocity detection method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105212961A CN105212961A (en) | 2016-01-06 |
| CN105212961B true CN105212961B (en) | 2018-08-31 |
Family
ID=54982505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510514501.6A Expired - Fee Related CN105212961B (en) | 2015-08-20 | 2015-08-20 | A kind of acoustic radiation shear-wave velocity detection method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105212961B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
| CN106645406B (en) * | 2016-12-02 | 2019-03-12 | 北京空间飞行器总体设计部 | The positioning system and localization method of a kind of spacecraft by impact from space debris |
| CN107616814A (en) * | 2017-08-25 | 2018-01-23 | 深圳中科乐普医疗技术有限公司 | A kind of biological tissue's shear-wave velocity measuring method and medical supersonic wave device |
| CN107753058B (en) * | 2017-11-22 | 2021-03-02 | 深圳中科乐普医疗技术有限公司 | Shear wave dynamic filtering method |
| CN109589138A (en) * | 2018-11-26 | 2019-04-09 | 深圳中科乐普医疗技术有限公司 | A kind of shear-wave velocity calculation method and elastogram equipment |
| CN109512465B (en) * | 2018-12-19 | 2021-04-06 | 深圳中科乐普医疗技术有限公司 | Acoustic radiation force bidirectional shear wave composite imaging method and device |
| CN110811689B (en) | 2019-10-31 | 2020-11-27 | 汕头市超声仪器研究所股份有限公司 | First-order estimation method for shear wave motion speed |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102551801A (en) * | 2010-10-06 | 2012-07-11 | 美国西门子医疗解决公司 | Shear wave velocity solution for medical ultrasound imaging |
| CN103269639A (en) * | 2010-12-22 | 2013-08-28 | 皇家飞利浦电子股份有限公司 | Shear wave velocity estimation using center of mass |
| CN103431874A (en) * | 2013-09-06 | 2013-12-11 | 中国科学院深圳先进技术研究院 | Method and system for estimating acoustic radiation force pulse imaging |
| CN104434216A (en) * | 2013-09-24 | 2015-03-25 | 美国西门子医疗解决公司 | Shear wave Estimation from Analytic Data |
| CN104510499A (en) * | 2013-09-30 | 2015-04-15 | 美国西门子医疗解决公司 | Shear wave detection in medical ultrasound imaging |
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
| CN104605890A (en) * | 2014-12-18 | 2015-05-13 | 深圳开立生物医疗科技股份有限公司 | Shear wave crest value waveform correction method, device and system and application thereof |
| CN104622507A (en) * | 2013-11-11 | 2015-05-20 | 中国科学院深圳先进技术研究院 | Elasticity modulus measurement method and system |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8147410B2 (en) * | 2009-03-23 | 2012-04-03 | The Hong Kong Polytechnic University | Method and apparatus for ultrasound imaging and elasticity measurement |
| US9468421B2 (en) * | 2012-02-16 | 2016-10-18 | Siemens Medical Solutions Usa, Inc. | Visualization of associated information in ultrasound shear wave imaging |
-
2015
- 2015-08-20 CN CN201510514501.6A patent/CN105212961B/en not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102551801A (en) * | 2010-10-06 | 2012-07-11 | 美国西门子医疗解决公司 | Shear wave velocity solution for medical ultrasound imaging |
| CN103269639A (en) * | 2010-12-22 | 2013-08-28 | 皇家飞利浦电子股份有限公司 | Shear wave velocity estimation using center of mass |
| CN103431874A (en) * | 2013-09-06 | 2013-12-11 | 中国科学院深圳先进技术研究院 | Method and system for estimating acoustic radiation force pulse imaging |
| CN104434216A (en) * | 2013-09-24 | 2015-03-25 | 美国西门子医疗解决公司 | Shear wave Estimation from Analytic Data |
| CN104510499A (en) * | 2013-09-30 | 2015-04-15 | 美国西门子医疗解决公司 | Shear wave detection in medical ultrasound imaging |
| CN104622507A (en) * | 2013-11-11 | 2015-05-20 | 中国科学院深圳先进技术研究院 | Elasticity modulus measurement method and system |
| CN104605890A (en) * | 2014-12-18 | 2015-05-13 | 深圳开立生物医疗科技股份有限公司 | Shear wave crest value waveform correction method, device and system and application thereof |
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
Non-Patent Citations (2)
| Title |
|---|
| J＇er＇emy Bercoff, Micka¨el Tanter, and Mathias Fink.Supersonic Shear Imaging: A New Technique for Soft Tissue Elasticity Mapping.《ieee transactions on ultrasonics, ferroelectrics, and frequency control》.2004,第51卷(第4期),第396-408页. * |
| Robust Estimation of Time-of-Flight Shear Wave Speed Using a Radon Sum Transformation;Ned C. Rouze，Mark L. Palmeri;《IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control》;20101231;第57卷(第12期);第2662-2670页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105212961A (en) | 2016-01-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105212961B (en) | A kind of acoustic radiation shear-wave velocity detection method and system | |
| CN109521422B (en) | Multi-target life detection method based on radar signals and detection radar | |
| CN104535655B (en) | A kind of ray tracing formula ultrasonic Lamb wave defect chromatography imaging method | |
| CN1696734A (en) | Determination of time-difference of arrival and angle of arrival | |
| CN101839893B (en) | Lamb wave virtual time reversal method with high spatial resolution | |
| JP2010511172A5 (en) | ||
| CN107942335A (en) | A kind of object identification method and equipment | |
| Draštík et al. | Fish hydroacoustic survey standardization: A step forward based on comparisons of methods and systems from vertical surveys of a large deep lake | |
| CN106019254A (en) | Separating and identifying method for multiple human body objects in distance direction of UWB impact biological radar | |
| CN106659471A (en) | Autocorrelation guided cross-correlation in ultrasound shear wave elastography | |
| CN105301114A (en) | Acoustic coating layer insertion loss measurement method based on multi-channel space-time inverse filtering technology | |
| CN109444896A (en) | Underwater sound array positioning system and its localization method | |
| CN104034801B (en) | Structural damage iteration focal imaging monitoring method anti-during based on synthesis | |
| KR101426226B1 (en) | Signal processing method of radar | |
| CN102764141A (en) | Elastography method, elastography system, and biological tissue displacement estimation method and biological tissue displacement estimation system in elastography | |
| CN110353729B (en) | Sound wave transit time detection method based on bidirectional long-short term memory network | |
| CN110507359B (en) | Shear wave imaging method and system | |
| CN106680823B (en) | Method and device for detecting target distance and speed by using sound pulse of sperm whale | |
| CN102636773A (en) | Single-element range ambiguity resistant method based on channel multipath characteristic | |
| CN105092022B (en) | Submarine target ranging pulse screening technique and device based on sperm whale cry | |
| CN111830484A (en) | Signal time location determination | |
| Qu et al. | Measuring the sound speed in deep-sea first sediment layer using a high-frequency submersible sub-bottom profiler: Method and sea trial application | |
| CN106108942B (en) | Parallel ARFI imaging method based on OpenCL | |
| CN104267413B (en) | Lifting Wavelet dual threshold Denoising Algorithm based on signal intensity self adaptation TABU search | |
| CN1284003C (en) | Method and system for measuring carrier to bottom velocity by correlated speed measuring sound radar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20161019 Address after: 518000 Guangdong city of Shenzhen province Nanshan District Guangdong streets Science Park Keyuan West Industrial Zone 25 Building 2 floor 205 Applicant after: Shenzhen red source Asset Management Co.,Ltd. Address before: West Industrial Zone, Nanshan District science and Technology Park Keyuan Shenzhen 518000 Guangdong 25 Building 2 floor No. 205 Applicant before: SHENZHEN HONGYING TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190410 Address after: 518000 Meilin Road North Ring Road, Futian District, Shenzhen City, Guangdong Province Patentee after: SHENZHEN JIAJUN INDUSTRY Co.,Ltd. Address before: 518000 Guangdong Province, Nanshan District, Shenzhen City, Guangdong Province, Guangdong Province, Yuehai Street Science and Technology Park, West Industrial Park, 25 buildings, 2 floors 205 Patentee before: Shenzhen red source Asset Management Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180831 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |