Next Article in Journal
On the Passivity of Hysteretic Systems with Double Hysteretic Loops
Next Article in Special Issue
Effect of Different Manufacturing Methods on the Conflict between Porosity and Mechanical Properties of Spiral and Porous Polyethylene Terephthalate/Sodium Alginate Bone Scaffolds
Previous Article in Journal
Removal of Chromium(VI) from Aqueous Solutions Using Fe3O4 Magnetic Polymer Microspheres Functionalized with Amino Groups
Previous Article in Special Issue
Reduction of Adipose Tissue Formation by the Controlled Release of BMP-2 Using a Hydroxyapatite-Coated Collagen Carrier System for Sinus-Augmentation/Extraction-Socket Grafting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultrasound Elasticity Imaging System with Chirp-Coded Excitation for Assessing Biomechanical Properties of Elasticity Phantom

1
Department of Biomedical Engineering, National Cheng Kung University, Tainan City 70101, Taiwan
2
Division of Infectious Diseases, Department of Medicine of Chi Mei Medical Center, Tainan City 71004, Taiwan
3
Department of Electrical Engineering, Southern Taiwan University of Science and Technology, Tainan City 71005, Taiwan
*
Author to whom correspondence should be addressed.
Materials 2015, 8(12), 8392-8413; https://doi.org/10.3390/ma8125458
Submission received: 31 August 2015 / Revised: 3 November 2015 / Accepted: 24 November 2015 / Published: 3 December 2015
(This article belongs to the Special Issue Selected Papers from ICBEI2015)

Abstract

:
The biomechanical properties of soft tissues vary with pathological phenomenon. Ultrasound elasticity imaging is a noninvasive method used to analyze the local biomechanical properties of soft tissues in clinical diagnosis. However, the echo signal-to-noise ratio (eSNR) is diminished because of the attenuation of ultrasonic energy by soft tissues. Therefore, to improve the quality of elastography, the eSNR and depth of ultrasound penetration must be increased using chirp-coded excitation. Moreover, the low axial resolution of ultrasound images generated by a chirp-coded pulse must be increased using an appropriate compression filter. The main aim of this study is to develop an ultrasound elasticity imaging system with chirp-coded excitation using a Tukey window for assessing the biomechanical properties of soft tissues. In this study, we propose an ultrasound elasticity imaging system equipped with a 7.5-MHz single-element transducer and polymethylpentene compression plate to measure strains in soft tissues. Soft tissue strains were analyzed using cross correlation (CC) and absolution difference (AD) algorithms. The optimal parameters of CC and AD algorithms used for the ultrasound elasticity imaging system with chirp-coded excitation were determined by measuring the elastographic signal-to-noise ratio (SNRe) of a homogeneous phantom. Moreover, chirp-coded excitation and short pulse excitation were used to measure the elasticity properties of the phantom. The elastographic qualities of the tissue-mimicking phantom were assessed in terms of Young’s modulus and elastographic contrast-to-noise ratio (CNRe). The results show that the developed ultrasound elasticity imaging system with chirp-coded excitation modulated by a Tukey window can acquire accurate, high-quality elastography images.

1. Introduction

The standard medical practice of soft tissue palpation is based on the qualitative assessment of stiffness of tissues. Human tissue lesions are generally correlated with the variation in elastic properties of tissues. In many cases, despite the difference in stiffness between a lesion and normal tissue, the detection and evaluation of a pathological lesion via palpation is difficult owing to its small size or its location. In general, lesions may or may not exhibit sonographic contrast, which would enable them to be ultrasonically detectable. For instance, tumors of the prostate or breast may be significantly stiffer than the embedding tissue, and yet be invisible or barely visible in standard ultrasound examinations. Diffuse diseases such as cirrhosis of the liver are known to significantly increase the stiffness of liver tissue; however, the tissue may appear normal in conventional ultrasound examination [1].
Therefore, we can obtain new information related to the biomechanical properties of tissues for differentiating normal tissues from abnormal tissues by imaging tissue stiffness or a related parameter, such as strain under stress. In 1991, Ophir et al. [1] proposed a method to quantify strain information in biological tissues. A method capable of quantitatively imaging the stiffness of tissues with good resolution, sensitivity, and diminished speckle was presented in [2]. Ultrasound elastography can generate numerous types of images referred to as elastograms. Ultrasound elastography is typically used to evaluate a mass as benign or malignant for lesions. Soft and hard materials exhibit large and low strain values, respectively, and large and low strain values are typically displayed as bright and dark regions, respectively. Ultrasound elastography is being increasingly used for assessing the biomechanical properties of tissues and has been applied to numerous organs and pathologies including many cancers, such as scirrhous carcinoma of the breast, liver cancer, and prostatic carcinoma [3,4,5]. The different types of ultrasound elastography methods include compression ultrasound elastography, acoustic radiation force impulse (ARFI) technique, and shear wave elastography [5,6,7].
Compression ultrasound elastography is the most commonly used method in which a constant stress is applied to the studied tissue [5,7,8]. The technique uses manual compression or a controlled stepper motor to move a transducer to generate stress on the tissue and measures tissue deformation to estimate the elasticity of the tissue. Displacement is estimated by comparing echoes before and after compression by correlation methods [9]. The measured displacement, generated strain (ε), and stress can be used to compute the Young’s modulus, which is a more objective parameter of stiffness. This strain map is often referred to as an elastogram because the applied stress is unknown and only strain is displayed. Young’s modulus is the ratio of stress to strain within the elastic limit. Nightingale proposed the ARFI method in which acoustic radiation force is employed in place of transducer compression to generate force [7]. Acoustic radiation force is a unidirectional force applied to absorbing or reflecting targets in the propagation path of an acoustic wave. Attenuation is a frequency-dependent phenomenon and is primarily caused by absorption in soft tissues. The momentum transfer from the acoustic wave generates a force that causes displacement of the tissue [10,11]. Sarvazyan et al. [12] proposed a method, which can be considered to be a precursor of elastography techniques, based on ultrasonic radiation pressure by combining radiation pressure or acoustic radiation force and the shear waves generated. This technique is referred to as shear wave elasticity imaging (SWEI). Thereafter, Nightingale et al. validated the clinical feasibility of radiation force induced shear-wave imaging [13]. When acoustic radiation force is applied to a given spatial volume for a short duration, transient shear waves that propagate away from the initial region of excitation and perpendicular to the push pulse are generated. In tissues, shear waves travel at a velocity of approximately 1–10 m/s and can be easily tracked using ultrasound [10].
Ultrasound elastography is being increasingly used by clinicians to detect lesions or cancers in soft tissues. However, poor quality of elasticity imaging may degrade its value in clinical applications [14]. Imaging quality is determined by the visibility of small lesions and is limited by decorrelation noise [15]. Elastographic signal-to-noise ratio (SNRe), one of the major quality metrics of elasticity imaging, is affected by numerous factors such as echo signal-to-noise ratio (eSNR), applied strain, attenuation, penetration depth, and axial resolution. In an ultrasound system, the energy of the trigger signal affects eSNR and penetration depth, whereas the bandwidth of the signal influences axial resolution. A low eSNR reduces the uniformity between continuous echo frames acquired to estimate correlation-based displacement. Therefore, large displacement errors are generated and appear in strain images as a high-intensity noise; this is referred to as the decorrelation noise. In this situation, longer correlation window lengths reduce decorrelation noise for small tissue deformations at the cost of slightly-reduced axial resolution [16]. Therefore, the performance of elasticity imaging can be improved by increasing the energy of the trigger signal. In ultrasound imaging systems, the energy of the signal can be increased in two ways. First, when the peak power of the trigger signal is increased, commercial ultrasound systems increase peak power of trigger signal to the maximum value that is close to the mechanical index (MI) specification. Second, the average power of the trigger signal can be increased by lengthening the duration of the trigger. Therefore, a coded excitation signal can increase the average power of the trigger signal without affecting the power amplitude to realize a better eSNR.
Since the elastic properties of soft tissues are an important diagnostic indicator in clinics, palpation is one of the most commonly used diagnostic methods. However, palpation is limited to the inspection of superficial tissue, and the diagnosis is subjective with low sensitivity. The disadvantages of assessing tissue elasticity via palpation can be overcome by using ultrasound elastography. However, the eSNR is diminished because of the attenuation of ultrasonic energy by soft tissues, thus resulting in reduction of elastography qualities such as SNRe and elastographic contrast-to-noise ratio (CNRe). Chirp-coded pulse excitation has been proposed for improving the SNRe and CNRe of ultrasound elastography. However, the elongated chirp-coded pulses decrease the axial resolution of an ultrasound image. Therefore, pulse compression technique can be employed to improve the axial resolution of images. However, the effects of chirp-coded pulse modulated with different window functions and different strain analysis algorithms on elastography qualities are still not clearly known.
The aim of this study is to develop an ultrasound elasticity imaging system with chirp-coded excitation for assessing the biomechanical properties of tissue-mimicking phantom. Further, we analyzed the correlation between the qualities of elastography and a chirp-coded pulse modulated with different window functions and the effects of different strain analysis algorithms on the qualities of elastography.

