Markerless Prediction of Diaphragm Displacement Based on Two-Step Subspace Mapping
-
摘要: 呼吸会引起体内器官和肿瘤的运动, 这会显著影响放射治疗的过程和效果. 人体内部膈肌和胸腹部外表面是当前两种与呼吸系统高度相关的结构, 本文对其进行系统研究, 提出了一种新的分步子空间映射(TSSM)算法, 通过对体外胸腹部表面的测量, 来预测体内膈肌的运动. 本文首先采用三维图像分割技术对4D CT图像进行分割, 在不使用标记物的情况下, 准确测量体内膈肌和体外胸腹部表面的位移. 为了解决跨空间的预测问题, TSSM首先构造特征子空间, 并将膈肌数据和胸腹外表面数据分别映射到各自的子空间中, 以减少数据的相关性和冗余信息; 然后通过线性岭回归优化过程, 对两个子空间进行二次映射, 从而有效地捕获跨空间数据之间的相关性. 根据训练得到的相关模型, 通过体外胸腹部外表面的运动情况, 对体内膈肌的运动情况进行准确的预测. 为了研究数据之间的非线性关系, 本文进一步将TSSM推广到了基于核的TSSM(kTSSM)算法. 实验表明, 该方法可以根据腹腔外表面的运动情况, 准确的对体内膈肌位移进行预测, 优于经典的线性模型和ANN模型. 本文给出了优化算法的解析解, 其运算速度快, 将有助于提高放射治疗中门控技术和跟踪技术的效率和精度.Abstract: Respiratory induced organ and tumor motion has great influence in radiation therapy. Two highly correlated structures with respiratory are comprehensively investigated including internal diaphragm and external thoracoabdominal surface in this paper. A novel two-step subspace mapping (TSSM) algorithm is proposed to predict the diaphragm displacement by markerless thoracoabdominal surface measurement. TSSM first incorporates 3D image segmentation to accurately measure the displacement of the diaphragm and the thoracoabdominal surface without markers on 4D CT images. To solve the cross-domain prediction problem, TSSM first constructs eigenspaces and then make a mapping in high-dimension subspace, which can effectively capture the correlation between cross-domain data. Experiments show the proposed method can accurately predict the displacement of the internal diaphragm by the external thoracoabdominal surface, and outperform the classical linear model and ANN model. The simple closed form for the optimization algorithm leads to an extremely fast algorithm, which has potential for improving the timing accuracy of surface-guided gating and tracking in radiotherapy.
-
Key words:
- respiratory motion /
- thoracoabdominal surface /
- diaphragm /
- subspace mapping /
- regression model /
- 4D CT
-
图 1 三维图像分割结果冠状平面展示图, 背景体素值为0, 身体区域的体素值为1, 肺部体素值为3.
Fig. 1 Segmentation images on the axial plane, where the background voxels is set to 0, the voxels of body area to 1, and the voxels of lungs to 3
图 2 三维图像分割结果, 及其对应的肺部区域和身体区域
Fig. 2 A 3D segmentation result and the corresponding separated 3D lungs and body
图 4 胸腹部表面, (a)初始胸腹部表面, (b)胸腹部掩膜, (c)最终胸腹部表面
Fig. 4 The thoracoabdominal surface. (a) is the initial thoracoabdominal surface obtained from 3D body mask, (b) is the produced thoracoabdominal mask, and (c) is the final thoracoabdominal surface.
图 6 膈肌和胸腹部表面的位移. 每一种颜色的曲线对应一个患者的数据.
Fig. 6 The displacement of diaphragm and thoracoabdominal surface, where each color corresponds to a specific patient’s data.
图 7 归一化后的膈肌和胸腹部表面的位移. 每一种颜色的曲线对应一个患者的数据.
Fig. 7 The displacement of diaphragm and thoracoabdominal surface after normalization, where each color corresponds to a specific patient’s data.
图 8 三个方向的预测结果. 该结果对应于100次运行MSE的中位数值.
Fig. 8 Prediction results corresponding the median MSE value in three directions after 100 independent runs.
图 9 三个方向各相位预测误差的统计箱图
Fig. 9 Boxplot of the prediction performance on every phase in three directions
图 10 主方向的预测结果. 该结果对应于100次运行MSE的中位数值.
Fig. 10 Prediction results corresponding the median MSE value in the principal direction after 100 independent runs.
图 11 主方向各相位预测误差的统计箱图
Fig. 11 Boxplot of the prediction performance on every phase in the principal direction
图 14 PCA中的阈值ε对TSSM预测模型的影响
Fig. 14 The influence of ε in PCA on the prediction performance of TSSM
表 1 预测结果在三个方向上的统计结果
Table 1 Statistic results of the prediction results in three directions
OLS模型[1] ANN TSSM(线性) kTSSM(多项式) kTSSM(高斯) 没有噪音 MSE 0.45 (0.23) 0.30 (0.16) 0.20 (0.09) 0.18 (0.08) 0.20 (0.10) R2 0.70 (0.14) 0.81 (0.07) 0.87 (0.05) 0.88 (0.04) 0.88 (0.04) MAPE 232.11 (156.44) 159.60 (112.37) 150.82 (83.98) 116.57 (71.28) 130.69 (87.26) 噪音($\sigma =0.1$) MSE 0.46 (0.22) 0.29 (0.15) 0.22 (0.10) 0.19 (0.09) 0.22 (0.09) R2 0.70 (0.15) 0.8207 (0.07) 0.86 (0.05) 0.88 (0.04) 0.86 (0.05) MAPE 276.86 (551.10) 214.82 (257.92) 157.28 (117.30) 137.25 (92.36) 170.47 (204.00) 噪音($\sigma =0.2$) MSE 0.54 (0.24) 0.39 (0.18) 0.25 (0.11) 0.23 (0.10) 0.24 (0.10) R2 0.69 (0.13) 0.76 (0.10) 0.84 (0.05) 0.86 (0.04) 0.85 (0.06) MAPE 277.42 (348.60) 211.14 (196.70) 172.08 (171.03) 131.98 (78.23) 166.12 (156.08) 噪音($\sigma =0.3$) MSE 0.65 (0.25) 0.42 (0.21) 0.30 (0.11) 0.28 (0.11) 0.29 (0.12) R2 0.59 (0.20) 0.75 (0.11) 0.82 (0.05) 0.82 (0.06) 0.82 (0.06) MAPE 426.45 (1371.37) 211.37 (228.00) 242.75 (274.21) 157.12 (94.87) 179.81 (115.10) 噪音($\sigma =0.4$) MSE 0.80 (0.28) 0.52 (0.33) 0.37 (0.13) 0.33 (0.13) 0.38 (0.13) R2 0.53 (0.21) 0.69 (0.19) 0.79 (0.07) 0.80 (0.06) 0.78 (0.07) MAPE 592.47 (1687.33) 262.36 (405.00) 198.89 (112.67) 174.14 (116.58) 185.34 (133.87) 噪音($\sigma =0.5$) MSE 0.98 (0.33) 0.65 (0.29) 0.46 (0.16) 0.41 (0.15) 0.45 (0.17) R2 0.47 (0.22) 0.64 (0.15) 0.75 (0.08) 0.76 (0.08) 0.76 (0.09) MAPE 365.33 (399.98) 254.35 (291.97) 234.93 (185.28) 194.18 (150.11) 214.82 (260.37) 
下载: 导出CSV
表 2 主方向预测结果的统计结果
Table 2 Statistic results of the prediction results in principal direction
OLS模型[1] ANN TSSM(线性) kTSSM(多项式) kTSSM(高斯) 没有噪音 MSE 0.47 (0.28) 0.27 (0.18) 0.14 (0.07) 0.13 (0.06) 0.15 (0.06) R2 0.73 (0.20) 0.86 (0.08) 0.93 (0.03) 0.93 (0.03) 0.92 (0.03) MAPE 91.07 (38.34) 63.62 (23.72) 47.85 (15.31) 45.17 (14.94) 46.47 (14.39) 噪音($\sigma =0.1$) MSE 0.40 (0.26) 0.25 (0.18) 0.15 (0.07) 0.14 (0.06) 0.15 (0.07) R2 0.78 (0.16) 0.86 (0.08) 0.92 (0.04) 0.92 (0.03) 0.91 (0.03) MAPE 143.90 (173.93) 100.38 (136.34) 116.46 (198.17) 68.29 (66.63) 107.08 (145.85) 噪音($\sigma =0.2$) MSE 0.43 (0.23) 0.29 (0.12) 0.18 (0.07) 0.17 (0.06) 0.18 (0.07) R2 0.76 (0.14) 0.84 (0.07) 0.90 (0.04) 0.91 (0.04) 0.90 (0.04) MAPE 188.52 (476.86) 138.61 (386.48) 133.32 (129.00) 79.96 (53.09) 82.96 (79.07) 噪音($\sigma =0.3$) MSE 0.61 (0.33) 0.38 (0.18) 0.23 (0.10) 0.20 (0.08) 0.24 (0.09) R2 0.64 (0.20) 0.79 (0.11) 0.87 (0.06) 0.87 (0.06) 0.88 (0.05) MAPE 192.18 (223.53) 136.84 (116.74) 134.21 (151.02) 86.07 (64.90) 101.09 (88.78) 噪音($\sigma =0.4$) MSE 0.70 (0.31) 0.45 (0.36) 0.30 (0.12) 0.26 (0.09) 0.29 (0.11) R2 0.65 (0.17) 0.77 (0.21) 0.83 (0.07) 0.87 (0.06) 0.85 (0.07) MAPE 221.12 (417.10) 175.31 (342.96) 142.09 (168.43) 94.42 (64.70) 109.11 (116.15) 噪音($\sigma =0.5$) MSE 1.01 (0.50) 0.52 (0.28) 0.39 (0.13) 0.32 (0.13) 0.36 (0.14) R2 0.47 (0.35) 0.75 (0.15) 0.80 (0.09) 0.84 (0.07) 0.81 (0.09) MAPE 173.44 (98.15) 193.44 (455.90) 146.14 (212.00) 102.30 (130.34) 94.49)
下载: 导出CSV
表 3 Parameter setting of TSSM和kTSSM for optimal results
模型 TSSM(线性) kTSSM(多项式) kTSSM(高斯) 参数 λ λ d c λ σkernel 三个方向的预测 没有噪音 1 100 6 9 0.1 0.2 噪音($\sigma =0.1$) 1 1000 3.5 9 0.1 0.2 噪音($\sigma =0.2$) 1 100 2.5 9 0.1 0.1 噪音($\sigma =0.3$) 1 1 0.5 4 0.1 0.1 噪音($\sigma =0.4$) 1 10 1.5 7 0.1 0.1 噪音($\sigma =0.5$) 1 10 1 4 1 0.1 主方向的预测 没有噪音 0.1 1 1 5 0.1 0.2 噪音($\sigma =0.1$) 0.1 0.1 1 7 0.1 0.2 噪音($\sigma =0.2$) 1 0.1 0.5 4 0.1 0.2 噪音($\sigma =0.3$) 1 1 0.5 3 0.1 0.1 噪音($\sigma =0.4$) 1 1 1 4 0.1 0.2 噪音($\sigma =0.5$) 1 0.1 0.5 7 0.1 0.1
下载: 导出CSV
-
[1] K. T. Malinowski, T. J. McAvoy, R. George, S. Dieterich, and W. D. D'Souza. Mitigating errors in external respiratory surrogate-based models of tumor position. International Journal of Radiation Oncology* Biology* Physics, 2012, 82(5): e709−e716 doi: 10.1016/j.ijrobp.2011.05.042 [2] H. Fayad et al., “A patient specific respiratory model based on 4D CT data and a time of flight camera (TOF),” in 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC), 2009, pp. 2594−2598: IEEE. [3] J. Schwaab et al. Ultrasound tracking for intra-fractional motion compensation in radiation therapy. Physica Medica, 2014, 30(5): 578−582 doi: 10.1016/j.ejmp.2014.03.003 [4] C. Ozhasoglu and M. J. Murphy. Issues in respiratory motion compensation during external-beam radiotherapy. International Journal of Radiation Oncology* Biology* Physics, 2002, 52(5): 1389−1399 doi: 10.1016/S0360-3016(01)02789-4 [5] S. Shimizu et al. Impact of respiratory movement on the computed tomographic images of small lung tumors in three-dimensional (3D) radiotherapy. International Journal of Radiation Oncology* Biology* Physics, 2000, 46(5): 1127−1133 doi: 10.1016/S0360-3016(99)00352-1 [6] J. R. McClelland et al. A continuous 4D motion model from multiple respiratory cycles for use in lung radiotherapy. Medical Physics, 2006, 33(9): 3348−3358 doi: 10.1118/1.2222079 [7] A. Torshabi, A. Pella, M. Riboldi, and G. Baroni. Targeting accuracy in real-time tumor tracking via external surrogates: a comparative study. Technology in Cancer Research & Treatment, 2010, 9(6): 551−561 [8] C. Paganelli et al. Magnetic resonance imaging–guided versus surrogate-based motion tracking in liver radiation therapy: A prospective comparative study. International Journal of Radiation Oncology* Biology* Physics, 2015, 91(4): 840−848 doi: 10.1016/j.ijrobp.2014.12.013 [9] H. Fayad, T. Pan, J. François Clement, and D. Visvikis. Correlation of respiratory motion between external patient surface and internal anatomical landmarks. Medical physics, 2011, 38(6Part1): 3157−3164 doi: 10.1118/1.3589131 [10] J. Hanley et al. Deep inspiration breath-hold technique for lung tumors: the potential value of target immobilization and reduced lung density in dose escalation. International Journal of Radiation Oncology* Biology* Physics, 1999, 45(3): 603−611 doi: 10.1016/S0360-3016(99)00154-6 [11] D. Mah et al. Technical aspects of the deep inspiration breath-hold technique in the treatment of thoracic cancer. International Journal of Radiation Oncology* Biology* Physics, 2000, 48(4): 1175−1185 doi: 10.1016/S0360-3016(00)00747-1 [12] R. Koshani, J. M. Balter, J. A. Hayman, G. T. Henning, and M. van Herk. Short-term and long-term reproducibility of lung tumor position using active breathing control (ABC). International Journal of Radiation Oncology* Biology* Physics, 2006, 65(5): 1553−1559 doi: 10.1016/j.ijrobp.2006.04.027 [13] W. Lu, K. J. Ruchala, M.-L. Chen, Q. Chen, and G. H. Olivera. Real-time respiration monitoring using the radiotherapy treatment beam and four-dimensional computed tomography (4DCT)—a conceptual study. Physics in Medicine & Biology, 2006, 51(18): 4469 [14] S. Minohara, T. Kanai, M. Endo, K. Noda, and M. Kanazawa. Respiratory gated irradiation system for heavy-ion radiotherapy. International Journal of Radiation Oncology* Biology* Physics, 2000, 47(4): 1097−1103 doi: 10.1016/S0360-3016(00)00524-1 [15] H. D. Kubo and B. C. Hill. Respiration gated radiotherapy treatment: a technical study. Physics in Medicine & Biology, 1996, 41(1): 83 [16] Y. Seppenwoolde et al. Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. International Journal of Radiation Oncology* Biology* Physics, 2002, 53(4): 822−834 doi: 10.1016/S0360-3016(02)02803-1 [17] C.-C. Shieh et al. A Bayesian approach for three-dimensional markerless tumor tracking using kV imaging during lung radiotherapy. Physics in Medicine & Biology, 2017, 62(8): 3065 [18] F. Preiswerk et al. Model-guided respiratory organ motion prediction of the liver from 2D ultrasound. Medical Image Analysis, 2014, 18(5): 740−751 doi: 10.1016/j.media.2014.03.006 [19] S. Shimizu et al. Detection of lung tumor movement in real-time tumor-tracking radiotherapy. International Journal of Radiation Oncology* Biology* Physics, 2001, 51(2): 304−310 doi: 10.1016/S0360-3016(01)01641-8 [20] S. Vedam, V. Kini, P. Keall, V. Ramakrishnan, H. Mostafavi, and R. Mohan. Quantifying the predictability of diaphragm motion during respiration with a noninvasive external marker. Medical Physics, 2003, 30(4): 505−513 doi: 10.1118/1.1558675 [21] R. I. Berbeco, S. B. Jiang, G. C. Sharp, G. T. Chen, H. Mostafavi, and H. Shirato. Integrated radiotherapy imaging system (IRIS): design considerations of tumour tracking with linac gantry-mounted diagnostic x-ray systems with flat-panel detectors. Physics in Medicine & Biology, 2004, 49(2): 243 [22] X. Tang, G. C. Sharp, and S. B. Jiang. Fluoroscopic tracking of multiple implanted fiducial markers using multiple object tracking. Physics in Medicine & Biology, 2007, 52(14): 4081 [23] G. C. Sharp, S. B. Jiang, S. Shimizu, and H. Shirato. Tracking errors in a prototype real-time tumour tracking system. Physics in Medicine & Biology, 2004, 49(23): 5347 [24] J. T. Booth et al. The first patient treatment of electromagnetic-guided real time adaptive radiotherapy using MLC tracking for lung SABR. Radiotherapy and Oncology, 2016, 121(1): 19−25 doi: 10.1016/j.radonc.2016.08.025 [25] P. Kupelian et al. Multi-institutional clinical experience with the Calypso System in localization and continuous, real-time monitoring of the prostate gland during external radiotherapy. International Journal of Radiation Oncology* Biology* Physics, 2007, 67(4): 1088−1098 doi: 10.1016/j.ijrobp.2006.10.026 [26] T. R. Willoughby et al. Target localization and real-time tracking using the Calypso 4D localization system in patients with localized prostate cancer. International Journal of Radiation Oncology* Biology* Physics, 2006, 65(2): 528−534 doi: 10.1016/j.ijrobp.2006.01.050 [27] J. D. Hoisak, K. E. Sixel, R. Tirona, P. C. Cheung, and J.-P. Pignol. Correlation of lung tumor motion with external surrogate indicators of respiration. International Journal of Radiation Oncology* Biology* Physics, 2004, 60(4): 1298−1306 doi: 10.1016/j.ijrobp.2004.07.681 [28] Q. Zhang et al. A patient‐specific respiratory model of anatomical motion for radiation treatment planning. Medical Physics, 2007, 34(12): 4772−4781 doi: 10.1118/1.2804576 [29] A. Fassi, J. Schaerer, M. Fernandes, M. Riboldi, D. Sarrut, and G. Baroni. Tumor tracking method based on a deformable 4D CT breathing motion model driven by an external surface surrogate. International Journal of Radiation Oncology* Biology* Physics, 2014, 88(1): 182−188 doi: 10.1016/j.ijrobp.2013.09.026 [30] J. Vandemeulebroucke, J. Kybic, P. Clarysse, and D. Sarrut, “Respiratory motion estimation from cone-beam projections using a prior model,” in International Conference on Medical Image Computing and Computer-Assisted Intervention, 2009, pp. 365−372: Springer. [31] P. J. Besl and N. D. McKay. A Method for Registration of 3-D Shapes. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1992, 14(2): 239−256 doi: 10.1109/34.121791 [32] M. Baktashmotlagh, M. T. Harandi, B. C. Lovell, and M. Salzmann, “Unsupervised domain adaptation by domain invariant projection,” in Proceedings of the IEEE International Conference on Computer Vision, 2013, pp. 769−776. [33] S. J. Pan, J. T. Kwok, and Q. Yang, “Transfer learning via dimensionality reduction,” in AAAI, 2008, vol. 8, pp. 677-682. [34] S. J. Pan, I. W. Tsang, J. T. Kwok, and Q. Yang. Domain adaptation via transfer component analysis. IEEE Transactions on Neural Networks, 2010, 22(2): 199−210 -
点击查看大图
计量
- 文章访问数: 26
- HTML全文浏览量: 10
- 被引次数: 0
