OSA's Digital Library

Journal of the Optical Society of America A

Journal of the Optical Society of America A

| OPTICS, IMAGE SCIENCE, AND VISION

  • Vol. 15, Iss. 10 — Oct. 1, 1998
  • pp: 2639–2650

Total least-squares reconstruction with wavelets for optical tomography

Wenwu Zhu, Yao Wang, and Jun Zhang  »View Author Affiliations


JOSA A, Vol. 15, Issue 10, pp. 2639-2650 (1998)
http://dx.doi.org/10.1364/JOSAA.15.002639


View Full Text Article

Enhanced HTML    Acrobat PDF (1899 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

In a previous paper [Zhu, J. Opt. Soc. Am. A 14, 799 (1997)] an iterative algorithm for obtaining the total least-squares (TLS) solution of a linear system based on the Rayleigh quotient formulation was presented. Here we derive what to our knowledge are the first statistical properties of this solution. It is shown that the Rayleigh-quotient-form TLS (RQF-TLS) estimator is equivalent to the maximum-likelihood estimator when noise terms in both data and operator elements are independent and identically distributed Gaussian. A perturbation analysis of the RQF-TLS solution is derived, and from it the mean square error of the RQF-TLS solution is obtained in closed form, which is valid at small noise levels. We then present a wavelet-based multiresolution scheme for obtaining the TLS solution. This method was employed with a multigrid algorithm to solve the linear perturbation equation encountered in optical tomography. Results from numerical simulations show that this method requires substantially less computation than the previously reported one-grid TLS algorithm. The method also allows one to identify regions of interest quickly from a coarse-level reconstruction and restrict the reconstruction in the following fine resolutions to those regions. Finally, the method is less sensitive to noise than the one-grid TLS and multigrid least-squares algorithms.

© 1998 Optical Society of America

OCIS Codes
(100.0100) Image processing : Image processing
(100.2000) Image processing : Digital image processing
(100.3010) Image processing : Image reconstruction techniques
(100.3190) Image processing : Inverse problems
(100.6950) Image processing : Tomographic image processing
(100.7410) Image processing : Wavelets

History
Original Manuscript: March 19, 1998
Revised Manuscript: June 10, 1998
Manuscript Accepted: May 11, 1998
Published: October 1, 1998

Citation
Wenwu Zhu, Yao Wang, and Jun Zhang, "Total least-squares reconstruction with wavelets for optical tomography," J. Opt. Soc. Am. A 15, 2639-2650 (1998)
http://www.opticsinfobase.org/josaa/abstract.cfm?URI=josaa-15-10-2639


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. R. L. Barbour, H. L. Graber, Y. Wang, J. Chang, R. Aronson, “A perturbation approach for optical diffusion tomography using continuous-wave and time-resolved data,” in Medical Optical Tomography: Functional Imaging and Monitoring, G. Mueller, B. Chance, R. Alfano, S. Arridge, J. Beuthan, E. Gratton, M. Kaschke, B. Masters, S. Svanberg, P. van der Zee, eds., Vol. IS11 of Institute Series (Society of Photo-optical Instrumentation Engineers, Bellingham, Washington, 1993), pp. 87–120.
  2. Y. Wang, J. Chang, R. Aronson, R. L. Barbour, H. L. Graber, J. Lubowsky, “Imaging scattering media by diffusion tomography: an iterative perturbation approach,” in Physiological Monitoring and Early Detection Diagnostic Methods, T. S. Mang, ed., Proc. SPIE1641, 58–71 (1992). [CrossRef]
  3. Y. Q. Yao, Y. Wang, Y. L. Pei, W. Zhu, J. H. Hu, R. L. Barbour, “Frequency domain optical tomography in human tissue,” in Experimental and Numerical Methods for Solving Ill-Posed Inverse Problems: Medical and Nonmedical Application, R. L. Barbour, M. J. Carvlin, M. A. Fiddy eds., Proc. SPIE2570, 254–266 (1995). [CrossRef]
  4. J. R. Singer, F. A. Grünbaum, P. Kohn, J. P. Zubelli, “Image reconstruction of the interior of bodies that diffuse radiation,” Science 248, 990–993 (1990). [CrossRef] [PubMed]
  5. S. R. Arridge, P. van der Zee, M. Cope, D. Delpy, “New results for the development of infrared-red absorption imaging,” in Biomedical Image Processing, A. C. Bovik, W. E. Higgins, eds., Proc. SPIE1245, 92–103 (1990). [CrossRef]
  6. S. R. Arridge, “The forward and inverse problems in time resolved infra-red imaging,” in Medical Optical Tomography: Functional Imaging and Monitoring, G. Mueller, B. Chance, R. Alfano, S. Arridge, J. Beuthan, E. Gratton, M. Kaschke, B. Masters, S. Svanberg, P. van der Zee, eds., Vol. IS11 of Institute Series (Society of Photo-optical Instrumentation Engineers, Bellingham, Washington, 1993), pp. 35–64.
  7. W. Zhu, Y. Wang, Y. Deng, Y. Yao, R. Barbour, “A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography,” IEEE Trans. Med. Imag. 16, 210–217 (1997). [CrossRef]
  8. W. Zhu, Y. Wang, Y. Yao, J. Chang, H. L. Graber, R. Barbour, “Iterative total least-squares image reconstruction algorithm for optical tomography by the conjugate gradient method,” J. Opt. Soc. Am. A 14, 799–807 (1997). [CrossRef]
  9. J. Chang, Y. Wang, R. Aronson, H. L. Graber, R. L. Barbour, “A layer-stripping approach for recovery of scattering media from time-resolved data,” in Inverse Problems in Scattering and Imaging, M. A. Fiddy, ed., Proc. SPIE1767, 384–395 (1992). [CrossRef]
  10. J. Chang, W. Zhu, Y. Wang, H. L. Graber, R. L. Barbour, “Regularized progressive expansion algorithm for recovery of scattering media from time-resolved data,” J. Opt. Soc. Am. A 14, 306–312 (1997). [CrossRef]
  11. G. H. Golub, C. F. Van Loan, “An analysis of the total least squares problem,” SIAM (Soc. Ind. Appl. Math.) J. Math. Anal. 17, 883–893 (1980).
  12. V. Z. Measarvić, N. P. Galatsanos, A. Katsaggelos, “Regularized constrained total least squares image restoration,” IEEE Trans. Image Process. 4, 1096–1108 (1995). [CrossRef]
  13. P. Li, S. W. Flax, E. S. Ebbini, M. O’Donnell, “Blocked element compensation in phased array imaging,” IEEE Trans. Ultrason. Frequencies Freq. Control 40, 282–292 (1993).
  14. G. Wang, J. Zhang, G. W. Pan, “Solution of inverse problems in image processing by wavelet expansion,” IEEE Trans. Image Process. 4, 579–593 (1995). [CrossRef] [PubMed]
  15. L. Blanc-Feraud, P. Charbonnier, P. Lobel, M. Barlaud, “A fast tomographic reconstruction algorithm in the 2D wavelet transform domain,” in Proceedings of the 1994 IEEE International Conference on Acoustics, Speech, and Signal Processing (Institute of Electrical and Electronics Engineers, Piscataway, N.J., 1994), pp. 305–308.
  16. A. H. Delaney, Y. Bresler, “Multiresolution tomographic reconstruction using wavelets,” IEEE Trans. Image Process. 4, 799–813 (1995). [CrossRef] [PubMed]
  17. M. R. Banham, Nikolas P. Galatsanos, H. L. Gonzalez, A. K. Katsaggelos, “Multichannel restoration of single channel images using a wavelet-based subband decomposition,” IEEE Trans. Image Process. 3, 821–833 (1994). [CrossRef] [PubMed]
  18. M. T. Silvia, E. C. Tacker, “Regularization of Marchenko’s integral equation by total least squares,” J. Acoust. Soc. Am. 72, 1202–1207 (1982). [CrossRef]
  19. G. H. Golub, “Some modified matrix eigenvalue problems,” SIAM (Soc. Ind. Appl. Math.) Rev. 15, 318–334 (1973).
  20. H. Chen, T. K. Sarkar, S. A. Dianat, J. D. Brule, “Adaptive spectral estimation by the conjugate gradient method,” IEEE Trans. Acoust. Speech Signal Process. 34, 272–284 (1986). [CrossRef]
  21. X. Yang, T. K. Sarkar, E. Arvas, “A survey of conjugate gradient algorithms for solution of extreme eigen-problem of a symmetric matrix,” IEEE Trans. Acoust. Speech Signal Process. 37, 1550–1556 (1989). [CrossRef]
  22. S. Van Huffel, J. Vandewalle, The Total Least Squares Problem: Computational Aspects and Analysis (Society for Industrial and Applied Mathematics, Philadelphia, Pa., 1991).
  23. T. J. Abatzoglou, J. M. Mendel, G. A. Harada, “The constrained total least squares technique and its applications to harmonic superresolution,” IEEE Trans. Signal Process. 39, 1070–1087 (1991). [CrossRef]
  24. V. Z. Measarvić, N. P. Galatsanos, “map and regularized constrained total least squares image restoration,” in Proceedings of the 1994 IEEE International Conference on Image Processing (Institute of Electrical and Electronics Engineers, Piscatway, N.J., 1994), Vol. 3, pp. 177–181.
  25. S. G. Mallat, “A theory for multiresolution signal decomposition: the wavelet representation,” IEEE Trans. Pattern. Anal. Mach. Intell. 11, 674–693 (1989). [CrossRef]
  26. W. Hackbusch, Multigrid Methods and Applications (Springer-Verlag, Berlin, 1985).
  27. W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery, Numerical Recipes in C: The Art of Scientific Computing, 2nd ed. (Cambridge U. Press, Cambridge, UK, 1992).
  28. Y. Q. Yao, Y. Wang, R. L. Barbour, H. L. Graber, J. W. Chang, “Scattering characteristics of photon density waves from an object in a spherically two-layer medium,” in Optical Tomography, Photon Migration, and Spectroscopy of Tissue and Model Media: Theory, Human Studies, and Instrumentation, B. Chance, R. R. Alfano, eds., Proc. SPIE2389, 291–303 (1995). [CrossRef]
  29. W. L. Briggs, A Multigrid Tutorial (Society for Industrial and Applied Mathematics, Philadelphia, Pa., 1987).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited