 Research
 Open Access
 Published:
A stochastic filtering approach to recover strain images from quasistatic ultrasound elastography
BioMedical Engineering OnLine volume 13, Article number: 15 (2014)
Abstract
Background
Modelbased reconstruction algorithms have shown potentials over conventional strainbased methods in quasistatic elastographic image by using realistic finite element (FE) or biomechanical model constraints. However, it is still difficult to properly handle the discrepancies between the model constraint and ultrasound data, and the measurement noise.
Methods
In this paper, we explore the usage of Kalman filtering algorithm for the estimation of strain imaging in quasistatic ultrasound elastography. The proposed strategy formulates the displacement distribution through biomechanical models, and the ultrasoundderived measurements through observation equations. Through this filtering strategy, the discrepancies are quantitatively modelled as one Gaussian white noise, and the measurement noise of ultrasound data is modelled as another independent Gaussian white noise. The optimal estimation of kinematic functions, i.e. the full displacement and velocity field, are computed through this Kalman filter. Then the strain images can be easily calculated from the estimated displacement field.
Results
The accuracy and robustness of our proposed framework is first evaluated in synthetic data in controlled conditions, and the performance of this framework is then evaluated in the real data collected from elastography phantoms and patients with favourable results.
Conclusions
The potential of our algorithm is to provide the distribution of mechanically meaningful strain under a proper biomechanical model constraint. We address the modeldata discrepancy and measurement noise by introducing process noise and measurement noise in our framework, and then the mechanically meaningful strain is estimated through the Kalman filter in the minimum mean square error (MMSE) sense.
Introduction
A tumour or a suspicious pathological growth is normally much stiffer than the background of normal soft tissue [1]. So when a mechanical compression or vibration is applied, the tumour deforms less than the surrounding tissue, i.e. the strain in the tumour is less than the surrounding tissue. Hence a strain image may, under particular simplifying assumptions, be interpreted as representative of the underlying Young’s modulus distribution. Ultrasound elastography [1] is a noninvasive technique in which stiffness or strain images of soft tissue are used to detect or classify tumours or cancers. Similar elastography techniques have been also developed using magnetic resonance imaging [2–4] and mammographic imaging [5]. As an emerging medical imaging modality, it has been broadly applied in the clinical applications, including improving the diagnostic accuracy of breast and prostate cancer [6–9], assessing plaque vulnerability [10–13], guiding minimally invasive therapy [14–16]. Currently, supersonic shear imaging (SSI) has become one cuttingedging technique for realtime visualisation of soft tissue mechanical properties. Using ultrasonic focused beams, SSI is able to remotely generate mechanical vibration sources radiating lowfrequency, shear waves inside tissues, and then reconstruct the viscoelastic properties of tissue from the propagated shear waves [17]. Among different forms of ultrasound elastography, quasistatic ultrasound elastography is one popular technique because of its simplicity: two sets of radiofrequency (RF) echo frames are first collected before and after deformation caused by external compressions, and then the displacement field, strain field and even the distribution of tissular elasticity are reconstructed from RF signals using different approaches [18]. Numerous techniques have been proposed to solve the reconstruction problem of strain images for quasistatic ultrasound elastography in the last decades [19], but the estimation of strain images in ultrasound elastography still remains as a challenging researching topic at present because of its illposed nature [20].
The general procedure of quasistatic elastography is first to recover internal tissue displacements from two ultrasound radiofrequency (RF) frames before and after tissue deformation due to a quasistatic external compression, and then strain image and even mechanical properties distribution, such as Young’s modulus, can be computed from the displacement field under certain assumptions. So the quality of tissue motion estimation determines the performance of quasistatic ultrasound elastography. The estimation of tissue displacements is inherently a threedimensional problem, which means the displacement vector components physically involve in all three directions (x, y, z axis) simultaneously and continuously. However, early methods only focus on axial displacement estimation. Among these techniques, time delay estimation (TDE) is an important and widely used displacement estimation method [1]. Typically, TDE method is to find the bestmatching segment in the delayed RF signal for a specific segment in the reference RF signal by computing the maximum or minimum of a pattern matching function. Crosscorrelation is mostly used as a pattern matching function in TDE method, and several other matching techniques have been also employed, such as correlation coefficient, hybrid sign correlations, sum absolute differences (SAD) and sum squared differences (SSD). TDE provides accurate estimation of axial displacement, but it is normally timeconsuming. In the quasistatic elastography, the motion of tissue caused by the compression from the probe normally occurs in the corresponding twodimensional scanning plane. However, the spatial resolution in the axial direction is much higher than that in the lateral direction in ultrasound imaging, that is why the estimation of the axial component of the motion has received the most attention [1, 21]. Recently, quite a few approaches have been proposed to recover the lateral and shear strain fields (i.e. reconstruct the strain tensor), Poisson’s ratio and Young’s modulus [21], but their performance is still compromised by the low resolution of lateral displacements and the noise from RF data.
Both for the calculation of strain images and the reconstruction of elastic parameters, accurate estimation of tissue displacements is the first important step that will critically affect the image qualities. Therefore, different motion tracking techniques have been developed to recover tissue displacements during the past two decades [21]. Different beamforming schemes have been proposed to improve the lateral resolution. For instance, using large beam steering angle can obtain better lateral displacement measurements, but this method can be only implemented by the phased array. Different postprocessing techniques have also been applied to improve the quality of lateral displacement, such as iteratively interpolation along the lateral direction or local affine transformation. In another work, a twodimensional displacement field is calculated using the analytic minimization (AM) of cost functions that incorporate both the similarity of the amplitudes of RF signals and the displacement continuity [21]. However, lateral displacement estimation is still of order of magnitude less accurate than the axial displacement estimation. Thus, lateral strain failed in identifying the ablation lesion in patient experiments [21]. Researchers have also tried to integrate biomechanical constraint to recover highresolution lateral displacements from the noisy data. With the assumption of a constant Poisson’s ratio (0.49), i.e. based on the biomechanical constraint of tissue incompressibility [22], lateral displacements were recovered from axialstrain measurements using the leastsquare technique. However, the leastsquare technique in this tissueincompressibilityassumption (TIA) method cannot perfectly eliminate measurement noise. In addition, errors will be introduced when the tissue incompressibility assumption is invalid. The recovery of strain and elasticity of tissue from the sparse displacement is one inverse problem [2, 23, 24], therefore, it is necessary to develop a robust framework with a more meaningful biomechanical constraint to recover a full displacement field from the ultrasonic measurements.
The rest of this paper is organised as follows: section Methodology describes the details of our elasticity reconstruction method. The linear FE model is first introduced, followed by the discussion about the integration of the stochastic state space strategy and the dynamic equation for the quasistatic ultrasound elastography. In section Experiments and results, the results obtained from the experiments using simulated data, real phantom data and invivo clinical data are presented and discussed respectively. In each set of data, our results are compared to the strainbased maps. As the ground truths are available in simulation and phantom study, the results generated by our strategy are compared to the ideal results. In the clinical experiment, the elastographic images are compared to CT images with lesions indicated by the doctor. Finally, the discussion of our framework, concluding remarks and future research endeavours are outlined in section Discussion and section Conclusion respectively.
Methodology
Linear elasticity
In order to construct a realistic, yet computationally feasible, analysis framework using the imaging data and other available physical measurements such as pressure induced by the probe, the structure and material of the biological tissue should be properly modelled. For computational simplicity, in our current 2D implementation, we adopt the linear isotropic continuum material for the quasistatic ultrasound, where the stress and strain relationship obeys the Hooke’s law. Hence, in our current two dimensional implementation the linear isotropic continuum material is adopted to describe the mechanical behaviour of the tissue in quasistatic elastography for the computational simplicity, and the relationship of the stress σ and the strain ε obeys the Hooke’s law [25]:
where S is the strainstress matrix.
Let u(x,y) and v(x,y) be the displacement along the x and yaxis of a point, the infinitesimal strain tensor of the point is:
Under the plane strain condition [25], strainstress matrix S can be derived as:
S is a materialdependent matrix, E stands for the Young’s modulus, and ν stands for the Poisson’s ratio. Here, the Young’s modulus E and the Poisson’s ratio ν are two materialspecific parameters. This fact is quite clear from these relationships that the internal stress caused by the deformation is a function of the displacement vector and the material parameters. Ideally, the problem could be tackled in three dimension to avoid the throughplane motion effect, but in this paper we only study the twodimensional problem because the image of quasistatic ultrasound elastography are two dimensional data collected by the linear array of ultrasonic probe. Further, more realistic biomechanical models can also be adopted into current framework in same way, but a linear material model is used to illustrate the fundamental ideas and rationales of our filtering estimation strategy to recover the distribution of strains for quasistatic ultrasound elastography.
Stochastic finite element method
In the past decades, the deterministic finite element method has been able to provide an effective and convenient platform for biomechanical studies in the past decades [26, 27]. However, it does not have the capability to process the uncertainties of material properties, and kinematic observations, i.e. measurements of the displacement. Especially the ultrasound imagingderived data are usually corrupted by noises of various nature. It is thus necessary to develop a strategy, the stochastic finite element method (SFEM), which has been used widely for structural dynamic analysis in probability analysis frameworks [28, 29]. In SFEM, structural material properties are described by random fields, possibly with known prior statistics.
In our implementation, a Delaunay triangulated finite element mesh is constructed at the first image frame before compression. An isoparametric formulation defined in a natural coordinate system is used, where, for trinodal linear element, the basis functions are linear functions of the nodal coordinates [25]. With mesh constructed, assuming that the material parameters E (Young’s modulus) and ν (Poisson’s ratio) are temporally constant but varying spatially, we can have the following dynamic governing equation [25]:
where U is the displacement vector, and R is the load vector. Mass matrix M here is a known function of the material density and is temporally and spatially constant. Stiffness matrix K is a function of Young’s modulus and Poisson’s ratio determined by the material constitutive law [25]. Damping matrix C is frequency dependent, and we assumed the presence of small proportional Rayleigh damping and so C = α M + β K in our implementation, where α and β are the the mass and stiffnessproportional Rayleigh damping weighting coefficients, respectively [25]. In practice it is difficult to determine the damping parameters because they are frequency dependent. Our assumption of Raleigh damping was based on the very low damping exhibited by biological tissues during quasistatic elastography, and fixed the two weighting coefficients at 1%.
One initialisation issue of Equation (4) is how to measure the external loading vector (R) during freehand elastography. Considering the object system dynamics embodied in Equation (4), if any knowledge of the displacement vector (U) is available, it can be used as essential boundary conditions to recover the motion parameters of all other nodes. The following experiments involving synthetic and real imaging data provided a set of displacements at nodal points of the boundary (e.g., axial displacements), and they are employed in the following fashion. Let U _{ b } = b be known from the imaging data at selected sampling nodes of the boundary, then the additional constraining equation μ U _{ b } = μ b is enforced on the system dynamics through
where weighting coefficient μ depends on the confidence of each displacement, with large μ values (1 × 10^{4} in this study) indicating highly trustworthy data points and small μ values for others. In this way it remains possible to describe the boundary condition without measuring the external force during freehand elastography. More details of this enforcement of boundary condition can be found elsewhere [25].
State space representation
In order to employ our filtering strategy to integrate biomechanical model, the dynamics equation (Equation (4)) needs to be transformed into a statespace representation of the continuoustime linear stochastic system. First the kinematic vector x(k) and the material parameter vector θ are defined as:
where the kinematic vector x(t) is consisted of displacement U(t) and velocity $\stackrel{\u0307}{U}(t)$ information, and the material parameter vector θ is consisted of Young’s modulus E and Poisson’s ratio ν. In general, tumour inclusion or tissue blocked from its blood nutrients is stiffer than normal tissue, which mostly reflects in the variation of E. Since benign and cancerous tumours usually have distinguishing elastic properties, i.e. different Young’s modulus value, the value of Poisson’s ratio can be fixed in the following implementation. For example, as one of incompressible materials, the Poisson’s ratio of tissue can be set close to 0.5.
The state space representation of Equation (4) becomes
where n(t) is the process noise and it is white noise (E[ n(t)] = 0; E[ n(t)n(s)^{′}] = Q _{ n }(t)δ _{ ts }, where Q _{ n } is the process noise covariance). The system matrices A _{ c }, B _{ c } and input forces W(t) are:
An associated measurement equation, which describes the observations y(t) provided by the ultrasonic RF data, can be expressed in the form:
where e(t) is the measurement noise and it is white noise (E[e(t)] = 0; E[e(t)e(s)^{′}] = R _{ e }(t)δ _{ ts }, and R _{ e } is the measurement noise covariance), independent of n(t). H is the measurement matrix which should be specified by the relation between state vector x(t) and measurement vector y(t).
In order to run the dynamic equation (Equation (4)) in computer, it should be transformed into a discrete statespace equation, typically seen in control and estimation literature [30, 31]. We discretize Equation (8) and (9) over a constant time interval T. Since the interval T is always a known constant, we can replace kT with k in following equations
where
The associated measurement equation can be:
where y(k) is the measurement vector contained the displacement extracted from ultrasound RF data. In most quasistatic elastography, only theaxial component of displacement vector is extracted. Therefore, the corresponding place of D will be set to 1 or 0 according to the available measurement data. n(k) and e(k) are process noise and the measurement noise respectively.
Experiments and results
Our filtering framework is validated in simulated data, real phantom data and clinical data respectively. All the programs in the following experiments are coded in Matlab 2009a (MathWorks Corporation, USA) and run in a desktop computer with a core 2 intel CPU and 12 G memory.
Simulated data
ABAQUS (DS Simulia Corporation, USA) software is used for generating simulation data, and our phantom model is shown in Figure 1. Because the displacements are collected by the ultrasound probe along axial direction and the lateral displacement is much less accurate than axial displacement, we only use the axial displacements as input data and estimate the full displacement field through our method. In this model, a 20×20 mm rectangular object (E _{ b } = 25 kPa and ν _{ b } = 0.49) is built as the background. A 4 mm hard circular inclusion with different elastic property (E _{ t } = 80 kPa and ν _{ t } = 0.49) is included in the centre of the model, as the target. The top side of the object is compressed to move 1.0 mm down along the vertical direction with a deformation ratio 5% and the bottom side is fixed. After the rectangular with circular is divided into a triangular mesh and the mechanical deformation is calculated by ABAQUS, the Gaussian noise is added into the axial displacements for generating simulated the measurements. The lateral displacement is recovered by our proposed Kalman filter method, then a linear interpolation function is applied to smooth the displacement image. The least square method is used to calculate the strain field. We design two experiments on the simulated data in different conditions:

