Compressed-sensing photoacoustic computed tomography in vivo with partially known support

: Compressed sensing (CS) can recover sparse signals from under-sampled measurements. In this work, we have developed an advanced CS framework for photoacoustic computed tomography (PACT). During the reconstruction, a small part of the nonzero signals’ locations in the transformed sparse domain is used as partially known support (PKS). PACT reconstructions have been performed with under-sampled in vivo image data of human hands and a rat. Compared to PACT with basic CS, PACT with CS-PKS can recover signals using fewer ultrasonic transducer elements and can improve convergence speed, which may ultimately enable high-speed, low-cost PACT for various biomedical applications.


Introduction
Photoacoustic imaging is a novel noninvasive imaging modality combining the advantages of both optical imaging and ultrasonic imaging [1].Compared with ultrasonic imaging, it provides high optical contrast sensitive to important physiological parameters (such as the oxygen saturation of hemoglobin and the metabolic rate of oxygen); compared with purely optical imaging, it offers better spatial resolution for deep imaging.Although young, It has attracted much attention with its wide applications in various preclinical imaging settings [2,3].Moreover, clinical photoacoustic imaging technologies are also emerging [4].By now, it has been explored for the early diagnosis of cancer, the imaging of tumor angiogenesis, and identification of cardiovascular vulnerable plaques, in both preclinical and clinical studies [4][5][6][7].Photoacoustic computed tomography (PACT) is one form of photoacoustic imaging [1], based on the reconstruction of the internal photoacoustic source distributions from the ultrasonic measurements over a plane or a volume.Generally, it uses either an unfocused ultrasonic transducer with mechanical scanning or an ultrasonic transducer array to acquire the image data [8][9][10]; then, the optical properties of the imaged objects can be recovered using a reconstruction algorithm.
For many biomedical applications, PACT with a high-frequency ultrasonic array is of particular interest, owing to its capability of both high-speed and high-resolution imaging [11].However, the cost of a high-frequency ultrasonic array can be very high.In addition, due to challenges in fabrication, a fair amount of today's high-frequency ultrasonic arrays are not dense enough to spatially fully sample the signals, which requires a pitch that is half of the center ultrasonic wavelength or smaller.For example, in one high-frequency linear-array PACT system, the pitch of the ultrasonic array is as large as twice of the center ultrasonic wavelength [12], resulting in reconstruction artifacts.Generally, to achieve high-quality PACT with conventional back-projection reconstruction, a number of densely packed ultrasonic transducers (or array elements) are needed, which can give rise to a high cost, and a long data acquisition and processing time.Compressed sensing (CS) is a novel information theory capable of recovering signals that are sparse from under-sampled measurements.To date, CS has been widely studied and used in MRI [13], diffuse optical tomography [14] and CT [15].Recently, explorations aiming to leverage the advantages of CS for PACT have also been reported.For example, J. Provost et al. have conducted numerical and phantom studies on CS-based PACT [16]; D. Liang et al. have numerically demonstrated CS-based photoacoustic imaging under random optical illumination [17]; Z. Guo et al. have reported CS-based photoacoustic reconstruction in time domain with images acquired in vivo [18]; M. Sun et al. have developed an arc-direction compressed-sensing PACT algorithm with numerical phantoms [19].All these studies have shown that CS-based reconstruction techniques can recover photoacoustic signals sampled at a (temporal and/or spatial) rate lower than that required by the Nyquist sampling theory.Therefore, both the number of ultrasonic transducers required for the PACT system and the amount of acquired data needed for reconstruction can be significantly reduced.
However, although CS has been used for PACT, this exploration is in its infancy.To our knowledge, the use of CS-based PACT in frequency domain has not been studied and validated with in vivo image data yet.More importantly, many advanced CS theories that may of great potential benefits for PACT still remain unexplored.For instance, recently, N. Vaswani et al. developed a modified-CS method, in which partially known support (PKS), i.e., some known nonzero locations in the transformed sparse domain, was used as an additional prior besides sparsity [20].Theoretical studies and applications in MRI have shown that PKS can further reduce the number of measurements needed for high-fidelity reconstruction [21], a significant benefit of great interest to PACT as well-including the potential to reduce the number of ultrasonic transducers (or array elements) and the image data acquisition time.
In this work, first, using in vivo photoacoustic image data acquired with a linear ultrasonic array, we have validated our basic CS reconstruction methods for PACT in both time and frequency domains.The CS-based PACT has produced significantly better images than those obtained with the conventional back-projection (BP) PACT.In addition, compared to CSbased PACT in time domain, its counterpart in frequency domain has produced images of a better contrast-to-noise ratio (CNR), with a significantly reduced computation time.Most importantly, for the first time to our knowledge, we have developed an advanced CS reconstruction model that incorporated PKS for PACT reconstruction.To solve the complex inverse problem, an iteratively reweighed conjugate gradient descent method (IR-CGD) has been developed.Results with in vivo image data have shown that, even in the case of heavily under-sampling, the developed CS-PKS strategy can recover the photoacoustic signals accurately, suppress the under-sampling artifacts effectively, and converge rapidly with a relatively small number of iterations, all of which are not attainable with previously developed basic CS approaches.These advantages of CS-PKS may ultimately enable highspeed, low-cost PACT to be constructed with less densely packed ultrasonic arrays or less number of ultrasonic transducers for various biomedical applications.

Photoacoustic computed tomography
According to the photoacoustic generation theory, the acoustic pressure ( , ) p t r at position r and time t satisfies the wave equation [22] where ' r denotes the position of the initial acoustic pressure.The above time dependent wave equation can be transformed into the temporal-frequency domain.Denoting the Fourier transform of p as p , we have [23] ( , ) ( , ) exp( ) , where , ω is the angular frequency equal to 2 f π , and f is the signal temporal frequency.Thus, according to Eq. ( 2) and Eq. ( 3), the forward problem in the temporalfrequency domain is expressed as

Compressed sensing
One important prerequisite of compressed sensing is that either the signal or its transform in a certain domain is sparse or compressible.A signal is called s-sparse when s of its coefficients are nonzero (usually s<<N).A signal is s-compressible if its ordered set of coefficients decays rapidly and x can be well approximated by the first s large coefficients.Fortunately, most medical images are sparse in a certain domain by finding an appropriate sparse transform : = ψ x ψθ , where θ is the original image, and x is the transformed one.It has been proven that photoacoustic images can be transformed into sparse domains by a variety of transforms, such as the numerical derivative (ND), the wavelet transform, etc [16][17][18][19].
In PACT, the goal of reconstruction is to recover the photoacoustic image signal θ through the measurement data y from ultrasonic transducers.Given that y is obtained with a measurement matrix K , we have = K y θ .Photoacoustic images with CS can usually be reconstructed by solving the following constrained optimization problem for x : To implement CS reconstruction for PACT, the derived discrete expressions of K in time and temporal-frequency domains for our PACT system are explicitly given below.According to Eq. ( 2) and the back projection principle, the measurement matrix in time domain can be designed as [18] , ( From Eq. ( 4), the measurement matrix in temporal-frequency domain can be written as [16] ( ) , .
In Eq. ( 6) and Eq. ( 7), , i j r represents the Cartesian coordinates of image pixels, m r represents the positions of transducers, p represents the number of transducers, s q and n q represent the number of sampling points in time and frequency domains, respectively.

Compressed sensing with partially known support
Recently, a new type of CS reconstruction framework-CS with partially known support (CS-PKS)-has been developed to incorporate certain prior support of a sparse signal into the CS reconstruction process, where the support is defined as the locations of non-zero signals in the sparse domain [20].Compared with basic CS, CS-PKS has been demonstrated in MRI to recover signals with fewer measurements or obtain higher quality images with the same number of measurements [21,24].
be the sparse coefficients of the desired image and denote the support 0 ( ) is the known part of T with a size of 0 T and is the unknown part of T with a size of ∆ .With the partially known support 0 T , the candidates for the s-sparse signal set x are restricted in a signal space smaller than that in basic CS, and thus requires fewer samples to yield an accurate reconstruction [25].
Since the signal x is known to be nonzero at some locations, CS-PKS allows us to minimize the number of nonzeros out of the known part during searching for a sparse solution to = Φx y .Using convex relaxation, this procedure can be formulated as [21] 1 2 min . ., where ∆ x denotes the signals outside the known support, ε the noise level, and

Reconstruction model
Based on the above theory of CS and PKS, the reconstruction models of basic CS and CS-PKS were developed for PACT, which are expressed by the Eq. ( 9) and Eq. ( 10), In Eq. ( 9) and Eq. ( 10), F is the objective function consisting of three parts-the first part represents the square error between the estimated measurements from the reconstructed signal and the experimentally acquired measurements, the second part represents the 1 l norm of signal in a sparse domain, and the third part represents the total variation (TV) penalty of the signal, and α and β are the regularization parameters determining the trade-off between the data consistency and the sparsity.The difference between the CS reconstruction model and the CS-PKS model is that the PKS prior information is incorporated into the second part of F in Eq. (10).Please note that it is often useful to include a TV item in CS even when other sparse transforms are used in the objective function [13,17,26].The definition of ∆ x in Eq. ( 10) is given by ( ) To solve Eq. ( 10), an iterative reweighed method is employed, and Eq. ( 10) is rewritten as where i represents the th i iteration, W is a diagonal matrix whose diagonal element is set to 1 when the corresponding signal element belongs to the unknown support ∆ in the ( 1)  T , which can be obtained by thresholding the reconstructed sparse coefficients [21] { } where ( ) i z x is the th z element of reconstructed signal x in the th i iteration and ( ) i τ is a threshold used to determine the set of partially known support ( ) 0 i T .It has been found that, the iterative reconstruction process can work well for fast-decaying signal ( ) , where δ is a hyper-parameter.

Optimization algorithm
In order to solve Eq. ( 12), an iteratively reweighed conjugate gradient descent method (IR-CGD) was developed.The result from the first few basic CS iterations, denoted as (0 ) x , is used as the initial value for subsequent CS-PKS reconstruction.The major steps of the algorithm are specified below.
Step 2: Do the inner loop 2.1 Initialize the inner loop variable x represents the reconstructed signal of the th l inner loop in the th i outer loop).

