CN111803803A - 基于医学影像数据的几何模型建立方法 - Google Patents

基于医学影像数据的几何模型建立方法 Download PDF

Info

Publication number
CN111803803A
CN111803803A CN202010467815.6A CN202010467815A CN111803803A CN 111803803 A CN111803803 A CN 111803803A CN 202010467815 A CN202010467815 A CN 202010467815A CN 111803803 A CN111803803 A CN 111803803A
Authority
CN
China
Prior art keywords
medical image
image data
tissue
geometric model
density
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.)
Granted
Application number
CN202010467815.6A
Other languages
English (en)
Other versions
CN111803803B (zh
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.)
Neuboron Medtech Ltd
Original Assignee
Neuboron Medtech 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 Neuboron Medtech Ltd filed Critical Neuboron Medtech Ltd
Publication of CN111803803A publication Critical patent/CN111803803A/zh
Application granted granted Critical
Publication of CN111803803B publication Critical patent/CN111803803B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30—Polynomial surface description
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/13—Edge detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/168—Segmentation; Edge detection involving transform domain methods
    • 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
    • A61N2005/1085—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
    • A61N2005/109—Neutrons
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30096—Tumor; Lesion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/32—Image data format

Abstract

本发明一方面提供一种基于医学影像数据的几何模型建立方法，包括：读取医学影像数据的步骤；通过医学影像数据与组织种类之间的转换关系定义组织种类的步骤；决定组织分群数的步骤；通过医学影像数据与密度之间的转换关系定义组织密度的步骤；建立带有组织和密度信息的3D编码矩阵的步骤；产生几何模型的步骤。另外，还可包括判断医学影像体素是否在ROI边界内的步骤。根据医学影像的数据与组织种类之间的转换关系，可以根据实际需要确定组织分群数，从而更加精确地提供组织种类、元素组成及密度，建立的几何模型更加匹配于医学影像数据反应出的真实情况。

Description

