Research Papers: Imaging

Comparison of photoacoustic image reconstruction algorithms using the channelized Hotelling observer

[+] Author Affiliations
Adam Petschke, Patrick J. La Rivière

The University of Chicago, Department of Radiology, 5841 South Maryland Avenue, Chicago, Illinois 60637

J. Biomed. Opt. 18(2), 026009 (Feb 05, 2013). doi:10.1117/1.JBO.18.2.026009
History: Received July 24, 2012; Revised January 14, 2013; Accepted January 16, 2013
Text Size: A A A

Open Access Open Access

Abstract.  We demonstrate the use of task-based image-quality metrics to compare various photoacoustic image-reconstruction algorithms, including a method based on the pseudoinverse of the system matrix, simple backprojection, filtered backprojection, and a method based on the Fourier transform. We use a three-dimensional forward model with a linear transducer array to simulate a photoacoustic imaging system. The reconstructed images correspond with two-dimensional slices of the object and are 128×128pixels. In order to compare the algorithms, we use channelized Hotelling observers that predict the detection ability of human observers. We use two sets of channels: constant Q and difference of Gaussian spatial frequency channels. We look at three tasks, identification of a point source in a uniform background, identification of a 0.5-mm cube in a uniform background, and identification of a point source in a lumpy background. For the lumpy background task, which is the most realistic of the tasks, the method based on the pseudoinverse performs best according to both sets of channels.

Figures in this Article

Photoacoustic imaging is based on the detection of acoustic waves that are generated by thermoelastic expansion induced by optical absorption of laser light. The majority of photoacoustic imaging systems can be classified as either microscopy systems or computed tomography (CT) systems.1 In general, the microscopy-based systems do not require a reconstruction step as the lateral resolution is set by either a focused transducer or focused laser illumination. CT systems can be further subdivided into scanning systems (also called synthetic aperture) and array systems. From a mathematical point of view, synthetic apertures and arrays require the same reconstruction algorithm, although arrays tend to have smaller individual transducers so the directivity can be more uniform over the imaging area. Transducer arrays have a speed advantage as they do not need to be mechanically scanned, which can be important in many applications. In addition, the array-based imaging systems have the potential to be integrated with ultrasound imaging systems, and indeed many of the arrays in use in photoacoustic systems are taken from ultrasound systems.24 In this paper we consider a linear transducer array that is similar to the ones used in current photoacoustic and ultrasound imaging systems.

There are several reconstruction algorithms used in linear array imaging systems. Filtered backprojection is based on an analytical inversion of the solution to the photoacoustic wave equation.5 It is theoretically exact for a continuous, infinite linear transducer array, and when the object to be imaged is acoustically homogeneous and has no acoustic absorption. Unfiltered (or simple) backprojection, although not theoretically exact, is numerically equivalent to the delay-and-sum beamformer method employed in ultrasound imaging systems, and hence is very readily implemented in combined ultrasound-photoacoustic imaging systems. The Fourier method of image reconstruction is also based on inverting the solution to the photoacoustic wave equation, but this inversion is performed in the frequency domain.6 For sufficiently small datasets or given sufficiently powerful computer hardware, all of these simple reconstruction algorithms can be used in real-time imaging systems. There are more sophisticated reconstruction algorithms that can correct for acoustic inhomogeneity and absorption such as iterative methods.7 In general, these iterative algorithms produce high-quality images; unfortunately, using current computers they are not capable of real-time imaging.