1.
Since we have the ground truth of simulated measurements, the first experiment is designed to examine the influence of measurement noise. In order to evaluate the performance quantitatively, we define the following two factors, the lateral displacement relative error (DRE) and contrasttonoise ratio (CNR):
$$\begin{array}{l}\mathit{\text{DRE}}=\frac{\overline{e}}{\overline{d}}\end{array}$$(14)$$\begin{array}{l}\mathit{\text{CNR}}=\sqrt{\frac{2{(\overline{{s}_{b}}\overline{{s}_{t}})}^{2}}{{\sigma}_{b}^{2}+{\sigma}_{t}^{2}}}\end{array}$$(15)
where $\overline{e}$ is the average nodal displacement error, $\overline{d}$ is the average nodal displacement value, $\overline{{s}_{t}}$ and $\overline{{s}_{b}}$ are the spatial strain average of the target and background, ${\sigma}_{t}^{2}$ and ${\sigma}_{b}^{2}$ are the spatial strain variance of the target and background. Figure 2 and Figure 3 show the lateral displacement and strains estimated results when the axial displacement is added by the Gaussian noise (SNR = 40 dB). Figure 4 and Figure 5 show the estimated results in different noise conditions (40 db, 20 db, 15 db). Table 1 shows the comparison of lateral displacement errors and CNR under different measurement noise levels.

2.
The second experiment is designed to evaluate the effect of initialisation of the model elastic property. In our method Young’s modulus distribution is assumed to be homogeneous when the whole displacement field is estimated. Figure 6 and Table 2 show the effect of Young’s modulus in lateral displacement estimation. The performance of our method is acceptable when Young’s modulus is greater than 20 kPa, but the CNR is in a higher level only from 20 kPa to 40 kPa.
Real phantom data
Two sets of phantom data are used to evaluate the performance of our method. The first set of phantom data is collected by us in one Elasticity QA Phantom (Model 049, CIRS Inc. USA) using a PCbased Ultrasonix RP ultrasound machine (Ultrasonix Medical Corporation, Burnaby, BC, Canada). One phaseshift method with prior estimates [32] is used to estimate the axial displacements from the RF signals, and the full displacement field is recovered from the axial displacement measurements. And then all the strain images are displayed in Figure 7. We also compare our method to the method developed by Johns Hopkins University (JHU) [21]. In this comparison, the second set of phantom data is provided by JHU and it is available online [21]. Figure 8 shows the recovered strain images from both methods and the strain profiles in the same depth.
Clinical data
The ultrasound experiments were performed with the approval of the Health Science Research Ethics Committee of JHU. The participants were provided written informed consent before beginning the experiment. Three sets of patients’ ultrasonic RF data are provided online in [21], and those data are collected in the patients undergoing open surgical RF thermal ablation for liver cancer enrolled between February 06, 2008 and July 28, 2009. The details of these patients’ clinical status have been extensively discussed in [21]. All the patients’ data are processed in the same way as the way for real phantom data. Since the axial and lateral displacement measurements have been provided in [21], the extraction of displacement is not necessary in this experiment. Figure 9 shows patients’ CT images scanned after RF ablation, strain images and our estimated results respectively.
Discussion
Our modelbased stochastic filtering strategy is able to provide a a highresolution displacement field (in both axial and lateral directions) in 2D elastography. In the following context the robustness, initialisation and computational cost of our filtering framework are discussed respectively.
Robustness of our filtering strategy
The kernel of our proposed approach is the Kalman filtering algorithm, therefore the tolerance of external noise is the one key performance factor of our filtering strategy. Our filtering framework can recursively absorb sparse measurements under the biomechanicalmodelbased constraint, and finally calculate physically meaningful optimal estimates, which is a full and smoother displacement field. In order to examine the robustness of our filtering strategy, we have designed several experiments to test the tolerance of the noise in our method.The tolerance of the proposed filtering strategy to external noise was evaluated by adding different levels of measurement noise to the simulated measurement. Our filtering framework was still able to recover the displacement/strain fields with a similar quality from measurements containing different levels of measurement noise after certain number of iterations, as indicated by the three pictures in Figure 4. But if we examine the tolerance of noise in much finer scale, we can see the quality of estimated displacement/strain field will not change when SNR is greater than 40 db as shown in Figure 5. This fact can be also verified by quantitative measures in Table 1.
Initialisation issues
The number of iterations and the initialisation of biomechanicalmodel constraint (Young’s modulus and Poisson’s ratio) are important to our filtering strategy. In order to generate a displacement/strain image of high quality, we simply performed 200 iterations of the filtering procedures as show in Figure 4(c). Since most biological tissues are incompressible materials, Poisson’s ratio used in our filtering strategy can be initialised to 0.49. When Young’s modulus is initialised greater than 20 kPa, the average position error will not increase anymore (see Figure 6(a) and Table 2), but a better CNR can only be obtained when Young’s modulus is from 20 kPa to 40 kPa (see Figure 6(b)). Therefore, the initial value of Young’s modulus should be in the range from 20 kPa to 40 kPa.
Computational cost
The computational cost of our method is the high in current implement scheme. In all the experiments performed in this study we can only handle measurements on 21 × 21 sample nodes. Although the spatial resolution of the recovered strain image is actually limited by the number of nodes, the quality of the estimated strain images generated by our modelbased filtering strategy is not yet compromised. Furthermore, we can take advantage from advanced computing power technologies, such as GPU, to accelerate our filtering strategy in future study.
Conclusion
In this paper, we have developed a biomechanically constrained filtering framework to extract a full displacement field from the measurements derived from ultrasound RF signals collected before and after deformation. Then, the strain tensor can be reconstructed from the full displacement field. The proposed framework have been validated by a series of experiments. However, the spatial resolution of recovered displacement field and strain images are restricted by the number of finite element nodes. In the future works, more efforts should be done on computational complexity reduction in order to increase spatial resolution.
Abbreviations
 MMSE:

Minimum mean square error
 RF:

Radio frequency
 SFEM:

Stochastic finite element method
 FEM:

Finite element method
 FE:

Finite element
 TIA:

Tissueincompressibilityassumption
 SNR:

Signaltonoise ratio
 DRE:

Displacement relative error
 CNR:

Contrasttonoise ratio
 JHU:

Johns Hopkins University.
References
 1.
Ophir J, Cespedes I, Ponnekanti H, Yazdi Y, Li X: Elastography: a quantitative method for imaging the elasticity of biological tissues. Ultrason Imaging 1991,13(2):111–134.
 2.
Sumi C, Nakayama K: A robust numerical solution to reconstruct a globally relative shear modulus distribution from strain measurements. IEEE Trans Med Imaging 1998, 17: 419–428. 10.1109/42.712131
 3.
Plewes D, Bishop J, Samani A, Sciarretta J: Visualization and quantification of breast cancer biomechanical properties with magnetic resonance elastography. Phys Med Biol 2000, 45: 1591–1610. 10.1088/00319155/45/6/314
 4.
McGee KP, Lake D, Mariappan Y, Hubmayr RD, Manduca A, Ansell K, Ehman RL: Calculation of shear stiffness in noise dominated magnetic resonance elastography data based on principal frequency estimation. Phys Med Biol 2011,56(14):4291–4309. 10.1088/00319155/56/14/006
 5.