基于医学影像数据的几何模型建立方法
技术领域
本发明涉及一种几何模型建立方法，尤其涉及一种基于医学影像数据的几何模型建立方法。
背景技术
随着原子科学的发展，例如钴六十、直线加速器、电子射束等放射线治疗已成为癌症治疗的主要手段之一。然而传统光子或电子治疗受到放射线本身物理条件的限制，在杀死肿瘤细胞的同时，也会对射束途径上大量的正常组织造成伤害；另外由于肿瘤细胞对放射线敏感程度的不同，传统放射治疗对于较具抗辐射性的恶性肿瘤(如：多行性胶质母细胞瘤(glioblastoma multiforme)、黑色素细胞瘤(melanoma))的治疗成效往往不佳。
为了减少肿瘤周边正常组织的辐射伤害，化学治疗(chemotherapy)中的标靶治疗概念便被应用于放射线治疗中；而针对高抗辐射性的肿瘤细胞，目前也积极发展具有高相对生物效应(relative biological effectiveness,RBE)的辐射源，如质子治疗、重粒子治疗、中子捕获治疗等。其中，中子捕获治疗便是结合上述两种概念，如硼中子捕获治疗，借由含硼药物在肿瘤细胞的特异性集聚，配合精准的中子射束调控，提供比传统放射线更好的癌症治疗选择。
硼中子捕获治疗(Boron Neutron Capture Therapy,BNCT)是利用含硼(10B)药物对热中子具有高捕获截面的特性，借由10B(n,α)7Li中子捕获及核分裂反应产生4He和7Li两个重荷电粒子，两粒子的总射程约相当于一个细胞大小，因此对于生物体造成的辐射伤害能局限在细胞层级，当含硼药物选择性地聚集在肿瘤细胞中，搭配适当的中子射源，便能在不对正常组织造成太大伤害的前提下，达到局部杀死肿瘤细胞的目的。
三维模型广泛应用于科学实验分析、科学实验模拟领域。比如在核辐射与防护领域，为了模拟人体在一定辐射条件下的吸收剂量，常常需要利用计算机技术对医学影像数据进行各种处理建立精确的MCNP需要的晶格模型，并结合MCNP(蒙特卡罗程序)进行模拟计算。
蒙特卡罗方法是目前能够对辐照目标内部三维空间核粒子碰撞轨迹和能量分布进行精确模拟的工具，蒙特卡罗方法与复杂的三维人体解剖模型相结合代表了模拟在计算机技术中的跃进。在诊断放射检查中，精确的人体器官剂量评估对于放射治疗是非常有益的。目前，国际上已经成功建立多种人体模型并结合蒙特卡罗模拟程序，对人体在辐射环境下的吸收剂量进行精确性的计算评估。人体三维解剖模型成功转换为蒙特卡罗程序所需要的几何描述是进行蒙特卡罗模拟计算的前提条件，也是目前国际上蒙特卡罗模拟研究的热点和难点。
核磁共振成像(Magnetic Resonance Imaging,MRI)或电子计算机断层扫描(Computed Tomography,CT)等医学影像数据能够针对人体体内特征提供较为详细的组织几何结构信息，为人体内部结构的实体建模提供了数据基础。而在中子捕获治疗领域如何根据医学影像数据建立MCNP所需的几何模型属于一个很重要的课题，换句话说，即如何根据医学影像数据建立MCNP软件输入档所需的晶格模型，从而提高治疗计划的精确性。
因此，有必要提出一种提高治疗计划精确性的根据医学影像数据建立MCNP所需的几何模型的方法。
发明内容
本发明的一个方面提供一种基于医学影像数据的几何模型建立方法，包括：
读取医学影像数据的步骤；
通过医学影像数据与组织种类之间的转换关系定义组织种类的步骤；
决定组织分群数的步骤；
通过医学影像数据与密度之间的转换关系定义组织密度的步骤；
医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码的步骤；
整合组织种类和组织密度信息，建立带有组织和密度信息的3D编码矩阵的步骤；
产生几何模型的步骤。
该几何模型建立方法根据医学影像的数据与组织种类之间的转换关系，可以根据实际需要确定组织分群数，从而更加精确地提供组织种类、元素组成及密度，建立的几何模型更加匹配于医学影像数据反应出的真实情况。
作为一种优选地，该几何模型建立方法应用于中子捕获治疗，其进一步包括给定B-10浓度的步骤和建立带有B-10浓度信息的3D编码矩阵的步骤。标记有B-10浓度信息的几何模型，便可清楚地知道，各个组织内的含硼药物浓度，然后进行中子捕获治疗照射模拟时，则更加真实地反应出实际情况。
组织分群数为用户手动定义的组织分群数加上数据库中已有的4种组织分群数或14种组织分群数。如果在已有数据库中并未建立有相对应的组织分群数，那么可以由用户自定义一个新的组织分群数。这样即避免了如果已有数据库中不能完全匹配相对应的组织分群数，只能近似选择的情况，从而有效地提高的建模的精确度。
更加优选地，该几何模型建立方法进一步包括建立3D组织编码矩阵的步骤和建立3D密度编码矩阵的步骤。根据医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码，从而建立起3D组织编码矩阵和3D密度编码矩阵。
几何模型包括MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡。通过医学影像数据最终获得MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡，从而为模拟计算提供理论依据并获得精确的模拟计算结果。
本发明的另一方面提供一种基于医学影像数据的几何模型建立方法，包括：
读取医学影像数据的步骤；
定义或读取ROI边界的步骤；
判断医学影像体素是否在ROI边界内的步骤：
如果是，则进入通过为每个ROI边界内体素给定一个特定组织与密度的方式进行用户手动定义组织种类与密度的步骤或进入通过医学影像数据与组织种类/密度之间的转换关系自动定义ROI组织种类与密度的步骤，
如果否，则进入通过医学影像数据与组织种类之间的转换关系自动定义组织种类的步骤并通过医学影像数据与密度之间的转换关系定义组织密度的步骤；
医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码的步骤；
整合组织种类和组织密度信息，建立带有组织和密度信息的3D编码矩阵的步骤；
产生几何模型的步骤。
所谓ROI是指感兴趣区域(下文统称ROI)，用户可以手动定义ROI的组织种类、元素组成以及密度。如果医学影像体素点不在ROI边界内，则根据医学影像的数据与组织种类之间的转换关系进行组织种类的定义，并根据实际需要确定组织分群数，从而更加精确地提供组织种类、元素组成及密度，建立的几何模型更加匹配于医学影像数据反应出的真实情况。
作为一种优选地，该几何模型建立方法应用于中子捕获治疗，几何模型建立方法包括给定B-10浓度的步骤和建立带有B-10浓度信息的3D编码矩阵的步骤。标记有B-10浓度信息的几何模型，便可清楚地知道，各个组织内的含硼药物浓度，然后进行中子捕获治疗照射模拟时，则更加真实地反应出实际情况。
组织分群数为用户手动定义的组织分群数加上数据库中已有的4种组织分群数或14种组织分群数。如果在已有数据库中并未建立有相对应的组织分群数，那么可以由用户自定义一个新的组织分群数。这样即避免了如果已有数据库中不能完全匹配相对应的组织分群数，只能近似选择的情况，从而有效地提高的建模的精确度。
更加优选地，该几何模型建立方法进一步包括建立3D组织编码矩阵的步骤和建立3D密度编码矩阵的步骤。根据医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码，从而建立起3D组织编码矩阵和3D密度编码矩阵。
几何模型包括MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡。通过医学影像数据最终获得MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡，从而为模拟计算提供理论依据并获得精确的模拟计算结果。
医学影像数据可以为核磁共振成像(Magnetic Resonance Imaging,MRI)或电子计算机断层扫描(Computed Tomography,CT)，下文实施例中将基于电子计算机断层扫描(CT)的数据来阐述，CT的文件格式通常为DICOM。但本领域技术人员熟知地，还可以使用其他的医学影像数据，只要该医学影像数据能够被转换成带有组织种类、密度及B-10浓度信息的MCNP晶格模型的，就能够应用于本发明揭示的基于医学影像数据的几何模型建立方法中。
本发明实施例中的有益效果和/或特点如下：
1、手动定义感兴趣区域(ROI)的组织种类、元素组成及密度；
2、针对非ROI的CT影像体素点，可自动进行组织种类匹配，依据CT值大小差异，已有数据库可区分出4种或14种不同元素组成的组织，也可以根据实际的实验结果确定其他数别的不同元素组成的组织；
3、对于CT值涵盖较广范围的ROI，如黏膜腔室，难以给定唯一的组织密度，用户可根据本发明实施例揭示的方法自动进行CT值与组织/密度的转换；
4、经输入正常血液含硼药物浓度、组织/肿瘤-血液硼浓度比值等参数后，本发明实施例揭示的方法自动将B-10元素编写入所有体素点中；
5、最终产生的三维MCNP晶格模型将带有组织种类(元素组成)、密度和B-10浓度等信息。
附图说明
图1是硼中子捕获反应示意图。
图2是10B(n,α)7Li中子捕获核反应方程式。
图3是本发明实施例中的基于医学影像数据的几何模型建立方法的逻辑框图。
图4是CT值(HU)与组织密度回归曲线公式及相关系数图。
具体实施方式
下面结合附图对本发明的实施例做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
作为一种优选地，以用于中子捕获治疗的基于医学影像数据的几何模型建立方法为本发明的实施例。下面将简单介绍一下中子捕获治疗，尤其是硼中子捕获治疗。
中子捕获治疗作为一种有效的治疗癌症的手段近年来的应用逐渐增加，其中以硼中子捕获治疗最为常见，供应硼中子捕获治疗的中子可以由核反应堆或加速器供应。本发明的实施例以加速器硼中子捕获治疗为例，加速器硼中子捕获治疗的基本组件通常包括用于对带电粒子(如质子、氘核等)进行加速的加速器、靶材与热移除系统和射束整形体，其中加速带电粒子与金属靶材作用产生中子，依据所需的中子产率与能量、可提供的加速带电粒子能量与电流大小、金属靶材的物化性等特性来挑选合适的核反应，常被讨论的核反应有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，两粒子的总射程约相当于一个细胞大小，因此对于生物体造成的辐射伤害能局限在细胞层级，当含硼药物选择性地聚集在肿瘤细胞中，搭配适当的中子射源，便能在不对正常组织造成太大伤害的前提下，达到局部杀死肿瘤细胞的目的。
本发明实施例的目的在于将医学影像数据转换为带有组织种类、密度及B-10浓度信息的MCNP晶格模型，以进行硼中子捕获治疗的组织剂量模拟计算。医学影像数据可以为核磁共振成像(Magnetic Resonance Imaging,MRI)或电子计算机断层扫描(ComputedTomography,CT)，下文实施例中将基于电子计算机断层扫描(CT)的数据来阐述，CT的文件格式通常为DICOM。但本领域技术人员熟知地，还可以使用其他的医学影像数据，只要该医学影像数据能够被转换成带有组织种类、密度及B-10浓度信息的MCNP晶格模型的，就能够应用于本发明揭示的基于医学影像数据的几何模型建立方法中。
简单来说，本发明实施例揭示的基于医学影像数据的几何模型建立方法主要包括如下流程：
1、输入计算机断层影像(DICOM格式)，CT影像将呈现于对应接口上；
2、自动读取DICOM档案中已定义的ROI边界，亦可另外新增ROI；
3、软件将对CT影像上每一体素点(voxel)的位置进行判别，将所有体素点分为「落于ROI边界内」及「落于ROI边界外」；
4、针对ROI边界内的体素点，可根据实际需求，由使用者手动对每一ROI定义相对的唯一组织种类及组织密度，或自动进行CT值与组织种类、密度的匹配，以避免因给定CT值涵盖范围较大的ROI(如黏膜腔室)唯一组织种类(元素组成)及密度，而造成剂量计算的误差；
5、而针对ROI边界外的体素点，则自动进行组织种类的定义，依据CT值大小差异，可区分为4或14种不同元素组成的组织，用户可根据实际判断，选择使用ICRU-46报告表列的4种组织(将在下文详述)，或使用Vanderstraeten等人于2007年发表文献中的14种不同元素组成的组织(将在下文详述)；
6、对于未手动定义密度的体素点，依据CT值大小差异，将自动给定密度，共可区分出96个密度分群；
7、使用者手动输入正常血液含硼药物浓度、肿瘤-血液硼浓度比值、组织-血液硼浓度比值等参数，将编写B-10元素于所有体素点中；
8、软件将整合组织种类(元素组成)、组织密度、组织B-10浓度等信息，产生三维MCNP晶格模型，并编写出MCNP输入档规定格式的晶格卡(lattice card)、栅元卡(cellcard)、曲面卡(surface card)及材料卡(material card)。
具体来说，请参照图3，本发明揭示的基于医学影像数据的几何模型建立方法包括两个实施例。
本发明的第一实施例提供一种基于医学影像数据的几何模型建立方法，包括：读取医学影像数据的步骤；通过医学影像数据与组织种类之间的转换关系定义组织种类的步骤；决定组织分群数的步骤；通过医学影像数据与密度之间的转换关系定义组织密度的步骤；建立带有组织和密度信息的3D编码矩阵的步骤；产生几何模型的步骤。
根据医学影像的数据与组织种类之间的转换关系，可以根据实际需要确定组织分群数，从而更加精确地提供组织种类、元素组成及密度，建立的几何模型更加匹配于医学影像数据反应出的真实情况。
几何模型建立方法进一步包括给定B-10浓度的步骤和建立带有B-10浓度信息的3D编码矩阵的步骤。标记有B-10浓度信息的几何模型，便可清楚地知道，各个体素点内的含硼药物浓度，然后进行中子捕获治疗照射模拟时，则更加真实地反应出实际情况。
组织分群数为用户手动定义的组织分群数加上数据库中已有的4种组织分群数或14种组织分群数。如果在已有数据库中并未建立有相对应的组织分群数，那么可以根据实验情况确定一个新的组织分群数。这样即避免了如果已有数据库中不能完全匹配相对应的组织分群数，只能近似选择的情况，从而有效地提高的建模的精确度。
几何模型建立方法进一步包括建立3D组织编码矩阵的步骤和建立3D密度编码矩阵的步骤。根据医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码，从而建立起3D组织编码矩阵和3D密度编码矩阵。
几何模型包括MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡。通过医学影像数据最终获得MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡，从而为模拟计算提供理论依据并获得精确的模拟计算结果。
本发明的第二实施例提供一种基于医学影像数据的几何模型建立方法，包括：读取医学影像数据的步骤；定义或读取ROI边界的步骤；判断医学影像体素是否在ROI边界内的步骤：如果是，则进入通过为每个ROI边界内体素给定一个特定组织与密度的方式进行用户手动定义组织种类与密度的步骤或或进入通过医学影像数据与组织种类/密度之间的转换关系自动定义ROI组织种类与密度的步骤，如果否，则进入通过医学影像数据与组织种类之间的转换关系自动定义组织种类的步骤并通过医学影像数据与密度之间的转换关系定义组织密度的步骤；建立带有组织和密度信息的3D编码矩阵的步骤；产生几何模型的步骤。
所谓ROI是指感兴趣区域，用户可以手动定义ROI的组织种类、元素组成以及密度。如果不在ROI边界内，根据医学影像的数据与组织种类/密度之间的转换关系进行组织种类的定义，并根据实际需要确定组织分群数，从而更加精确地提供组织种类、元素组成及密度，建立的几何模型更加匹配于医学影像数据反应出的真实情况。
几何模型建立方法包括给定B-10浓度的步骤和建立带有B-10浓度信息的3D编码矩阵的步骤。标记有B-10浓度信息的几何模型，便可清楚地知道，各个体素点内的含硼药物浓度，然后进行中子捕获治疗照射模拟时，则更加真实地反应出实际情况。
组织分群数为用户手动定义的组织分群数加上数据库中已有的4种组织分群数或14种组织分群数。如果单一ROI边界内涵盖CT值范围较大(如黏膜腔室)，则可自动进行CT值与组织种类、密度的匹配，以避免因给定唯一组织种类(元素组成)及密度，而造成剂量计算的误差。如果在已有数据库中并未建立有相对应的组织分群数，那么可以根据实验情况确定一个新的组织分群数。这样即避免了如果已有数据库中不能完全匹配相对应的组织分群数，只能近似选择的情况，从而有效地提高的建模的精确度。
几何模型建立方法进一步包括建立3D组织编码矩阵的步骤和建立3D密度编码矩阵的步骤。根据医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码，从而建立起3D组织编码矩阵和3D密度编码矩阵。
几何模型包括MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡。通过医学影像数据最终获得MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡，从而为模拟计算提供理论依据并获得精确的模拟计算结果。
请参照图4，下面将阐述本发明揭示的基于医学影像数据的几何模型建立方法中已有数据库中的CT值-组织种类及CT值-组织密度转换关系图表。
CT值，又称为Hounsfield Unit(HU)，为反应光衰减系数之单位，其定义如公式一：
Figure BDA0002513265400000081
引用Vanderstraeten等人于2007年发表之文献(Barbara Vanderstraeten etal,“Convension of CT numbers into tissue parametersfor Monte Carlo dosecalculations:a multi-centre study”,Phys.Med.Biol.52(2007)539–562.)，依据CT值大小不同，可转换出1种空气(air)、1种肺组织(lung)、2种软组织(脂肪(adipose)及肌肉软组织(muscle))、10种骨组织(bone)，换言之，不同CT值共对应出14种不同元素组成的组织，如下表一。
表一.CT值对应14种不同元素组成的组织(元素重量百分比)
Figure BDA0002513265400000082
Figure BDA0002513265400000091
引用ICRU-46号报告(International Commission on Radiation Units andMeasurements,Photon,electron,proton and neutron interaction data for bodytissues,ICRU-46,Tech.Rep.,1992.)，撷取4种人体脑部主要组织，包含空气(air)、脑组织(adult brain)、皮肤(adult skin)、头骨(cranium)，其对应之密度与元素组成如表二。
表二.CT值对应4种不同元素组成的组织(元素重量百分比)
Figure BDA0002513265400000092
同样引用Vanderstraeten等人的文献，该份文献统整了医院真实的实验值，整理出CT值对应组织密度的关系公式，如图4所示；本发明揭示的基于医学影像数据的几何模型建立方法使用图4的三组回归公式，将CT值(-1000～2000)区分为96组密度分群，如表三。
表三.CT值与质量密度的转换.HUL和HUU分别为CT值的下限和上限
Figure BDA0002513265400000093
Figure BDA0002513265400000101
本发明实施例中的有益效果和/或特点如下：
1、手动定义感兴趣区域(ROI)的组织种类、元素组成及密度；
2、针对非ROI的CT影像体素点，可自动进行组织种类匹配，依据CT值大小差异，已有数据库可区分出4种或14种不同元素组成的组织，也可以根据实际的实验结果确定其他数别的不同元素组成的组织；
3、对于CT值涵盖较广范围的ROI，如黏膜腔室，难以给定唯一的组织元素组成与密度，用户可根据本发明实施例揭示的方法自动进行CT值与组织/密度的转换；
4、经输入正常血液含硼药物浓度、组织/肿瘤-血液硼浓度比值等参数后，本发明实施例揭示的方法自动将B-10元素编写入所有体素点中；
5、最终产生的三维MCNP晶格模型将带有组织种类(元素组成)、密度和B-10浓度等信息。
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，都在本发明要求保护的范围之内。

Claims (5)

1.一种基于医学影像数据的几何模型建立方法，包括：
读取医学影像数据的步骤，所述医学影像数据通过采集人体的医学影像获得；
通过医学影像数据与组织种类之间的转换关系定义组织种类的步骤；
决定组织分群数的步骤；
通过医学影像数据与密度之间的转换关系定义组织密度的步骤；
医学影像数据的切片通过相对应的转换关系，每一张切片建立相应的组织编码和密度编码的步骤；
整合组织种类和组织密度信息，建立带有组织和密度信息的3D编码矩阵的步骤；
产生几何模型的步骤，所述几何模型用于中子捕获治疗系统的组织剂量模拟计算，所述中子捕获治疗系统包括对带电粒子进行加速的加速器、靶材与热移除系统和射束整形体，所述加速器加速带电粒子与所述靶材作用产生中子。
2.根据权利要求1所述的基于医学影像数据的几何模型建立方法，其特征在于：所述几何模型建立方法应用于中子捕获治疗，所述几何模型建立方法进一步包括给定B-10浓度的步骤和建立带有B-10浓度信息的3D编码矩阵的步骤。
3.根据权利要求1所述的基于医学影像数据的几何模型建立方法，其特征在于：所述组织分群数为用户手动定义的组织分群数加上数据库中已有的4种组织分群数或14种组织分群数。
4.根据权利要求3所述的基于医学影像数据的几何模型建立方法，其特征在于：所述几何模型建立方法进一步包括建立3D组织编码矩阵的步骤和建立3D密度编码矩阵的步骤。
5.根据权利要求1所述的基于医学影像数据的几何模型建立方法，其特征在于：所述几何模型包括MCNP软件输入档所需的晶格卡、栅元卡、曲面卡和材料卡。
CN202010467815.6A 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法 Active CN111803803B (zh)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN2015107902487 2015-11-17
CN201510790248.7A CN106474634A (zh) 2015-11-17 2015-11-17 基于医学影像数据的几何模型建立方法
CN201680017272.XA CN107427692B (zh) 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法
PCT/CN2016/102335 WO2017084459A1 (zh) 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201680017272.XA Division CN107427692B (zh) 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法

Publications (2)

Publication Number Publication Date
CN111803803A true CN111803803A (zh) 2020-10-23
CN111803803B CN111803803B (zh) 2023-02-10

Family

ID=58238342

Family Applications (4)

Application Number Title Priority Date Filing Date
CN201510790248.7A Pending CN106474634A (zh) 2015-11-17 2015-11-17 基于医学影像数据的几何模型建立方法
CN201810305219.0A Pending CN108310683A (zh) 2015-11-17 2015-11-17 基于医学影像数据的几何模型建立方法
CN201680017272.XA Active CN107427692B (zh) 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法
CN202010467815.6A Active CN111803803B (zh) 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法

Family Applications Before (3)

Application Number Title Priority Date Filing Date
CN201510790248.7A Pending CN106474634A (zh) 2015-11-17 2015-11-17 基于医学影像数据的几何模型建立方法
CN201810305219.0A Pending CN108310683A (zh) 2015-11-17 2015-11-17 基于医学影像数据的几何模型建立方法
CN201680017272.XA Active CN107427692B (zh) 2015-11-17 2016-10-18 基于医学影像数据的几何模型建立方法

Country Status (6)

Country Link
US (1) US10692283B2 (zh)
EP (1) EP3357537B1 (zh)
JP (1) JP6754841B2 (zh)
CN (4) CN106474634A (zh)
TW (1) TWI639134B (zh)
WO (1) WO2017084459A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016179381A1 (en) 2015-05-06 2016-11-10 Neutron Therapeuutics Inc. Neutron target for boron neutron capture therapy
CN108295384B (zh) * 2017-01-11 2020-02-28 南京中硼联康医疗科技有限公司 基于医学影像的组织元素质量比例解构方法及几何模型建立方法
US10462893B2 (en) 2017-06-05 2019-10-29 Neutron Therapeutics, Inc. Method and system for surface modification of substrate for ion beam target
CN107290774A (zh) * 2017-07-21 2017-10-24 四川瑶天纳米科技有限责任公司 中子剂量快速确定方法
CN109308733A (zh) * 2017-07-27 2019-02-05 南京中硼联康医疗科技有限公司 基于医学影像数据的几何模型建立方法及剂量计算方法
JP7175964B2 (ja) * 2017-08-24 2022-11-21 南京中硼▲聯▼康医▲療▼科技有限公司 中性子捕捉療法システム
CN110335281A (zh) * 2018-03-28 2019-10-15 北京连心医疗科技有限公司 一种肿瘤边界确定方法、设备和存储介质
CN108852511B (zh) * 2018-04-26 2019-06-07 北京罗森博特科技有限公司 骨内植入物最优位置的自动规划方法及装置
CN109344459B (zh) * 2018-09-13 2022-11-29 北京应用物理与计算数学研究所 一种蒙特卡罗程序的可视建模与转换方法及系统
CN110013613A (zh) * 2019-04-16 2019-07-16 东莞东阳光高能医疗设备有限公司 一种硼中子俘获治疗计划系统
CN111128317A (zh) * 2019-11-20 2020-05-08 中国辐射防护研究院 一种电离辐射组织等效材料配方设计方法及系统
CN113797447A (zh) 2020-06-11 2021-12-17 中硼(厦门)医疗器械有限公司 放射治疗系统及其治疗计划生成方法
CN114073827B (zh) 2020-08-15 2023-08-04 中硼(厦门)医疗器械有限公司 放射线照射系统及其控制方法
CN114367061B (zh) * 2020-10-14 2023-06-23 中硼(厦门)医疗器械有限公司 硼中子捕获治疗系统及其治疗计划生成方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5341292A (en) * 1992-06-04 1994-08-23 New England Medical Center Hospitals, Inc. Monte Carlo based treatment planning for neutron capture therapy
WO2001045556A1 (en) * 1999-12-21 2001-06-28 Bechtel Bwxt Idaho, Llc Monte carlo simulation of neutron transport for use in radiotherapy
US6285969B1 (en) * 1998-05-22 2001-09-04 The Regents Of The University Of California Use of single scatter electron monte carlo transport for medical radiation sciences
EP1658878A1 (en) * 2004-11-17 2006-05-24 The European Community, represented by the European Commission BNCT treatment planning
US20060160157A1 (en) * 2005-01-19 2006-07-20 Zuckerman Mathew M Method, compositions and classification for tumor diagnostics and treatment
CN103549953A (zh) * 2013-10-25 2014-02-05 天津大学 一种基于医学核磁共振影像提取微波检测乳房模型的方法
CN103985122A (zh) * 2014-05-17 2014-08-13 清华大学深圳研究生院 基于心脏ct图像的全心脏提取方法
CN104224277A (zh) * 2013-06-18 2014-12-24 三星电子株式会社 用于产生超声的方法、设备和系统
CN104267425A (zh) * 2014-10-16 2015-01-07 中国科学院合肥物质科学研究院 一种基于CT数据的内照射HPGe探测器探测效率确定方法
US20150131883A1 (en) * 2013-10-25 2015-05-14 The Johns Hopkins University Joint estimation of tissue types and linear attenuation coefficients for computed tomography
CN104899891A (zh) * 2015-06-24 2015-09-09 重庆金山科技（集团）有限公司 一种识别孕囊组织的方法、装置及宫腔吸引装置
WO2015144540A1 (en) * 2014-03-28 2015-10-01 Koninklijke Philips N.V. Method and device for generating one or more computer tomography images based on magnetic resonance images with the help of tissue class separation

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5873830A (en) * 1997-08-22 1999-02-23 Acuson Corporation Ultrasound imaging system and method for improving resolution and operation
CA2491703A1 (en) * 2002-07-09 2004-01-15 Medispectra, Inc. Methods and apparatus for characterization of tissue samples
JP5171251B2 (ja) * 2004-05-04 2013-03-27 スティフテルセン ウニヴァズィテッフォースクニング ベルゲン 灌流磁気共鳴画像化装置の作動方法、灌流磁気共鳴画像化データ操作方法及びこれらの方法に用いるデータ処理装置、コンピュータソフトウェア製品
US8125813B2 (en) * 2005-06-16 2012-02-28 Best Medical International, Inc. Variance reduction simulation system, program product, and related methods
CN101458826B (zh) * 2008-11-25 2012-05-30 中国科学院等离子体物理研究所 利用ct值赋予密度、组成成分的数字人体建模方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5341292A (en) * 1992-06-04 1994-08-23 New England Medical Center Hospitals, Inc. Monte Carlo based treatment planning for neutron capture therapy
US6285969B1 (en) * 1998-05-22 2001-09-04 The Regents Of The University Of California Use of single scatter electron monte carlo transport for medical radiation sciences
WO2001045556A1 (en) * 1999-12-21 2001-06-28 Bechtel Bwxt Idaho, Llc Monte carlo simulation of neutron transport for use in radiotherapy
EP1658878A1 (en) * 2004-11-17 2006-05-24 The European Community, represented by the European Commission BNCT treatment planning
US20060160157A1 (en) * 2005-01-19 2006-07-20 Zuckerman Mathew M Method, compositions and classification for tumor diagnostics and treatment
CN104224277A (zh) * 2013-06-18 2014-12-24 三星电子株式会社 用于产生超声的方法、设备和系统
CN103549953A (zh) * 2013-10-25 2014-02-05 天津大学 一种基于医学核磁共振影像提取微波检测乳房模型的方法
US20150131883A1 (en) * 2013-10-25 2015-05-14 The Johns Hopkins University Joint estimation of tissue types and linear attenuation coefficients for computed tomography
WO2015144540A1 (en) * 2014-03-28 2015-10-01 Koninklijke Philips N.V. Method and device for generating one or more computer tomography images based on magnetic resonance images with the help of tissue class separation
CN103985122A (zh) * 2014-05-17 2014-08-13 清华大学深圳研究生院 基于心脏ct图像的全心脏提取方法
CN104267425A (zh) * 2014-10-16 2015-01-07 中国科学院合肥物质科学研究院 一种基于CT数据的内照射HPGe探测器探测效率确定方法
CN104899891A (zh) * 2015-06-24 2015-09-09 重庆金山科技（集团）有限公司 一种识别孕囊组织的方法、装置及宫腔吸引装置

