title: Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation.
abstract: OBJECTIVE
  Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction.
  
  
  METHODS
  Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according to Fessler ["Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography," IEEE Trans. Image Process. 5(3), 493-506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map.
  
  
  RESULTS
  Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization "map" designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter.
  
  
  CONCLUSIONS
  Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.
authors.affiliation: Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, Ontario M5G 2M9, Canada and Department of Biomedical Engineering, Johns Hopkins University, Baltimore Maryland 21205.
authors.forename: Grace J
authors.initials: GJ
authors.lastname: Gang
authors.affiliation: Department of Biomedical Engineering, Johns Hopkins University, Baltimore Maryland 21205.
authors.forename: J Webster
authors.initials: JW
authors.lastname: Stayman
authors.affiliation: Department of Biomedical Engineering, Johns Hopkins University, Baltimore Maryland 21205.
authors.forename: Wojciech
authors.initials: W
authors.lastname: Zbijewski
authors.affiliation: Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, Ontario M5G 2M9, Canada and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland 21205.
authors.forename: Jeffrey H
authors.initials: JH
authors.lastname: Siewerdsen
mesh_terms: D016880:Anisotropy; D016477:Artifacts; D003198:Computer Simulation; D006801:Humans; D016013:Likelihood Functions; D008954:Models, Biological; D009038:Motion; D019047:Phantoms, Imaging; D014057:Tomography, X-Ray Computed / Q000295:instrumentation / Q000379:methods*
references.citation: J Opt Soc Am A Opt Image Sci Vis. 1997 Sep;14(9):2420-42
references.pmid: 9291611
references.citation: Med Phys. 2005 May;32(5):1397-413
references.pmid: 15984691
references.citation: Med Phys. 1979 Sep-Oct;6(5):441-51
references.pmid: 492079
references.citation: Med Phys. 1998 May;25(5):614-28
references.pmid: 9608470
references.citation: Comput Biol Med. 1976 Oct;6(4):307-323
references.pmid: 1000957
references.citation: J Med Imaging (Bellingham). 2014 Oct;1(3):031013
references.pmid: 26158053
references.citation: Med Phys. 2012 Jun;39(6):3240-52
references.pmid: 22755707
references.citation: J Opt Soc Am A. 1985 Oct;2(10):1752-9
references.pmid: 4056949
references.citation: IEEE Trans Nucl Sci. 2013 Jan 11;60(1):182-193
references.pmid: 24436497
references.citation: Med Phys. 2011 Oct;38(10):5612-29
references.pmid: 21992379
references.citation: Phys Med Biol. 2013 Oct 21;58(20):7419-45
references.pmid: 24080994
references.citation: IEEE Trans Med Imaging. 2004 Dec;23(12):1543-56
references.pmid: 15575411
references.citation: Phys Med Biol. 2013 Jun 21;58(12):R63-96
references.pmid: 23739261
references.citation: Med Phys. 2010 May;37(5):2074-81
references.pmid: 20527540
references.citation: Phys Med Biol. 1987 May;32(5):565-75
references.pmid: 3588670
references.citation: Proc SPIE Int Soc Opt Eng. 2013 Feb;8668:
references.pmid: 34248250
references.citation: Med Phys. 2008 Dec;35(12):5510-29
references.pmid: 19175110
references.citation: IEEE Trans Image Process. 1996;5(3):493-506
references.pmid: 18285134
references.citation: IEEE Trans Med Imaging. 2000 Jun;19(6):601-15
references.pmid: 11026463
references.citation: Med Phys. 2009 Feb;36(2):500-12
references.pmid: 19291989
references.citation: IEEE Trans Med Imaging. 2006 Jan;25(1):28-41
references.pmid: 16398412
references.citation: Med Phys. 2003 Apr;30(4):590-600
references.pmid: 12722811
references.citation: IEEE Trans Med Imaging. 2001 Aug;20(8):815-22
references.pmid: 11513032
references.citation: Med Phys. 2011 Mar;38(3):1444-58
references.pmid: 21520856
references.citation: Phys Med Biol. 1999 Nov;44(11):2835-51
references.pmid: 10588288
references.citation: Med Phys. 2012 Aug;39(8):5145-56
references.pmid: 22894440
references.citation: J Opt Soc Am A Opt Image Sci Vis. 2006 Dec;23(12):2989-96
references.pmid: 17106455
references.citation: J Nucl Med. 1992 Mar;33(3):451-7
references.pmid: 1740719
references.citation: Med Phys. 2011 Nov;38(11):6188
references.pmid: 22047384
references.citation: J Opt Soc Am A Opt Image Sci Vis. 2007 Dec;24(12):B99-B109
references.pmid: 18059919
references.citation: Phys Med Biol. 1984 Nov;29(11):1343-52
references.pmid: 6505016
references.citation: Med Phys. 2011 Jul;38 Suppl 1:S2
references.pmid: 21978114
references.citation: Proc SPIE Int Soc Opt Eng. 2014 Mar 19;9033:903335
references.pmid: 25328638
references.citation: Med Phys. 2013 Mar;40(3):031104
references.pmid: 23464284
references.citation: Int J Radiat Oncol Biol Phys. 1993 Oct 20;27(3):697-705
references.pmid: 8226167
references.citation: Med Phys. 1997 Jan;24(1):71-89
references.pmid: 9029542
references.citation: Med Phys. 2011 Apr;38(4):2122-31
references.pmid: 21626945
references.citation: Phys Med Biol. 2011 Oct 7;56(19):6223-42
references.pmid: 21896963
references.citation: IEEE Trans Med Imaging. 2002 Sep;21(9):1040-50
references.pmid: 12564872
references.citation: Med Phys. 2008 Nov;35(11):5043-53
references.pmid: 19070238
references.citation: Med Phys. 2013 Mar;40(3):031908
references.pmid: 23464323
references.citation: Theranostics. 2013 Oct 04;3(10):741-56
references.pmid: 24312148
references.citation: J Nucl Med. 1993 Jul;34(7):1198-203
references.pmid: 8315502
references.citation: Med Phys. 2000 Aug;27(8):1903-14
references.pmid: 10984236
references.citation: Phys Med Biol. 2006 Aug 21;51(16):4017-29
references.pmid: 16885621
references.citation: Med Phys. 2012 Jun;39(6):3214-28
references.pmid: 22755705
references.citation: Med Phys. 2010 Apr;37(4):1591-600
references.pmid: 20443480
references.citation: IEEE Trans Image Process. 1996;5(9):1346-58
references.pmid: 18285223
references.citation: Med Phys. 2002 Nov;29(11):2655-71
references.pmid: 12462733
references.citation: IEEE Trans Med Imaging. 2010 May;29(5):1097-113
references.pmid: 20236879
references.citation: IEEE Trans Med Imaging. 2004 Mar;23(3):269-84
references.pmid: 15027520
references.citation: Proc Natl Acad Sci U S A. 1993 Nov 1;90(21):9758-65
references.pmid: 8234311
references.citation: Med Phys. 2012 Jul;39(7):4115-22
references.pmid: 22830744
references.citation: IEEE Trans Med Imaging. 2010 Apr;29(4):1050-8
references.pmid: 20335088
references.citation: Phys Med Biol. 2008 May 21;53(10):2471-93
references.pmid: 18424878
references.citation: Med Phys. 2011 Apr;38(4):1754-68
references.pmid: 21626910