A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography

Wenwu Zhu, Yao Wang, Yining Deng, Yuqi Yao, Randall L. Barbour

Research output: Contribution to journalArticle

Abstract

In this paper we present a wavelet-based multigrid approach to solve the perturbation equation encountered in optical tomography. With this scheme the unknown image the data as well as the weight matrix are all represented by wavelet expansions thus yielding a multiresolution representation of the original perturbation equation in the wavelet domain. This transformed equation is then solved using a multigrid scheme by which an increasing portion of wavelet coefficients of the unknown image are solved in successive approximations. One can also quickly identify regions of interest (ROI's) from a coarse level reconstruction and restrict the reconstruction in the following fine resolutions to those regions. At each resolution level a regularized least squares solution is obtained using the conjugate gradient descent method. This approach has been applied to continuous wave data calculated based on the diffusion approximation of several two-dimensional (2-D) test media. Compared to a previously reported one grid algorithm the multigrid method requires substantially shorter computation time under the same reconstruction quality criterion.

Original languageEnglish (US)
Pages (from-to)210-217
Number of pages8
JournalIEEE Transactions on Medical Imaging
Volume16
Issue number2
StatePublished - 1997

Fingerprint

Optical Tomography
Optical tomography
Least-Squares Analysis
Weights and Measures

Keywords

  • Image reconstruction
  • Multigrid method
  • Optical tomography
  • Wavelet transform

ASJC Scopus subject areas

  • Biomedical Engineering
  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology
  • Electrical and Electronic Engineering
  • Computer Science Applications
  • Computational Theory and Mathematics

Cite this

A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography. / Zhu, Wenwu; Wang, Yao; Deng, Yining; Yao, Yuqi; Barbour, Randall L.

In: IEEE Transactions on Medical Imaging, Vol. 16, No. 2, 1997, p. 210-217.

Research output: Contribution to journalArticle

Zhu, Wenwu ; Wang, Yao ; Deng, Yining ; Yao, Yuqi ; Barbour, Randall L. / A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography. In: IEEE Transactions on Medical Imaging. 1997 ; Vol. 16, No. 2. pp. 210-217.
@article{28441a8924f245e0b3b9afabd6506b05,
title = "A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography",
abstract = "In this paper we present a wavelet-based multigrid approach to solve the perturbation equation encountered in optical tomography. With this scheme the unknown image the data as well as the weight matrix are all represented by wavelet expansions thus yielding a multiresolution representation of the original perturbation equation in the wavelet domain. This transformed equation is then solved using a multigrid scheme by which an increasing portion of wavelet coefficients of the unknown image are solved in successive approximations. One can also quickly identify regions of interest (ROI's) from a coarse level reconstruction and restrict the reconstruction in the following fine resolutions to those regions. At each resolution level a regularized least squares solution is obtained using the conjugate gradient descent method. This approach has been applied to continuous wave data calculated based on the diffusion approximation of several two-dimensional (2-D) test media. Compared to a previously reported one grid algorithm the multigrid method requires substantially shorter computation time under the same reconstruction quality criterion.",
keywords = "Image reconstruction, Multigrid method, Optical tomography, Wavelet transform",
author = "Wenwu Zhu and Yao Wang and Yining Deng and Yuqi Yao and Barbour, {Randall L.}",
year = "1997",
language = "English (US)",
volume = "16",
pages = "210--217",
journal = "IEEE Transactions on Medical Imaging",
issn = "0278-0062",
publisher = "Institute of Electrical and Electronics Engineers Inc.",
number = "2",

}

TY - JOUR

T1 - A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography

AU - Zhu, Wenwu

AU - Wang, Yao

AU - Deng, Yining

AU - Yao, Yuqi

AU - Barbour, Randall L.

PY - 1997

Y1 - 1997

N2 - In this paper we present a wavelet-based multigrid approach to solve the perturbation equation encountered in optical tomography. With this scheme the unknown image the data as well as the weight matrix are all represented by wavelet expansions thus yielding a multiresolution representation of the original perturbation equation in the wavelet domain. This transformed equation is then solved using a multigrid scheme by which an increasing portion of wavelet coefficients of the unknown image are solved in successive approximations. One can also quickly identify regions of interest (ROI's) from a coarse level reconstruction and restrict the reconstruction in the following fine resolutions to those regions. At each resolution level a regularized least squares solution is obtained using the conjugate gradient descent method. This approach has been applied to continuous wave data calculated based on the diffusion approximation of several two-dimensional (2-D) test media. Compared to a previously reported one grid algorithm the multigrid method requires substantially shorter computation time under the same reconstruction quality criterion.

AB - In this paper we present a wavelet-based multigrid approach to solve the perturbation equation encountered in optical tomography. With this scheme the unknown image the data as well as the weight matrix are all represented by wavelet expansions thus yielding a multiresolution representation of the original perturbation equation in the wavelet domain. This transformed equation is then solved using a multigrid scheme by which an increasing portion of wavelet coefficients of the unknown image are solved in successive approximations. One can also quickly identify regions of interest (ROI's) from a coarse level reconstruction and restrict the reconstruction in the following fine resolutions to those regions. At each resolution level a regularized least squares solution is obtained using the conjugate gradient descent method. This approach has been applied to continuous wave data calculated based on the diffusion approximation of several two-dimensional (2-D) test media. Compared to a previously reported one grid algorithm the multigrid method requires substantially shorter computation time under the same reconstruction quality criterion.

KW - Image reconstruction

KW - Multigrid method

KW - Optical tomography

KW - Wavelet transform

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

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

M3 - Article

VL - 16

SP - 210

EP - 217

JO - IEEE Transactions on Medical Imaging

JF - IEEE Transactions on Medical Imaging

SN - 0278-0062

IS - 2

ER -