The minimum-norm least-squares (MNLS) solution to reconstructing a photoacoustic image is given by the pseudoinverse of the system matrix, which relates the initial pressure distribution to the measured pressure wave. In this paper, we calculate the pseudoinverse using a singular-value decomposition (SVD). The pseudoinverse calculation involves dividing by the nonzero singular values of the system matrix. In practice, we have to regularize the problem by throwing out small singular values that would otherwise amplify the noise. Where to set the regularization threshold depends on the specific signal and noise levels. The noise in an imaging system can generally be divided into object-dependent and object-independent noise. In many imaging systems, such as x-ray CT, the object-dependent noise varies strongly across detectors. In photoacoustics, however, the object-dependent noise is relatively uniform across transducer elements (as long as the object is not too close to the transducer). The object-independent noise (such as electronic noise) is also uniform across the transducer elements. The regularization level still depends on the relative strength of the signal and noise, but it is possible to calculate a few pseudoinverse matrices based on different singular-value truncations, store these in memory, and then decide in real time which one to use based on image quality. Due to the more uniform noise properties, we expect this method to work relatively better in photoacoustics than, say, x-ray CT, which would benefit from a more advanced regularization. Using this method the reconstruction step consists only of a matrix-vector multiplication, which is very fast. In a recent publication, the system matrix of a photoacoustic imaging system was measured experimentally, and the corresponding pseudoinverse was calculated ahead of time and used for real-time imaging.8 This method has the advantage of including all the relevant acoustic properties in the system matrix, but the disadvantages of requiring the experimental measurement of the system matrix and the inevitable inclusion of system noise in the measurement. In this paper we develop a model for a two-dimensional (2-D) system matrix and use that to calculate the pseudoinverse. We also develop a three-dimensional (3-D) forward model to simulate a photoacoustic imaging system.

A natural question then arises as how best to select an algorithm to optimize image quality. In a previous work,9 we performed a preliminary comparison between the filtered backprojection and pseudoinverse reconstruction algorithms. In that paper, we looked at various metrics such as resolution, streak length, and ideal observer signal-to-noise ration (SNR) for a point source object; however, these metrics can be misleading. The resolution and streak length tell nothing of the noise properties of the algorithms, and the ideal observer SNR only measures how much information is lost for each algorithm in going from the raw data to the reconstructed image, which does not necessarily correlate with image quality as judged by human observer task performance.

In this paper we demonstrate the use in photoacoustics of so-called task-based metrics of image quality that seek to define image quality based not on physical measures but on the utility of an image for performing a predefined task. Such task-based image-quality assessment has gained widespread attention in nuclear medicine and x-ray imaging modalities.10 While the gold standard of task-based assessment is human observer studies using the receiver operating characteristic (ROC) methodology, these are tedious and labor-intensive. As an alternative, various model observers have been proposed that seek to emulate human observer performance in well-defined tasks, such as the detection of a known signal in a known or random background. The channelized Hotelling observer (CHO) has been shown to predict human performance particularly well.10 It entails integrating images over a set of frequency-domain channels that mimic the frequency response of the human visual system and then applying a linear discriminant to the output. To assess the image quality of the reconstruction algorithms in this study, we use the CHO with two sets of channels that have been shown to predict the performance of human observers: constant Q11 and difference of Gaussian (DOG) channels.12

The photoacoustic equation can be written as1Display Formula

(21vs22t2)p(r,t)=βCPH(r,t)t,(1)
where vs is the speed of sound, p(r,t) is the pressure, β is the thermal coefficient of volume expansion, CP is the specific heat at constant pressure, and H(r,t) is the heating function. The Green’s function solution for the pressure spectrum can be written as13Display Formula
p˜(r,ω)=iωβ4πCPV0d3r0H˜(r0,ω)eik|rr0||rr0|,(2)
where k=ω/vs. The geometry we are considering is shown in Fig. 1 and consists of a linear transducer array located on the x axis with a cylindrical lens on its face and consisting of N individual transducers, an object to be imaged in the half-space z0>0, and a light source above the object pointing in the positive y0 direction. The cylindrical lens rejects acoustic waves that do not originate near the imaging (x0z0) plane. We assume uniform light illumination on the imaging plane to avoid modeling photon propagation. With this geometry, we can separate the heating function spectrum into H˜(r0,ω)=μ(r0)I˜(ω), where μ(r0) is the absorption coefficient and I˜(ω) is the irradiance spectrum. The average pressure on each individual transducer is the pressure integrated over x and y divided by the area, where the width in the x direction is w. The cylindrical acoustic lens imparts a phase of eikdeiky2/(2f) to the wave, where d is the acoustic path length of the lens (refractive index times total thickness of the lens) and f is the focal length. The lens has a height of h and the transducer has a transfer function denoted by T˜(ω). This gives the measured voltage spectrum as Display Formula
V˜n(ω)=iωβI˜(ω)T˜(ω)eikd4πhwCPw/2+nww/2+nwh/2h/2dxdyeiky2/(2f)V0d3r0μ(r0)z0eik|rr0||rr0|2,(3)
where n denotes each individual transducer and runs from (N1)/2 to (N1)/2. The middle transducer has n=0, and the center of the middle transducer is located at x=0. The extra z0/|rr0| is an obliquity factor.

