CN116832345B - 电子设备、放射治疗系统及相关装置 - Google Patents
电子设备、放射治疗系统及相关装置 Download PDFInfo
- Publication number
- CN116832345B CN116832345B CN202310809495.1A CN202310809495A CN116832345B CN 116832345 B CN116832345 B CN 116832345B CN 202310809495 A CN202310809495 A CN 202310809495A CN 116832345 B CN116832345 B CN 116832345B
- Authority
- CN
- China
- Prior art keywords
- dose
- particle
- radiation
- treatment plan
- patient
- 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
- 238000001959 radiotherapy Methods 0.000 title claims abstract description 63
- 239000002245 particle Substances 0.000 claims abstract description 411
- 238000011282 treatment Methods 0.000 claims abstract description 287
- 238000009826 distribution Methods 0.000 claims abstract description 199
- 238000000275 quality assurance Methods 0.000 claims abstract description 157
- 230000005855 radiation Effects 0.000 claims abstract description 151
- 238000005259 measurement Methods 0.000 claims abstract description 148
- 230000005466 cherenkov radiation Effects 0.000 claims abstract description 100
- 230000003287 optical effect Effects 0.000 claims abstract description 88
- 238000000034 method Methods 0.000 claims abstract description 87
- 238000004364 calculation method Methods 0.000 claims abstract description 62
- 238000004590 computer program Methods 0.000 claims abstract description 35
- 238000003860 storage Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims description 48
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 30
- 238000006243 chemical reaction Methods 0.000 claims description 29
- 230000006872 improvement Effects 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000012797 qualification Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 12
- 210000001519 tissue Anatomy 0.000 description 35
- 230000000694 effects Effects 0.000 description 28
- 230000008901 benefit Effects 0.000 description 14
- 206010028980 Neoplasm Diseases 0.000 description 13
- 238000004422 calculation algorithm Methods 0.000 description 13
- 210000000056 organ Anatomy 0.000 description 11
- 238000002661 proton therapy Methods 0.000 description 11
- 230000008021 deposition Effects 0.000 description 10
- 238000000151 deposition Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 239000000463 material Substances 0.000 description 7
- 239000000126 substance Substances 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 238000013135 deep learning Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000002560 therapeutic procedure Methods 0.000 description 6
- 239000012634 fragment Substances 0.000 description 5
- 230000003993 interaction Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000004020 luminiscence type Methods 0.000 description 4
- 239000011164 primary particle Substances 0.000 description 4
- 230000001225 therapeutic effect Effects 0.000 description 4
- 238000004435 EPR spectroscopy Methods 0.000 description 3
- 238000013136 deep learning model Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 150000002500 ions Chemical class 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000010399 physical interaction Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000001950 radioprotection Effects 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 208000002454 Nasopharyngeal Carcinoma Diseases 0.000 description 2
- 206010061306 Nasopharyngeal cancer Diseases 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 208000012191 childhood neoplasm Diseases 0.000 description 2
- 238000002648 combination therapy Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 210000003414 extremity Anatomy 0.000 description 2
- 210000003141 lower extremity Anatomy 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 201000001441 melanoma Diseases 0.000 description 2
- 201000011216 nasopharynx carcinoma Diseases 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 231100000331 toxic Toxicity 0.000 description 2
- 230000002588 toxic effect Effects 0.000 description 2
- 210000001364 upper extremity Anatomy 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 208000006678 Abdominal Neoplasms Diseases 0.000 description 1
- 208000018084 Bone neoplasm Diseases 0.000 description 1
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 206010068771 Soft tissue neoplasm Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 210000003423 ankle Anatomy 0.000 description 1
- 210000000617 arm Anatomy 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 244000309466 calf Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 235000019800 disodium phosphate Nutrition 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000008713 feedback mechanism Effects 0.000 description 1
- 210000002683 foot Anatomy 0.000 description 1
- 210000000245 forearm Anatomy 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000011866 long-term treatment Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 208000025440 neoplasm of neck Diseases 0.000 description 1
- 208000025426 neoplasm of thorax Diseases 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000004224 protection Effects 0.000 description 1
- 238000004549 pulsed laser deposition Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 208000011581 secondary neoplasm Diseases 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 201000003957 thoracic cancer Diseases 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
- 210000000707 wrist Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N2005/1085—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
- A61N2005/1087—Ions; Protons
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- Algebra (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Animal Behavior & Ethology (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
本申请提供了用于放射治疗的质量保证方法、电子设备、计算机可读存储介质、计算机程序产品及放射治疗系统，所述方法包括：获取切伦科夫辐射‑粒子剂量对应关系；使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射；计算得到患者光学测量信息对应的所述患者表面的粒子剂量；计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果。本申请通过对切伦科夫辐射的检测，获取剂量分布的实际测量数据并与理论计算数据进行对比，从而判断治疗计划的质量保证结果，提高放射治疗的准确性和安全性。
Description
技术领域
本申请涉及放射治疗的技术领域，尤其涉及用于放射治疗的质量保证方法、电子设备、计算机可读存储介质、计算机程序产品及放射治疗系统。
背景技术
随着技术日新月异的发展，放射治疗及放射治疗参与的联合治疗在肿瘤治疗中的作用日益彰显。质子治疗在鼻咽癌、黑色素瘤、儿童肿瘤治疗方面有明显的治疗优势。由于质子与物质相互作用时，剂量沉积表现出的Bragg峰（布拉格峰）曲线分布，使得质子治疗在取得肿瘤治疗显著疗效的同时，能够更好的保护危机器官，减少放射治疗的毒副作用。
质量保证是质子精准治疗的重要手段，目前的质量保证方法采用水模体进行模拟，无法保证验证结果的准确性。
基于此，本申请提供了用于放射治疗的质量保证方法、电子设备、计算机可读存储介质、计算机程序产品及放射治疗系统，以改进相关技术。
发明内容
本申请的目的在于提供用于放射治疗的质量保证方法、电子设备、计算机可读存储介质、计算机程序产品及放射治疗系统，通过光学测量设备对切伦科夫辐射的检测，准确获取患者表面的粒子剂量，用于获取剂量分布的实际测量数据并与理论计算数据进行对比，从而判断治疗计划的质量保证结果，提高放射治疗的准确性和安全性，为辐射治疗带来更广泛的应用前景。
本申请的目的采用以下技术方案实现：
本申请提供了一种用于放射治疗的质量保证方法，所述方法包括：
获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系；
在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；
根据所述患者光学测量信息和所述切伦科夫辐射-粒子剂量对应关系，计算得到所述患者光学测量信息对应的所述患者表面的粒子剂量；
根据所述患者表面的粒子剂量，计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格。
在一些可能的实现方式中，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在每种粒子辐射条件下执行以下处理：使用所述光学测量设备测量模体受到所述粒子辐射而发射的切伦科夫辐射，以得到模体光学测量信息；使用剂量测量设备获取所述模体表面的粒子剂量；其中，所述模体为人体模拟组织模体；
根据多种粒子辐射条件下的模体光学测量信息和模体表面的粒子剂量，建立所述模体光学测量信息和所述模体表面的粒子剂量之间的对应关系，作为所述切伦科夫辐射-粒子剂量对应关系。
在一些可能的实现方式中，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在多种粒子辐射条件下，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到人体表面的粒子剂量和次级电子剂量之间的粒子剂量-电子剂量对应关系；
根据所述粒子剂量-电子剂量对应关系和预设的电子剂量-切伦科夫辐射对应关系，获取所述切伦科夫辐射-粒子剂量对应关系；
其中，所述电子剂量-切伦科夫辐射对应关系用于指示切伦科夫辐射强度与电子剂量之间的正比关系。
在一些可能的实现方式中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
根据所述治疗计划对应的粒子辐射条件，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息；
根据所述粒子在人体内部的剂量分布信息，生成所述粒子辐射的参考深度剂量分布曲线；
根据所述患者表面的粒子剂量和所述粒子辐射的参考深度剂量分布曲线，计算出所述粒子在所述患者内部的剂量分布信息。
在一些可能的实现方式中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
将所述治疗计划对应的粒子辐射条件和所述患者表面的粒子剂量输入基于深度学习的剂量分布模型，以得到所述粒子在所述患者内部的剂量分布信息。
在一些可能的实现方式中，所述根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，包括：
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到所述实际测量数据与所述理论计算数据对应的伽马通过率；
当所述伽马通过率满足预设的通过率条件时，将所述治疗计划的质量保证结果设置为合格；
当所述伽马通过率不满足所述通过率条件时，将所述治疗计划的质量保证结果设置为不合格。
在一些可能的实现方式中，所述方法还包括：
当所述质量保证结果为不合格时，根据所述治疗计划和所述伽马通过率，获取对应的质量改进方案；
生成质量保证报告，在所述质量保证报告中包括所述治疗计划的质量保证结果、所述伽马通过率和所述质量改进方案。
在一些可能的实现方式中，每种对应关系采用以下任意一种方式表示：转换系数表；转换模型；转换矩阵；转换公式。
在一些可能的实现方式中，所述粒子为质子；或者，
所述粒子为质子治疗对应的任意一种粒子。
第二方面，本申请提供了一种电子设备，所述电子设备包括存储器和至少一个处理器，所述存储器存储有计算机程序，所述至少一个处理器被配置成执行所述计算机程序时实现以下步骤：
获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系；
在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；
根据所述患者光学测量信息和所述切伦科夫辐射-粒子剂量对应关系，计算得到所述患者光学测量信息对应的所述患者表面的粒子剂量；
根据所述患者表面的粒子剂量，计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格。
第三方面，本申请提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被至少一个处理器执行时实现上述任一项方法的步骤或者实现上述任一项电子设备的功能。
第四方面，本申请提供了一种计算机程序产品，所述计算机程序产品包括计算机程序，所述计算机程序被至少一个处理器执行时实现上述任一项方法的步骤或者实现上述任一项电子设备的功能。
第五方面，本申请提供了一种放射治疗系统，包括：
质量保证装置，用于采用上述任一项用于放射治疗的质量保证方法，获取治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格；
确定人体剂量的装置，用于在所述质量保证结果为合格后，根据所述治疗计划确定放射治疗中使用的剂量。
附图说明
下面结合说明书附图和具体实施方式进一步说明本申请。
图1是本申请实施例提供的一种用于放射治疗的质量保证方法的流程示意图。
图2是本申请实施例提供的一种放射治疗系统的结构示意图。
图3是本申请实施例提供的一种电子设备的结构框图。
图4是本申请实施例提供的一种计算机程序产品的结构示意图。
具体实施方式
下面将结合本申请的说明书附图以及具体实施方式，对本申请中的技术方案进行描述，需要说明的是，在不相冲突的前提下，以下描述的各实施方式之间或各技术特征之间可以任意组合形成新的实施方式。
本申请实施例中，“示例性的”或者“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性的”或者“例如”的任何实施方式或设计方案不应被解释为比其他实施方式或设计方案更优选或更具优势。确切而言，使用“示例性的”或者“例如”等词旨在以具体方式呈现相关概念。
本申请实施例中出现的第一、第二等描述，仅作示意与区分描述对象之用，没有次序之分，也不表示本申请实施例中对数量的特别限定，不能构成对本申请实施例的任何限制。
下面对本申请实施例的技术领域和相关术语进行简单说明。
质子治疗（Proton Therapy）是一种放射治疗技术，它利用高能质子束来精确治疗肿瘤。相比传统的X射线放疗，质子治疗可以更好地控制放射剂量的传递，减少对正常组织的伤害，提高治疗效果。
质子治疗的原理是利用质子的物理特性，即在进入人体后，质子束会在一定深度达到最大的剂量（布拉格峰），然后急剧减小直至停止。这种特性使得质子束能够在肿瘤内释放最大剂量，而减少在肿瘤后方正常组织中的剂量沉积，从而降低了治疗引起的副作用。质子治疗适用于多种类型的肿瘤，包括儿童肿瘤、颅脑肿瘤、头颈部肿瘤、胸部肿瘤、腹部肿瘤、骨和软组织肿瘤等。质子治疗尤其适用于在关键器官周围或对放射敏感的肿瘤。质子治疗相比传统放疗可以更好地保护正常组织和器官，减少治疗引起的副作用。特别是对于儿童患者，质子治疗可以减少长期的治疗后遗症，降低二次肿瘤的风险。一些研究表明，质子治疗在某些情况下可以提供与传统放疗相当的治疗效果，同时减少了不良反应。
蒙特卡罗（又译为蒙特卡洛，Monte Carlo, MC）方法被誉为粒子剂量计算的金标准。相较于传统的解析算法，其具有计算精度高等优点。临床放疗中，MC计算粒子在人体感兴趣区域（Regions of interest, ROI）内的剂量沉积时，会根据治疗部位及CT值等信息，将CT图像上的每个体素替换为对应的真实人体器官材料，从而实现从虚拟图像到真实人体的映射，对粒子在人体组织中进行输运和模拟。由于具有了真实的器官材料属性，MC对粒子的输运和模拟，粒子的相互作用都认为是在人体组织中发生的，产生的剂量沉积属于器官（人体组织）剂量。
与MC对应的另一种剂量计算方法是解析算法。由于算法属性的差异，解析算法进行剂量计算时，直接使用不同密度的水对CT图像上的体素进行替换，密度值由对应位置上的CT值决定。因此其计算的剂量认为是水中的剂量。解析算法计算效率高，对计算机硬件要求低，目前是剂量计算中使用最广泛、普遍的算法。
随着技术日新月异的发展，放射治疗及放射治疗参与的联合治疗在肿瘤治疗中的作用日益彰显。质子治疗在鼻咽癌、黑色素瘤、儿童肿瘤治疗方面有明显的治疗优势。由于质子与物质相互作用时，剂量沉积表现出的Bragg峰（布拉格峰）曲线分布，使得质子治疗在取得肿瘤治疗显著疗效的同时，能够更好的保护危机器官，减少放射治疗的毒副作用。随着技术的快速发展，将会有更多患者接受质子治疗这种精准的治疗方式。
质量保证是质子精准治疗的重要手段。相关质量保证的实现方式是将治疗计划复制在水模体上，利用外部的剂量测量设备测量到的水模体剂量分布与治疗计划系统计算的水模体剂量分布进行比较，从而验证治疗计划剂量配送的准确性。但水模体中的剂量分布和人体内的剂量分布存在差异。目前缺乏对人体内剂量分布进行有效准确测量的有效手段。
基于此，本申请提供了用于放射治疗的质量保证方法、电子设备、计算机可读存储介质、计算机程序产品及放射治疗系统，以改进相关技术。
本申请实施例提供的方案涉及放射治疗、模拟仿真等技术，具体通过如下实施例进行说明。需要说明的是，以下实施例的描述顺序不作为对实施例优选顺序的限定。
（用于放射治疗的质量保证方法）
参见图1，图1是本申请实施例提供的一种用于放射治疗的质量保证方法的流程示意图。
本申请实施例提供了一种用于放射治疗的质量保证方法，所述方法包括：
步骤S101：获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系；
步骤S102：在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；
步骤S103：根据所述患者光学测量信息和所述切伦科夫辐射-粒子剂量对应关系，计算得到所述患者光学测量信息对应的所述患者表面的粒子剂量；
步骤S104：根据所述患者表面的粒子剂量，计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；
步骤S105：根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格。
本申请实施例对用于放射治疗的粒子不作限定，其例如可以采用质子、重离子、电子、X射线、γ射线等。这些粒子例如具有不同的能量和质量，可用于针对肿瘤的定向辐射。
切伦科夫辐射(Cherenkov Radiation，CR)是指带电粒子在介质中的运动速度超过光速时产生的一种发光现象，其本质是带电粒子穿过介质时引起的介质中原子或分子的极化。例如，人体受到质子辐射后，在治疗过程中会产生可见光的辐射现象，该光场可以被探测到。
本申请实施例中，粒子本身为带电粒子（如质子等），或者，粒子辐射人体时可以产生带电粒子（如次级电子等）。
用于探测切伦科夫辐射光场的光学测量设备例如可以采用光检测器，用于检测来自患者和/或透光电介质的切伦科夫辐射。透光电介质可以是气体、液体和固体。人体和透光电介质的厚度或体积不同，对应的切伦科夫辐射也随之变化。光检测器的数量可以是一个或多个。多个光检测器可以间隔开，形成一个光检测器阵列。在患者可能遮挡单个光检测器的视野（例如粒子辐射源和透光电介质位于患者上方，该光检测器位于患者下方）的情况下，其他光检测器仍然能够检测到切伦科夫辐射。透光电介质可以定位成使得自身产生的任何切伦科夫辐射不会被患者遮挡，例如可以将透光电介质和光检测器布置成比患者定位更靠近粒子辐射源。光检测器可以是光学照相机、闭路电视、摄像机等，在一些实施例中，它们也用于拍摄、记录、监测患者或者其他对象。其他对象例如可以是人体模拟组织模体、水模体、治疗床、治疗仪器等。
粒子辐射源可以安装在可旋转的机架上，机架可围绕患者或者治疗床旋转。在一些实施例中，光检测器可以安装在机架上，并且可以与粒子辐射源一起旋转。在放射治疗或标定过程期间，粒子辐射源围绕患者或者治疗床旋转。
人体表面的粒子剂量指治疗过程中人体表面的粒子剂量沉积，粒子剂量的单位例如是戈瑞（Gy）、毫戈瑞（mGy）等。
切伦科夫辐射-粒子剂量对应关系是指经过一系列测量或模拟分析，确定的切伦科夫辐射强度与患者表面粒子剂量之间的关系，对应关系例如可以采用模型、函数、表格、矩阵、公式等来表示，本申请不对此设限。
患者光学测量信息是指使用光学测量设备对患者受到质子辐射后发射的切伦科夫辐射进行测量，记录下的光学测量信息。作为示例，光学测量信息是切伦科夫辐射强度（的数值）。
患者表面的粒子剂量是根据患者光学测量信息和切伦科夫辐射-粒子剂量对应关系得出的计算结果。作为一个示例，计算结果为50 Gy。
粒子在患者内部的剂量分布信息是基于患者表面的粒子剂量，计算得到的粒子在患者内部的剂量分布，描述了粒子剂量随深度和位置的变化情况。剂量分布可以采用剂量分布曲线、三维剂量分布图、剂量分布表、等剂量线、剂量体积直方图等来表示。该粒子在患者内部的剂量分布信息作为治疗计划的剂量分布的实际测量数据。
剂量分布可以通过不同的表现形式来呈现。剂量分布曲线以剂量值为纵轴，位置（例如深度或距离）为横轴，展示了在特定位置上的剂量分布情况。剂量分布于曲线可以帮助理解剂量随深度或距离的变化。
等剂量线是连接相同剂量值的线条，通常以不同颜色或线型表示。等剂量线的密集程度显示了剂量在空间中的分布情况，可以直观地展示剂量的高低和分布范围。
三维剂量分布图使用三维空间来展示剂量的分布情况，通常以彩色或灰度表示。它提供了更全面的剂量信息，可以清晰地展示剂量在不同位置的强度和空间分布。
剂量体积直方图以剂量值为横轴，体积或体积百分比为纵轴，展示了不同剂量值下的体积或体积百分比，可以提供剂量分布的统计信息，比如特定体积接收到的剂量。
治疗计划的剂量分布的理论计算数据是指放射治疗系统根据治疗计划中设定的参数、计算模型和计算方法，通过数值模拟或计算算法得到的理论剂量分布数据。
根据实际测量数据与理论计算数据的比较，判断治疗计划的剂量分布是否符合预期要求，得到质量保证结果，如治疗计划合格或不合格。质量保证结果例如可以采用中文、字母、数字、符号中的一种或多种来表示。如使用“合格”、“Y”、“1”、“√”、“A1”等来表示合格，使用“不合格”、“N”、“0”、“×”、“B2”等来表示不合格。
由此，首先，通过实验或模拟研究，获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的对应关系。在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，从中获取患者光学测量信息。利用患者光学测量信息和切伦科夫辐射-粒子剂量对应关系，进行计算，得到患者表面的粒子剂量。根据患者表面的粒子剂量，通过相关算法或模型，计算出粒子在患者内部的剂量分布信息，这提供了治疗计划对应的剂量分布的实际测量数据。将剂量分布的实际测量数据与治疗计划对应的剂量分布的理论计算数据进行比较，以评估治疗计划的质量。根据实际测量数据与理论计算数据的对比，判断治疗计划的质量保证结果，该结果将指示质量保证是否合格或不合格。
这样做的好处是，基于实际测量数据，确保粒子治疗的准确性和精确性，通过比较实际测量数据与理论计算数据，可以评估和保证治疗计划的质量；通过获取患者表面的粒子剂量和粒子在患者内部的剂量分布信息，提供了更准确的剂量分布数据，从而提高治疗计划的精确性；通过实际测量数据与理论计算数据的比较，可以及时发现和纠正治疗计划中的潜在问题，从而提高治疗的安全性和有效性；该方法提供了一种基于实际测量数据的反馈机制，可以用于优化和改进治疗计划，从而进一步提高治疗的质量和效果。
具体而言，通过获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系，结合光学测量设备测量患者受到的切伦科夫辐射，可以计算得到患者表面的粒子剂量，这种基于切伦科夫辐射的测量方法能够提供实际的粒子剂量数据，从而实现对质子治疗质量的准确保证。通过患者表面粒子剂量计算得到治疗计划对应的剂量分布的实际测量数据，可以将实际测量数据与理论计算数据进行对比。这种对比能够评估治疗计划的准确性和精确度，为粒子精准治疗的质量保证提供可靠的依据。基于实际测量数据与理论计算数据的对比，可以获取治疗计划的质量保证结果。根据所得结果，可以判断治疗计划是合格还是不合格，能够确保质子治疗过程中的准确性和安全性，提高肿瘤治疗的效果。
综上所述，该方法提供了一种能够实现放射治疗质量保证的方法，通过光学测量设备对切伦科夫辐射的检测，该方法能够准确获取患者表面的粒子剂量，用于获取剂量分布的实际测量数据并与理论计算数据进行对比，从而判断治疗计划的质量保证结果，这将提高放射治疗的准确性和安全性，为肿瘤辐射治疗带来更广泛的应用前景。
在一些实施例中，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在每种粒子辐射条件下执行以下处理：使用所述光学测量设备测量模体受到所述粒子辐射而发射的切伦科夫辐射，以得到模体光学测量信息；使用剂量测量设备获取所述模体表面的粒子剂量；其中，所述模体为人体模拟组织模体；
根据多种粒子辐射条件下的模体光学测量信息和模体表面的粒子剂量，建立所述模体光学测量信息和所述模体表面的粒子剂量之间的对应关系，作为所述切伦科夫辐射-粒子剂量对应关系。
获取模体（人体模拟组织模体）的切伦科夫辐射测量结果（模体光学测量信息）与模体表面的粒子剂量之间的关系，作为步骤S101中的切伦科夫辐射-粒子剂量对应关系。通过建立对应关系，可以根据例如是切伦科夫辐射强度的模体光学测量信息推导出相应的粒子剂量（数值）。举例来说，在粒子辐射实验中，使用光学测量设备测量人体模拟组织模体受到粒子辐射后发射的切伦科夫辐射的强度，同时使用剂量测量设备获取模体表面的粒子剂量。通过分析多组测量数据，建立模体光学测量信息和模体表面粒子剂量之间的对应关系，例如，假设根据实验结果得出切伦科夫辐射强度每增加1个单位，模体表面的粒子剂量增加2 Gy，则该对应关系可以采用函数或者公式来表示。
模体光学测量信息指的是通过光学测量设备测量得到的模体（人体模拟组织模体）受到粒子辐射后发射的切伦科夫辐射的光场探测信息，例如是切伦科夫辐射强度。举例来说，使用光学测量设备对人体模拟组织模体受粒子辐射而发射的切伦科夫辐射进行测量，记录下切伦科夫辐射的强度数值。例如，测量结果显示，在某种粒子辐射条件下，模体的切伦科夫辐射强度为5个单位。
模体表面的粒子剂量指的是通过剂量测量设备获取到的模体（人体模拟组织模体）表面的粒子辐射剂量沉积。举例来说，使用剂量测量设备测量人体模拟组织模体表面的粒子辐射剂量，并以Gray（Gy）为单位表示。例如，测量结果显示，在某种粒子辐射条件下，模体表面的粒子剂量为60 Gy。
在一些实施例中，所述剂量测量设备包括电离室、闪烁体探测器、热释电探测器、光刺激发光剂量计和电子自旋共振剂量计中的一种或多种。
电离室是一种常用的剂量测量设备，通过测量粒子辐射通过电离室时产生的电荷量来确定粒子剂量。
闪烁体探测器利用闪烁效应来测量粒子辐射，当粒子通过闪烁体时，闪烁体会发出光信号，该信号与粒子剂量相关。
热释电探测器使用热释电材料来测量粒子辐射剂量，粒子辐射会使材料中的电子被激发，当材料加热时，被激发的电子会重新回到基态，释放出可测量的热量。
光刺激发光剂量计使用材料中的固有晶格缺陷来测量粒子辐射剂量，通过光激发这些缺陷并测量释放的发光强度来确定剂量。
电子自旋共振剂量计利用材料中的自由基或离子缺陷来测量粒子辐射剂量，通过测量自旋共振信号的光强度来确定剂量。
由此，在每种粒子辐射条件下进行以下两种处理：使用光学测量设备测量模体的切伦科夫辐射，使用剂量测量设备获取模体表面的粒子剂量。之后，建立模体光学测量信息和模体表面粒子剂量之间的对应关系，这个对应关系可以作为切伦科夫辐射-粒子剂量的参考关系（或者说标定关系）。
这样做的好处是，通过测量模体受到粒子辐射而发射的切伦科夫辐射和模体表面的粒子剂量，建立了切伦科夫辐射-粒子剂量的对应关系，这个对应关系可以用于后续的粒子治疗过程中的质量保证；通过光学测量设备测量模体的切伦科夫辐射，可以准确获取模体光学测量信息，为后续的对应关系建立提供准确的数据基础；通过建立切伦科夫辐射-粒子剂量对应关系，有助于实现粒子（辐射）治疗的质量保证方法，提高治疗计划的准确性和治疗效果，同时保护周围正常组织。综上所述，该方法通过测量模体的切伦科夫辐射和表面粒子剂量，建立了切伦科夫辐射-粒子剂量的对应关系，这个方法能够提供精确的模体光学测量信息，为粒子治疗的质量保证提供参考或者标定依据，并进一步完善治疗计划的准确性和治疗效果。
在一些实施例中，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在多种粒子辐射条件下，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到人体表面的粒子剂量和次级电子剂量之间的粒子剂量-电子剂量对应关系；
根据所述粒子剂量-电子剂量对应关系和预设的电子剂量-切伦科夫辐射对应关系，获取所述切伦科夫辐射-粒子剂量对应关系；
其中，所述电子剂量-切伦科夫辐射对应关系用于指示切伦科夫辐射强度与电子剂量之间的正比关系。
本申请实施例中，粒子辐射条件用于指示一个或多个粒子辐射参数，包括粒子类型、粒子能量、辐射方向、粒子辐射强度、辐射时间等。
粒子类型是指粒子的种类，例如质子、重离子、电子等。
粒子能量是指粒子的能量水平，例如以兆电子伏（MeV）表示。
辐射方向是指粒子辐射的入射方向和角度。
粒子辐射强度是指粒子辐射（照射）的强度，例如以剂量率（例如每分钟或每秒辐射剂量）或总剂量（例如总辐射剂量）表示。
辐射时间是指粒子辐射的持续时间，例如以秒、分钟、小时或天数表示。辐射时间可以是多个分段辐射时间之和。
上述粒子辐射参数的具体数值的选择因治疗计划和患者的个性化情况而有所不同。例如，在质子治疗中，常见的粒子能量范围可以从70 MeV到250 MeV，剂量率可以在每分钟1 Gy到10 Gy之间，辐射时间可以从几分钟到几周。
粒子剂量-电子剂量对应关系指的是人体表面的粒子剂量与次级电子剂量之间的关系。通过建立该对应关系，可以根据人体表面的粒子剂量值推导出相应的次级电子剂量值。通过使用蒙特卡罗方法模拟多种粒子辐射条件下粒子在人体中的输运过程，得到人体表面的粒子剂量和次级电子剂量之间的对应关系。例如，在特定粒子辐射条件下，当人体表面的粒子剂量为50 Gy时，相应的次级电子剂量为30 Gy。
电子剂量-切伦科夫辐射对应关系指的是电子剂量与切伦科夫辐射强度之间的关系。通过建立该对应关系，可以根据次级电子剂量值推导出相应的切伦科夫辐射强度值。例如，根据预设的电子剂量-切伦科夫辐射对应关系，当次级电子剂量为30 Gy时，相应的切伦科夫辐射强度为40个单位。
在本实施例中，切伦科夫辐射-粒子剂量对应关系的获取过程是通过使用蒙特卡罗方法模拟粒子在人体中的输运过程，得到人体表面的粒子剂量和次级电子剂量之间的关系，再根据电子剂量-切伦科夫辐射对应关系，获得切伦科夫辐射-粒子剂量对应关系，这样就可以根据切伦科夫辐射的光学测量信息推导出相应的粒子剂量值。
由此，使用蒙特卡罗方法，在多种粒子辐射条件下，模拟粒子在人体中的输运过程。根据模拟结果可以得到人体表面的粒子剂量和次级电子剂量之间的粒子剂量-电子剂量对应关系，这个对应关系用于衔接粒子剂量和切伦科夫辐射之间的关系。在预设的电子剂量-切伦科夫辐射对应关系的基础上，进一步获取切伦科夫辐射-粒子剂量对应关系，该对应关系可以用于后续的质量保证过程。
这样做的好处是，通过蒙特卡罗模拟粒子在人体中的输运过程，可以获得粒子剂量和次级电子剂量之间的对应关系，这种对应关系能够提供准确的粒子剂量测量和计算的基础；通过预设的电子剂量-切伦科夫辐射对应关系和粒子剂量-电子剂量对应关系，可以获取切伦科夫辐射-粒子剂量的对应关系；通过粒子剂量和切伦科夫辐射之间的对应关系，可以对治疗计划中的粒子剂量进行精确评估。综上所述，通过蒙特卡罗模拟，实现了粒子剂量-电子剂量-切伦科夫辐射对应关系的获取，这有助于提高治疗计划的准确性，优化质子治疗的剂量测量和质量保证过程。
在一些实施例中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
根据所述治疗计划对应的粒子辐射条件，获取所述粒子辐射的参考深度剂量分布曲线；
根据所述患者表面的粒子剂量和所述粒子辐射的参考深度剂量分布曲线，计算出所述粒子在所述患者内部的剂量分布信息。
其中，所述根据所述治疗计划对应的粒子辐射条件，获取所述粒子辐射的参考深度剂量分布曲线，包括：
根据所述治疗计划对应的粒子辐射条件，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息；根据所述粒子在人体内部的剂量分布信息，生成所述粒子辐射的参考深度剂量分布曲线；或者，
根据所述治疗计划对应的粒子辐射条件，使用剂量测量设备测量得到所述粒子辐射的参考深度剂量分布曲线。
在一些实施例中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
根据所述治疗计划对应的粒子辐射条件，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息；
根据所述粒子在人体内部的剂量分布信息，生成所述粒子辐射的参考深度剂量分布曲线；
根据所述患者表面的粒子剂量和所述粒子辐射的参考深度剂量分布曲线，计算出所述粒子在所述患者内部的剂量分布信息。
剂量分布信息是指粒子辐射在患者内部各部位（或者说各器官、各人体组织）的剂量分布情况。它描述了粒子辐射在患者组织中沉积的能量分布，用于评估治疗计划的效果和确定治疗区域。例如，在质子治疗中，使用蒙特卡罗方法模拟质子在人体组织中的运动和相互作用，以得到患者内部的剂量分布信息。在蒙特卡罗模拟中，所使用的人体例如是参考人。
在一些实施例中，根据所述治疗计划对应的粒子辐射条件以及所述患者对应的参考人信息，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息。所述参考人信息包括参考人的性别信息、部位信息和年龄信息中的一种或多种。
参考人例如可以是ICRP参考人或者特定国家、特定地区的参考人。作为一个示例，参考人可以是中国参考人。ICRP是国际放射防护委员会（International Commission onRadiological Protection）的缩写。ICRP参考人是国际放射防护委员会定义的一个标准人类模型，用于辐射防护计算中，作为一个参考对象。ICRP参考人最初定义为一名20-30岁、体重70千克、身高170厘米、生活在平均气温在10-20摄氏度的气候中的男性。但随着时间的推移，ICRP已经扩展了这一定义，包括不同年龄段和性别的参考人。目前，ICRP提供了一系列参考人的数据，包括男性和女性，每种性别分别包括新生儿、1岁、5岁、10岁、15岁和成年人。其中，性别信息包括男性、女性。人体由许多不同的部位组成，包括头部、颈部、躯干、四肢等。躯干包括胸部、腹部和背部。四肢包括上肢和下肢。上肢包括肩膀、上臂、前臂、手腕和手。下肢包括髋部、大腿、小腿、脚踝和脚。每种部位包括一种或多种器官，每种器官包括多种（人体）组织。
参考深度剂量分布曲线描述粒子辐射在人体内部不同深度处的剂量分布情况的参考信息。参考深度剂量分布曲线例如是蒙特卡罗模拟得到的剂量分布的IDD（IntegratedDepth Dose，积分深度剂量分布）曲线。在质子治疗中，根据患者体表的粒子剂量和模拟得到的参考深度剂量分布曲线，可以推断出质子在患者组织深度处的剂量分布情况。
作为示例，粒子在患者内部的剂量分布信息的获取过程是通过使用蒙特卡罗方法模拟粒子在人体中的输运过程，得到治疗计划对应的粒子在人体内部的剂量分布信息，进一步得到剂量分布的IDD曲线，IDD曲线是蒙特卡罗模拟结果，而非实际测量结果。然后，根据IDD曲线和患者表面的粒子剂量（由探测出的切伦科夫辐射对应得到），计算出粒子在患者内部的剂量分布信息。
由此，根据治疗计划对应的粒子辐射条件，使用蒙特卡罗方法进行模拟，模拟粒子在人体中的输运过程，得到治疗计划对应的粒子在人体内部的剂量分布信息；基于蒙特卡罗模拟结果，生成粒子辐射的参考深度剂量分布曲线，参考深度剂量分布曲线显示了粒子在人体内部不同位置的剂量分布情况；利用患者表面的粒子剂量和粒子辐射的参考深度剂量分布曲线进行计算，得到粒子在患者内部的剂量分布信息，这提供了治疗计划对应的剂量分布的实际测量数据。
这样做的好处是，通过蒙特卡罗模拟粒子的输运过程，可以准确模拟和计算粒子在患者内部的剂量分布，这提供了治疗计划对应的粒子剂量分布信息，为精确的剂量计划和质量保证提供了依据；通过生成粒子辐射的参考深度剂量分布曲线，可以对粒子在患者内部的剂量分布进行可视化和分析，这有助于评估治疗计划的准确性和剂量分布的一致性；通过计算粒子在患者内部的剂量分布信息，可以更准确地评估和调整治疗计划，以确保粒子的剂量分布符合治疗目标，提高治疗的准确性和安全性。
在一些实施例中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
将所述治疗计划对应的粒子辐射条件和所述患者表面的粒子剂量输入基于深度学习的剂量分布模型，以得到所述粒子在所述患者内部的剂量分布信息。
治疗计划是为了治疗患者的肿瘤或疾病而制定的放射治疗方案。治疗计划可以包括粒子辐射条件，例如粒子类型、粒子能量、辐射方向等。例如，在一种质子治疗计划中，使用质子粒子，质子能量为150 MeV，入射角度为45度，辐射时间为5分钟。
剂量分布模型是一种用于预测和模拟粒子在人体内部的剂量分布的数学模型，基于已知的治疗计划和人体的解剖结构，通过算法或统计方法来估计粒子在人体内部各处的剂量分布。举例来说，使用多个治疗计划和相应的人体模拟组织模体表面的粒子剂量、模体内部的剂量分布信息训练初始的深度学习模型，当满足预设的训练结束条件时停止训练。该训练好的剂量分布模型可以针对任意输入预测出剂量分布信息，例如以剂量分布曲线或者三维剂量分布图来表示，显示出粒子在患者不同部位的剂量分布情况。
在本实施例中，粒子在患者内部的剂量分布信息的获取过程是通过将治疗计划对应的粒子辐射条件和患者表面的粒子剂量输入基于深度学习的剂量分布模型，来预测得到粒子在患者内部的剂量分布信息。
由此，将治疗计划对应的粒子辐射条件和患者表面的粒子剂量作为输入，利用基于深度学习的剂量分布模型，对输入的数据进行处理和分析。剂量分布模型经过训练，可以预测和生成粒子在患者内部的剂量分布信息。通过剂量分布模型输出粒子在患者内部的剂量分布信息，提供了实际治疗计划对应的剂量分布数据。
这样做的好处是，基于深度学习的剂量分布模型能够通过输入治疗计划和患者表面的粒子剂量，快速而准确地获取粒子在患者内部的剂量分布信息，这提高了数据处理的效率和精确度；通过深度学习模型的训练和优化，可以提高剂量分布模型的准确性，例如可以根据每个患者的个性化信息，对剂量分布模型的模型参数进行调整，以得到适用于每个患者的剂量分布模型，这有助于更准确地预测和生成粒子在特定患者内部的剂量分布信息；通过获取粒子在患者内部的剂量分布信息，可以对治疗计划进行优化和调整，这有助于提高治疗计划的准确性和个体化，以更好地满足患者的治疗需求；准确获取粒子在患者内部的剂量分布信息有助于提高治疗的精确性和安全性，这可以帮助放射治疗专家更好地控制剂量分布，减少对健康组织的副作用，提高治疗的效果。
在一些实施例中，所述根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，包括：
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到所述实际测量数据与所述理论计算数据对应的伽马通过率；
当所述伽马通过率满足预设的通过率条件时，将所述治疗计划的质量保证结果设置为合格；
当所述伽马通过率不满足所述通过率条件时，将所述治疗计划的质量保证结果设置为不合格。
实际测量数据是指在治疗计划对应的放射治疗过程中，使用光学测量设备对切伦科夫辐射进行测量得到的光学测量信息所对应的粒子在患者体内的剂量分布的计算结果。由于切伦科夫辐射是实际测量得到的，因此其所对应的剂量分布的计算结果看作是实际测量数据。这些数据反映了治疗计划对应的剂量分布在患者体内的实际情况。
理论计算数据是指在放射治疗之前，基于治疗计划，使用数学计算或模拟算法预测的剂量分布数据。这些数据是根据理论模型和参数计算得到的，用于比较和评估实际测量数据。例如，使用计算算法或蒙特卡罗模拟方法，预测治疗计划对应的粒子在患者体内各位置的剂量分布数据。
伽马通过率（Gamma通过率，即γ通过率）是用于评估实际测量数据与理论计算数据之间的一致性的指标。它表示实际测量数据与理论计算数据之间的差异程度，用于质量保证过程中的剂量验证。例如，假设通过统计分析比较实际测量数据和理论计算数据，确定伽马通过率为95%，表示95%的剂量数据在实际测量和理论计算之间具有一致性。所采用的基准例如是5%/3mm、3%/3mm、2%/2mm等。
预设的通过率条件是指在质量保证过程中，设定的用于判断治疗计划质量的通过率条件。它表示实际测量数据与理论计算数据之间的一致性要求。例如设定预设的通过率条件为：3%/3mm基准下伽马通过率大于90%，意味着实际测量数据与理论计算数据的一致性需达到90%以上。
根据伽马通过率的评估，判断治疗计划的质量保证结果是合格或不合格。当伽马通过率满足预设的通过率条件时，质量保证结果设置为合格，反之则设置为不合格。假设预设的通过率条件为3%/3mm基准下伽马通过率大于90%，而伽马通过率为92%，则质量保证结果为合格；如果伽马通过率为85%，则质量保证结果为不合格。
由此，根据治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算实际测量数据与理论计算数据对应的伽马通过率。根据预设的通过率条件，确定伽马通过率的阈值，这个阈值可以根据治疗计划的要求和标准进行设置。根据计算得到的伽马通过率和通过率条件的比较结果，设定治疗计划的质量保证结果。如果伽马通过率满足预设的通过率条件，将质量保证结果设置为合格。如果伽马通过率不满足通过率条件，将质量保证结果设置为不合格。
这样做的好处是，通过比较实际测量数据与理论计算数据的伽马通过率，可以准确评估治疗计划的质量，这有助于确定治疗计划是否符合预期的剂量分布标准；通过使用实际测量数据与理论计算数据的伽马通过率作为指标，可以实时监测（或者说在线监测）治疗计划的实际测量剂量分布与预期目标（理论计算剂量分布）的一致性，这有助于及时发现和纠正潜在的问题，确保治疗计划的准确性和安全性；通过设置通过率条件和判断质量保证结果，可以确保治疗计划的剂量分布符合预期要求，这有助于提高治疗计划的准确性和安全性，保证治疗效果。
在一些实施例中，所述方法还包括：
生成质量保证报告，在所述质量保证报告中包括所述治疗计划的质量保证结果和所述伽马通过率。
在一些实施例中，所述方法还包括：
当所述质量保证结果为不合格时，根据所述治疗计划和所述伽马通过率，获取对应的质量改进方案；
生成质量保证报告，在所述质量保证报告中包括所述治疗计划的质量保证结果、所述伽马通过率和所述质量改进方案。
质量改进方案是指在质量保证过程中，当治疗计划的质量保证结果为不合格时，根据伽马通过率和治疗计划的特定情况，制定相应的改进方案。这些方案旨在纠正质量问题、提高治疗计划的准确性和效果。举例来说，如果伽马通过率未达到预设的通过率条件，质量改进方案可以包括对剂量计算算法进行调整、优化治疗计划参数或校准光学测量设备。
质量保证报告是指在质量保证过程完成后生成的报告，用于记录和总结治疗计划的质量保证结果以及相关信息。该报告可以包括质量保证结果、伽马通过率和质量改进方案等内容。例如，患者A的质量保证报告包括患者A的姓名、治疗计划信息、质量保证结果（不合格）、伽马通过率（例如85%）以及针对不合格结果提出的质量改进方案。
在本实施例中，当质量保证结果为不合格时，根据治疗计划和伽马通过率获取相应的质量改进方案。随后，生成质量保证报告，其中包括治疗计划的质量保证结果、伽马通过率和质量改进方案等信息，以便记录和评估治疗计划的质量。
由此，根据实际测量剂量分布和理论计算剂量分布之间的比较结果，判断治疗计划的质量保证结果是否合格。如果质量保证结果为不合格，则基于质量保证结果和伽马通过率的分析，获取治疗计划对应的质量改进方案。这些方案可以包括调整治疗计划参数、优化剂量分布等措施，旨在提高治疗计划的质量和效果。根据质量保证结果、伽马通过率和质量改进方案，生成质量保证报告。该报告包括治疗计划的质量保证结果（合格或不合格）、伽马通过率的数值和质量改进方案的详细描述。
这样做的好处是，通过分析质量保证结果和伽马通过率，获取针对不合格结果的质量改进方案，这有助于指导医生进行必要的调整和优化，以提高治疗计划的质量和效果；通过生成质量保证报告，将治疗计划的质量保证结果、伽马通过率和质量改进方案整合在一起，这提供了一个系统化的评估和记录工具，方便医学专家和团队之间的交流和参考；通过周期性的质量保证报告生成和质量改进方案的实施，实现持续的质量改进，这有助于不断优化治疗计划的质量、精确性和安全性；质量保证报告中记录了质量保证结果、伽马通过率和质量改进方案的详细信息，增强了治疗计划的可追溯性，有助于质量管理和质量控制。
在一些实施例中，每种对应关系采用以下任意一种方式表示：转换系数表；转换模型；转换矩阵；转换公式。
由此，每种对应关系（包括切伦科夫辐射-粒子剂量对应关系、粒子剂量-电子剂量对应关系、电子剂量-切伦科夫辐射对应关系）均可以采用转换系数表、转换模型、转换矩阵、转换公式、转换函数的方式进行表示。这些方式用于建立不同参数之间的关系，以便进行剂量计算、剂量转换或剂量校准（即剂量标定）等操作。
转换系数表通过建立一个表格，将不同参数之间的转换系数列出，以实现参数之间的转换。每个参数的取值对应于特定的转换系数，可以直接从表格中查找和应用。
转换模型利用数学模型或函数，建立不同参数之间的转换关系。这个模型可以是线性的、非线性的、统计学模型等，根据具体情况选择合适的模型来描述参数之间的转换关系。
转换矩阵使用矩阵表示不同参数之间的线性转换关系。每个矩阵元素代表不同参数之间的转换系数，通过矩阵乘法操作，可以实现参数之间的转换。
转换公式通过数学公式或方程式，描述不同参数之间的转换关系。这些公式可以基于物理原理、经验规律或统计分析等进行推导和建立。
这样做的好处是，每种对应关系可以采用相同或不同的表示方式，根据具体需要和应用场景选择最适合的转换方法，这提供了灵活性和适应性，使得参数之间的转换更加方便和可靠；通过建立准确的转换关系，可以进行剂量计算、剂量转换或剂量校准等操作，这有助于确保治疗计划的准确性和精确性，提高放射治疗的效果；使用统一的转换方式，可以实现不同实验室、机构或设备之间的参数一致性和标准化，这有助于提高数据的比较性和可重复性，促进放射治疗的质量管理和质量保证；采用适当的转换方式，可以简化剂量计算和转换的过程，提高工作效率。
在一些实施例中，所述粒子为质子；或者，
所述粒子为质子治疗对应的任意一种粒子。
质子治疗对应的粒子例如可以包括初级粒子、次级电子、次级碎片等。初级粒子主要与物质发生电离相互作用损失能量；次级电子来源于初级粒子、次级碎片与物质的电离相互作用，通过与物质发生电离损失能量；次级碎片由初级粒子与物质发生核反应产生，其中带电的核碎片（如H-2, H-3, He-3等）通过电离相互作用损失能量，不带电的碎片如中子主要通过与物质发生弹性、非弹性碰撞损失能量。
由此，如果粒子为质子，则将针对质子的特性和物理相互作用进行处理。这包括质子的输运和剂量沉积特性等。根据这些特性，可以进行剂量计算、剂量分布分析和质量保证等操作。对于其他质子治疗所涉及的粒子，该方法同样适用。虽然粒子的物理特性和相互作用可能有所不同，但可以通过相应的物理模型和转换关系进行处理，以实现剂量计算和剂量分布分析。
这样做的好处是，该质量保证方法兼容不同类型的粒子，如质子和其他质子治疗所涉及的粒子，这增强了治疗的多样性和适应性，使放射治疗能够针对不同病例和治疗需求进行个性化处理；通过针对不同类型的粒子进行处理，可以实现客观和准确的剂量计算和评估，这有助于优化治疗计划，确保粒子的剂量分布满足治疗要求，提高治疗的精确性和效果；根据不同类型的粒子进行处理，可以实现治疗的个性化和优化，对于每种粒子，可以针对其特定的物理相互作用和剂量沉积特性进行调整和优化，以提高治疗计划的效果和患者的治疗结果；通过适用于不同粒子的治疗计划，提高放射治疗的准确性、效率和安全性，进一步提升放射治疗的应用前景。
在一个具体应用场景中，本申请实施例还提供了一种用于放射治疗的质量保证方法，所述方法包括：
获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系；所述粒子为质子；
在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；
根据所述患者光学测量信息和所述切伦科夫辐射-粒子剂量对应关系，计算得到所述患者光学测量信息对应的所述患者表面的粒子剂量；
根据所述治疗计划对应的粒子辐射条件，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息；
根据所述粒子在人体内部的剂量分布信息，生成所述粒子辐射的参考深度剂量分布曲线（IDD曲线）；
根据所述患者表面的粒子剂量和所述粒子辐射的参考深度剂量分布曲线，计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到所述实际测量数据与所述理论计算数据对应的伽马通过率；
当所述伽马通过率满足预设的通过率条件时，将所述治疗计划的质量保证结果设置为合格；
当所述伽马通过率不满足所述通过率条件时，将所述治疗计划的质量保证结果设置为不合格；
当所述质量保证结果为不合格时，根据所述治疗计划和所述伽马通过率，获取对应的质量改进方案；
生成质量保证报告，在所述质量保证报告中包括所述治疗计划的质量保证结果、所述伽马通过率和所述质量改进方案。
其中，采用以下任意一种方式获取所述切伦科夫辐射-粒子剂量对应关系：
第一种方式：
在每种粒子辐射条件下执行以下处理：使用所述光学测量设备测量模体受到所述粒子辐射而发射的切伦科夫辐射，以得到模体光学测量信息；使用剂量测量设备获取所述模体表面的粒子剂量；其中，所述模体为人体模拟组织模体；
根据多种粒子辐射条件下的模体光学测量信息和模体表面的粒子剂量，建立所述模体光学测量信息和所述模体表面的粒子剂量之间的对应关系，作为所述切伦科夫辐射-粒子剂量对应关系。
第二种方式：
在多种粒子辐射条件下，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到人体表面的粒子剂量和次级电子剂量之间的粒子剂量-电子剂量对应关系；
根据所述粒子剂量-电子剂量对应关系和预设的电子剂量-切伦科夫辐射对应关系，获取所述切伦科夫辐射-粒子剂量对应关系；
其中，所述电子剂量-切伦科夫辐射对应关系用于指示切伦科夫辐射强度与电子剂量之间的正比关系。
上述实施例中，每种对应关系采用以下任意一种方式表示：转换系数表；转换模型；转换矩阵；转换公式。
作为一个示例，首先，获取切伦科夫辐射-粒子剂量对应关系。例如，如果使用质子作为粒子，可以在模体（人体模拟组织模体）上进行一系列的实验。在这些实验中，可以通过不同能量的质子辐射照射模体，然后使用光学测量设备测量模体发射的切伦科夫辐射，同时也使用剂量测量设备获取模体表面的粒子剂量。例如，当模体接收到100 Gy的质子剂量时，测量到的切伦科夫辐射强度可能是1000个光子每平方毫米。
在进行放射治疗过程中，测量患者受到的切伦科夫辐射。例如，如果在实际治疗过程中，测量到患者发射的切伦科夫辐射强度为800个光子每平方毫米，根据之前获取的切伦科夫辐射-粒子剂量对应关系，推算出患者表面的质子剂量为80 Gy。
利用蒙特卡罗方法模拟质子在人体中的输运过程。在模拟过程中，可以得到粒子在人体内部的剂量分布信息。比如，在模拟过程中，发现质子主要在深度为5-10厘米的组织中停止，造成的剂量为90 Gy。根据这些信息，可以生成粒子辐射的参考深度剂量分布曲线。
根据患者表面的粒子剂量和蒙特卡洛模拟得到的参考深度剂量分布曲线，推算出粒子在患者内部的剂量分布信息。比如，根据测量的患者表面剂量（80 Gy），推算出患者体内深度为5-10厘米的组织的剂量可能为72 Gy。
根据治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到实际测量数据与理论计算数据对应的伽马通过率。比如，如果实际测量数据（72 Gy）与理论计算数据（90 Gy）对应的伽马通过率为85%，不满足预设的通过率条件（比如3%/3mm基准下伽马通过率大于90%），则治疗计划的质量保证结果被设置为不合格，此时需要获取对应的质量改进方案，如调整治疗设备的参数或更换设备。
最后，生成质量保证报告，包括治疗计划的质量保证结果、伽马通过率和质量改进方案。这份报告可以为医生和患者提供重要的信息，帮助他们了解治疗效果，并在必要时做出调整。
本申请实施例中，人体模拟组织模体可以是完整人体的模体，也可以是人体局部组织的模体，例如各部位、各器官、各人体组织对应的模体。在获取切伦科夫辐射-粒子剂量对应关系时，可以使用多种不同类型的人体模拟组织模体，如肺、骨、脂肪和肌肉等对应的模体，以更准确地模拟实际的人体组织。由于蒙特卡罗模拟需要大量计算，可能会导致放射治疗计划的延迟。在应用中，可以部分替代性地使用深度学习模型学习质子在人体中的输运特性，能更快地生成剂量分布信息。在治疗过程中，可以使用光检测器阵列来测量切伦科夫辐射，以便更准确地计算患者的剂量分布信息。
对于预设的通过率条件，可以根据患者的具体状况或治疗过程的特定阶段动态调整通过率阈值，以提供更精确的质量保证。在治疗过程中，如果伽马通过率不满足预设的通过率条件，可以设置实时反馈系统，以人工方式或者智能方式立即调整治疗设备的参数，减少不必要的辐射暴露。对于不合格的治疗计划，除了调整设备参数或更换设备外，还可以考虑其他可能的改进方案，如增加或减少治疗次数等。
需要说明的是，在能够实现的前提下，上述各实施例中的各步骤的顺序可以进行调整，本申请不对此设限。
（放射治疗系统）
参见图2，图2是本申请实施例提供的一种放射治疗系统的结构示意图。
本申请实施例还提供一种放射治疗系统，包括：
质量保证装置，用于采用上述任一项用于放射治疗的质量保证方法，获取治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格；
确定人体剂量的装置，用于在所述质量保证结果为合格后，根据所述治疗计划确定放射治疗中使用的剂量。
由此，在放射治疗系统中，质量保证装置使用上面任一种放射治疗的质量保证方法来评估治疗计划的质量保证结果，判断是否合格或不合格。质量保证方法可以涉及剂量测量、光学测量、模拟计算等步骤，确保治疗计划的质量和精确性。在质量保证结果为合格后，使用确定人体剂量的装置来确定实际进行放射治疗时所使用的剂量。根据治疗计划和质量保证结果，该装置将计算和确定在治疗过程中要使用的剂量。这可以包括根据治疗计划参数进行剂量计算、调整和校准等步骤，以确保患者接受到正确的剂量。
这样做的好处是，通过使用质量保证装置，可以对治疗计划进行质量保证评估，确保其符合预期的质量标准，这有助于提高治疗计划的准确性、精确性和安全性；通过确定人体剂量的装置，可以根据治疗计划和质量保证结果，计算和确定患者实际接受的剂量，这有助于优化患者的剂量分布，确保患者接受到精确、适当的放射治疗剂量；通过有效的质量保证和剂量确定，放射治疗系统可以提高治疗的效果，质量保证保证了治疗计划的质量和一致性，而正确的剂量确定确保了患者获得准确的剂量分布，从而提高了治疗的精确性和有效性；通过使用质量保证装置和确定人体剂量的装置，可以增强放射治疗的安全性，质量保证过程确保了治疗计划的准确性和合规性，而剂量确定过程确保了患者接受的剂量是安全和适当的。综上所述，该放射治疗系统通过质量保证装置和确定人体剂量的装置，实现了治疗计划的质量保证和剂量确定，提高了治疗的效果和安全性。
在一个具体应用场景中，假设该放射治疗系统被用于对一个患有脑瘤的患者进行治疗。首先，使用质量保证装置进行质量保证流程。质量保证装置将在人体模拟组织模体上执行一系列的实验，例如，通过不同能量的质子辐射照射模体，然后使用光学测量设备测量模体发射的切伦科夫辐射，同时也使用剂量测量设备获取模体表面的粒子剂量。这些数据将用于建立切伦科夫辐射-质子剂量对应关系。
然后，放射治疗系统将按照治疗计划对患者进行照射。在治疗过程中，质量保证装置会实时测量患者发射的切伦科夫辐射，例如，测量结果是800个光子每平方毫米。根据切伦科夫辐射-质子剂量对应关系，计算出患者表面的质子剂量为80 Gy。然后，根据这个剂量，结合蒙特卡罗模拟得到的IDD曲线，并推算出质子在患者内部的剂量分布信息。
此后，根据治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到实际测量数据与理论计算数据对应的伽马通过率。比如，如果实际测量数据（72 Gy）与理论计算数据（75 Gy）对应的伽马通过率为97%，满足预设的通过率条件（比如95%），则治疗计划的质量保证结果被设置为合格。
在质量保证结果为合格后，确定人体剂量的装置将根据治疗计划确定放射治疗中使用的剂量。在这个例子中，因为治疗计划是合格的，所以确定人体剂量的装置将确定剂量为72 Gy，作为治疗的实际剂量。
这种放射治疗系统具有高精度和实时反馈（或者说在线反馈）的优点，可以有效地提高放射治疗的质量和患者的治疗效果。
在一些实施例中，所述放射治疗系统还包括以下一种或多种：
光学测量设备，用于测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；测量模体受到所述粒子辐射而发射的切伦科夫辐射，以得到模体光学测量信息，其中所述模体为人体模拟组织模体；
剂量测量设备，用于获取所述模体表面的粒子剂量。
在一些实施例中，所述光学测量设备包括一个或多个光检测器。当采用多个光检测器时，可以形成光检测器阵列，当部分光检测器的视野被遮挡时，其他光检测器仍然能够检测切伦科夫辐射。由于切伦科夫辐射的产生时间很短，可以认为粒子辐射源照射患者或者模体后立即产生，因此，光检测器的检测窗口可以和粒子辐射源的发射窗口相匹配，即二者至少部分重合，以实现对切伦科夫辐射的检测。作为一个示例，可以使用控制模块开启和关闭光检测器的检测功能，光检测器被操作为具有与粒子辐射源的一系列关闭时段同相的一系列关闭时段。例如，可以向光检测器提供电子选通信号，使得其检测脉冲打开和关闭以形成一系列检测窗口。此外，可以使用快门来控制光检测器的检测功能的开启和关闭，该快门通过移动以阻止光在检测窗口之间到达光检测器，以及通过移动以允许包括切伦科夫辐射在内的光在检测窗口期间到达光检测器。快门的切换效果可以通过电子、机械、光学或软件实现的方式来实现。在实际应用中，光检测器的检测窗口可以由控制模块和/或快门控制。
在一些实施例中，所述剂量测量设备包括电离室、闪烁体探测器、热释电探测器、光刺激发光剂量计和电子自旋共振剂量计中的一种或多种。
（电子设备）
本申请实施例还提供了一种电子设备，其具体实施例与上述方法实施例中记载的实施例、所达到的技术效果一致，部分内容不再赘述。
所述电子设备包括存储器和至少一个处理器，所述存储器存储有计算机程序，所述至少一个处理器被配置成执行所述计算机程序时实现以下步骤：
获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系；
在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；
根据所述患者光学测量信息和所述切伦科夫辐射-粒子剂量对应关系，计算得到所述患者光学测量信息对应的所述患者表面的粒子剂量；
根据所述患者表面的粒子剂量，计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格。
在一些实施例中，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在每种粒子辐射条件下执行以下处理：使用所述光学测量设备测量模体受到所述粒子辐射而发射的切伦科夫辐射，以得到模体光学测量信息；使用剂量测量设备获取所述模体表面的粒子剂量；其中，所述模体为人体模拟组织模体；
根据多种粒子辐射条件下的模体光学测量信息和模体表面的粒子剂量，建立所述模体光学测量信息和所述模体表面的粒子剂量之间的对应关系，作为所述切伦科夫辐射-粒子剂量对应关系。
在一些实施例中，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在多种粒子辐射条件下，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到人体表面的粒子剂量和次级电子剂量之间的粒子剂量-电子剂量对应关系；
根据所述粒子剂量-电子剂量对应关系和预设的电子剂量-切伦科夫辐射对应关系，获取所述切伦科夫辐射-粒子剂量对应关系；
其中，所述电子剂量-切伦科夫辐射对应关系用于指示切伦科夫辐射强度与电子剂量之间的正比关系。
在一些实施例中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
根据所述治疗计划对应的粒子辐射条件，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息；
根据所述粒子在人体内部的剂量分布信息，生成所述粒子辐射的参考深度剂量分布曲线；
根据所述患者表面的粒子剂量和所述粒子辐射的参考深度剂量分布曲线，计算出所述粒子在所述患者内部的剂量分布信息。
在一些实施例中，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
将所述治疗计划对应的粒子辐射条件和所述患者表面的粒子剂量输入基于深度学习的剂量分布模型，以得到所述粒子在所述患者内部的剂量分布信息。
在一些实施例中，所述至少一个处理器被配置成执行所述计算机程序时采用以下方式根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果：
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到所述实际测量数据与所述理论计算数据对应的伽马通过率；
当所述伽马通过率满足预设的通过率条件时，将所述治疗计划的质量保证结果设置为合格；
当所述伽马通过率不满足所述通过率条件时，将所述治疗计划的质量保证结果设置为不合格。
在一些实施例中，所述至少一个处理器被配置成执行所述计算机程序时还实现以下步骤：
当所述质量保证结果为不合格时，根据所述治疗计划和所述伽马通过率，获取对应的质量改进方案；
生成质量保证报告，在所述质量保证报告中包括所述治疗计划的质量保证结果、所述伽马通过率和所述质量改进方案。
在一些实施例中，每种对应关系采用以下任意一种方式表示：转换系数表；转换模型；转换矩阵；转换公式。
在一些实施例中，所述粒子为质子；或者，
所述粒子为质子治疗对应的任意一种粒子。
参见图3，图3是本申请实施例提供的一种电子设备10的结构框图。
电子设备10例如可以包括至少一个存储器11、至少一个处理器12以及连接不同平台系统的总线13。
存储器11可以包括易失性存储器形式的（计算机）可读介质，例如随机存取存储器(RAM)111和/或高速缓存存储器112，还可以进一步包括只读存储器(ROM)113。其中，存储器11还存储有计算机程序，计算机程序可以被处理器12执行，使得处理器12实现上述任一项方法的步骤。存储器11还可以包括具有至少一个程序模块115的实用工具114，这样的程序模块115包括但不限于：操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例的每一个或某种组合中可能包括网络环境的实现。
相应的，处理器12可以执行上述计算机程序，以及可以执行实用工具114。处理器12可以采用一个或多个应用专用集成电路（ASIC，Application Specific IntegratedCircuit）、DSP、可编程逻辑器件（PLD，Programmable Logic Device）、复杂可编程逻辑器件（CPLD，Complex Programmable Logic Device）、现场可编程门阵列（FPGA，Field-Programmable Gate Array）或其他电子元件。
总线13可以为表示几类总线结构的一种或多种，包括存储器总线或者存储器控制器、外围总线、图形加速端口、处理器或者使用多种总线结构的任意总线结构的局域总线。
电子设备10也可以与一个或多个外部设备例如键盘、指向设备、蓝牙设备等通信，还可与一个或者多个能够与该电子设备10交互的设备通信，和/或与使得该电子设备10能与一个或多个其它计算设备进行通信的任何设备(例如路由器、调制解调器等)通信。这种通信可以通过输入输出接口14进行。并且，电子设备10还可以通过网络适配器15与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。网络适配器15可以通过总线13与电子设备10的其它模块通信。应当明白，尽管图中未示出，但在实际应用中可以结合电子设备10使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储平台等。
（计算机可读存储介质）
本申请实施例还提供了一种计算机可读存储介质，其具体实施例与上述方法实施例中记载的实施例、所达到的技术效果一致，部分内容不再赘述。
所述计算机可读存储介质存储有计算机程序，所述计算机程序被至少一个处理器执行时实现上述任一项方法的步骤或者实现上述任一项电子设备的功能。
计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。在本申请实施例中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。计算机可读存储介质例如可以为但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
计算机可读存储介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了可读程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读存储介质还可以是任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读存储介质上包含的程序代码可以用任何适当的介质传输，包括但不限于无线、有线、光缆、RF等，或者上述的任意合适的组合。可以以一种或多种程序设计语言的任意组合来编写用于执行本发明操作的程序代码，程序设计语言包括Java、C++、Python、C#、JavaScript、PHP、Ruby、Swift、Go、Kotlin等。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中，远程计算设备可以通过任意种类的网络，包括局域网(LAN)或广域网(WAN)，连接到用户设备，或者，可以连接到外部计算设备(例如利用因特网服务提供商来通过因特网连接)。
（计算机程序产品）
本申请实施例还提供了一种计算机程序产品，其具体实施例与上述方法实施例中记载的实施例、所达到的技术效果一致，部分内容不再赘述。
本申请提供了一种计算机程序产品，所述计算机程序产品包括计算机程序，所述计算机程序被至少一个处理器执行时实现上述任一项方法的步骤或者实现上述任一项电子设备的功能。
参见图4，图4是本申请实施例提供的一种计算机程序产品的结构示意图。
所述计算机程序产品用于实现上述任一项方法的步骤或者实现上述任一项电子设备的功能。计算机程序产品可以采用便携式紧凑盘只读存储器(CD-ROM)并包括程序代码，并可以在终端设备，例如个人电脑上运行。然而，本发明的计算机程序产品不限于此，计算机程序产品可以采用一个或多个计算机可读介质的任意组合。
本申请从使用目的上，效能上，进步及新颖性等观点进行阐述，已符合专利法所强调的功能增进及使用要件，本申请以上的说明书及说明书附图，仅为本申请的较佳实施例而已，并非以此局限本申请，因此，凡一切与本申请构造，装置，特征等近似、雷同的，即凡依本申请专利申请范围所作的等同替换或修饰等，皆应属本申请的专利申请保护的范围之内。
Claims (10)
1.一种电子设备，其特征在于，所述电子设备用于放射治疗的质量保证，所述电子设备包括存储器和至少一个处理器，所述存储器存储有计算机程序，所述至少一个处理器被配置成执行所述计算机程序时实现以下步骤：
获取人体受到粒子辐射而发射的切伦科夫辐射和人体表面的粒子剂量之间的切伦科夫辐射-粒子剂量对应关系；
在治疗计划对应的放射治疗过程中，使用光学测量设备测量患者受到粒子辐射而发射的切伦科夫辐射，以得到患者光学测量信息；
根据所述患者光学测量信息和所述切伦科夫辐射-粒子剂量对应关系，计算得到所述患者光学测量信息对应的所述患者表面的粒子剂量；
根据所述患者表面的粒子剂量，计算出所述粒子在所述患者内部的剂量分布信息，作为所述治疗计划对应的剂量分布的实际测量数据；
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格；
所述至少一个处理器被配置成执行所述计算机程序时采用以下方式根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，获取所述治疗计划的质量保证结果：
根据所述治疗计划对应的剂量分布的实际测量数据与理论计算数据，计算得到所述实际测量数据与所述理论计算数据对应的伽马通过率；
当所述伽马通过率满足预设的通过率条件时，将所述治疗计划的质量保证结果设置为合格；
当所述伽马通过率不满足所述通过率条件时，将所述治疗计划的质量保证结果设置为不合格。
2.根据权利要求1所述的电子设备，其特征在于，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在每种粒子辐射条件下执行以下处理：使用所述光学测量设备测量模体受到所述粒子辐射而发射的切伦科夫辐射，以得到模体光学测量信息；使用剂量测量设备获取所述模体表面的粒子剂量；其中，所述模体为人体模拟组织模体；
根据多种粒子辐射条件下的模体光学测量信息和模体表面的粒子剂量，建立所述模体光学测量信息和所述模体表面的粒子剂量之间的对应关系，作为所述切伦科夫辐射-粒子剂量对应关系。
3.根据权利要求1所述的电子设备，其特征在于，所述切伦科夫辐射-粒子剂量对应关系的获取过程包括：
在多种粒子辐射条件下，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到人体表面的粒子剂量和次级电子剂量之间的粒子剂量-电子剂量对应关系；
根据所述粒子剂量-电子剂量对应关系和预设的电子剂量-切伦科夫辐射对应关系，获取所述切伦科夫辐射-粒子剂量对应关系；
其中，所述电子剂量-切伦科夫辐射对应关系用于指示切伦科夫辐射强度与电子剂量之间的正比关系。
4.根据权利要求1所述的电子设备，其特征在于，所述粒子在所述患者内部的剂量分布信息的获取过程包括：
根据所述治疗计划对应的粒子辐射条件，使用蒙特卡罗方法模拟所述粒子在人体中的输运过程，以得到所述治疗计划对应的所述粒子在人体内部的剂量分布信息；
根据所述粒子在人体内部的剂量分布信息，生成所述粒子辐射的参考深度剂量分布曲线；
根据所述患者表面的粒子剂量和所述粒子辐射的参考深度剂量分布曲线，计算出所述粒子在所述患者内部的剂量分布信息。
5.根据权利要求1所述的电子设备，其特征在于，所述至少一个处理器被配置成执行所述计算机程序时还实现以下步骤：
当所述质量保证结果为不合格时，根据所述治疗计划和所述伽马通过率，获取对应的质量改进方案；
生成质量保证报告，在所述质量保证报告中包括所述治疗计划的质量保证结果、所述伽马通过率和所述质量改进方案。
6.根据权利要求1-5任一项所述的电子设备，其特征在于，每种对应关系采用以下任意一种方式表示：转换系数表；转换模型；转换矩阵；转换公式。
7.根据权利要求6所述的电子设备，其特征在于，所述粒子为质子；或者，
所述粒子为质子治疗对应的任意一种粒子。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被至少一个处理器执行时实现权利要求1-7任一项所述电子设备的功能。
9.一种计算机程序产品，其特征在于，所述计算机程序产品包括计算机程序，所述计算机程序被至少一个处理器执行时实现权利要求1-7任一项所述电子设备的功能。
10.一种放射治疗系统，其特征在于，包括：
质量保证装置，用于采用权利要求1-7任一项所述的电子设备被执行的用于放射治疗的质量保证方法，获取治疗计划的质量保证结果，所述质量保证结果用于指示合格或者不合格；
确定人体剂量的装置，用于在所述质量保证结果为合格后，根据所述治疗计划确定放射治疗中使用的剂量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310809495.1A CN116832345B (zh) | 2023-07-04 | 2023-07-04 | 电子设备、放射治疗系统及相关装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310809495.1A CN116832345B (zh) | 2023-07-04 | 2023-07-04 | 电子设备、放射治疗系统及相关装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116832345A CN116832345A (zh) | 2023-10-03 |
| CN116832345B true CN116832345B (zh) | 2024-02-06 |
Family
ID=88161206
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310809495.1A Active CN116832345B (zh) | 2023-07-04 | 2023-07-04 | 电子设备、放射治疗系统及相关装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116832345B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117323584B (zh) * | 2023-10-18 | 2024-03-29 | 迈胜医疗设备有限公司 | 用于放射治疗计划调整的电子设备、放射治疗系统及相关装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104548372A (zh) * | 2015-01-07 | 2015-04-29 | 上海联影医疗科技有限公司 | 放射治疗计划、剂量确定、质量保证的方法及装置 |
| WO2016181140A1 (en) * | 2015-05-12 | 2016-11-17 | Vision Rt Limited | Radiation dosage monitoring system |
| CN109893778A (zh) * | 2019-03-15 | 2019-06-18 | 南京航空航天大学 | 一种粒子放射治疗射束实时监测方法 |
| CN111991710A (zh) * | 2020-08-24 | 2020-11-27 | 西安交通大学医学院第一附属医院 | 放疗位置和剂量实时监测及定位装置、肿瘤放射治疗系统 |
| CN114887238A (zh) * | 2022-04-15 | 2022-08-12 | 南京航空航天大学 | 一种应用切伦科夫光的外照射放射治疗在线监测方法与装置 |
| CN114949630A (zh) * | 2022-04-15 | 2022-08-30 | 南京航空航天大学 | 一种碳量子点薄膜制备方法及其在切伦科夫剂量监测中的应用 |
| CN115393256A (zh) * | 2022-04-15 | 2022-11-25 | 南京航空航天大学 | 一种基于神经网络和切伦科夫光的辐照剂量测量方法与装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018003179A1 (ja) * | 2016-06-28 | 2018-01-04 | 株式会社日立製作所 | 放射線治療装置 |
-
2023
- 2023-07-04 CN CN202310809495.1A patent/CN116832345B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104548372A (zh) * | 2015-01-07 | 2015-04-29 | 上海联影医疗科技有限公司 | 放射治疗计划、剂量确定、质量保证的方法及装置 |
| WO2016181140A1 (en) * | 2015-05-12 | 2016-11-17 | Vision Rt Limited | Radiation dosage monitoring system |
| CN107636492A (zh) * | 2015-05-12 | 2018-01-26 | 维申Rt有限公司 | 辐射剂量监测系统 |
| CN109893778A (zh) * | 2019-03-15 | 2019-06-18 | 南京航空航天大学 | 一种粒子放射治疗射束实时监测方法 |
| CN111991710A (zh) * | 2020-08-24 | 2020-11-27 | 西安交通大学医学院第一附属医院 | 放疗位置和剂量实时监测及定位装置、肿瘤放射治疗系统 |
| CN114887238A (zh) * | 2022-04-15 | 2022-08-12 | 南京航空航天大学 | 一种应用切伦科夫光的外照射放射治疗在线监测方法与装置 |
| CN114949630A (zh) * | 2022-04-15 | 2022-08-30 | 南京航空航天大学 | 一种碳量子点薄膜制备方法及其在切伦科夫剂量监测中的应用 |
| CN115393256A (zh) * | 2022-04-15 | 2022-11-25 | 南京航空航天大学 | 一种基于神经网络和切伦科夫光的辐照剂量测量方法与装置 |
Non-Patent Citations (1)
| Title |
|---|
| 基于切伦科夫辐射的核素治疗剂量测量可行性分析;艾遥;舒迪昀;龚春慧;耿长冉;张旭东;汤晓斌;;强激光与粒子束(第12期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116832345A (zh) | 2023-10-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kry et al. | AAPM TG 158: measurement and calculation of doses outside the treated volume from external‐beam radiation therapy | |
| Hälg et al. | Neutron dose and its measurement in proton therapy—current state of knowledge | |
| Sadagopan et al. | Characterization and clinical evaluation of a novel IMRT quality assurance system | |
| Ojala et al. | The accuracy of Acuros XB algorithm for radiation beams traversing a metallic hip implant—comparison with measurements and Monte Carlo calculations | |
| Lansonneur et al. | Simulation and experimental validation of a prototype electron beam linear accelerator for preclinical studies | |
| CN116832345B (zh) | 电子设备、放射治疗系统及相关装置 | |
| Jones et al. | Characterization of Compton-scatter imaging with an analytical simulation method | |
| CN111001097B (zh) | 一种放疗剂量评估系统、装置及存储介质 | |
| Mijnheer | EPID-based dosimetry and its relation to other 2D and 3D dose measurement techniques in radiation therapy | |
| Jayakody et al. | Dosimetry procedure to verify dose in High Dose Rate (HDR) brachytherapy treatment of cancer patients: A systematic review | |
| Jung et al. | Fast transit portal dosimetry using density‐scaled layer modeling of aSi‐based electronic portal imaging device and Monte Carlo method | |
| Juste et al. | Comparison of experimental 3D dose curves in a heterogeneous phantom with results obtained by MCNP5 simulation and those extracted from a commercial treatment planning system | |
| Cui et al. | A comparison between GATE and MCNPX for photon dose calculations in radiation protection using a male voxel phantom | |
| CN117323584B (zh) | 用于放射治疗计划调整的电子设备、放射治疗系统及相关装置 | |
| Gilling et al. | Organ dose from Varian XI and Varian OBI systems are clinically comparable for pelvic CBCT imaging | |
| Gilling | A GATE Monte Carlo dose analysis from Varian XI conebeam computed tomography | |
| Cornelius et al. | Commissioning of a Geant4 based treatment plan simulation tool: linac model and dicom-rt interface | |
| Moreno | Monte Carlo simulations for dosimetric verification in photon and electron beam radiotherapy | |
| Tan | 2D transit dosimetry using electronic portal imaging device | |
| Farrag et al. | Photo-neutron production of medical linear accelerator Elekta, precise model: A study with the Geant4 MC toolkit | |
| Nabha | Biophysical characterization of primary and secondary radiation fields in proton therapy | |
| Scarboro | Understanding the Influence of Photon Energy on 6MV Non Reference Dosimetry and CT Dosimetry using TLD and OSLD | |
| Ba Sunbul | Novel Developments and Applications of Dosimetry in Radiation Therapy Treatment | |
| Dukov et al. | Realistic breast phantoms with segmented real tumour formations from tomographic images | |
| Botha | The evaluation of an algorithmic model, created for the image guided radiotherapy treatment couch for integration into the Pinnacle Treatment Planning System |
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 |