1
|
Dong W, Zhu C, Xie D, Zhang Y, Tao S, Tian C. Image restoration for ring-array photoacoustic tomography system based on blind spatially rotational deconvolution. PHOTOACOUSTICS 2024; 38:100607. [PMID: 38665365 PMCID: PMC11044036 DOI: 10.1016/j.pacs.2024.100607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 03/17/2024] [Accepted: 04/09/2024] [Indexed: 04/28/2024]
Abstract
Ring-array photoacoustic tomography (PAT) system has been widely used in noninvasive biomedical imaging. However, the reconstructed image usually suffers from spatially rotational blur and streak artifacts due to the non-ideal imaging conditions. To improve the reconstructed image towards higher quality, we propose a concept of spatially rotational convolution to formulate the image blur process, then we build a regularized restoration problem model accordingly and design an alternating minimization algorithm which is called blind spatially rotational deconvolution to achieve the restored image. Besides, we also present an image preprocessing method based on the proposed algorithm to remove the streak artifacts. We take experiments on phantoms and in vivo biological tissues for evaluation, the results show that our approach can significantly enhance the resolution of the image obtained from ring-array PAT system and remove the streak artifacts effectively.
Collapse
Affiliation(s)
- Wende Dong
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 211106, China
- Key Laboratory of Space Photoelectric Detection and Perception (Nanjing University of Aeronautics and Astronautics), Ministry of Industry and Information Technology, Nanjing, Jiangsu 211106, China
| | - Chenlong Zhu
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 211106, China
- Key Laboratory of Space Photoelectric Detection and Perception (Nanjing University of Aeronautics and Astronautics), Ministry of Industry and Information Technology, Nanjing, Jiangsu 211106, China
| | - Dan Xie
- School of Engineering Science, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yanli Zhang
- College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 211106, China
- Key Laboratory of Space Photoelectric Detection and Perception (Nanjing University of Aeronautics and Astronautics), Ministry of Industry and Information Technology, Nanjing, Jiangsu 211106, China
| | - Shuyin Tao
- School of Computer Science and Engineering, Nanjing University of Science and Technology, Nanjing, Jiangsu 210094, China
| | - Chao Tian
- School of Engineering Science, University of Science and Technology of China, Hefei, Anhui 230026, China
- Institute of Artificial Intelligence, Hefei Comprehensive National Science Center, Hefei, Anhui 230088, China
- Anhui Province Key Laboratory of Biomedical Imaging and Intelligent Processing, Institute of Artificial Intelligence, Hefei Comprehensive National Science Center, Hefei, Anhui 230088, China
| |
Collapse
|
2
|
Qin Z, Ma Y, Ma L, Liu G, Sun M. Convolutional sparse coding for compressed sensing photoacoustic CT reconstruction with partially known support. BIOMEDICAL OPTICS EXPRESS 2024; 15:524-539. [PMID: 38404320 PMCID: PMC10890869 DOI: 10.1364/boe.507831] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Revised: 12/01/2023] [Accepted: 12/14/2023] [Indexed: 02/27/2024]
Abstract
In photoacoustic tomography (PAT), imaging speed is an essential metric that is restricted by the pulse laser repetition rate and the number of channels on the data acquisition card (DAQ). Reconstructing the initial sound pressure distribution with fewer elements can significantly reduce hardware costs and back-end acquisition pressure. However, undersampling will result in artefacts in the photoacoustic image, degrading its quality. Dictionary learning (DL) has been utilised for various image reconstruction techniques, but they disregard the uniformity of pixels in overlapping blocks. Therefore, we propose a compressive sensing (CS) reconstruction algorithm for circular array PAT based on gradient domain convolutional sparse coding (CSCGR). A small number of non-zero signal positions in the sparsely encoded feature map are used as partially known support (PKS) in the reconstruction procedure. The CS-CSCGR-PKS-based reconstruction algorithm can use fewer ultrasound transducers for signal acquisition while maintaining image fidelity. We demonstrated the effectiveness of this algorithm in sparse imaging through imaging experiments on the mouse torso, brain, and human fingers. Reducing the number of array elements while ensuring imaging quality effectively reduces equipment hardware costs and improves imaging speed.
Collapse
Affiliation(s)
- Zezheng Qin
- School of Astronautics, Harbin Institute of Technology, Harbin 150000, China
| | - Yiming Ma
- School of Astronautics, Harbin Institute of Technology, Harbin 150000, China
| | - Lingyu Ma
- School of Astronautics, Harbin Institute of Technology, Harbin 150000, China
| | - Guangxing Liu
- Suzhou Institute of Biomedical Engineering Technology, Chinese Academy of Sciences, Suzhou 215163, China
- College of Biomedical Engineering, University of Science and Technology of China, Hefei 230026, China
| | - Mingjian Sun
- School of Astronautics, Harbin Institute of Technology, Harbin 150000, China
| |
Collapse
|
3
|
Xie D, Dong W, Zheng J, Tian C. Spatially-variant image deconvolution for photoacoustic tomography. OPTICS EXPRESS 2023; 31:21641-21657. [PMID: 37381257 DOI: 10.1364/oe.486846] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 05/25/2023] [Indexed: 06/30/2023]
Abstract
Photoacoustic tomography (PAT) system can reconstruct images of biological tissues with high resolution and contrast. However, in practice, the PAT images are usually degraded by spatially variant blur and streak artifacts due to the non-ideal imaging conditions and chosen reconstruction algorithms. Therefore, in this paper, we propose a two-phase restoration method to progressively improve the image quality. In the first phase, we design a precise device and measuring method to obtain spatially variant point spread function samples at preset positions of the PAT system in image domain, then we adopt principal component analysis and radial basis function interpolation to model the entire spatially variant point spread function. Afterwards, we propose a sparse logarithmic gradient regularized Richardson-Lucy (SLG-RL) algorithm to deblur the reconstructed PAT images. In the second phase, we present a novel method called deringing which is also based on SLG-RL to remove the streak artifacts. Finally, we evaluate our method with simulation, phantom and in vivo experiments, respectively. All the results show that our method can significantly improve the quality of PAT images.
Collapse
|
4
|
Zheng S, Yingsa H, Meichen S, Qi M. Quantitative photoacoustic tomography with light fluence compensation based on radiance Monte Carlo model. Phys Med Biol 2023; 68. [PMID: 36821863 DOI: 10.1088/1361-6560/acbe90] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Accepted: 02/23/2023] [Indexed: 02/25/2023]
Abstract
Objective. Photoacoustic tomography (PAT) is a rapidly evolving imaging modality that provides images with high contrast and spatial resolution showing the optical properties of biological tissues. The photoacoustic pressure is proportional to the product of the optical absorption coefficient and the local light fluence. The essential challenge in reconstructing quantitative images representing spatially varying absorption coefficients is the unknown light fluence. In addition, optical attenuation induces spatial variations in the light fluence, and the heterogeneity of the fluence determines the limits of reconstruction quality and depth.Approach.In this work, a reconstruction enhancement scheme is proposed to compensate for the variation in the light fluence in the absorption coefficient recovery. The inverse problem of the radiance Monte Carlo model describing light transport through the tissue is solved by using an alternating optimization strategy. In the iteration, the absorption coefficients and photon weights are alternately updated.Main results.The method provides highly accurate quantitative images of absorption coefficients in simulations, phantoms, andin vivostudies. The results show that the method has great potential for improving the accuracy of absorption coefficient recovery compared to conventional reconstruction methods that ignore light fluence variations. Comparison with state-of-the-art fluence compensation methods shows significant improvements in root mean square error, normalized mean square absolute distance, and structural similarity metrics.Significance.This method achieves high precision quantitative imaging by compensating for nonuniform light fluence without increasing the complexity and operation of the imaging system.
Collapse
Affiliation(s)
- Sun Zheng
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
- Hebei Key Laboratory of Power Internet of Things Technology, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| | - Hou Yingsa
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| | - Sun Meichen
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| | - Meng Qi
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| |
Collapse
|
5
|
Yip LCM, Omidi P, Rascevska E, Carson JJL. Approaching closed spherical, full-view detection for photoacoustic tomography. JOURNAL OF BIOMEDICAL OPTICS 2022; 27:JBO-220034GRR. [PMID: 36042544 PMCID: PMC9424748 DOI: 10.1117/1.jbo.27.8.086004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Accepted: 07/01/2022] [Indexed: 05/28/2023]
Abstract
SIGNIFICANCE Photoacoustic tomography (PAT) is a widely explored imaging modality and has excellent potential for clinical applications. On the acoustic detection side, limited-view angle and limited-bandwidth are common key issues in PAT systems that result in unwanted artifacts. While analytical and simulation studies of limited-view artifacts are fairly extensive, experimental setups capable of comparing limited-view to an ideal full-view case are lacking. AIMS A custom ring-shaped detector array was assembled and mounted to a 6-axis robot, then rotated and translated to achieve up to 3.8π steradian view angle coverage of an imaged object. APPROACH Minimization of negativity artifacts and phantom imaging were used to optimize the system, followed by demonstrative imaging of a star contrast phantom, a synthetic breast tumor specimen phantom, and a vascular phantom. RESULTS Optimization of the angular/rotation scans found ≈212 effective detectors were needed for high-quality images, while 15-mm steps were used to increase the field of view as required depending on the size of the imaged object. Example phantoms were clearly imaged with all discerning features visible and minimal artifacts. CONCLUSIONS A near full-view closed spherical system has been developed, paving the way for future work demonstrating experimentally the significant advantages of using a full-view PAT setup.
Collapse
Affiliation(s)
- Lawrence C. M. Yip
- Lawson Health Research Institute, Imaging Program, London, Ontario, Canada
- Western University, Schulich School of Medicine and Dentistry, Department of Medical Biophysics, London, Ontario, Canada
| | - Parsa Omidi
- Lawson Health Research Institute, Imaging Program, London, Ontario, Canada
- Western University, School of Biomedical Engineering, London, Ontario, Canada
| | - Elina Rascevska
- Lawson Health Research Institute, Imaging Program, London, Ontario, Canada
- Western University, School of Biomedical Engineering, London, Ontario, Canada
| | - Jeffrey J. L. Carson
- Lawson Health Research Institute, Imaging Program, London, Ontario, Canada
- Western University, Schulich School of Medicine and Dentistry, Department of Medical Biophysics, London, Ontario, Canada
- Western University, School of Biomedical Engineering, London, Ontario, Canada
- Western University, Schulich School of Medicine and Dentistry, Department of Surgery, London, Ontario, Canada
| |
Collapse
|
6
|
Sun Z, Sun H. Image reconstruction for endoscopic photoacoustic tomography including effects of detector responses. Exp Biol Med (Maywood) 2022; 247:881-897. [PMID: 35232296 DOI: 10.1177/15353702221079570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In photoacoustic tomography (PAT), conventional image reconstruction methods are generally based on the assumption of an ideal point-like ultrasonic detector. This assumption is appropriate when the receiving surface of the detector is sufficiently small and/or the distance between the imaged object and the detector is large enough. However, it does not hold in endoscopic applications of PAT. In this study, we propose a model-based image reconstruction method for endoscopic photoacoustic tomography (EPAT), considering the effect of detector responses on image quality. We construct a forward model to physically describe the imaging process of EPAT, including the generation of the initial pressure due to optical absorption and thermoelastic expansion, the propagation of photoacoustic waves in tissues, and the acoustic measurement. The model outputs the theoretical sampling voltage signal, which is the response of the ultrasonic detector to the acoustic pressure reaching its receiving surface. The images representing the distribution map of the optical absorption energy density on cross-sections of the imaged luminal structures are reconstructed from the sampling voltage signals output by the detector through iterative inversion of the forward model. Compared with the conventional approaches based on back-projection and other imaging models, our method improved the quality and spatial resolution of the resulting images.
Collapse
Affiliation(s)
- Zheng Sun
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, China.,Hebei Key Laboratory of Power Internet of Things Technology, North China Electric Power University, Baoding 071003, China
| | - Huifeng Sun
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, China.,Hebei Key Laboratory of Power Internet of Things Technology, North China Electric Power University, Baoding 071003, China
| |
Collapse
|
7
|
Fast Correction of “Finite Aperture Effect” in Photoacoustic Tomography Based on Spatial Impulse Response. PHOTONICS 2021. [DOI: 10.3390/photonics8090356] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Photoacoustic computed tomography (PACT) is a fast-developing imaging technique, which can provide structural and functional information in biological tissues with high-resolution beyond the depth of the optical diffusion limit. However, the most current PACT reconstruction method generally employs a point detector assumption, whereas in most PAT systems with circular or spherical scanning modes, the transducer is mostly flat and with a finite size. This model mismatch leads to a notable deterioration in the lateral direction in regions far from the rotation center, which is known as the “finite aperture effect”. In this work, we propose to compensate a novel Back-projection (BP) method based on the transducer’s spatial impulse response (SIR) for fast correction of the “finite aperture effect”. The SIR accounts for the waveform change of the transducer for an arbitrary point source due to the geometry of the detection surface. Simulation results showed that the proposed SIR-BP method can effectively improve the lateral resolution and signal to noise ratio (SNR) in the off-center regions. For a target 4.5 mm far from the rotation center, this new method improved the lateral resolution about five times along with a 7 dB increase in the SNR. Experimental results also showed that this SIR-BP method can well restore the image angular blur to recover small structures, as demonstrated by the imaging of leaf veins. This new method offers a valuable alternative to the conventional BP method, and can guide the design of PAT systems based on circular/spherical scan.
Collapse
|
8
|
Melchert O, Wollweber M, Roth B. Optoacoustic inversion via convolution kernel reconstruction in the paraxial approximation and beyond. PHOTOACOUSTICS 2019; 13:1-5. [PMID: 30510898 PMCID: PMC6257913 DOI: 10.1016/j.pacs.2018.10.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/23/2018] [Revised: 10/05/2018] [Accepted: 10/23/2018] [Indexed: 06/09/2023]
Abstract
In this article we address the numeric inversion of optoacoustic signals to initial stress profiles. Therefore we study a Volterra integral equation of the second kind that describes the shape transformation of propagating stress waves in the paraxial approximation of the underlying wave-equation. Expanding the optoacoustic convolution kernel in terms of a Fourier-series, a best fit to a pair of observed near-field and far-field signals allows to obtain a sequence of expansion coefficients that describe a given "apparative" setup. The resulting effective kernel is used to solve the optoacoustic source reconstruction problem using a Picard-Lindelöf correction scheme. We verify the validity of the proposed inversion protocol for synthetic input signals and explore the feasibility of our approach to also account for the shape transformation of signals beyond the paraxial approximation including the inversion of experimental data stemming from measurements on melanin doped PVA hydrogel tissue phantoms.
Collapse
Affiliation(s)
- O. Melchert
- Hannover Centre for Optical Technologies (HOT), Interdisciplinary Research Centre of the Leibniz Universität Hannover, Nienburger Str. 17, D-30167 Hannover, Germany
| | | | | |
Collapse
|
9
|
Wang J, Wang Y. Photoacoustic imaging reconstruction using combined nonlocal patch and total-variation regularization for straight-line scanning. Biomed Eng Online 2018; 17:105. [PMID: 30075784 PMCID: PMC6076421 DOI: 10.1186/s12938-018-0537-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2018] [Accepted: 07/30/2018] [Indexed: 01/06/2023] Open
Abstract
Background For practical straight-line scanning in photoacoustic imaging (PAI), serious artifacts caused by missing data will occur. Traditional total variation (TV)-based algorithms fail to obtain satisfactory results, with an over-smoothed and blurred geometric structure. Therefore, it is important to develop a new algorithm to improve the quality of practical straight-line reconstructed images. Methods In this paper, a combined nonlocal patch and TV-based regularization model for PAI reconstruction is proposed to solve these problems. A modified adaptive nonlocal weight function is adopted to provide more reliable estimations for the similarities between patches. Similar patches are searched for throughout the entire image; thus, this model realizes adaptive search for the neighborhood of the patch. The optimization problem is simplified to a common iterative PAI reconstruction problem. Results and conclusion The proposed algorithm is validated by a series of numerical simulations and an in vitro experiment for straight-line scanning. The results of patch-TV are compared to those of two mainstream TV-based algorithms as well as the iterative algorithm only with patch-based regularization. Moreover, the peak signal-to-noise ratio, the noise robustness, and the convergence and calculation speeds are compared and discussed. The results show that the proposed patch-TV yields significant improvement over the other three algorithms qualitatively and quantitatively. These simulations and experiment indicate that the patch-TV algorithm successfully solves the problems of PAI reconstruction and is highly effective in practical PAI applications.
Collapse
Affiliation(s)
- Jin Wang
- Department of Electronic Engineering, Fudan University, No. 220 Handan Road, Shanghai, 200433, China
| | - Yuanyuan Wang
- Department of Electronic Engineering, Fudan University, No. 220 Handan Road, Shanghai, 200433, China. .,Key Laboratory of Medical Imaging Computing and Computer Assisted Intervention (MICCAI) of Shanghai, Fudan University, Shanghai, 200433, China.
| |
Collapse
|
10
|
Wang J, Zhang C, Wang Y. A photoacoustic imaging reconstruction method based on directional total variation with adaptive directivity. Biomed Eng Online 2017; 16:64. [PMID: 28558769 PMCID: PMC5450113 DOI: 10.1186/s12938-017-0366-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2017] [Accepted: 05/24/2017] [Indexed: 01/05/2023] Open
Abstract
BACKGROUND In photoacoustic tomography (PAT), total variation (TV) based iteration algorithm is reported to have a good performance in PAT image reconstruction. However, classical TV based algorithm fails to preserve the edges and texture details of the image because it is not sensitive to the direction of the image. Therefore, it is of great significance to develop a new PAT reconstruction algorithm to effectively solve the drawback of TV. METHODS In this paper, a directional total variation with adaptive directivity (DDTV) model-based PAT image reconstruction algorithm, which weightedly sums the image gradients based on the spatially varying directivity pattern of the image is proposed to overcome the shortcomings of TV. The orientation field of the image is adaptively estimated through a gradient-based approach. The image gradients are weighted at every pixel based on both its anisotropic direction and another parameter, which evaluates the estimated orientation field reliability. An efficient algorithm is derived to solve the iteration problem associated with DDTV and possessing directivity of the image adaptively updated for each iteration step. RESULTS AND CONCLUSION Several texture images with various directivity patterns are chosen as the phantoms for the numerical simulations. The 180-, 90- and 30-view circular scans are conducted. Results obtained show that the DDTV-based PAT reconstructed algorithm outperforms the filtered back-projection method (FBP) and TV algorithms in the quality of reconstructed images with the peak signal-to-noise rations (PSNR) exceeding those of TV and FBP by about 10 and 18 dB, respectively, for all cases. The Shepp-Logan phantom is studied with further discussion of multimode scanning, convergence speed, robustness and universality aspects. In-vitro experiments are performed for both the sparse-view circular scanning and linear scanning. The results further prove the effectiveness of the DDTV, which shows better results than that of the TV with sharper image edges and clearer texture details. Both numerical simulation and in vitro experiments confirm that the DDTV provides a significant quality improvement of PAT reconstructed images for various directivity patterns.
Collapse
Affiliation(s)
- Jin Wang
- Department of Electronic Engineering, Fudan University, Shanghai, 200433 China
| | - Chen Zhang
- Department of Electronic Engineering, Fudan University, Shanghai, 200433 China
| | - Yuanyuan Wang
- Department of Electronic Engineering, Fudan University, Shanghai, 200433 China
- Key Laboratory of Medical Imaging Computing and Computer-Assisted Intervention of Shanghai, Shanghai, 200433 China
| |
Collapse
|
11
|
An Efficient Compensation Method for Limited-View Photoacoustic Imaging Reconstruction Based on Gerchberg–Papoulis Extrapolation. APPLIED SCIENCES-BASEL 2017. [DOI: 10.3390/app7050505] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
12
|
Deán-Ben XL, Razansky D. On the link between the speckle free nature of optoacoustics and visibility of structures in limited-view tomography. PHOTOACOUSTICS 2016; 4:133-140. [PMID: 28066714 PMCID: PMC5200938 DOI: 10.1016/j.pacs.2016.10.001] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Revised: 09/30/2016] [Accepted: 10/13/2016] [Indexed: 05/08/2023]
Abstract
Similar to pulse-echo ultrasound, optoacoustic imaging encodes the location of optical absorbers by the time-of-flight of ultrasound waves. Yet, signal generation mechanisms are fundamentally different for the two modalities, leading to significant distinction between the optimum image formation strategies. While interference of back-scattered ultrasound waves with random phases causes speckle noise in ultrasound images, speckle formation is hindered by the strong correlation between the optoacoustic responses corresponding to individual sources. However, visibility of structures is severely hampered when attempting to acquire optoacoustic images under limited-view tomographic geometries. In this tutorial article, we systematically describe the basic principles of optoacoustic signal generation and image formation for objects ranging from individual sub-resolution absorbers to a continuous absorption distribution. The results are of relevance for the proper interpretation of optoacoustic images acquired under limited-view scenarios and may also serve as a basis for optimal design of tomographic acquisition geometries and image formation strategies.
Collapse
|
13
|
Pulkkinen A, Cox BT, Arridge SR, Goh H, Kaipio JP, Tarvainen T. Direct Estimation of Optical Parameters From Photoacoustic Time Series in Quantitative Photoacoustic Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2016; 35:2497-2508. [PMID: 27323361 DOI: 10.1109/tmi.2016.2581211] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Estimation of optical absorption and scattering of a target is an inverse problem associated with quantitative photoacoustic tomography. Conventionally, the problem is expressed as two folded. First, images of initial pressure distribution created by absorption of a light pulse are formed based on acoustic boundary measurements. Then, the optical properties are determined based on these photoacoustic images. The optical stage of the inverse problem can thus suffer from, for example, artefacts caused by the acoustic stage. These could be caused by imperfections in the acoustic measurement setting, of which an example is a limited view acoustic measurement geometry. In this work, the forward model of quantitative photoacoustic tomography is treated as a coupled acoustic and optical model and the inverse problem is solved by using a Bayesian approach. Spatial distribution of the optical properties of the imaged target are estimated directly from the photoacoustic time series in varying acoustic detection and optical illumination configurations. It is numerically demonstrated, that estimation of optical properties of the imaged target is feasible in limited view acoustic detection setting.
Collapse
|
14
|
Efficient block-sparse model-based algorithm for photoacoustic image reconstruction. Biomed Signal Process Control 2016. [DOI: 10.1016/j.bspc.2015.12.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
15
|
Zhang C, Zhang Y, Wang Y. A photoacoustic image reconstruction method using total variation and nonconvex optimization. Biomed Eng Online 2014; 13:117. [PMID: 25129644 PMCID: PMC4148921 DOI: 10.1186/1475-925x-13-117] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2014] [Accepted: 08/12/2014] [Indexed: 11/25/2022] Open
Abstract
Background In photoacoustic imaging (PAI), the reduction of scanning time is a major concern for PAI in practice. A popular strategy is to reconstruct the image from the sparse-view sampling data. However, the insufficient data leads to reconstruction quality deteriorating. Therefore, it is very important to enhance the quality of the sparse-view reconstructed images. Method In this paper, we proposed a joint total variation and Lp-norm (TV-Lp) based image reconstruction algorithm for PAI. In this algorithm, the reconstructed image is updated by calculating its total variation value and Lp-norm value. Along with the iteration, an operator-splitting framework is utilized to reduce the computational cost and the Barzilai-Borwein step size selection method is adopted to obtain the faster convergence. Results and conclusion Through the numerical simulation, the proposed algorithm is validated and compared with other widely used PAI reconstruction algorithms. It is revealed in the simulation result that the proposed algorithm may be more accurate than the other algorithms. Moreover, the computational cost, the convergence, the robustness to noises and the tunable parameters of the algorithm are all discussed respectively. We also implement the TV-Lp algorithm in the in-vitro experiments to verify its performance in practice. Through the numerical simulations and in-vitro experiments, it is demonstrated that the proposed algorithm enhances the quality of the reconstructed images with faster calculation speed and convergence.
Collapse
Affiliation(s)
| | | | - Yuanyuan Wang
- Department of Electronic Engineering, Fudan University, Shanghai 200433, China.
| |
Collapse
|
16
|
Zhang Y, Wang Y, Zhang C. Total variation based gradient descent algorithm for sparse-view photoacoustic image reconstruction. ULTRASONICS 2012; 52:1046-55. [PMID: 22986153 DOI: 10.1016/j.ultras.2012.08.012] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2012] [Revised: 08/20/2012] [Accepted: 08/20/2012] [Indexed: 05/06/2023]
Abstract
In photoacoustic imaging (PAI), reconstruction from sparse-view sampling data is a remaining challenge in the cases of fast or real-time imaging. In this paper, we present our study on a total variation based gradient descent (TV-GD) algorithm for sparse-view PAI reconstruction. This algorithm involves the total variation (TV) method in compressed sensing (CS) theory. The objective function of the algorithm is modified by adding the TV value of the reconstructed image. With this modification, the reconstructed image could be closer to the real optical energy distribution map. Additionally in the proposed algorithm, the photoacoustic data is processed and the image is updated individually at each detection point. In this way, the calculation with large matrix can be avoided and a more frequent image update can be obtained. Through the numerical simulations, the proposed algorithm is verified and compared with other reconstruction algorithms which have been widely used in PAI. The peak signal-to-noise ratio (PSNR) of the image reconstructed by this algorithm is higher than those by the other algorithms. Additionally, the convergence of the algorithm, the robustness to noise and the tunable parameter are further discussed. The TV-based algorithm is also implemented in the in vitro experiment. The better performance of the proposed method is revealed in the experiments results. From the results, it is seen that the TV-GD algorithm may be a practical and efficient algorithm for sparse-view PAI reconstruction.
Collapse
Affiliation(s)
- Yan Zhang
- Department of Electronics Engineering, Fudan University, Shanghai 200437, China
| | | | | |
Collapse
|
17
|
Wang X, Roberts WW, Carson PL, Wood DP, Fowlkes JB. Photoacoustic tomography: a potential new tool for prostate cancer. BIOMEDICAL OPTICS EXPRESS 2010; 1:1117-1126. [PMID: 21258534 PMCID: PMC3018085 DOI: 10.1364/boe.1.001117] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2010] [Revised: 10/01/2010] [Accepted: 10/05/2010] [Indexed: 05/18/2023]
Abstract
The feasibility of photoacoustic tomography (PAT) for noninvasive imaging of prostate cancer was explored through the study on a canine model in vivo. Imaging of blood-rich lesions mimicking prostate tumors was achieved using a commercial medical ultrasound (US) system without affecting its original imaging functions. Based on the optical contrast between hemoglobin and other tissues, PAT has demonstrated good sensitivity and high contrast-to-noise ratio in visualizing deep lesions; while US has presented the morphological features including the boundary and the urethral of the prostate. PAT of prostate cancer may facilitate improved tumor localization, staging of disease, and detection of recurrences.
Collapse
Affiliation(s)
- Xueding Wang
- Department of Radiology, University of Michigan School of Medicine, Ann Arbor, MI 48109, USA
| | - William W. Roberts
- Department of Urology, University of Michigan School of Medicine, Ann Arbor, MI 48109, USA
| | - Paul L. Carson
- Department of Radiology, University of Michigan School of Medicine, Ann Arbor, MI 48109, USA
| | - David P. Wood
- Department of Urology, University of Michigan School of Medicine, Ann Arbor, MI 48109, USA
| | - J. Brian Fowlkes
- Department of Radiology, University of Michigan School of Medicine, Ann Arbor, MI 48109, USA
| |
Collapse
|
18
|
Zhang C, Li C, Wang LV. Fast and robust deconvolution-based image reconstruction for photoacoustic tomography in circular geometry: experimental validation. IEEE PHOTONICS JOURNAL 2010; 2:57-66. [PMID: 20953277 PMCID: PMC2954598 DOI: 10.1109/jphot.2010.2042801] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Photoacoustic tomography (PAT) is a fast-developing biomedical imaging technology suitable for in vivo imaging. PAT in spherical or circular geometry gives good image resolution yet is slow or expensive in signal acquisition and image formation. Reducing the number of detection angles can ameliorate such issues usually at the expense of image quality. This paper introduces a deconvolution-based algorithm that models the imaging process as a linear and shift-invariant system. As demonstrated by the in vivo experiment, this algorithm not only runs much faster than the back-projection algorithm, but also shows stronger robustness in that it provides better image quality when detection angles are sparse. Therefore, this algorithm promises to enable real-time PAT in circular geometry.
Collapse
Affiliation(s)
- Chi Zhang
- Optical Imaging Laboratory, Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130-4899, USA
| | | | | |
Collapse
|