- Research
- Open Access

# A regularization-free elasticity reconstruction method for ultrasound elastography with freehand scan

- Xiaochang Pan
^{1}, - Ke Liu
^{2}, - Jing Bai
^{1}and - Jianwen Luo
^{1, 3}Email author

**13**:132

https://doi.org/10.1186/1475-925X-13-132

© Pan et al.; licensee BioMed Central Ltd. 2014

**Received:**8 August 2014**Accepted:**2 September 2014**Published:**7 September 2014

## Abstract

### Background

In ultrasound elastography, reconstruction of tissue elasticity (e.g., Young’s modulus) requires regularization and known information of forces and/or displacements on tissue boundaries. In practice, it is challenging to choose an appropriate regularization parameter; and the boundary conditions are difficult to obtain *in vivo*. The purpose of this study is to develop a more applicable algorithm that does not need any regularization or boundary force/displacement information.

### Methods

The proposed method adopts the bicubic B-spline as the tissue motion model to estimate the displacement fields. Then the estimated displacements are input to the finite element inversion scheme to reconstruct the Young’s modulus of each element. In the inversion, a modulus boundary condition is used instead of force/displacement boundary conditions. Simulation and experiments on tissue-mimicking phantoms are carried out to test the proposed method.

### Results

The simulation results demonstrate that Young’s modulus reconstruction of the proposed method has a relative error of −3.43 ± 0.43% and root-squared-mean error of 16.94 ± 0.25%. The phantom experimental results show that the target hardening artifacts in the strain images are significantly reduced in the Young’s modulus images. In both simulation and phantom studies, the size and position of inclusions can be accurately depicted in the modulus images.

### Conclusions

The proposed method can reconstruct tissue Young’s modulus distribution with a high accuracy. It can reduce the artifacts shown in the strain image and correctly delineate the locations and sizes of inclusions. Unlike most modulus reconstruction methods, it does not need any regularization during the inversion procedure. Furthermore, it does not need to measure the boundary conditions of displacement or force. Thus this method can be used with a freehand scan, which facilitates its usage in the clinic.

## Keywords

- Optical Flow
- Axial Displacement
- Displacement Estimation
- Ultrasound Elastography
- Theoretical Displacement

## Background

Breast cancer is the most common cancer among women, and the second-leading cause of cancer deaths in women in the United States [1]. The pathological state of the breast cancer highly correlates with their mechanical properties, such as Young’s modulus (or shear modulus) and viscoelasticity [2]. This lays the foundation of manual palpation routinely used in breast cancer detection. Palpation is especially helpful in the detection and localization of breast lesions [3]. However, it is limited in the cases when the tumor is small or deep beneath the skin surface [4]. Besides, palpation is subjective and only provides qualitative information. Based on the concept of palpation, quasi-static elastography (or compression elastography) is proposed to non-invasively estimate the mechanical property of soft tissues using ultrasound imaging [5].

In quasi-static elastography, the axial strain field is interpreted as relative Young’s modulus distribution within the tissue with the assumption of constant stress distribution [5]. The true Young’s modulus distribution could be computed from the internal strain and stress field [6]. However, the internal stress distribution cannot be measured *in vivo*
[7]. Due to the non-uniform stress distribution within the tissue, several mechanical artifacts could exist in the axial strain image and may compromise the diagnosis in the clinic [8]. For examples, stress decay with depth causes the target-hardening artifact [5]. To overcome this limitation, many researchers are devoted to reconstructing the Young’s or shear modulus within the tissue by using certain constraints and the estimated displacements or strains [8–15].

The reconstruction of Young’s modulus is an inverse problem [8, 16]. Most inversion approaches assume that the tissue is linearly elastic, isotropic, continuous and incompressible [8, 11, 13–15]. In addition, the three-dimensional (3D) elasticity problem is usually simplified to a two-dimensional (2D) problem by assuming plane-strain [8] or plane-stress [14] conditions. The approaches for this inverse problem can be mainly divided into two categories, i.e., direct inversion approaches [14, 15, 17, 18] and iterative inversion approaches [8, 10, 16, 19–21].

The direct approaches compute the Young’s or shear modulus by solving the partial differential equations of equilibrium [7, 16]. These approaches rearrange the equilibrium equation used in the forward problem, and the Young’s modulus could be directly reconstructed by using special boundary conditions and estimated strain fields [16]. However, the boundary conditions are the known modulus values or pressures on the boundaries of the region of interest (ROI), which are difficult to measure *in vivo* [7]. In addition, the errors in the estimated strain field can be greatly amplified in the direct inversion, and thus may make the reconstruction result unreasonable [16].

The iterative approaches treat the inverse problem as a nonlinear optimization problem. The shear modulus distribution is calculated from minimization of the difference between the measured displacements and predicted displacements computed in the forward problem (i.e., the objective function) [7, 22]. Generally, the iterative inversion approaches are more robust than the direct inversion approaches [16]. However, these approaches also have some limitations. Firstly, they need to solve the forward problem several times, and hence are computationally expensive [16, 21]. Besides, if the objective function contains multiple minima, it may converge to a wrong solution rather than the true one [12]. Furthermore, the iterative inversion requires a regularization term. And choosing an appropriate regularization parameter is challenging [20]. The purpose of regularization is to suppress the reconstruction noise. The value of the regularization parameter can affect the size of inclusion shown in the reconstructed modulus image and the contrast between the inclusion and background [20]. An appropriate regularization parameter should be used to reduce the noise and preserve the contrast of the reconstructed modulus image simultaneously [20]. However, it is difficult to obtain the optimal regularization parameter in practice.

To the authors’ knowledge, current elasticity inversion schemes require complex equipment to obtain the boundary conditions of displacement or force [23–26]. These complex equipment limit the Young’s modulus reconstruction technique to be used in a freehand ultrasound scan for breast elastography. The purpose of this study is to develop an elasticity reconstruction method that does not need too much manual intervention or measurements of boundary conditions, which can be implemented in common medical ultrasound systems with a freehand scan. To overcome the limitations of the existing inversion approaches, a direct inversion scheme based on finite element method (FEM) is developed. It combines the advantages of direct and iterative inversion approaches. A novel feature of the proposed method is the utilization of the bicubic B-spline function as the tissue motion model to suppress the noise in the inversion, meanwhile it does not need any regularization in the inversion. Besides, a more practical modulus boundary condition is used in this method instead of force or displacement boundary condition. Simulations and phantom experiments are conducted to validate the effectiveness of the proposed method.

## Methods

### Inversion method

In the proposed inversion method of Young’s modulus reconstruction, the tissue was modelled as linearly elastic, isotropic, continuous and incompressible [8, 11, 13–15]. Besides, the commonly used plane-strain approximation is adopted in the proposed method. As the boundary forces of the reconstruction ROI are unknown, the reconstructed Young’s moduli are relative values (i.e., contrast) rather than absolute values for there is no known stress information [12]. The proposed method is comprised of two algorithms. The first algorithm is a bicubic B-spline fitting-based technique which estimates the 2D axial and lateral displacement fields from a pre-estimated axial displacement field, the plane-strain constraints and incompressibility assumption. The second algorithm is FEM-based direct inversion with modulus boundary condition. These algorithms are described in the following sections.

#### Displacement estimation

The displacements within the tissue are first estimated through a two-step optical flow method using ultrasound radiofrequency (RF) signals [27]. It has been proved that the displacements estimated from this method have a rather high accuracy. However, they cannot be directly used in the FEM-based inversion method, since the inversion results are very sensitive to the errors in the displacement measurements. The modulus reconstruction accuracy degrades rapidly when the standard deviation of the displacement error exceeds 10^{−4} mm [11]. However, the displacement errors are typically larger than this threshold [28]. Hence, the estimated displacements need to be processed before the inversion.

*a*

_{ i,j }with size of

*n*

_{ x }

*×n*

_{ y }as axial displacement parameters of bicubic B-spline function, and

*b*

_{ i,j }as lateral displacement parameters with the same size.

*n*

_{ x }and

*n*

_{ y }are the number of uniformly-spaced knots in the lateral and axial directions, respectively. The axial displacement field

*V*and lateral displacement field

*U*can be presented as

where *i* = [*x*/*n*
_{
x
}] − 1, *j* = [*y*/*n*
_{
y
}] − 1, *p* = *x*/*n*
_{
x
} − [*x*/*n*
_{
x
}], *q* = *y*/*n*
_{
y
} − [*y*/*n*
_{
y
}] and where *B*
_{
m
} and *B*
_{
n
} stands for the *m*-th and *n*-th basis function of the B-spline [29, 30], respectively.

*a*

_{ i,j }are calculated by fitting the axial displacements

*V*(x, y) of B-spline and estimated axial displacements

*V*

_{ 0 }(

*x*,

*y*) from optical flow using the least square technique [31]. Hence, the axial displacements estimated with this method could preserve a rather high precision and smoothness. However, the lateral displacements estimated from optical flow are much noisier than axial displacements, due to lower resolution, lower sampling rate and lack of phase information in the lateral direction [32]. It is inappropriate to obtain the lateral displacement parameters

*b*

_{ i,j }in a similar way as the axial parameters. To overcome this limitation, the constraints of 2D plane strain and the incompressibility of the tissue are utilized to compute the lateral displacement. These constraints could be formulated as.

the lateral displacement parameters *b*
_{
i,j
} can also be calculated by the least square technique [31]. Finally, the axial and lateral displacements of B-spline model could be computed by substituting *a*
_{
i,j
} and *b*
_{
i,j
} into Eqs. (1) and (2), respectively.

#### Elasticity reconstruction

where **K**
_{
boundary
}, **d**
_{
boundary
} and **f** are the global stiffness matrix, the global nodal displacement vector and the global nodal force vector of the boundary, respectively. Node displacement **d**
_{
boundary
} can be obtained from the estimated displacement field. The global stiffness matrix **K**
_{
boundary
} is assembled from the elements with the known node position, Young’s modulus and Poisson ratio [34]. The global nodal force vector of boundary **f** determines the forces of the boundary nodes (white circles, Figure 1) and interior nodes (black circles, Figure 1).

**E**is the variable to be solved. Now, consider all the elements within the ROI as an object. Denote the number of elements and nodes as

*N*

_{ element }and

*N*

_{ node }, respectively. As described in [11], the Young’s modulus vector

**E**can be extracted from the multiplication of the global stiffness matrix

**K**

_{ roi }and the displacement vector

**d**. The right term of Equation (5) becomes [11].

**D**is a 2

*N*

_{ node }×

*N*

_{ element }matrix which is related to the elements and Poisson ratio, while the size of

**E**is

*N*

_{ element }× l. For more details about the procedure of assembly of matrix

**D**, the readers are referred to [11]. For the quadrilateral element used in this method,

*N*

_{ node }is greater than

*N*

_{ element }. The Young’s modulus vector

**E**is solved by the least square method

The boundary nodal force vector **f**
_{
boundary
} is a subset of **f** which has been solved in (5).

### Simulations

^{2}(axially × laterally, Figure 2). The inclusions had the same diameter of 6.5 mm. The Young’s moduli of the inclusions and background were 75 kPa and 25 kPa, respectively. Compressional axial strains between 0.5% and 2.0% were applied to the model, with a perfect slip boundary.

A 320-element linear array transducer with a center frequency of 6 MHz and a −3 dB bandwidth of 50% was simulated with Field_II software (http://field-ii.dk/) [35, 36]. The speed of sound was assumed to be 1540 m/s, and the sampling frequency of the RF signals was 32 MHz. The pitch of the transducer was 0.12 mm. The focal depth was 19 mm and the F number was 0.5. These transducer parameters were chosen to correspond to the actual transducer used in the phantom experiments. 320 lines were simulated in each ultrasound image, and the distance between adjacent lines was equal to the pitch. The imaging zone had a depth of 38 mm and a width of 38 mm. Randomly distributed scatterers were simulated within the model. The scatterer density was 200/mm^{2}, which fulfilled the requirement for fully developed speckle in the simulated ultrasound image [37]. The pre-deformed RF signals of the model were simulated from Field_II using the parameters described above. Finite element analysis was performed with FEMLAB 2.3 software (Comsol Inc. Burlington, MA, USA) and MATLAB 6.5 (The MathWorks Inc., Natick, MA, USA) to calculate the theoretical displacements of the model. The scatterers were moved according to finite element solutions with different applied strains. The post-deformed RF signals were then simulated from Field_II using the new scatterers’ distribution. Gaussian white noises were added in the simulated RF signals, with a signal-to-noise ratio (SNR) of 30 dB.

### Phantom experiments

### Quantitative evaluation

In the optical flow-based displacement estimation, a window size of 2.2 × 1.8 mm^{2} (height × width) was used. In the following B-spline fitting, the axial and lateral knot spacing were 4.82 mm and 4.75 mm, respectively. For the modulus reconstruction, the quadrilateral element used in the inversion was 0.48 × 0.48 mm^{2} (height × width). These parameters were used in the data processing of both simulation and phantom studies. The displacement estimation and elasticity reconstruction methods described above were used in both the simulation and phantom studies.

## Results

### Simulation results

^{−3}and 2.6 × 10

^{−3}mm, respectively. The lateral displacements obtained from B-spline fitting (Figure 4(f)) are smoother than those estimated from optical flow (Figure 4(e)), and closer to the theoretical displacement (Figure 4(d)). The standard deviations of lateral displacement errors of optical flow and B-spline fitting are 13.4 × 10

^{−3}and 6.9 × 10

^{−3}mm, respectively. This suggests that the lateral displacements estimated from B-spline fitting are more precise than those from optical flow.

^{2}(height × width). Due to the noise in the lateral displacements estimated from optical flow (Figure 4(e)), the modulus could not be correctly reconstructed (Figure 5(c)). The Young’s moduli of the inclusion reconstructed from both the theoretical and estimated displacements are lower than the theoretical value. The Young’s moduli reconstructed from the theoretical displacement were relatively uniform within the inclusions, with a rapid transition from the inclusion to the background (Figure 5(b)). In contrast, the Young’s moduli reconstructed from the estimated displacement have the highest values in the center of the inclusions and gradually decrease with the distance from the center. A smooth transition from the inclusion to the background is shown (Figure 5(d)).To quantitatively investigate the error of Young’s modulus reconstruction, the mean relative error and RMSE between the theoretical and reconstructed modulus values are calculated with different applied strains (Figure 6). The mean relative errors of reconstructed Young’s modulus from the estimated displacements of B-spline model is −3.43 ± 0.43% and is about 22% smaller than those from the theoretical displacements (−4.40%). The RMSE of the modulus estimates of B-spline displacements is 16.94 ± 0.25% and is about 25% larger than those of theoretical displacements (13.52%). These can be explained by the small errors in the estimated displacements.

### Phantom results

^{2}(height × width, Figure 10(a)), 23.1 × 28.5 mm

^{2}(Figure 10(b)) and 22.2 × 23.7 mm

^{2}(Figure 10(c)), respectively. The numbers of elements used in are 50 × 70, 48 × 60 and 46 × 50, respectively, and the computational time are 2.2 s, 1.2 s and 0.75 s, respectively. All the programs were performed on the platform of MATLAB 6.5 (The MathWorks Inc., Natick, MA, USA). The reconstructed modulus values of a circular ROI in the inclusion center with a radius of 2 mm are 1.87 ± 0.16, 1.98 ± 0.17 and 2.02 ± 0.18 for Figures 10(a), 10(b) and 10(c), respectively. Although the reconstruction ROIs are different in size, the estimated Young’s modulus distributions of the inclusion are similar.

## Discussion

The goal of this study is to develop a robust elasticity reconstruction method for ultrasound elastography with freehand scan. Conventional elastography typically shows axial strain images [3, 38], which suffer from mechanical artifacts (such as target hardening and stress concentration) and lack of quantitative information [7, 8]. The reconstruction of Young’s modulus or shear modulus (i.e., solving the inverse problem) can overcome these drawbacks. However, the inverse problem generally has an ill-posed nature [9, 25, 33]. In this work, we proposed a B-spline function-based displacement estimation method, and employed a modulus boundary condition to convert the ill-posed problem to a well-posed one. Simulations and phantom experiments were conducted to validate the performance of the proposed method. In the simulation study, the mean relative error and RMSE were used to evaluate the accuracy and precision of Young’s modulus reconstruction. In the phantom study, the reconstructed modulus images were compared with the corresponding strain images.

The axial and lateral displacement fields used in the inversion scheme are defined as bicubic B-spline function. This tissue motion model is also called free-form deformation [39], and has been used in ultrasound cardiac motion estimation [40, 41]. The conventional motion estimation methods are usually based on local information (Figures 4(b) and (e)), while the proposed B-spline fitting method takes the global information into consideration, and can obtain much smoother displacement estimation (Figures 4(c) and (f)). Although the axial displacement estimated from the optical method is more precise than that of the B-spline fitting method, they cannot be directly used in the inversion since the noise in the displacement will be amplified without the restriction of continuity. The B-spline displacement model has the implicit continuous constraint, and it is helpful to reduce the ill-condition of the Young’s modulus inversion. The tissue incompressibility and plane strain assumptions are used to estimate the lateral displacement, which is a necessary part for modulus reconstruction. The B-spline motion model could be implemented by the proposed displacement fitting or non-rigid registration of the RF signals [42]. However, the non-rigid registration is much more computationally expensive than displacement fitting, and may present local minima during the iteration process, due to the highly oscillatory nature of the RF signals [20]. Thus, the displacement fitting technique is used in this study.

Conventional ultrasound elastography typically shows strain images, and can provide improved ability to determine the lesions’ location and shape when compared to the corresponding B-mode image [3]. Nevertheless, under the assumption of stress uniformity, the strain image suffers from several artifacts and therefore do not reveal the true tissue elastic properties [8, 9]. In reality, the stress together with strain attenuates with depth, hence the tissue on the bottom part of the tissue seems stiffer than that on the top (Figures 7(c)-(h)). These are the so-called target hardening artifacts [5, 8, 9]. In the Young’s modulus images, these artifacts are eliminated (Figures 7(i)-(l), and Figures 9(c)-(d)), which may facilitate the diagnosis of the breast tumors in the clinic. It would have a high contrast of modulus image for the inclusion by fixing the scale of the image colour bar around 2, since the modulus contrast between normal and cancerous tissues is usually great than 100% [43].

The boundary conditions currently used in the iterative inversion schemes include displacement boundary conditions [13, 25], force boundary conditions [8] and a combination of both [16, 20, 23]. However, the force distribution within the tissue is difficult to measure *in vivo*
[7]. The combination boundary conditions of displacement and force need special experimental setup [20, 23], and hence limit the usage in the clinic. The boundary condition used in the proposed method is the modulus boundary condition which assumes the moduli of the boundary elements around the reconstruction ROI are the same [13, 15, 44]. It may induce some artifacts if the moduli on the ROI boundary are different. However, the ROI around the inclusion could be chosen according to the B-mode and/or strain images in order to avoid this situation. Besides, the size or shape of the ROI can be changed. If the moduli of the inclusion do not have much difference when different ROIs are used, as shown in Figure 10, the modulus reconstruction result can be deemed reliable.

In the simulation study, the RMSE of the Young’s moduli reconstructed from the estimated displacements are about 25% higher than those from theoretical displacements, while the mean error of estimated Young’s moduli reconstructed from the estimated displacements are about 22% lower than those from the theoretical displacements (Figure 6). These results suggest that the precision of the Young’s modulus reconstructed by the proposed method is comparable to those from the theoretical displacements. As shown, the modulus distributions reconstructed from the theoretical displacements are not the same as the theoretical modulus distribution (Figures 5(a) and (c)). This can be explained by the measurement noise (i.e., displacement estimation noise) and process noise caused by the mismatch between the models used in the forward and inverse problems, which are the two main factors that affect the modulus reconstruction results [45]. In the forward problem, the object was meshed with non-uniform size elements, and the meshes were finer near the inclusion boundary (Figure 2). In the inverse problem, as the internal geometry of the tissue is usually unknown, the uniform quadrilateral elements were used (Figure 1). The discrepancy of the meshes used in the forward and inverse problems leads to reconstruction error even if the theoretical displacements are used as the input. The reconstruction errors in the simulation study are similar to the results of Zhu et al. [11].

Another limitation of the proposed method is the smoothing effect on the interface between the inclusion and background. This is mainly attribute to the smoothness of displacements by the bicubic B-spline fitting in the displacement estimation. The high-frequency estimation noise together with high-frequency displacement information is filtered out, and hence the high-frequency components in the modulus distribution are also eliminated. The smoothness of displacements used in the proposed method is similar to the regularization used in the iterative inversion schemes, which also tends to smooth over the sharp dips in the reconstructed modulus image [2, 16, 24, 25]. The lower reconstructed Young’s modulus values in the inclusions (or contrasts) are consistent with the results of a previous study [24]. As shown in Figure 7, the smaller the size of the inclusion is, the lager the influence of the smoothing effect will be. Hence, the result of the reconstructed Young’s modulus becomes worse when the inclusion gets smaller (Figure 7). To overcome these limitations, the mesh adaptation algorithm has been proposed to improve the accuracy of Young’s modulus reconstruction [46]. Nevertheless, the adaptive mesh generation requires image segmentation using the internal tissue geometry or axial strain information.

Tissue deformation is typically 3D in nature, but ultrasound RF data are often collected from a 2D imaging plane. The 2D plane strain simplification is utilized in the proposed method. Some errors may be incurred when the 3D elasticity problem is modelled by either plane strain or plane stress approximation [7, 20]. However, we believe that plane strain is a reasonable assumption for breast elastography in the clinic because it has proved that the motion perpendicular to the image plane is small in the breast *in vivo*
[11]. And when 3D data are available, the FEM-based inversion scheme could be extended to 3D case. As a final note, the displacement fields could also be estimated from other imaging modalities such as magnetic resonance imaging (MRI), and the Young’s modulus distribution of tissue could be inferred by using the proposed method.

## Conclusions

In this paper, a finite element based Young’s modulus reconstruction method is proposed to facilitate the quantitative investigation of the elasticity of soft tissues. The B-spline function fitting method is proposed to reduce the displacement estimation noise. Then the Young’s modulus distribution within the tissue is reconstructed based on the framework of FEM with the input of displacements and modulus boundary condition. The simulation and phantom studies demonstrated that the proposed method yields the inclusion’s relative modulus with reasonably high accuracy, and correctly delineates the locations and sizes of inclusions. Besides, the computational time of the proposed reconstruction method ranges from 0.75 to 2.2 s, which is much faster than the conventional iterative reconstruction methods (e.g., 18 hours) [21]. More importantly, the proposed method does not need any regularization or boundary force information, thus avoiding too much manual intervention and special equipment. With these advantages, the proposed method could be implemented in conventional ultrasound systems with freehand scan.

## Declarations

### Acknowledgements

This work was supported in part by the National Basic Research Program of China (973) under Grant No. 2011CB707701, National Major Scientific Instrument and Equipment Development Project under Grant No. 2011YQ03011402, National Natural Science Foundation of China under Grant Nos. 61271131, 61322101, 81471665 and 81227901, National Key Technology R&D Program of China under Grant No. 2011BAI02B03, Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant No. 20130002110061, Tsinghua University Initiative Scientific Research Program, Tsinghua National Laboratory for Information Science and Technology (TNList) Cross-discipline Foundation, and Beijing Science and Technology Plan Project under Grant No. Z131100005213001. The authors would like to thank Dr. Lingyun Huang, with Philips Research China, Shanghai, China, for his assistance in data acquisition and helpful discussions.

## Authors’ Affiliations

## References

- Siegel R, Naishadham D, Jemal A:
**Cancer statistics, 2013.***CA Cancer J Clin*2013,**63**(1):11–30. 10.3322/caac.21166View ArticleGoogle Scholar - Goenezen S, Dord J-F, Sink Z, Barbone PE, Jiang J, Hall TJ, Oberai AA:
**Linear and nonlinear elastic modulus imaging: an application to breast cancer diagnosis.***IEEE Trans Med Imaging*2012,**31**(8):1628–1637.View ArticleGoogle Scholar - Hiltawsky KM, Krüger M, Starke C, Heuser L, Ermert H, Jensen A:
**Freehand ultrasound elastography of breast lesions: clinical results.***Ultrasound Med Biol*2001,**27**(11):1461–1469. 10.1016/S0301-5629(01)00434-3View ArticleGoogle Scholar - Ophir J, Alam SK, Garra B, Kallel F, Konofagou E, Krouskop T, Varghese T:
**Elastography: Ultrasonic estimation and imaging of the elastic properties of tissues.***Proc Inst Mech Eng H J Eng Med*1999,**213**(3):203–233. 10.1243/0954411991534933View ArticleGoogle Scholar - 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. 10.1177/016173469101300201View ArticleGoogle Scholar - Malvern LE:
*Introduction to the Mechanics of a Continuous Medium*. 1969.Google Scholar - Doyley MM:
**Model-based elastography: a survey of approaches to the inverse elasticity problem.***Phys Med Biol*2012,**57**(3):R35-R73. 10.1088/0031-9155/57/3/R35View ArticleGoogle Scholar - Kallel F, Bertrand M:
**Tissue elasticity reconstruction using linear perturbation method.***IEEE Trans Med Imaging*1996,**15**(3):299–313. 10.1109/42.500139View ArticleGoogle Scholar - Doyley MM, Srinivasan S, Pendergrass SA, Wu Z, Ophir J:
**Comparative evaluation of strain-based and model-based modulus elastography.***Ultrasound Med Biol*2005,**31**(6):787–802. 10.1016/j.ultrasmedbio.2005.02.005View ArticleGoogle Scholar - Oberai AA, Gokhale NH, Doyley MM, Bamber JC:
**Evaluation of the adjoint equation based algorithm for elasticity imaging.***Phys Med Biol*2004,**49**(13):2955–2974. 10.1088/0031-9155/49/13/013View ArticleGoogle Scholar - Zhu Y, Hall TJ, Jiang J:
**A finite-element approach for Young’s modulus reconstruction.***IEEE Trans Med Imaging*2003,**22**(7):890–901. 10.1109/TMI.2003.815065View ArticleGoogle Scholar - Doyley M, Meaney P, Bamber J:
**Evaluation of an iterative reconstruction method for quantitative elastography.***Phys Med Biol*2000,**45**(6):1521–1540. 10.1088/0031-9155/45/6/309View ArticleGoogle Scholar - Bishop J, Samani A, Sciarretta J, Plewes DB:
**Two-dimensional MR elastography with linear inversion reconstruction: methodology and noise analysis.***Phys Med Biol*2000,**45**(8):2081–2091. 10.1088/0031-9155/45/8/302View ArticleGoogle Scholar - Sumi C, Suzuki A, Nakayama K:
**Estimation of shear modulus distribution in soft tissue from strain distribution.***IEEE Trans Biomed Eng*1995,**42**(2):193–202. 10.1109/10.341832View ArticleGoogle Scholar - Skovoroda A, Emelianov S, O’Donnell M:
**Tissue elasticity reconstruction based on ultrasonic displacement and strain images.***IEEE Trans Ultrason Ferroelectr Freq Control*1995,**42**(4):747–765.View ArticleGoogle Scholar - Oberai AA, Gokhale NH, Feijóo GR:
**Solution of inverse problems in elasticity imaging using the adjoint method.***Inverse Probl*2003,**19**(2):297–313. 10.1088/0266-5611/19/2/304View ArticleGoogle Scholar - Floc’h L, Ohayon J, Tracqui P, Finet G, Gharib AM, Maurice RL, Cloutier G, Pettigrew RI:
**Vulnerable atherosclerotic plaque elasticity reconstruction based on a segmentation-driven optimization procedure using strain measurements: theoretical framework.***IEEE Trans Med Imaging*2009,**28**(7):1126–1137.View ArticleGoogle Scholar - Raghavan K, Yagle AE:
**Forward and inverse problems in elasticity imaging of soft tissues.***IEEE Trans Nucl Sci*1994,**41**(4):1639–1648. 10.1109/23.322961View ArticleGoogle Scholar - Mousavi SR, Sadeghi-Naini A, Czarnota GJ, Samani A:
**Towards clinical prostate ultrasound elastography using full inversion approach.***Med Phys*2014,**41**(3):1–12.View ArticleGoogle Scholar - Richards MS, Barbone PE, Oberai AA:
**Quantitative three-dimensional elasticity imaging from quasi-static deformation: a phantom study.***Phys Med Biol*2009,**54**(3):757–779. 10.1088/0031-9155/54/3/019View ArticleGoogle Scholar - Samani A, Bishop J, Plewes DB:
**A constrained modulus reconstruction technique for breast cancer assessment.***IEEE Trans Med Imaging*2001,**20**(9):877–885. 10.1109/42.952726View ArticleGoogle Scholar - Doyley MM, Parker KJ:
**Elastography: General Principles and Clinical Applications.***Ultrasound Clinics*2014,**9**(1):1–11. 10.1016/j.cult.2013.09.006View ArticleGoogle Scholar - Karimi H, Fenster A, Samani A:
**A novel fast full inversion based breast ultrasound elastography technique.***Phys Med Biol*2013,**58**(7):2219–2233. 10.1088/0031-9155/58/7/2219View ArticleGoogle Scholar - Richards M, Doyley M:
**Investigating the impact of spatial priors on the performance of model-based IVUS elastography.***Phys Med Biol*2011,**56**(22):7223–7246. 10.1088/0031-9155/56/22/014View ArticleGoogle Scholar - Jiang J, Varghese T, Brace CL, Madsen EL, Hall TJ, Bharat S, Hobson MA, Zagzebski JA, Lee FT:
**Young’s modulus reconstruction for radio-frequency ablation electrode-induced displacement fields: a feasibility study.***IEEE Trans Med Imaging*2009,**28**(8):1325–1334.View ArticleGoogle Scholar - Plewes DB, Bishop J, Samani A, Sciarretta J:
**Visualization and quantification of breast cancer biomechanical properties with magnetic resonance elastography.***Phys Med Biol*2000,**45**(6):1591–1610. 10.1088/0031-9155/45/6/314View ArticleGoogle Scholar - Pan X, Gao J, Tao S, Liu K, Bai J, Luo J:
**A two-step optical flow method for strain estimation in elastography: Simulation and phantom study.***Ultrasonics*2014,**54**(4):990–996. 10.1016/j.ultras.2013.11.010View ArticleGoogle Scholar - Lopata RG, Nillesen MM, Hansen HH, Gerrits IH, Thijssen JM, de Korte CL:
**Performance evaluation of methods for two-dimensional displacement and strain estimation using ultrasound radio frequency data.***Ultrasound Med Biol*2009,**35**(5):796–812. 10.1016/j.ultrasmedbio.2008.11.002View ArticleGoogle Scholar - Lee S, Wolberg G, Shin SY:
**Scattered data interpolation with multilevel B-splines.***IEEE Trans Vis Comput Graph*1997,**3**(3):228–244. 10.1109/2945.620490View ArticleGoogle Scholar - Lee S, Wolberg G, Chwa KY, Shin SY:
**Image metamorphosis with scattered feature constraints.***IEEE Trans Vis Comput Graph*1996,**2**(4):337–354. 10.1109/2945.556502View ArticleGoogle Scholar - Press WH, Teukolsky SA, Vetterling WT, Flannery BP:
*Numerical Recipes: The art of Scientific Computing.*3rd edition. New York: Cambridge University Press; 2007:776–780.Google Scholar - Luo J, Konofagou EE:
**Effects of various parameters on lateral displacement estimation in ultrasound elastography.***Ultrasound Med Biol*2009,**35**(8):1352–1366. 10.1016/j.ultrasmedbio.2009.03.001View ArticleGoogle Scholar - Barbone PE, Bamber JC:
**Quantitative elasticity imaging: what can and cannot be inferred from strain images.***Phys Med Biol*2002,**47**(12):2147–2164. 10.1088/0031-9155/47/12/310View ArticleGoogle Scholar - Bang H, Kwon YW:
*The Finite Element Method Using MATLAB.*2nd edition. Edited by: Frank K. CRC Press; 2000:51–56.Google Scholar - Jensen JA, Svendsen NB:
**Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers.***EEE Trans Ultrason Ferroelectr Freq Control*1992,**39**(2):262–267.View ArticleGoogle Scholar - Jensen JA:
**Field: A program for simulating ultrasound systems.***Med Biol Eng Comput*1996,**34:**351–352. 10.1007/BF02520003View ArticleGoogle Scholar - Oosterveld B, Thijssen J, Verhoef W:
**Texture of B-mode echograms: 3-D simulations and experiments of the effects of diffraction and scatterer density.***Ultrason Imaging*1985,**7**(2):142–160. 10.1177/016173468500700204View ArticleGoogle Scholar - Garra BS, Cespedes EI, Ophir J, Spratt SR, Zuurbier RA, Magnant CM, Pennanen MF:
**Elastography of breast lesions: initial clinical results.***Radiology*1997,**202**(1):79–86. 10.1148/radiology.202.1.8988195View ArticleGoogle Scholar - Rueckert D, Sonoda LI, Hayes C, Hill DLG, Leach MO, Hawkes DJ:
**Nonrigid registration using free-form deformations: application to breast MR images.***IEEE Trans Med Imaging*1999,**18**(8):712–721. 10.1109/42.796284View ArticleGoogle Scholar - Heyde B, Jasaityte R, Barbosa D, Robesyn V, Bouchez S, Wouters P, Maes F, Claus P, D’Hooge J:
**Elastic image registration versus speckle tracking for 2-D myocardial motion estimation: a direct comparison in vivo.***IEEE Trans Med Imaging*2013,**32**(2):449–459.View ArticleGoogle Scholar - Ledesma-Carbayo MJ, Kybic J, Desco M, Santos A, Suhling M, Hunziker P, Unser M:
**Spatio-temporal nonrigid registration for ultrasound cardiac motion estimation.***IEEE Trans Med Imaging*2005,**24**(9):1113–1126.View ArticleGoogle Scholar - Xiaochang P, Jing G, Jinhua S, Jianwen L, Jing B:
**A regularization-free Young’s modulus reconstruction algorithm for ultrasound elasticity imaging.***Engineering in Medicine and Biology Society (EMBC), 2013 35th Annual International Conference of the IEEE; 3–7 July 2013*2013, 1132–1135.View ArticleGoogle Scholar - Krouskop TA, Wheeler TM, Kallel F, Garra BS, Hall T:
**Elastic moduli of breast and prostate tissues under compression.***Ultrason Imaging*1998,**20**(4):260–274. 10.1177/016173469802000403View ArticleGoogle Scholar - Steele DD, Chenevert TL, Skovoroda AR, Emelianov SY:
**Three-dimensional static displacement, stimulated echo NMR elasticity imaging.***Phys Med Biol*2000,**45**(6):1633–1648. 10.1088/0031-9155/45/6/316View ArticleGoogle Scholar - Lu M, Zhang H, Wang J, Yuan J, Hu Z, Liu H:
**Reconstruction of elasticity: a stochastic model-based approach in ultrasound elastography.***Biomed Eng Online*2013,**12**(1):79. 10.1186/1475-925X-12-79View ArticleGoogle Scholar - Goksel O, Eskandari H, Salcudean SE:
**Mesh adaptation for improving elasticity reconstruction using the FEM inverse problem.***IEEE Trans Med Imaging*2013,**32**(2):408–418.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.