TomoTherapy Inc., Madison, WI 53717, USA.
Phys Med Biol. 2010 Dec 7;55(23):7211-29. doi: 10.1088/0031-9155/55/23/003. Epub 2010 Nov 16.
IMRT optimization requires a fast yet relatively accurate algorithm to calculate the iteration dose with small memory demand. In this paper, we present a dose calculation algorithm that approaches these goals. By decomposing the infinitesimal pencil beam (IPB) kernel into the central axis (CAX) component and lateral spread function (LSF) and taking the beam's eye view (BEV), we established a non-voxel and non-beamlet-based dose calculation formula. Both LSF and CAX are determined by a commissioning procedure using the collapsed-cone convolution/superposition (CCCS) method as the standard dose engine. The proposed dose calculation involves a 2D convolution of a fluence map with LSF followed by ray tracing based on the CAX lookup table with radiological distance and divergence correction, resulting in complexity of O(N(3)) both spatially and temporally. This simple algorithm is orders of magnitude faster than the CCCS method. Without pre-calculation of beamlets, its implementation is also orders of magnitude smaller than the conventional voxel-based beamlet-superposition (VBS) approach. We compared the presented algorithm with the CCCS method using simulated and clinical cases. The agreement was generally within 3% for a homogeneous phantom and 5% for heterogeneous and clinical cases. Combined with the 'adaptive full dose correction', the algorithm is well suitable for calculating the iteration dose during IMRT optimization.
调强放疗优化需要一个快速且相对准确的算法来计算迭代剂量,同时对内存的需求较小。在本文中,我们提出了一种满足这些目标的剂量计算算法。通过将微射束(IPB)核分解为中心轴(CAX)分量和侧向扩展函数(LSF),并采用束视野(BEV),我们建立了一种无体素和无束块的剂量计算公式。LSF 和 CAX 均由使用坍缩锥卷积/叠加(CCCS)方法作为标准剂量引擎的调试过程确定。所提出的剂量计算涉及用 LSF 对通量图进行二维卷积,然后根据 CAX 查找表进行基于射线追踪的计算,其中包括辐射距离和发散校正,因此在空间和时间上的复杂度均为 O(N(3))。这个简单的算法比 CCCS 方法快几个数量级。由于不需要预先计算束块,因此其实现也比传统的基于体素的束块叠加(VBS)方法小几个数量级。我们使用模拟和临床病例比较了所提出的算法和 CCCS 方法。对于均匀体模,一致性通常在 3%以内,对于非均匀体模和临床病例,一致性在 5%以内。与“自适应全剂量校正”相结合,该算法非常适合在调强放疗优化过程中计算迭代剂量。