[go: up one dir, main page]

CN116603176A - 硼中子捕获治疗系统及其照射参数选取方法 - Google Patents

硼中子捕获治疗系统及其照射参数选取方法 Download PDF

Info

Publication number
CN116603176A
CN116603176A CN202310597333.6A CN202310597333A CN116603176A CN 116603176 A CN116603176 A CN 116603176A CN 202310597333 A CN202310597333 A CN 202310597333A CN 116603176 A CN116603176 A CN 116603176A
Authority
CN
China
Prior art keywords
boron
voxel
dose
boron concentration
tissue model
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.)
Pending
Application number
CN202310597333.6A
Other languages
English (en)
Inventor
邓逸樵
陈江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Boron Xiamen Medical Equipment Co ltd
Original Assignee
China Boron Xiamen Medical Equipment Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Boron Xiamen Medical Equipment Co ltd filed Critical China Boron Xiamen Medical Equipment Co ltd
Priority to CN202310597333.6A priority Critical patent/CN116603176A/zh
Publication of CN116603176A publication Critical patent/CN116603176A/zh
Pending legal-status Critical Current

Links

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
    • 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
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Urology & Nephrology (AREA)
  • Surgery (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Radiation-Therapy Devices (AREA)

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)、放射性核素半衰期(HalfTime)和影像晶格强度(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)、放射性核素半衰期(HalfTime)和影像晶格强度(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)、放射性核素半衰期(HalfTime)，可以理解，还可以包括药物种类(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 (11)

1.一种硼中子捕获治疗系统，其特征在于，包括：
中子束照射装置，所述中子束照射装置用于在照射治疗时产生治疗用中子束并照射到摄入含硼药物的被照射体形成被照射部位；
治疗计划模块，所述治疗计划模块根据所述被照射部位的医学影像数据定义组织种类并建立带有组织种类数据的三维体素假体组织模型，所述被照射部位的医学影像数据包括组织相关信息和硼浓度相关信息，所述治疗计划模块根据所述组织相关信息建立对应的带有组织种类数据的三维体素假体组织模型，所述治疗计划模块对所述三维体素假体组织模型的不同种类的组织给定硼浓度数据并能够对同一种类的组织给定不同的硼浓度数据，所述治疗计划模块根据所述三维体素假体组织模型、所述硼浓度数据和所述中子束照射装置产生的所述治疗用中子束的参数进行剂量模拟计算并生成治疗计划；
控制模块，所述控制模块从所述治疗计划模块调取与所述被照射体对应的所述治疗计划，并控制所述中子束照射装置按照所述治疗计划对所述被照射体进行照射治疗。
2.根据权利要求1所述的硼中子捕获治疗系统，其特征在于：所述治疗计划模块根据所述硼浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。
3.根据权利要求2所述的硼中子捕获治疗系统，其特征在于：所述组织相关信息通过所述被照射部位的非放射性核素医学影像获得，所述治疗计划模块根据所述非放射性核素医学影像数据与组织种类之间的转换关系自动定义或手动定义所述三维体素假体组织模型中每一体素单元的组织种类。
4.根据权利要求2所述的硼中子捕获治疗系统，其特征在于：所述硼浓度相关信息通过所述被照射部位的放射性核素医学影像获得，所述被照射体摄入放射性标记的含硼药物或与含硼药物的肿瘤细胞亲和性近似的非含硼药物用于进行所述放射性核素医学影像扫描，所述治疗计划模块根据所述放射性核素医学影像数据与硼浓度之间的转换关系自动定义或手动定义所述三维体素假体组织模型中每一体素单元的硼浓度。
5.根据权利要求2所述的硼中子捕获治疗系统，其特征在于：所述硼浓度相关信息包括被照射体体重(Body Weight)、施打药物药量(Injection Dose)、药物活性测量时间(Measure Time)、造影时间(Scan Time)、放射性核素半衰期(Half Time)和影像晶格强度(Image Pixel Intensitypixel)；
所述治疗计划模块根据所述硼浓度相关信息计算所述三维体素假体组织模型中每一体素单元的硼浓度，包括：
采用公式一计算近似的血液影像晶格强度Image Pixel Intensityblood：
其中，SUVblood为血液的标准摄取值，Calibration Factor为医学影像扫描设备的校正值；
采用公式二计算所述三维体素假体组织模型中每一体素单元的硼浓度对血液硼浓度的比值：
其中，Bpixel为所述三维体素假体组织模型中每一体素单元的硼浓度，Bblood为血液硼浓度，SUVpixel为所述三维体素假体组织模型中每一体素单元的标准摄取值；
所述治疗计划模块根据计算结果给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。
6.根据权利要求5所述的硼中子捕获治疗系统，其特征在于：所述治疗计划模块通过蒙特卡罗模拟程序模拟所述三维体素假体组织模型单位时间的硼剂量(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γ为光子的相对生物效应。
7.根据权利要求2所述的硼中子捕获治疗系统，其特征在于：所述治疗计划模块根据所述中子束照射装置产生的所述治疗用中子束的参数和带有所述组织种类、硼浓度数据的所述三维体素假体组织模型通过蒙特卡罗模拟程序模拟所述被照射部位进行所述治疗用中子束的照射治疗时的物理剂量率分布。
8.根据权利要求7所述的硼中子捕获治疗系统，其特征在于：所述治疗计划模块根据不同照射角度取样模拟计算的等效剂量率分布进行评估，并选择出至少一个照射角度。
9.一种硼中子捕获治疗系统的照射参数选取方法，其特征在于，包括：
建立带有组织种类数据的三维体素假体组织模型；
根据硼浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼浓度数据；
在蒙特卡罗模拟程序中定义射束参数，通过对不同照射角度取样进行剂量模拟计算；
根据计算结果选择最佳的照射角度。
10.根据权利要求9所述的照射参数选取方法，其特征在于：所述根据被照射部位的医学影像数据的硼浓度相关信息给定所述三维体素假体组织模型中每一体素单元的硼浓度数据的步骤，包括：
解读被照射部位的医学影像数据的硼浓度相关信息，获得被照射体体重(BodyWeight)、施打药物药量(Injection Dose)、药物活性测量时间(Measure Time)、造影时间(Scan Time)、放射性核素半衰期(HalfTime)和影像晶格强度(Image PixelIntensitypixel)；
计算所述三维体素假体组织模型中每一体素单元的硼浓度，包括：
采用公式一计算近似的血液影像晶格强度Image Pixel Intensityblood：
其中，SUVblood为血液的标准摄取值，Calibration Factor为医学影像扫描设备的校正值；
采用公式二计算所述三维体素假体组织模型中每一体素单元的硼浓度Bpixel对血液硼浓度Bblood的比值：
其中，SUVpixel为所述三维体素假体组织模型中每一体素单元的标准摄取值；
根据计算结果给定所述三维体素假体组织模型中每一体素单元的硼浓度数据。
11.根据权利要求10所述的照射参数选取方法，其特征在于：所述在蒙特卡罗模拟程序中定义射束参数，通过对不同照射角度取样进行剂量模拟计算的步骤包括：
模拟在取样的照射角度下，所述三维体素假体组织模型中每一体素单元在所定义的射束照射下单位时间所受到的物理剂量，包括硼剂量(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γ为光子的相对生物效应。
CN202310597333.6A 2020-10-14 2020-10-14 硼中子捕获治疗系统及其照射参数选取方法 Pending CN116603176A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310597333.6A CN116603176A (zh) 2020-10-14 2020-10-14 硼中子捕获治疗系统及其照射参数选取方法

Applications Claiming Priority (2)

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 硼中子捕获治疗系统及其照射参数选取方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN202011094034.3A Division CN114367061B (zh) 2020-10-14 2020-10-14 硼中子捕获治疗系统及其治疗计划生成方法

Publications (1)

Publication Number Publication Date
CN116603176A true CN116603176A (zh) 2023-08-18

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 After (1)

Application Number Title Priority Date Filing Date
CN202011094034.3A Active CN114367061B (zh) 2020-10-14 2020-10-14 硼中子捕获治疗系统及其治疗计划生成方法

Country Status (6)

Country Link
US (1) US20230241414A1 (zh)
EP (1) EP4230259A4 (zh)
JP (2) JP2023540347A (zh)
CN (2) CN116603176A (zh)
TW (1) TWI762424B (zh)
WO (1) WO2022078175A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119028424A (zh) * 2024-01-03 2024-11-26 华硼中子科技(杭州)有限公司 一种硼中子俘获治疗相对生物学效应预测方法
CN119049581A (zh) * 2024-10-31 2024-11-29 华硼中子科技(杭州)有限公司 用于bnct的中子源适用性评估方法和装置、电子设备

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113877073B (zh) 2020-07-03 2023-09-12 中硼(厦门)医疗器械有限公司 放射治疗系统及其治疗计划生成方法
JP2025534781A (ja) * 2022-10-24 2025-10-17 中硼（厦▲門▼）医▲療▼器械有限公司 ホウ素中性子捕捉療法システム及びその動作方法
CN117982811A (zh) * 2022-11-07 2024-05-07 中硼(厦门)医疗器械有限公司 硼中子捕获治疗系统及照射剂量修正方法
CN118022199A (zh) * 2022-11-07 2024-05-14 中硼(厦门)医疗器械有限公司 硼中子捕获治疗系统及其工作方法
CN117982810A (zh) * 2022-11-07 2024-05-07 中硼(厦门)医疗器械有限公司 中子捕获治疗系统及剂量评估方法
CN115831319A (zh) * 2022-11-30 2023-03-21 西安交通大学医学院第一附属医院 基于核素显像确定分化型甲癌131碘治疗剂量的方法及系统
CN116598010B (zh) * 2023-05-16 2024-01-09 兰州大学 硫中子俘获治疗中载能粒子对正常细胞的影响分析方法
CN118136279A (zh) * 2024-03-13 2024-06-04 华硼中子科技(杭州)有限公司 一种微剂量学模型与蒙特卡罗中子输运模型耦合方法
CN118211421B (zh) * 2024-04-12 2025-07-11 华硼中子科技(杭州)有限公司 硼中子俘获多角度照射时间系数获取方法及装置、存储介质、终端和计算机程序产品
CN118483404B (zh) * 2024-04-29 2025-03-18 中山大学肿瘤防治中心（中山大学附属肿瘤医院、中山大学肿瘤研究所） 一种放疗射线的相对生物效应的评估方法
CN118892625B (zh) * 2024-10-09 2025-03-28 华硼中子科技(杭州)有限公司 Bnct放射剂量获取方法、放射中止时刻的获取方法、终端及介质
CN119993514A (zh) * 2025-01-16 2025-05-13 华硼中子科技(杭州)有限公司 基于黑色素瘤的硼浓度分布模型获取方法、终端及介质
CN120087166B (zh) * 2025-05-06 2025-07-15 华硼中子科技(杭州)有限公司 中子运输模拟方法、装置、介质及设备与bnct治疗系统
CN120393316B (zh) * 2025-06-25 2025-10-03 华硼中子科技(杭州)有限公司 基于屏障算法的bnct治疗计划生成方法、装置、介质与终端
CN120376048B (zh) * 2025-06-27 2025-09-26 华硼中子科技(杭州)有限公司 Bnct治疗计划优化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1658878A1 (en) * 2004-11-17 2006-05-24 The European Community, represented by the European Commission BNCT treatment planning
CN106853272A (zh) * 2015-12-08 2017-06-16 南京中硼联康医疗科技有限公司 射束的照射角度评价方法
CN109011221A (zh) * 2018-09-04 2018-12-18 东莞东阳光高能医疗设备有限公司 一种剂量引导的中子俘获治疗系统及其操作方法
CN109308733A (zh) * 2017-07-27 2019-02-05 南京中硼联康医疗科技有限公司 基于医学影像数据的几何模型建立方法及剂量计算方法

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU3320995A (en) * 1994-08-11 1996-03-07 University Of Utah Research Foundation A method for quantitation of boron concentration by magnetic resonance imaging
JP5641503B2 (ja) * 2010-10-15 2014-12-17 国立大学法人 筑波大学 マルチステップ・ラティス・ボクセル法
EP3032927B1 (en) * 2014-12-08 2016-11-09 Neuboron Medtech Ltd. A beam shaping assembly for neutron capture therapy
JP6538423B2 (ja) * 2015-05-25 2019-07-03 株式会社東芝 ホウ素中性子捕捉療法用治療装置及びその制御方法
JP6490535B2 (ja) * 2015-08-12 2019-03-27 株式会社東芝 照射線量計算誤差推定装置および照射線量計算誤差推定方法
CN108310683A (zh) * 2015-11-17 2018-07-24 南京中硼联康医疗科技有限公司 基于医学影像数据的几何模型建立方法
EP3395407B9 (en) * 2016-01-15 2020-04-08 Neuboron Medtech Ltd. Neutron capture therapy system
JP6602714B2 (ja) * 2016-03-29 2019-11-06 住友重機械工業株式会社 中性子捕捉療法システム
KR101862645B1 (ko) * 2016-10-31 2018-06-01 가천대학교 산학협력단 Mri를 이용한 bnct를 위한 붕소의 농도분포 측정방법, 및 bnct의 치료 계획방법
CN108295384B (zh) * 2017-01-11 2020-02-28 南京中硼联康医疗科技有限公司 基于医学影像的组织元素质量比例解构方法及几何模型建立方法
WO2018168713A1 (ja) * 2017-03-13 2018-09-20 住友重機械工業株式会社 中性子捕捉療法システム及び制御装置
WO2019051557A1 (en) * 2017-09-14 2019-03-21 Australian Nuclear Science And Technology Organisation METHOD AND SYSTEM OF IRRADIATION
JP2019130077A (ja) * 2018-01-31 2019-08-08 住友重機械工業株式会社 中性子捕捉療法システム
CN209253964U (zh) * 2018-09-03 2019-08-16 东莞东阳光高能医疗设备有限公司 一种磁共振引导的硼中子俘获治疗系统
JPWO2020111085A1 (ja) * 2018-11-26 2021-10-21 国立大学法人 筑波大学 ホウ素中性子捕捉療法（ｂｎｃｔ）用治療計画システム、ｂｎｃｔシステム及びｂｎｃｔ
CN110013613A (zh) * 2019-04-16 2019-07-16 东莞东阳光高能医疗设备有限公司 一种硼中子俘获治疗计划系统
CN113877073B (zh) * 2020-07-03 2023-09-12 中硼(厦门)医疗器械有限公司 放射治疗系统及其治疗计划生成方法
CN114073827B (zh) * 2020-08-15 2023-08-04 中硼(厦门)医疗器械有限公司 放射线照射系统及其控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1658878A1 (en) * 2004-11-17 2006-05-24 The European Community, represented by the European Commission BNCT treatment planning
CN106853272A (zh) * 2015-12-08 2017-06-16 南京中硼联康医疗科技有限公司 射束的照射角度评价方法
CN109308733A (zh) * 2017-07-27 2019-02-05 南京中硼联康医疗科技有限公司 基于医学影像数据的几何模型建立方法及剂量计算方法
CN109011221A (zh) * 2018-09-04 2018-12-18 东莞东阳光高能医疗设备有限公司 一种剂量引导的中子俘获治疗系统及其操作方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119028424A (zh) * 2024-01-03 2024-11-26 华硼中子科技(杭州)有限公司 一种硼中子俘获治疗相对生物学效应预测方法
CN119028424B (zh) * 2024-01-03 2025-01-21 华硼中子科技(杭州)有限公司 一种硼中子俘获治疗相对生物学效应预测方法
CN119049581A (zh) * 2024-10-31 2024-11-29 华硼中子科技(杭州)有限公司 用于bnct的中子源适用性评估方法和装置、电子设备

Also Published As

Publication number Publication date
TWI762424B (zh) 2022-04-21
WO2022078175A1 (zh) 2022-04-21
EP4230259A4 (en) 2024-03-20
CN114367061A (zh) 2022-04-19
CN114367061B (zh) 2023-06-23
JP2023540347A (ja) 2023-09-22
US20230241414A1 (en) 2023-08-03
JP2025114748A (ja) 2025-08-05
EP4230259A1 (en) 2023-08-23
TW202214319A (zh) 2022-04-16

Similar Documents

Publication Publication Date Title
CN114367061B (zh) 硼中子捕获治疗系统及其治疗计划生成方法
EP3357537A1 (en) Geometric model establishment method based on medical image data
JP7605866B2 (ja) 放射線療法システム及びその治療計画の生成方法
US12172031B2 (en) Radiotherapy system and therapy plan generation method therefor
US20250269203A1 (en) Boron neutron capture therapy system and treatment plan generation method thereof
EP1658878A1 (en) BNCT treatment planning
RU2824926C1 (ru) Система бор-нейтронозахватной терапии и способ создания плана лечения для нее
Amato et al. Applications of the monte carlo method in medical physics
RU2850118C1 (ru) Способ планирования бор-нейтронозахватной терапии и оценки результата ее проведения
TWI857841B (zh) 硼中子捕獲治療系統及其工作方法
TWI903287B (zh) 蒙特卡羅計算的優化方法及中子捕獲治療系統
Ssentongo Simulation of prompt gamma production in particle therapy and beam trajectory imaging
Kumada et al. Treatment planning
Stichelbaut et al. Radiation protection studies for proton therapy centres
Robinson Verification of direct brachytherapy dosimetry for a single seed implant
TW202423499A (zh) 蒙特卡羅計算的優化方法及中子捕獲治療系統
Stabin Dosimetry of lonizing

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