一种术前剂量验证方法、装置及放射治疗设备与流程

1.本发明属于放射治疗技术领域，具体涉及一种术前剂量验证方法、装置及放射治疗设备。


背景技术：

2.目前剂量验证通过水箱，二维矩阵(电离室或者半导体探头)或者胶片等验证工具进行。水箱虽然能较为精确反应绝对剂量的分布，但是架水箱的步骤繁琐耗时，不利于日常计划的qa工作。二维矩阵虽然操作简单，但是受限于精度限制，不能很好的反映计划中迅速跌落的剂量梯度。胶片虽然空间精度较高，但是后续的处理较为复杂，线性响应范围较窄，也不适合目前的治疗的验证工作。其次，目前三维验证设备在旋转调强的验证中，需要额外的设备和工具来保证和机架旋转的同步。而epid由于其安装，天然满足和机架的同步，且操作简单(几乎不需要额外的工作)精度高，是理想的加速器自带的验证设备。目前，基于epid的剂量验证主要包括以下几类：基于光栅的logfile重建，基于通量(透射)和剂量算法重建，基于测量进行重建。基于logfile的方法，在不同加速器重可能会遇到权限问题，因此有通用性问题(对于不同的厂商不一定能获取该数据)。基于测量重建的方法，虽然能标定每个探头对绝对(相对)剂量的响应，但是epid探头本身的不稳定和标定的过程复杂而且漫长，且不同加速器之间可能需要重新标定(架水箱)和不断的校准，因此成本较高。基于算法的重建，对剂量算法本身的精确性提出了更高的要求，其中各项参数拟合过程复杂且繁琐。


技术实现要素：