Miga M: A new approach to elastography using mutual information and finite elements. Phys Med Biol 2003, 48: 467–480. 10.1088/00319155/48/4/304
 6.
Garra B, Cespedes E, Ophir J: Elastography of breast leisons: initial clinical results. Radiology 1997, 202: 79–86.
 7.
Souchon R, Rouviere O, Gelet A: Visualisation of HIFU lesions using elastography of the human prostate in vivo: preliminary results. Ultrasound Med Biol 2003,29(7):1007–1015. 10.1016/S03015629(03)000656
 8.
Hiltawsky KM, Kruger M, Starke C: Freehand ultrasound elastography of breast lesions: clinical results. Ultrasound Med Biol 2001,27(11):1461–1469. 10.1016/S03015629(01)004343
 9.
Samani A, Bishop J, Plewes D: A constrained modulus reconstruction technique for breast cancer assessment. IEEE Trans Med Imaging 2001,20(9):877–885. 10.1109/42.952726
 10.
Brusseau E, Fromageau J, Finet G, Delachartre P, Vray D: Axial strain imaging of intravascular data: results on polyvinyl alcohol cryogel phantoms and carotid artery. Ultrasound Med Biol 2001,27(12):1631–1642. 10.1016/S03015629(01)004513
 11.
Korte C, Pasterkamp G, van der Steen A, Woutman H, Bom N: Characterization of plaque components with intravascular ultrasound elastography in human femoral and coronary arteries in vitro. Circulation 2000,102(6):617–623. 10.1161/01.CIR.102.6.617
 12.
de Korte C, Carlier S, Mastik F: Morphological and mechanical information of coronary arteries obtained with intravascular elastography. Feasibility study in vivo. Eur Heart J 2002,23(5):405–413. 10.1053/euhj.2001.2806
 13.
Doyley M, Mastik F, de Korte C: Advancing intravascular ultrasonic palpation toward clinical applications. Ultrasound Med Biol 2001b,27(11):1471–1480. 10.1016/S03015629(01)004574
 14.
Kallel F, Stafford R, Pric R: The feasibility of elastographic visualization of HIFUinduced thermal lesions in soft tissues. Imageguided highintensity focused ultrasound. Ultrasound Med Biol 1999,25(4):641–647. 10.1016/S03015629(98)001847
 15.
Righetti R, Kallel F, Stafford R: Elastographic characterization of HIFUinduced lesions in canine livers. Ultrasound Med Biol 1999,25(7):1099–1113. 10.1016/S03015629(99)000447
 16.
Varghese T, Techavipoo U, Liu W: Elastographic measurement of the area and volume of thermal lesions resulting from radiofrequency ablation: pathologic correlation. Am J Roentgenol 2003,181(3):701–707. 10.2214/ajr.181.3.1810701
 17.
Bercoff J, Tanter M, Fink M: Supersonic shear imaging: a new technique for soft tissue elasticity mapping. IEEE Trans Ultrason Ferroelectr Freq Control 2004,51(4):396–409.
 18.
Shiina T, Nitta N, Ueno E, Bamber J: Real time tissue elasticity imaging using the combined autocorrelation mehtod. J Med Ultrason 2002,29(Auturmn):119–128.
 19.
Ophir J, Alam S, Garra B, Kallel F, Konofagou E, Krouskop T, Merritt C, Righetti R, Souchon R, Srinivasan S, Varghese T: Elastography: imaging the elastic properties of soft tissues with ultrasound. J Med Ultrasound 2002, 29: 155–171. 10.1007/BF02480847
 20.
Doyley M, Srinivasan S, Pendergrass S, Wu Z: Comparative elvaluation of strainbased and modelbased modulus elastography. Ultrasound Med Biol 2005,31(6):787–802. 10.1016/j.ultrasmedbio.2005.02.005
 21.
Rivaz H, Boctor E, Choti M, Hager G: Realtime regularized ultrasound elastography. IEEE Trans Med Imaging 2010. in press
 22.
Lubinski M, Emelianov S, Raghavan K, Yagle A, Skovoroda A, O’Donnell M: Lateral displacement estimation using tissue incompressibility. IEEE Trans Ultrason Ferroelectr Freq Control 1996,43(2):247–256.
 23.
Samani A, Plewes D: An inverse problem solution for measuring the elastic modulus of intact ex vivo breast tissue tumours. Phys Med Biol 2007,52(5):1247–1260. 10.1088/00319155/52/5/003
 24.
Luo J, Ying K, Bai J: Elasticity reconstruction for ultrasound elastography using a radial compression: an inverse approach. Ultrasonics 2006, 44: e195e198.
 25.
Bathe K: Finite Element Procedures in Engineering Analysis.. Prentice Hall Inc.; 1982.
 26.
Benayoun S, Ayache N: Dense nonrigid motion estimation in sequences of medical images using differential constraints. Int J Comput Vis 1998, 26: 25–40. 10.1023/A:1007932523504
 27.
Creswell LL, Moulton MJ, Wyers SG, Pirolo JS, Fishman DS, Perman WH, Myers KW, Actis RL, Vannier MW, Szabo BA, Pasque MK: An experimental method for evaluating constitutive models of myocardium in vivo hearts. Am J Physiol 1994, 267: H853H863.
 28.
Contreras H: The stochastic finite element method. Comput Struct 1980, 12: 314–348.
 29.
Liu W, Belytschko T, Man A: Random field finite element. Int J Numerical Methods Eng 1980, 23: 1831–1845.
 30.
Glad T, Ljung L: Control Theory. London: Taylor and Francis; 2000.
 31.
Shi P, Liu H: Stochastic finite element framework for simultaneous estimation of cardiac kinematic functions and material parameters. Med Image Anal 2003, 7: 445–464. 10.1016/S13618415(03)000665
 32.
Yuan JW, Zhang HY, Lu MH, Chen SP, Liu HF: A strainbased ultrasound elastography using phase shift with prior estimates and meshfree shape function. IEEE ISBI. 2011, 532–535.
Acknowledgements
This work was supported in part by National Basic Research Program 973 (2010CB732606), 863 Program of China (No.2012AA02A603), Guangdong Innovation Research Team Fund for LowCost Healthcare Technologies in China, External Cooperation Program of the Chinese Academy of Sciences (GJHZ1212), Key Lab for Health Informatics of Chinese Academy of Sciences, Distinguished Youth Talents Fund of Shenzhen (JC201006020025A), the National Natural Science Foundation of China (81101120, 61233012, 61031003, 61170076), and Guangdong Province Science and technology plan project (No. 2011B031800104).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
Dr. ML and Miss DW contributed equally to this work and the manuscript writing as they both conceived of the study, designed the whole work, and drafted the manuscript. Mr. WL and Miss WhL carried out the experiments on phantom and synthetic data. Dr. HZ participated in the design of this study and coordination and helped to draft the manuscript. Dr. WH provided the clinical advices into this manuscript. All authors read and approved the final manuscript.
Minhua Lu, Dan Wu contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Lu, M., Wu, D., Lin, W. et al. A stochastic filtering approach to recover strain images from quasistatic ultrasound elastography. BioMed Eng OnLine 13, 15 (2014). https://doi.org/10.1186/1475925X1315
Received:
Accepted:
Published:
Keywords
 Lateral Displacement
 Axial Displacement
 Strain Image
 Time Delay Estimation
 Ultrasound Elastography