Signal-Processing Framework for Ultrasound Compressed Sensing Data: Envelope Detection and Spectral Analysis

: Acquisition times and storage requirements have become increasingly important in signal-processing applications, as the sizes of datasets have increased. Hence, compressed sensing (CS) has emerged as an alternative processing technique, as original signals can be reconstructed using fewer data samples collected at frequencies below the Nyquist sampling rate. However, further analysis of CS data in both time and frequency domains requires the reconstruction of the original form of the time-domain data, as traditional signal-processing techniques are designed for uncompressed data. In this paper, we propose a signal-processing framework that extracts spectral properties for frequency-domain analysis directly from under-sampled ultrasound CS data, using an appropriate basis matrix, and e ﬃ ciently converts this into the envelope of a time-domain signal, avoiding full reconstruction. The technique generates more accurate results than the traditional framework in both time- and frequency-domain analyses, and is simpler and faster in execution than full reconstruction, without any loss of information. Hence, the proposed framework o ﬀ ers a new standard for signal processing using ultrasound CS data, especially for small and portable systems handling large datasets.


Introduction
"Big Data" is an increasingly ubiquitous term that refers to the growing size and complexity of datasets due to the variety of sources from which they can be acquired. Big Data is transforming the fields of science, engineering, medicine, healthcare, finance, business, and society in general. As the sizes of datasets increase, efficient compression techniques are desperately required for acquisition, storage, and reliable retrieval of data. Compressed sensing (CS) has emerged as an alternative signal-processing solution that addresses the requirements for reductions in data storage and acquisition time, since, assuming they are sufficiently sparse, original signals can be reconstructed using fewer data samples acquired at frequencies lower than those dictated by the Shannon-Nyquist theory [1][2][3]. Accordingly, the CS technique has been applied in many fields, including imaging, communications, pattern recognition, sound processing, and video processing [4][5][6][7][8][9][10]. The CS technique also offers promising benefits for medical image processing applications [7,[11][12][13]. For instance, the relatively long data acquisition times for MRI systems, which make patients uncomfortable as they are required to remain motionless for prolonged periods, is dramatically reduced with CS [14]. Similarly, many recent studies have focused on reducing the size of the hardware and data storage of portable systems used for medical ultrasound, by lowering the sampling rate using the CS theory [15]. 2 of 16 In general, when CS signals are processed for additional analysis, reconstruction is required to convert them to their original forms, as most traditional methods for signal analysis utilize uncompressed signals acquired in the time domain. In particular, with medical imaging analysis, the envelope of the received signal in the time domain contains important anatomical information, while, in the frequency domain, spectral information provides fundamental background for quantitative analysis of tissue properties. For the applications of nuclear magnetic resonance spectroscopy, the CS reconstruction contains the spectral information since the sparse transformation is the Fourier transform itself [16,17]. Similarly, Shim et al. showed that the traditional processing framework for frequency-domain analysis of compressed medical ultrasound data contains a redundant transformation from the time domain to the frequency domain [18]. We have observed that the same redundancy exists with time-domain analyses such as envelope detection for time-varying signals for B-mode ultrasound images.
In this paper, we propose a signal-processing framework for ultrasound CS data that is simpler and faster in execution than the traditional framework, and produces more accurate results in both time-and frequency-domain analyses. With the traditional framework, CS signals first need to be recovered to the uncompressed forms and then transformed from the time domain into the frequency domain using a technique such as the fast Fourier transform (FFT), to estimate the spectral information required for frequency-domain analysis. Here, when an appropriate basis matrix is used, this spectral information is directly extracted from the CS data, avoiding reconstruction of the original signal and the redundant transformation of the reconstructed time-domain signal to the frequency domain. In addition, time-domain properties like the envelope of the original signal are computed during the reconstruction process without fully converting them into their original forms. More specifically, because the envelope of a signal is defined by its magnitude in a complex analytic form, the traditional framework uses the Hilbert transform to obtain analytic time-domain signals for envelope detection. In contrast, with the proposed method, the analytic signal is calculated directly from the frequency information obtained during the reconstruction process, using the inverse Fourier transform (IFFT).
To validate our framework, we performed envelope detection on ultrasound CS data to form a B-mode image, for time-domain analysis. Similarly, for frequency-domain analysis, we estimated ultrasonic attenuation coefficients from the same CS data using the centroid downshift method (CDM) [19]. Additionally, we compared the estimation performances of the discrete Fourier transform (DFT) and discrete cosine transform (DCT), to select the proper basis matrix for the measurement and reconstruction processes. We used several numerical tissue-mimicking (TM) phantoms with different attenuation and backscatterer properties to simulate the propagation of radiofrequency (RF) medical ultrasound signals, for further validation of the operation of the proposed framework in practical experiments.
The rest of the paper is organized as follows: in the next section, the CS technique and the CDM are briefly summarized, as well as the theoretical background of envelope detection using Fourier transform coefficients. We then introduce details about the proposed framework, experimental procedures, and numerical phantoms used in the experiments. In Section 3, we discuss the results of envelope detection, comparing the performances of the two basis matrices at a variety of compression rates, and then estimate attenuation coefficients using the DCT basis matrix. We also compare the traditional and proposed methods with respect to their computational times. In Section 4, we discuss several aspects of the proposed framework and introduce a complete signal-processing framework for medical ultrasound CS data, introducing the new concept of a Q-mode image that simultaneously presents qualitative and quantitative information from the time and frequency domains. The contributions of this paper to the literature are summarized in the final section.

Compressed Sensing Technique
CS enables the reconstruction of under-sampled discrete-time signals when the original signals satisfy the relevant conditions for sparsity and incoherence [1]. For the former condition, a signal is considered to be K-sparse if the number of non-zero values is K [1].
The core of CS theory is finding a solution to the following equation: where y is an (M × 1) vector containing the compressive measurements from the original signal, x is the original (N × 1) discrete-time signal, and A is the (M × N) measurement matrix, in which random sampling is generally modeled with the Gaussian random distribution or the Bernoulli random distribution. If x is K-sparse, there is a high probability that it can be completely recovered from y if M , which is the number of compressive measurements, satisfies the condition that M ≥ O(K log(N/K)) . While x is not necessarily sparse in its original form, specific transformations exist to make it sparse in other domains. When these are applied, x is expressed as follows: where B is an (N × N) basis (transformation) matrix, and the (N × 1) vector, f , is the coefficient for the specific basis matrix. To apply CS theory, a basis matrix should be chosen such that f satisfies the K-sparse condition. In this case, with the matrix T = AB −1 , Equation (1) can be rewritten as follows: In reconstruction, f is first estimated from y , before x is recovered using Equation (2). In general, Equation (3) has an infinite number of solutions, since this problem is underdetermined. However, if T satisfies the restricted isometry property (RIP) condition, it becomes possible to determine a unique solution using optimization techniques. This condition can be expressed as follows: where δ 2K is the RIP constant, located in the range 0 < δ 2K < 1 . Several different optimization approaches exist for determining a unique solution to Equation (3), including l 0 -minimization and l 1 -minimization. While l 0 -minimization is not practical, since it is a nondeterministic polynomial-time hard problem [20], l 1 -minimization becomes a convex optimization problem by solving the linear programming in Equation (5) [21,22]:

Centroid Downshift Method
Of the various medical imaging modalities, ultrasound pulse-echo systems have been widely utilized for visualizing the human body's internal organs, as well as for evaluating the pathological state of soft tissue. To analyze the status of scanned tissue, several ultrasonic parameters, including speed of sound [23,24], integrated backscatter [25,26] and scatterer size [27], are extracted from RF data. Ultrasound attenuation is another critical parameter for evaluation of soft tissue, used to distinguish malignant from benign tumors, as well as provide fundamental information for quantitative ultrasound (QUS) [28].
To simplify representation in the frequency domain, most attenuation estimation methods for ultrasound pulse-echo systems assume constant sound speed and linear frequency dependence in a region of interest (ROI). After applying the Born approximation, which ignores multiple weak scattering of propagating acoustic waves in soft tissue, the intensity of the reflected RF signal in the frequency domain, R( f , z) , can be represented as a product of pulse characteristics (G( f )), diffraction effects (D( f , z)), and backscatterer (B( f )) and attenuation terms (A( f , z)) [29]: where z denotes the distance between the face of a transducer and the center of the ROI. Here, G( f ) , represents the transmit pulse characteristics which are closely related to the geometrical design and beam shaping parameters of the transducer. Insonifying pulses are typically modeled as Gaussian in the spectral domain, which is an approximation that has shown to be accurate within 10%, from root-mean-square-error analysis [30]. Most array transducers in modern ultrasound systems utilize beam focusing techniques to improve the spatial resolution at the ROI. As a result, the propagating acoustic waves experience the diffraction effects, expressed as D( f , z) , along the depth of the tissue, and reflected RF signals acquired by linear array transducers contain diffraction effects around the pre-and post-focal sections. The spectral properties of the power spectra are also distorted before and/or after the beam focus. Fink et al. [31] showed, using classical sound theory and focused transducer models, that diffraction effects could underestimate (or overestimate) attenuation coefficients at pre-(or post-) focal regions. To compensate for the diffraction effects, we utilize a simple correction method using two reference TM phantoms whose attenuation properties are known a priori [28]. Since the reflected RF signals for reference and sample phantoms acquired with the same transducer settings (i.e., the same focus, TGC, and beam pattern) exhibit identical diffraction effects, these can be eliminated by comparing the changes to the center frequency at the same depth.
Because of the large number of randomly distributed small scatterers in the ROI, the frequency-domain characteristics for backscattered echo signals, B( f ) , are expressed as a stochastic model [32]. In general, this backscatterer model is typically in a form that considers only nonlinear frequency dependence, which is incorporated as a power of frequency.
The cumulative attenuation through the beam path is modeled as being dependent on frequency and the distance. In a medium with uniform attenuation characteristics, A( f , z) , is given as, where z is the distance from the transducer to the ROI and α( f ) is the frequency-dependent attenuation expressed in Neper/cm (or dB/cm). Since, α( f ) is known to exhibit linear dependence on frequency [33,34], this function can be written as α( f ) = β· f , where β is the attenuation slope, also called the attenuation coefficient, expressed in Neper/cm/MHz (or dB/cm/MHz). In medical ultrasound systems, estimation of the attenuation coefficients is based on the frequency-dependent attenuation in soft tissue. The power spectrum of RF signals reflected at a given depth demonstrate a downshift to lower frequencies with respect to the position of the center frequency (or centroid), since higher-frequency components attenuate faster than lower-frequency components. Previous research in the literature has found that the attenuation coefficient in soft tissue is almost linearly proportional at frequencies up to 10 MHz. Due to this linearity of the attenuation at higher frequencies, the downshift in center frequency with respect to depth is proportional to the attenuation coefficient, β . Assuming that a Gaussian-shaped transmit pulse is propagated, the relationship between the slope of this downshift with respect to depth and the attenuation coefficient is given by [19]: where z is the distance from the transducer to the ROI, σ 2 is the variance of the transmit pulse, and f c (z) denotes the center frequency of the power spectrum at depth z .

Envelope Detection from Fourier Transform Coefficients
The envelope of a signal in the time domain represents the instantaneous amplitude of a time-varying signal. In medical imaging applications, this envelope conveys structural or anatomical information along the signal path. In general, an analytic expression of a signal is first derived to obtain the envelope of a time-domain signal. These analytic expressions, a widely used concept in signal processing, are alternative complex-valued representations of real-valued signals that make various processes related to a signal's time dependence conceptually simpler [35]. In the time domain, the real part of the analytic signal, x a (t) , is simply the original signal, x(t) , and the new imaginary part, x h (t) , is calculated using the Hilbert transform (x h (t) = HT x(t) ), which suppresses the negative frequency components, i.e., x The mathematical definition of the Hilbert transform of x(t) , x h (t) , is given by: where p.v. is the Cauchy principal value of the integral. Figure 1 shows the spectral properties of the original and analytic signals in the frequency domain. Then, the envelope of a signal, E , is defined as the magnitude of the analytic signal as below, Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 16

Envelope Detection from Fourier Transform Coefficients
The envelope of a signal in the time domain represents the instantaneous amplitude of a timevarying signal. In medical imaging applications, this envelope conveys structural or anatomical information along the signal path. In general, an analytic expression of a signal is first derived to obtain the envelope of a time-domain signal. These analytic expressions, a widely used concept in signal processing, are alternative complex-valued representations of real-valued signals that make various processes related to a signal's time dependence conceptually simpler [35]. In the time domain, the real part of the analytic signal, x ( ) , is simply the original signal, x( ) , and the new imaginary part, x ( ) , is calculated using the Hilbert transform ( x ( ) = HT{x( )}), which suppresses the negative frequency components, i.e., The mathematical definition of the Hilbert transform of x( ) , x ( ) , is given by: where p.v. is the Cauchy principal value of the integral. Figure 1 shows the spectral properties of the original and analytic signals in the frequency domain. Then, the envelope of a signal, , is defined as the magnitude of the analytic signal as below, An alternative approach that is valid for most engineering applications has been proposed for calculating discrete-time domain analytic signals, based on the spectral properties shown in Figure 1 [36]. Since, from the spectrum, an analytic signal has no negative frequency components, and its positive frequency components are double the size of that of the original signal, it can be estimated An alternative approach that is valid for most engineering applications has been proposed for calculating discrete-time domain analytic signals, based on the spectral properties shown in Figure 1 [36].
Since, from the spectrum, an analytic signal has no negative frequency components, and its positive frequency components are double the size of that of the original signal, it can be estimated based on a DFT of the original signal. In this procedure, X[m] represents the coefficients of the N-point DFT of the original signal, x[n] , defined as, The coefficients of the DFT of the analytic signal, Z[m] , can thus be obtained from X[m] as follows: The analytic signal in the discrete-time domain can subsequently be recovered from Z[m] using the inverse DFT. Finally, the envelope is also directly calculated from the real and imaginary parts of the complex discrete-time analytic signal, x[n] and x h [n] , respectively, as follows:

Signal-Processing Framework for Ultrasound Compressed Sensing Data
Although CS signals are advantageous for acquisition and storage processes in terms of the time and the space they save, most signal-processing frameworks used for temporal and spectral analysis require an uncompressed time-domain signal, and cannot be directly applied to CS data. The signal-processing framework for ultrasound CS data we propose in this paper involves the direct computation of the analytic representation of a measured signal, and estimation of its spectral properties without full reconstruction in the time domain. This is possible due to CS theory, which dictates that a basis matrix can be used to convert original signals into sparse signals. If an appropriate basis matrix is selected, the sparse signal's coefficients will contain the spectral information necessary for analysis, avoiding the need for reconstruction and transformation [18].
In Figure 2, two signal-processing procedures for CS data are summarized, with the blue dotted line representing the traditional framework, and the red solid line representing the proposed framework. For simplicity, a signal is recovered from sparse data using the l 1 -minimization algorithm. For envelope detection with the traditional method, the original time-domain signal ( ∼ x ) is first reconstructed from CS data (y). The Hilbert transform is then applied to calculate the analytic signal ( x . In contrast, with our proposed framework, spectral information (f ) obtained by a proper basis matrix during the signal recovery process is directly used for frequency-domain analysis, where attenuation coefficients are identified using the CDM, while a simple IFFT is performed to obtain the analytic signal (x a ) for detection of the envelope of the original time-domain signal. As there are no redundant transformations between time and frequency domains, shorter processing times for both temporal and spectral analyses are expected. In addition, we expect more accurate signal-processing results due to a reduced accumulation of round-off errors.

Numerical Phantoms
To verify the framework, we generated numerical RF signals with differing attenuation properties, for stimulation of the reference and sample TM phantoms, using a frequency-domain model based on the classical linear diffraction theory of continuous waves [38]. In this study, transmitters and receivers were modeled as linear array transducers consisting of 128 rectangular 0.2 × 10 mm elements (azimuthal direction × elevation direction), with a 0.2-mm center-to-center separation. Each element in the array generates an individual beam line. Hence, the array's collective field is a superposition of simultaneous excitation from single elements. Dynamic focusing was applied for reflected RF signals in combination with a dynamic aperture, and the F-number was fixed at 2 during the reception process. Each transmit pulse was modeled as a Gaussian function with a center frequency of 5 MHz and a bandwidth of 80%. In this simulation, a complex wave number is used to calculate parameters related to attenuation and sound speed during wave propagation.
The numerical phantoms used in this study are shown in Figure 3. All phantoms were 34 mm wide, 80 mm deep, and 10 mm thick, and the interval between simulated beam lines was set to 0.2 mm. We assumed that the face of the transducer was in direct contact with the phantoms, and set both the axial transmit focus and elevation focus to 80 mm. We assumed a uniform sound speed of 1540 m/s through the entire depth of the phantoms, while distributions of 50-μm polystyrene beads were included to model backscatter. All simulations were performed in 0.1-MHz increments in the 1 to 15 MHz range. Figure 3a,b depict the sample phantoms used in the experiments. Each phantom includes a cylindrical region with contrasting backscatter properties to the background medium, generated by In our CS framework, the Gaussian random matrix was utilized as a measurement matrix to randomly sample the signal. The Basis Pursuit was a principle for decomposing a signal into a sub-optimal superposition of dictionary elements which have the smallest l 1 -norm of coefficients among all decompositions, and used as the reconstruction algorithm in the proposed framework [37]. For correct operation of the devised framework, the coefficients of the basis matrix should contain physically relevant spectral information, a fact that restricts the possible candidates for selection. To verify the effect of the basis matrix on the framework's operation, we compare the performances of the DFT and DCT matrices when applied for measurement and reconstruction.

Numerical Phantoms
To verify the framework, we generated numerical RF signals with differing attenuation properties, for stimulation of the reference and sample TM phantoms, using a frequency-domain model based on the classical linear diffraction theory of continuous waves [38]. In this study, transmitters and receivers were modeled as linear array transducers consisting of 128 rectangular 0.2 × 10 mm elements (azimuthal direction × elevation direction), with a 0.2-mm center-to-center separation. Each element in the array generates an individual beam line. Hence, the array's collective field is a superposition of simultaneous excitation from single elements. Dynamic focusing was applied for reflected RF signals in combination with a dynamic aperture, and the F-number was fixed at 2 during the reception process. Each transmit pulse was modeled as a Gaussian function with a center frequency of 5 MHz and a bandwidth of 80%. In this simulation, a complex wave number is used to calculate parameters related to attenuation and sound speed during wave propagation.
The numerical phantoms used in this study are shown in Figure 3. All phantoms were 34 mm wide, 80 mm deep, and 10 mm thick, and the interval between simulated beam lines was set to 0.2 mm. We assumed that the face of the transducer was in direct contact with the phantoms, and set both the axial transmit focus and elevation focus to 80 mm. We assumed a uniform sound speed of 1540 m/s through the entire depth of the phantoms, while distributions of 50-µm polystyrene beads were included to model backscatter. All simulations were performed in 0.1-MHz increments in the 1 to 15 MHz range.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 16 adjusting the number density of polystyrene scatterers in the cylindrical inclusion. Hence, in Figure  3a (p12) the backscatter intensity is 12 dB higher in the cylindrical inclusion than it is in the background medium, while it is 12 dB lower in Figure 3b (m12). In contrast to the backscatter property, the attenuation coefficients of the inclusion and the background were same and set to 5 dB/cm/MHz in both phantoms. The reference phantoms used for diffraction compensation (DC) are shown in Figure 3c,d. With these phantoms, the backscatter and attenuation properties in the entire medium are uniform. The scatterer number density was set as 9.7/mm 3 for each phantom, while the attenuation coefficient was 0.3 dB/cm/MHz for Figure 3c, and 0.7 dB/cm/MHz for Figure 3d.

Comparison of the DFT and DCT Basis Matrices
The DFT and DCT basis matrices for converting the original signal to a sparse form were chosen because they both contain the frequency information that is necessary to make the proposed framework viable. To identify the preferred option, we conducted a comparison of the performances of both matrices, based on the resulting sparsity after each was applied to the original signal. Here, the original signal is a random ultrasound pulse, 6 mm long in the axial direction, that was transformed by the DFT and DCT. The transformation coefficients for both matrices were sorted in descending order and normalized with respect to their maximum. Coefficients larger the 1% of maximum relative weight ( Figure 4) were counted, and sparsity was subsequently compared based on the number of significant coefficients. From Figure 4, it can be observed that significant information was encoded in 128 coefficients with the DCT, compared to the 211 required with the DFT. This indicates that the DCT exhibits better sparsity than the DFT in contrasting inclusion images, suggesting that DCT basis matrix can achieve better spectral estimation performances under the same compression ratio. Note that although the transformation coefficients of the DCT do not strictly match the physical frequency components of the original signal, following averaging, this should not affect the accuracy of most spectral analysis processes.  Figure 3a,b depict the sample phantoms used in the experiments. Each phantom includes a cylindrical region with contrasting backscatter properties to the background medium, generated by adjusting the number density of polystyrene scatterers in the cylindrical inclusion. Hence, in Figure 3a (p12) the backscatter intensity is 12 dB higher in the cylindrical inclusion than it is in the background medium, while it is 12 dB lower in Figure 3b (m12). In contrast to the backscatter property, the attenuation coefficients of the inclusion and the background were same and set to 5 dB/cm/MHz in both phantoms.
The reference phantoms used for diffraction compensation (DC) are shown in Figure 3c,d. With these phantoms, the backscatter and attenuation properties in the entire medium are uniform. The scatterer number density was set as 9.7/mm 3 for each phantom, while the attenuation coefficient was 0.3 dB/cm/MHz for Figure 3c, and 0.7 dB/cm/MHz for Figure 3d.

Comparison of the DFT and DCT Basis Matrices
The DFT and DCT basis matrices for converting the original signal to a sparse form were chosen because they both contain the frequency information that is necessary to make the proposed framework viable. To identify the preferred option, we conducted a comparison of the performances of both matrices, based on the resulting sparsity after each was applied to the original signal. Here, the original signal is a random ultrasound pulse, 6 mm long in the axial direction, that was transformed by the DFT and DCT. The transformation coefficients for both matrices were sorted in descending order and normalized with respect to their maximum. Coefficients larger the 1% of maximum relative weight ( Figure 4) were counted, and sparsity was subsequently compared based on the number of significant coefficients. From Figure 4, it can be observed that significant information was encoded in 128 coefficients with the DCT, compared to the 211 required with the DFT. This indicates that the DCT exhibits better sparsity than the DFT in contrasting inclusion images, suggesting that DCT basis matrix can achieve better spectral estimation performances under the same compression ratio. Note that although the transformation coefficients of the DCT do not strictly match the physical frequency components of the original signal, following averaging, this should not affect the accuracy of most spectral analysis processes.   Figure 5 shows the results of envelope detection of CS ultrasound data using the p12 numerical phantom depicted in Figure 5a, where the backscatter intensity in the cylindrical inclusion is 12 dB higher than the background medium. Here, CS data were randomly sampled at a 50% ratio and reconstructed using the DCT basis matrix. For comparison, original RF data (dotted line) is plotted in the same graph as CS data processed using the traditional method (blue line) and the proposed framework (red line). Although both signal-processing methods produce representative envelopes of the original RF signal, the result obtained with the proposed framework is much smoother than that obtained with the traditional framework, as seen in Figure 5b (the magnified view of the region highlighted in Figure 5a).   Figure 5 shows the results of envelope detection of CS ultrasound data using the p12 numerical phantom depicted in Figure 5a, where the backscatter intensity in the cylindrical inclusion is 12 dB higher than the background medium. Here, CS data were randomly sampled at a 50% ratio and reconstructed using the DCT basis matrix. For comparison, original RF data (dotted line) is plotted in the same graph as CS data processed using the traditional method (blue line) and the proposed framework (red line). Although both signal-processing methods produce representative envelopes of the original RF signal, the result obtained with the proposed framework is much smoother than that obtained with the traditional framework, as seen in Figure 5b (the magnified view of the region highlighted in Figure 5a). reconstructed using the DCT basis matrix. For comparison, original RF data (dotted line) is plotted in the same graph as CS data processed using the traditional method (blue line) and the proposed framework (red line). Although both signal-processing methods produce representative envelopes of the original RF signal, the result obtained with the proposed framework is much smoother than that obtained with the traditional framework, as seen in Figure 5b (the magnified view of the region highlighted in Figure 5a). For further quantitative analysis, quality of B-mode images was compared by the weighted mean and standard deviation for original data and reconstructed CS data with traditional and the For further quantitative analysis, quality of B-mode images was compared by the weighted mean and standard deviation for original data and reconstructed CS data with traditional and the proposed methods [39]. The ROIs were selected inside of the inclusion region of +12dB (p12), 0dB (uniform), and −12dB (m12) contrast TM phantoms of the 70% compression. As shown in Figure 6, the mean pixel values in B-mode images of two reconstructed data using traditional and the proposed methods were close to the original data for all three phantoms. As for the estimation variances, two reconstructed data showed larger than the original, but the proposed method exhibited smaller standard deviations than the traditional method for all three cases as expected in Figure 5b. proposed methods [39]. The ROIs were selected inside of the inclusion region of +12dB (p12), 0dB (uniform), and −12dB (m12) contrast TM phantoms of the 70% compression. As shown in Figure 6, the mean pixel values in B-mode images of two reconstructed data using traditional and the proposed methods were close to the original data for all three phantoms. As for the estimation variances, two reconstructed data showed larger than the original, but the proposed method exhibited smaller standard deviations than the traditional method for all three cases as expected in Figure 5b. We also calculated the root mean square error (RMSE) of envelopes detected using the traditional framework and the proposed framework, with DFT and DCT basis matrices, at compression rates varying from 50% to 90% to present the overall performance. The resulting RMSE values were normalized with respect to the maximum, i.e., the results of envelope detection using DFT at 90% compression. A summary of these results is given in Table 1, which indicates that DCT produces more accurate results than the DFT, particularly at low compression rates, due to its effectiveness in producing a sparse signal. In addition, the proposed framework produces better We also calculated the root mean square error (RMSE) of envelopes detected using the traditional framework and the proposed framework, with DFT and DCT basis matrices, at compression rates varying from 50% to 90% to present the overall performance. The resulting RMSE values were normalized with respect to the maximum, i.e., the results of envelope detection using DFT at 90% compression. A summary of these results is given in Table 1, which indicates that DCT produces more accurate results than the DFT, particularly at low compression rates, due to its effectiveness in producing a sparse signal. In addition, the proposed framework produces better results than the traditional framework for all compression rates. The combination of qualitative and quantitative evidence shown in Figure 5 and Table 1, respectively, demonstrates that the proposed framework can suitably replace existing frameworks for time-domain analysis. We also compared the relative computational speeds of the traditional and proposed frameworks with DCT basis matrix, noting that the proposed framework completed computations in 63.6% of the time required for the traditional framework, based on the average for all compression rates. This shorter computational time was expected, as the lengthy time-domain conversion and Hilbert transform processes used by the traditional framework to obtain the analytic signal were replaced by the IFFT, as shown in Figure 2.

Attenuation Estimation
Since the DCT performs better than the DFT, as demonstrated in the previous two sections, in estimating the attenuation coefficients from CS ultrasound data, we selected a DCT basis matrix. In these simulations, CS data was randomly sampled from original RF data at a 50% ratio. Diffraction compensation (DC) was provided by the uniform attenuating phantoms as references depicted in Figure 3c,d. Figure 7a shows the variation of the estimated spectral centroid with respect to the depth dimension of the ROI, for the p12 numerical phantom. The beam focusing applied alleviates the original signal's decay in the axial direction, as shown by the blue dotted line, making estimation of attenuation properties less accurate. A more representative relationship between the centroid and depth, shown as the red solid line, is obtained after applying DC, enabling accurate estimation of attenuation coefficients, as shown in Figure 7b. In spite of the fact that the ROI in the p12 phantom is brighter than the background in the B-mode image (due to its larger backscatter intensity), at 0.5 dB/cm/MHz, the attenuation coefficients estimated are, on average, the same as the coefficients estimated for the background.
A similar simulation was conducted with the m12 numerical phantom. The results are shown in Figure 8. Similar to the behavior noted for the p12 phantom, the attenuation coefficient estimated for the ROI following the application of DC (Figure 8b) is the same as that of the background, despite its darker appearance. These two cases demonstrate that spectral analysis of CS data can provide important signal characteristics that are unavailable from time-domain representations.
As before, for a quantitative comparison of the performance of the traditional and proposed frameworks in frequency-domain analysis, we calculated the RMSE of the attenuation coefficients estimated for the p12 phantom at compression rates varying from 50% to 90%. All resulting RMSE values were normalized with respect to the largest value, which was again obtained at a 90% compression rate, using the traditional framework. The results of this analysis are summarized in Table 2, which shows that the proposed framework produces much smaller RMSEs than the traditional method for lower compression rates, with the difference decreasing as the compression rate is increased. We posit that this improved performance is due to the ability of the proposed framework to avoid cumulative estimation errors, since its direct use of spectral information from CS data removes redundant transformation steps. As a further consequence, the speed of computation is also improved; the analysis using the proposed framework required only 24.6% of the time taken using the traditional method, based on the average time for all the compression rates.
original signal's decay in the axial direction, as shown by the blue dotted line, making estimation of attenuation properties less accurate. A more representative relationship between the centroid and depth, shown as the red solid line, is obtained after applying DC, enabling accurate estimation of attenuation coefficients, as shown in Figure 7b. In spite of the fact that the ROI in the p12 phantom is brighter than the background in the B-mode image (due to its larger backscatter intensity), at 0.5 dB/cm/MHz, the attenuation coefficients estimated are, on average, the same as the coefficients estimated for the background.
A similar simulation was conducted with the m12 numerical phantom. The results are shown in Figure 8. Similar to the behavior noted for the p12 phantom, the attenuation coefficient estimated for the ROI following the application of DC (Figure 8b) is the same as that of the background, despite its darker appearance. These two cases demonstrate that spectral analysis of CS data can provide important signal characteristics that are unavailable from time-domain representations.    As before, for a quantitative comparison of the performance of the traditional and proposed frameworks in frequency-domain analysis, we calculated the RMSE of the attenuation coefficients estimated for the p12 phantom at compression rates varying from 50% to 90%. All resulting RMSE values were normalized with respect to the largest value, which was again obtained at a 90% compression rate, using the traditional framework. The results of this analysis are summarized in Table 2, which shows that the proposed framework produces much smaller RMSEs than the traditional method for lower compression rates, with the difference decreasing as the compression rate is increased. We posit that this improved performance is due to the ability of the proposed framework to avoid cumulative estimation errors, since its direct use of spectral information from CS data removes redundant transformation steps. As a further consequence, the speed of computation is also improved; the analysis using the proposed framework required only 24.6% of the time taken using the traditional method, based on the average time for all the compression rates.

Discussion
In this paper, we proposed a complete signal-processing framework for temporal and spectral analysis of ultrasound CS data. With our method, relevant spectral information is directly extracted from a compressed measurement. This framework is designed for ultrasound CS data that avoids a redundant reconstruction process. This is enabled by selecting an appropriate basis matrix for use in the data acquisition and reconstruction processes. Hence, we compared the performances of two fundamental basis matrices, determining that the DCT matrix offered the best performance for time and frequency domain analysis.
From our simulations, the proposed framework demonstrated good performance in both temporal and spectral analysis. Based on this, we introduce a novel imaging mode for medical ultrasound applications combining the results of time and frequency domains analysis, since each domain provides information that is inaccessible from the other. We term the resulting combination of a B-mode image and an attenuation coefficient map, shown in Figure 9, a QUS image (Q-mode image). As proposed, the Q-mode image offers a comprehensive representation of the quantitative measurements from frequency-domain medical ultrasound, and the spatial (or anatomical for in vivo case) relationship from the surrounding B-mode image. This thus makes it possible to conclude professional pathological judgments that cannot be provided from exclusive time-domain analysis (i.e., analysis of traditional B-mode images). With the careful selection of ROI, in addition, the combination of spatial information in the time domain and quantitative results in the frequency domain would aid the development of various deep learning applications for medical ultrasound data. As the deep learning application with a combined MRI images exhibited good performance [40], the new category of comprehensive dataset such as the Q-mode image, also makes it possible to obtain promising results for recognition of pathological characteristics which is not shown in a B-mode image alone.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 16 case) relationship from the surrounding B-mode image. This thus makes it possible to conclude professional pathological judgments that cannot be provided from exclusive time-domain analysis (i.e., analysis of traditional B-mode images). With the careful selection of ROI, in addition, the combination of spatial information in the time domain and quantitative results in the frequency domain would aid the development of various deep learning applications for medical ultrasound data. As the deep learning application with a combined MRI images exhibited good performance [40], the new category of comprehensive dataset such as the Q-mode image, also makes it possible to obtain promising results for recognition of pathological characteristics which is not shown in a Bmode image alone. The concept of the Q-mode image is also flexible, as any combination of results of quantitative analysis can be included with the B-mode image, which is recognized as the most important representation of time-domain information in medical ultrasound applications. Although our demonstration features the attenuation coefficient map as the result of frequency-domain analysis, QUS information such as backscatterer properties, speed of sound, or elasticity, resulting from frequency-domain analysis could be combined with a B-mode image for a new kind of Q-mode image.

Conclusions
The CS technique has been widely adopted in signal-processing applications due to the faster data acquisition times and smaller storage spaces it offers. However, its adoption in some applications has been hindered by the fact that compressed data requires reconstruction before The concept of the Q-mode image is also flexible, as any combination of results of quantitative analysis can be included with the B-mode image, which is recognized as the most important representation of time-domain information in medical ultrasound applications. Although our demonstration features the attenuation coefficient map as the result of frequency-domain analysis, QUS information such as backscatterer properties, speed of sound, or elasticity, resulting from frequency-domain analysis could be combined with a B-mode image for a new kind of Q-mode image.

Conclusions
The CS technique has been widely adopted in signal-processing applications due to the faster data acquisition times and smaller storage spaces it offers. However, its adoption in some applications has been hindered by the fact that compressed data requires reconstruction before further signal-processing for analysis in either the time or the frequency domains, a process that features redundant transformations between the two domains. In this paper, we proposed a simpler framework that avoids full reconstruction of ultrasound CS data. Here, using an appropriate basis matrix to transform the signal to its sparse form, spectral properties can be directly extracted from compressed data, and the conversion to a complex analytic signal for detection of the envelope of the original measurement in the time domain is made more efficient. Simulations conducted with numerical TM phantoms confirmed that this framework is simpler and faster in execution times than the traditional framework, average 63.6% of the time required for the envelope detection and average 24.6% for the attenuation estimation at compression rates varying from 50% to 90%. For the estimation accuracy, the RMSEs of envelope detection and attenuation estimation were averagely reduced by 10% and 38% compared to traditional method at all compression rates. Hence, we suggest that the proposed framework offers a new standard for signal processing using ultrasound CS data, especially for small and portable systems, and applications that handle large datasets.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: