CN105212961A - 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
- CN105212961A CN105212961A CN201510514501.6A CN201510514501A CN105212961A CN 105212961 A CN105212961 A CN 105212961A CN 201510514501 A CN201510514501 A CN 201510514501A CN 105212961 A CN105212961 A CN 105212961A
- Authority
- CN
- China
- Prior art keywords
- detecting position
- shearing wave
- transposed matrix
- depth
- acoustic radiation
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
This application discloses a kind of acoustic radiation shear-wave velocity detection method and system, comprise and select impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact is carried out to impact position; Respectively the degree of depth-shearing wave transposed matrix m-when two groups is added up along depth direction, obtain two groups of cumulative displacement arrays be made up of cumulative sum; Respectively maximum value search is carried out, using the time of maximum as shearing wave crest location place to two groups of displacement arrays that add up; According to the shearing wave crest location time of the distance between described first detecting position and described second detecting position and described first detecting position and described second detecting position, calculate shear wave velocity.The application is without the need to finding the focal position in transposed matrix, but first transposed matrix is added up along depth direction, the shearing wave crest information of comprehensive utilization different depth, find out the time location there is maximum shearing wave crest passing through, obtain the lateral attitude at maximum place, can shear wave velocity be calculated.
Description
Technical field
The application relates to ultrasonic imaging technique, particularly relates to a kind of acoustic radiation shear-wave velocity detection method and system.
Background technology
The core of acoustic radiation force imaging carries out quantitative analysis by catching shearing wave to organizing hardness.It is the key that Young's modulus (hardness) calculates that shear-wave velocity is estimated, is also the basis of imaging.Seminar of Duke University proposes TOF (TimeofFlight, flight time) method and carries out shear wave velocity estimation in paper " AcousticRadiationForceImpulseImaging:inVivoDemonstration ofClinicalFeasibility ".As shown in Figure 1, utilizing after acoustic radiation force impacts impact 30 position, tissue 10 is chosen two location points: the first detecting position and the second detecting position, the distance in the middle of them is denoted as Δ R, the displacement caused 2 positions by cross-correlation method calculating shearing wave 20.Shearing wave 20 crest is designated as t respectively by the time of 2
1and t
2, shearing wave 20 speed is
tOF method realizes very simple, and computational efficiency is high.
But TOF method also has following defect: TOF method requires the shearing wave depth location (i.e. focal position) first finding out the strongest displacement, then draws flight curve, then calculate shear wave velocity by the crest location found in flight curve.But find focal position accurately and have following requirement: the precision that computing cross-correlation 1, must be increased, namely must select fenestella and small step progress row computing cross-correlation.Greatly can increase add time like this, and computing Duplication and redundancy very large.2, require high to radiant force Signal-to-Noise, otherwise be easy to noise spot is thought by mistake be focus.TOF also has some shortcomings, such as, only make use of the shearing wave information of focal position, and the shearing wave data of other degree of depth of comprehensive utilization useless, like this shearing wave is assumed to be the perfect condition of single ripple.This hypothesis is being imitated inside body can also accept with medium uniformly, if but having uneven, propagate in inhomogeneous object, shearing wave has the characteristic of more complexity, such as, divide, refraction etc.The shearing wave of single calculating focal position is easy to make mistakes.TOF method only make use of the information of focal position in transposed matrix, does not well utilize the shearing wave displacement in other degree of depth.When radiant force is effective, TOF method can estimate shear wave velocity exactly.But in the application of actual human body or animal, radiant force effect can be subject to a lot of interference, and the focal position in transposed matrix is not easy to find, and when acoustic radiation force effect is more weak, the displacement of tissue that shearing wave causes is not obvious, uses TOF method be difficult to find crest location accurately and carry out wave velocity estimation.
Summary of the invention
The application provides a kind of acoustic radiation shear-wave velocity detection method and system.
According to the first aspect of the application, the application provides a kind of acoustic radiation shear-wave velocity detection method, it is characterized in that, comprising:
Select impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact is carried out to described impact position;
The data gathering described first detecting position and described second detecting position process respectively, the m-degree of depth-shearing wave transposed matrix when drawing two groups;
Respectively the m-degree of depth-shearing wave transposed matrix of time described in two groups is added up along depth direction, obtain two groups of cumulative displacement arrays be made up of cumulative sum;
Respectively maximum value search is carried out, using the time of maximum as shearing wave crest location place to displacement array cumulative described in two groups;
According to the shearing wave crest location time of the distance between described first detecting position and described second detecting position and described first detecting position and described second detecting position, calculate shear wave velocity.
Said method, describedly adds up to the m-degree of depth-shearing wave transposed matrix of time described in two groups respectively along depth direction, comprising: choose local data and add up.
Said method, described along depth direction respectively to the m-degree of depth-shearing wave transposed matrix of time described in two groups carry out cumulative before, also comprise: lateral interpolation is carried out to the degree of depth-shearing wave transposed matrix m-time described.
Said method, described along depth direction respectively to the m-degree of depth-shearing wave transposed matrix of time described in two groups carry out cumulative before, also comprise: horizontal filtering is carried out to the degree of depth-shearing wave transposed matrix m-time described.
Said method, the data of described first detecting position of described collection and described second detecting position process respectively, specifically comprise: utilize cross-correlation analysis to process the data in described first detecting position and described second detecting position collection respectively.
According to the second aspect of the application, the application provides a kind of acoustic radiation shear-wave velocity detection system, comprising:
Selecting module, for selecting impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact being carried out to described impact position;
Acquisition module, processes respectively for the data gathering described first detecting position and described second detecting position, the m-degree of depth-shearing wave transposed matrix when drawing two groups;
Accumulator module, for adding up to the m-degree of depth-shearing wave transposed matrix of time described in two groups respectively along depth direction, obtains two groups of cumulative displacement arrays be made up of cumulative sum;
Search module, for carrying out maximum value search, using the time of maximum as shearing wave crest location place to displacement array cumulative described in two groups respectively;
Processing module, for the shearing wave crest location time according to the distance between described first detecting position and described second detecting position and described first detecting position and described second detecting position, calculates shear wave velocity.
Said system, described accumulator module, also adds up for choosing local data.
Said system, described accumulator module, also for carrying out lateral interpolation to the degree of depth-shearing wave transposed matrix m-time described.
Said system, described accumulator module, also for carrying out horizontal filtering to the degree of depth-shearing wave transposed matrix m-time described.
Said system, described acquisition module, also for utilizing cross-correlation analysis to process the data in described first detecting position and described second detecting position collection respectively.
Owing to have employed above technical scheme, the beneficial effect that the application is possessed is:
In the detailed description of the invention of the application, respectively the degree of depth-shearing wave transposed matrix m-when two groups is added up along depth direction, obtain two groups of cumulative displacement arrays be made up of cumulative sum, respectively maximum value search is carried out to two groups of displacement arrays that add up, using the time of maximum as shearing wave crest location place, according to the shearing wave crest location time of the distance between the first detecting position and the second detecting position and the first detecting position and the second detecting position, calculate shear wave velocity.The application is without the need to finding the focal position in transposed matrix, but first transposed matrix is added up along depth direction, the shearing wave crest information of comprehensive utilization different depth, find out the time location there is maximum shearing wave crest passing through, obtain the lateral attitude at maximum place, can shear wave velocity be calculated.
Accompanying drawing explanation
Fig. 1 is shear-wave velocity estimation principle figure;
Fig. 2 is the acoustic radiation shear-wave velocity detection method flow chart in one embodiment of the application;
Fig. 3 is the acoustic radiation shear-wave velocity detection system structural representation in one embodiment of the application.
Detailed description of the invention
By reference to the accompanying drawings the application is described in further detail below by detailed description of the invention.
Embodiment one:
As shown in Figure 1, the acoustic radiation shear-wave velocity detection method of the application, its a kind of embodiment, comprises the following steps:
Step 102: select impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact is carried out to impact position.
Step 104: the data gathering the first detecting position and the second detecting position process respectively, the m-degree of depth-shearing wave transposed matrix when drawing two groups.
In one embodiment, the data gathering the first detecting position and the second detecting position process respectively, specifically can comprise: utilize cross-correlation analysis to process the data in the first detecting position and the second detecting position collection respectively.Do not need too little window length and step-length in cross-correlation calculation, reduce computing consuming time, can system resource be saved.
Step 106: add up to the degree of depth-shearing wave transposed matrix m-when two groups respectively along depth direction, obtains two groups of cumulative displacement arrays be made up of cumulative sum.
In one embodiment, respectively the degree of depth-shearing wave transposed matrix m-when two groups is added up along depth direction, can comprise and choose local data and add up.Wherein choosing of local data position comprises: when setting, the m-degree of depth-shearing wave transposed matrix size is as M capable N row, and centered by longitudinal center's point, the data length longitudinally respectively getting M/4 up and down adds up.Because the position of focus in transposed matrix) probably near longitudinal center, so the data near longitudinal center have the strongest shearing wave displacement.
In one embodiment, along depth direction respectively to the m-degree of depth-shearing wave transposed matrix of time described in two groups carry out cumulative before, when can also to comprise pair, the m-degree of depth-shearing wave transposed matrix carries out the step of lateral interpolation.In another embodiment, along depth direction respectively to the m-degree of depth-shearing wave transposed matrix of time described in two groups carry out cumulative before, can also comprise and horizontal filtering is carried out to the degree of depth-shearing wave transposed matrix m-time described.
Step 108: respectively maximum value search is carried out, using the time of maximum as shearing wave crest location place to two groups of displacement arrays that add up.
The shearing wave crest information of comprehensive utilization different depth, finds out the time location having maximum shearing wave crest and pass through.This way has more statistical significance, when radiant force poor signal or medium uneven make final shear wave velocity estimated result more stable.
Step 110: according to the shearing wave crest location time of the distance between the first detecting position and the second detecting position and the first detecting position and the second detecting position, calculate shear wave velocity.Specifically calculate by following formula:
wherein △ R is the displacement between the first detecting position and the second detecting position, and △ t is the time difference of crest by the first detecting position and the second detecting position of shearing wave.
Because the first detecting position is the same with the echo data processing mode of the second detecting position, so describe in detail for the echo data of the first detecting position here:
1. suppose that the first detecting position image data size is M*N, can be regarded as and acquires at the first detecting position the data that longitudinal length is M, acquire N time.
2. pair data matrix adjacent column carries out the calculating of cross-correlation Displacement Estimation, if step-length is Step, window is long is WinLength, finally show that cross-correlation transposed matrix size is K* (N-1), wherein K=[(M-WinLength)/Step]+1.The m-degree of depth-shearing wave transposed matrix when this matrix is.
3. time pair, the m-degree of depth-shearing wave transposed matrix is along depth direction, and namely column direction carries out data accumulation, obtains the cumulative sum one-dimension array that length is N-1.
4. in cumulative sum array, travel through maximizing, and be designated as T1 under recording maximum.
5. above method is run the second detecting position data, draws T2.According to formula
try to achieve shear wave velocity, wherein
Δr is the spacing of the first detecting position and the second detecting position, determines by during chosen position.
In the ideal situation, the position of carrying out radiant force impact is the focus in transposed matrix; But reality is impossible, due to the uncertainty of shearing wave wave source position excited, the focus (the propagation position that namely wave source energy is the strongest) of transposed matrix is unknown, needing the transposed matrix according to calculating to calculate, estimating focal position.Common method directly finds out the position of transposed matrix maximum, be namely judged as focal position, and the application first adds up transposed matrix along depth direction, obtain maximum lateral attitude, directly can calculate shear wave velocity, settle at one go.The shearing wave crest information of the acoustic radiation shear-wave velocity detection method comprehensive utilization different depth of the application, finds out the time location having maximum shearing wave crest and pass through.This way has more statistical significance, when radiant force poor signal or medium uneven make final shear wave velocity estimated result more stable.Present application addresses acoustic radiation force effect more weak time, the displacement of tissue that shearing wave causes is not obvious, uses TOF method be difficult to find crest location accurately and carry out the problem of wave velocity estimation.
Embodiment two:
As shown in Figure 3, the acoustic radiation shear-wave velocity detection system of the application, its a kind of embodiment, comprises and selects module, acquisition module and accumulator module.Selecting module, for selecting impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact being carried out to impact position.Acquisition module, processes respectively for the data gathering the first detecting position and the second detecting position, the m-degree of depth-shearing wave transposed matrix when drawing two groups.Accumulator module, for adding up to the degree of depth-shearing wave transposed matrix m-when two groups respectively along depth direction, obtains two groups of cumulative displacement arrays be made up of cumulative sum.Search module, for carrying out maximum value search, using the time of maximum as shearing wave crest location place to displacement array cumulative described in two groups respectively.Processing module, for the shearing wave crest location time according to the distance between the first detecting position and the second detecting position and the first detecting position and the second detecting position, calculates shear wave velocity.
The acoustic radiation shear-wave velocity detection system of the application, accumulator module, can also be used for choosing local data and add up.Wherein choosing of local data position comprises: when setting, the m-degree of depth-shearing wave transposed matrix size is as M capable N row, and centered by longitudinal center's point, the data length longitudinally respectively getting M/4 up and down adds up.Because the position of focus in transposed matrix) probably near longitudinal center, so the data near longitudinal center have the strongest shearing wave displacement.
In one embodiment, accumulator module, when can also to be used for pair, the m-degree of depth-shearing wave transposed matrix carries out lateral interpolation.In another embodiment, accumulator module, can also be used for carrying out horizontal filtering to the degree of depth-shearing wave transposed matrix m-time described.
The acoustic radiation shear-wave velocity detection system of the application, acquisition module, can also be used for utilizing cross-correlation analysis to process the data in the first detecting position and the second detecting position collection respectively.Do not need too little window length and step-length in cross-correlation calculation, reduce computing consuming time, can system resource be saved.
Because the first detecting position is the same with the echo data processing mode of the second detecting position, so describe in detail for the echo data of the first detecting position here:
1. suppose that the first detecting position image data size is M*N, can be regarded as and acquires at the first detecting position the data that longitudinal length is M, acquire N time.
2. pair data matrix adjacent column carries out the calculating of cross-correlation Displacement Estimation, if step-length is Step, window is long is WinLength, finally show that cross-correlation transposed matrix size is K* (N-1), wherein K=[(M-WinLength)/Step]+1.The m-degree of depth-shearing wave transposed matrix when this matrix is.
3. time pair, the m-degree of depth-shearing wave transposed matrix is along depth direction, and namely column direction carries out data accumulation, obtains the cumulative sum one-dimension array that length is N-1.
4. in cumulative sum array, travel through maximizing, and be designated as T1 under recording maximum.
5. above method is run the second detecting position data, draws T2.According to formula
try to achieve shear wave velocity, wherein
Δr is the spacing of the first detecting position and the second detecting position, determines by during chosen position.
In the ideal situation, the position of carrying out radiant force impact is the focus in transposed matrix; But reality is impossible, due to the uncertainty of shearing wave wave source position excited, the focus (the propagation position that namely wave source energy is the strongest) of transposed matrix is unknown, needing the transposed matrix according to calculating to calculate, estimating focal position.Common method directly finds out the position of transposed matrix maximum, be namely judged as focal position, and the application first adds up transposed matrix along depth direction, obtain maximum lateral attitude, directly can calculate shear wave velocity, settle at one go.The shearing wave crest information of the acoustic radiation shear-wave velocity detection method comprehensive utilization different depth of the application, finds out the time location having maximum shearing wave crest and pass through.This way has more statistical significance, when radiant force poor signal or medium uneven make final shear wave velocity estimated result more stable.Present application addresses acoustic radiation force effect more weak time, the displacement of tissue that shearing wave causes is not obvious, uses TOF method be difficult to find crest location accurately and carry out the problem of wave velocity estimation.
Above content is the further description done the application in conjunction with concrete embodiment, can not assert that the concrete enforcement of the application is confined to these explanations.For the application person of an ordinary skill in the technical field, under the prerequisite not departing from the application's design, some simple deduction or replace can also be made.
Claims (10)
1. an acoustic radiation shear-wave velocity detection method, is characterized in that, comprising:
Select impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact is carried out to described position;
The data gathering described first detecting position and described second detecting position process respectively, the m-degree of depth-shearing wave transposed matrix when drawing two groups;
Respectively the m-degree of depth-shearing wave transposed matrix of time described in two groups is added up along depth direction, obtain two groups of cumulative displacement arrays be made up of cumulative sum;
Respectively maximum value search is carried out, using the time of maximum as shearing wave crest location place to displacement array cumulative described in two groups;
According to the shearing wave crest location time of the distance between described first detecting position and described second detecting position and described first detecting position and described second detecting position, calculate shear wave velocity.
2. acoustic radiation shear-wave velocity detection method as claimed in claim 1, is characterized in that, describedly adds up to the m-degree of depth-shearing wave transposed matrix of time described in two groups respectively along depth direction, comprising:
Choose local data to add up.
3. acoustic radiation shear-wave velocity detection method as claimed in claim 1, is characterized in that, described along depth direction respectively to the m-degree of depth-shearing wave transposed matrix of time described in two groups carry out cumulative before, also comprise:
Lateral interpolation is carried out to the degree of depth-shearing wave transposed matrix m-time described.
4. acoustic radiation shear-wave velocity detection method as claimed in claim 1, is characterized in that, described along depth direction respectively to the m-degree of depth-shearing wave transposed matrix of time described in two groups carry out cumulative before, also comprise:
Horizontal filtering is carried out to the degree of depth-shearing wave transposed matrix m-time described.
5. the acoustic radiation shear-wave velocity detection method according to any one of Claims 1-4, is characterized in that, the data of described first detecting position of described collection and described second detecting position process respectively, specifically comprise:
Cross-correlation analysis is utilized to process the data in described first detecting position and described second detecting position collection respectively.
6. an acoustic radiation shear-wave velocity detection system, is characterized in that, comprising:
Selecting module, for selecting impact position, the first detecting position and the second detecting position respectively, acoustic radiation force impact being carried out to described impact position;
Acquisition module, processes respectively for the data gathering described first detecting position and described second detecting position, the m-degree of depth-shearing wave transposed matrix when drawing two groups;
Accumulator module, for adding up to the m-degree of depth-shearing wave transposed matrix of time described in two groups respectively along depth direction, obtains two groups of cumulative displacement arrays be made up of cumulative sum;
Search module, for carrying out maximum value search, using the time of maximum as shearing wave crest location place to displacement array cumulative described in two groups respectively;
Processing module, for the shearing wave crest location time according to the distance between described first detecting position and described second detecting position and described first detecting position and described second detecting position, calculates shear wave velocity.
7. acoustic radiation shear-wave velocity detection system as claimed in claim 6, is characterized in that, described accumulator module, also adding up for choosing local data.
8. acoustic radiation shear-wave velocity detection system as claimed in claim 6, is characterized in that, described accumulator module, also for carrying out lateral interpolation to the degree of depth-shearing wave transposed matrix m-time described.
9. acoustic radiation shear-wave velocity detection system as claimed in claim 6, is characterized in that, described accumulator module, also for carrying out horizontal filtering to the degree of depth-shearing wave transposed matrix m-time described.
10. the acoustic radiation shear-wave velocity detection system according to any one of claim 6 to 9, is characterized in that, described acquisition module, also for utilizing cross-correlation analysis to process the data in described first detecting position and described second detecting position collection respectively.
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 true CN105212961A (en) | 2016-01-06 |
| CN105212961B 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) |
Cited By (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 |
| CN106645406A (en) * | 2016-12-02 | 2017-05-10 | 北京空间飞行器总体设计部 | Positioning system and method of spacecraft subjected to space junk collision |
| CN107616814A (en) * | 2017-08-25 | 2018-01-23 | 深圳中科乐普医疗技术有限公司 | A kind of biological tissue's shear-wave velocity measuring method and medical supersonic wave device |
| CN107753058A (en) * | 2017-11-22 | 2018-03-06 | 深圳中科乐普医疗技术有限公司 | A kind of shearing wave dynamic filter method |
| CN109512465A (en) * | 2018-12-19 | 2019-03-26 | 深圳中科乐普医疗技术有限公司 | A kind of acoustic radiation force double direction shear wave composite imaging method and device |
| CN109589138A (en) * | 2018-11-26 | 2019-04-09 | 深圳中科乐普医疗技术有限公司 | A kind of shear-wave velocity calculation method and elastogram equipment |
| WO2021082043A1 (en) * | 2019-10-31 | 2021-05-06 | 汕头市超声仪器研究所股份有限公司 | First-order estimation method for shear wave motion velocity |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100240994A1 (en) * | 2009-03-23 | 2010-09-23 | The Hong Kong Polytechnic University | Method and apparatus for ultrasound imaging and elasticity measurement |
| CN102551801A (en) * | 2010-10-06 | 2012-07-11 | 美国西门子医疗解决公司 | Shear wave velocity solution for medical ultrasound imaging |
| US20130218011A1 (en) * | 2012-02-16 | 2013-08-22 | Siemens Medical Solutions Usa, Inc. | Visualization of Associated Information in Ultrasound Shear Wave 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 |
-
2015
- 2015-08-20 CN CN201510514501.6A patent/CN105212961B/en not_active Expired - Fee Related
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100240994A1 (en) * | 2009-03-23 | 2010-09-23 | The Hong Kong Polytechnic University | Method and apparatus for ultrasound imaging and elasticity measurement |
| 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 |
| US20130218011A1 (en) * | 2012-02-16 | 2013-08-22 | Siemens Medical Solutions Usa, Inc. | Visualization of Associated Information in Ultrasound Shear Wave Imaging |
| 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》 * |
| NED C. ROUZE，MARK L. PALMERI: "Robust Estimation of Time-of-Flight Shear Wave Speed Using a Radon Sum Transformation", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》 * |
Cited By (10)
| 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 |
| CN106645406A (en) * | 2016-12-02 | 2017-05-10 | 北京空间飞行器总体设计部 | Positioning system and method of spacecraft subjected to space junk collision |
| 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 |
| CN107753058A (en) * | 2017-11-22 | 2018-03-06 | 深圳中科乐普医疗技术有限公司 | A kind of shearing wave dynamic filter method |
| CN109589138A (en) * | 2018-11-26 | 2019-04-09 | 深圳中科乐普医疗技术有限公司 | A kind of shear-wave velocity calculation method and elastogram equipment |
| CN109512465A (en) * | 2018-12-19 | 2019-03-26 | 深圳中科乐普医疗技术有限公司 | A kind of acoustic radiation force double direction shear wave composite imaging method and device |
| CN109512465B (en) * | 2018-12-19 | 2021-04-06 | 深圳中科乐普医疗技术有限公司 | Acoustic radiation force bidirectional shear wave composite imaging method and device |
| WO2021082043A1 (en) * | 2019-10-31 | 2021-05-06 | 汕头市超声仪器研究所股份有限公司 | First-order estimation method for shear wave motion velocity |
| US11903769B2 (en) | 2019-10-31 | 2024-02-20 | Shantou Institute Of Ultrasonic Instruments Co., Ltd. | First-order estimation method of shear wave velocity |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105212961B (en) | 2018-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105212961A (en) | A kind of acoustic radiation shear-wave velocity detection method and system | |
| CN104535655B (en) | A kind of ray tracing formula ultrasonic Lamb wave defect chromatography imaging method | |
| EP2317336B1 (en) | Method for estimating target range error and sonar system thereof | |
| CN1696734A (en) | Determination of time-difference of arrival and angle of arrival | |
| Nosal | Methods for tracking multiple marine mammals with wide-baseline passive acoustic arrays | |
| JP2010511172A5 (en) | ||
| US11399805B2 (en) | Ultrasound diagnostic device and ultrasound signal processing method | |
| CN104363836A (en) | System and method for shear wave elastography by transmitting ultrasound with subgroups of ultrasound transducer elements | |
| CN108169327B (en) | MUSIC corrosion monitoring process based on excitation beam forming and weighted image fusion | |
| CN107942335A (en) | A kind of object identification method and equipment | |
| US10194889B2 (en) | Methods, systems and computer program products for multi-resolution imaging and analysis | |
| CN103885057A (en) | Self-adaptation variable-sliding-window multi-target tracking method | |
| CN105488795A (en) | Composite material damage identification method | |
| CN104765064A (en) | Microseism interference imaging method | |
| CN112650300B (en) | Unmanned aerial vehicle obstacle avoidance method and device | |
| CN101839893A (en) | Lamb wave virtual time reversal method with high spatial resolution | |
| CN105676205A (en) | Airborne LiDAR waveform data Gaussian decomposition method | |
| CN105241459A (en) | Delay estimation method and device used for indoor underwater target positioning | |
| CN106572838A (en) | Elasticity measurement and detection method and system | |
| CN103487796B (en) | A kind of method utilizing underwater acoustic channel Statistically invariant feature to realize passive ranging | |
| CN107003403A (en) | The method and apparatus of acoustic imaging | |
| CN110507360B (en) | Shear wave imaging method and system | |
| CN102636773A (en) | Single-element range ambiguity resistant method based on channel multipath characteristic | |
| CN103824082A (en) | Pedestrian detection method and detection system thereof | |
| US11635538B2 (en) | Equivalent linear velocity for first arrival picking of seismic refraction |
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 |