Compute the gradient F
∇ of the objective function F with respect to signal x in Eq. ( 11).

Obtain ( ) 1 i l +
x with the CGD method.
, go to step3; otherwise, let 1 l l = + , and go to step 2.2.

In vivo imaging
Noninvasive photoacoustic imaging of a rat and human hands was performed in vivo using a linear-array PACT system.The major components of the system include: (1) a tunable dye laser (Cobra, Sirah Laser-und Plasmatechnik GmbH, Germany) pumped by a Q-switched Nd:YLF laser (INNOSLAB, Edgewave GmbH, Germany), (2) a custom-built linear ultrasonic array with a center frequency of 30 MHZ, and (3) an 8-core PC (Dell Precision 490) equipped with an 8-channel high-speed PCI data acquisition (DAQ) card (Octopus CompuScope 8389, GaGe, USA).Specifically, the high-frequency ultrasonic array was fabricated from a 2-2-piezo-composite by the NIH Resource Center for Medical Ultrasonic Transducer Technology at the University of Southern California.The array had 48 elements (82 m µ × 2 mm) with 100 m µ spacing.The elements were elevationally focused at 8.2 mm, providing an elevational resolution of 200 -300 m µ within a ~3.5-mm focal zone.The mean fractional bandwidth was ~70% for receiving-only operation, as used in our PACT system.The axial and lateral resolutions of the system at 8 mm normal depth from the transducer surface were ~25 and ~80 m µ , respectively.Further details of the imaging system can be found in our previous publications [10,12].Note that, (1) the name ""microscopy" was used in references [10,12] only because high microscopic spatial resolution was achieved, the imaging system (which was also used for this study) is essentially a linear-detection geometry (2) because a linear array was used, our imaging system is considered as a limitedview PACT; (3) as the pitch of the ultrasonic array was 100 m µ , about twice of the center ultrasonic wavelength, the array is considered as not spatially fully sampled [27].
For in vivo imaging experiments, the dye laser output was tuned to 584 nm-an isosbestic point at which the oxy-and deoxy-hemoglobin absorb equally.For small animal imaging, a Sprague Dawley rat (Harlan Laboratories, Inc., USA) was used.Intra-dermal injection of a mixture of ketamine (85 mg/kg) and xylazine (15 mg/kg) was implemented for anesthesia.Before photoacoustic imaging, the hair in the imaged region was removed with commercial depilatory lotion.For human hand imaging, the optical fluence on the skin surface was set to ~0.5 mJ/cm 2 per pulse, well below the ANSI [28] recommended Maximum Permissible Exposure (MPE) of 20 mJ/cm 2 for a single pulse in the visible spectral range.The timeaveraged light intensity during image acquisition was ~150 mW/cm 2 , also below the ANSI recommended MPE of 196 mW/cm 2 , calculated by 1.1t 0.25 W/cm 2 .As the ANSI safety limit for this pulse width region is dominantly based on the thermal mechanism, our compliance with the ANSI standards guarantees no thermal damage to the tissue.All animal and human experiments described here were carried out in compliance with Washington University approved protocols.

