AUTOREGRESSIVE SIGNAL PROCESSING APPLIED TO HIGH-FREQUENCY ACOUSTIC MICROSCOPY OF SOFT TISSUES

Abstract

A method to create a parameter map depicting acoustical and mechanical properties of biological tissue at microscopic resolutions to identify potential health related issues. The method including mounting the biological tissue on a substrate, raster scanning the biological tissue with an RF frequency, recovering RF echo signals from said substrate and from a plurality of locations on said biological tissue, wherein each of the plurality of locations corresponds to a specific pixel comprising the parameter map, the recovered RF echo signals including a reference signal recovered from the substrate at a point devoid of tissue, a first sample signal recovered from an interface between the biological tissue and water, and a second sample signal recovered from an interface between said biological tissue and said substrate, repeatedly applying a plurality of computer-generated calculation steps based on the reference signal, the first sample signal and the second sample signal to generate estimated values for a plurality of parameters associated with each of the specific pixels in the parameter map. The plurality of computer-generated calculation steps includes a denoising step, and using the generated estimated values to create said parameter map depicting parameters including, but not limited, to acoustic impedance, speed of sound, ultrasound attenuation, mass density, bulk modulus and nonlinear attenuation.

Claims

1. A method to create a parameter map depicting acoustical and mechanical properties of biological tissue at microscopic resolutions to identify potential health related issues, comprising, mounting said biological tissue on a substrate, raster scanning the biological tissue with an RF frequency, recovering RF echo signals from said substrate and from a plurality of locations on said biological tissue, wherein each of the plurality of locations corresponds to a specific pixel comprising the parameter map, the recovered RF echo signals including: a reference signal recovered from the substrate at a point devoid of tissue, a first sample signal recovered from an interface between the biological tissue and water, and a second sample signal recovered from an interface between said biological tissue and said substrate, repeatedly applying a plurality of computer-generated calculation steps based on the reference signal, the first sample signal and the second sample signal to generate estimated values for a plurality of parameters associated with each of the specific pixels in the parameter map; wherein the plurality of computer-generated calculation steps including a denoising step, and using the generated estimated values to create said parameter map depicting parameters including, but not limited, to acoustic impedance, speed of sound, ultrasound attenuation, mass density, bulk modulus and nonlinear attenuation.

2. The method as recited in claim 1 wherein the biological tissue has a thickness less than 6 m.

3. The method as recited in claim 1 the RF echo signals are two or more signals.

4. The method as recited in claim 1 wherein the acoustic impedance is less than 1.56 MRayl.

5. The method as recited in claim 1 further comprising a perturbation signal, wherein the perturbation signal is equal to or greater than 0.2.

6. The method as recited in claim 1 wherein an autoregressive model is implemented.

7. The method as recited in claim 1 further comprising estimating a non linear attenuation based on a power law model.

8. The method as recited in claim 3, wherein the signals overlap.

9. The method as recited in claim 6, wherein the autoregressive model uses an entire normalized spectra.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0010] The patent or application filed contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.

[0011] Further objects, features and advantages of the invention will become apparent from the following detailed description taken in conjunction with the accompanying figures showing illustrative embodiments of the invention, in which:

[0012] FIG. 1 is a diagram of the QAM approach;

[0013] FIG. 2 shows the results of parameter variation and resulting error in estimating acoustic-properties values;

[0014] FIGS. 3a-3e provides simulation results showing estimation error and standard deviation of exponent for varying components using the NLAR model;

[0015] FIGS. 4a-4d provides simulation results showing estimation error and standard deviation for varying components for varying amplitude of the perturbation signal;

[0016] FIGS. 5a-5j shows simulated and measured signals with fit signals using the Hozumi model and the AR model;

[0017] FIGS. 6a-6h show a QAM amplitude map and parameter maps generated using the Hozumi model and the AR model;

[0018] FIGS. 7a-7f show parameter maps from RO12 of FIG. 6a generated using the Hozumi model and the AR model;

[0019] FIGS. 8a-8d show parameter maps of a human lymph node derived using the NLAR model;

[0020] FIGS. 9a-9h show an amplitude and histology map of a pig comes and parameter maps generated using the Hozumi model and the AR model; and

[0021] FIGS. 10a-10c show normalized and varied spectra of a signal.

DETAILED DESCRIPTION

[0022] 2. Theory

[0023] A. Forward Model

[0024] FIG. 1 depicts the experimental approach used to collect data using transducers 40 and tissue sample 15 on glass plate substrate 10. The sample is raster scanned in 2D, and RF echo signals are acquired at each scan location. The RF echo signal acquired from a location devoid of tissue is referred to as the reference signal, 20, and is symbolized by s.sub.0(tt.sub.0). This notation indicates that the reference signal is composed of only one echo at the glass-water interface. Other scanned locations are referred as sample signals 25. 30, and symbolized with s.sub.1(t) being the reflection from the sample-water interface and s.sub.2(t) being the reflection from the sample-substrate interface. d denotes sample thickness of tissue sample 15. Table 1 below lists all symbols used to describe the theory.

[0025] Similarly, we refer to signals derived at all other scanned locations as sample signals, symbolized by s(t). Our forward model is described using the following expression:

[00001]$\begin{array}{cc}s\left(t\right)=s_1\left(t\right)+s_2\left(t\right)+\mathrm{.Math.}+s_n\left(t\right)& \left(1\right)\\ =C_1.Math.s_0^{*}\left(t-t_1\right)+C_2.Math.s_0^{*}\left(t-t_2\right)+\mathrm{.Math.}+C_n.Math.s_0^{*}\left(t-t_n\right),& \left(2\right)\end{array}$

[0026] which means that the sample signals are the sum of n weighted and delayed versions of the reference signal and the * symbol represents frequency-dependent attenuation effects. Among these n signals, two (i.e., s.sub.k.sub.1, and s.sub.k.sub.2 with k.sub.1k.sub.2 and t.sub.k.sub.1<t.sub.k.sub.2) are the echoes from the water-tissue and tissue-glass interfaces, respectively. These two signals are the ones needed to obtain quantitative tissue properties. The remaining n2 signals (i.e., s.sub.k.sub.3, . . . , s.sub.k.sub.n) account for potential multiple reflections, scattering, or noise.

[0027] By taking the Fourier transform of Eq. (2), we obtain:

S(f)=S.sub.0(f)[C.sub.1 exp(f(.sub.1+j2(t.sub.1t.sub.0)))+ . . . +C.sub.n exp(f(.sub.n+j2(t.sub.nt.sub.0)))],(3)

[0028] where .sub.k is the attenuation coefficient of the k.sup.th signal in Np/Hz and S.sub.0(f) is the Fourier transform of the reference signal (i.e., S.sub.0(f)=FFT(s.sub.0)). By dividing S(f) by s.sub.0(f) eq. (3) yields the normalized spectrum:

N(f)=S(f)/S.sub.0(f)

=.sub.k=1.sup.nC.sub.k(f[.sub.k+j2(t.sub.kt.sub.0)])(4)

TABLE-US-00001 TABLE 1 Used symbols Symbol Explanation A.sub.M amplitude in mV and determined by the maximum of the envelope signal =.sub.2/2d .sub.k attenuation coefficient of the k.sup.th signal in Np/Hz C.sub.k signal amplitude c, c.sub.w speed of sound in tissue and coupling medium in m/s d tissue thickness non linear attenuation exponent in (f) = f.sup. .sub.i error term F discrete frequency index f, f.sub.m frequency vector, center frequency in MHz f =fs/(2M) where M is total duration of the sampled signal in points and fs is the sampling frequency f.sub.min, f.sub.m frequencies of the extrema of |N|.sup.2 column vector of unknown C.sub.q H.sub.q, p Hankel matrix = N.sub.q+pi j complex number (i.e., j = {square root over (1)}) k signal index Vandermonde-like matrix of size n by F.sub.2 F.sub.1 + 1 whose element q, p is .sub.q.sup.p, .sub.k =exp(f[.sub.k + j2(t.sub.k t.sub.0)]) N(f) normalized spectrum n number of signals .sub.u(N) unwrapped phase of N p, q used as integer index S.sub.k(f) Fourier Transform of s.sub.k(t) s.sub.k(t), s signal and reference signal in the time domain attenuation coefficient ((f) = f.sup.) in dB/cm TOF time of flight in s and determined by the maximum of the envelope signal t time vector in s x.sub.k AR coefficients Z, Z.sub.w, Z acoustic impedance in MRayl of tissue, coupling, and glass plate, respectively indicates data missing or illegible when filed

[0029] B. Hozumi Inverse Model

[0030] The Hozumi model is compared here to the present invention AR model. Briefly, the Hozumi inverse model assumes n=2 in the forward model and implementation starts by assessing the extrema of the squared magnitude (|N|.sup.2) of the normalized spectra. From eq. (4), using n=2 and assuming .sub.1=0 (i.e., the first signal suffers no attenuation), following expression can be derived:

[00002]$\begin{array}{cc}{.Math.N.Math.}^2={.Math.\frac{\mathrm{FFT}.Math.\left\{s_1\left(t\right)+s_2\left(t\right)\right\}}{\mathrm{FFT}.Math.\left\{s_0\left(t\right)\right\}}.Math.}^2& \left(5\right)\\ =C_1^2+2.Math.C_1.Math.C_2.Math.e^{f.Math..Math.{}_2}.Math.\cos \left(2.Math..Math..Math.f\left(t_2-t_1\right)\right)+C_2^2.Math.e^{2.Math.f.Math..Math.{}_2}.& \left(6\right)\end{array}$

[0031] Then,

[00003]$\begin{array}{cc}{}_u\left(N\right)=2.Math.\left(f_{\mathrm{ma}.Math..Math.x}.Math.\frac{2.Math.d}{c_w}-p\right),p& \left(7\right)\\ \mathrm{and}& \\ {}_u\left(N\right)=4.Math..Math..Math.f_{m.Math..Math.i.Math..Math.n}.Math.\frac{d}{c_w}-\left(2.Math.p-1\right),p& \left(8\right)\end{array}$

[0032] where d is the tissue thickness, f.sub.min and f.sub.max are the frequencies of the extrema of |N|.sup.2, and .sub.u(N) is the unwrapped phase of N. Eq. (7), and (8) are used to estimate the thickness of the specimen. The speed of sound (c) can be estimated using

[00004]$\begin{array}{cc}d=\frac{{}_u\left(N\right)}{4.Math..Math..Math.f_{\mathrm{ma}.Math..Math.x}\left(\frac{1}{c_w}-\frac{1}{c}\right)}& \left(9\right)\end{array}$

[0033] where c.sub.w is the speed of sound in the coupling medium. C.sub.1, C.sub.2 and =.sub.2/2d can be found from the amplitudes of |N|.sup.2. The acoustic impedance of the sample (Z) was estimated from C.sub.1 using first principles. To estimate we used a dichotomy method applied to |N|.sup.2.

[0034] C. Autoregressive Inverse Model (AR Model)

[0035] Initially, the AR model assumes that the signals are composed of more than one signal, i.e., n is assumed to be 2, because that value provides robustness and stability. Most of the n decomposed signals are related to noise and estimation artifacts. The aim of this chapter is to find the two signals that are related to the water-sample and sample-substrate interfaces.

[0036] The inverse model consists of rewriting Eq. (4) at discreet frequencies denoted by f.sub.i=if. The step size (f) is related to the total duration of the sampled signal in points (M) and the sampling frequency (fs) and is defined by fs/(2M). (Note that zero-padding the signal will make f smaller as it is typically done, but in this application it does not provide new information and is therefore unnecessary and was not done.) Discretization yields the following equation for N.sub.i=N(if):

[00005]$\begin{array}{cc}{N}_i={.Math.}_{k=1}^n.Math.{C_k\left[\exp \left(.Math..Math.f\left({}_k+j.Math..Math.2.Math..Math.\left(t_k-t_0\right)\right]\right)\right]}^i& \left(10\right)\\ ={.Math.}_{k=1}^n.Math.C_k.Math.{}_k^i,& \left(11\right)\end{array}$

[0037] where .sub.k=exp(f[.sub.k+j2(t.sub.kt.sub.0)]).

[0038] Then, the AR process is introduced by assuming that N.sub.i can be estimated using a linear combination of the f.sub.i-previous frequencies:

N.sub.i=.sub.k=1.sup.{circumflex over (n)}x.sub.kN.sub.ik+.sub.i,(12)

[0039] where .sub.i is an error term and the x.sub.k are AR coefficients and we choose {circumflex over (n)}=n.

[0040] Based on Eq. (12), the AR inverse model consists of the following four steps: 1) estimating the C.sub.k and .sub.k, 2) Cadzow denoising, 3) determining k.sub.1 and k.sub.2, and 4) estimating all acoustic parameters from C.sub.k.sub.1, C.sub.k.sub.2, .sub.k.sub.1, and .sub.k.sub.2.

[0041] 1) Estimation of C.sub.k and .sub.k

[0042] Equation (12) is rewritten in matrix notation and only is used between frequencies f.sub.1=F.sub.1f and f.sub.2=F.sub.2f, which are determined by the 20-dB bandwidth of the transducer:

N=RX+,(13)

[0043] where N is the column vector of length F.sub.2F.sub.1n+1 and is composed of the values of N.sub.i from F.sub.1+n to F.sub.2, R is the n by F.sub.2F.sub.1n+1 matrix whose element (q,p) is N.sub.F.sub.1.sub.+nq+p1, and is the column vector of length F.sub.2F.sub.1n+1 composed of the values Eq from q=F.sub.1+n to q=F.sub.2.

[0044] Equation (13) is solved for X using a least-squares minimization the sum of the magnitude-squared errors (i.e., .sup.t, where the superscript t symbolizes the transposition operation):

X=(R.sup.tR).sup.1(RN).(14)

[0045] Equations (10) and (12) are combined to establish a relationship between X and the .sub.k:

N.sub.i=.sub.l=1.sup.nx.sub.i.sub.k=1.sup.nC.sub.k.sub.k.sup.il(15)

.sub.k=1.sup.nx.sub.kN.sub.ik=.sub.l=1.sup.nx.sub.i.sub.k=1.sup.nC.sub.k.sub.k.sup.il,(16)

[0046] which yields:

[00006]$\begin{array}{cc}{.Math.}_{k=1}^n.Math.C_k.Math.{}_k^i.Math.P\left(\frac{1}{{}_k}\right)=0,& \left(17\right)\end{array}$

[0047] where P is a polynomial of degree n defined by:

P(z)=1+.sub.l=1.sup.nx.sub.lz.sup.l.(18)

[0048] Equation (17) must be true at all frequencies, thus,

[00007]$P.Math.\left(\frac{1}{{}_k}\right)=0$

for all k. Therefore, the .sub.k are the reciprocal of the n zeros of P. P(z)=0 can be easily solved because the coefficients x.sub.l are known through eq. (14).

[0049] To find the C.sub.k, eq. (10) is written for all i between f.sub.1 a f.sub.2, yielding the following matrix equation:

F=N.sub.2,(19)

[0050] where is a Vandermonde-like matrix of size n by F.sub.2F.sub.1+1 whose element q, p is .sub.q.sup.p, is a column vector of unknown C.sub.q, and N.sub.2 is a column vector whose q.sup.th element is N.sub.F.sub.1.sub.+q1. This equation is solved using a least-squares approach to yield:

=(.sup.t).sup.1(N.sub.2),(20)

[0051] 2) Cadzow Denoising

[0052] To improve the signal decomposition by further taking advantage of the low expected rank of R, Cadzow de-noising was applied to preprocess N and to yield a new vector N.sup.D. The method makes iterative use of singular-value decomposition (SVD) and anti-diagonal averaging to reduce the rank of the Hankel matrix (H.sub.q,p=N.sub.q+p1) as described briefly in the following.

[0053] The first step consists of calculating the SVD such that

H=U*S*V.(21)

[0054] where S is a diagonal matrix with nonnegative diagonal elements (i.e., singular values) in decreasing order and U and V are unitary matrices. In the next step (.sub.i) is reconstructed by keeping only the n largest singular values of S

.sub.i=U*S.sub.n*V,(22)

[0055] where S.sub.n(q,p)=S(q,p) for q, p<n and S.sub.n=0 otherwise. A de-noised version (N.sup.D) of N can be reconstructed by taking the average of all anti-diagonals of .sub.i

N.sub.i.sup.D=mean.sub.q+p=i+1().(23)

[0056] Eqs. (21) to (23) are repeated iteratively (i.e., five times in QAM applications). Although we choose n>2, the aim is to find the two signals, s.sub.k.sub.1, and s.sub.k.sub.2 (k.sub.1k.sub.2 and t.sub.k.sub.1<t.sub.k.sub.2, see eq. (2)), that represent the water-sample and sample-substrate reflections, respectively. Other signal components (i.e., s.sub.k.sub.3 to s.sub.k.sub.n) either are related to noise, estimation artefacts or are actual signals originating from structures inside the tissue (e.g., caused by scattering). To simplify the parameter estimation further, we eliminated signals related to noise and estimation artifacts. To do this, we tested the following two strategies: we set n to a fixed predefined constant (e.g., n=2) or we estimated it based on a threshold procedure. We estimated n to be the K singular values that contribute more than 10% to the overall signal. During the first iteration of the Cadzow de-noising, we calculated

[00008]$\begin{array}{cc}n=1+.Math.\left\{q.Math..Math.s.t..Math.q1.Math..Math.\mathrm{and}.Math..Math.\frac{{S\left(q,q\right)}^2}{{.Math.}_{p1}.Math.{\left(S\left(p,p\right)\right)}^2}>0.1\right\}.Math.,& \left(24\right)\end{array}$

[0057] where represents the cardinal symbol that expresses the number of elements in the enclosed set ({ }). This approach can produce signals with higher orders (i.e., n>2). To find s.sub.k.sub.1 and s.sub.k.sub.2 to estimate the acoustic parameters, we applied Cadzow denoising followed by replacing N by N.sup.D in eq. (10).

[0058] 3) Determining k.sub.1 and k.sub.2

[0059] The next step of the inverse model consists of finding the indices (k.sub.1 and k.sub.2) corresponding to the water-tissue reflection defined by s.sub.k.sub.1 and the tissue-glass reflection defined by s.sub.k.sub.2. The AR inverse model can find signals that are related to noise or artifacts and are not related to the original signal. Some of these unreliable signals easily can be identified and excluded by thresholding the signal amplitude or phase shift (i.e., t.sub.kt.sub.0). For QAM applications, we implemented a simple, but fast and efficient, algorithm to find s.sub.k.sub.1 and s.sub.k.sub.2. The first step of the algorithm consists of calculating the amplitude (A.sub.M in milliVolts) and time of flight (TOF.sub.M in microseconds) of the maximum of the Hilbert-transformed signals. Assuming that the amplitude (C.sub.p=1/{circumflex over (R)}.sub.p where {circumflex over (R)}.sub.p is the reflection coefficient of signal s.sub.p with p{1, . . . , n} and is estimated using eq. (40)) and time of flight (TOF.sub.p=t.sub.pt.sub.0, see eq. (26)) of at least one of the signals s.sub.k.sub.1 or s.sub.k.sub.2 must be close to A.sub.M and TOF.sub.M of the measured signal, the algorithm first finds the s.sub.p, with

p=arg min.sub.p({square root over ((A.sub.MC.sub.p).sup.2+(TOF.sub.MTOF.sub.p).sup.2)}).(25)

[0060] The second signal is then selected by finding the C.sub.q with q=arg min.sub.qp|A.sub.MC.sub.q|. The final step consists of sorting s.sub.p and s.sub.q so that s.sub.p=s.sub.k.sub.1; s.sub.q=s.sub.k.sub.2 if TOF.sub.pTOF.sub.q and s.sub.q=s.sub.k.sub.1; s.sub.p=s.sub.k.sub.2 otherwise. This approach assumes that the two amplitudes from s.sub.k.sub.1 and s.sub.k.sub.2 are the two largest among the n signals of the forward model. If the amplitude from a third signal, originating from scattering or reflections from within the sample, is stronger, then C.sub.k.sub.1 or C.sub.k.sub.2 in the above algorithm will misidentify the desired signals. However, we never encountered this situation in any of our experiments in soft tissue with a glass plate as a substrate.

[0061] 4) Acoustic-Parameter Estimation

[0062] The final step of the AR inverse model consists of estimating, acoustic impedance, speed of sound and attenuation from C.sub.k.sub.1, C.sub.k.sub.2, .sub.k.sub.1, and .sub.k.sub.2.

[0063] The definition of .sub.k (see eq. (10)) directly yields:

[00009]$\begin{array}{cc}{t}_{k_1}-t_0=\frac{\mathrm{imag}\left(\log \left({}_{k_1}\right)\right)}{.Math..Math.f}& \left(26\right)\\ t_{k_2}-t_0=\frac{\mathrm{imag}\left(\log \left({}_{k_2}\right)\right)}{.Math..Math.f},& \left(27\right)\\ \mathrm{and}& \\ {}_{k_1}=\frac{\mathrm{real}\left(\log \left({}_{k_1}\right)\right)}{.Math..Math.f}& \left(28\right)\\ {}_{k_2}=\frac{\mathrm{real}\left(\log \left({}_{k_2}\right)\right)}{.Math..Math.f},& \left(29\right)\end{array}$

[0064] where imag and real symbolize functions that take the imaginary and real part of the argument.

[0065] First principles also lead to the following expressions:

[00010]$\begin{array}{cc}{t}_{k_1}-t_0=\frac{2.Math.d}{c_w}& \left(30\right)\\ t_{k_2}-t_0=\frac{2.Math.d}{c}-\frac{2.Math.d}{c_w},& \left(31\right)\end{array}$

[0066] which can be simultaneously solved to yield:

[00011]$\begin{array}{cc}d=\frac{c_w}{2}.Math.\frac{\mathrm{imag}\left(\log \left({}_{k_1}\right)\right)}{.Math..Math.f}& \left(32\right)\\ c=c_w.Math.\frac{\mathrm{imag}\left(\log \left({}_{k_1}\right)\right)}{\mathrm{imag}\left(\log \left({}_{k_1}\right)\right)+\mathrm{imag}\left(\log \left({}_{k_2}\right)\right)},& \left(33\right)\end{array}$

[0067] To estimate Z, we exploit the fact that C.sub.k.sub.1 is the ratio between the pressure reflection coefficients at the water-tissue and at the water-glass interfaces (i.e., R.sub.wt and R.sub.wg, respectively).

[00012]$\begin{array}{cc}{C}_{k_1}=\frac{R_{\mathrm{wt}}}{R_{\mathrm{wg}}}=\frac{Z-Z_w}{Z+Z_w}.Math.\frac{Z_g+Z_w}{Z_g-Z_w}& \left(34\right)\end{array}$

[0068] which yields

[00013]$\begin{array}{cc}Z=Z_w.Math.\frac{1+\frac{C_{k_1}}{R_{\mathrm{wg}}}}{1-\frac{C_{k_1}}{R_{\mathrm{wg}}}}.& \left(35\right)\end{array}$

[0069] In eq. (34) and eq. (35), Z.sub.w and Z.sub.g stands for the known acoustic impedance of water and glass, respectively.

[0070] To estimate the attenuation, we note that .sub.k.sub.2 is the attenuation coefficient suffered by the signal through round-trip travel into a tissue that has a thickness of d, therefore,