Also Published As

Publication number Publication date
WO2017084459A1 (zh) 2017-05-26
EP3357537B1 (en) 2020-09-30
US10692283B2 (en) 2020-06-23
TW201719580A (zh) 2017-06-01
CN107427692B (zh) 2020-08-07
CN111803803B (zh) 2023-02-10
CN108310683A (zh) 2018-07-24
EP3357537A4 (en) 2018-10-24
EP3357537A1 (en) 2018-08-08
TWI639134B (zh) 2018-10-21
US20180247452A1 (en) 2018-08-30
CN106474634A (zh) 2017-03-08
JP2019502505A (ja) 2019-01-31
JP6754841B2 (ja) 2020-09-16
CN107427692A (zh) 2017-12-01

Similar Documents

Publication Publication Date Title
CN107427692B (zh) 基于医学影像数据的几何模型建立方法
US10643761B2 (en) Method for evaluating irradiation angle of beam
CN108295384B (zh) 基于医学影像的组织元素质量比例解构方法及几何模型建立方法
US11087524B2 (en) Method for establishing smooth geometric model based on data of medical image
CN116603176A (zh) 硼中子捕获治疗系统及其照射参数选取方法
EP1658878A1 (en) BNCT treatment planning
CN107292075A (zh) 增进放射治疗系统计算效益的方法
TWI791390B (zh) 中子捕獲治療設備及照射參數選取裝置的使用方法
JP2023531798A (ja) 放射線療法システム及びその治療計画生成方法
Deng et al. The dosimetry calculation for boron neutron capture therapy
CN117982812A (zh) 硼中子捕获治疗系统及其治疗计划生成方法

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