CN118416412B - 探测器阵列校正方法、模块、磁共振加速器及存储介质 - Google Patents
探测器阵列校正方法、模块、磁共振加速器及存储介质 Download PDFInfo
- Publication number
- CN118416412B CN118416412B CN202410885122.7A CN202410885122A CN118416412B CN 118416412 B CN118416412 B CN 118416412B CN 202410885122 A CN202410885122 A CN 202410885122A CN 118416412 B CN118416412 B CN 118416412B
- Authority
- CN
- China
- Prior art keywords
- detector array
- data
- gold standard
- accelerator
- correction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 145
- 238000000034 method Methods 0.000 title claims abstract description 101
- 230000005855 radiation Effects 0.000 claims abstract description 21
- 230000011218 segmentation Effects 0.000 claims abstract description 8
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 68
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000000523 sample Substances 0.000 claims description 7
- BTCSSZJGUNDROE-UHFFFAOYSA-N gamma-aminobutyric acid Chemical compound NCCCC(O)=O BTCSSZJGUNDROE-UHFFFAOYSA-N 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 2
- 238000001959 radiotherapy Methods 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 4
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000007408 cone-beam computed tomography Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000275 quality assurance Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 235000003332 Ilex aquifolium Nutrition 0.000 description 1
- 241000209027 Ilex aquifolium Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 231100000673 dose–response relationship Toxicity 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011347 external beam therapy Methods 0.000 description 1
- 238000002786 image-guided radiation therapy Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
本发明属于放疗设备技术领域，具体涉及一种探测器阵列校正方法、模块、磁共振加速器及存储介质。本发明提供进行磁共振兼容探测器阵列校正的方法，通过判断加速器最大射野大小是否可覆盖探测器阵列，分别采用基于中心位置的单野校正方法，或基于偏中心位置的多野分段校正方法。本发明还进一步提供了实现该方法的模块及集成有该模块的磁共振加速器。本发明具有简单和准确的优势，能够更好地改善磁共振兼容探测器阵列在磁共振加速器中的应用效果，精确测量加速器束流轮廓数据，从而保证QA的精确性，提高工作效率，具有很好的应用前景。
Description
技术领域
本发明属于放疗设备技术领域，具体涉及一种探测器阵列校正方法、模块、磁共振加速器及存储介质。
背景技术
放射治疗，简称放疗，是治疗肿瘤的主要手段之一，而直线加速器是目前放疗中最常用的设备。为了保证治疗的精确性，常规加速器通常配备锥形束CT（Cone Beam ComputedTomography，CBCT）、电子射野影像装置（Electronic Portal Imaging Device，EPID）和兆伏CT（Megavoltage Computed Tomography）等设备用于图像引导放射治疗。磁共振加速器是将磁共振成像（Magnetic Resonance Imaging，MRI）与直线加速器的功能整合与一体的革命性技术，它在精准放疗领域引起了极大关注，其优势在于将MRI提供的高软组织对比度、实时成像等优点与加速器的精确剂量传递相结合，从而实现更高水平的放射剂量传递与实时引导治疗等目标。
高精度的放射治疗设备更需要高精度的质量保证（Quality Assurance，QA），在常规加速器的QA中，二维探测器阵列的应用非常广泛，其在加速器束流调试验证与患者计划验证中发挥着重要作用。二维探测器阵列由于其易用性与精确性，已经在加速器束流轮廓（profile）调试与日常QA中得到了广泛应用，甚至在部分情况下可以代替传统的水箱扫描系统进行加速器调试及建模数据采集。
然而，磁共振加速器系统中的高强度磁场对于加速器束流性能造成了本质影响，光子束流与介质相互作用产生的次级电子会受到磁场的影响，从而导致探测器阵列的剂量响应在磁场中发生变化。
因此，需要针对磁共振兼容的探测器阵列进行校正。阵列校正因子通常由每个探测器的灵敏度数值组成，这些灵敏度数值可以应用于探测器的输出测量，从而使得校正后的剂量分布独立于探测器的敏感度。为了避免探测器的能量与角度依赖性，目前磁共振兼容探测器阵列的校正方法通常是为磁共振加速器修改后的基于宽野的校正方法（WideField Calibration）。该方法中，整个阵列校正的过程包括了10次照射曝光，每次曝光前均需将探测器阵列放置于射野中不同的位置，且每次曝光均使用40 cm×22 cm的射野与200MU的照射剂量，每一步骤均被记录为独立的测量。最终，将十组测量文件进行分析处理，生成相应的校正文件。该方法目前主要存在如下问题：1.使用该方法校正后的束流轮廓（profile）数据与水箱测量得到的金标准数据仍然存在不可忽略的偏差；2.使用该方法校正后的束流轮廓数据无法反映加速器束流真实情况，会引起临床QA结果的误判；3.使用该方法校正后的束流轮廓数据存在异常区域（异常点、曲线不平滑），影响束流特征指标（如对称性）的准确计算；4.该方法的校正过程较为繁琐，重复性差，增加临床工作量。
因此，本领域仍然亟需开发更加简单和准确的磁共振兼容探测器阵列校正方法。
发明内容
针对现有技术的问题，本发明提供一种探测器阵列校正方法、模块、磁共振加速器及存储介质。
一种探测器阵列校正方法，加速器最大射野大小无法覆盖探测器阵列时，进行探测器阵列校正，包括如下步骤：
步骤1，将探测器阵列中心与射野中心对齐，将探测器阵列设置在与金标准数据相同深度处，收集加速器最大射野的束流轮廓数据，得到中心位置处探测器阵列数据；
步骤2，计算金标准数据与所述中心位置处探测器阵列数据的比值，得到中心位置处探测器阵列的校正因子Ic(r)，其中，r是指探测器阵列中的探测器；
步骤3，将探测器阵列维持在与金标准数据相同深度处，以所述中心位置为起点，设计探测器阵列的至少一条移动路线，所述探测器阵列在一条移动路线上移动n次，n为大于等于1的整数，不同移动路线上n的取值相同或不同，经过所有移动后探测器阵列的所有探测器至少一次位于加速器最大射野范围内；
在同一移动路线中，对于第j次移动，保证至少有一个共同的探测器rj，所述探测器rj满足在第j次移动前后，均位于加速器最大射野范围内；收集每次移动后加速器最大射野的束流轮廓数据，得到n组偏中心位置处探测器阵列数据；其中，j的取值为1至n的整数；
步骤4，在同一移动路线中，计算金标准数据与所述偏中心位置处探测器阵列数据的比值，得到偏中心位置处探测器阵列的校正因子ij(r)；
步骤5，在同一移动路线中，选择步骤3中所述共同的探测器rj作为参考点，按照如下公式计算修正因子uj：
当j=1，uj=Ic(rj)/ij(rj)；
当j大于1，uj= ij-1(rj)uj-1/ij(rj)；
步骤6，按照如下公式对偏中心位置处探测器阵列的校正因子进行修正：
Ij(r)=ij(r)×uj(r)，其中，Ij(r)即为修正后的偏中心位置校正因子；
步骤7，将Ic(r)和所有移动路线中得到的Ij(r)进行组合，得到最终的探测器阵列校正因子，将所述最终的探测器阵列校正因子替换到原始校正文件中。
优选的，步骤2和步骤4中，所述金标准数据通过水箱测量法采集。
优选的，步骤2中，所述金标准数据使用前进行如下处理：
步骤a1，将所述金标准数据与所述中心位置处探测器阵列数据按照射野中心轴进行探头数据归一；
步骤a2，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据；
步骤4中，所述金标准数据使用前进行如下处理：
步骤b1，根据步骤3中的移动距离调整所述金标准数据的坐标；
步骤b2，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据；
步骤b3，将所述金标准数据与所述偏中心位置处探测器阵列数据按照射野中心轴进行探头数据归一。
优选的，步骤3中，保证每一次移动后探测器阵列的中心不在射野区域内。
优选的，步骤7中，当同一探测器在不同移动路线中计算得到多个Ij(r)时，则最终取值按照如下方式中的一种进行选择:
1)选择初次获得的值；
2)任意选取一个取值；
3)计算平均值。
本发明还提供一种探测器阵列校正方法，进行探测器阵列校正，包括如下步骤：
步骤A，获取加速器束流轮廓金标准数据；
步骤B，判断加速器最大射野大小是否可覆盖探测器阵列；
若无法覆盖则按照权利要求1-5任一项所述的探测器阵列校正方法执行，进行基于偏中心位置的多野分段校正方法；
若可以覆盖则进行基于中心位置的单野校正方法，所述基于中心位置的单野校正方法按照如下步骤执行：
步骤C 1，将探测器阵列中心与射野中心对齐，将探测器阵列设置在与金标准数据相同深度处，收集加速器最大射野的束流轮廓数据，得到中心位置处探测器阵列数据；
步骤C2，计算金标准数据与所述中心位置处探测器阵列数据的比值，得到校正因子，将所得校正因子替换到原始校正文件中。
优选的，步骤C2中，所述金标准数据使用前进行如下处理：
步骤C21，将所述金标准数据与所述中心位置处探测器阵列数据按照射野中心轴进行探头数据归一；
步骤C22，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据。
本发明还提供一种探测器阵列校正模块，包括：
控制器，被配置为执行上述探测器阵列校正方法；
存储器，被配置为存储程序代码和中间数据；
数据接口，被配置为与外部设备连接，执行数据的输入和输出。
本发明还提供一种磁共振加速器，包括：放射源、加速器组件、探测器阵列和控制系统；所述控制系统中集成有上述探测器阵列校正模块。
本发明还提供一种计算机可读存储介质，其上存储有：用于实现上述探测器阵列校正方法的计算机程序。
采用本发明提供的方法、模块或装置进行磁共振兼容探测器阵列校正，具有如下有益的技术效果：
1、校正后的束流轮廓（profile）数据与水箱测量得到的金标准数据一致性好、吻合度高，可以反映加速器束流真实情况。
2、校正后的束流轮廓中无异常区域，且校正方法简单，可重复性高。
3、本发明方法具有非常高的普适性，适用于探测器阵列的校正，尤其是现有方法存在重大缺陷的情况，可推广到一维、二维、三维甚至多维探测器阵列，且在有磁场和无磁场下均适用。
综上所述，本发明提供的方法具有简单和准确的优势，具有很好的应用前景。
显然，根据本发明的上述内容，按照本领域的普通技术知识和惯用手段，在不脱离本发明上述基本技术思想前提下，还可以做出其它多种形式的修改、替换或变更。
以下通过实施例形式的具体实施方式，对本发明的上述内容再作进一步的详细说明。但不应将此理解为本发明上述主题的范围仅限于以下的实例。凡基于本发明上述内容所实现的技术均属于本发明的范围。
附图说明
图1（a）为基于宽野的校正方法的校正结果中X轴束流轮廓数据比较结果图；
图1（b）为基于中心位置的单野校正方法的校正结果中X轴束流轮廓数据比较结果图；
图2（a）为基于宽野的校正方法的校正结果中Y轴束流轮廓数据比较结果图；
图2（b）为基于中心位置的单野校正方法的校正结果中Y轴束流轮廓数据比较结果图；
图3（a）为基于宽野的校正方法的校正结果中正对角线轴束流轮廓数据比较结果图；
图3（b）为基于中心位置的单野校正方法的校正结果中正对角线轴束流轮廓数据比较结果图；
图4（a）为基于宽野的校正方法的校正结果中负对角线轴束流轮廓数据比较结果图；
图4（b）为基于中心位置的单野校正方法的校正结果中负对角线轴束流轮廓数据比较结果图；
图5（a）为基于宽野的校正方法的校正结果中X轴束流轮廓数据比较结果图；
图5（b）为基于偏中心位置的多野分段校正方法的校正结果中X轴束流轮廓数据比较结果图；
图6（a）为基于宽野的校正方法的校正结果中Y轴束流轮廓数据比较结果图；
图6（b）为基于偏中心位置的多野分段校正方法的校正结果中Y轴束流轮廓数据比较结果图；
图7（a）为基于宽野的校正方法的校正结果中正对角线轴束流轮廓数据比较结果图；
图7（b）为基于偏中心位置的多野分段校正方法的校正结果中正对角线轴束流轮廓数据比较结果图；
图8（a）为基于宽野的校正方法的校正结果中负对角线轴束流轮廓数据比较结果图；
图8（b）为基于偏中心位置的多野分段校正方法的校正结果中负对角线轴束流轮廓数据比较结果图。
具体实施方式
需要特别说明的是，实施例中未具体说明的数据采集、传输、储存和处理等步骤的算法，以及未具体说明的硬件结构、电路连接等均可通过现有技术已公开的内容实现。
以下实施例和实验例中，部分术语的含义如下：
“磁共振加速器”是将磁共振成像系统与直线加速器整合一体的设备，该设备基于磁共振成像引导，使用X射线进行放射治疗，因此设备中存在高强度磁场；此外由于高强度磁场的存在，需要应用专门的磁共振兼容探测器对其进行QA。
“加速器束流轮廓”用于描述外放射治疗设备加速器产生的X线束流特征，指不同方向（通常包括横向（IEC-X轴）、纵向（IEC-Y轴）、正负对角线轴等）、不同深度的X线剂量分布形状，通常可使用探测器阵列、三维水箱扫描系统测量得到。
“射野”代表直线加速器照射范围的形状、大小。
“射野区域”定义为实际射野大小80%以内的区域。
“最大射野”定义为加速器在准直器（限束系统）开到最大时的照射野。
“插值”是一种利用函数对数据范围内的某一点的数值进行预测的常用方法。
“伽马通过率”指在一定剂量阈值和伽马分析的标准下，绝对或者相对剂量的通过率，用于描述测量剂量与标准剂量之间的差异。
实施例1 磁共振兼容探测器阵列校正方法
本实施例的方法包括如下步骤：
步骤A，将探测器阵列的校正因子重置为1；获取磁共振加速器束流轮廓金标准数据；所述金标准数据可采用水箱测量法采集；
步骤B，判断加速器最大射野大小是否可覆盖探测器阵列；
若无法覆盖则进行基于偏中心位置的多野分段校正方法；
若可以覆盖则进行基于中心位置的单野校正方法。
其中，基于偏中心位置的多野分段校正方法的具体步骤如下：
步骤1，将探测器阵列中心与射野中心对齐，将探测器阵列设置在与金标准数据相同深度处，收集加速器最大射野的束流轮廓数据，得到中心位置处探测器阵列数据。
步骤2，计算金标准数据与所述中心位置处探测器阵列数据的比值，得到中心位置处探测器阵列的校正因子Ic(r)，其中，r是指探测器阵列中的探测器。
上述步骤中，所述金标准数据使用前进行如下处理：
步骤a1，将所述金标准数据与所述中心位置处探测器阵列数据按照射野中心轴进行探头数据归一；
步骤a2，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据。
步骤3，
将探测器阵列维持在与金标准数据相同深度处，以所述中心位置为起点，设计探测器阵列的至少一条移动路线，所述探测器阵列在一条移动路线上移动n次，n为大于等于1的整数，不同移动路线上n的取值相同或不同。
对于移动路线的数量、移动路线的移动方向、同一移动路线内的移动次数n均可根据探测器阵列的形状和探测器分布进行适应性调整。目的在于，经过所有移动后探测器阵列的所有探测器至少一次位于加速器最大射野范围内。
在同一移动路线中，对于第j次移动，保证至少有一个共同的探测器rj，所述探测器rj满足在第j次移动前后，均位于加速器最大射野范围内；收集每次移动后加速器最大射野的束流轮廓数据，得到n组偏中心位置处探测器阵列数据；其中，j的取值为1至n的整数。
步骤4，在同一移动路线中，计算金标准数据与所述偏中心位置处探测器阵列数据的比值，得到偏中心位置处探测器阵列的校正因子ij(r)。
上述步骤中，所述金标准数据使用前进行如下处理：
步骤b1，根据步骤3中的移动距离调整所述金标准数据的坐标；
步骤b2，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据；
步骤b3，将所述金标准数据与所述偏中心位置处探测器阵列数据按照射野中心轴进行探头数据归一。
步骤5，在同一移动路线中，选择步骤3中所述共同的探测器rj作为参考点，按照如下公式计算修正因子uj：
当j=1，uj=Ic(rj)/ij(rj)；
当j大于1，uj= ij-1(rj) uj-1/ij(rj)。
步骤6，按照如下公式对偏中心位置处探测器阵列的校正因子进行修正：
Ij(r)=ij(r)×uj(r)，其中，Ij(r)即为修正后的偏中心位置校正因子。
步骤7，将Ic(r)和所有移动路线中得到的Ij(r)进行组合，得到最终的探测器阵列校正因子，将所述最终的探测器阵列校正因子替换到原始校正文件中。
上述步骤中，对于一些实施例，可能出现同一探测器在不同的移动路线中多次被最大射野覆盖的情况。当同一探测器在不同移动路线中计算得到多个Ij(r)时，考虑到越靠近中心位置进行修正的结果越可靠，因此本实施例中最终的Ij(r)取值选取初次修正的值。
基于中心位置的单野校正方法的具体步骤如下：
步骤C 1，将探测器阵列中心与射野中心对齐，将探测器阵列设置在与金标准数据相同深度处，收集加速器最大射野的束流轮廓数据，得到中心位置处探测器阵列数据；
步骤C2，计算金标准数据与所述中心位置处探测器阵列数据的比值，得到校正因子，将所得校正因子替换到原始校正文件中。
上述步骤中，所述金标准数据使用前进行如下处理：
步骤C21，将所述金标准数据与所述中心位置处探测器阵列数据按照射野中心轴进行探头数据归一；
步骤C22，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据。
实施例2 磁共振加速器
本实施例提供一种磁共振加速器，包括：放射源、加速器组件、探测器阵列和控制系统。所述控制系统中集成有探测器阵列校正模块。
探测器阵列校正模块包括：
控制器，被配置为执行实施例1所述的探测器阵列校正方法；
存储器，被配置为存储程序代码和中间数据；
数据接口，被配置为与外部设备连接，执行数据的输入和输出。
本实施例中，除了所述探测器阵列校正模块之外，所述磁共振加速器的其他硬件结构和软件可根据现有技术实现。
所述探测器阵列校正模块的硬件可通过现有技术实现，所述探测器阵列校正模块被配置为执行实施例1所述的探测器阵列校正方法。
下面通过实验对本发明的技术方案做进一步说明。
实验例1 校正准确性验证
一、实验方法
本实验例对实施例1中所述方法的磁共振兼容探测器阵列校正准确性进行验证。本实验例中，加速器束流轮廓金标准数据使用PTW 31021 Semiflex 3D指型电离室、UnidosTango静电计、BeamScan MR水箱测量系统（PTW，Frieberg，Germany），在医科达Unity 1.5TMR-Linac（Elekta，Stockholm，Sweden）上测量得到。使用的探测器阵列为IC PROFILER-MR电离室探测器阵列（Sun Nuclear，Melbourne，USA），该设备在Y轴方向上均匀排列65个探测器，在X轴和正负对角线上分别有63个探测器。
实验组包括：实施例1中所述基于中心位置的单野校正方法、基于偏中心位置的多野分段校正方法，以及现有技术基于宽野的校正方法。
实验组具体说明如下：
1、基于中心位置的单野校正方法
加速器最大射野为57 cm×22 cm，探测器阵列尺寸为32 cm×32 cm。使用ICPROFILER-MR采集深度为4.9 cm、射野大小为40 cm×22 cm的束流轮廓数据，得到中心位置处探测器阵列数据，相同条件下水箱测量得到的金标准数据。其他操作按照实施例1的记载执行。本实施例中探测器阵列Y方向范围超出了最大射野大小，执行中心位置的单野校正方法时，忽略超出最大射野范围的探测器。
2、基于偏中心位置的多野分段校正方法
使用IC PROFILER-MR采集深度为4.9 厘米、射野大小为40 cm×22 cm的束流轮廓数据，相同条件下水箱测量得到的金标准数据。探测器阵列的移动方式为：将IC PROFILER-MR在X轴和Y轴方向上各偏移±8 cm，在正对角线轴和负对角线轴方向上各偏移±8√2 cm，得到偏中心位置处探测器阵列数据。其他操作按照实施例1的记载执行。
3、基于宽野的校正方法
具体方法参考如下文献执行：Simon T A, Simon W E, Kahler D, et al. Widefield array calibration dependence on the stability of measureddosedistributions[J]. Medical Physics, 2010, 37(7): 3501-3509。
二、实验结果
基于中心位置的单野校正方法与基于宽野的校正方法的校正结果比较如图1（a）至图4（b）所示，虚线为金标准束流轮廓数据，实线为探测器阵列测量并校正后得到的束流轮廓数据；左图为使用基于宽野的校正方法得到的结果，右图为使用基于中心位置的单野校正方法得到的结果。可见，现有技术基于宽野校正得到的结果与金标准数据差距较大，且部分数据区域存在异常点、数据不平滑，影响临床质控判断。而使用本发明中提到的基于中心位置的单野校正方法校正后的数据与金标准数据非常一致。如表1所示，基于中心射野校正方法得到的数据与金标准的平均伽马通过率（1%/1mm的gamma分析标准）为99.97%，基于宽野的校正方法的平均伽马通过率为80.62%，两者间存在统计学差异。
基于偏中心位置的多野分段校正方法与基于宽野的校正方法的校正结果比较如图5（a）至图8（b）所示，虚线为金标准束流轮廓数据，实线为探测器阵列测量并校正后得到的束流轮廓数据；左图为使用基于宽野的校正方法得到的结果，右图为使用基于偏中心位置的多野分段校正方法得到的结果。可见，现有技术基于宽野校正得到的结果与金标准数据差距较大，且部分数据区域存在异常点、数据不平滑，影响临床质控判断；而使用本发明中提供的基于偏中心位置的多野分段校正方法校正后的数据与金标准数据更加一致。如表1所示，基于偏中心位置的多野分段校正方法得到的数据与金标准的平均伽马通过率（1%/1mm的gamma分析标准）为93.13%，与基于宽野的校正方法的伽马通过率（平均80.62%）存在统计学差异。
表1 采用不同校正方法在3种不同深度下（4.9cm, 2.9cm, 0.9cm）测得的束流轮廓数据与金标准数据进行gamma分析得到的gamma通过率（gamma passing rate, GPR）的比较
可见，本发明的磁共振兼容探测器阵列校正方法校正结果的准确性相比于现有技术的方法具有显著提升。
通过上述实施例和实验例可以看到，本发明提供了一种更加简单和准确的磁共振兼容探测器阵列校正方法，并进一步提供了实现该方法的模块及集成有该方法的磁共振加速器。本发明能够更好地改善磁共振兼容探测器阵列在磁共振加速器中的应用效果，精确测量加速器束流轮廓数据，从而保证QA的精确性，提高工作效率，具有很好的应用前景。
Claims (10)
1.一种探测器阵列校正方法，其特征在于，加速器最大射野大小无法覆盖探测器阵列时，进行探测器阵列校正，包括如下步骤：
步骤1，将探测器阵列中心与射野中心对齐，将探测器阵列设置在与金标准数据相同深度处，收集加速器最大射野的束流轮廓数据，得到中心位置处探测器阵列数据；
步骤2，计算金标准数据与所述中心位置处探测器阵列数据的比值，得到中心位置处探测器阵列的校正因子Ic(r)，其中，r是指探测器阵列中的探测器；
步骤3，将探测器阵列维持在与金标准数据相同深度处，以所述中心位置为起点，设计探测器阵列的至少一条移动路线，所述探测器阵列在一条移动路线上移动n次，n为大于等于1的整数，不同移动路线上n的取值相同或不同，经过所有移动后探测器阵列的所有探测器至少一次位于加速器最大射野范围内；
在同一移动路线中，对于第j次移动，保证至少有一个共同的探测器rj，所述探测器rj满足在第j次移动前后，均位于加速器最大射野范围内；收集每次移动后加速器最大射野的束流轮廓数据，得到n组偏中心位置处探测器阵列数据；其中，j的取值为1至n的整数；
步骤4，在同一移动路线中，计算金标准数据与所述偏中心位置处探测器阵列数据的比值，得到偏中心位置处探测器阵列的校正因子ij(r)；
步骤5，在同一移动路线中，选择步骤3中所述共同的探测器rj作为参考点，按照如下公式计算修正因子uj：
当j=1，uj=Ic(rj)/ij(rj)；
当j大于1，uj= ij-1(rj)uj-1/ij(rj)；
步骤6，按照如下公式对偏中心位置处探测器阵列的校正因子进行修正：
Ij(r)=ij(r)×uj(r)，其中，Ij(r)即为修正后的偏中心位置校正因子；
步骤7，将Ic(r)和所有移动路线中得到的Ij(r)进行组合，得到最终的探测器阵列校正因子，将所述最终的探测器阵列校正因子替换到原始校正文件中。
2.按照权利要求1所述的探测器阵列校正方法，其特征在于：步骤2和步骤4中，所述金标准数据通过水箱测量法采集。
3.按照权利要求2所述的探测器阵列校正方法，其特征在于：步骤2中，所述金标准数据使用前进行如下处理：
步骤a1，将所述金标准数据与所述中心位置处探测器阵列数据按照射野中心轴进行探头数据归一；
步骤a2，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据；
步骤4中，所述金标准数据使用前进行如下处理：
步骤b1，根据步骤3中的移动距离调整所述金标准数据的坐标；
步骤b2，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据；
步骤b3，将所述金标准数据与所述偏中心位置处探测器阵列数据按照射野中心轴进行探头数据归一。
4.按照权利要求1所述的探测器阵列校正方法，其特征在于：步骤3中，保证每一次移动后探测器阵列的中心不在射野区域内。
5.按照权利要求1所述的探测器阵列校正方法，其特征在于：步骤7中，当同一探测器在不同移动路线中计算得到多个Ij(r)时，则最终取值按照如下方式中的一种进行选择:
1)选择初次获得的值；
2)任意选取一个取值；
3)计算平均值。
6.一种探测器阵列校正方法，其特征在于，进行探测器阵列校正，包括如下步骤：
步骤A，获取加速器束流轮廓金标准数据；
步骤B，判断加速器最大射野大小是否可覆盖探测器阵列；
若无法覆盖则按照权利要求1-5任一项所述的探测器阵列校正方法执行，进行基于偏中心位置的多野分段校正方法；
若可以覆盖则进行基于中心位置的单野校正方法，所述基于中心位置的单野校正方法按照如下步骤执行：
步骤C 1，将探测器阵列中心与射野中心对齐，将探测器阵列设置在与金标准数据相同深度处，收集加速器最大射野的束流轮廓数据，得到中心位置处探测器阵列数据；
步骤C2，计算金标准数据与所述中心位置处探测器阵列数据的比值，得到校正因子，将所得校正因子替换到原始校正文件中。
7.按照权利要求6所述的探测器阵列校正方法，其特征在于：步骤C2中，所述金标准数据使用前进行如下处理：
步骤C21，将所述金标准数据与所述中心位置处探测器阵列数据按照射野中心轴进行探头数据归一；
步骤C22，对所述金标准数据进行插值，获取探测器阵列中每个探测器位置坐标处对应的数据。
8.一种探测器阵列校正模块，其特征在于，包括：
控制器，被配置为执行权利要求1-7任一项所述的探测器阵列校正方法；
存储器，被配置为存储程序代码和中间数据；
数据接口，被配置为与外部设备连接，执行数据的输入和输出。
9.一种磁共振加速器，其特征在于，包括：放射源、加速器组件、探测器阵列和控制系统；所述控制系统中集成有权利要求8所述的探测器阵列校正模块。
10.一种计算机可读存储介质，其特征在于，其上存储有：用于实现权利要求1-7任一项所述的探测器阵列校正方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410885122.7A CN118416412B (zh) | 2024-07-03 | 2024-07-03 | 探测器阵列校正方法、模块、磁共振加速器及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410885122.7A CN118416412B (zh) | 2024-07-03 | 2024-07-03 | 探测器阵列校正方法、模块、磁共振加速器及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118416412A CN118416412A (zh) | 2024-08-02 |
| CN118416412B true CN118416412B (zh) | 2024-08-30 |
Family
ID=92323391
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410885122.7A Active CN118416412B (zh) | 2024-07-03 | 2024-07-03 | 探测器阵列校正方法、模块、磁共振加速器及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118416412B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103562740A (zh) * | 2011-05-31 | 2014-02-05 | 皇家飞利浦有限公司 | Mri放射治疗装置的静态磁场的校正 |
| CN105854191A (zh) * | 2016-04-26 | 2016-08-17 | 中国科学院合肥物质科学研究院 | 一种放射治疗中三维剂量验证系统及验证方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6807446B2 (en) * | 2002-09-03 | 2004-10-19 | Celsion Corporation | Monopole phased array thermotherapy applicator for deep tumor therapy |
| CA2572019A1 (en) * | 2004-06-24 | 2006-01-05 | Calypso Medical Technologies, Inc. | Systems and methods for treating a lung of a patient using guided radiation therapy or surgery |
| WO2007035775A2 (en) * | 2005-09-19 | 2007-03-29 | Feng Ma | Imaging system and method utilizing primary radiation |
| CN102473471A (zh) * | 2009-07-23 | 2012-05-23 | 太阳核公司 | 用于测量来自辐射源的辐射的多轴扫描系统和方法 |
| US10022564B2 (en) * | 2016-02-05 | 2018-07-17 | Varian Medical Systems International Ag | Systems, methods, and devices for radiation beam alignment and radiation beam measurements using electronic portal imaging devices |
| CN112635024A (zh) * | 2021-03-10 | 2021-04-09 | 四川大学 | 一种放射治疗自动计划设计系统及其构建方法 |
| CN118105638A (zh) * | 2024-01-09 | 2024-05-31 | 江苏瑞尔医疗科技有限公司 | 一种全球面体探测器阵列-模体装置及其使用方法 |
-
2024
- 2024-07-03 CN CN202410885122.7A patent/CN118416412B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103562740A (zh) * | 2011-05-31 | 2014-02-05 | 皇家飞利浦有限公司 | Mri放射治疗装置的静态磁场的校正 |
| CN105854191A (zh) * | 2016-04-26 | 2016-08-17 | 中国科学院合肥物质科学研究院 | 一种放射治疗中三维剂量验证系统及验证方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118416412A (zh) | 2024-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11691031B2 (en) | Systems, methods, and devices for radiation beam asymmetry measurements using electronic portal imaging devices | |
| US20240042241A1 (en) | Calibration of a stereoscopic camera system for use with a radio therapy treatment apparatus | |
| US20210138272A1 (en) | Systems, methods, and devices for radiation beam alignment and radiation beam measurements using electronic portal imaging devices | |
| US11058393B2 (en) | Imaging-based self-adjusting radiation therapy systems, devices, and methods | |
| Clivio et al. | Evaluation of the machine performance check application for TrueBeam Linac | |
| US20190175949A1 (en) | Time-resolved pre-treatment portal dosimetry systems, devices, and methods | |
| Low et al. | Evaluation of polymer gels and MRI as a 3‐D dosimeter for intensity‐modulated radiation therapy | |
| AU767922B2 (en) | Roentgen device for localization and radiation therapy of cancellation cancers | |
| Kanazawa et al. | Application of an RI-beam for cancer therapy: In-vivo verification of the ion-beam range by means of positron imaging | |
| EP1813308B1 (en) | Relative and absolute calibration for dosimetric devices | |
| Treuer et al. | On isocentre adjustment and quality control in linear accelerator basedradiosurgery with circular collimators and room lasers | |
| Quino et al. | Patient specific pre-treatment QA verification using an EPID approach | |
| Winkler et al. | Introducing a system for automated control of rotation axes, collimator and laser adjustment for a medical linear accelerator | |
| Dong et al. | Verification of radiosurgery target point alignment with an electronic portal imaging device (EPID) | |
| Dunscombe et al. | The investigation and rectification of field placement errors in the delivery of complex head and neck fields | |
| CN118416412B (zh) | 探测器阵列校正方法、模块、磁共振加速器及存储介质 | |
| US12109438B2 (en) | Positioning method realized by computer, and radiotherapy system | |
| JP2018000379A (ja) | 線量計算装置、線量計算方法、及び線量計算プログラム | |
| WO2013119887A1 (en) | Method and statistical validation technique for detecting differences between radiation therapy images with application to the detection and control of radiation therapy treatment delivery errors | |
| Roed et al. | The potential of polymer gel dosimeters for 3D MR-IGRT quality assurance | |
| Yang et al. | Using the volumetric effect of a finite‐sized detector for routine quality assurance of multileaf collimator leaf positioning | |
| JP2021069882A (ja) | 粒子線治療システム、粒子線照射判定装置、および、粒子線照射判定方法 | |
| EP3338860A1 (en) | Registration of particle beam radiography data | |
| Schippers et al. | Applications in radiation therapy of a scintillating screen viewed by a CCD camera | |
| KR101435497B1 (ko) | 2차원 전리함 선량계를 이용한 토모테라피용 환자별 정도 관리를 위한 방법 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |