螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法
|
摘要
螺旋锥束计算机断层成像(CT)作为常用的临床诊断工具,如何尽可能地减少其辐射剂量是热点研究领域之一.局部成像利用准直器减小射线直照区域,能够有效降低CT辐射剂量.然而,局部成像会造成投影数据横向截断,产生局部重建问题.现有螺旋反投影滤波(BPF)算法只能实现局部曲面重建,难以实现局部体区域重建.在圆轨迹扇束BPF算法的基础上,通过加权修正和坐标扩展,提出了螺旋锥束CT倾斜扇束反投影滤波(TFB-BPF)重建算法.该算法把重建区域按层划分,对各层构建倾斜的扇形束几何,并使用经过加权修正的TFB-BPF算法逐层进行重建.该算法最大的优势是滤波线(与原始螺旋锥束BPF算法中PI线等价)在二维平面内选择,算法更加简洁高效,并且能够应用于局部体区域的重建.实验结果表明,算法能够有效重建物体局部体区域,并且重建图像质量较好,没有明显的截断伪影.
Dose reduction becomes one of the hot research fields in the most commonly used helical computed tomography(CT) for clinical diagnostic. Local imaging using a collimator can effectively lower the CT radiation dose by reducing the direct irradiation area. Due to the limitation of the exposing area, the projection data used for local imaging reconstruction are usually truncated, resulting in local reconstruction problems. The key in local image reconstruction is how to deal with the horizontal truncation of the projection data. The helical cone beam back-projection filtration(BPF) algorithm only needs to ensure the integrity of the projections of PI line to realize the reconstruction of the entire PI line, making local reconstruction possible. Due to the complexity and irregularity of the spatial distribution of PI lines, the existing helical BPF algorithms can only realize local surface reconstruction, whereas the local volume reconstruction is difficult. For the BPF algorithm in designing the PI line and the sampling points in helical cone beam CT and the difficulty in local volume reconstruction,the tilted fan-beam back-projection filtration(TFB-BPF) reconstruction algorithm is proposed by utilizing the weighted correction and coordinate expansion, based on the circular fan beam BPF. The algorithm divides the reconstruction area into several slices, constructs the inclined fan beam geometry for each layer, and slice-byslice reconstruction is conducted by using the weighted modified tilted fan beam BPF algorithm. The most powerful feature of the algorithm is that the filter line, equivalent to the PI line in the original helical BPF algorithm, is selected in a two-dimensional plane. Therefore, it is more concise and efficient and can be applied to the reconstruction of local volume regions. In this paper, the helical cone beam CT imaging geometry and the original helical cone beam BPF algorithm are introduced. Then, the TFB-BPF reconstruction algorithm is deduced. Experimental results show that the algorithm can effectively realize the local volume reconstruction and overall improved image quality without obvious truncation artifacts.
Dose reduction becomes one of the hot research fields in the most commonly used helical computed tomography(CT) for clinical diagnostic. Local imaging using a collimator can effectively lower the CT radiation dose by reducing the direct irradiation area. Due to the limitation of the exposing area, the projection data used for local imaging reconstruction are usually truncated, resulting in local reconstruction problems. The key in local image reconstruction is how to deal with the horizontal truncation of the projection data. The helical cone beam back-projection filtration(BPF) algorithm only needs to ensure the integrity of the projections of PI line to realize the reconstruction of the entire PI line, making local reconstruction possible. Due to the complexity and irregularity of the spatial distribution of PI lines, the existing helical BPF algorithms can only realize local surface reconstruction, whereas the local volume reconstruction is difficult. For the BPF algorithm in designing the PI line and the sampling points in helical cone beam CT and the difficulty in local volume reconstruction,the tilted fan-beam back-projection filtration(TFB-BPF) reconstruction algorithm is proposed by utilizing the weighted correction and coordinate expansion, based on the circular fan beam BPF. The algorithm divides the reconstruction area into several slices, constructs the inclined fan beam geometry for each layer, and slice-byslice reconstruction is conducted by using the weighted modified tilted fan beam BPF algorithm. The most powerful feature of the algorithm is that the filter line, equivalent to the PI line in the original helical BPF algorithm, is selected in a two-dimensional plane. Therefore, it is more concise and efficient and can be applied to the reconstruction of local volume regions. In this paper, the helical cone beam CT imaging geometry and the original helical cone beam BPF algorithm are introduced. Then, the TFB-BPF reconstruction algorithm is deduced. Experimental results show that the algorithm can effectively realize the local volume reconstruction and overall improved image quality without obvious truncation artifacts.
引文
[1] Huda W 2015 Curr. Raio. Rpt. 3 80
[2] Wang J, Liang Z G, Lu H B, Xing L 2010 Curr. Med.Imaging Rev. 6 72
[3] Ohno Y, Yaguc,i A, Okazaki T, Aoyagi K, Yamagata H,Sugihara N, Koyama H, Yoshikawa T, Sugimira K 2016 Eur.J. Radio. 85 1375
[4] Zheng X, Ravishankar S, Long Y, Fessler J A 2018 IEEE Trans. Med. Imaging 37 1498
[5] Shi Y Y, Yu H Y, Zhang Y B, LiI R, Kalra M, Wang G,Mou X Q 2017 the 14th International Meeting on Fully ThreeDimensional Image Reconstruction. in Radiology and Nuclear Medicine Xi'An,China, June 18-23, 2017 p837
[6] Maruyama S, Fukushima Y, Miyamae Y, Koizumi K 2018Radio. Phys. Tech. 11 235
[7] Kataria B, Althen J N, SmedbyO, Persson A, Sokjer H,Sandborg M 2018 Eur. Radiol. 28 2464
[8] Yang Q, Yan P, Zhang Y, Yu H, Shi Y, Mou X Q, Kalra MK, Zhang Y, Sun L, Wang G 2018 IEEE Trans. Med.Imaging 37 1348
[9] Yang X, Andrade V D, Scullin W, Dyer E L, Kasthuri N,Carlo F D 2018 Sci. Rep. 8 2575
[10] Yan B, Deng L, Han Y, Zhang F, Wang X Ch, Li L 2014Chin. Phys. C38 108201
[11] Han Y, Li L, Yan B, Xi X Q, Hu G E 2015 Acta. Phys. Sin.64 058704(in Chinese)[韩玉,李磊,闫镔,席晓琦,胡国恩2015物理学报64 058704]
[12] Han Y, Yan B, Li L, Xi X Q, Hu G E 2014 IEEE Trans.Nucl. Sci. 61 2753
[13] Hashemi S M, Beheshti S, Patrick R G, Narinder S P,Richard S C 2015 Comp. Math. Meth. Med. 2015 161797
[14] Sidky E Y, Kao C, Pan X C 2006 J. X-ray Sci. Tech. 14 119
[15] Sidky E Y, Pan X C 2008 Phys. Med. Biol. 53 4777
[16] Sidky E Y, Chartrand R, Boone J M, Pan X C 2014 IEEE J.Trans. Eng. Heal. Med. 2 1
[17] Courdurier M, Noo F, Defrise M, Kudo H 2008 Inverse.Probl. 24 065001
[18] Kudo H, Courdurier M, Noo F, Defrise M 2008 Phys. Med.Biol. 53 2207
[19] Ye Y, Yu H, Wei Y, Wang G 2007 J. Biol. Imaging 1 2
[20] Yu H, Ye Y, Wang G 2008 J. X-Ray Sci. Technol. 16 243
[21] Sidky E Y, Kraemer D N, Roth E G, Ullberg C, Reiser I S,Pan X 2014 J. Med. Imaging 1 031007
[22] Shi Y, Mou X 2016 Proc. SPIE 9967 99671N
[23] Zhao Y, Brun E, Coan P, Huang Z, Sztrokay A, Diemoz P C,Liebhardt S, Mittone A, Gasilov S, Miao J, Bravin A 2012PNAS 109 18290
[24] Zanette I, Bech M, Rack A, Due G L, Tafforeau P, David C,Mohr J, Pfeiffer F. Weitkamp T 2012 Natl. Acad. Sci. 10910199
[25] Turbell H, Danielsson P E 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference Toronto,Canada November 8-14, 1998 p8
[26] Noo F, Defrise M, Clackdoyle R 1999 Phys. Med. Biol. 44 561
[27] Tang X Y, Hsieh J, Roy A N, Dutta S, Samsonov D,Hagiwara A 2006 Phys. Med. Biol. 5 855
[28] Han Y, Yan B, Yu Ch Q, Li L, Li J X, Bao S L 2012 Chi7n.Phys. B 21 068701
[29] Katsevich A 2002 SIAM J. Appl. Math. 62 2012
[30] Zou Y, Pan X 2004 Phys. Med. Biol. 49 941
[31] Zou Y, Pan X 2004 Phys. Med. Biol. 49 383
[32] Zou Y, Pan X 2004 Phys. Med. Biol. 49 2717
[33] Noo F, Clackdoyle R, Pack J D 2004 Phys. Med. Biol. 493903
[34] Yu L F, Zou Y, Sidky E Y 2006 IEEE Trans. Med. Imag. 25869
[35] Wang G, Lin T H, Cheng P C, Shinozaki D M, Kim H G1991 Proc. SPIE 1556 99
[36] Shepp L A, Logan B F 1974 IEEE Trans. Nucl. Sci. 21 21
[2] Wang J, Liang Z G, Lu H B, Xing L 2010 Curr. Med.Imaging Rev. 6 72
[3] Ohno Y, Yaguc,i A, Okazaki T, Aoyagi K, Yamagata H,Sugihara N, Koyama H, Yoshikawa T, Sugimira K 2016 Eur.J. Radio. 85 1375
[4] Zheng X, Ravishankar S, Long Y, Fessler J A 2018 IEEE Trans. Med. Imaging 37 1498
[5] Shi Y Y, Yu H Y, Zhang Y B, LiI R, Kalra M, Wang G,Mou X Q 2017 the 14th International Meeting on Fully ThreeDimensional Image Reconstruction. in Radiology and Nuclear Medicine Xi'An,China, June 18-23, 2017 p837
[6] Maruyama S, Fukushima Y, Miyamae Y, Koizumi K 2018Radio. Phys. Tech. 11 235
[7] Kataria B, Althen J N, SmedbyO, Persson A, Sokjer H,Sandborg M 2018 Eur. Radiol. 28 2464
[8] Yang Q, Yan P, Zhang Y, Yu H, Shi Y, Mou X Q, Kalra MK, Zhang Y, Sun L, Wang G 2018 IEEE Trans. Med.Imaging 37 1348
[9] Yang X, Andrade V D, Scullin W, Dyer E L, Kasthuri N,Carlo F D 2018 Sci. Rep. 8 2575
[10] Yan B, Deng L, Han Y, Zhang F, Wang X Ch, Li L 2014Chin. Phys. C38 108201
[11] Han Y, Li L, Yan B, Xi X Q, Hu G E 2015 Acta. Phys. Sin.64 058704(in Chinese)[韩玉,李磊,闫镔,席晓琦,胡国恩2015物理学报64 058704]
[12] Han Y, Yan B, Li L, Xi X Q, Hu G E 2014 IEEE Trans.Nucl. Sci. 61 2753
[13] Hashemi S M, Beheshti S, Patrick R G, Narinder S P,Richard S C 2015 Comp. Math. Meth. Med. 2015 161797
[14] Sidky E Y, Kao C, Pan X C 2006 J. X-ray Sci. Tech. 14 119
[15] Sidky E Y, Pan X C 2008 Phys. Med. Biol. 53 4777
[16] Sidky E Y, Chartrand R, Boone J M, Pan X C 2014 IEEE J.Trans. Eng. Heal. Med. 2 1
[17] Courdurier M, Noo F, Defrise M, Kudo H 2008 Inverse.Probl. 24 065001
[18] Kudo H, Courdurier M, Noo F, Defrise M 2008 Phys. Med.Biol. 53 2207
[19] Ye Y, Yu H, Wei Y, Wang G 2007 J. Biol. Imaging 1 2
[20] Yu H, Ye Y, Wang G 2008 J. X-Ray Sci. Technol. 16 243
[21] Sidky E Y, Kraemer D N, Roth E G, Ullberg C, Reiser I S,Pan X 2014 J. Med. Imaging 1 031007
[22] Shi Y, Mou X 2016 Proc. SPIE 9967 99671N
[23] Zhao Y, Brun E, Coan P, Huang Z, Sztrokay A, Diemoz P C,Liebhardt S, Mittone A, Gasilov S, Miao J, Bravin A 2012PNAS 109 18290
[24] Zanette I, Bech M, Rack A, Due G L, Tafforeau P, David C,Mohr J, Pfeiffer F. Weitkamp T 2012 Natl. Acad. Sci. 10910199
[25] Turbell H, Danielsson P E 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference Toronto,Canada November 8-14, 1998 p8
[26] Noo F, Defrise M, Clackdoyle R 1999 Phys. Med. Biol. 44 561
[27] Tang X Y, Hsieh J, Roy A N, Dutta S, Samsonov D,Hagiwara A 2006 Phys. Med. Biol. 5 855
[28] Han Y, Yan B, Yu Ch Q, Li L, Li J X, Bao S L 2012 Chi7n.Phys. B 21 068701
[29] Katsevich A 2002 SIAM J. Appl. Math. 62 2012
[30] Zou Y, Pan X 2004 Phys. Med. Biol. 49 941
[31] Zou Y, Pan X 2004 Phys. Med. Biol. 49 383
[32] Zou Y, Pan X 2004 Phys. Med. Biol. 49 2717
[33] Noo F, Clackdoyle R, Pack J D 2004 Phys. Med. Biol. 493903
[34] Yu L F, Zou Y, Sidky E Y 2006 IEEE Trans. Med. Imag. 25869
[35] Wang G, Lin T H, Cheng P C, Shinozaki D M, Kim H G1991 Proc. SPIE 1556 99
[36] Shepp L A, Logan B F 1974 IEEE Trans. Nucl. Sci. 21 21