2. System Architecture

In this study, an ultrasound elastography system was developed for quantitatively assessing the elasticity properties of tissues and phantoms. The schematic diagram of the proposed ultrasound elastography system with two types of excitation signals is shown in Figure 1. First, a 7.5-MHz single-element focused transducer (Model V320; Panametrics, Waltham, MA, USA) was excited by an ultrasonic pulser (Model 5058PR; Panametrics). Subsequently, a chirp waveform was generated using the MATLAB program and used as the trigger signal input for the arbitrary waveform generator (AWG; Model AFG3252; Tektronix Inc., Beaverton, OR, USA) connected to a radio-frequency power amplifier (Model 25A250A; Amplifier Research, Souderton, PA, USA). The echoes were then amplified by the receiver (Model 5073PR; Panametrics).
Figure 1. Block diagram of ultrasound elasticity imaging system.
Figure 1. Block diagram of ultrasound elasticity imaging system.
Materials 08 05458 g001
The received signals were digitized by an 8-bit analog-to-digital converter (Model PCI5152; National Instruments, Austin, TX, USA) at 2 GS/s sampling rate housed in a personal computer. The LabVIEW program was used to control the motor, acquire the data, and process the signals. Additionally, a three-axis step motor was used to move the ultrasonic probe and the indentation system containing the load cell (Model SMT S-Type; Interface Inc., Scottsdale, AZ, USA) and the compression plate (Polymethylpentene (PMP/TPX; Plastics Industry Development Center, Tainan City, Taiwan)), as shown in Figure 2a. The load cell was connected in series with a compression plate to record the corresponding force response. The indentation system was combined with a fixture and fixed on a Z-axis motor, and the ultrasonic probe was fixed on an X-axis motor, as shown in Figure 2b.
Figure 2. (a) Load cell and compression plate and (b) indentation system.
Figure 2. (a) Load cell and compression plate and (b) indentation system.
Materials 08 05458 g002

3. Methodology Description

The location of indentation system was adjusted by motor control interface which ensured the compression plate to press the object lightly. The position of the probe was controlled by step motor system for lateral scanning to construct two-dimensional B-mode image, as shown in Figure 3. A pre-compression radio frequency (RF)-signal was acquired during uniaxial loading and unloading of the phantom, respectively. Deformations of 1% of the phantom initial height were employed [2,17,18]. Finally, the software of MATLAB R2009a (The MathWorks, Natick, MA, USA) was used to process signal, and the elastic information of the object can be obtained.
Figure 3. Measurement method.
Figure 3. Measurement method.
Materials 08 05458 g003

3.1. Chirp Pulse

A typical chirp signal is a linear frequency modulation waveform that means its instantaneous frequency is known as linear change with time [19], and linearly spans a frequency bandwidth B = f2f1, where f1 and f2 are the starting and ending frequencies, respectively. If the chirp sweeps from f1 to f2 over a time, T, then chirp-code excitation is described by formula [14]:
C ( t ) = w ( t ) cos 2 π ( f 1 + B 2 T t ) t ,           0 < t < T
where w(t) is a tapering window function. Signal duration is adjusted by compression ratio. When compression ratio increases, the signal duration is elongated. On the other hand, when compression ratio decreases, the signal duration is shortened. There are many window functions have be used for code pulse excitation. Such as Tukey, Hanning, Gaussian [14,15,20] et al., in order to make the amplitude of the waveform is more prominent but the trigger signal maintains the same bandwidth. The window functions can be used as tapering or filtering window functions and many window functions have several parameters, it is nearly impossible to test all window functions to find the optimal chirp scheme for strain imaging.
The chirp pulse was programed using the MATLAB program and input to the AWG that set the center frequency to be 7.5 MHz, as shown in Figure 4. In order to investigate the effect of window function and compression filter on strain imaging quality, two chirp schemes were designed, as shown in Table 1.
Table 1. Chirp schemes.
Table 1. Chirp schemes.
Name of SchemeWindow FunctionType of Compression Filter
Chirp 1TukeyMatched
Chirp 2GaussianMatched
Figure 4. Chirp without window (w(t) = 1).
Figure 4. Chirp without window (w(t) = 1).
Materials 08 05458 g004

3.1.1. Tukey Window

The Tukey window, often referred to as the cosine-tapered window, is shown in Figure 5. The window is defined as:
w ( t ) { 1.0 , 0 | t |     β T 2 1 2 ( 1 + cos ( π   | t | −  β ( T / 2 ) 2 ( 1 β ) ( T / 2 ) ) ) , β T 2   | t |   T 2
where T is the window length, and β ∈ (0, 1) such that β = 0 is a rectangular window and β = 1 is a Hanning window. Liu et al. reported that a chirp pulse excitation with a 40% Tukey window and matched compression filter exhibits better strain imaging performance [20].
Figure 5. (a) Tukey window function; (b) chirp with Tukey window; and (c) frequency responses.
Figure 5. (a) Tukey window function; (b) chirp with Tukey window; and (c) frequency responses.
Materials 08 05458 g005

3.1.2. Gaussian Window

In mathematics, a Gaussian function, shown in Figure 6 and often simply referred to as a Gaussian, is a function of the form:
w ( t ) = e 1 2 ( α 2 t T 1 ) 2
Figure 6. (a) Gaussian window function; (b) chirp with Gaussian window; and (c) frequency responses.
Figure 6. (a) Gaussian window function; (b) chirp with Gaussian window; and (c) frequency responses.
Materials 08 05458 g006

3.1.3. Matched Filter

A matched filter is a common compression filter that eliminates phase difference to achieve compressed signals in response to extend pulse duration in the time domain. A matched filter can be defined as follows [21]:
h ( n ) = s * n
From Equation (3), the conjugate complex of s(n) shows that the phase of the matched filter is contrary to that of an echo signal, *, complex conjugate. Therefore, the phase difference can be eliminated by computing the convolution of h(n) and s(n), thus focusing the echo signal energy to complete pulse compression , shown as Figure 7.
Figure 7. (a) Frequency responses before matched filter; and (b) frequency responses of Gaussian and Tukey after matched filter.
Figure 7. (a) Frequency responses before matched filter; and (b) frequency responses of Gaussian and Tukey after matched filter.
Materials 08 05458 g007

3.2. Tissue-Mimicking Phantom

Timothy et al. employed agar- and gelatin-based materials for elastography experimental validation [22]. The Young’s modulus values of these agar and gelatin dispersions are in the range of 5–135 kPa, a range representative of soft tissue, while the acoustic velocity and attenuation values are in the range of 1492–1575 m/s and 0.1–0.52 dB/(cm·MHz), respectively, again representative of soft tissue [23].
The block-shaped phantom used in this study was fabricated using a mold. A rod with 6 mm diameter was placed in a cubic mold with 40 mm height, 40 mm width, and 80 mm length, as shown in Figure 8. The phantom was fabricated using two different types of tissues to mimic normal and tumor tissues. In this study, the background material is an agar-gelatin mixture containing 7 wt % gelatin powder (SI-G2500, Uni-Onward Corp., New Taipei City, Taiwan) and 1 wt % agar powder (FLU-05093, Uni-Onward Corp.) dissolved in hot water. Tissue scattering was mimicked by graphite powder (AL-282863, Uni-Onward Corp.) of which 3 wt % was added to the gelatin/agar solution. A cylindrical tumor was fabricated using a mixture of 7 wt % gelatin powder, 2.5 wt % agar powder, and 3 wt % graphite powder. The agar-gelatin materials were produced according to the method described by Lopata et al. [22,24].
Figure 8. (a) A mold for constructing the block shaped phantom; (b) Phantom.
Figure 8. (a) A mold for constructing the block shaped phantom; (b) Phantom.
Materials 08 05458 g008

3.3. Signal Processor

Objects with different elasticities exhibit different deformations under the same compression force. The ultrasound waves emitted from the transducer were used to measure the original tissue thickness and the variation in thickness of the soft tissue layer in the indented site on the basis of time-of-flight and sound speed. The average speed of sound in soft tissues of the human body is 1540 m/s. For instance, let us assume that the boundaries A, B are located on the path of the ultrasound waves emitted from the ultrasound transducer F, as shown in Figure 9. The distance between A and B is A B ¯ and A B ¯ = F A ¯ F A ¯ . Let us assume that the boundary A moves to A′ and B moves to B′ when pressure is exerted on the path of the ultrasound waves emitted from the transducer. The distance between A′ and B′ is A B ¯ and A B ¯ = A B ¯ A B ¯ . The waveform of echo signal varies when the object undergoes deformation under certain compression, as shown in Figure 10. Congruent A-lines were divided into temporal segments, and the corresponding time shifts of the segments were measured using cross correlation (CC) and absolution difference (AD) algorithms. Moreover, the Kalman filter was employed to estimate statistical noise such as white noise, as shown in Figure 11.
Figure 9. The schematic diagram of tissue before and after compression. RF: radio frequency.
Figure 9. The schematic diagram of tissue before and after compression. RF: radio frequency.
Materials 08 05458 g009
Figure 10. The signals of pre- and post-compression with time shift.
Figure 10. The signals of pre- and post-compression with time shift.
Materials 08 05458 g010
Figure 11. Flow chart of algorithms. (a) Cross correlation (CC); (b) Absolution difference (AD).
Figure 11. Flow chart of algorithms. (a) Cross correlation (CC); (b) Absolution difference (AD).
Materials 08 05458 g011

3.3.1. Cross Correlation Algorithm, CC

Before applying the cross correlation algorithm, the companding technique was used to correct the lateral displacement. Strain images are formed by analyzing pre-compression and post-compression echo fields U and C, respectively. Both the echo fields are Z × Y dimensional matrices of samples from continuous echo waveforms. For example, U = ( U z , U y ); z = 1… Z; y = 1…Y, for Y waveforms each consisting of Z digitized echo values.
The sum of absolute differences (SAD) algorithm is implemented by selecting the corresponding data kernels from U and C of size   L × M . For all (i, j) in a P × Q search region in C where P > L and Q > M, the SAD coefficients ( i , j ) are computed as follows:
i , j = l = 1 L m = 1 M | C l , m U l + i , m + j |
Let i , j = min { i , j } be the minimum SAD coefficient in the search region. The location of i , j indicates the position ( Z C , Y C ) in C that corresponds to ( Z U , Y U ) in U.
Sixteen non-overlapping SAD kernels were equally distributed over U at points ( Z U , Y U ). From the resultant displacement vectors, 16 corresponding points ( Z C , Y C ) were estimated. Using linear regression analysis, the companding parameters m y and m z   and the shift parameters   b y   and   b z   were solved as follows:
Z c = m z Z U + b z ,         Y c = m y Y U + b y
We can use the CC algorithm to estimate the axial shift of the ultrasound echo signal. For continuous functions, the cross-correlation function is defined as follows [1]:
R X Y ( t ) = def X * ( τ ) Y ( t + τ ) d τ
where R X Y ( t ) is the result of cross-correlation analysis, and X ( τ ) and Y ( τ ) are the precompression and post-compression radio frequency signals, respectively. When the functions match, the value of R X Y ( t ) is maximized.
Cross-correlation is easy to implement. The only reason why cross-correlation is not used in real time applications is its computing load. In addition, cross-correlation can produce a “false peak” [22], that occurs when noise, finite window size, and signal decorrelation increase the amplitude of a secondary correlation peak above that of the primary correlation peak. False peaks are relatively large in magnitude and appear as discontinuities in the displacement vector. Thus, false peak correction is needful. Calculating the slope between two points, when slope exceeds a preset threshold, the previous point replaces the first point.

3.3.2. Absolute Difference Algorithm, AD

Dynamic programming technique is used to estimate the displacement [13]. Let us consider two echo signals g(i) and g′(i) corresponding to the two A-lines acquired before and after compression. The difference between the two signals can be quantified using absolute differences as below:
Δ ( i , d ) = | g ( i ) g ( i + d ) |
where d min d d max is the displacement at the sample, and d min and d max   specify the allowed displacement. In ultrasound machines, the gains of radio frequency data can be altered to improve visualization. In order to reduce the effect of these variations on Δ, both pre-compression and post-compression ultrasound images are divided by the maximum value of one of the images. The smoothness of the displacements (S) is given as:
S ( d i , d i 1 ) = ( d i d i 1 ) k
where d i is the displacement at sample i, and d i 1 is the displacement at sample i1 of g(i). The cost function C at the ith point and the associated displacement d i is defined as a recursive function as follows:
C ( i , d i ) = min d i 1 { C ( i 1 , d i 1 ) + w S ( d i , d i 1 ) } + Δ ( i , d i )
where w is a regularization weight which governs smoothness.
The minimum of a scalar function of several variables must then be determined. The value of d i 1 that minimizes is also memorized in a function for later use [19]:
M ( i , d i ) = arg  min d i 1 { C ( i 1 , d i 1 ) + w S ( d i , d i 1 ) }
nonlinear:
{ D ( i ) = arg  min d i { C ( i , d i ) } ,   i = m D ( i ) = M ( i + 1 , D ( i + 1 ) ) ,   i = 1 m 1
The displacements map of all the A-lines are independently calculated using the same procedure.

4. Experimental Results

4.1. Trigger Single Trial

In this thesis, a tissue-mimicking phantom, which is 50 mm long × 40 mm wide × 15 mm high is used. The tissue-mimicking phantom includes gelatin as a scattering factor to simulate the cell tissue of the human body. Figure 12 shows the photograph of the tissue-mimicking phantom. A 0.3-V chirp and short pulse with center frequency of 25 MHz as trigger signal is used to scan the tissue-mimicking phantom.
Figure 12. The photograph of a tissue-mimicking phantom.
Figure 12. The photograph of a tissue-mimicking phantom.
Materials 08 05458 g012
Figure 13 shows the ultrasound images of the tissue-mimicking phantom by the short pulse imaging and chirp-coded excitation imaging. The image of the chirp-coded excitation is brighter than the image of the short pulse. For quantitative comparison, Figure 14 shows the signal-to-noise ratio as a function of depth from the middle A-line of both of the images. The average SNR of the chirp-coded excitation is about 40 dB and the average SNR of the short pulse is about 25 dB. There is a 15 dB SNR improvement by chirp-coded excitation.
Figure 13. The images of tissue-mimicking phantom (a) short pulse and (b) chirp-coded excitation.
Figure 13. The images of tissue-mimicking phantom (a) short pulse and (b) chirp-coded excitation.
Materials 08 05458 g013
Figure 14. The signal-to-noise ratio as a function of depth.
Figure 14. The signal-to-noise ratio as a function of depth.
Materials 08 05458 g014
The average power of the chirp signal is higher so the higher signal energy of the chirp signal is able to penetrate deeper into the tissue with the same tissue attenuation. Then pulse compression can recover the axial resolution and focus the energy of the signal to get stronger signal. Figure 15 shows the images of the tissue-mimicking phantom by unipolar pulse, bipolar pulse and chirp-coded excitation. There is a 1–2 mm penetration depth improvement by the chirp-coded excitation.
Figure 15. The images of the tissue-mimicking (a) unipolar pulse; (b) bipolar pulse; or (c) chirp-coded excitation.
Figure 15. The images of the tissue-mimicking (a) unipolar pulse; (b) bipolar pulse; or (c) chirp-coded excitation.
Materials 08 05458 g015

4.2. Strain-Stress Curve

Human soft tissue can be seen as a continuous material; therefore, it can be formulated as a sequential mechanics model. Hence, the stress-strain relations of human tissue can be modeled using the Hooke’s law. The relationship between stress and strain of a tendon is known to be in the longitudinal direction [25]. Moreover, the curve can be divided into four sections: toe area, linear elastic region, nonlinear elastic region, and inelastic region. In this study, we used scales and a three-axis step motor to record the force and compression distance under a pressure of 0.1 mm for each record reading of scale.
The percentage strain values for each were plotted against the applied stress values. Figure 16 shows the experimental strain-stress loading curves of the phantoms. These curves are very similar to those mentioned above. Thereafter, we calculated the slope of the curve to obtain the Young’s modulus of the phantoms. The Young’s modulus values of the hard and soft phantoms are 26.86 and 13.27 kPa, respectively.
Figure 16. The strain-stress loading curves of two phantoms of different hardness.
Figure 16. The strain-stress loading curves of two phantoms of different hardness.
Materials 08 05458 g016

4.3. Optimal Parameters of Algorithm in Elastography

In order to accurately test the performance of chirp-coded pulse excitation for strain imaging, we employed two phantoms. One uniform phantom without any lesions was used for investigating the effect of pulse length, applied strain, and window length on strain imaging. Another lesion phantom having one lesion was used for testing the CNRe of the chirp pulse. Sections 4.3.1–4.3.3 describe the experimental method, parameters, and the results obtained using the uniform phantom, respectively. Sections 4.3.4–4.3.7 describe the experimental method, parameters, and the results obtained with the lesion phantom. A short pulse is a sine wave with a center frequency (f0) of 7.5 MHz. For a chirp pulse, the default center frequency is 7.5 MHz, and the starting and ending frequencies are 4 and 11 MHz, respectively.

4.3.1. Effect of Pulse Length

A longer chirp pulse must be applied to achieve a larger eSNR gain. However, a longer pulse length does not necessarily ensure better SNRe. Figure 17 shows the optimal chirp pulse length for strain imaging. The longer the chirp pulse, the more distorted the received coded waveform owing to dynamic focusing [20,26]. When the pulse length is greater than 20 cycles, SNRe decreases with pulse length. Moreover, we can see that chirp 1 performs better than chirp 2.
Figure 17. Elastographic signal-to-noise ratio (SNRe) for chirp pulse as a function of pulse length. Algorithm for (a) AD and (b) CC. The correlation window length is 0.77 mm.
Figure 17. Elastographic signal-to-noise ratio (SNRe) for chirp pulse as a function of pulse length. Algorithm for (a) AD and (b) CC. The correlation window length is 0.77 mm.
Materials 08 05458 g017

4.3.2. Effect of Applied Strain

Since the post-deformation echo signal is a function of displacements of tissue scatters, the performance of chirp pulse is also a function of applied strain. In order to determine the optimal applied strain for a chirp pulse, the applied strains were varied from 0.25% to 1.25%. Figure 18 shows that an optimal applied strain exists for chirp-coded strain imaging. The SNRe is the maximum when the strain is 1%, which is in accordance with previous investigation. In elastography, the applied compression is typically up to 1% of the total depth of the tissue [2] because a large applied strain will cause a large decorrelation noise of echo signal. However, for calculating strains, a large applied strain is less sensitive to displacement errors caused by system noise, finite window length, and signal decorrelation noise when compared with a small strain. Chirp 1 provides the greatest improvement over the short pulse. Conversely, the effect of chirp 2 is worse than that of the other two pulses.
Figure 18. SNRe for all pulse as a function of applied strain. Algorithm for (a) AD and (b) CC. The correlation window length is 0.77 mm. The pulse length of chirp is 20 cycle.
Figure 18. SNRe for all pulse as a function of applied strain. Algorithm for (a) AD and (b) CC. The correlation window length is 0.77 mm. The pulse length of chirp is 20 cycle.
Materials 08 05458 g018

4.3.3. Effect of Correlation Window Length

For strain imaging, attainable axial strain resolution is largely determined by correlation window length [27]. A shorter window length will lead to better axial strain resolution than a long window length. In order to compare the performance of chirp and short pulses in terms of window length, we analyzed the effects of different window lengths on strain imaging for both short and chirp pulses. Figure 19 shows that the chirp 1 pulse requires a shorter window length to realize the same SNRe as that obtained with the short pulse. For example, to achieve a SNRe of 15 dB, the window length for short pulse is 1 mm, whereas it is only 0.77 mm for chirp 1 pulse. Therefore, the chirp pulse achieves a better axial strain resolution than the short pulse because the increased SNRe with the chirp pulse permits the use of a shorter window length. We selected a window length of 0.77 mm because the other window lengths are relatively unstable. Moreover, the smaller the window length, the higher is the axial resolution. However, chirp 2 performs poorly when compared with the other pulses. Although a long window length typically results in a higher SNRe, it also lowers the axial strain resolution. Therefore, a tradeoff exists between window length and axial strain resolution.
Figure 19. SNRe as a function of correlation window length.
Figure 19. SNRe as a function of correlation window length.
Materials 08 05458 g019

4.3.4. CNRe

In order to investigate the CNRe of chirp and short pulses, we use the lesion phantom containing one lesion. Figure 20 shows the measured CNRe values for the lesion phantom in Figure 21. The CNRe values of the strain images obtained using AD algorithm are 49.76 ± 6.17 dB and 51.89 ± 4.87 dB for the short and chirp pulses, respectively. The CNRe values of the strain images obtained using the CC algorithm are 53.05 ± 6.03 dB and 57.18 ± 3.5 dB for the short and chirp pulses, respectively. Therefore, the chirp pulse achieves a slightly higher CNRe than the short pulse. Moreover, the short pulse is unstable when compared with the chirp pulse.
Figure 20. A strain image of short pulse and chirp pulse by CC and AD, respectively. The pulse length of chirp is 20 cycle. The correlation window length is 0.77 mm and applied strain is 1% are fixed.
Figure 20. A strain image of short pulse and chirp pulse by CC and AD, respectively. The pulse length of chirp is 20 cycle. The correlation window length is 0.77 mm and applied strain is 1% are fixed.
Materials 08 05458 g020
Figure 21. Elastographic contrast-to-noise ratio (CNRe) of phantom.
Figure 21. Elastographic contrast-to-noise ratio (CNRe) of phantom.
Materials 08 05458 g021

4.3.5. Correlation Window Length

To validate the conclusion that a chirp pulse exhibits better axial strain resolution than a short pulse, we compare the axial strain resolutions of the chirp and short pulses on the lesion phantom. The SNRe values of the two backgrounds and targets in Figure 22 are nearly same; however, the window length for the short pulse in Figure 23a is 1 mm and that for the chirp pulse in Figure 22b is 0.77 mm. In order to compare the target in Figure 22, the target in Figure 22a is blurred more than the target in Figure 22b. Therefore, Figure 22 shows that the chirp pulse has a better axial strain resolution than the short pulse because the increased SNRe with the chirp pulse permits the use of a shorter pulse length.
Figure 23 shows a comparison of the strain images of the lesion phantom with 1% applied strain for short and chirp pulses for correlation window lengths of 0.385, 0.77, 1.16 and 1.54 mm. Using shorter correlation window, the strain image has higher resolution, it has more calculating operation time. The measured CNRe of the Figure 23 results are shown in Table 2. For the same correlation window length, the chirp pulses provide greater CNRe values than short pulses. Therefore, to achieve the same level of lesion detectability, the chirp pulse requires a shorter correlation window.
Table 2. Elastographic contrast-to-noise ratio (CNRe) measured for lesion phantom on Figure 23.
Table 2. Elastographic contrast-to-noise ratio (CNRe) measured for lesion phantom on Figure 23.
Windows Length (mm)Short PulseChirp Pulse
0.3812.1614.71
0.7727.2829.83
1.1645.0647.03
1.5448.4849.32
Figure 22. Strain images for two pulse. (a) The window length is 1 mm for short pulse, and the SNRe of the background and target is 15.99 and 19.69, respectively; and (b) the window length is 0.77 mm for chirp pulse, and the SNRe of the background and target is 16.03 and 20.12, respectively.
Figure 22. Strain images for two pulse. (a) The window length is 1 mm for short pulse, and the SNRe of the background and target is 15.99 and 19.69, respectively; and (b) the window length is 0.77 mm for chirp pulse, and the SNRe of the background and target is 16.03 and 20.12, respectively.
Materials 08 05458 g022
Figure 23. Strain imaging with different correlation window length. The correlation windows are 0.38, 0.77, 1.16 and 1.54 mm from top to bottom.
Figure 23. Strain imaging with different correlation window length. The correlation windows are 0.38, 0.77, 1.16 and 1.54 mm from top to bottom.
Materials 08 05458 g023

4.3.6. Differential Strain

In elastography, the applied strain is theoretically up to 1% of the total depth of the tissue [2]. When the strain is greater than 1%, decorrelation noise will appear in strain images. Figure 24 shows the strain images of the lesion phantom for applied strains of 1%, 2% and 3% using short and chirp pulses. The square boxes in Figure 24 indicate the regions from which the CNRe values shown in Table 3 are computed. As can be clearly seen, the chirp pulse successfully suppresses decorrelation strain noise.
Table 3. CNRe measured for lesion phantom on Figure 24.
Table 3. CNRe measured for lesion phantom on Figure 24.
PhantomCCAD
StrainShortChirpShortChirp
1%52.156.748.550.2
2%40.842.313.216.1
3%0.229.15.25.3
Figure 24. Strain images of a lesion phantom with 1%, 2% and 3% strain applied from top. (a) CC algorithm; and (b) AD algorithm.
Figure 24. Strain images of a lesion phantom with 1%, 2% and 3% strain applied from top. (a) CC algorithm; and (b) AD algorithm.
Materials 08 05458 g024

4.3.7. Lateral Resolution

In order to evaluate the accuracy of the system, we employed a hard cylinder embedded in the phantom at a depth of 20 mm at diameter of 6, 5, 4 and 3 mm as the scanning target. Figure 25a shows the strain images for the chirp pulse obtained using the two algorithms. We selected a few regions including the hard cylinder, whose locations and sizes are shown in Figure 25a, to compute the average strain profile of the elastogram. Figure 25b shows that the mean values of the estimates obtained using the CC algorithm are similar to those obtained using the AD algorithm. Further, the contrast of the strain image obtained using the AD algorithm is lower than that of the strain image obtained using the CC algorithm.
Figure 25. (a) Image of strain; and (b) curve of strain of strain image at 22 mm of depth.
Figure 25. (a) Image of strain; and (b) curve of strain of strain image at 22 mm of depth.
Materials 08 05458 g025

4.4. Young’s Modulus

According to above results, we used chirp pulse length, applied strain, and window length values of 20 cycles, 0.1%, and 0.77 mm, respectively, for all the measurements. The peak eSNR measured for our laboratory system is 30 dB for the lesion phantom. Considering the performances of chirp 1 and chirp 2 pulses, the chirp 1 scheme is applied to the ultrasonic strain imaging system for the phantom experiment.
Young’s modulus indicates the quality of stiffness. The phantoms were fabricated using different concentrations of agar. One phantom contained a cylinder softer than the background, and another phantom contained a cylinder harder than the background. Figure 26 shows the strain images of the two phantoms, and we can clearly see the different elasticities of the phantoms. In Figure 26, the red and yellow square boxes indicate the regions from which the Young’s modulus values shown in Table 4 and Figure 27 are computed, respectively.
Table 4. Young’s modulus of phantom.
Table 4. Young’s modulus of phantom.
AlgorithmPhantomShort Pulse (kPa)Chirp (kPa)
CCSoft14.50 ± 0.1614.97 ± 0.66
Hard25.03 ± 0.7725.52 ± 1.22
ADSoft14.66 ± 0.5314.43 ± 0.08
Hard22.52 ± 0.7022.72 ± 0.32
Figure 26. Strain image of two phantoms by CC and AD algorithm.
Figure 26. Strain image of two phantoms by CC and AD algorithm.
Materials 08 05458 g026
Figure 27. Young’s modulus of phantom including hard substance for two pulses by two algorithms.
Figure 27. Young’s modulus of phantom including hard substance for two pulses by two algorithms.
Materials 08 05458 g027

4.5. Discussion

The chirp pulse increases eSNR (30 dB) owing to its higher energy, and the increased eSNR improves the SNRe (15.7 dB). Therefore, the chirp pulse realizes a significantly higher SNRe than a short pulse (Figure 14 and Figure 15). The results of [17] show that some chirp schemes such as a chirp pulse excitation with 40% tapered Tukey window and matched compression filter performs better than short pulse scheme. Theoretically, the increased eSNR with chirp pulse should improve SNRe; however, the reason for the poor performance of the chirp 2 scheme when compared with that of the short pulse is unknown. Chirp pulse excitation with Gaussian window and matched compression filter is the most commonly used scheme for B mode imaging. However, Figure 14 and Figure 15 show that the common chirp scheme for B mode imaging is not the optimal scheme for elastograms.
These results can be explained in terms of range side lobe level (RSLL) and main lobe width. The main lobe of the chirp scheme blurs pre-compressed and post-compressed echo signals; therefore, the broader the main lobe, the greater is the decorrelation noise caused by the main lobe. Figure 28 shows that the chirp 1 scheme is closer to the main lobe than the short pulse, and the chirp 1 scheme has a broader main lobe when compared with those of the other pulses. However, the chirp 1 scheme exhibits a higher RSLL than the short pulse. However, for the chirp 1 scheme, the normalized amplitudes of the side lobes except those close to the main lobe are below −20 dB. Therefore, the decorrelation noise due to the RSLL of chirp 1 scheme is small when compared with that due to the main lobe and can be neglected. On the other hand, the chirp 1 scheme exhibits greater eSNR than the short pulse and the chirp 2 scheme, as shown in Table 5. However, the main lobe width was not considered for evaluating the optimal chirp scheme for B mode imaging in [24]. Therefore, the performance of chirp pulse primarily depends on the main lobe width and energy. In strain imaging, an appropriate chirp scheme must be selected; otherwise, the chirp scheme may perform worse than the short pulse.
Table 5. The eSNR of three pulses.
Table 5. The eSNR of three pulses.
PulseShort PulseChirp 1 (Tukey)Chirp 2 (Gaussian)
eSNR (dB)20.2233.0930.36
Figure 28. (a) Comparison of envelopes of short pulse, chirp 1 scheme and chirp 2 scheme with 20 cycles of pulse length; and (b) the detail in (a) from 27 to 29 MHz.
Figure 28. (a) Comparison of envelopes of short pulse, chirp 1 scheme and chirp 2 scheme with 20 cycles of pulse length; and (b) the detail in (a) from 27 to 29 MHz.
Materials 08 05458 g028
A long chirp pulse penetrates deeper and achieves a higher eSNR gain than a short chirp pulse. However, the chirp pulse length is limited by dynamic receiving focusing, which distorts the received coded waveform prior to beam sum [26]. Figure 17 shows that an optimal pulse length exists for chirp-coded strain imaging. Therefore, we must select a short chirp pulse (approximately 20 cycles) for the excitation.
In the CC algorithm, correlation window length affects the attainable axial strain resolution. Since the chirp pulse leads to increased SNRe, it can use a shorter window length than short pulse to achieve the same SNRe, as shown in Figure 19. Therefore, the chirp pulse exhibits a better axial strain resolution than the short pulse. Moreover, Figure 18 and Figure 19 show that chirp pulse has better lesion detectability than the short pulse because of the higher CNRe, which is consistent with the results of [20]. Moreover, both the short and chirp pulses must have the same the lateral resolution, which is proportional to the beamwidth of the ultrasound system, because chirp coding and decoding procedures have no effect on the transducer aperture size [20]. Figure 13 and Figure 14 show that the optimal chirp pulse length is 20 cycles and the optimal strain is 1% of the total depth of the tissue for both the CC and AD algorithms. As can be seen in Figure 20, the CC algorithm performs better than the AD algorithm in elastograms because it can achieve a higher CNRe. Figure 24 shows that the optimal strain imaging can be achieved with a chirp pulse for excitation signal using the CC algorithm and over 1% strain. In order to compare the performances of the two algorithms, we plot the strain curve for the lesion phantom. Figure 25b shows that the strain curve corresponding to the AD algorithm is less obvious than that of the CC algorithm for the hard phantom. On the other hand, Figure 27 shows the elasticity distribution along the axial centerline of the phantom in Figure 26. The results of the phantom experiment indicate that the estimated and real values are significantly different, in accordance with the results in [28]. Makoto et al. reported that elasticity distribution can be reconstructed more quantitatively using a 3D model than with a two-dimensional (2D) model because these methods assume a 2D state, namely, a plane strain or plane stress state, which is not always satisfied in actual biological measurements [28]. The conventional method only considers axial and lateral displacement and the lack of positive displacement amount.

5. Conclusions

In this study, we proposed an ultrasound elastography system with chirp-coded excitation for evaluating the elasticity of a phantom (Young’s modulus of the background materials and cylindrical inclusion were 13.27 and 26.86 kPa, respectively) and investigating the strain performance of CC and AD algorithms. In order to minimize side lobes without degrading axial resolution, a mismatched compression filter was used to decode the echo signals.
Development of a chirp coded excitation ultrasound system is developed by using chirp signal as the ultrasound trigger signal in order to achieve better signal-to-noise ratio and deeper penetration depth. Table 6 shows that the echo signal of short pulse is less than the echo signal of chirp-coded excitation. According to experiment results and discussion, there are a 15 dB SNR improvement and a 1–2 mm penetration depth improvement by chirp-coded excitation, pulse compression, the appropriate receiver and the handmade expander of noise reduction and SNR enhancement.
Table 6. The relationship between the pulse type and echo peak amplitude.
Table 6. The relationship between the pulse type and echo peak amplitude.
Pulse TypePeak Frequency (MHz)Echo Peak Amplitude (V)
Unipolar Pulse250.96
Bipolar Pulse251.04
Coded excitation (Chirp)252.08
We compared the performance of different chirp schemes with the conventional short pulse. The chirp with Gaussian window, which is common and suitable for B-mode imaging, was found to exhibit poor performance of SNRe than both the short pulse and the chirp with Tukey window because the width of the main lobe was not been considered for evaluating the optimal chirp scheme for B-mode imaging. Thereafter, the effects of different factors such as chirp pulse length, applied strain, and correlation window length on strain imaging were investigated. Optimal parameters of the three pulses were the same for both the CC and AD algorithms.
The homogenous phantom experimental results show that the SNRe of elastography measured using a short pulse is 11 dB. The SNRe values measured using a 20-cycle chirp-coded ultrasound system modulated using CC and AD algorithms were 15 and 13 dB, respectively. The CNRe of the image obtained using the chirp-coded pulse can be improved by 4.1 dB when compared with that obtained using the short pulse. The results validate that a chirp with Tukey window has better lesion detectability than a short pulse. Additionally, the Young’s modulus values of the cylindrical inclusion analyzed using the CC and AD algorithms were 25.52 and 22.72 kPa, respectively. These results show that the ultrasound elasticity imaging system with chirp-coded excitation modulated by a Tukey window can acquire highly accurate and high quality elastography images. In the future, the ultrasound elasticity imaging system will be used on human subjects to validate the feasibility of collecting in vivo data [29,30,31].

Acknowledgments

This study is supported in part by the Ministry of Science and Technology, Taiwan, under contract number: MOST 102-2221-E-006-012-MY3, duration: 1 August 2013–31 July 2016.

Author Contributions

G.-C.C., H.-J.C. and K.-H.L. conceived and designed the experiments; H.-J.C. performed the experiments; H.-J.C. and C.-M. L. analyzed the data; P.-J. C. and T. C. contributed reagents/materials/analysis tools; G.-C.C. and T. C. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ophir, J.; Alam, S.K.; Garra, B.S.; Kallel, F.; Konofagou, E.E.; Krouskop, T.; Merritt, C.R.B.; Righetti, R.; Souchon, R.; Srinivasan, S. Elastography: Imaging the elastic properties of soft tissues with ultrasound. J. Med. Ultrason. 2002, 29, 155–171. [Google Scholar] [CrossRef]
  2. 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, 111–134. [Google Scholar] [CrossRef] [PubMed]
  3. Zhai, L.; Madden, J.; Foo, W.C.; Palmeri, M.L.; Mouraviev, V.; Polascik, T.J.; Nightingale, K.R. Acoustic radiation force impulse imaging of human prostates ex vivo. Ultrasound Med. Biol. 2010, 36, 576–588. [Google Scholar] [CrossRef] [PubMed]
  4. Gallotti, A.; D’onofrio, M.; Romanini, L.; Cantisani, V.; Mucelli, R.P. Acoustic Radiation Force Impulse (ARFI) ultrasound imaging of solid focal liver lesions. Eur. J. Radiol. 2012, 81, 451–455. [Google Scholar] [CrossRef] [PubMed]
  5. Balleyguier, C.; Canale, S.; Hassen, W.B.; Vielh, P.; Bayou, E.H.; Mathieu, M.C.; Uzan, C.; Bourgier, C.; Dromain, C. Breast elasticity: Principles, technique, results: An update and overview of commercially available software. Eur. J. Radiol. 2013, 82, 427–434. [Google Scholar] [CrossRef] [PubMed]
  6. Destounis, S.; Gruttadauria, J.L. Elasticity imaging 101. J. Radiol. Nurs. 2013, 32, 124–130. [Google Scholar] [CrossRef]
  7. Gennisson, J.L.; Deffieux, T.; Fink, M.; Tanter, M. Ultrasound elastography: Principles and techniques. Diagn. Interv. Imaging 2013, 94, 487–495. [Google Scholar] [CrossRef] [PubMed]
  8. Varghese, T. Quasi-static ultrasound elastography. Ultrasound Clin. 2009, 4, 323–338. [Google Scholar] [CrossRef] [PubMed]
  9. Gao, L.; Parker, K.J.; Lerner, R.M.; Levinson, S.F. Imaging of the elastic properties of tissue—A review. Ultrasound Med. Biol. 1996, 22, 959–977. [Google Scholar] [CrossRef]
  10. Benson, J.; Fan, L. Tissue Strain Analytics—A Complete Ultrasound Solution for Elastography. 2012. Available online: http://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0ahUKEwjRqMarqKvJAhUKJ5QKHZ8CDd0QFggjMAA&url=http%3A%2F%2Fsonoworld.com%2FCommon%2FDownloadFile.aspx%3FModuleDocumentsId%3D66&usg=AFQjCNFvP5Ymuj2jd1BYncxUJoXKBSF9Xg&bvm=bv.108194040,d.dGo&cad=rja (accessed on 25 November 2015).
  11. Nightingale, K.R.; Palmeri, M.L.; Nightingale, R.W.; Trahey, G.E. On the feasibility of remote palpation using acoustic radiation force. J. Acoust. Soc. Am. 2001, 110, 625–634. [Google Scholar] [CrossRef] [PubMed]
  12. Sarvazyan, A.P.; Rudenko, O.V.; Swanson, S.D.; Fowlkes, J.B.; Emelianov, S.Y. Shear wave elasticity imaging: A new ultrasonic technology of medical diagnostics. Ultrasound Med. Biol. 1998, 24, 1419–1435. [Google Scholar] [CrossRef]
  13. Nightingale, K.; McAleavey, S.; Trahey, G. Shear-wave generation using acoustic radiation force: In vivo and ex vivo results. Ultrasound Med. Biol. 2003, 29, 1715–1723. [Google Scholar] [CrossRef] [PubMed]
  14. Peng, H.; Liu, D.C. Chirp-coded pulse excitation for ultrasound elasticity imaging. In Proceedings of the 2010 4th International Conference on Bioinformatics and Biomedical Engineering (iCBBE), Chengdu, China, 18–20 June 2010; pp. 1–4.
  15. Liu, J.; Insana, M.F. Coded pulse excitation for ultrasonic strain imaging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 231–240. [Google Scholar] [PubMed]
  16. Chaturvedi, P.; Insana, M.F.; Hall, T.J. 2-D companding for noise reduction in strain imaging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1998, 45, 179–191. [Google Scholar] [CrossRef] [PubMed]
  17. Cespedes, I.; Ophir, J.; Ponnekanti, H.; Maklad, N. Elastography: Elasticity imaging using ultrasound with application to muscle and breast in vivo. Ultrason. Imaging 1993, 15, 73–88. [Google Scholar] [CrossRef] [PubMed]
  18. Fung, Y.C. Biomechanics: Mechanical Properties of Living Tissues; Springer: Berlin/Heidelberg, Germany, 1993. [Google Scholar]
  19. Qiu, W.; Yu, Y.; Tsang, F.K.; Zheng, H.; Sun, L. A novel modulated excitation imaging system for microultrasound. IEEE Trans. Biomed. Eng. 2013, 60, 1884–1890. [Google Scholar] [PubMed]
  20. Peng, H.; Liu, D.C. Enhanced ultrasound strain imaging using chirp-coded pulse excitation. Biomed. Signal Process. Control 2013, 8, 130–141. [Google Scholar] [CrossRef]
  21. Turin, G. An introduction to matched filters. IRE Trans. Inf. Theory 1960, 6, 311–329. [Google Scholar] [CrossRef]
  22. Zahiri-Azar, R.; Salcudean, S.E. Motion estimation in ultrasound images using time domain cross correlation with prior estimates. IEEE Trans. Biomed. Eng. 2006, 53, 1990–2000. [Google Scholar] [CrossRef] [PubMed]
  23. Chaturvedi, P.; Insana, M.F.; Hall, T.J. Testing the limitations of 2-D companding for strain imaging using phantoms. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1998, 45, 1022–1031. [Google Scholar] [CrossRef] [PubMed]
  24. Varghese, T.; Ophir, J. Enhancement of echo-signal correlation in elastography using temporal stretching. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1997, 44, 173–180. [Google Scholar] [CrossRef] [PubMed]
  25. Price, J.M.; Patitucci, P.; Fung, Y.C. Biomechanics. Mechanical Properties of Living Tissues; Springer: Berlin/Heidelberg, Germany, 1981. [Google Scholar]
  26. O’Donnell, M. Coded excitation system for improving the penetration of real-time phased-array imaging systems. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1992, 39, 341–351. [Google Scholar] [CrossRef] [PubMed]
  27. Liu, J.; Abbey, C.K.; Insana, M.F. Linear approach to axial resolution in elasticity imaging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2004, 51, 716–725. [Google Scholar] [PubMed]
  28. Yamakawa, M.; Shiina, T. Tissue elasticity reconstruction based on 3-dimensional finite-element model. Jpn. J. Appl. Phys. 1999, 38, 3393. [Google Scholar] [CrossRef]
  29. Hu, C.H.; Liu, R.; Zhou, Q.; Yen, J.; Shung, K.K. Coded excitation using biphase-coded pulse with mismatched filters for high-frequency ultrasound imaging. Ultrasonics 2006, 44, 330–336. [Google Scholar] [CrossRef] [PubMed]
  30. Santos, S.D.; Domenjoud, M.; Prevorovsky, Z. Ultrasonic imaging of human tooth using chirp-coded nonlinear time reversal acoustics. Phys. Procedia 2010, 3, 913–918. [Google Scholar] [CrossRef]
  31. Pedersen, M.H.; Misaridis, T.X.; Jensen, J.A. Clinical evaluation of chirp-coded excitation in medical ultrasound. Ultrasound Med. Biol. 2003, 29, 895–905. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Chun, G.-C.; Chiang, H.-J.; Lin, K.-H.; Li, C.-M.; Chen, P.-J.; Chen, T. Ultrasound Elasticity Imaging System with Chirp-Coded Excitation for Assessing Biomechanical Properties of Elasticity Phantom. Materials 2015, 8, 8392-8413. https://doi.org/10.3390/ma8125458

AMA Style

Chun G-C, Chiang H-J, Lin K-H, Li C-M, Chen P-J, Chen T. Ultrasound Elasticity Imaging System with Chirp-Coded Excitation for Assessing Biomechanical Properties of Elasticity Phantom. Materials. 2015; 8(12):8392-8413. https://doi.org/10.3390/ma8125458

Chicago/Turabian Style

Chun, Guan-Chun, Hsing-Jung Chiang, Kuan-Hung Lin, Chien-Ming Li, Pei-Jarn Chen, and Tainsong Chen. 2015. "Ultrasound Elasticity Imaging System with Chirp-Coded Excitation for Assessing Biomechanical Properties of Elasticity Phantom" Materials 8, no. 12: 8392-8413. https://doi.org/10.3390/ma8125458

APA Style

Chun, G. -C., Chiang, H. -J., Lin, K. -H., Li, C. -M., Chen, P. -J., & Chen, T. (2015). Ultrasound Elasticity Imaging System with Chirp-Coded Excitation for Assessing Biomechanical Properties of Elasticity Phantom. Materials, 8(12), 8392-8413. https://doi.org/10.3390/ma8125458

Article Metrics

Back to TopTop