Results in time and frequency domains
In this section, Eq. ( 9) was adopted as the basic CS reconstruction model, and a four-level Daubechies wavelet was used as the sparse transform.A photograph of the imaged hand (human hand-1) is shown in Fig. 1(A), in which a dashed box indicates the imaged area (of ~6.4 mm × 9.6 mm).For each 3D image, 160 B-scan frames were acquired, with each reconstructed B-scan image consisting of 128 × 128 pixels (corresponding to a cross-section of ~6.4 mm × 1.6 mm).Note that, in all B-scans, although Hilbert transform was used for envelop-detection to display the images, it was applied only after reconstructions.The rest of Fig. 1 shows the reconstructed images with data from all 48 ultrasonic transducer elements.Figure 1(B) is the maximum amplitude projection image (MAP) -the maximum photoacoustic amplitudes projected along the depth direction to the surface of the handreconstructed by the conventional photoacoustic BP method.Figure 1 With data from all 48 transducer elements, it can be seen that, although the linear-array is not considered as spatially fully sampled, the photoacoustic images can be reconstructed well using any of the three reconstruction methods (note that the imaging system using BP for reconstruction-with data from all 48 transducer elements-has already been validated with optical microscopy of excised tissue in small animals [12]).In addition, compared with the TD-CS reconstruction, the FD-CS reconstruction has produced images of a slightly higher contrast-to-noise ratio (CNR), yet has consumed ~50% less computation time (for each B-scan reconstruction, it takes about 40 s using FD-CS vs. 75 s using TD-CS on a PC with a Intel Core i5-2400 CPU of 3.1 GHz).still retain the major features seen in the images reconstructed with data from all 48 transducer elements (Fig. 1).Similarly, compared with TD-CS, the FD-CS approach has shown a higher CNR with ~50% less reconstruction time.
The results suggest that, in the case of under-sampling, the FD-CS reconstruction can produce images of a higher CNR, compared to those obtained with BP and TD-CS.However, when using data from only 16 transducer elements (i.e., heavily undersampling), significant under-sampling artifacts emerge in all reconstructed images (Fig. 2).

