CN114367061B - 硼中子捕获治疗系统及其治疗计划生成方法 - Google Patents
硼中子捕获治疗系统及其治疗计划生成方法 Download PDFInfo
- Publication number
- CN114367061B CN114367061B CN202011094034.3A CN202011094034A CN114367061B CN 114367061 B CN114367061 B CN 114367061B CN 202011094034 A CN202011094034 A CN 202011094034A CN 114367061 B CN114367061 B CN 114367061B
- Authority
- CN
- China
- Prior art keywords
- boron
- neutron
- voxel
- dose
- boron concentration
- 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
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 title claims abstract description 177
- 229910052796 boron Inorganic materials 0.000 title claims abstract description 177
- 238000011282 treatment Methods 0.000 title claims abstract description 100
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000003814 drug Substances 0.000 claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 229940079593 drug Drugs 0.000 claims abstract description 24
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 230000001225 therapeutic effect Effects 0.000 claims abstract description 19
- 230000001678 irradiating effect Effects 0.000 claims abstract description 7
- 238000013334 tissue model Methods 0.000 claims description 64
- 239000008280 blood Substances 0.000 claims description 38
- 210000004369 blood Anatomy 0.000 claims description 38
- 238000002560 therapeutic procedure Methods 0.000 claims description 37
- 230000004071 biological effect Effects 0.000 claims description 23
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 238000005070 sampling Methods 0.000 claims description 13
- 210000004881 tumor cell Anatomy 0.000 claims description 11
- 230000037396 body weight Effects 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 2
- 239000002245 particle Substances 0.000 description 13
- 230000005855 radiation Effects 0.000 description 10
- 206010028980 Neoplasm Diseases 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 6
- 238000001959 radiotherapy Methods 0.000 description 6
- 238000002595 magnetic resonance imaging Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000004215 lattice model Methods 0.000 description 3
- 229910052751 metal Inorganic materials 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 201000010915 Glioblastoma multiforme Diseases 0.000 description 2
- 229910052790 beryllium Inorganic materials 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 238000002512 chemotherapy Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 208000005017 glioblastoma Diseases 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 229910052744 lithium Inorganic materials 0.000 description 2
- 201000001441 melanoma Diseases 0.000 description 2
- 238000011275 oncology therapy Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- 230000002105 relative biological effectiveness Effects 0.000 description 2
- 230000009885 systemic effect Effects 0.000 description 2
- 208000024827 Alzheimer disease Diseases 0.000 description 1
- ZOXJGFHDIHLPTG-BJUDXGSMSA-N Boron-10 Chemical compound [10B] ZOXJGFHDIHLPTG-BJUDXGSMSA-N 0.000 description 1
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 1
- 231100000987 absorbed dose Toxicity 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- ATBAMAFKBVZNFJ-UHFFFAOYSA-N beryllium atom Chemical compound [Be] ATBAMAFKBVZNFJ-UHFFFAOYSA-N 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- -1 carbon ions Chemical class 0.000 description 1
- 229910017052 cobalt Inorganic materials 0.000 description 1
- 239000010941 cobalt Substances 0.000 description 1
- GUTLYIVDDKVIGB-UHFFFAOYSA-N cobalt atom Chemical compound [Co] GUTLYIVDDKVIGB-UHFFFAOYSA-N 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 230000004992 fission Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 239000012216 imaging agent Substances 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000002727 particle therapy Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000002661 proton therapy Methods 0.000 description 1
- 206010039073 rheumatoid arthritis Diseases 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 239000013077 target material Substances 0.000 description 1
Images
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/103—Treatment planning systems
-
- 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/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
-
- 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/103—Treatment planning systems
- A61N5/1039—Treatment planning systems using functional images, e.g. PET or MRI
-
- 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/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
-
- 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
-
- 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/1077—Beam delivery systems
-
- 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
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- 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
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
- A61N2005/1034—Monte Carlo type methods; particle tracking
-
- 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
- A61N2005/1074—Details of the control system, e.g. user interfaces
-
- 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/109—Neutrons
-
- 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/1092—Details
- A61N2005/1098—Enhancing the effect of the particle by an injected agent or implanted device
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Abstract
一种硼中子捕获治疗系统及其治疗计划生成方法，硼中子捕获治疗系统包括中子束照射装置、治疗计划模块和控制模块。中子束照射装置用于在照射治疗时产生治疗用中子束并照射到摄入含硼(10B)药物的被照射体形成被照射部位。治疗计划模块根据被照射部位的医学影像数据和中子束照射装置产生的治疗用中子束的参数进行剂量模拟计算并生成治疗计划，被照射部位的医学影像数据包括组织相关信息和硼(10B)浓度相关信息。控制模块从治疗计划模块调取与被照射体对应的治疗计划，并控制中子束照射装置按照治疗计划对被照射体进行照射治疗。本发明的硼中子捕获治疗系统及其治疗计划生成方法，可以提高模型建立及剂量计算的精确性。
Description
技术领域
本发明一方面涉及一种放射治疗系统，尤其涉及一种硼中子捕获治疗系统；本发明另一方面涉及一种治疗计划生成方法，尤其涉及一种硼中子捕获治疗系统的治疗计划生成方法。
背景技术
随着原子科学的发展，例如钴六十、直线加速器、电子射束等放射线治疗已成为癌症治疗的主要手段之一。然而传统光子或电子治疗受到放射线本身物理条件的限制，在杀死肿瘤细胞的同时，也会对射束途径上大量的正常组织造成伤害；另外由于肿瘤细胞对放射线敏感程度的不同，传统放射治疗对于较具抗辐射性的恶性肿瘤(如：多行性胶质母细胞瘤(glioblastoma multiforme)、黑色素细胞瘤(melanoma))的治疗成效往往不佳。
为了减少肿瘤周边正常组织的辐射伤害，化学治疗(chemotherapy)中的标靶治疗概念便被应用于放射线治疗中；而针对高抗辐射性的肿瘤细胞，目前也积极发展具有高相对生物效应(relative biological effectiveness,RBE)的辐射源，如质子治疗、重粒子治疗、中子捕获治疗等。其中，中子捕获治疗便是结合上述两种概念，如硼中子捕获治疗(Boron Neutron Capture Therapy,BNCT)，借由含硼药物在肿瘤细胞的特异性集聚，配合精准的射束调控，提供比传统放射线更好的癌症治疗选择。
三维模型广泛应用于科学实验分析、科学实验模拟领域。比如在核辐射与防护领域，为了模拟人体在一定辐射条件下的吸收剂量以帮助医生制定治疗计划，常常需要利用计算机技术对医学影像数据进行各种处理建立精确的蒙特卡罗软件需要的晶格模型，并结合蒙特卡罗软件进行模拟计算。在中子捕获治疗领域根据医学影像数据建立蒙特卡罗软件所需的晶格模型并进行剂量计算和评估时，需要在模型中定义每个晶格反映的生物体基本信息，如组织种类、硼浓度信息等，信息的准确性和精确度决定了剂量计算结果的可靠度。而通常硼浓度信息是根据血样检验或切片检验获得样本的硼浓度数据，以此推算相应的组织和肿瘤硼浓度，从而在对应的模型区域中给定区域硼浓度值，这样给定的硼浓度信息没有考虑生物体内真实的硼药分布及硼药随时间的代谢情况，进而影响剂量计算结果的可靠度。
因此，有必要提出一种硼中子捕获治疗系统及其治疗计划生成方法。
发明内容
为了克服现有技术的缺陷，本发明的一个方面提供一种硼中子捕获治疗系统，包括中子束照射装置、治疗计划模块和控制模块。中子束照射装置用于在照射治疗时产生治疗用中子束并照射到摄入含硼(10B)药物的被照射体形成被照射部位。治疗计划模块根据所述被照射部位的医学影像数据和所述中子束照射装置产生的所述治疗用中子束的参数进行剂量模拟计算并生成治疗计划，所述被照射部位的医学影像数据包括组织相关信息和硼(10B)浓度相关信息。控制模块从所述治疗计划模块调取与所述被照射体对应的所述治疗计划，并控制所中子束照射装置按照所述治疗计划对所述被照射体进行照射治疗。根据被照射部位的医学影像数据的硼浓度相关信息进行剂量模拟和制定治疗计划，可以提高模型建立及剂量计算的精确性。
进一步地，所述治疗计划模块根据所述组织相关信息建立对应的带有组织种类数据的三维体素假体组织模型，根据所述硼浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。根据被照射部位的医学影像数据的硼浓度相关信息给定硼浓度数据，更加符合实际情况。
作为一种优选地，所述组织相关信息通过所述被照射部位的非放射性核素医学影像获得，所述治疗计划模块根据所述非放射性核素医学影像数据与组织种类之间的转换关系自动定义或手动定义所述三维体素假体组织模型中每一体素单元的组织种类。三维体素假体组织模型是根据医学影像数据与组织种类之间的转换关系建立的，更加精确地提供组织种类(元素组成)，建立的几何模型更加匹配于医学影像数据反应出的真实情况。进一步地，所述非放射性核素医学影像为CT。
作为一种优选地，所述硼浓度相关信息通过所述被照射部位的放射性核素医学影像获得，所述被照射体摄入放射性标记的含硼药物或与含硼药物的肿瘤细胞亲和性近似的非含硼药物用于进行所述放射性核素医学影像扫描，所述治疗计划模块根据所述放射性核素医学影像数据与硼浓度之间的转换关系自动定义或手动定义所述三维体素假体组织模型中每一体素单元的硼浓度。进一步地，所述放射性核素医学影像为PET，所述放射性标记的含硼药物为18F-BPA。
作为一种优选地，所述硼浓度相关信息包括被照射体体重(Body Weight)、施打药物药量(Injection Dose)、药物活性测量时间(Measure Time)、造影时间(Scan Time)、放射性核素半衰期(Half Time)和影像晶格强度(Image Pixel Intensitypixel)，所述治疗计划模块根据所述硼浓度相关信息计算所述三维体素假体组织模型中每一体素单元的硼浓度，所述治疗计划模块根据计算结果给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。进一步地，所述治疗计划模块根据所述硼浓度相关信息计算所述三维体素假体组织模型中每一体素单元的硼浓度包括：
采用公式一计算近似的血液影像晶格强度Image Pixel Intensityblood：
其中，SUVblood为血液的标准摄取值，Calibration Factor为医学影像扫描设备的校正值；
采用公式二计算所述三维体素假体组织模型中每一体素单元的硼浓度对血液硼浓度的比值：
其中，Bpixel为所述三维体素假体组织模型中每一体素单元的硼浓度，Bblood为血液硼浓度，SUVpixel为所述三维体素假体组织模型中每一体素单元的标准摄取值。
进一步地，所述治疗计划模块通过蒙特卡罗模拟程序模拟所述三维体素假体组织模型单位时间的硼剂量(DB)、快中子剂量(Df)、超热中子剂量(Depi)、热中子剂量(Dth)及光子剂量(Dγ)，并采用公式三计算所述三维体素假体组织模型的等效剂量率D：
D(Gy-Eq)＝CBE×Bpixel(ppm)×DB(Gy/ppm)+RBEf×Df(Gy)+RBEepi×Depi(Gy)+RBEth×Dth(Gy)+RBEγ×Dγ(Gy) (公式三)
其中，CBE为含硼药物的单位浓度的化合物生物效应，RBEf为快中子的相对生物效应，RBEepi为超热中子的相对生物效应，RBEth为热中子的相对生物效应，RBEγ为光子的相对生物效应。
作为一种优选地，所述治疗计划模块根据所述中子束照射装置产生的所述治疗用中子束的参数和带有所述组织种类、硼浓度数据的所述三维体素假体组织模型通过蒙特卡罗模拟程序模拟所述被照射部位进行所述治疗用中子束的照射治疗时的物理剂量率分布。
进一步地，所述治疗计划模块根据不同照射角度取样模拟计算的所述等效剂量率分布进行优选，并选择出至少一个照射角度。
本发明的另一个方面提供一种硼中子捕获治疗系统的治疗计划生成方法，其特征在于，包括：根据被照射部位的医学影像数据的组织相关信息建立对应的带有组织种类数据的三维体素假体组织模型；根据被照射部位的医学影像数据的硼(10B)浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼(10B)浓度数据；在蒙特卡罗模拟程序中定义射束参数，通过对不同照射角度取样进行剂量模拟计算；根据计算结果对照射角度进行优选生成治疗计划。根据被照射部位的医学影像数据的硼浓度相关信息给定硼浓度数据，并进行剂量模拟和制定治疗计划，更加符合实际情况，可以提高模型建立及剂量计算的精确性。
作为一种优选地，所述根据被照射部位的医学影像数据的硼浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼浓度数据的步骤，包括：解读被照射部位的医学影像数据的硼浓度相关信息，获得被照射体体重(Body Weight)、施打药物药量(Injection Dose)、药物活性测量时间(Measure Time)、造影时间(Scan Time)、放射性核素半衰期(Half Time)和影像晶格强度(Image Pixel Intensitypixel)；计算所述三维体素假体组织模型中每一体素单元的硼浓度；根据计算结果给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。进一步地，计算所述三维体素假体组织模型中每一体素单元的硼浓度包括：
采用公式一计算近似的血液影像晶格强度Image Pixel Intensityblood：
其中，SUVblood为血液的标准摄取值，Calibration Factor为医学影像扫描设备的校正值；
采用公式二计算所述三维体素假体组织模型中每一体素单元的硼浓度Bpixel对血液硼浓度Bblood的比值：
其中，SUVpixel为所述三维体素假体组织模型中每一体素单元的标准摄取值。
进一步地，所述在蒙特卡罗模拟程序中定义射束参数，通过对不同照射角度取样进行剂量模拟计算的步骤包括：模拟在取样的照射角度下，所述三维体素假体组织模型中每一体素单元在所定义的射束照射下单位时间所受到的物理剂量，包括硼剂量(DB)、快中子剂量(Df)、超热中子剂量(Depi)、热中子剂量(Dth)和光子剂量(Dγ)；采用公式三计算所述三维体素假体组织模型中每一体素单元在所定义的射束照射下单位时间的等效剂量率D：
D(Gy-Eq)＝CBE×Bpixel(ppm)×DB(Gy/ppm)+RBEf×Df(Gy)+RBEepi×Depi(Gy)+RBEth×Dth(Gy)+RBEγ×Dγ(Gy) (公式三)
其中，CBE为含硼药物的单位浓度的化合物生物效应，RBEf为快中子的相对生物效应，RBEepi为超热中子的相对生物效应，RBEth为热中子的相对生物效应，RBEγ为光子的相对生物效应。
本发明的再一个方面提供一种硼中子捕获治疗系统，包括中子束照射装置、治疗计划模块和控制模块。中子束照射装置用于在照射治疗时产生治疗用中子束并照射到摄入含硼(10B)药物的被照射体形成被照射部位。治疗计划模块根据所述被照射部位的医学影像数据定义组织种类并建立带有组织种类数据的三维体素假体组织模型，所述治疗计划模块对所述三维体素假体组织模型的不同种类的组织给定硼浓度数据并能够对同一种类的组织给定不同的硼浓度数据，所述治疗计划模块根据所述三维体素假体组织模型、所述硼浓度数据和所述中子束照射装置产生的所述治疗用中子束的参数进行剂量模拟计算并生成治疗计划。控制模块从所述治疗计划模块调取与所述被照射体对应的所述治疗计划，并控制所中子束照射装置按照所述治疗计划对所述被照射体进行照射治疗。根据实际情况给定硼浓度数据并进行剂量模拟和制定治疗计划，尤其是对同一种类的组织可以给定不同的硼浓度数据，可以提高模型建立及剂量计算的精确性。
本发明的硼中子捕获治疗系统及其治疗计划生成方法，可以提高模型建立及剂量计算的精确性。
附图说明
图1是硼中子捕获反应示意图。
图2是10B(n,α)7Li中子捕获核反应方程式。
图3是本发明实施例的中子捕获治疗系统的框图。
图4是本发明实施例中的治疗计划模块生成治疗计划的方法的流程图。
图5是本发明实施例中的建立三维体素假体组织模型的方法的流程图。
图6是本发明实施例中的给定三维体素假体组织模型中每一体素单元的硼浓度数据的方法的流程图。
图7是本发明实施例中的剂量模拟计算的方法的流程图。
具体实施方式
下面结合附图对本发明的实施例做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
作为一种优选地，以中子捕获治疗系统及其治疗计划生成方法为本发明的实施例。下面将简单介绍一下中子捕获治疗，尤其是硼中子捕获治疗。
中子捕获治疗作为一种有效的治疗癌症的手段近年来的应用逐渐增加，其中以硼中子捕获治疗最为常见，供应硼中子捕获治疗的中子可以由核反应堆或加速器供应。本发明的实施例以加速器硼中子捕获治疗为例，加速器硼中子捕获治疗的基本组件通常包括用于对带电粒子(如质子、氘核等)进行加速的加速器、靶材与热移除系统和射束整形体，其中加速带电粒子与金属靶材作用产生中子，依据所需的中子产率与能量、可提供的加速带电粒子能量与电流大小、金属靶材的物化性等特性来挑选合适的核反应，常被讨论的核反应有7Li(p,n)7Be及9Be(p,n)9B，这两种反应皆为吸热反应。两种核反应的能量阀值分别为1.881MeV和2.055MeV，由于硼中子捕获治疗的理想中子源为keV能量等级的超热中子，理论上若使用能量仅稍高于阀值的质子轰击金属锂靶材，可产生相对低能的中子，不须太多的缓速处理便可用于临床，然而锂金属(Li)和铍金属(Be)两种靶材与阀值能量的质子作用截面不高，为产生足够大的中子通量，通常选用较高能量的质子来引发核反应。
硼中子捕获治疗(Boron Neutron Capture Therapy,BNCT)是利用含硼(10B)药物对热中子具有高捕获截面的特性，借由10B(n,α)7Li中子捕获及核分裂反应产生4He和7Li两个重荷电粒子。参照图1和图2，其分别示出了硼中子捕获反应的示意图和10B(n,α)7Li中子捕获核反应方程式，两荷电粒子的平均能量约为2.33MeV，具有高线性转移(Linear EnergyTransfer,LET)、短射程特征，α粒子的线性能量转移与射程分别为150keV/μm、8μm，而7Li重荷粒子则为175keV/μm、5μm，两粒子的总射程约相当于一个细胞大小，因此对于生物体造成的辐射伤害能局限在细胞层级，当含硼药物选择性地聚集在肿瘤细胞中，搭配适当的中子射源，便能在不对正常组织造成太大伤害的前提下，达到局部杀死肿瘤细胞的目的。
蒙特卡罗方法能够对辐照目标内部三维空间核粒子碰撞轨迹和能量分布进行精确模拟，在中子捕获治疗中，为了模拟人体在一定辐射条件下的吸收剂量以帮助医生制定治疗计划，常常需要利用计算机技术对医学影像进行各种处理建立精确的蒙特卡罗软件需要的晶格模型，并结合蒙特卡罗软件进行模拟计算。医学影像数据可以为核磁共振成像(Magnetic Resonance Imaging,MRI)、电子计算机断层扫描(Computed Tomography,CT)、正电子发射型计算机断层扫描(Positron Emission Tomography,PET)、PET-CT或X射线成像(X-Ray imaging)。但本领域技术人员熟知地，还可以使用其他的医学影像数据，只要该医学影像数据能够被转换成三维体素假体组织模型，就能够应用于本发明揭示的放射治疗系统及其治疗计划生成方法中。
参阅图3，本实施例的放射治疗系统优选为硼中子捕获治疗系统100，包括中子束照射装置10、治疗计划模块20和控制模块30。中子束照射装置10包括中子产生装置11和治疗台12，中子产生装置12用于在照射治疗时产生治疗用中子束N并照射到治疗台12上摄入含硼(10B)药物的患者形成被照射部位。治疗前，治疗计划模块20根据患者被照射部位的医学影像数据和中子产生装置11产生的治疗用中子束N的参数生成治疗计划。照射治疗时，控制模块30从治疗计划模块20调取当前患者对应的治疗计划，并根据治疗计划控制中子束照射装置10的照射。
一实施例中，患者在治疗前需摄入放射性标记的含硼(10B)药物并通过放射性核素医学影像(如PET)获得硼(10B)浓度相关信息，即患者被照射部位的医学影像数据包括组织相关信息和硼浓度相关信息，其中，组织相关信息可以是同时从放射性核素医学影像获得，也可以是从其他非放射性核素医学影像(如CT)获得。可以理解，放射性核素医学影像还可以为PET-CT等，非放射性核素医学影像还可以为MRI等，本发明对此不做具体限定。治疗计划模块20根据被照射部位的医学影像数据的组织相关信息建立对应的带有组织种类、组织密度数据三维体素假体组织模型，根据被照射部位的医学影像数据的硼浓度相关信息给定三维体素假体组织模型中每一体素单元的硼浓度数据，并根据中子产生装置11产生的治疗用中子束N的参数和带有组织种类、组织密度、硼浓度数据的三维体素假体组织模型通过蒙特卡罗模拟程序模拟计算患者进行照射治疗时的剂量分布并生成治疗计划。可以理解，也可以不定义组织密度。根据被照射部位的医学影像数据的硼浓度相关信息给定硼浓度数据，并进行剂量模拟和制定治疗计划，更加符合实际情况，可以提高模型建立及剂量计算的精确性。
请参照图4，本实施例的治疗计划模块20生成治疗计划的方法，具体包括如下步骤：
S410：根据被照射部位的医学影像数据的组织相关信息建立对应的带有组织种类、组织密度数据的三维体素假体组织模型；
S420：根据被照射部位的医学影像数据的硼浓度相关信息给定三维体素假体组织模型中每一体素单元的硼浓度数据；
S430：在蒙特卡罗模拟程序(如MCNP，Monte Carlo N Particle Transport Code)中定义射束参数，通过对不同照射角度取样进行剂量模拟计算；
S440：根据计算结果对照射角度进行优选生成治疗计划。
请参照图5，一实施例中，根据医学影像数据建立三维体素假体组织模型的步骤S410，可以进一步包括：
S510：读取医学影像数据；
S520：建立三维医学影像体素模型；
S530：定义或读取感兴趣区域的边界；
S540：定义每一体素单元的组织种类(元素组成)、组织密度，可以是根据CT影像数据与组织种类、组织密度之间的转换关系自动定义；也可以是用户手动定义，如为每个感兴趣区域的边界内的体素单元给定一个特定的组织种类和组织密度。
S550：建立三维体素假体组织模型。
三维体素假体组织模型是根据医学影像数据与组织种类、组织密度之间的转换关系建立的，更加精确地提供组织种类(元素组成)及组织密度，建立的几何模型更加匹配于医学影像数据反应出的真实情况。根据医学影像数据建立三维体素假体组织模型的详细过程，可以参照2017年03月08日公开的、公开号为CN106474634A、发明名称为“基于医学影像数据的几何模型建立方法”的专利申请，在此全文引入。
步骤S420根据被照射部位的医学影像数据的硼浓度相关信息给定三维体素假体组织模型中每一体素单元的硼浓度数据，标记有组织硼浓度信息的几何模型，便可清楚地知道，各个组织内的含硼药物浓度，然后进行中子照射模拟时，则更加真实地反应出实际情况。硼浓度相关信息通过被照射部位的放射性核素医学影像获得，被照射体摄入放射性标记的含硼药物用于进行放射性核素医学影像扫描，治疗计划模块根据放射性核素医学影像数据与硼浓度之间的转换关系自动定义或手动定义三维体素假体组织模型中每一体素单元的硼浓度。一实施例中，放射性核素医学影像为PET，被照射体摄入的放射性标记的含硼药物为18F-BPA。可以理解，18F-BPA还可以替换为其他放射性标记或其他含硼药物；也可以为放射性标记的与含硼药物的肿瘤细胞亲和性近似的非含硼药物，如18F-FDG。由于目前BPA等可以作为硼中子捕获治疗用的含硼药物价格普遍较高，通过非含硼药物代替含硼药物做治疗计划的模拟，能够极大的降低成本。
请参照图6，一实施例中，步骤S420，可以进一步包括：
S610：解读被照射部位的医学影像数据的硼浓度相关信息，即放射性核素医学影像扫描获得的硼浓度相关的医学影像数据信息IOD(Information Object Definition)。
医学影像数据通常采用DICOM(Digital Imaging and Communications inMedicine)格式，DICOM数据中与硼浓度相关的IOD包括被照射体体重(Body Weight)、施打药物药量(Injection Dose)、药物活性测量时间(Measure Time)、造影时间(Scan Time)、放射性核素半衰期(Half Time)，可以理解，还可以包括药物种类(Radiopharmaceutical)等，这些信息在开始进行放射性核素医学影像扫描时即可以确定，信息来源可以是操作者手动输入的，也可以是自动获得或调取的，如表一，列出了DICOM数据中与硼浓度相关信息的数据标签(Tag)及对应的数据名称(Description)：
表一.DICOM数据中与硼浓度相关的IOD
DICOM数据中与硼浓度相关的IOD还包括影像晶格强度Image Pixel Intensity，放射性核素医学影像扫描，如PET，利用放射性核素衰变产生的正电子与组织中的电子相遇发生正负电子对湮灭反应，此时用光电倍增管(PMT)等检测器便可以检测到由此发射出来的伽马射线，之后计算机就会形成反映正电子放射性同位素分布的横截面图像。一实施例中，采用18F-BPA-PET扫描，影像每一晶格的原始数据为18F衰变释出的正子与电子互毁作用产生光子的计数率转换为PET影像每一晶格的影像晶格强度Image Pixel Intensitypixel作为医学影像数据的输出，又18F标志于BPA，因此PET影像上的原始数据可做为10B量化的依据。可以理解，如采用放射性标记的与含硼药物的肿瘤细胞亲和性近似的非含硼药物进行放射性核素医学影像扫描，影像上的原始数据也可同样做为10B量化的依据。
S620：根据解读的硼浓度相关信息计算三维体素假体组织模型中每一体素单元的硼浓度。
标准摄取值(Standard Uptake Value,SUV)为放射性核素医学影像采集的半定量指标，是指局部组织摄取的显像剂(放射性标记的药物)的放射性活度与全身平均注射活度的比值，如公式一，通过步骤S610获得的与硼浓度相关的IOD可以得到相应晶格的SUV，放射性核素医学影像扫描的每一晶格与三维体素假体组织模型中每一体素单元一一对应。
设定血液的SUV值为1，即近似认为血液的放射性活度等于全身平均注射放射性活度。由公式二获得SUVblood＝1时的放射性核素医学影像原始数据，即近似的血液影像晶格强度Image Pixel Intensityblood。
其中，SUVblood为血液的标准摄取值，Calibration Factor为医学影像扫描设备的校正值。
同时，通过每个晶格和血液标准摄取值比值，获得每个晶格和血液硼浓度比值，即采用公式三计算三维体素假体组织模型中每一体素单元的硼浓度对血液硼浓度的比值。
其中，Bpixel为三维体素假体组织模型中每一体素单元的硼浓度，Bblood为血液硼浓度，SUVpixel为三维体素假体组织模型中每一体素单元的标准摄取值。
将三维体素假体组织模型中每一体素单元的硼浓度以血液硼浓度作为基准值做换算SUVpixel/SUVblood，即可获得每一体素单元的硼浓度对血液硼浓度比值，建立非均匀硼浓度分布。通过抽血检测等获得血硼浓度值后，即可以获得每一体素单元具体的硼浓度值。
S630：根据计算结果给定每一体素单元的硼浓度数据，可以是根据计算结果自动给定的；也可以是用户手动给定的，如根据计算结果为每个感兴趣区域的边界内的体素单元给定一个特定的硼浓度数据。
通过计算更加贴合实际情况，治疗计划模块对三维体素假体组织模型的不同种类的组织给定硼浓度数据并能够对同一种类的组织给定不同的硼浓度数据，可以提高模型建立及剂量计算的精确性。建立带有组织种类和组织硼浓度的三维体素假体组织模型后，便可以通过蒙特卡罗模拟程序模拟BNCT中患者被中子束照射时内部三维空间核粒子碰撞轨迹和能量分布，即可以模拟出物理剂量率分布，根据物理剂量率分布计算等效剂量率分布，再根据剂量指标筛选治疗计划方案。请参照图7，具体来说，步骤S430包括：
S710：在蒙特卡罗模拟程序中定义射束参数(如射束能量、强度、半径等)，通过对不同照射角度取样，模拟三维体素假体组织模型在不同照射角度的物理剂量率分布，即分别模拟在取样的照射角度下，三维体素假体组织模型中每一体素单元在所定义的射束照射下单位时间所受到的物理剂量。
取样时需要决定计算射束的起始位置及射束角度，计算中起始位置与角度的决定，可以是顺向算法或逆向算法，顺向算法中是将起始位置决定于体外位置并可依固定角度或距离间隔依序取样计算，也可以透过随机取样的方式进行；射束角度的部分，则可设定为照射点至肿瘤质心或肿瘤最深处的向量方向，具体的肿瘤端点位置可视用户需求而调整；逆向算法中，则是将起始位置决定于肿瘤范围内，其起始位置可以是肿瘤质心、最深处或肿瘤范围内随机取点，而射束角度则可以利用随机取样或依指定间隔取样的方式进行。
取样时还可以对射束角度进行筛选，如进行射束角度评价，根据评价结果选择进行后续计算的射束角度；或者在取样计算后再进行射束角度的筛选，如根据辐射剂量分布的结果或射束角度评价的结果进行筛选。射束角度的评价方法在此不再详述，可以参照2017年06月16日公开的、公开号为CN106853272A、发明名称为“射束的照射角度评价方法”的专利申请，在此全文引入。
硼中子捕获治疗中会造成剂量贡献的主要因子有三部分：
1)硼剂量DB：剂量来自组织、肿瘤中的含硼药物与中子发生10B(n,α)7Li中子捕获反应而产生的高线性能量转移α、7Li粒子。
2)中子剂量：依中子能量分为快中子剂量Df、超热中子剂量Depi及热中子剂量Dth，其中对于组织器官组成元素主要为1H(n,n’)p中子与氢的弹性散射作用产生的回跳质子造成的剂量、14N(n,p)14C作用产生的质子与回跳碳离子造成的剂量以及中子与其他元素作用的微剂量。
3)光子剂量Dγ：包含中子与屏蔽结构发生作用、与人体组织产生捕获反应而诱发的光子，而后者主要为热中子发生1H(n,γ)2H反应而产生2.22MeV之光子。
一实施例，步骤S710中，单位时间所受到的物理剂量包括通过蒙特卡罗模拟程序模拟得到的单位时间的硼剂量(DB)、快中子剂量(Df)、超热中子剂量(Depi)、热中子剂量(Dth)及光子剂量(Dγ)。
S720：计算三维体素假体组织模型在不同照射角度的等效剂率量分布。
硼中子捕获治疗中，由于光子、中子会造成的生物效应不同，所以针对快中子、超热中子、热中子和光子剂量项，分别乘上不同组织的相对生物效应(Relative biologicaleffectiveness,RBE)以求得等效剂量(equivalent dose)；而硼剂量方面，因α、7Li粒子的射程短，伤害通常局限在硼中子捕获作用发生位置，又不同种类的细胞吸收含硼药物的能力不同，使得硼浓度在体内分布不均匀，为求得等效剂量，此项剂量须乘上各组织的化合物生物效应(Compound biological effectiveness,CBE)及相应的硼浓度(在步骤S630中已定义)，即一体素单元的等效剂量率D采用公式四进行计算：
D(Gy-Eq)＝CBE×Bpixel(ppm)×DB(Gy/ppm)+RBEf×Df(Gy)+RBEepi×Depi(Gy)+RBEth×Dth(Gy)+RBEγ×Dγ(Gy) (公式四)
其中，CBE为含硼药物的单位浓度的化合物生物效应，RBEf为快中子的相对生物效应，RBEepi为超热中子的相对生物效应，RBEth为热中子的相对生物效应，RBEγ为光子的相对生物效应。
步骤S440，根据S720计算得到的等效剂量率分布，结合数学算法对不同的治疗计划方案(照射角度及对应的照射时间)进行评估或优选。如选择其中一个感兴趣区域作为目标，以该感兴趣区域的最大剂量、平均剂量或处方剂量等作为约束条件求得取样的照射角度对应的照射时间，以该照射时间获得取样的照射角度下三维体素假体组织模型的等效剂量分布，然后采用剂量体积直方图(DVH)、等剂量曲线、剂量表格等对不同照射角度及对应的照射时间在三维体素假体组织模型模拟计算得出的等效剂量分布进行评估或优选，还可以进行上文所述的照射角度评价来进行评估以供医师等操作者选出更符合需求的治疗计划方案。可以理解，还可以通过优化算法等选择至少2个照射角度依次对患者进行照射，可以手动设置照射角度具体的数量或通过算法自动获得或采用弧形连续调控，照射角度的取样可以是在患者同侧或对侧。
可以理解，上述公式一至公式四及剂量计算、评估、优选方法的一些简单变换，仍然在本发明要求保护的范围之内。
可以理解，本发明还可以应用于其他能够用硼中子捕获治疗进行治疗的疾病，如阿尔兹海默症、类风湿关节炎，则肿瘤细胞为其他病变组织；患者也可以为其他被照射体。
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，都在本发明要求保护的范围之内。
Claims (6)
1.一种硼中子捕获治疗系统，其特征在于，包括：
中子束照射装置，所述中子束照射装置用于在照射治疗时产生治疗用中子束并照射到摄入含硼药物的被照射体形成被照射部位；
治疗计划模块，所述治疗计划模块根据所述被照射部位的医学影像数据和所述中子束照射装置产生的所述治疗用中子束的参数进行剂量模拟计算并生成治疗计划，所述被照射部位的医学影像数据包括组织相关信息和硼浓度相关信息，所述治疗计划模块根据所述组织相关信息建立对应的带有组织种类的数据的三维体素假体组织模型，根据所述硼浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼浓度数据，所述硼浓度相关信息包括被照射体体重Body Weight、施打药物药量Injection Dose、药物活性测量时间Measure Time、造影时间Scan Time、放射性核素半衰期Half Time和影像晶格强度ImagePixel Indensitypixel；
控制模块，所述控制模块从所述治疗计划模块调取与所述被照射体对应的所述治疗计划，并控制所述中子束照射装置按照所述治疗计划对所述被照射体进行照射治疗；
其中，所述治疗计划模块根据所述硼浓度相关信息计算所述三维体素假体组织模型中每一体素单元的硼浓度，包括：
采用公式一计算近似的血液影像晶格强度Image Pixel Intensityblood：
其中，SUVblood为血液的标准摄取值，Calibration Factor为医学影像扫描设备的校正值；
采用公式二计算所述三维体素假体组织模型中每一体素单元的硼浓度对血液硼浓度的比值：
其中，Bpixel为所述三维体素假体组织模型中每一体素单元的硼浓度，Bblood为血液硼浓度，SUVpixel为所述三维体素假体组织模型中每一体素单元的标准摄取值；
所述治疗计划模块根据计算结果给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。
2.根据权利要求1所述的硼中子捕获治疗系统，其特征在于，所述组织相关信息通过所述被照射部位的非放射性核素医学影像获得，所述治疗计划模块根据所述非放射性核素医学影像的数据与组织种类之间的转换关系自动定义或手动定义所述三维体素假体组织模型中每一体素单元的组织种类。
3.根据权利要求1所述的硼中子捕获治疗系统，其特征在于，所述硼浓度相关信息通过所述被照射部位的放射性核素医学影像获得，所述被照射体摄入放射性标记的含硼药物或与含硼药物的肿瘤细胞亲和性近似的非含硼药物用于进行所述放射性核素医学影像扫描，所述治疗计划模块根据所述放射性核素医学影像的数据与硼浓度之间的转换关系自动定义或手动定义所述三维体素假体组织模型中每一体素单元的硼浓度。
4.根据权利要求1所述的硼中子捕获治疗系统，其特征在于，所述治疗计划模块通过蒙特卡罗模拟程序模拟所述三维体素假体组织模型单位时间的硼剂量DB、快中子剂量Df、超热中子剂量Depi、热中子剂量Dth及光子剂量Dγ，并采用公式三计算所述三维体素假体组织模型的等效剂量率D：
D(Gy-Eq)=CBE×Bpixel(ppm)×DB(Gy/ppm)+RBEf×Df(Gy)+RBEepi×Depi(Gy)+RBEth×Dth(Gy)+RBEγ×Dγ(Gy)（公式三）
其中，CBE为含硼药物的单位浓度的化合物生物效应， RBEf为快中子的相对生物效应，RBEepi为超热中子的相对生物效应，RBEth为热中子的相对生物效应，RBEγ为光子的相对生物效应。
5.根据权利要求1所述的硼中子捕获治疗系统，其特征在于，所述治疗计划模块根据所述中子束照射装置产生的所述治疗用中子束的参数和带有所述组织种类、硼浓度数据的所述三维体素假体组织模型通过蒙特卡罗模拟程序模拟所述被照射部位进行所述治疗用中子束的照射治疗时的物理剂量率分布。
6.根据权利要求5所述的硼中子捕获治疗系统，其特征在于，所述治疗计划模块根据不同照射角度取样模拟计算的等效剂量率分布进行评估，并选择出至少一个照射角度。
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011094034.3A CN114367061B (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其治疗计划生成方法 |
| CN202310597333.6A CN116603176A (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其照射参数选取方法 |
| PCT/CN2021/119903 WO2022078175A1 (zh) | 2020-10-14 | 2021-09-23 | 硼中子捕获治疗系统及其治疗计划生成方法 |
| JP2023515083A JP2023540347A (ja) | 2020-10-14 | 2021-09-23 | ホウ素中性子捕捉療法システム及びその治療計画生成方法 |
| EP21879221.6A EP4230259A4 (en) | 2020-10-14 | 2021-09-23 | BORN NEUTRON CAPTURE THERAPY SYSTEM AND TREATMENT PLAN GENERATION METHOD THEREOF |
| TW110136802A TWI762424B (zh) | 2020-10-14 | 2021-09-30 | 硼中子捕獲治療系統及其治療計畫生成方法 |
| US18/131,416 US20230241414A1 (en) | 2020-10-14 | 2023-04-06 | Boron neutron capture therapy system and treatment plan generation method therefor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011094034.3A CN114367061B (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其治疗计划生成方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310597333.6A Division CN116603176A (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其照射参数选取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114367061A CN114367061A (zh) | 2022-04-19 |
| CN114367061B true CN114367061B (zh) | 2023-06-23 |
Family
ID=81137970
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310597333.6A Pending CN116603176A (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其照射参数选取方法 |
| CN202011094034.3A Active CN114367061B (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其治疗计划生成方法 |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310597333.6A Pending CN116603176A (zh) | 2020-10-14 | 2020-10-14 | 硼中子捕获治疗系统及其照射参数选取方法 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20230241414A1 (zh) |
| EP (1) | EP4230259A4 (zh) |
| JP (1) | JP2023540347A (zh) |
| CN (2) | CN116603176A (zh) |
| TW (1) | TWI762424B (zh) |
| WO (1) | WO2022078175A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113877073B (zh) * | 2020-07-03 | 2023-09-12 | 中硼(厦门)医疗器械有限公司 | 放射治疗系统及其治疗计划生成方法 |
| CN116598010B (zh) * | 2023-05-16 | 2024-01-09 | 兰州大学 | 硫中子俘获治疗中载能粒子对正常细胞的影响分析方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1996004844A1 (en) * | 1994-08-11 | 1996-02-22 | University Of Utah Research Foundation | A method for quantitation of boron concentration by magnetic resonance imaging |
| EP1658878A1 (en) * | 2004-11-17 | 2006-05-24 | The European Community, represented by the European Commission | BNCT treatment planning |
| JP2012088771A (ja) * | 2010-10-15 | 2012-05-10 | Univ Of Tsukuba | マルチステップ・ラティス・ボクセル法 |
| JP2016214760A (ja) * | 2015-05-25 | 2016-12-22 | 株式会社東芝 | ホウ素中性子捕捉療法用治療装置及びその制御方法 |
| JP2017176354A (ja) * | 2016-03-29 | 2017-10-05 | 住友重機械工業株式会社 | 中性子捕捉療法システム |
| WO2018080290A1 (ko) * | 2016-10-31 | 2018-05-03 | 가천대학교 산학협력단 | Mri를 이용한 bnct를 위한 붕소의 농도분포 측정방법, 및 bnct의 치료 계획방법 |
| CN109011221A (zh) * | 2018-09-04 | 2018-12-18 | 东莞东阳光高能医疗设备有限公司 | 一种剂量引导的中子俘获治疗系统及其操作方法 |
| CN109308733A (zh) * | 2017-07-27 | 2019-02-05 | 南京中硼联康医疗科技有限公司 | 基于医学影像数据的几何模型建立方法及剂量计算方法 |
| CN209253964U (zh) * | 2018-09-03 | 2019-08-16 | 东莞东阳光高能医疗设备有限公司 | 一种磁共振引导的硼中子俘获治疗系统 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108310683A (zh) * | 2015-11-17 | 2018-07-24 | 南京中硼联康医疗科技有限公司 | 基于医学影像数据的几何模型建立方法 |
| CN106853272A (zh) * | 2015-12-08 | 2017-06-16 | 南京中硼联康医疗科技有限公司 | 射束的照射角度评价方法 |
| CN108295384B (zh) * | 2017-01-11 | 2020-02-28 | 南京中硼联康医疗科技有限公司 | 基于医学影像的组织元素质量比例解构方法及几何模型建立方法 |
| WO2018168713A1 (ja) * | 2017-03-13 | 2018-09-20 | 住友重機械工業株式会社 | 中性子捕捉療法システム及び制御装置 |
| CN117731964A (zh) * | 2017-09-14 | 2024-03-22 | 澳大利亚核科学和技术组织 | 辐照方法及系统 |
| CN110013613A (zh) * | 2019-04-16 | 2019-07-16 | 东莞东阳光高能医疗设备有限公司 | 一种硼中子俘获治疗计划系统 |
| CN113877073B (zh) * | 2020-07-03 | 2023-09-12 | 中硼(厦门)医疗器械有限公司 | 放射治疗系统及其治疗计划生成方法 |
| CN114073827B (zh) * | 2020-08-15 | 2023-08-04 | 中硼(厦门)医疗器械有限公司 | 放射线照射系统及其控制方法 |
-
2020
- 2020-10-14 CN CN202310597333.6A patent/CN116603176A/zh active Pending
- 2020-10-14 CN CN202011094034.3A patent/CN114367061B/zh active Active
-
2021
- 2021-09-23 WO PCT/CN2021/119903 patent/WO2022078175A1/zh unknown
- 2021-09-23 JP JP2023515083A patent/JP2023540347A/ja active Pending
- 2021-09-23 EP EP21879221.6A patent/EP4230259A4/en active Pending
- 2021-09-30 TW TW110136802A patent/TWI762424B/zh active
-
2023
- 2023-04-06 US US18/131,416 patent/US20230241414A1/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1996004844A1 (en) * | 1994-08-11 | 1996-02-22 | University Of Utah Research Foundation | A method for quantitation of boron concentration by magnetic resonance imaging |
| EP1658878A1 (en) * | 2004-11-17 | 2006-05-24 | The European Community, represented by the European Commission | BNCT treatment planning |
| JP2012088771A (ja) * | 2010-10-15 | 2012-05-10 | Univ Of Tsukuba | マルチステップ・ラティス・ボクセル法 |
| JP2016214760A (ja) * | 2015-05-25 | 2016-12-22 | 株式会社東芝 | ホウ素中性子捕捉療法用治療装置及びその制御方法 |
| JP2017176354A (ja) * | 2016-03-29 | 2017-10-05 | 住友重機械工業株式会社 | 中性子捕捉療法システム |
| WO2018080290A1 (ko) * | 2016-10-31 | 2018-05-03 | 가천대학교 산학협력단 | Mri를 이용한 bnct를 위한 붕소의 농도분포 측정방법, 및 bnct의 치료 계획방법 |
| CN109308733A (zh) * | 2017-07-27 | 2019-02-05 | 南京中硼联康医疗科技有限公司 | 基于医学影像数据的几何模型建立方法及剂量计算方法 |
| CN209253964U (zh) * | 2018-09-03 | 2019-08-16 | 东莞东阳光高能医疗设备有限公司 | 一种磁共振引导的硼中子俘获治疗系统 |
| CN109011221A (zh) * | 2018-09-04 | 2018-12-18 | 东莞东阳光高能医疗设备有限公司 | 一种剂量引导的中子俘获治疗系统及其操作方法 |
Non-Patent Citations (1)
| Title |
|---|
| BNCT治疗计划系统综述;肖刚等;《中国医学物理学杂志》;20060228(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EP4230259A4 (en) | 2024-03-20 |
| WO2022078175A1 (zh) | 2022-04-21 |
| CN116603176A (zh) | 2023-08-18 |
| TWI762424B (zh) | 2022-04-21 |
| CN114367061A (zh) | 2022-04-19 |
| EP4230259A1 (en) | 2023-08-23 |
| TW202214319A (zh) | 2022-04-16 |
| JP2023540347A (ja) | 2023-09-22 |
| US20230241414A1 (en) | 2023-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20230241414A1 (en) | Boron neutron capture therapy system and treatment plan generation method therefor | |
| EP1658878A1 (en) | BNCT treatment planning | |
| CN113877073B (zh) | 放射治疗系统及其治疗计划生成方法 | |
| Martínez‐Rovira et al. | Implementation of biological washout processes within GATE/Geant4—A Monte Carlo study in the case of carbon therapy treatments | |
| CN117982812A (zh) | 硼中子捕获治疗系统及其治疗计划生成方法 | |
| Amato et al. | Applications of the monte carlo method in medical physics | |
| Singh | Oral | |
| Ssentongo | Simulation of prompt gamma production in particle therapy and beam trajectory imaging | |
| Li et al. | PET/SPECT/Spectral‐CT/CBCT imaging in a small‐animal radiation therapy platform: A Monte Carlo study—Part II: Biologically guided radiotherapy | |
| Carriço | Monte Carlo investigations for a novel small-animal proton radiographic setup based on a pixelated Silicon detector and energy stacking | |
| Safavi-Naeini et al. | Radioactive beams for ion therapy: Monte Carlo simulations and experimental verifications | |
| Papoutsis | From dose prescription to dose delivery-can dose painting by numbers be accurately delivered? | |
| Esen et al. | Oral | |
| Ramanathan | Oral | |
| Sweetlove et al. | Oral | |
| Mishra et al. | Oral | |
| Das et al. | Oral | |
| Sha et al. | Oral | |
| Holla et al. | Oral | |
| Siwach et al. | Oral | |
| Nair et al. | Oral | |
| Cruz et al. | Oral | |
| Deshpande et al. | Oral | |
| Patel et al. | Oral | |
| Perumal et al. | Oral |
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 |