[00014]$\begin{array}{cc}=\frac{{}_{k_2}}{2.Math.d}=\frac{\mathrm{real}\left(\log \left({}_{k_2}\right)\right)}{2.Math.d.Math..Math..Math..Math.f},& \left(36\right)\end{array}$

[0071] which is then converted to dB/cm/MHz.

[0072] First principles also provide an estimate of Z based on C.sub.k.sub.2 using

[00015]$\begin{array}{cc}{C}_{k_2}=\frac{2.Math.Z}{Z_w+Z}.Math.\frac{Z_g-Z}{Z_g+Z}.Math.\frac{2.Math.Z_w}{Z_w+Z}.Math.\frac{Z_w+Z_g}{Z_w-Z_g},& \left(37\right)\end{array}$

[0073] which results because the tissue-glass echo propagates from water to tissue, reflects at the glass-tissue interface, and propagates back from tissue to water. The last term in Eq. (37) also appears in Eq. (34) from the division by S.sub.0 and is the reciprocal of the water-glass pressure reflection coefficient.

[0074] Equation (37) can be rewritten as

Z.sup.3+[(2Z.sub.w+Z.sub.g)+1]Z2+[(Z.sup.2.sub.w+2Z.sub.wZ.sub.g)Z.sub.g]Z+Z.sub.w.sup.2Z.sub.g=0,(38)

[0075] where

[00016]$\begin{array}{cc}={\mathrm{Ck}}_2.Math.\frac{Z_g-Z_w}{4.Math.Z_w\left(Z_g+Z_w\right)}.& \left(39\right)\end{array}$

[0076] Equation (38) is a third-degree polynomial which is solved using closed-form equations to provide three roots. Typically, finding the correct root is straightforward because the correct root is close to the acoustic impedance value of water. One of the other two roots is usually smaller than one MRayl and the last one greater than three MRayl.

[0077] Initial tests with measured signals showed in some cases that the attenuation .sub.k.sub.1 for the first signal (s.sub.k.sub.1), which originates only from ultrasound paths in water (see FIG. 1), is not negligible (i.e., .sub.k.sub.10). Hence, assuming that attenuation in water is negligible at frequencies of 250-MHz and higher can lead to an inaccurate estimate of the acoustic impedance (Z). To obtain an attenuation-corrected estimate of the first signal amplitude (i.e., .sub.k.sub.1), we used

[00017]$\begin{array}{cc}{\hat{C}}_{k_1}=C_{k_1}.Math.{.Math.{}_{k_1}.Math.}^{\frac{\mathrm{fm}}{.Math..Math.f}}=C_{k_1}.Math.\exp \left({}_{k_1}.Math.f_m\right),& \left(40\right)\end{array}$

[0078] where f.sub.m is the center frequency of the transducer. Equation (40) is used to correct the amplitude of the first signal for a potentially negative value of .sub.k.sub.1 by providing a corrected value at the center frequency. .sub.k.sub.1 is then used instead of C.sub.k.sub.1 to estimate Z in eq. (35) (Note that because R.sub.k.sub.1=C.sub.k.sub.1 when .sub.k.sub.1=0, eq. (40) is always used.)

[0079] D. Nonlinear Autoregressive (NLAR) Inverse Model

[0080] This describes a refinement of the AR inverse model in the case of acoustic attenuation that is not linear with frequency.

[0081] 1) Affine Attenuation

[0082] We start from the following approximation for attenuation:

(f)=.sub.0+.sub.1f,(41)

[0083] where .sub.0 is not equal to zero as assumed above. In this case, simple algebraic manipulations can be used to estimate .sub.0 and to establish that .sub.1 can be estimated using eq. (36) because

[00018]$\begin{array}{cc}\frac{C_{k_2}.Math.S_{k_2}\left(f\right)}{S_0\left(f\right)}.Math..Math.\begin{array}{c}=.Math.C_{k_2}.Math.\exp \left(2.Math..Math..Math.f\left({}_{k_2}+j\left(t_{k_2}-t_0\right)\right)\right)\\ =.Math.C_{k_2}.Math.\exp \left(2.Math.\left[2.Math.d\left({}_0+{}_1.Math.f\right)\right]+2.Math..Math..Math.\mathrm{fj}\left(t_{k_2}-t_0\right)\right)\\ =.Math.C_{k_2}.Math.\exp \left(2.Math.\left[2.Math.d.Math..Math.{}_0\right]\right).Math.\exp \left(2.Math..Math..Math.f\left(2.Math..Math.d.Math..Math.{}_1+2.Math..Math..Math.\mathrm{fj}\left(t_{k_2}-t_0\right)\right)\right)\\ =.Math.C_{k_2}^{*}.Math.\exp \left(2.Math..Math..Math.f\left({}_{k_2}^{*}+j\left(t_{k_2}-t_0\right)\right)\right),\end{array}& \left(42\right)\end{array}$

[0084] where

C.sub.i.sub.2*=C.sub.i.sub.2 exp(2[2d.sub.0])(43)

.sub.i.sub.2*=2d.sub.1.(44)

[0085] Therefore, Eq. (44) confirms that .sub.1 can be directly obtained from Eq. (36). To obtain .sub.0, we use Eq. (37) with the value of Z found from Eq. (35) to obtain C.sub.i.sub.2. Because of the identity of Eq. (42), Eq. (20) yields C.sub.i.sub.2*. Equation (43) then can estimate .sub.0.

[0086] 2) Power-Law Attenuation

[0087] Finally, although the affine law of Eq. (41) is often used to obtain an attenuation approximation when attenuation values are known over a finite bandwidth, in this final section, we consider the power-law model of attenuation as expressed by:

(f)=f.sup.,(45)

[0088] To estimate and , we use the affine approach over m distinct frequency bandwidths (i.e., BW.sub.1{BW.sub.1, BW.sub.2, . . . , BW.sub.m}) that are composed of frequencies BW.sub.l=f.sub.1.sup.1, . . . , f.sub.N.sub.1.sup.1, BW.sub.2=f.sub.1.sup.2 . . . , f.sub.N.sub.2.sup.2, . . . , BW.sub.m=f.sub.1.sup.m, . . . , f.sub.N.sub.m.sup.m. This procedure yields .sub.0.sup.l and .sub.1.sup.l where the superscript l specifies which bandwidth was used to obtain these estimates. Then, we derive the following equations for all f.sup.l included in each BW:

log()+ log(f.sup.l)=log(.sub.0.sup.l+.sub.1.sup.lf.sup.l),(46)

[0089] which can be written in the following matrix form:

M[log(),]=T,(47)

[0090] where

[00019]$\begin{array}{cc}M=\left(\begin{array}{cc}1& \log \left(f_1^1\right)\\ \mathrm{.Math.}& \mathrm{.Math.}\\ \mathrm{.Math.}& \log \left(f_{N_1}^1\right)\\ \mathrm{.Math.}& \mathrm{.Math.}\\ \mathrm{.Math.}& \log \left(f_1^m\right)\\ \mathrm{.Math.}& \mathrm{.Math.}\\ 1& \log \left(f_{N_m}^m\right)\end{array}\right);T=\left(\begin{array}{c}\log \left({}_0^1+{}_1^1.Math.f_1^1\right)\\ \mathrm{.Math.}\\ \log \left({}_0^1+{}_1^1.Math.f_{N_1}^1\right)\\ \mathrm{.Math.}\\ \log \left({}_0^m+{}_1^m.Math.f_1^m\right)\\ \mathrm{.Math.}\\ \log \left({}_0^m+{}_1^m.Math.f_{N_m}^m\right)\end{array}\right);& \left(48\right)\end{array}$

[0091] which is solved using a least-squares matrix formulation to obtain (from log()) and .

[0092] In the current implementation, we used three bandwidths (i.e., m=3): BW1 was the full 20-dB bandwidth of the transducer extending from f1 to f2 (i.e., the same as used in the AR approach), BW2 extended from f1 to f21/11(f2f1), and BW3 extended from f1 to f22/11(f2f1).

[0093] We limit our presentation below of results of the NLAR approach only to results obtained using the power-law attenuation model. Therefore, the NLAR approach always is discussed in terms of the the power-law attenuation method described in the present section.

[0094] 3. Material and Methods

[0095] A working QAM system has been designed, fabricated, tested, and applied operating at frequencies ranging from 100 to 500 MHz and, although all the methods described in the Section II above are applicable over a wide range of frequency, the remaining discussion pertains to the QAM system discussed and equipped with a broadband transducer operating a center frequency of 250 MHz.

TABLE-US-00002 TABLE 2 Simulated parameter variation Parameter Lower value Step size Upper value SNR [dB] 40 10 100 d [m] 3 1 8 Z [MRayl] 1.51 0.01 1.6 [dB/MHz/cm] 8 1 13 1 0.1 1.5

[0096] A. Simulations

[0097] Simulations have been conducted to evaluate performance of the inverse models. To mimic our experiments closely, we used a measured reference signal obtained using the apparatus described in section 111.2 below. Simulations consisted of setting values for d, c, Z, , and and reconstructing the simulated signal s.sup.Sim(t) accordingly:

s.sup.Sim(t)=C.sub.1s.sub.1.sup.Sim(t)+C.sub.2s.sub.2.sup.Sim(t)+C.sub.perts.sub.0*(t.sub.t.sub.p.sub.ert)+(t),(49)

[0098] where the first two terms are the simulated reflections from the water-tissue and tissue-glass interfaces whose defining parameters are obtained directly from the set parameters. The third term is a perturbation term having a shape similar to the first two terms, which can be turned off by selecting C.sub.pert=0. The last term is a white Gaussian noise of power .sup.2 with

=10.sup.(20 log .sup.10.sup.(((s.sup.0.sup.)SNR)/20));(50)

[0099] where SNR is the signal-to-noise ratio expressing the difference in dB between the amplitudes of the reference signal and noise. The symbol (expresses the maximum of the hilbert transform of s.sub.0. Equation (49) depends on a very large number of parameters; therefore, we limited the range of parameter variations to experimentally-relevant values. We tested the effects of decreasing the SNR, decreasing the signal separation (i.e., sample thickness, d), decreasing the amplitude of the first signal (i.e., acoustic impedance, Z), increasing the attenuation coefficient (i.e., , see eq. (45)), and increasing the frequency exponent (). Table 2 gives a summary of all parameter-value ranges used in the simulations. These ranges were selected to be representative of realistic scenarios for QAM-applications and were based on preliminary tests to find the optimal range between easily separable cases (i.e., large SNR, d, Z, small and =1) and difficult cases (i.e., small SNR, d, Z, large and >1). In each simulation scenario, the value of the parameter under investigation was varied and all the remaining parameter values were kept constant with SNR=60 dB, d=8 m, Z=1.63 MRayl, c=1600 m/s =10 dB/cm at 250 MHz, C.sub.pert=0 and =1. To assess statistical variations, each case was performed for 200 realizations of (t). This procedure required (7+6+10+6+6)*200=7000 simulations.

[0100] To investigate the impact of a third signal, we varied C.sub.pert from 0.Math.C.sub.s.sub.1 to 0.9.Math.C.sub.s.sub.1 (increment 0.1.Math.C.sub.s.sub.1) and randomly placed the third signal between s.sub.1 and s.sub.2 (i.e., t.sub.1t.sub.pertt.sub.2). This procedure ensured that C.sub.pert was always smaller then C.sub.s.sub.1 and C.sub.s.sub.2. The other acoustic parameter values were kept constant at SNR=60 dB, d=8 m, Z=1.63 MRayl, c=1600 m/s =10 dB/cm.

[0101] To assess the performance of the inverse models, we calculated the mean error and standard deviation of estimated parameters (i.e., c, Z, , and ) with the simulated parameter. In addition, we used the Grubb's test to detect and report outliers in terms of percentages to provide another metric to compare the performance of the AR- and Hozumi-model approaches in the simulation experiments.

[0102] B. Experiments

[0103] The inverse models were tested using experimental data, which were collected with our previously described 250-MHz QAM system. The device was equipped with a 250-MHz transducer and signals were digitized at 2.5 GHz with 12-bit accuracy. The data were acquired from a 12-m thin human lymph node and a 6-m thin section of a human cornea. The thickness of the samples was chosen to be thin enough that scattering effects are strongly mitigated. The specimens were raster scanned in 2D with a 2-m step size in both directions. RF were data acquired at each location and processed individually using the inverse models. 2D maps of d, c, Z, and attenuation were generated. Adjacent 3-m thin section were stained with H&E and digitally imaged at 20 to provide a reference for tissue microstructure. Outliers in experimental data were detected by using absolute thresholds for values of c and Z parameters because the Grubb's test could not be applied. Thresholds were selected based on results of our previous studies. Specifically, estimates were rejected if c<1500, c>2200, Z<1.48, or Z>2.2.

[0104] 4. Results

[0105] A. Simulations

[0106] 1) General Simulations

[0107] FIG. 2 shows the results of parameter variation and resulting error in estimating acoustic-properties values, which are c (i.e., column 1), Z (i.e., column 2), and (i.e., column 3). Column 4 shows the percentage of outliers for each model.

[0108] Columns one to three of FIG. 2, show error of the estimated acoustic property (i.e., Speed of Sound (c), Acoustic Impedance (Z), and Attenuation coefficient ()). Each row shows different parameter-set variations for signal to noise ratio (SNR), separation (i.e., d, thickness of the specimen), amplitude of the first signal (i.e., Z), attenuation coefficient (i.e., ), and attenuation exponent (i.e., ). Blue error bars show error and std of Hozumi-based estimation and black bars show error and std of AR-based estimation. Gray bars show results if a is estimated using the non linear estimation procedure (see eq. 42). The fourth column shows the outliers that in % that were removed from the parameter estimation. Mean and standard deviations were obtained after outliers were removed. Blue and gray bars show results for Hozumi and AR model, respectively.

[0109] The inverse AR and Hozumi models performed well (i.e., average errors were approximately zero for all three properties and small number of outliers) in the easy cases (i.e., large SNR, d, and Z and small and =1). Both methods were stable in estimating c, Z, and a down to 50-dB SNR, which is lower than values found in our measurements (i.e., an SNR60 dB was obtained in an experimental water-glass reflection signal).

[0110] However, in the harder cases, variance, bias, and the number of outliers, of the Hozumi model were significantly greater than those of the AR model. This difference in the performance of the two inverse models is particularly, apparent for less-separated signals (i.e., variation in d, second row FIG. 2) and for small amplitudes of the first signal (i.e., variation of Z, third row FIG. 2). The Hozumi model is unable to estimate c and Z when d is smaller than 6 m and Z is smaller than 1.56 MRayl. (In FIGS. 2-4 the absence of a data-point and error bars indicate 100% outlier as shown in column 4, i.e., the two parameters could not be estimated and the two signals could not be detected in all 200 cases). The AR model remains stable (i.e., c has an average error of 0 m/s and Z has an average error of 0 MRayl and the variation in c is <10 m/s and in Z is <0.05 MRayl) down to a d of 4 m and a Z of 1.52 MRayl.

[0111] Simulated variations in a had a limited impact on estimates of c and Z (fourth row of FIG. 2), although the Hozumi model showed a trend toward increasing negative bias in c and Z when increases. In addition, results clearly indicate that the NLAR and AR models significantly outperform the Hozumi model in estimating . The simulation results also indicate that the NLAR model yielded smaller variance (even though the simulated attenuation was linear). Furthermore, the AR and NLAR models did not produce any outliers, whereas the Hozumi model produced an average outlier percentage of approximately 7%.

[0112] 2) Nonlinear Attenuation Results

[0113] Simulating nonlinear attenuation (i.e., >1) had no impact on estimates of c and Z (last row of FIG. 2) obtained using the AR and Hozumi models. However, an important impact on estimating was visible with the Hozumi and AR models showing a strongly increasing positive bias when increased. However, when was estimated using a power-law model for attenuation estimation (see eq. (45)) the NLAR model produced satisfactory estimates of up to =1.4.

[0114] To investigate the NLAR model performance further, FIGS. 3a-d shows errors in estimates as a function of simulated SNR, d, Z, and , respectively =1 (black) and =1.2 (red). The results obtained for =1 indicate that the NLAR model was capable of performing satisfactorily of estimating up to an SNR of 50 dB, d=5 m, Z=1.52 MRayl, =20 dB/MHz/cm. Similarly, results obtained for =1.2 were nearly identical in terms of variance, but a small positive bias of 0.05 was observed in most cases.

[0115] FIG. 3 provides simulation results showing estimation error and std of exponent for varying parameters using the NLAR model. FIG. 3a shows SNR, FIG. 3b shows d, FIG. 3c shows Z, FIG. 3d shows , and using the NLAR-model.

[0116] FIG. 3e shows errors in estimates of as a function of simulated from 1 to 1.5. Results indicate very small variances (i.e., <0.02); however, a slowly increasing positive bias occurs as the simulated value of increases up to a maximum bias of 0.13 when the simulated was 1.5.

[0117] 3) Perturbation-Signal Results

[0118] Estimation results based on simulations that included a perturbation signal showed that, as the amplitude of the perturbation signal increased, estimation errors increased for all four parameters and for all models as shown in FIGS. 4a-4d. (Estimation error and std of c (4a), Z (4b), (4c), and exponent (4d) for varying amplitude of the perturbation signal (C.sub.pert)).

[0119] However, the AR-model outperformed the Hozumi model. Nevertheless, the impact of a perturbation signal on estimates of c and Z was small (FIGS. 4a and 4b) with nearly no bias and with standard deviations smaller than 5 m/s and 0.05 MRayl for all models. The error in Z was nearly zero for all tested amplitudes of C.sub.pert when the AR model was used, whereas it increased with C.sub.pert when the Hozumi model is used.

[0120] The existence of a perturbation signal had an important effect on estimates of and (FIG. 4c-4d). For a, the AR models (i.e., NLAR and AR) significantly outperformed the Hozumi model. In fact, the error bars shown in FIG. 4c in the case of the Hozumi model were greater than 0.5 dB/MHz/cm if C.sub.pert was greater than 0.2, essentially making the estimates unreliable. Finally, estimates obtained using the NLAR model showed an overall trend toward increasing variance as C.sub.pert increased, but the variances remained below 0.1 up to C.sub.pert=0.7 except for the outlier result obtained for C.sub.pert=0.6.

[0121] In summary, the existence of a perturbation signal has limited effects on estimates of c and Z, but strong effects on attenuation parameters. Overall, the AR and NLAR models perform much better than the Hozumi model by producing much smaller variances and more-reliable estimates for nearly all values of C.sub.pert

[0122] 4) Illustrative Simulation Fits

[0123] FIGS. 5a-5j show representative simulated and fitted signals (FIG. 5a-f) using the Hozumi and AR-inverse models, respectively. The depicted signals are for representative cases of the parameter-variation experiments with SNR=50 dB (FIGS. 5a-5b), d=6 m (FIGS. 5c-5d) and Z=1.56 MRayl (FIGS. 5e-5f). These cases were selected at thresholds where the Hozumi or AR models start to fail more frequently (as illustrated in FIG. 2 rows one to three). Both models show moderately good fits even in the low SNR case (FIGS. 5a-5b). However, the AR model shows a smaller estimation error, in particular for the more-challenging case with a small sample thickness (i.e., d=6 m, FIGS. 5c-5d). Although the fit looks reasonable, the errors in estimating Z and c are 25% and 16% of an assumed natural variations of 1.5 to 1.9 MRayl and 1500 to 1800 m/s in soft biological tissues, respectively. Similar results were observed for the variations in signal amplitude (FIGS. 5e-5f) with an underestimation of Z.

[0124] FIGS. 5a-5f show simulated signals and FIGS. 5g-5h show measured signals (i.e., depicted in black) together with the fitted signals (i.e., depicted by red dotted lines) using Hozumi model (i.e., FIGS. 5a,5c,5e,5g and 5i) and AR model (i.e., FIGS. 5b,5d,5f,5h and 5j), respectively. FIGS. 5a and 5b show a signal with a low SNR (i.e., simulated SNR=40 dB, c=1600 m/s, Z=1.63 MRayl, d=8 m), FIGS. 5c and 5d show a signal with small d (i.e., simulated d=6 m, c=1600 m/s, Z=1.63 MRayl), FIGS. 5e and 5f show a signal with small Z (i.e., Z=1.56 MRayl, c=1600 m/s, d=8 m), FIGS. 5g and 5h show a measured signal of the Lymph node sample, and FIGS. 5i and 5j show a measured signal from the cornea sample. Estimated parameters of the signals are given in the figures upper right cornea.

[0125] B. Experiments

[0126] 1) Lymph-Node Results

[0127] FIG. 6-8 depict parameter maps of a 12-m lymph-node section. Parameter maps obtained with both models show similar features and allow distinguishing between different tissue types (i.e., fibrous tissue of the capsule, Ca, and lymph-node parenchyma tissue, LT). For example, Z, c, and values were significantly higher in Ca compared to LT (T-test, p<0.01). Table 3 below summarizes average values for Ca and LT estimated using Hozumi and AR models. Average values of c and Z estimated with the Hozumi model are slightly lower than those estimated using the AR-model. On average, values are higher when estimated with the AR model. However, all differences are within the standard deviations of the parameter variations. Differences are apparent in the structural features. In particular, the variation of c (FIGS. 6e-6f and FIGS. 7c-7d) is more noticeable in images based on the AR model than those based on the Hozumi model, and the images based on the AR model more clearly depict anatomic structures matching those visible in the histology images (FIG. 6b). Note that QAM parameter maps and H&E stained images are obtained from adjacent but different sides, which can result in dissimilarities between H&E and QAM images. The color maps in FIG. 6a-6h are optimized to highlight features in the fibrous tissue of the Ca and FIGS. 7a-7f optimizes the color maps to compare Hozumi and AR-model in LT. No difference in the number of outliers was observed between Hozumi and AR-model. Parameter maps of and generated using the NLAR model are shown in FIGS. 8a-8d. The contrast between Ca and LT is weak for . Values for show no significant difference between Ca and LT.

TABLE-US-00003 TABLE 3 Acoustic parameters of lymphocyte tissue (LT) and capsule (Ca) of lymph nodes and epithelial (Ep) and stroma (St) tissue of cornea samples estimated using the Hozumi and the AR model, respectively Parameter Hozumi AR/NLAR Z.sub.LT [MRayl] 1.71 0.09 1.76 0.10 c.sub.LT [m/s] 1541 12 1547 15 .sub.LT [dB/MHz/cm] 6.4 1.4 5.7 1.0 .sub.LT [dB/MHz] 3.9 1.6 .sub.LT 1.32 0.22 Z.sub.Ca [MRayl] 1.78 0.08 1.84 0.08 c.sub.Ca [m/s] 1598 21 1610 30 .sub.Ca [dB/MHz/cm] 9.5 2.7 8.5 2.5 .sub.Ca [dB/MHz] 5.9 2.8 .sub.ca 1.37. 0.21 Z.sub.Ep [MRayl] 1.58 0.03 1.59 0.03 c.sub.Ep [m/s] 1526 20 1548 22 .sub.Ep [dB/MHz/cm] 2.4 1.5 3.0 1.1 .sub.Ep [dB/MHz] 2.3 1.1 .sub.Ep 1.25 0.15 Z.sub.St [MRayl] 1.55 0.01 1.55 0.03 c.sub.St [m/s] 1604 38 1684 60 .sub.St [dB/MHz/cm] 3.8 1.8 6.6 2.8 .sub.St [dB/MHz] 7.6 3.4 .sub.St 1.0 0.5

[0128] FIGS. 6a-6h show QAM-amplitude maps derived from the maximum amplitude of the envelope signals (FIG. 6a) and H&E stained Histology map (FIG. 6b) of a human lymph node. FIGS. 6c,6e, and 6g show parameter maps generated using the Hozumi model and FIGS. 6d,6f, and 6h show parameter maps generated using the AR model. FIGS. 6c and 6d plot Z, FIGS. 6e and 6f plot c, and FIGS. 6g and 6h show . Two tissue types are identified on the histology section, which are capsule (Ca) and Lymphocyte tissue (LT). FIGS. 6c to 6h are from ROI1.

[0129] FIGS. 7a, 7c and 7e show parameter maps from ROI2 of FIG. 6a generated using the Hozumi model and FIGS. 7b,7d, and 7f show parameter maps generated using the AR model. FIGS. 7a and 7b plot Z, FIGS. 7c and 7d plot c, and FIGS. 7e and 7f show . The depicted tissue type is Lymphocyte tissue.

[0130] FIGS. 8a-8d show parameter maps of (FIGS. 8a and 8c) and (FIGS. 8b and 8d) of a human lymph node of ROI1 (FIGS. 8a and 8b) and ROI2 (FIGS. 8c and 8d) of FIG. 6a derived using the NLAR-model.

[0131] 2) Cornea Results

[0132] The parameter maps of the 6-m cornea sample are shown in FIGS. 9a-9h and the average values are summarized in Table 3 above. Two tissue types (i.e., epithelium or Ep, and stroma or St) were identified with both models. Differences in acoustic properties and structural contrast between Hozumi and AR-model are evident. They are particularly apparent in St. The AR model shows higher c, Z, and values. Furthermore, the features follow the collagen-fiber orientation as seen in the H&E stained Histology image (FIG. 9b) and the Hozumi model failed in additional locations (i.e., white areas in the parameter map) compared to the AR model. The values of c, Z, and in the epithelium are only slightly lower in the Hozumi-based estimates.

[0133] FIG. 9a shows an amplitude map and FIG. 9b shows a histology map of a pig cornea. FIGS. 9c,9e, and 9g show parameter maps generated using the Hozumi model and FIGS. 9d,9f, and 9h show parameter maps generated using the AR model. FIGS. 9c and 9d plot Z, FIGS. 9e and 9f plot c, and FIGS. 9g and 9h show . Two tissue types are identified on the histology section, which are epithelium (Ep) and stroma tissue (St).

[0134] 3) Illustrative Experimental Fits

[0135] Above findings are confirmed by comparing the fitted models with the measured signals as shown in FIGS. 5g-j. The fitted amplitudes of the first signals (i.e., reflections from the tissue-water interface) are underestimated when the Hozumi model is used, which leads to lower estimated Z values. Lower Z values were also observed in experimental results (FIG. 6 e-f and FIG. 7 c-d), and the number of outliers was significantly higher in the thin cornea samples (white pixels in FIGS. 9c, e, and g). Furthermore, simulations suggest that non-linear attenuation (i.e., >1) can cause an overestimation of a if the Hozumi model is used, as shown in FIG. 2 column 3 and row 5. The experimental data show -values larger than 1 for Ca and LT (see Table 3). The Hozumi model shows higher attenuation values compared to the average values obtained using the AR model. However, the St tissue showed an value of approximately 1 and in this case, attenuation estimated with the AR model was almost twice as high when compared to the Hozumi model.

[0136] C. Results Summary

[0137] Overall, experiment and simulation results are consistent. For challenging cases (e.g., small d, small Z, or large ), the Hozumi model underestimates Z (FIG. 2 column 2 and rows 2 to 3) and the percentage of outliers is larger than that of the AR-model (FIG. 2 column 4, rows 2 to 3). Simulations and experimental results show that both methods perform equally well in well separated signals with low, linear attenuation and high Z contrast between the tissue and the coupling medium. However, the AR model outperforms the Hozumi model in terms of parameter-estimation accuracy and stability (i.e., parameter error variation and number of outliers) if (i) the signals overlap (e.g., because of small sample thickness or high c), (ii) the contrast in Z between the tissue and the coupling medium decreases (i.e., smaller amplitude of the first signal), (iii) the attenuation increases, and (iv) in case of non-linear attenuation (i.e., >1).

[0138] V. Discussion and Conclusion

[0139] In this study, we investigated the application of an autoregressive inverse model for estimating acoustical properties of thin, soft-tissue sections at microscopic resolutions. In addition, our extension to better deal with power-law attenuation and the elegant use of Cadzow denoising have never been investigated and provide significant improvements in QAM imaging. The proposed AR model is similar to Prony's method, which was successfully applied in ultrasound research to separate fast and slow waves as observed in ultrasound through-transmission experiments of bone specimens. Furthermore, the inversion method used to obtain the modes of the AR process based on the zeros of a polynomial (see. Eq. (18)), also exists in signal-processing literature and is often termed annihilating filter or polynomial and has applications in sparsity-based methods.

[0140] The AR and Hozumi model are ultrasound-frequency-spectrum approaches. However, time-domain methods were also suggested in SAM and QAM applications and work well when the two signals are completely separated in time. However, the spectral methods (e.g., Hozumi, AR, cepstrum) are the only methods available to separate the two signals when they overlap in time. In fact, our results (FIG. 2, and FIGS. 5a-5j) demonstrate that the AR approach is far superior when the signals overlap significantly.

[0141] The most widely-used inverse methods for QAM-based thin-section assessment remain those based on the approach suggested by Hozumi et al. However, the Hozumi approach has limitations. For example, it only models order-two signals, strongly depends on the transducer bandwidth, relies on peak detection in the frequency domain, and does not allow estimating non-linear attenuation. Therefore, the motivation of the present invention was to improve performance for tissues with (i) acoustic impedance close to water (i.e., signals with small amplitudes and low SNR), (ii) strongly overlapping signals, (iii) high attenuation, and (iv) non-linear attenuation. Our simulation results show better performance (e.g., a smaller parameter-estimation error) and better stability (e.g., a smaller variation of estimation error and fewer outliers) for our AR-inverse model approach in all four investigated cases (FIG. 2) when compared to the Hozumi approach. Particularly relevant to soft-tissue measurements is the better performance of the AR model when Z is close to the impedance of water (i.e., when the amplitude of the first signal is small and the SNR is low). Furthermore, in very thin sections that produce overlapping signals, the AR-model is more stable. This feature is a very important attribute for soft-tissue measurements because it allows performing acoustic microscopy and optical microscopy on the same sample without sacrificing histology-image quality caused by thick samples.

[0142] Interestingly, the present invention provides that non-linear attenuation in ranges typical for soft tissues has only a small impact on c and Z estimates for the Hozumi and AR models. The Hozumi-model shows only a small bias for Z estimates with increasing attenuation. Nevertheless, the Hozumi-inverse model has difficulties estimating with increasing attenuation and (FIG. 2). This occurs because the Hozumi-model estimates Z and values only from the resonance-peak amplitudes in the frequency domain, rather than from the normalized spectrum over the entire bandwidth as it is done for the AR model. The peak locations and amplitudes also are affected by the attenuation which can lead to estimation errors (FIGS. 10a-10c).

[0143] FIG. 10 shows normalized spectra (N(f)) of the standard signal (i.e., Easy, black spectra with SNR=60 dB, d=8 m, Z=1.63 MRayl, c=1600 m/s =10 dB/cm at 250 MHz, C.sub.pert=0 and =1) and spectra with (FIG. 10a) varied d=6 t (blue curve), and Z=1.53 MRayl (red curve); (FIG. 10b) varied =30 dB/cm (blue curve), =1.5 (red curve); and (FIG. 10c) varied SNR=40 dB (blue) and Cadzow de-noised spectra (red).

[0144] The results from ex-vivo-tissue experiments are consistent with those obtained from simulations and demonstrate a better performance of the AR-model as follows:

[0145] 1) Analysis of the parameter-maps of the very thin sections (e.g., the cornea sample, FIG. 9) that were derived using the AR-model revealed anatomical structures that match histology more closely than parameter maps generated by the Hozumi model.

[0146] 2) The Hozumi model failed to separate the two signals in more cases than the AR-model as indicated by the white areas (i.e., outliers) in FIG. 9)

[0147] 3) In the thicker, 12-m lymph-node sections, the parameter maps of the AR-model showed enhanced structural features in the Z and c parameter maps, which is a result of greater robustness and sensitivity of the AR-model to small parameter variations, as shown in the simulation results. The simulations indicate that, in low-Z and -d conditions, the AR-model produces lower variation in parameter-estimation errors and still provides reliable results when the Hozumi-model completely fails (see FIG. 2).

[0148] 4) If a third signal is introduced (FIG. 4), then the AR-model performs better than the Hozumi model, which may be one reason for the better contrast of the parameter maps in the 12-m sections (FIGS. 6-8) obtained from the experimental data. Perturbations can arise from scattering or strong reflection inside the tissue. The impact of scattering is assumed to be small because echoes from the scattering signals would be about one order of magnitude smaller than the sample-tissue and tissue-substrate interface signals, which are planar reflections. Scattering is also mitigated by the use of relatively thin samples (i.e., around the wavelength). In fact, our extensive experimental data from our previous studies did not show any larger signal between the signals of the two interface reflections, which led to the assumption of our current approach that assumes the largest amplitudes as the reflections from the water-tissue and tissue-substrate reflections. The AR approach can be tailored to pick the right two signals out of several signals by using another criterion than amplitude (e.g., time of flight, attenuation, etc.), but the existence of a prominent middle signal would also mean that our forward model (i.e., single layer) is most likely violated and estimated parameters would most likely need to be excluded. Multilayer approaches could also be taken care of using a different version of AR model.

[0149] A striking advantage of the proposed AR approach is its ability to estimate nonlinear attenuation based on a power-law model, which has not been reported to the best of the authors' knowledge. If the power-law model (i.e., the NLAR model) is used, estimating was significantly improved compared to the Hozumi model, and the AR model shows stable results (i.e., average error in estimating 0 dB/MHz/cm and error-variation<0.1 dB/MHz/cm) over the entire range of simulated values and also can estimate the exponent correctly (i.e. error variation of 0.4) over the simulated parameter range (FIG. 3). Biological tissues typically exhibit non-linear attenuation, but assessing ultrasound attenuation at microscopic resolutions has been challenging. However, our experimental results demonstrate the high value of the non-linear attenuation parameter of the lymph-node sample (FIG. 8). These results demonstrate the usefulness of the AR method. Accurate measurements of the non-linear attenuation also allow to assess the dispersion effects that can occur in soft tissues and can lead to speed of sound values that vary slightly with ultrasound frequency. Using causality relationships, the increase can be exactly quantified using the derivative of the frequency-dependent attenuation.

[0150] In addition, higher-quality information provided by the AR-model approach is important for modeling needs and for ultrasound applications at lower frequencies. Currently, QAM is the only method that can provide multiple acoustical and mechanical properties at fine-resolutions and over large-scale areas as is required for numerical modeling of sound propagation or finite element modeling. Such properties cannot be assessed using conventional optical microscopy methods. Use of advanced computer simulations is rising, which allows investigating complex phenomena that are difficult or cannot be examined experimentally. However, the results of these simulations will only be as accurate as the underlying models whose accuracy in turn depends on realistic input data. Furthermore, in many quantitative-ultrasound (QUS) applications, assumptions are made about the acoustic attenuation of tissue to correct QUS-parameter estimation. We hope that QAM will provide better estimates of common tissues to improve novel QUS-methods using acoustic-impedance-map methods, for example. The methods provided in the present invention to separate ultrasound signals and to estimate acoustical-parameter values also are suitable for applications at lower frequencies.

[0151] The present invention, a new AR-model, is compared with the current standard method (i.e., Hozumi-model). Although some of the improved performance of the AR-model may be directly related to implementation details, the AR model has some general advantages that are illustrated in FIG. 10, which shows representative normalized spectra for different simulated signals. The Hozumi model relies on the detection of resonance peaks in the frequency domain. In well-separated signals (FIG. 10a, Easy) with strong amplitudes, the normalized spectra exhibit several peaks that are easily detected. In theory, one maximum-minimum pair is sufficient to estimate d and c values, but additional peaks increase robustness (i.e., decrease bias and variance) to noise by averaging estimates obtained from all possible extremum pairs. Therefore, the Hozumi model is sensitive to the bandwidth of the QAM system (FIG. 10a, d), which determines the number of peaks that can be analyzed. In contrast, the AR model uses the entire normalized spectra and not only the peaks. Thus, the AR-model is more robust at the bandwidth limits (i.e., small d), as indicated by the simulation results. Similarly, the Hozumi model depends strongly on the resonance-peak amplitudes. If the amplitudes are too small, peak-detection methods may fail to detect the resonance peaks accurately. This is the case if contrast in Z between tissue and coupling medium is very small (FIG. 10a, Z). Strong attenuation (linear and non-linear) produce similar effects on the peak detection (FIG. 10b) and explain the better performance of the AR-model as observed in the simulation results. Furthermore, although not discussed in great detail here in the interest of space, the performance of the AR-model is strongly enhanced by the Cadzow denoising. FIG. 10c shows an example of a noisy signal before (FIG. 10c, SNR) and after de-noising (FIG. 10c, Cadzow) compared to a noise-free signal.

[0152] To summarize, our AR-model approach for QAM-based parameter estimation showed better performance for four highly-relevant scenarios: (i) tissues with acoustic impedance close to water, (ii) tissue samples yielding overlapping signals and (iii) tissues with nonlinear attenuation. We demonstrated in experiments and simulations the improved robustness and precision of acoustic-parameter estimation of the AR-model (e.g., smaller variation and bias of errors for samples with Z1.56 MRayl, d6 m, 5 dB/MHz/cm, and 1). Furthermore, the AR method is easily implemented and allows direct estimation of all acoustic properties including those related to non-linear attenuation. Another advantage is the unique ability of the AR model to remove spurious signals, such as perturbation signals (FIG. 4), and still provide accurate estimates of acoustic properties. Our ex-vivo experiments demonstrate that average acoustic-parameter values are similar for both methods in 12-m sections. However, the AR model seems to show more contrast for thinner tissue sections. This would allow to compare QAM measurements at different frequencies (e.g., 250 MHz and 500 MHz) on the same tissue samples. Furthermore, we believe that our new AR-model approach in general can markedly improve current QAM technology.

[0153] A computer system suitable for storing and/or executing program code for the present invention includes at least one processor coupled directly or indirectly to memory elements through a system bus. The memory elements include local memory employed during actual execution of the program code, bulk storage, and cache memories that provide temporary storage of at least some program code to reduce the number of times code is retrieved from bulk storage during execution. Input/output (I/O) devices (including but not limited to keyboards, displays, pointing devices, etc.) can be coupled to the computer system either directly or through intervening I/O controllers. Network adapters may also be coupled to the computer system in order to enable the computer system to become coupled to other computer systems or remote printers or storage devices through intervening private or public networks. Modems, cable modems, and Ethernet cards are just a few of the currently available types of network adapters. The computer system can also include an operating system and a compute file-system.

[0154] Although the present invention has been described in conjunction with specific embodiments, those of ordinary skill in the art will appreciate the modifications and variations that can be made without departing from the scope and the spirit of the present invention. Such modifications and variations are envisioned to be within the scope of the appended claims.