www.69555.net利用独立准直器开展调强放疗算法研究戴建荣胡逸民 提要研究利用独立准直器(简称IC)调整射野强度分布。设计了两种递推的办法来计算射野片序列，采用模拟煺火算法优化射野片的照射顺序。并以在Varian 600C加速器上实施三个临床调强射野为例，估计IC调强放疗的照射时间。照射时间一般在5至10min。它近似与机器剂量率成反比，并受调强算法、矩阵元素大小以及强度分级数的影响。 关键词独立准直器算法调强放疗 Intensity-Modulation Radiotherapy Using Independent Collimator:Algorithm Study Dai JianrongHu Yimin Cancer Institute (Hospital),Chinese Academy of Medical Science Peking Union Medical College ABSTRACTTo investigate algorith ms for independent colli mator (IC) to modulate beam intensity distributions.Field-setting sequences wer e determined for three clinical intensity distributions with two forward algorit hms.The delivery order of all fields in the sequences was optimized in order to minimize the overall jaw-moving time.Suppose treatments were performed with Var ian 600C,delivery time was calculated.The delivery time was between 5min to 10mi n,It was approximately inversely proportional to machine dose rate,and also rela ted to algorithm,number of intensity levels and metrix element size. KEY WORDSIndependent collimator AlgorithmIntensity-modulation radiotherapy 适形放疗(Conformal Radiotherapy、简称 CR)是指射线照射形成的剂量分布适合靶区形状，从而可以在减少或不增加周围正常组织受照剂量的同时、增加靶区剂量、达到提高治疗增益比的目的。从理论上讲，实施CR的任何一种方法必须能够调整射野内的射线强度分布，故适形放疗又被称为调强放疗(Intensity_modulationRadiotherapy,简称IMRT)。目前临床开展调强放疗的方法归纳起来有四种1：(1)射束修整器(BeamMo difier),包括挡块、楔形板和补偿器。其缺点是制做费时、摆位困难。(2)断层治疗(Tomotherapy)，以Peacock公司的NOMOS系统为代表2。其优点是将工作人员从制做、搬动挡块和补偿器的体力劳动中解放出来；缺点是治疗时间长，并且由于床步进运动的误差，可能造成薄层相邻区受到超剂量或欠剂量照射。(3)多叶准直器(Multileaf Collimator，简称MLC)，采用分片方式(Segmental mode)3或动态方式(Dynamic mode)调强4-5。MLC具有NOMOS系统的优点，并可缩短治疗时间至数分钟，故成为实施调强放疗的主要方法。其缺点是MLC加工精度高，计算机控制复杂；价格昂贵(约30万美元)；叶片间有漏射线，形成的射野半影大。(4)窄束扫描(Scanning Beam)，典型代表是Scanditronix MM50回旋加速器6。其优点是不仅可以做X线调强，且可调整射线能量以及做电子线调强。但由于价格极其昂贵(约800万美元)，至今全世界仅有约十家放疗中心应用。独立准直器(Independant Collimator，简称IC)已成为国外医用直线加速器准直系统的标准配置，国内北京医疗器械研究所也已研制成功。IC由两对叶片构成，四个电机分别单独驱动每个叶片，可形成射野中心偏移准直器轴线的矩形野。国外已利用IC各叶片独立运动特点实现动态楔形板功能7，这实际是沿一个方向调整射野强度分布。我们希望利用IC调整射野内2维的强度分布，以期建立一种新的既劳动强度轻又经济的调强放疗实施方法，国内外文献尚未见类似报道。本文介绍IC调强放疗的射野片序列算法，并以三个临床调强射野为例估计所需的照射时间。 1.材料与方法 1.1计算IC射野片序列 与MLC类似，IC调强可以采用分片方式，且一组射野相继射来实现所需的强度分布，每个射野片用形成照射范围的准直器叶片位置和照射强度来定义。射野内初始的强度分布可用一个两维矩阵ID0(m,n)表示，其中m和n分别表示矩阵沿准直器X和Y两个方向的大小，下标0表示强度矩阵是初始矩阵。矩阵中每一个元素的强度值可用0和MaxI之间的一个整数值表示，MaxI表示最大强度级。例如，在图1所示的强度矩阵中，mn5，MaxI10。由于每个元素的实际大小和位置不影响计算射野片序列，因此在计算射野片序列时，本文没有直接采用准直器坐标系，而是采用矩阵坐标系，两者之间是平移缩放关系。矩阵坐标系原点位于矩阵元素(0，0)位置，坐标轴刻度是矩阵元素编号。在计算实施照射的计算机控制文件时，每个射野片的准直器叶片位置将转换到准直器坐标系中。 图1随机产生的55强度矩阵，强度分级数为11。 对于任何已知的强度矩阵，可以用许多个IC射野片序列来实现。显然，效率最低的办法是每次只照射一个矩阵元素；效率稍高的一种办法是将IC视为MLC的一对叶片，也就是将所有射野片的宽或长限定为一个元素的宽或长；效率最高的办法应该保证照射时间最短，因此计算射野片序列实际上是一个组合最优化问题。因直接求解很困难，作者设计了两种递推的办法。 方法1基本思路是在强度矩阵中找出最大可以照射的矩形范围作为一个射野片，其计算过程分为三步：(1)搜索初始强度矩阵，找出最大可以照射的矩形区域作为射野片，1，照射强度等于此区域内所有元素的强度最小值。(2)将射野片1照射对ID0(m,n)的强度贡献减去，得到强度矩阵ID1(m,n)。(3)对ID1做第一和第二步处理可求得射野片2的参数；依次类推求其它射野片直至强度矩阵为零。 方法2基本思路是在强度矩阵中找出一个矩形范围作为一个射野以保证剩余的强度矩阵复杂度最低。强度分布的复杂程度可用它所包含的块数来衡量。所谓“块”是指强度大于零且相等的相邻元素所能构成的最大矩形范围，一个非零元素不能同时属于两个或两个以上的块中。例如，在图1所示的强度矩阵中，强度2、5、8、9均有两个相邻元素构成块，而其它元素均单独成块，故此强度矩阵总块数为21。方法2的计算过程分为四步：(1)从初始强度矩阵ID0(m,n)找出一个强度最大的元素。(2)围绕这个元素。根据可使强度矩阵总块数减少最多的原则确定射野片1的叶片位置和照射强度。(3)将射野片1照射对ID0(m,n)的强度贡献减去，得到强度矩阵ID1(m,n)。(4)对ID1(m,n)做第一和第二步处理可求得射野片2的参数；依次类推求其它射野片直至强度矩阵为零。 .2优化IC射野片的照射顺序 当射野片序列确定后，需要优化射野片的照射顺序以减少照射过程中因准直器叶片移动所需的时间。从一个射野片位置变化到下一个射野片位置时准直器四个叶片同时开始移动，但可能不会同时到达各自的指定位置，因此叶片总移动时间TM等于准直器四个叶片移动所需的最长时间之和， 其中N是射野片数目；max是求最大值函数；fabs是求绝对值函数；i-1和i是两个相邻射野片的编号；X1i、X2i、Y1i和Y2i是形成第i个射野片的准直器四个叶片在准直器坐标系中的坐标；Vx、Vy分别是准直器X和Y方向叶片运动速度。 求TM最小值的问题类似与“旅行推销员”问题，后者是指有一个推销员要去N个分别位于(Xi，Yi)的城市进行推销，并于最后返回他出发的城市，要求每个城市只能去一次，而且所经过的路径尽可能短，求推销员途径城市的顺序。这类问题可以用模拟煺火算法求解8。 1.3估计IC调强放疗的照射时间 尽管由于目前还没有实现IC调强放疗的计算机控制功能，IC调强放疗的照射时间不能根据实验确定，但可以根据治疗机参数来估计。射野一次放疗所需时间(即治疗时间)由患者摆位时间和照射时间两部分组成。摆位时间主要取决于治疗部位和是否采用了辅助摆位装置，与具体的调强方法无关，因此可以仅根据照射时间来比较两种IC射野片序列计算方法。对于传统放序，照射时间就是出束时间；对于IC调强放疗，因为采用的是分片技术，照射过程中需要多次中断照射、调整射野位置，因此照射时间TI等于总出束时间TB与准直器叶片总移动时间TM之和， TITBTM(2) (3) 其中Ii是第i个射野片的照射强度；MUPI是单位强度对应的机器跳数；MUR是机器每分钟照射的跳数，即机器剂量率；其它符号参见公式(1)。 剂量率和准直器叶片运动速度与加速器型号有关，同样的强度分布在不同加速器上实施所需的照射时间会有差别，本文以Varian600C加速器为例，计算用两种射野片序列计算方法实施三个临床调强射野所需的照射时间。Varian600C最高剂量率是250MU/min;准直器叶片运动速度卡表测量结果是X方向120cm/min，Y方向72cm/min；请注意，这里未考虑准直器叶片运动范围限制。射野强度分布数据来自纽约MemorialSloanKettering医院的治疗计划系统。系统采用笔束卷积的物理模型计算剂量，采用迭代优化法计算射野强度分布9。因为该院采用MLC动态技术实施调强放疗，计算时笔束Y方向大小只能取MLC叶片宽度1cm，X方向大小取0.2cm。野1和野2是鼻咽癌射野，野3是前列腺癌射野。它们的强度分布如图1所示，在离散为强度矩阵时采用这样一些参数；矩阵元素X方向大小分别取0.2、0.5和1.0cm，Y方向大小固定为1.0cm；最大强度级分别取5、10和20，也就是强度分级数分取6、11和21。 图2(a)鼻咽癌射野强度分布 图2(b)鼻咽癌射野强度分布 图2(c)前列腺癌射野强度分布 2结果 2.1优化射野片照射顺序的作用 经优化后准直器叶片总移动时间明显下降。例如当采用方法1计算，强度矩阵元素大小为0.5cm1.0cm，强度分级数为11时，三个临床调强射野总移动时间从优化前的2至3min减少到优化后的约1min，减少量接近70%(表1)。 表1优化前后准直器叶片总移动时间 (单位：min)。 野号 1 2 3 矩阵大小 1411 188 188 优化前 3.01 2.12 2.56 优化后 0.89 0.66 0.82 减少量(%) 70.43 68.87 67.97 2.2计算方法对照射时间的影响 由表可知，对于同一个强度分布，实现它所需的照射时间采用的计算方法有关。在组成照射时间的两项中，出束时间的贡献是主要的；称动时间贡献小，无论是采用哪种方法，其长短都在1min左右。当强度