Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter September 1, 2020

Sampling in thermoacoustic tomography

Chase Mathison


We explore the effect of sampling rates when measuring data given by Mf for special operators M arising in Thermoacoustic Tomography. We start with sampling requirements on Mf given f satisfying certain conditions. After this we discuss the resolution limit on f posed by the sampling rate of Mf without assuming any conditions on these sampling rates. Next we discuss aliasing artifacts when Mf is known to be under sampled in one or more of its variables. Finally, we discuss averaging of measurement data and resulting aliasing and artifacts, along with a scheme for anti-aliasing.

Funding source: National Science Foundation

Award Identifier / Grant number: DMS-1600327

Funding statement: The author was partly supported by NSF Grant DMS-1600327.


The author would like to thank Dr. Plamen Stefanov for suggesting this problem and for his guidance in the analysis of this problem.


[1] S. Bartels, Numerical Approximation of Partial Differential Equations, Texts Appl. Math. 64, Springer, Cham, 2016. 10.1007/978-3-319-32354-1Search in Google Scholar

[2] P. Burgholzer, C. Hofer, G. J. Matt, G. Paltauf, M. Haltmeier and O. Scherzer, Thermoacoustic tomography using a fiber-based fabry-perot interferometer as an integrating line detector, Proc. SPIE 6086 (2006), 434–442. 10.1117/12.644047Search in Google Scholar

[3] R. Courant, K. Friedrichs and H. Lewy, On the partial difference equations of mathematical physics, Technical report, 1959. Search in Google Scholar

[4] M. J. Grote and I. Sim, Efficient pml for the wave equation, preprint (2010), Search in Google Scholar

[5] H. Grün, G. Paltauf, M. Haltmeier and P. Burgholzer, Photoacoustic tomography using a fiber based fabry-perot interferometer as an integrating line detector and image reconstruction by model-based time reversal method, Novel Optical Instrumentation for Biomedical Applications. III, Optical Society of America, Washington (2007), 6631–6636. 10.1364/ECBO.2007.6631_6Search in Google Scholar

[6] V. Guillemin and S. Sternberg, Semi-Classical Analysis, International Press, Boston, 2013. Search in Google Scholar

[7] M. Haltmeier, Sampling conditions for the circular Radon transform, IEEE Trans. Image Process. 25 (2016), no. 6, 2910–2919. 10.1109/TIP.2016.2551364Search in Google Scholar PubMed

[8] M. Haltmeier, O. Scherzer, P. Burgholzer, R. Nuster and G. Paltauf, Thermoacoustic tomography and the circular Radon transform: Exact inversion formula, Math. Models Methods Appl. Sci. 17 (2007), no. 4, 635–655. 10.1142/S0218202507002054Search in Google Scholar

[9] M. Haltmeier, O. Scherzer, P. Burgholzer and G. Paltauf, Thermoacoustic computed tomography with large planar receivers, Inverse Problems 20 (2004), no. 5, 1663–1673. 10.1088/0266-5611/20/5/021Search in Google Scholar

[10] Y. Hristova, Time reversal in thermoacoustic tomography—an error estimate, Inverse Problems 25 (2009), no. 5, Article ID 055008. 10.1088/0266-5611/25/5/055008Search in Google Scholar

[11] Y. Hristova, P. Kuchment and L. Nguyen, Reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media, Inverse Problems 24 (2008), no. 5, Article ID 055006. 10.1088/0266-5611/24/5/055006Search in Google Scholar

[12] R. A. Kruger, W. L. Kiser, K. D. Miller, H. E. Reynolds, D. R. Reinecke, G. A. Kruger and P. J. Hofacker, Thermoacoustic ct, 2000 IEEE MTT-S International Microwave Symposium Digest (Cat. No. 00CH37017). Volume 2, IEEE Press, Piscataway (2000), 933–936. 10.1109/MWSYM.2000.863510Search in Google Scholar

[13] R. A. Kruger, D. R. Reinecke and G. A. Kruger, Thermoacoustic computed tomography-technical considerations, Med. Phys. 26 (1999), no. 9, 1832–1837. 10.1118/1.598688Search in Google Scholar PubMed

[14] P. Kuchment, The Radon Transform and Medical Imaging, CBMS-NSF Regional Conf. Ser. in Appl. Math. 85, SIAM, Philadelphia, 2014. 10.1137/1.9781611973297Search in Google Scholar

[15] L. A. Kunyansky, Thermoacoustic tomography with detectors on an open curve: An efficient reconstruction algorithm, Inverse Problems 24 (2008), no. 5, Article ID 055021. 10.1088/0266-5611/24/5/055021Search in Google Scholar

[16] R. J. Marks, II, Introduction to Shannon Sampling and Interpolation Theory, Springer Texts Elect. Eng., Springer, New York, 1991. 10.1007/978-1-4613-9708-3Search in Google Scholar

[17] A. Martinez, An Introduction to Semiclassical and Microlocal Analysis, Universitext, Springer, New York, 2002. 10.1007/978-1-4757-4495-8Search in Google Scholar

[18] A. A. Oraevsky, S. L. Jacques, R. O. Esenaliev and F. K. Tittel, Laser-based optoacoustic imaging in biological tissues, Laser-Tissue Interaction V; and Ultraviolet Radiation Hazards 2134, International Society for Optics and Photonics, Bellingham (1994), 122–129. 10.1117/12.182927Search in Google Scholar

[19] D. P. Petersen and D. Middleton, Sampling and reconstruction of wave-number-limited functions in N-dimensional Euclidean spaces, Information Control 5 (1962), 279–323. 10.1016/S0019-9958(62)90633-2Search in Google Scholar

[20] P. Stefanov, Semiclassical sampling and discretization of certain linear inverse problems, preprint (2018), 10.1137/19M123868XSearch in Google Scholar

[21] P. Stefanov and G. Uhlmann, Thermoacoustic tomography with variable sound speed, Inverse Problems 25 (2009), no. 7, Article ID 075011. 10.1088/0266-5611/25/7/075011Search in Google Scholar

[22] P. Stefanov and Y. Yang, Thermo- and photoacoustic tomography with variable speed and planar detectors, SIAM J. Math. Anal. 49 (2017), no. 1, 297–310. 10.1137/16M1073716Search in Google Scholar

[23] G. Strang, The wave equation and staggered leapfrog, 2006. Search in Google Scholar

[24] M. E. Taylor, Geometrical Optics and Fourier Integral Operators, Princeton University, Princeton, 1981. Search in Google Scholar

[25] K. Wang and M. A. Anastasio, Photoacoustic and thermoacoustic tomography: Image formation principles, Handbook of Mathematical Methods in Imaging. Vol. 1, 2, 3, Springer, New York (2015), 1081–1116. 10.1007/978-1-4939-0790-8_50Search in Google Scholar

[26] G. Zangerl, O. Scherzer and M. Haltmeier, Circular integrating detectors in photo and thermoacoustic tomography, Inverse Probl. Sci. Eng. 17 (2009), no. 1, 133–142. 10.1080/17415970802166782Search in Google Scholar

[27] M. Zworski, Semiclassical Analysis, Grad. Stud. Math. 138, American Mathematical Society, Providence, 2012. 10.1090/gsm/138Search in Google Scholar

Received: 2020-01-08
Revised: 2020-06-09
Accepted: 2020-07-29
Published Online: 2020-09-01
Published in Print: 2020-12-01

© 2020 Walter de Gruyter GmbH, Berlin/Boston