Graphic Jump LocationF1 :

Geometry of photoacoustic imaging system.

In this paper we consider a 32-element transducer with cylindrical lens focal length of 6 mm, height of 3 mm, individual element width of 0.15 mm, and center-to-center element spacing of 0.15 mm. The transducer transfer function is Gaussian with a center frequency of 3.5 MHz and 50% bandwidth, and we assume the laser pulse is sufficiently short so that I˜(ω)T˜(ω)=I0T˜(ω), where I0 is a constant. Note that the transducer spacing satisfies the Nyquist sampling criterion. For the forward model, we use the full 3-D expression given by Eq. (3). For the system matrix used to calculate the pseudoinverse, we use a 2-D approximation of Eq. (3) with y=y0=0. For the discrete version, we use the notation g=Hf, where g is a vector representing the voltage measured by the transducers, H is the system matrix, and f is a vector representing the initial pressure distribution of the object being imaged. We use a 128×128pixel basis with a pixel size of 37.5 μm. The temporal sampling rate is 100 ns with 71 samples per transducer, which produces a system matrix of size 2272×16384.

Pseudoinverse

We performed an SVD on the system matrix, which was then used to calculate the pseudoinverse. The SVD factors the system matrix into H=USV, where U and V are unitary matrices, the dagger represents conjugate transpose, and S is a diagonal matrix consisting of the singular values in decreasing order. The pseudoinverse is then given by H+=VS+U, where S+ is a diagonal matrix consisting of the reciprocal of singular values above a threshold. The thresholding of the singular values serves as a regularization. Note that this method is also referred to as truncated SVD inversion. Figure 2 shows a plot of the magnitude of the singular values. The image, which is an estimation of the initial pressure f, is then given by f^=H+g.

Graphic Jump LocationF2 :

Magnitude of singular values of system matrix.

Filtered Backprojection and Simple Backprojection

The filtered backprojection formula5 can be written as Display Formula

μ(x0,z0)=z0ngn(τn)*R(fcτn),(4)
where gn(τn) is the time integrated pressure (along with other constants), the convolution is with respect to τn=[(nwx0)2+z02]1/2/vs, R(u)=4sinc(2u)2sinc2(u), sinc(u)=sin(πu)/(πu), and fc is approximately the highest temporal frequency in the data. Note that in the Fourier domain, the filter is the Ram-Lak filter (ramp filter with cutoff fc). In order to construct the filtered backprojection matrix, we used impulses for gn(τn) that were time-integrated and interpolated to increase the sampling rate to prevent quantization errors.14 A Hamming window was applied along the transducer elements for apodization. The simple backprojection algorithm is equivalent to the filtered backprojection algorithm except the pressure is not time-integrated or filtered with the Ram-Lak filter. Note that this filtered backprojection formula is built on a foundation developed by Norton in the context of reflection tomography.15 An alternative algorithm has been developed by Burgholzer et al. but has not been implemented in the present study.16

Fourier Method

The Fourier transform-based formula6 can be written as Display Formula

μ˜(kx,kz)=kzkV˜(kx,ω),(5)
where k2=kx2+kz2. V˜(kx,ω) is constructed by taking the real part of the Fourier transform with respect to time, and taking the Fourier transform with respect to x. The image is formed by employing bilinear interpolation and taking the inverse Fourier transform of μ˜(kx,kz). Again, a Hamming window is used for apodization. We are ignoring multiplicative constants that would be required for quantitative results, as we did for the backprojection algorithms.

Deconvolution

For the filtered backprojection, simple backprojection, and Fourier based methods we deconvolve the transducer response using Display Formula

D˜(ω)=T˜*(ω)|T˜(ω)|2+ζ,(6)
where the asterisk indicates complex conjugation and ζ is a regularization parameter that varies from 106 to 10. For small values of ζ the deconvolution is the inverse of the transducer response and noise amplification occurs. When ζ is 10, D˜(ω)T˜*(ω), which is a matched filter. The matched filter lowers the noise at the expense of resolution. For all of the reconstructed images we perform envelope detection by taking the absolute value of the analytic function of the image.

Figures 3 and 4 show the image reconstruction of a point source and a 0.5-mm cube, respectively, for various singular-value truncations as well as the backprojection, filtered backprojection, and Fourier-based algorithms with different regularization parameters. Uncorrelated Gaussian noise is added to g before reconstruction to simulate electronic noise. The standard deviation of the noise is set to be equal to the peak signal level divided by 10. There is no background for these images. Figure 5 shows the image reconstruction of a lumpy background17 for the different algorithms. For each lumpy background image, there is a corresponding image with lumpy background and point source signal. All of the images in Fig. 5 contain the same background; in the next section the number of lumps and their locations are treated as random variables. The lumpy background is only present in a volume represented by a 0.5 mm cube at the center of the object. Uncorrelated Gaussian noise is also added to g before reconstruction, where, again, the standard deviation is set to the peak signal level divided by 10.

Graphic Jump LocationF3 :

Image reconstruction of a point source for singular-value truncations of (a) 105; (b) 106; (c) 107; (d) 108; filtered backprojection with regularization parameters of (e) 105; (f) 103; (g) 101; (h) 101; backprojection with regularization parameters of (i) 105; (j) 103; (k) 101; (l) 101; and the Fourier-based method with regularization parameters of (m) 105; (n) 103; (o) 101; (p) 101.

Graphic Jump LocationF4 :

Image reconstruction of a 0.5-mm cube for singular-value truncations of (a) 105; (b) 106; (c) 107; (d) 108; filtered backprojection with regularization parameters of (e) 105; (f) 103; (g) 101; (h) 101; backprojection with regularization parameters of (i) 105; (j) 103; (k) 101; (l) 101; and the Fourier-based method with regularization parameters of (m) 105; (n) 103; (o) 101; (p) 101.

Graphic Jump LocationF5 :

Image reconstruction of a lumpy background for singular-value truncations of (a) 107; (b) 107; (c) 108; (d) 108; filtered backprojection with regularization parameters of (e) 102; (f) 102; (g) 101; (h) 101; backprojection with regularization parameters of (i) 102; (j) 102; (k) 101; (l) 101; and the Fourier-based method with regularization parameters of (m) 102; (n) 102; (o) 101; (p) 101. A point source signal is also included in (b); (d); (f); (h); (j); (l); (n);, and (p). The lumpy background is only present in a volume represented by a 0.5-mm cube at the center of the object.

As mentioned in the Introduction, the CHO has been shown to predict human observer performance for various tasks. The choice of channels depends on the specific task to be performed. In this paper we use two sets of channels: constant Q11 and DOG channels.12 Constant Q channels have been shown to work well at predicting human performance when the noise is additive and Gaussian. DOG channels work well for additive Gaussian noise as well as for a lumpy background. The channel profiles are shown in Fig. 6.

Graphic Jump LocationF6 :

Channel profiles for (a) constant Q channels and (b) DOG channels.

In order to quantify image quality, we compute the channelized Hotelling SNR given by11Display Formula

SNRch2=f^chKch1f^ch,(7)
where the brackets indicate expectation values, f^ch=Uchf^, Uch is a matrix that represents the spatial frequency filtering and summing of the channels, and Kch is the covariance matrix of f^ch. Note that the covariance matrix is only a 5×5 matrix for the constant Q channels and a 3×3 matrix for the DOG channels, which means the number of simulations required to make the matrix full rank is greatly reduced as compared to the full image covariance matrix. In this study we use 500 simulations to construct Eq. (7). For the uniform background tasks, uncorrelated Gaussian noise is added to g before reconstruction, as discussed in the previous section. These tasks are referred to as signal known exactly, background known exactly (SKE/BKE). For the lumpy background17 task, the number of lumps and their locations are treated as random variables. The number of lumps, L, is modeled as a Gaussian random variable with a mean and variance of 1250. That number was selected using L/N=1/100, where N=1.25×105 is the number of discrete object points in a 0.5-mm cube in our simulations. The location of the lumps is modeled as a uniform random variable over the volume represented by a 0.5-mm cube at the center of the object. Each lump is modeled as a point source with a magnitude that is 100 times lower than the point source signal magnitude. Uncorrelated Gaussian noise is also added to g before reconstruction for this task. This task is referred to as SKE in a lumpy background and is a good model of tumor contrast agent imaging where the signal represents uptake of the contrast agent by the tumor and the lumpy background represents contrast agent uptake by the surrounding tissue.

We define the CHO efficiency as Display Formula

η=SNRch2SNRdata2,(8)
where SNRdata2 is given by Display Formula
SNRdata2=gKg1g.(9)

In Figs. 78 through 9, we show the CHO efficiency for the three tasks. For the SKE/BKE tasks, the efficiency is highest when the regularization is chosen so that the noise is lowest even though the resolution is worse as seen in Figs. 3 and 4. The performance of the algorithms is similar for the SKE/BKE tasks except that the Fourier-based algorithm performs best for the identification of the cube as determined by the difference of Gaussian channels. However, the SKE/BKE task is generally not a realistic model of actual imaging systems. For the SKE in a lumpy background task, which is more realistic, the method based on the pseudoinverse performs best according to both sets of channels. Also notice that for this task the optimal regularization is not necessarily to choose the lowest possible noise amplification. It is also interesting that for the three tasks studied in this paper, the simple backprojection performs as well as the filtered backprojection.

Graphic Jump LocationF7 :

CHO efficiency for the (a) constant Q channels and (b) DOG channels for the point source SKE/BKE task.

Graphic Jump LocationF8 :

CHO efficiency for the (a) constant Q channels and (b) DOG channels for the 0.5-mm cube SKE/BKE task.

Graphic Jump LocationF9 :

CHO efficiency for the (a) constant Q channels and (b) DOG channels for the point source SKE in a lumpy background task.

In this paper, we have used task-based image quality metrics to compare image reconstruction algorithms in photoacoustic tomography. The study was not intended to be the final word on algorithm selection but rather to introduce the use of task-based assessment and model observers to the field. As the field grows and matures, it could be valuable to deploy such tools for optimizing systems and algorithms for various tasks.

We have shown that for the case of a linear transducer array and an SKE in a lumpy background task, the image reconstruction algorithm based on the pseudoinverse performs best according to the CHOs we look at in this paper. Future work is planned to study the validity of the channels used in this paper by performing human observer experiments. Future work is also planned to study tasks that include signal variability. It should be noted that in this paper we assume uniform light illumination on the imaging plane. Depending on the specific system and task to be performed this may not be a realistic assumption, and in this case the laser irradiation geometry and photon propagation should be included in the forward model.

Another interesting conclusion of this study is that the simple backprojection algorithm performs as well as the filtered backprojection algorithm for the tasks we looked at. Since the simple backprojection algorithm is numerically equivalent to the delay-and-sum beamformer reconstruction employed in existing ultrasound imaging systems, combined ultrasound-photoacoustic imaging systems may be able to use this algorithm without loss of performance. This would be advantageous in terms of system design. The type of analysis presented in this paper could be used to determine if this conclusion is valid for the specific task the combined ultrasound-photoacoustic imaging system is being built for.

This work was supported by ACS Research Scholar Grant RSG-08-117-01-CCE to Patrick La Rivière.

Li  C., Wang  L. V., “Photoacoustic tomography and sensing in biomedicine,” Phys. Med. Biol.. 54, (19 ), R59 –R96 (2009). 0031-9155 CrossRef
Niederhauser  J. J. et al., “Combined ultrasound and optoacoustic system for real-time high-contrast vascular imaging in vivo,” IEEE Trans. Med. Imaging. 24, (4 ), 436 –440 (2005). 0278-0062 CrossRef
Yang  D. W. et al., “Fast full-view photoacoustic imaging by combined scanning with a linear transducer array,” Opt. Express. 15, (23 ), 15566 –15575 (2007). 1094-4087 CrossRef
Olafsson  R. et al., “Real-time, contrast enhanced photoacoustic imaging of cancer in a mouse window chamber,” Opt. Express. 18, (18 ), 18625 –18632 (2010). 1094-4087 CrossRef
Modgil  D., La Rivière  P. J., “Implementation and comparison of reconstruction algorithms for two-dimensional optoacoustic tomography using a linear array,” J. Biomed. Opt.. 14, (4 ), 044023  (2009). 1083-3668 CrossRef
Köstli  K. P., Beard  P. C., “Two-dimensional photoacoustic imaging by use of Fourier-transform image reconstruction and a detector with an anisotropic response,” Appl. Opt.. 42, (10 ), 1899 –1908 (2003). 0003-6935 CrossRef
Zhang  J. et al., “Simultaneous reconstruction of speed-of-sound and optical absorption properties in photoacoustic tomography via a time-domain iterative algorithm,” Proc. SPIE. 6856, , 68561F  (2008). 0277-786X CrossRef
Roumeliotis  M. B. et al., “Singular value decomposition analysis of a photoacoustic imaging system and 3D imaging at 0.7 FPS,” Opt. Express. 19, (14 ), 13405 –13417 (2011). 1094-4087 CrossRef
Petschke  A., La Rivière  P. J., “Photoacoustic image reconstruction using the pseudoinverse of the system matrix with the potential for real time imaging,” Proc. SPIE. 8223, , 82233J  (2012). 0277-786X CrossRef
Barrett  H. H., Myers  K. J., Foundations of Image Science. , pp. 913 –1000,  Wiley ,  Hoboken, NJ  (2004).
Myers  K. J., Barrett  H. H., “Addition of a channel mechanism to the ideal-observer model,” J. Opt. Soc. Am. A. 4, (12 ), 2447 –2457 (1987). 0740-3232 CrossRef
Abbey  C. K., Barrett  H. H., Wilson  D. W., “Observer signal-to-noise ratios for the ML-EM algorithm,” Proc. SPIE. 2712, , 47 –58 (1996). 0277-786X CrossRef
Petschke  A., La Rivière  P. J., “Comparison of intensity-modulated continuous-wave lasers with a chirped modulation frequency to pulsed lasers for photoacoustic imaging applications,” Biomed. Opt. Express. 1, (4 ), 1188 –1195 (2010). 2156-7085 CrossRef
Cobbold  R. S. C., Foundations of Biomedical Ultrasound. , pp. 492 –607,  Oxford ,  New York , (2006).
Norton  S. J., “Reconstruction of a reflectivity field from line integrals over circular paths,” J. Acoust. Soc. Am.. 67, (3 ), 853 –863 (1980). 0001-4966 CrossRef
Burgholzer  P. et al., “Temporal back-projection algorithms for photoacoustic tomography with integrating line detectors,” Inverse Probl.. 23, (6 ), S65 –S80 (2007). 0266-5611 CrossRef
Rolland  J. P., Barrett  H. H., “Effect of random background inhomogeneity on observer detection performance,” J. Opt. Soc. Am. A. 9, (5 ), 649 –658 (1992). 0740-3232 CrossRef
© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Adam Petschke and Patrick J. La Rivière
"Comparison of photoacoustic image reconstruction algorithms using the channelized Hotelling observer", J. Biomed. Opt. 18(2), 026009 (Feb 05, 2013). ; http://dx.doi.org/10.1117/1.JBO.18.2.026009


Figures

Graphic Jump LocationF1 :

Geometry of photoacoustic imaging system.

Graphic Jump LocationF4 :

Image reconstruction of a 0.5-mm cube for singular-value truncations of (a) 105; (b) 106; (c) 107; (d) 108; filtered backprojection with regularization parameters of (e) 105; (f) 103; (g) 101; (h) 101; backprojection with regularization parameters of (i) 105; (j) 103; (k) 101; (l) 101; and the Fourier-based method with regularization parameters of (m) 105; (n) 103; (o) 101; (p) 101.

Graphic Jump LocationF3 :

Image reconstruction of a point source for singular-value truncations of (a) 105; (b) 106; (c) 107; (d) 108; filtered backprojection with regularization parameters of (e) 105; (f) 103; (g) 101; (h) 101; backprojection with regularization parameters of (i) 105; (j) 103; (k) 101; (l) 101; and the Fourier-based method with regularization parameters of (m) 105; (n) 103; (o) 101; (p) 101.

Graphic Jump LocationF2 :

Magnitude of singular values of system matrix.

Graphic Jump LocationF9 :

CHO efficiency for the (a) constant Q channels and (b) DOG channels for the point source SKE in a lumpy background task.

Graphic Jump LocationF8 :

CHO efficiency for the (a) constant Q channels and (b) DOG channels for the 0.5-mm cube SKE/BKE task.

Graphic Jump LocationF7 :

CHO efficiency for the (a) constant Q channels and (b) DOG channels for the point source SKE/BKE task.

Graphic Jump LocationF6 :

Channel profiles for (a) constant Q channels and (b) DOG channels.

Graphic Jump LocationF5 :

Image reconstruction of a lumpy background for singular-value truncations of (a) 107; (b) 107; (c) 108; (d) 108; filtered backprojection with regularization parameters of (e) 102; (f) 102; (g) 101; (h) 101; backprojection with regularization parameters of (i) 102; (j) 102; (k) 101; (l) 101; and the Fourier-based method with regularization parameters of (m) 102; (n) 102; (o) 101; (p) 101. A point source signal is also included in (b); (d); (f); (h); (j); (l); (n);, and (p). The lumpy background is only present in a volume represented by a 0.5-mm cube at the center of the object.

Tables

References

Li  C., Wang  L. V., “Photoacoustic tomography and sensing in biomedicine,” Phys. Med. Biol.. 54, (19 ), R59 –R96 (2009). 0031-9155 CrossRef
Niederhauser  J. J. et al., “Combined ultrasound and optoacoustic system for real-time high-contrast vascular imaging in vivo,” IEEE Trans. Med. Imaging. 24, (4 ), 436 –440 (2005). 0278-0062 CrossRef
Yang  D. W. et al., “Fast full-view photoacoustic imaging by combined scanning with a linear transducer array,” Opt. Express. 15, (23 ), 15566 –15575 (2007). 1094-4087 CrossRef
Olafsson  R. et al., “Real-time, contrast enhanced photoacoustic imaging of cancer in a mouse window chamber,” Opt. Express. 18, (18 ), 18625 –18632 (2010). 1094-4087 CrossRef
Modgil  D., La Rivière  P. J., “Implementation and comparison of reconstruction algorithms for two-dimensional optoacoustic tomography using a linear array,” J. Biomed. Opt.. 14, (4 ), 044023  (2009). 1083-3668 CrossRef
Köstli  K. P., Beard  P. C., “Two-dimensional photoacoustic imaging by use of Fourier-transform image reconstruction and a detector with an anisotropic response,” Appl. Opt.. 42, (10 ), 1899 –1908 (2003). 0003-6935 CrossRef
Zhang  J. et al., “Simultaneous reconstruction of speed-of-sound and optical absorption properties in photoacoustic tomography via a time-domain iterative algorithm,” Proc. SPIE. 6856, , 68561F  (2008). 0277-786X CrossRef
Roumeliotis  M. B. et al., “Singular value decomposition analysis of a photoacoustic imaging system and 3D imaging at 0.7 FPS,” Opt. Express. 19, (14 ), 13405 –13417 (2011). 1094-4087 CrossRef
Petschke  A., La Rivière  P. J., “Photoacoustic image reconstruction using the pseudoinverse of the system matrix with the potential for real time imaging,” Proc. SPIE. 8223, , 82233J  (2012). 0277-786X CrossRef
Barrett  H. H., Myers  K. J., Foundations of Image Science. , pp. 913 –1000,  Wiley ,  Hoboken, NJ  (2004).
Myers  K. J., Barrett  H. H., “Addition of a channel mechanism to the ideal-observer model,” J. Opt. Soc. Am. A. 4, (12 ), 2447 –2457 (1987). 0740-3232 CrossRef
Abbey  C. K., Barrett  H. H., Wilson  D. W., “Observer signal-to-noise ratios for the ML-EM algorithm,” Proc. SPIE. 2712, , 47 –58 (1996). 0277-786X CrossRef
Petschke  A., La Rivière  P. J., “Comparison of intensity-modulated continuous-wave lasers with a chirped modulation frequency to pulsed lasers for photoacoustic imaging applications,” Biomed. Opt. Express. 1, (4 ), 1188 –1195 (2010). 2156-7085 CrossRef
Cobbold  R. S. C., Foundations of Biomedical Ultrasound. , pp. 492 –607,  Oxford ,  New York , (2006).
Norton  S. J., “Reconstruction of a reflectivity field from line integrals over circular paths,” J. Acoust. Soc. Am.. 67, (3 ), 853 –863 (1980). 0001-4966 CrossRef
Burgholzer  P. et al., “Temporal back-projection algorithms for photoacoustic tomography with integrating line detectors,” Inverse Probl.. 23, (6 ), S65 –S80 (2007). 0266-5611 CrossRef
Rolland  J. P., Barrett  H. H., “Effect of random background inhomogeneity on observer detection performance,” J. Opt. Soc. Am. A. 9, (5 ), 649 –658 (1992). 0740-3232 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.