Total least-squares reconstruction with wavelets for optical tomography

Wenwu Zhu, Yao Wang, Jun Zhang

Research output: Contribution to journalArticle

Abstract

In a previous paper [Zhu et al., 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.

Original languageEnglish (US)
Pages (from-to)2639-2650
Number of pages12
JournalJournal of the Optical Society of America A: Optics and Image Science, and Vision
Volume15
Issue number10
StatePublished - 1998

Fingerprint

Optical Tomography
Optical tomography
Least-Squares Analysis
Noise
Mean square error
Maximum likelihood
Linear systems
Computer simulation

ASJC Scopus subject areas

  • Electronic, Optical and Magnetic Materials
  • Computer Vision and Pattern Recognition

Cite this

Total least-squares reconstruction with wavelets for optical tomography. / Zhu, Wenwu; Wang, Yao; Zhang, Jun.

In: Journal of the Optical Society of America A: Optics and Image Science, and Vision, Vol. 15, No. 10, 1998, p. 2639-2650.

Research output: Contribution to journalArticle

@article{e76d49869be14c4d90f934eeb3bdc1c8,
title = "Total least-squares reconstruction with wavelets for optical tomography",
abstract = "In a previous paper [Zhu et al., 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.",
author = "Wenwu Zhu and Yao Wang and Jun Zhang",
year = "1998",
language = "English (US)",
volume = "15",
pages = "2639--2650",
journal = "Journal of the Optical Society of America A: Optics and Image Science, and Vision",
issn = "0740-3232",
publisher = "The Optical Society",
number = "10",

}

TY - JOUR

T1 - Total least-squares reconstruction with wavelets for optical tomography

AU - Zhu, Wenwu

AU - Wang, Yao

AU - Zhang, Jun

PY - 1998

Y1 - 1998

N2 - In a previous paper [Zhu et al., 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.

AB - In a previous paper [Zhu et al., 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.

UR - http://www.scopus.com/inward/record.url?scp=0032186422&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0032186422&partnerID=8YFLogxK

M3 - Article

VL - 15

SP - 2639

EP - 2650

JO - Journal of the Optical Society of America A: Optics and Image Science, and Vision

JF - Journal of the Optical Society of America A: Optics and Image Science, and Vision

SN - 0740-3232

IS - 10

ER -