Results with CS-PKS method
To further improve the image quality, PKS was incorporated into our FD-CS photoacoustic reconstruction process (section 3).Specifically, the result after a few iterations of traditional CS was used as the initial image.To start the iterative reconstruction process, it was transformed into a sparse domain using a four-level Daubechies wavelet.According to Eq. ( 13), in each iteration, the locations (in the sparse domain) of which the coefficients are larger than a threshold were used as the partially known support.In vivo imaging of a human hand (human hand-2) and a rat were performed to validate and demonstrate the advantages of the developed CS-PKS strategy.All CS-PKS reconstructions were performed based on the optimization model defined by Eq. ( 12).The imaged volume is about 6.4 mm × 10 mm × 1.6 mm (in x, y, and z directions, respectively).For each of the following experiments, 166 Bscan frames were acquired, with each reconstructed B-scan image consisting of 128 × 128 pixels.
Figure 3 shows the reconstructed photoacoustic images of human hand-2.Figures 3(A) and (B) are the MAP images reconstructed by BP with data from 48 and 16 transducer elements, respectively.Figures 3(C) and (D) are the MAP images reconstructed by the basic FD-CS and the FD-CS-PKS methods, respectively, with data from only 16 transducer elements.Figures 3(a) -(h) are the corresponding B-scan images along the vertical dashed lines shown in Figs.3(A) -(D).Compared with the reconstructed images using data from all 48 transducer elements (Figs.3(A almost to a degree obtainable only with data from all 48 transducer elements using the conventional BP reconstruction.In addition, our experiments have also shown that the CS-PKS method is able to improve the convergence speed.For each B-scan reconstruction, it takes 7 s using FD-CS-PKS vs. 12 s using FD-CS, on a PC with an Intel Core i5-2400 CPU of 3.1 GHz.x is the (lateral) direction of the transducer array, y is the mechanical scanning direction, and z is the depth direction.MAP, maximum amplitude projection; BP, back-projection; FD-CS, compressed sensing in frequency domain; FD-CS-PKS, compressed sensing with partially known support in frequency domain.The color-scale represents relative optical absorption.
In vivo imaging of a rat was also performed to further demonstrate the advantages of the developed FD-CS-PKS strategy.Figure 4 shows the reconstructed images of the subcutaneous vasculature of the back of the rat.Figures 4(A) and (B) are the MAP images reconstructed using BP with data from 48 and 16 transducer elements, respectively.Figures 4(C) and (D) are the MAP images reconstructed using the basic CS and the FD-CS-PKS methods, respectively, with data from 16 transducer elements.Figures 4(a) -(h) are the corresponding B-scan images along the dashed lines shown in Figs.4(A) -(D).In this case, similar to the results obtained with human hand-2, BP recovers only a small part of useful features.Furthermore, while FD-CS clearly improves the image quality, FD-CS-PKS produces the best images of a high CNR and significantly reduced artifacts.x is the (lateral) direction of the transducer array, y is the mechanical scanning direction, and z is the depth direction.MAP, maximum amplitude projection; BP, back-projection; FD-CS, compressed sensing in frequency domain; FD-CS-PKS, compressed sensing with partially known support in frequency domain.The color-scale represents relative optical absorption.

Quantitative comparisons
From the above results, the advantages of our CS-PKS reconstruction strategy can be seen visually.In this section, to further compare the results, three quantitative analyses are performed (in which all images were normalized by dividing their corresponding maximum values).
The first quantitative approach for comparison is to compute the histograms of the difference images-images obtained by subtracting the control from the reconstructed images and then taking the absolute value (Fig. 5).In our calculation, the FD-CS image reconstructed with data from all 48 transducer elements was used as the control (using the image reconstructed by BP with data from all 48 transducer elements as the control leads to the same conclusions).Figures 5(A) -(C) are the difference histograms between the control and the images reconstructed by BP, TD-CS, and FD-CS, respectively, with data from 16 transducer elements.From Fig. 5(A) to Fig. 5(C), the number of counts lying in the lower ranges of the difference amplitudes increases significantly, suggesting that the images reconstructed using the CS-based methods (in particular, FD-CS) possessing higher fidelity.From Figs. 5(D) -(F), it is obvious that, for the images reconstructed with CS, the majority of the difference amplitudes lies in the lower ranges.In particular, in the case of FD-CS-PKS, the number of counts in [0, 0.1] represents >90% of all pixel numbers.Figures 5(G) -(I), obtained with the rat model, show very similar results.Further, localized comparisons were performed by plotting the photoacoustic amplitudes along chosen lines in the MAP images (Fig. 6).Figures 6(A) and (B) are the plots for human hand-2 and the rat along the pink lines shown in Fig. 3(A) and Fig. 4(A), respectively.In order to have a better quantitative comparison, the CNRs of selected signal peaks were also computed (Fig. 6, insets) (in our computation, the range from 0.0 mm to 1.5 mm along the xaxis is used as an estimation of the background).It can be seen that, (1) images from CSbased reconstructions generally have higher CNRs compared to those by the BP method, and (2) the images reconstructed by CS-PKS with only 16 transducers have shown equivalent CNRs to that of BP images with data from all 48 ultrasonic elements.Fig. 6.Photoacoustic amplitudes (relative optical absorption) along the chosen pink lines in MAP images of (A) human hand-2 (Fig. 3 (A)), and (B) the rat (Fig. 4 (A)).Insets, the contrast to noise ratios (CNRs) of selected signal peaks.BP, back-projection; CS, compressed sensing; PKS, compressed sensing with partially known support; the number 48 or 16 after each reconstruction method name indicates that the reconstruction is performed with data from 48 or 16 transducer elements.
The third quantitative approach for comparison is based on the computation of the crosscorrelation coefficient (corr), as defined by Eq. ( 14): where o I is the control image, r I is the reconstructed image, i and j represent the row and column of the images.r I and o I are the mean values of the reconstructed images and the control, respectively.For all models ("Hand-1", "Hand-2", and "Rat"), using the images reconstructed by FD-CS48 as the control, the computed cross-correlation coefficients are shown in Table 1.It can be seen clearly that the images reconstructed using CS have larger cross-correlation coefficients than those reconstructed by BP.Furthermore, the crosscorrelation coefficients of the images reconstructed using our CS-PKS strategy are among the largest of all approaches.

Discussion and conclusions
The above results have clearly demonstrated the superiority of our advanced CS-PKS strategy over the basic CS approaches.Here, we discuss a few steps during the reconstruction process that are critical for achieving optimal results.First, in the optimization process with the IR-CGD method, two regularization parameters α and β -representing the weight of PKS and the ND regularization item, respectively-need to be determined appropriately.Overweighed values will result in the distortion of the reconstruction.On the contrary, if their proportions to the objective function are too small, their support would become ineffective.What's more, the weight distribution between PKS and ND, i.e., the relative magnitude of α and β , also matters.Unfortunately, at present, there is no general method to effectively compute the optimal values of these regularization parameters, which vary with different reconstruction algorithms and raw data sets.In our experiments, they are empirically determined by trying different combinations and choosing the best one.Second, the parameter δ , which is used to determine the set of the detection support, is another key factor affecting the final reconstruction results.According to the relationship between δ and τ (section 3.1), a large δ will result in too many false locations in the detection support, and a small δ will result in too few locations in the detection support.Finally, several other factors may also influence the reconstructed image quality, including the step for updating signals in each iteration, the sparsity of the signals, and the choice of the sparse transform basis.The best values of the reconstruction parameters used in our experiments are listed in Table 2.In summary, we have developed CS-based reconstruction methods for PACT in time and frequency domains.Supported by in vivo image data, we have found that, in the case of under-sampling, the CS-based photoacoustic reconstructions are superior to the conventional photoacoustic BP reconstruction.Furthermore, the CS-PKS reconstruction framework developed by us has demonstrated accurate reconstruction results even with heavily undersampled data from only 16 ultrasonic transducer elements.Compared to PACT with basic CS, PACT with the developed CS-PKS strategy has demonstrated significant advantages, including a higher CNR, a faster convergence, and a more effective suppression of undersampling artifacts.The promising results demonstrated in this work suggest that the development of novel advanced CS reconstruction strategies has the potential to enable highquality PACT with a significantly reduced number of ultrasonic transducers (or array elements).In return, this will reduce the cost and the amount of image data required for accurate reconstruction, which may ultimately have a profound impact in promoting various preclinical and clinical applications of PACT.

#
167093 -$15.00USD Received 19 Apr 2012; revised 22 May 2012; accepted 22 May 2012; published 6 Jul 2012 (C) 2012 OSA PACT; (C) is the MAP image reconstructed by CS in time domain (TD), and Fig. 1(D) is the MAP image reconstructed by CS in frequency domain (FD), in which a band-pass filter of 0 -40 MHz is applied.Figures 1(a) -(f) are the B-scan images along the positions represented by the vertical dash lines in Figs.1(B) -(D).

Fig. 1 .
Fig. 1.Reconstructed images of the subcutaneous vasculature of human hand-1, with data from all 48 transducer elements.(A) Photograph of the hand; (B) -(D) MAP images reconstructed by the BP, TD-CS, and FD-CS methods, respectively; (a) -(f) B-scan images along the dashed lines in (B) -(D).x is the (lateral) direction of the transducer array, y is the mechanical scanning direction, and z is the depth direction.MAP, maximum amplitude projection; BP, back-projection; TD-CS, compressed sensing in time domain; FD-CS, compressed sensing in frequency domain.The color-scale represents relative optical absorption.

Figure 2
Figure 2 shows the reconstructed images in the case of more heavily undersamplingwith data from only 16 uniformly distributed transducer elements of the array.Similar to the approaches above, while all 128 time points per transducer element are used for BP and TD-CS reconstructions, only signals representing frequencies between 0 -40 MHz (corresponding to 40 points in frequency domain) are used for the FD-CS reconstruction.In this case, the images reconstructed by BP (Figs. 2(A), (a) and (b)) are becoming significantly worse, while the images reconstructed by TD-CS and FD-CS (Figs. 2(B) -(C) and (c) -(f))still retain the major features seen in the images reconstructed with data from all 48 transducer elements (Fig.1).Similarly, compared with TD-CS, the FD-CS approach has shown a higher CNR with ~50% less reconstruction time.The results suggest that, in the case of under-sampling, the FD-CS reconstruction can produce images of a higher CNR, compared to those obtained with BP and TD-CS.However, when using data from only 16 transducer elements (i.e., heavily undersampling), significant under-sampling artifacts emerge in all reconstructed images (Fig.2).

Fig. 2 .
Fig. 2. Reconstructed images of the subcutaneous vasculature of human hand-1, with data from 16 transducer elements.(A) -(C) MAP images reconstructed by BP, TD-CS, and FD-CS respectively; (a) -(f) B-scan images along the dashed lines in (A) -(C).x is the (lateral) direction of the transducer array, y is the mechanical scanning direction, and z is the depth direction.MAP, maximum amplitude projection; BP, back-projection; TD-CS, compressed sensing in time domain; FD-CS, compressed sensing in frequency domain.The color-scale represents relative optical absorption.
), (a), and (b)), the images reconstructed by BP with data from only 16 transducer elements are much worse (Figs.3(B), (c) and (d))-the features are almost submerged in noise.With the use of regular CS, the image quality improves significantly (Figs.3(C)), (e) and (f)), although some under-sampling artifacts are still notable.Surprisingly, when the PKS prior information was exploited, the under-sampling artifacts are greatly suppressed (Figs.3(D), (g) and (h)), and the image can be recovered #167093 -$15.00USD Received 19 Apr 2012; revised 22 May 2012; accepted 22 May 2012; published 6 Jul 2012 (C) 2012 OSA

Fig. 3 .
Fig. 3. Reconstructed images of the subcutaneous vasculature of human hand-2.(A,B) MAP images reconstructed by BP with data from 48 and 16 transducer elements, respectively; (C,D) MAP images reconstructed by the FD-CS and the FD-CS-PKS methods, respectively, with data from 16 transducer elements; (a) -(h) B-scan images along the dashed lines in (A) -(D).x is the (lateral) direction of the transducer array, y is the mechanical scanning direction, and z is the depth direction.MAP, maximum amplitude projection; BP, back-projection; FD-CS, compressed sensing in frequency domain; FD-CS-PKS, compressed sensing with partially known support in frequency domain.The color-scale represents relative optical absorption.

Fig. 4 .
Fig. 4. Reconstructed images of the subcutaneous vasculature of the back of a rat.(A,B) MAP images reconstructed by BP with data from 48 and 16 transducer elements, respectively; (C,D) MAP images reconstructed by the FD-CS and the FD-CS-PKS methods, respectively, with data from 16 transducer elements; (a) -(h) B-scan images along the dashed lines in (A) -(D).x is the (lateral) direction of the transducer array, y is the mechanical scanning direction, and z is the depth direction.MAP, maximum amplitude projection; BP, back-projection; FD-CS, compressed sensing in frequency domain; FD-CS-PKS, compressed sensing with partially known support in frequency domain.The color-scale represents relative optical absorption.

Fig. 5 .
Fig. 5. Histograms of the amplitudes of the difference images using the image reconstructed by FD-CS with data from all 48 transducer elements as the control.(A) -(C) Histograms of the amplitudes of the difference images for Hand-1; (D) -(F) Histograms of the amplitudes of the difference images for Hand-2; (G) -(I) Histograms of the amplitudes of the difference images for the Rat.BP, back-projection; CS, compressed sensing; TD-CS, compressed sensing in time domain; FD-CS, compressed sensing in frequency domain; FD-CS-PKS, compressed sensing with partially known support in frequency domain; the number 16 after each reconstruction method name indicates that the reconstruction is performed with data from only 16 transducer elements.