3.为了解决上述技术问题，本发明提出了一种术前剂量验证方法、装置及放射治疗设备。
4.为了达到上述目的，本发明的技术方案如下：
5.一方面，本发明公开一种术前剂量验证方法，具体包括以下步骤：
6.s1：epid连续采图，在所有图像中提取有效图像；
7.s2：对s1得到的每个有效图像进行图像修正；
8.s3：对经过s2修正后的每个图像进行标定，叠加处理后的所有图像和权重修正，得到修正的通量图；
9.s4：通过s3得到修正的通量图进行蒙特卡洛剂量计算，得到重建的剂量分布；
10.s5：将s4得到的重建的剂量分布与tps计算比对，根据gamma分析计算通过率，输出验证结果。
11.在上述技术方案的基础上，还可做如下改进：
12.作为优选的方案，s1中提取有效图像通过以下内容实现：
13.对每个图像的所有像素值从大到小排序确定每个图像的初始最大值max1；
14.若最大值max1大于一定像素阈值y1，则该图像为有效图像。
15.作为优选的方案，s2中，对s1得到的每个有效图像进行去坏点和暗电流处理、去条
纹处理及均匀处理中的一种或多种。
16.作为优选的方案，s3具体包括以下步骤：
17.s3.1：对经过s2修正后的每个图像进行大小野标定和离轴标定；
18.s3.2：对每个图像，以最大值max1为基准，对整个平板的像素数据逐点计算分析，建立主射线区域的像素映射；
19.s3.3：根据主射线区域计算出该区域的面积、周长和离轴位置；
20.s3.4：根据周长和面积计算出大小野权重因子，该因子和最大值max1值相乘，进行权重修正；
21.s3.5：根据离轴位置计算出离轴因子，该因子和最大值max1值相乘，进行权重修正；
22.s3.6：叠加处理后的所有图像，进行mu标定；
23.s3.7：图像进行分辨率的调整和裁剪，得到修正的通量图。
24.作为优选的方案，s3.1中大小野标定具体包括以下内容：
25.采集不同开野大小在同一跳数下epid的连续图像，对中心处像素累加和平均，然后用像素值和开野大小做函数拟合进行大小野的标定；
26.离轴标定具体包括以下内容：
27.采集其中一开野下，不同离轴距离下epid的连续图像，对中心处像素累加和平均，然后用像素值和对应离轴距离做函数拟合进行离轴标定；
28.mu标定具体包括以下内容：
29.采集其中一开野下，不同跳数mu下epid的连续图像，对中心处像素累加和平均，然后用像素值和对应的跳数mu做函数拟合进行mu标定。
30.作为优选的方案，s4具体包括以下内容：
31.根据修正的通量图在均匀水模体中进行蒙特卡洛剂量计算，得到重建的剂量分布。
32.作为优选的方案，s5具体包括以下内容：
33.提取重建的剂量分布中任一层的剂量作为最终的二维验证结果和tps在均匀水模体中的计算进行比对，根据gamma分析计算通过率，输出验证结果。
34.另一方面，本发明还公开一种术前剂量验证装置，包括：
35.提取模块，用于在epid连续采图后的所有图像中提取有效图像；
36.图像修正模块，用于对提取模块提取的每个有效图像进行图像修正；
37.图像转通量模块，用于对图像修正模块修正后的每个图像进行标定，叠加处理后的所有图像和权重修正，得到修正的通量图；
38.剂量分布模块，用于对图像转通量模块得到的修正的通量图进行蒙特卡洛剂量计算，得到重建的剂量分布；
39.验证模块，用于将剂量分布模块得到的重建的剂量分布与tps计算比对，根据gamma分析计算通过率，输出验证结果。
40.作为优选的方案，提取模块通过以下内容提取有效图像；
41.对每个图像的所有像素值从大到小排序确定每个图像的初始最大值max1；
42.若最大值max1大于一定像素阈值y1，则该图像为有效图像。
43.另一方面，本发明公开一种放射治疗设备，利用上述任一种术前剂量验证方法进行验证；或，包括上述任一种术前剂量验证装置。
44.本发明一种术前剂量验证方法、装置及放射治疗设备，具有以下有益效果：
45.第一，采用蒙特卡洛剂量计算作为像素到剂量的最主要映射手段，准确性高。
46.第二，通过逐帧对每个图像进行修正处理，最大限度避免了由于采集过程导致的细节损失。
47.第三，通过水箱等验证工具以及epid本身的采图特性，进行标定，最大限度保证了重建的精确性。
48.第四，基于逐帧标定和修正的连续采图方式减小了对采集的数据的数量要求、降低了对epid探头稳定性的要求和一定程度上降低了剂量算法的复杂性。
附图说明
49.为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
50.图1为本发明实施例提供的术前剂量验证方法的流程图。
51.图2为本发明实施例提供的术前剂量验证方法的流程图(具有预步骤)。
具体实施方式
52.下面结合附图详细说明本发明的优选实施方式。
53.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0054]“包括”元件的表述是“开放式”表述，该“开放式”表述仅仅是指存在对应的部件，不应当解释为排除附加的部件。
[0055]
为了达到本发明的目的，一种术前剂量验证方法、装置及放射治疗设备的其中一些实施例中，术前剂量验证方法具体包括以下步骤：
[0056]
s1：epid连续采图，在所有图像中提取有效图像；
[0057]
s2：对s1得到的每个有效图像进行图像修正；
[0058]
s3：对经过s2修正后的每个图像进行标定，叠加处理后的所有图像和权重修正，得到修正的通量图；
[0059]
s4：通过s3得到修正的通量图进行蒙特卡洛剂量计算，得到重建的剂量分布；
[0060]
s5：将s4得到的重建的剂量分布与tps计算比对，根据gamma分析计算通过率，输出验证结果。
[0061]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，s1中提取有效图像通过以下内容实现：
[0062]
对每个图像的所有像素值从大到小排序确定每个图像的初始最大值max1；
[0063]
若最大值max1大于一定像素阈值y1，则该图像为有效图像。
[0064]
根据探头的响应，设计算法确定有效采集的开始和终止，确保程序积累有效的像素信息，对每帧图像的所有像素值从大到小排序(一般选取第若干个最大值)确定每一帧图像的初始最大值max1，选取的数目可以根据采图的像素决定。
[0065]
比如1024x1024像素，一般选取第1024个作为最大值max1。最大值max1需要大于一定像素阈值y1，则作为有效图像。
[0066]
用类似的方法，选取第二个最大值max2，如果max2大于一定像素阈值y2，则该图像为稳定图像数据，否则为过渡图像数据。
[0067]
对于有效图像数据的不同阶段，调整部分设置的参数。
[0068]
进一步，在上述实施例的基础上，s2中，对s1得到的每个有效图像进行去坏点和暗电流处理、去条纹处理及均匀处理。
[0069]
对epid采集的图像进行影像学的图像修正，该图像修正需要对采集的每帧图像进行单独的修正。
[0070]
根据采集的暗场图像，确定坏点和暗电流的影像，建立暗场的坏点数据和暗电流数据，记录到暗场测试文件；在实际操作中，对每个点的像素值和周围的平均值做比较，若果大于一定比例，那么确定为坏点，记录其像素坐标(i，j)；对整个暗场做像素点的平均，作为本次暗场的基准值。
[0071]
根据采集的图像，设计消除条纹和不均性的算法，最大限度消除采图本身的影响。
[0072]
消除条纹：以加速器床板的长度方向为平板epid的y方向，对于每个像素点做y方向值卷积两次，分别为包含该点和不包含该点，记为value1和value2。如果value1和value2的值差距过大，那么该点的像素值将会被取代。
[0073]
不均性处理：选取一定的卷积板对数据做平滑或者锐化处理。
[0074]
将上述方法写入数据的预处理类。
[0075]
进一步，在上述实施例的基础上，s3具体包括以下步骤：
[0076]
s3.1：对经过s2修正后的每个图像进行大小野标定和离轴标定；
[0077]
s3.2：对每个图像，以最大值max1为基准，对整个平板的像素数据逐点计算分析，建立主射线区域的像素映射；
[0078]
s3.3：根据主射线区域计算出该区域的面积、周长和离轴位置；
[0079]
s3.4：根据周长和面积计算出大小野权重因子，该因子和最大值max1值相乘，进行权重修正；
[0080]
s3.5：根据离轴位置计算出离轴因子，该因子和最大值max1值相乘，进行权重修正；
[0081]
s3.6：叠加处理后的所有图像，进行mu标定；
[0082]
s3.7：图像进行分辨率的调整和裁剪，得到修正的通量图。
[0083]
进一步，在上述实施例的基础上，s3.1中大小野标定具体包括以下内容：
[0084]
采集不同开野大小，取3x3、5x5、8x8、10x10、12x12、15x15、20x20开野在同一跳数(100mu或者其他大跳数)下epid的连续图像，对中心处像素累加和平均，然后用像素值和开野大小做函数拟合进行大小野的标定。
[0085]
在程序中，需要粗略计算每一帧野的面积和周长。该大小野标定在累加图像前，逐
帧修正。
[0086]
离轴标定具体包括以下内容：
[0087]
采集10x10开野下，不同离轴距离，取0cm、2cm、4cm、6cm、8cm下epid的连续图像，对中心处像素累加和平均，然后用像素值和对应离轴距离做函数拟合进行离轴标定。
[0088]
该离轴标定在累加图像前，逐帧修正。
[0089]
mu标定具体包括以下内容：
[0090]
采集10x10开野下，不同跳数mu，取2mu、5mu、10mu、20mu、50mu、100mu、200mu、300mu下epid的连续图像，对中心处像素累加和平均，然后用像素值和对应的跳数mu做函数拟合进行mu标定。
[0091]
该mu标定在累加图像后，非逐帧修正。
[0092]
进一步，在上述实施例的基础上，s3.3中，根据最大值max1为基准，对于整个平板的像素数据逐点计算分析，建立主射线区域、半影区和非透射区域的像素映射。此外，阈值的响应标准可以通过参数调整，为了得到更好的验证效果。整个阈值响应通过逐帧处理。
[0093]
进一步，在上述实施例的基础上，s4具体包括以下内容：
[0094]
根据修正的通量图在均匀水模体中进行蒙特卡洛剂量计算，得到重建的剂量分布。
[0095]
进一步，在上述实施例的基础上，s5具体包括以下内容：
[0096]
提取重建的剂量分布中任一层的剂量作为最终的二维验证结果和tps在均匀水模体中的计算进行比对，根据gamma分析计算通过率，输出验证结果。
[0097]
本发明采集规则野和若干计划的epid连续图像和实测剂量，作为检验和测试。
[0098]
为了保证本发明的有效性，本发明还具有以下预步骤：
[0099]
s01：将epid采图计算和规则野实测数据比对，测试需要达到一定的标准，3310需要达到100％，225需要达到95％以上，作为加速器和epid稳定性的保证。
[0100]
s02：测试暗场数据，建立本次测量的坏点数据和暗电流数据。
[0101]
本发明实施例还公开一种术前剂量验证装置，包括：
[0102]
提取模块，用于在epid连续采图后的所有图像中提取有效图像；
[0103]
图像修正模块，用于对提取模块提取的每个有效图像进行图像修正；
[0104]
图像转通量模块，用于对图像修正模块修正后的每个图像进行标定，叠加处理后的所有图像和权重修正，得到修正的通量图；
[0105]
剂量分布模块，用于对图像转通量模块得到的修正的通量图进行蒙特卡洛剂量计算，得到重建的剂量分布；
[0106]
验证模块，用于将剂量分布模块得到的重建的剂量分布与tps计算比对，根据gamma分析计算通过率，输出验证结果。
[0107]
进一步，在上述实施例的基础上，提取模块通过以下内容提取有效图像；
[0108]
对每个图像的所有像素值从大到小排序确定每个图像的初始最大值max1；
[0109]
若最大值max1大于一定像素阈值y1，则该图像为有效图像。
[0110]
进一步，在上述实施例的基础上，图像修正模块对得到的每个有效图像进行去坏点和暗电流处理、去条纹处理及均匀处理。
[0111]
进一步，在上述实施例的基础上，图像转通量模块通过以下步骤得到修正的通量
图：
[0112]
s3.1：对经过s2修正后的每个图像进行大小野标定和离轴标定；
[0113]
s3.2：对每个图像，以最大值max1为基准，对整个平板的像素数据逐点计算分析，建立主射线区域的像素映射；
[0114]
s3.3：根据主射线区域计算出该区域的面积、周长和离轴位置；
[0115]
s3.4：根据周长和面积计算出大小野权重因子，该因子和最大值max1值相乘，进行权重修正；
[0116]
s3.5：根据离轴位置计算出离轴因子，该因子和最大值max1值相乘，进行权重修正；
[0117]
s3.6：叠加处理后的所有图像，进行mu标定；
[0118]
s3.7：图像进行分辨率的调整和裁剪，得到修正的通量图。
[0119]
进一步，在上述实施例的基础上，s3.1中大小野标定具体包括以下内容：
[0120]
采集不同开野大小在同一跳数下epid的连续图像，对中心处像素累加和平均，然后用像素值和开野大小做函数拟合进行大小野的标定；
[0121]
离轴标定具体包括以下内容：
[0122]
采集其中一开野下，不同离轴距离下epid的连续图像，对中心处像素累加和平均，然后用像素值和对应离轴距离做函数拟合进行离轴标定；
[0123]
mu标定具体包括以下内容：
[0124]
采集其中一开野下，不同跳数mu下epid的连续图像，对中心处像素累加和平均，然后用像素值和对应的跳数mu做函数拟合进行mu标定。
[0125]
进一步，在上述实施例的基础上，剂量分布模块通过以下内容得到重建的剂量分布：
[0126]
根据修正的通量图在均匀水模体中进行蒙特卡洛剂量计算，得到重建的剂量分布。
[0127]
进一步，在上述实施例的基础上，验证模块通过以下内容进行验证：
[0128]
提取重建的剂量分布中任一层的剂量作为最终的二维验证结果和tps在均匀水模体中的计算进行比对，根据gamma分析计算通过率，输出验证结果。
[0129]
本发明实施例还公开一种放射治疗设备，利用上述任一实施例公开的术前剂量验证方法进行验证；或，包括上述任一实施例公开的术前剂量验证装置。
[0130]
本发明根据蒙特洛剂量计算、epid影像分析技术和基于epid采图的剂量标定进行术前的剂量验证工作。在整个放射治疗之前，执行预定的计划，并在此期间，通过epid进行连续采图，通过逐帧对每个图像进行标定和修正处理，避免采集导致的细节损失，减少对剂量算法的依赖，重建的精确性高。基于逐帧标定和修正的连续采图和蒙卡计算相结合的剂量重建方式一方面和像素标定的重建方法相比，减小了对采集的数据量要求和降低对epid探头稳定性的要求，另一方面和剂量算法为基础的重建方法相比，降低了剂量计算的复杂度。该方法可以用在静态、动态和旋转调强的术前的剂量验证工作。
[0131]
本发明一种术前剂量验证方法、装置及放射治疗设备，具有以下有益效果：
[0132]
第一，采用蒙特卡洛剂量计算作为像素到剂量的最主要映射手段，准确性高。
[0133]
第二，通过逐帧对每个图像进行修正处理，最大限度避免了由于采集过程导致的
细节损失。
[0134]
第三，通过水箱等验证工具以及epid本身的采图特性，进行标定，最大限度保证了重建的精确性。
[0135]
第四，基于逐帧标定和修正的连续采图方式减小了对采集的数据的数量要求、降低了对epid探头稳定性的要求和一定程度上降低了剂量算法的复杂性。
[0136]
应当理解，上述描述的各种技术可结合硬件或软件，或者它们的组合一起实现。从而，本发明的方法和设备，或者本发明的方法和设备的某些方面或部分可采取嵌入有形媒介，例如软盘、cd
‑
rom、硬盘驱动器或者其它任意机器可读的存储介质中的程序代码(即指令)的形式，其中当程序被载入诸如计算机之类的机器，并被该机器执行时，该机器变成实践本发明的设备。
[0137]
上述实施例只为说明本发明的技术构思及特点，其目的在于让本领域普通技术人员能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围内。