Low-dose 4DCT reconstruction via temporal nonlocal means

Zhen Tian, Xun Jia, Bin Dong, Yifei Lou, Steve B. Jiang

Research output: Contribution to journalArticle

47 Citations (Scopus)

Abstract

Purpose: Four-dimensional computed tomography (4DCT) has been widely used in cancer radiotherapy for accurate target delineation and motion measurement for tumors in the thorax and upper abdomen areas. However, its prolonged scanning duration causes a considerable increase of radiation dose compared to conventional CT, which is a major concern in its clinical application. This work is to develop a new algorithm to reconstruct 4DCT images from undersampled projections acquired at low mA s levels in order to reduce the imaging dose. Methods: Conventionally, each phase of 4DCT is reconstructed independently using the filtered backprojection (FBP) algorithm. The basic idea of the authors' new algorithm is that by utilizing the common information among different phases, the input information required to reconstruct the image of high quality, and thus the imaging dose, can be reduced. The authors proposed a temporal nonlocal means (TNLM) method to explore the interphase similarity. All phases of the 4DCT images are reconstructed simultaneously by minimizing a cost function consisting of a data fidelity term and a TNLM regularization term. The authors utilized a modified forward-backward splitting algorithm and a Gauss-Jacobi iteration method to efficiently solve the minimization problem. The algorithm was also implemented on a graphics processing unit (GPU) to improve the computational speed. The authors' reconstruction algorithm has been tested on a digital NCAT thorax phantom in three low dose scenarios: All projections with low mA s level, undersampled projections with high mA s level, and undersampled projections with low mA s level. Results: In all three low dose scenarios, the new algorithm generates visually much better CT images containing less image noise and streaking artifacts compared to the standard FBP algorithm. Quantitative analysis shows that by comparing the authors' TNLM algorithm to the standard FBP algorithm, the contrast-to-noise ratio has been improved by a factor of 3.9-10.2 and the signal-to-noise ratio has been improved by a factor of 2.1-5.9, depending on the cases. In the situation of undersampled projection data, the majority of the streaks in the images reconstructed by FBP can be suppressed using the authors' algorithm. The total reconstruction time for all ten phases of a slice ranges from 40 to 90 s on an NVIDIA Tesla C1060 GPU card. Conclusions: The experimental results indicate that the authors' new algorithm outperforms the conventional FBP algorithm in effectively reducing the image artifacts due to undersampling and suppressing the image noise due to the low mA s level.

Original languageEnglish (US)
Pages (from-to)1359-1365
Number of pages7
JournalMedical Physics
Volume38
Issue number3
DOIs
StatePublished - Mar 2011

Fingerprint

Artifacts
Thorax
Four-Dimensional Computed Tomography
Interphase
Signal-To-Noise Ratio
Abdomen
Noise
Neoplasms
Radiotherapy
Radiation
Costs and Cost Analysis

Keywords

  • 4DCT reconstruction
  • dose reduction
  • GPU
  • temporal nonlocal means

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

Low-dose 4DCT reconstruction via temporal nonlocal means. / Tian, Zhen; Jia, Xun; Dong, Bin; Lou, Yifei; Jiang, Steve B.

In: Medical Physics, Vol. 38, No. 3, 03.2011, p. 1359-1365.

Research output: Contribution to journalArticle

Tian, Zhen ; Jia, Xun ; Dong, Bin ; Lou, Yifei ; Jiang, Steve B. / Low-dose 4DCT reconstruction via temporal nonlocal means. In: Medical Physics. 2011 ; Vol. 38, No. 3. pp. 1359-1365.
@article{45416c649dcc4e958ea89795f98ac6ce,
title = "Low-dose 4DCT reconstruction via temporal nonlocal means",
abstract = "Purpose: Four-dimensional computed tomography (4DCT) has been widely used in cancer radiotherapy for accurate target delineation and motion measurement for tumors in the thorax and upper abdomen areas. However, its prolonged scanning duration causes a considerable increase of radiation dose compared to conventional CT, which is a major concern in its clinical application. This work is to develop a new algorithm to reconstruct 4DCT images from undersampled projections acquired at low mA s levels in order to reduce the imaging dose. Methods: Conventionally, each phase of 4DCT is reconstructed independently using the filtered backprojection (FBP) algorithm. The basic idea of the authors' new algorithm is that by utilizing the common information among different phases, the input information required to reconstruct the image of high quality, and thus the imaging dose, can be reduced. The authors proposed a temporal nonlocal means (TNLM) method to explore the interphase similarity. All phases of the 4DCT images are reconstructed simultaneously by minimizing a cost function consisting of a data fidelity term and a TNLM regularization term. The authors utilized a modified forward-backward splitting algorithm and a Gauss-Jacobi iteration method to efficiently solve the minimization problem. The algorithm was also implemented on a graphics processing unit (GPU) to improve the computational speed. The authors' reconstruction algorithm has been tested on a digital NCAT thorax phantom in three low dose scenarios: All projections with low mA s level, undersampled projections with high mA s level, and undersampled projections with low mA s level. Results: In all three low dose scenarios, the new algorithm generates visually much better CT images containing less image noise and streaking artifacts compared to the standard FBP algorithm. Quantitative analysis shows that by comparing the authors' TNLM algorithm to the standard FBP algorithm, the contrast-to-noise ratio has been improved by a factor of 3.9-10.2 and the signal-to-noise ratio has been improved by a factor of 2.1-5.9, depending on the cases. In the situation of undersampled projection data, the majority of the streaks in the images reconstructed by FBP can be suppressed using the authors' algorithm. The total reconstruction time for all ten phases of a slice ranges from 40 to 90 s on an NVIDIA Tesla C1060 GPU card. Conclusions: The experimental results indicate that the authors' new algorithm outperforms the conventional FBP algorithm in effectively reducing the image artifacts due to undersampling and suppressing the image noise due to the low mA s level.",
keywords = "4DCT reconstruction, dose reduction, GPU, temporal nonlocal means",
author = "Zhen Tian and Xun Jia and Bin Dong and Yifei Lou and Jiang, {Steve B.}",
year = "2011",
month = "3",
doi = "10.1118/1.3547724",
language = "English (US)",
volume = "38",
pages = "1359--1365",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "3",

}

TY - JOUR

T1 - Low-dose 4DCT reconstruction via temporal nonlocal means

AU - Tian, Zhen

AU - Jia, Xun

AU - Dong, Bin

AU - Lou, Yifei

AU - Jiang, Steve B.

PY - 2011/3

Y1 - 2011/3

N2 - Purpose: Four-dimensional computed tomography (4DCT) has been widely used in cancer radiotherapy for accurate target delineation and motion measurement for tumors in the thorax and upper abdomen areas. However, its prolonged scanning duration causes a considerable increase of radiation dose compared to conventional CT, which is a major concern in its clinical application. This work is to develop a new algorithm to reconstruct 4DCT images from undersampled projections acquired at low mA s levels in order to reduce the imaging dose. Methods: Conventionally, each phase of 4DCT is reconstructed independently using the filtered backprojection (FBP) algorithm. The basic idea of the authors' new algorithm is that by utilizing the common information among different phases, the input information required to reconstruct the image of high quality, and thus the imaging dose, can be reduced. The authors proposed a temporal nonlocal means (TNLM) method to explore the interphase similarity. All phases of the 4DCT images are reconstructed simultaneously by minimizing a cost function consisting of a data fidelity term and a TNLM regularization term. The authors utilized a modified forward-backward splitting algorithm and a Gauss-Jacobi iteration method to efficiently solve the minimization problem. The algorithm was also implemented on a graphics processing unit (GPU) to improve the computational speed. The authors' reconstruction algorithm has been tested on a digital NCAT thorax phantom in three low dose scenarios: All projections with low mA s level, undersampled projections with high mA s level, and undersampled projections with low mA s level. Results: In all three low dose scenarios, the new algorithm generates visually much better CT images containing less image noise and streaking artifacts compared to the standard FBP algorithm. Quantitative analysis shows that by comparing the authors' TNLM algorithm to the standard FBP algorithm, the contrast-to-noise ratio has been improved by a factor of 3.9-10.2 and the signal-to-noise ratio has been improved by a factor of 2.1-5.9, depending on the cases. In the situation of undersampled projection data, the majority of the streaks in the images reconstructed by FBP can be suppressed using the authors' algorithm. The total reconstruction time for all ten phases of a slice ranges from 40 to 90 s on an NVIDIA Tesla C1060 GPU card. Conclusions: The experimental results indicate that the authors' new algorithm outperforms the conventional FBP algorithm in effectively reducing the image artifacts due to undersampling and suppressing the image noise due to the low mA s level.

AB - Purpose: Four-dimensional computed tomography (4DCT) has been widely used in cancer radiotherapy for accurate target delineation and motion measurement for tumors in the thorax and upper abdomen areas. However, its prolonged scanning duration causes a considerable increase of radiation dose compared to conventional CT, which is a major concern in its clinical application. This work is to develop a new algorithm to reconstruct 4DCT images from undersampled projections acquired at low mA s levels in order to reduce the imaging dose. Methods: Conventionally, each phase of 4DCT is reconstructed independently using the filtered backprojection (FBP) algorithm. The basic idea of the authors' new algorithm is that by utilizing the common information among different phases, the input information required to reconstruct the image of high quality, and thus the imaging dose, can be reduced. The authors proposed a temporal nonlocal means (TNLM) method to explore the interphase similarity. All phases of the 4DCT images are reconstructed simultaneously by minimizing a cost function consisting of a data fidelity term and a TNLM regularization term. The authors utilized a modified forward-backward splitting algorithm and a Gauss-Jacobi iteration method to efficiently solve the minimization problem. The algorithm was also implemented on a graphics processing unit (GPU) to improve the computational speed. The authors' reconstruction algorithm has been tested on a digital NCAT thorax phantom in three low dose scenarios: All projections with low mA s level, undersampled projections with high mA s level, and undersampled projections with low mA s level. Results: In all three low dose scenarios, the new algorithm generates visually much better CT images containing less image noise and streaking artifacts compared to the standard FBP algorithm. Quantitative analysis shows that by comparing the authors' TNLM algorithm to the standard FBP algorithm, the contrast-to-noise ratio has been improved by a factor of 3.9-10.2 and the signal-to-noise ratio has been improved by a factor of 2.1-5.9, depending on the cases. In the situation of undersampled projection data, the majority of the streaks in the images reconstructed by FBP can be suppressed using the authors' algorithm. The total reconstruction time for all ten phases of a slice ranges from 40 to 90 s on an NVIDIA Tesla C1060 GPU card. Conclusions: The experimental results indicate that the authors' new algorithm outperforms the conventional FBP algorithm in effectively reducing the image artifacts due to undersampling and suppressing the image noise due to the low mA s level.

KW - 4DCT reconstruction

KW - dose reduction

KW - GPU

KW - temporal nonlocal means

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

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

U2 - 10.1118/1.3547724

DO - 10.1118/1.3547724

M3 - Article

C2 - 21520846

AN - SCOPUS:79952131914

VL - 38

SP - 1359

EP - 1365

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 3

ER -