CN102920481A - Position estimation method of ultrasonic elastography one-dimensional axial displacement estimation window - Google Patents
Position estimation method of ultrasonic elastography one-dimensional axial displacement estimation window Download PDFInfo
- Publication number
- CN102920481A CN102920481A CN2012104861180A CN201210486118A CN102920481A CN 102920481 A CN102920481 A CN 102920481A CN 2012104861180 A CN2012104861180 A CN 2012104861180A CN 201210486118 A CN201210486118 A CN 201210486118A CN 102920481 A CN102920481 A CN 102920481A
- Authority
- CN
- China
- Prior art keywords
- window
- estimation
- displacement
- coordinate
- displacement estimation
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention discloses a position estimation method of an ultrasonic elastography one-dimensional axial displacement estimation window, which comprises the following steps of: 1, setting an initial position of the center of the one-dimensional axial displacement estimation window as a displacement estimation sampling point; 2, generating a position coordinate of the one-dimensional axial displacement estimation window by using an iteration method; 3, using a coordinate generated by iteration as a coordinate of a center point of the displacement estimation window, and taking the one-dimensional window with the point as the center for carrying out axial displacement estimation; and 4, carrying out gradient operation on a displacement estimation value to generate an elastic image. According to the invention, by adopting a new position estimation algorithm of the one-dimensional estimation window, the displacement estimation is carried out on the two-dimensional estimation window so that a more accurate displacement of a preset sampling point can be generated, and the influence of signal amplitude random fluctuation and signal compression to displacement estimation is avoided, therefore, amplitude modulation noise of elastography is effectively inhibited, strain texture fluctuation is avoided being caused, and the elastography quality and quality are improved.
Description
Technical field
The invention belongs to the medical ultrasonic engineering technical field, the method for estimation of one-dimensional estimated window position when being specifically related to carry out the axial displacement estimation in a kind of ultrasonic elastograph imaging.
Background technology
Ultrasonic elastograph imaging is a kind of novel imaging pattern that produces the tissue elasticity characteristic information and show with image mode by the correlation analysis to ultrasonic echo before and after the compression.Owing to can reflect the variation of tissue elasticity parameter, this imaging pattern has become clinically various diseases detection and the important means of diagnosing, such as breast carcinoma, carcinoma of prostate, thyroid nodule and hepatic fibrosis etc.
The ultrasonic elastograph imaging process generally includes Displacement Estimation, strain is estimated and show post processing three large links.Displacement Estimation is core link in the elastogram, and its estimated accuracy directly affects the elastogram quality.Although image-based motion tracking displacement estimation method is feasible, but because radiofrequency signal comprises abundanter original structure information and phase information (and image usually passed through go phase place and down-sampling to process), Displacement Estimation is done the more accurate estimated value of generations in the signal aspect.
Current, based on the displacement estimation method of signal be: preset displacement is estimated sampled point in the compression front signal, and then fetch bit moves estimating window centered by these points; Carry out motion tracking with this window, seek the signal window that matches in the echo-signal frame after compression; Calculate the time shift between two signal windows, be the estimation displacement of the time delay form of preset.As shown in Figure 1, " * " expression Displacement Estimation sampled point, " " expression estimating window central point; Fig. 1 adopts traditional algorithm, and getting the estimating window central point is the Displacement Estimation sampled point, and namely window center point overlaps with the Displacement Estimation sampled point.Yet this simply method that fetch bit moves estimating window centered by the Displacement Estimation sampled point has the following disadvantages and defective:
(1) owing to affected by the signal amplitude random fluctuation, the displacement that current displacement algorithm for estimating estimates not is the precise displacement of estimating window central point, namely is not the real displacement of preset.The strain pattern that still uses the gradient operation to produce with this displacement will produce obvious noise;
(2) the current displaced position algorithm for estimating that has proposed at first estimates position corresponding to displacement, then carries out strain with displacement difference divided by alternate position spike and calculates.The method can improve the quality of imaging to a certain extent, but the strain that this quadrat method produces not is the strain of preset, therefore can produce the image texture fluctuation problem.
Therefore, expect the precise displacement of preset rules point, reduce the picture amplitude zoop, prevent the image texture fluctuation, improve image quality, need to estimate the position of displacement estimating window.
Summary of the invention
For above shortcomings in the prior art, the invention provides a kind of location estimation method in order to the ultrasonic elastograph imaging one dimension axial displacement estimating window that improves elastogram performance and quality.
In order to solve the problems of the technologies described above, the present invention has adopted following technical scheme:
The location estimation method of ultrasonic elastograph imaging one dimension axial displacement estimating window, the method comprises the steps:
(1), use the ultrasound probe compress tissue, backscattering echo signal frame before and after the pressure contracting carries out the location estimation of estimating window to compression front signal frame;
(2), require to determine the long L of window and the displacement sampled point of one-dimensional estimated window according to signal to noise ratio and spatial resolution;
(3), the axial coordinate of Displacement Estimation sampled point is made as the initial coordinate of one-dimensional estimated window center point, that is:
, y here
0The axial coordinate of displacement sampled point, y
1It is the initial coordinate of one-dimensional estimated window center point;
(4), take signal amplitude as weight term, use iterative method to produce the axial location coordinate of one-dimensional estimated window center point, that is:
In the following formula, y
0The axial coordinate of Displacement Estimation sampled point, y
kBe the coordinate of the Displacement Estimation window of the k time iteration generation, k is iterations, k=1, and 2,3 ..., N(N is that iteration finishes number of times, is determined by the algorithm expectation quality), L is that window is long, H (y) is the weight term based on signal amplitude,
, wherein f (y) is signal amplitude;
(5), the coordinate y that produces with iteration
NAs the estimating window center point coordinate, centered by this point, get long window as L and carry out Displacement Estimation;
(6), carry out the gradient operation with this Displacement Estimation value and produce strain pattern.
The invention has the beneficial effects as follows: the present invention adopts the location-estimation algorithm of new one dimension displacement estimating window, carry out the more precise displacement that Displacement Estimation can produce default sampled point with this window, avoid signal amplitude random fluctuation and Signal Compression on the impact of Displacement Estimation, thereby the amplitude modulation(PAM) noise of establishment elastogram, also can avoid producing strain texture wave phenomenon, improve elastogram quality and quality.
Description of drawings
Fig. 1 is that traditional one-dimensional estimated window position is estimated sketch map;
Fig. 2 is that one-dimensional estimated window position of the present invention is estimated sketch map.
The specific embodiment
The invention will be further described below in conjunction with the drawings and specific embodiments.
The location estimation method of ultrasonic elastograph imaging one dimension axial displacement estimating window, the method comprises the steps:
(1), use the ultrasound probe compress tissue, backscattering echo signal frame before and after the pressure contracting carries out the location estimation of estimating window to compression front signal frame;
(2), require to determine the long L of window and the displacement sampled point of one-dimensional estimated window according to signal to noise ratio and spatial resolution;
(3), the axial coordinate of Displacement Estimation sampled point is made as the initial coordinate of one-dimensional estimated window center point, that is:
, y here
0The axial coordinate of displacement sampled point, y
1It is the initial coordinate of one-dimensional estimated window center point;
(4), take signal amplitude as weight term, use iterative method to produce the axial location coordinate of one-dimensional estimated window center point, that is:
In the following formula, y
0The axial coordinate of Displacement Estimation sampled point, y
kBe the coordinate of the Displacement Estimation window of the k time iteration generation, k is iterations, k=1, and 2,3 ..., N(N is that iteration finishes number of times, is determined by the algorithm expectation quality), L is that window is long, H (y) is the weight term based on signal amplitude,
, wherein f (y) is signal amplitude;
(5), the coordinate y that produces with iteration
NAs the estimating window center point coordinate, centered by this point, get long window as L and carry out Displacement Estimation;
(6), carry out the gradient operation with this Displacement Estimation value and produce strain pattern.
The present invention has adopted the location-estimation algorithm of new one dimension displacement estimating window, as shown in Figure 2, and " * " expression Displacement Estimation sampled point, " " expression estimating window central point.The estimating window central point that estimates with new method might not overlap with the Displacement Estimation sampled point.Owing to be subjected to the impact of signal amplitude random fluctuation, carry out the shift value of Displacement Estimation generation closer to the real displacement of Displacement Estimation sampled point with this window.The elastic image that produces with the method not only can effectively be avoided the impact of signal zoop, also can avoid the fluctuation of strain texture, thereby improve elastogram quality and quality.
Explanation is at last, above embodiment is only unrestricted in order to technical scheme of the present invention to be described, although with reference to preferred embodiment the present invention is had been described in detail, those of ordinary skill in the art is to be understood that, can make amendment or be equal to replacement technical scheme of the present invention, and not breaking away from aim and the scope of technical solution of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.
Claims (1)
1. the location estimation method of ultrasonic elastograph imaging one dimension axial displacement estimating window is characterized in that the method comprises the steps:
(1), use the ultrasound probe compress tissue, backscattering echo signal frame before and after the pressure contracting carries out the location estimation of Displacement Estimation window to compression front signal frame;
(2), require to determine the long L of window and the Displacement Estimation sampled point of one-dimensional estimated window according to signal to noise ratio and spatial resolution;
(3), the axial coordinate of Displacement Estimation sampled point is made as the initial coordinate of one-dimensional estimated window center point, that is:
, y here
0The axial coordinate of displacement sampled point, y
1It is the initial coordinate of one-dimensional estimated window center point;
(4), take signal amplitude as weight term, use iterative method to produce the axial location coordinate of one-dimensional estimated window center point, that is:
In the following formula, y
0The axial coordinate of Displacement Estimation sampled point, y
kBe the coordinate of the Displacement Estimation window of the k time iteration generation, k is iterations, k=1, and 2,3 ..., N, N are that iteration finishes number of times, determine that by the algorithm expectation quality L is that window is long, H (y) is the weight term based on signal amplitude,
, wherein f (y) is signal amplitude;
(5), the coordinate y that produces with iteration
NAs the estimating window center point coordinate, centered by this point, get long window as L and carry out Displacement Estimation;
(6), carry out the gradient operation with this Displacement Estimation value and produce strain pattern.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210486118.0A CN102920481B (en) | 2012-11-26 | 2012-11-26 | Position estimation method of ultrasonic elastography one-dimensional axial displacement estimation window |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210486118.0A CN102920481B (en) | 2012-11-26 | 2012-11-26 | Position estimation method of ultrasonic elastography one-dimensional axial displacement estimation window |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102920481A true CN102920481A (en) | 2013-02-13 |
| CN102920481B CN102920481B (en) | 2014-05-07 |
Family
ID=47635519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210486118.0A Expired - Fee Related CN102920481B (en) | 2012-11-26 | 2012-11-26 | Position estimation method of ultrasonic elastography one-dimensional axial displacement estimation window |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102920481B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070154069A1 (en) * | 2005-12-30 | 2007-07-05 | Pixart Imaging Inc., R.O.C. | Displacement estimation device and method for the same |
| US20080019609A1 (en) * | 2006-07-20 | 2008-01-24 | James Hamilton | Method of tracking speckle displacement between two images |
| CN101569543A (en) * | 2008-04-29 | 2009-11-04 | 香港理工大学 | Two-dimension displacement estimation method of elasticity imaging |
| CN102626327A (en) * | 2012-04-26 | 2012-08-08 | 声泰特（成都）科技有限公司 | Ultrasonic elastography and pressure feedback method based on receive-side spatial-compounding |
-
2012
- 2012-11-26 CN CN201210486118.0A patent/CN102920481B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070154069A1 (en) * | 2005-12-30 | 2007-07-05 | Pixart Imaging Inc., R.O.C. | Displacement estimation device and method for the same |
| US20080019609A1 (en) * | 2006-07-20 | 2008-01-24 | James Hamilton | Method of tracking speckle displacement between two images |
| CN101569543A (en) * | 2008-04-29 | 2009-11-04 | 香港理工大学 | Two-dimension displacement estimation method of elasticity imaging |
| CN102626327A (en) * | 2012-04-26 | 2012-08-08 | 声泰特（成都）科技有限公司 | Ultrasonic elastography and pressure feedback method based on receive-side spatial-compounding |
Non-Patent Citations (4)
| Title |
|---|
| JOEL E. LINDOP ET AL: "Phase-based ultrasonic deformation estimation", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》 * |
| LINDOP J.E ET AL: "Estimation of displacement location for enhanced strain imaging", 《ULTRASONICS, FERROELECTRICS AND FREQUENCY CONTROL, IEEE TRANSACTIONS ON》 * |
| 崔少国等: "基于相位的2D超声弹性成像形变估计", 《北京生物医学工程》 * |
| 崔少国等: "弹性成像中不同频率子图像的相关性与频率复合", 《中国生物医学工程学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102920481B (en) | 2014-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zahiri-Azar et al. | Motion estimation in ultrasound images using time domain cross correlation with prior estimates | |
| Rivaz et al. | Ultrasound elastography: a dynamic programming approach | |
| CN103040488B (en) | System and method for real-time ultrasonic elastography displacement estimation | |
| RU2689174C2 (en) | Autocorrelation-induced cross correlation in shear wave ultrasonic elastography | |
| CN104688266A (en) | Motion correction in three-dimensional elasticity ultrasound imaging | |
| CN102217953B (en) | Image tracking method and device based on multi-neighborhood-aided two-dimensional ultrasonic deformed microstructure | |
| Yuan et al. | Analytical phase-tracking-based strain estimation for ultrasound elasticity | |
| CN104739442B (en) | Compressive resilience imaging displacement detection method, device and supersonic imaging apparatus | |
| CN102860842A (en) | Real-time quasi-static ultrasonic elastography method | |
| Jalilian et al. | Increasing frame rate of echocardiography based on a novel 2d spatio-temporal meshless interpolation | |
| KR20130051241A (en) | Method for generating diagnosis image, apparatus and medical imaging system for performing the same | |
| CN110353729A (en) | A kind of sound wave transition time detection method based on two-way shot and long term memory network | |
| CN102920481B (en) | Position estimation method of ultrasonic elastography one-dimensional axial displacement estimation window | |
| KR101893383B1 (en) | Apparatus and method for generating ultrasonic image | |
| CN110554393B (en) | High-contrast minimum variance imaging method based on deep learning | |
| CN102920479B (en) | Position estimation method of ultrasonic elastography two-dimensional axial displacement estimation window | |
| Yuan et al. | A strain-based ultrasound elastography using phase shift with prior estimates and meshfree shape function | |
| CN112731520A (en) | Full waveform inversion method and system based on structure tensor diffusion filtering | |
| CN102920480B (en) | Ultrasonic elastography property enhancement method | |
| CN106108942A (en) | Parallel ARFI formation method based on OpenCL | |
| Ekroll et al. | Quantitative vascular blood flow imaging: A comparison of vector velocity estimation schemes | |
| Zhang et al. | Fast axial and lateral displacement estimation in myocardial elastography based on RF signals with predictions | |
| Revell et al. | Ultrasound speckle tracking for strain estimation | |
| KR101085220B1 (en) | Apparatus and method for embodying elastic image | |
| Revell et al. | Variable sized block matching for in vivo musculoskeletal motion analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140507 Termination date: 20161126 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |