SPECT reconstruction using DCT-induced tight framelet regularization

Jiahan Zhang, Si Li, Yuesheng Xu, C. R. Schmidtlein, Edward D. Lipson, David H. Feiglin, Andrzej Krol

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Wavelet transforms have been successfully applied in many fields of image processing. Yet, to our knowledge, they have never been directly incorporated to the objective function in Emission Computed Tomography (ECT) image reconstruction. Our aim has been to investigate if the 1-norm of non-decimated discrete cosine transform (DCT) coefficients of the estimated radiotracer distribution could be effectively used as the regularization term for the penalized-likelihood (PL) reconstruction, where a regularizer is used to enforce the image smoothness in the reconstruction. In this study, the 1-norm of 2D DCT wavelet decomposition was used as a regularization term. The Preconditioned Alternating Projection Algorithm (PAPA), which we proposed in earlier work to solve penalized likelihood (PL) reconstruction with non-differentiable regularizers, was used to solve this optimization problem. The DCT wavelet decompositions were performed on the transaxial reconstructed images. We reconstructed Monte Carlo simulated SPECT data obtained for a numerical phantom with Gaussian blobs as hot lesions and with a warm random lumpy background. Reconstructed images using the proposed method exhibited better noise suppression and improved lesion conspicuity, compared with images reconstructed using expectation maximization (EM) algorithm with Gaussian post filter (GPF). Also, the mean square error (MSE) was smaller, compared with EM-GPF. A critical and challenging aspect of this method was selection of optimal parameters. In summary, our numerical experiments demonstrated that the 1-norm of discrete cosine transform (DCT) wavelet frame transform DCT regularizer shows promise for SPECT image reconstruction using PAPA method.

Original languageEnglish (US)
Title of host publicationProgress in Biomedical Optics and Imaging - Proceedings of SPIE
PublisherSPIE
Volume9412
ISBN (Print)9781628415025
DOIs
StatePublished - 2015
EventMedical Imaging 2015: Physics of Medical Imaging - Orlando, United States
Duration: Feb 22 2015Feb 25 2015

Other

OtherMedical Imaging 2015: Physics of Medical Imaging
CountryUnited States
CityOrlando
Period2/22/152/25/15

Fingerprint

Wavelet Analysis
discrete cosine transform
Discrete cosine transforms
Single-Photon Emission-Computed Tomography
Computer-Assisted Image Processing
norms
Wavelet decomposition
image reconstruction
Image reconstruction
lesions
Emission-Computed Tomography
projection
decomposition
filters
Noise
wavelet analysis
Mean square error
Wavelet transforms
Tomography
image processing

Keywords

  • DCT
  • Framelet
  • Image reconstruction
  • Noise reduction
  • Re gularization

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics
  • Electronic, Optical and Magnetic Materials
  • Biomaterials
  • Radiology Nuclear Medicine and imaging

Cite this

Zhang, J., Li, S., Xu, Y., Schmidtlein, C. R., Lipson, E. D., Feiglin, D. H., & Krol, A. (2015). SPECT reconstruction using DCT-induced tight framelet regularization. In Progress in Biomedical Optics and Imaging - Proceedings of SPIE (Vol. 9412). [94123H] SPIE. https://doi.org/10.1117/12.2082118

SPECT reconstruction using DCT-induced tight framelet regularization. / Zhang, Jiahan; Li, Si; Xu, Yuesheng; Schmidtlein, C. R.; Lipson, Edward D.; Feiglin, David H.; Krol, Andrzej.

Progress in Biomedical Optics and Imaging - Proceedings of SPIE. Vol. 9412 SPIE, 2015. 94123H.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Zhang, J, Li, S, Xu, Y, Schmidtlein, CR, Lipson, ED, Feiglin, DH & Krol, A 2015, SPECT reconstruction using DCT-induced tight framelet regularization. in Progress in Biomedical Optics and Imaging - Proceedings of SPIE. vol. 9412, 94123H, SPIE, Medical Imaging 2015: Physics of Medical Imaging, Orlando, United States, 2/22/15. https://doi.org/10.1117/12.2082118
Zhang J, Li S, Xu Y, Schmidtlein CR, Lipson ED, Feiglin DH et al. SPECT reconstruction using DCT-induced tight framelet regularization. In Progress in Biomedical Optics and Imaging - Proceedings of SPIE. Vol. 9412. SPIE. 2015. 94123H https://doi.org/10.1117/12.2082118
Zhang, Jiahan ; Li, Si ; Xu, Yuesheng ; Schmidtlein, C. R. ; Lipson, Edward D. ; Feiglin, David H. ; Krol, Andrzej. / SPECT reconstruction using DCT-induced tight framelet regularization. Progress in Biomedical Optics and Imaging - Proceedings of SPIE. Vol. 9412 SPIE, 2015.
@inproceedings{e2386b9883b94f639a1e1332961ae873,
title = "SPECT reconstruction using DCT-induced tight framelet regularization",
abstract = "Wavelet transforms have been successfully applied in many fields of image processing. Yet, to our knowledge, they have never been directly incorporated to the objective function in Emission Computed Tomography (ECT) image reconstruction. Our aim has been to investigate if the 1-norm of non-decimated discrete cosine transform (DCT) coefficients of the estimated radiotracer distribution could be effectively used as the regularization term for the penalized-likelihood (PL) reconstruction, where a regularizer is used to enforce the image smoothness in the reconstruction. In this study, the 1-norm of 2D DCT wavelet decomposition was used as a regularization term. The Preconditioned Alternating Projection Algorithm (PAPA), which we proposed in earlier work to solve penalized likelihood (PL) reconstruction with non-differentiable regularizers, was used to solve this optimization problem. The DCT wavelet decompositions were performed on the transaxial reconstructed images. We reconstructed Monte Carlo simulated SPECT data obtained for a numerical phantom with Gaussian blobs as hot lesions and with a warm random lumpy background. Reconstructed images using the proposed method exhibited better noise suppression and improved lesion conspicuity, compared with images reconstructed using expectation maximization (EM) algorithm with Gaussian post filter (GPF). Also, the mean square error (MSE) was smaller, compared with EM-GPF. A critical and challenging aspect of this method was selection of optimal parameters. In summary, our numerical experiments demonstrated that the 1-norm of discrete cosine transform (DCT) wavelet frame transform DCT regularizer shows promise for SPECT image reconstruction using PAPA method.",
keywords = "DCT, Framelet, Image reconstruction, Noise reduction, Re gularization",
author = "Jiahan Zhang and Si Li and Yuesheng Xu and Schmidtlein, {C. R.} and Lipson, {Edward D.} and Feiglin, {David H.} and Andrzej Krol",
year = "2015",
doi = "10.1117/12.2082118",
language = "English (US)",
isbn = "9781628415025",
volume = "9412",
booktitle = "Progress in Biomedical Optics and Imaging - Proceedings of SPIE",
publisher = "SPIE",

}

TY - GEN

T1 - SPECT reconstruction using DCT-induced tight framelet regularization

AU - Zhang, Jiahan

AU - Li, Si

AU - Xu, Yuesheng

AU - Schmidtlein, C. R.

AU - Lipson, Edward D.

AU - Feiglin, David H.

AU - Krol, Andrzej

PY - 2015

Y1 - 2015

N2 - Wavelet transforms have been successfully applied in many fields of image processing. Yet, to our knowledge, they have never been directly incorporated to the objective function in Emission Computed Tomography (ECT) image reconstruction. Our aim has been to investigate if the 1-norm of non-decimated discrete cosine transform (DCT) coefficients of the estimated radiotracer distribution could be effectively used as the regularization term for the penalized-likelihood (PL) reconstruction, where a regularizer is used to enforce the image smoothness in the reconstruction. In this study, the 1-norm of 2D DCT wavelet decomposition was used as a regularization term. The Preconditioned Alternating Projection Algorithm (PAPA), which we proposed in earlier work to solve penalized likelihood (PL) reconstruction with non-differentiable regularizers, was used to solve this optimization problem. The DCT wavelet decompositions were performed on the transaxial reconstructed images. We reconstructed Monte Carlo simulated SPECT data obtained for a numerical phantom with Gaussian blobs as hot lesions and with a warm random lumpy background. Reconstructed images using the proposed method exhibited better noise suppression and improved lesion conspicuity, compared with images reconstructed using expectation maximization (EM) algorithm with Gaussian post filter (GPF). Also, the mean square error (MSE) was smaller, compared with EM-GPF. A critical and challenging aspect of this method was selection of optimal parameters. In summary, our numerical experiments demonstrated that the 1-norm of discrete cosine transform (DCT) wavelet frame transform DCT regularizer shows promise for SPECT image reconstruction using PAPA method.

AB - Wavelet transforms have been successfully applied in many fields of image processing. Yet, to our knowledge, they have never been directly incorporated to the objective function in Emission Computed Tomography (ECT) image reconstruction. Our aim has been to investigate if the 1-norm of non-decimated discrete cosine transform (DCT) coefficients of the estimated radiotracer distribution could be effectively used as the regularization term for the penalized-likelihood (PL) reconstruction, where a regularizer is used to enforce the image smoothness in the reconstruction. In this study, the 1-norm of 2D DCT wavelet decomposition was used as a regularization term. The Preconditioned Alternating Projection Algorithm (PAPA), which we proposed in earlier work to solve penalized likelihood (PL) reconstruction with non-differentiable regularizers, was used to solve this optimization problem. The DCT wavelet decompositions were performed on the transaxial reconstructed images. We reconstructed Monte Carlo simulated SPECT data obtained for a numerical phantom with Gaussian blobs as hot lesions and with a warm random lumpy background. Reconstructed images using the proposed method exhibited better noise suppression and improved lesion conspicuity, compared with images reconstructed using expectation maximization (EM) algorithm with Gaussian post filter (GPF). Also, the mean square error (MSE) was smaller, compared with EM-GPF. A critical and challenging aspect of this method was selection of optimal parameters. In summary, our numerical experiments demonstrated that the 1-norm of discrete cosine transform (DCT) wavelet frame transform DCT regularizer shows promise for SPECT image reconstruction using PAPA method.

KW - DCT

KW - Framelet

KW - Image reconstruction

KW - Noise reduction

KW - Re gularization

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

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

U2 - 10.1117/12.2082118

DO - 10.1117/12.2082118

M3 - Conference contribution

AN - SCOPUS:84943328557

SN - 9781628415025

VL - 9412

BT - Progress in Biomedical Optics and Imaging - Proceedings of SPIE

PB - SPIE

ER -