Christensen Gary E, Song Joo Hyun, Lu Wei, El Naqa Issam, Low Daniel A
Department of Electrical and Computer Engineering and Department of Radiation Oncology, The University of Iowa, Iowa City, Iowa 52242, USA.
Med Phys. 2007 Jun;34(6):2155-63. doi: 10.1118/1.2731029.
Breathing motion is one of the major limiting factors for reducing dose and irradiation of normal tissue for conventional conformal radiotherapy. This paper describes a relationship between tracking lung motion using spirometry data and image registration of consecutive CT image volumes collected from a multislice CT scanner over multiple breathing periods. Temporal CT sequences from 5 individuals were analyzed in this study. The couch was moved from 11 to 14 different positions to image the entire lung. At each couch position, 15 image volumes were collected over approximately 3 breathing periods. It is assumed that the expansion and contraction of lung tissue can be modeled as an elastic material. Furthermore, it is assumed that the deformation of the lung is small over one-fifth of a breathing period and therefore the motion of the lung can be adequately modeled using a small deformation linear elastic model. The small deformation inverse consistent linear elastic image registration algorithm is therefore well suited for this problem and was used to register consecutive image scans. The pointwise expansion and compression of lung tissue was measured by computing the Jacobian of the transformations used to register the images. The logarithm of the Jacobian was computed so that expansion and compression of the lung were scaled equally. The log-Jacobian was computed at each voxel in the volume to produce a map of the local expansion and compression of the lung during the breathing period. These log-Jacobian images demonstrate that the lung does not expand uniformly during the breathing period, but rather expands and contracts locally at different rates during inhalation and exhalation. The log-Jacobian numbers were averaged over a cross section of the lung to produce an estimate of the average expansion or compression from one time point to the next and compared to the air flow rate measured by spirometry. In four out of five individuals, the average log-Jacobian value and the air flow rate correlated well (R2 = 0.858 on average for the entire lung). The correlation for the fifth individual was not as good (R2 = 0.377 on average for the entire lung) and can be explained by the small variation in tidal volume for this individual. The correlation of the average log-Jacobian value and the air flow rate for images near the diaphragm correlated well in all five individuals (R2 = 0.943 on average). These preliminary results indicate a strong correlation between the expansion/compression of the lung measured by image registration and the air flow rate measured by spirometry. Predicting the location, motion, and compression/expansion of the tumor and normal tissue using image registration and spirometry could have many important benefits for radiotherapy treatment. These benefits include reducing radiation dose to normal tissue, maximizing dose to the tumor, improving patient care, reducing treatment cost, and increasing patient throughput.
呼吸运动是传统适形放疗中减少正常组织剂量和照射的主要限制因素之一。本文描述了利用肺活量测定数据跟踪肺部运动与在多个呼吸周期内从多层CT扫描仪采集的连续CT图像体积的图像配准之间的关系。本研究分析了5名个体的时间CT序列。治疗床从11个不同位置移动到14个不同位置以对整个肺部进行成像。在每个治疗床位置,在大约3个呼吸周期内采集15个图像体积。假定肺组织的扩张和收缩可建模为弹性材料。此外,假定肺在五分之一呼吸周期内的变形较小,因此可使用小变形线性弹性模型对肺的运动进行充分建模。因此,小变形反向一致线性弹性图像配准算法非常适合此问题,并用于配准连续的图像扫描。通过计算用于配准图像的变换的雅可比行列式来测量肺组织的逐点扩张和压缩。计算雅可比行列式的对数,以便肺的扩张和压缩得到同等缩放。在体积中的每个体素处计算对数雅可比行列式,以生成呼吸周期内肺局部扩张和压缩的映射图。这些对数雅可比行列式图像表明,肺在呼吸周期内并非均匀扩张,而是在吸气和呼气期间以不同速率局部扩张和收缩。在肺的一个横截面上对对数雅可比行列式数值求平均值,以得出从一个时间点到下一个时间点的平均扩张或压缩估计值,并与肺活量测定法测量的气流速率进行比较。在5名个体中的4名中,平均对数雅可比行列式值与气流速率相关性良好(整个肺部平均R2 = 0.858)。第5名个体的相关性不太好(整个肺部平均R2 = 0.377),这可以用该个体潮气量变化较小来解释。在所有5名个体中,靠近膈肌的图像的平均对数雅可比行列式值与气流速率的相关性都很好(平均R2 = 0.943)。这些初步结果表明,通过图像配准测量的肺扩张/压缩与通过肺活量测定法测量的气流速率之间存在很强的相关性。利用图像配准和肺活量测定法预测肿瘤和正常组织的位置、运动以及压缩/扩张情况,对于放射治疗可能具有许多重要益处。这些益处包括减少对正常组织的辐射剂量、使肿瘤剂量最大化、改善患者护理、降低治疗成本以及提高患者通量。