Optimization Online


An Improved Method of Total Variation Superiorization Applied to Reconstruction in Proton Computed Tomography

Blake Schultze(Blake_Schultze***at***baylor.edu)
Yair Censor(yair***at***math.haifa.ac.il)
Paniz Karbasi(Paniz_Karbasi***at***baylor.edu)
Keith E. Schubert(Keith_Schubert***at***baylor.edu)
Reinhard W. Schulte(rschulte***at***llu.edu)

Abstract: Previous work showed that total variation superiorization (TVS) improves reconstructed image quality in proton computed tomography (pCT). The structure of the TVS algorithm has evolved since then and this work investigated if this new algorithmic structure provides additional benefits to pCT image quality. Structural and parametric changes introduced to the original TVS algorithm included: (1) inclusion or exclusion of TV reduction requirement, (2) a variable number, $N$, of TV perturbation steps per feasibility-seeking iteration, and (3) introduction of a perturbation kernel $0<\alpha<1$. The structural change of excluding the TV reduction requirement check tended to have a beneficial effect for $3\le N\le 6$ and allows full parallelization of the TVS algorithm. Repeated perturbations per feasibility-seeking iterations reduced total variation (TV) and material dependent standard deviations for $3\le N\le 6$. The perturbation kernel $\alpha$, equivalent to $\alpha=0.5$ in the original TVS algorithm, reduced TV and standard deviations as $\alpha$ was increased beyond $\alpha=0.5$, but negatively impacted reconstructed relative stopping power (RSP) values for $\alpha>0.75$. The reductions in TV and standard deviations allowed feasibility-seeking with a larger relaxation parameter $\lambda$ than previously used, without the corresponding increases in standard deviations experienced with the original TVS algorithm. This work demonstrates that the modifications related to the evolution of the original TVS algorithm provide benefits in terms of both pCT image quality and computational efficiency for appropriately chosen parameter values.

Keywords: Feasibility-seeking algorithms, image reconstruction, perturbations, proton computed tomography (pCT), superiorization, total variation superiorization (TVS)

Category 1: Applications -- Science and Engineering (Biomedical Applications )

Category 2: Nonlinear Optimization (Constrained Nonlinear Optimization )

Category 3: Convex and Nonsmooth Optimization (Convex Optimization )

Citation: Technical Report, March 3, 2018.

Download: [PDF]

Entry Submitted: 03/06/2018
Entry Accepted: 03/06/2018
Entry Last Modified: 03/06/2018

Modify/Update this entry

  Visitors Authors More about us Links
  Subscribe, Unsubscribe
Digest Archive
Search, Browse the Repository


Coordinator's Board
Classification Scheme
Give us feedback
Optimization Journals, Sites, Societies
Mathematical Optimization Society