Lawson Joshua D, Schreibmann Eduard, Jani Ashesh B, Fox Tim
Department of Radiation Oncology, Emory University, Atlanta, Georgia, U.S.A.
J Appl Clin Med Phys. 2007 Nov 5;8(4):96-113. doi: 10.1120/jacmp.v8i4.2432.
Deformable (non-rigid) registration is an essential tool in both adaptive radiation therapy and image-guided radiation therapy to account for soft-tissue changes during the course of treatment. The evaluation method most commonly used to assess the accuracy of deformable image registration is qualitative human evaluation. Here,we propose a method for systematically measuring the accuracy of an algorithm in recovering artificially introduced deformations in cases of rigid geometry, and we use that method to quantify the ability of a modified basis spline (B-Spline) registration algorithm to recover artificially introduced deformations. The evaluation method is entirely computer-driven and eliminates biased interpretation associated with human evaluation; it can be applied to any chosen method of image registration. Our method involves using planning computed tomography (PCT) images acquired with a conventional CT simulator and cone-beam computed tomography (CBCT) images acquired daily by a linear accelerator-mounted kilovoltage image system in the treatment delivery room. The deformation that occurs between the PCT and daily CBCT images is obtained using a modified version of the B-Spline deformable model designed to overcome the low soft-tissue contrast and the artifacts and distortions observed in CBCT images. Clinical CBCT images and contours of phantom and central nervous system cases were deformed (warped) with known random deformations. In registering the deformed with the non-deformed image sets, we tracked the algorithm's ability to recover the original, non-deformed set. Registration error was measured as the mean and maximum difference between the original and the registered surface contours from outlined structures. Using this approach, two sets of tests can be devised. To measure the residual error related to the optimizer's convergence performance, the warped CBCT image is registered to the unwarped version of itself, eliminating unknown factors such as noise and positioning errors. To study additional errors introduced by artifacts and noise in the CBCT image, the warped CBCT image is registered to the original PCT image. Using a B-Spline deformable image registration algorithm, mean residual error introduced by the algorithm's performance on noise-free images was less than 1 mm, with a maximum of 2 mm. The chosen deformable image registration model was capable of accommodating significant variability in structures over time, because the artificially introduced deformation magnitude did not significantly influence the residual error. On the second type of test, noise and artifacts reduced registration accuracy to a mean of 1.33 mm and a maximum of 4.86 mm.The accuracy of deformable image registration can be easily and consistently measured by evaluating the algorithm's ability to recover artificially introduced deformations in rigid cases in which the true solution is known a priori. The method is completely automated, applicable to any chosen registration algorithm, and does not require user interaction of any kind.
可变形(非刚性)配准是适形放射治疗和图像引导放射治疗中的一项重要工具,用于处理治疗过程中软组织的变化。评估可变形图像配准准确性最常用的方法是定性的人工评估。在此,我们提出一种系统地测量算法在刚性几何情况下恢复人工引入变形准确性的方法,并使用该方法量化改进的基样条(B样条)配准算法恢复人工引入变形的能力。该评估方法完全由计算机驱动,消除了与人工评估相关的有偏差的解释;它可应用于任何选定的图像配准方法。我们的方法涉及使用通过传统CT模拟器获取的计划计算机断层扫描(PCT)图像以及在治疗交付室中由安装在线性加速器上的千伏图像系统每天获取的锥束计算机断层扫描(CBCT)图像。使用为克服CBCT图像中观察到的低软组织对比度、伪影和畸变而设计的B样条可变形模型的改进版本,来获取PCT图像和每日CBCT图像之间发生的变形。对临床CBCT图像以及体模和中枢神经系统病例的轮廓施加已知的随机变形(扭曲)。在将变形后的图像集与未变形的图像集进行配准时,我们跟踪算法恢复原始未变形集的能力。配准误差通过轮廓结构的原始表面轮廓和配准后的表面轮廓之间的平均差异和最大差异来衡量。使用这种方法,可以设计两组测试。为了测量与优化器收敛性能相关的残余误差,将扭曲的CBCT图像与自身的未扭曲版本进行配准,消除噪声和定位误差等未知因素。为了研究CBCT图像中的伪影和噪声引入的额外误差,将扭曲的CBCT图像与原始PCT图像进行配准。使用B样条可变形图像配准算法,该算法在无噪声图像上的性能引入的平均残余误差小于1毫米,最大为2毫米。所选的可变形图像配准模型能够适应结构随时间的显著变化,因为人工引入的变形幅度并未显著影响残余误差。在第二类测试中,噪声和伪影将配准精度降低到平均1.33毫米,最大4.86毫米。通过评估算法在已知先验真实解的刚性情况下恢复人工引入变形的能力,可以轻松且一致地测量可变形图像配准的准确性。该方法完全自动化,适用于任何选定的配准算法,并且不需要任何类型的用户交互。