本发明涉及计算机辅助诊断领域，尤其是一种能谱计算机断层扫描剂量的蒙特卡罗模拟方法及系统。
背景技术：
：名词解释：ct：计算机断层扫描。能谱ct成像：根据物质在不同x射线能量下产生的不同的吸收原理，采用两种或以上x射线管电压，提供比常规ct更多的影像信息。ctdi：ct剂量指数，用于反映ct扫描的剂量水平。ctdi100：ct剂量指数100，在剂量平面前后各50mm范围内的剂量曲线的积分，具体计算公式为：其中d(z)为z方向(即ct轴线方向)上的剂量曲线，n为扫描层数，t为层厚。ctdiw：加权ct剂量指数，反应扫描平面中的平均剂量，与扫描平面中心的ct剂量指数100ctdi100(中心)和扫描平面中心外周的ct剂量指数100ctdi100(外周)有关，具体计算公式为：ctdivol：容积ct剂量指数，用于反映整个螺旋扫描容积中的层面平均剂量，ctdivol＝ctdiw/螺距因子。ctdiair：扫描腔内无其它物体，仅有空气和扫描床时测得的扫描平面中心处的ct剂量指数100；dlp：评价受检者接受一次ct曝光扫描后接受的总的辐射剂量，dlp＝ctdivol×l，ctdivol为容积ct剂量指数，l为沿z轴的扫描长度。tld：热释光剂量计，用于直接测量剂量。tasmip：tungstenanodespectralmodelusinginterpolatingpolynomials，利用插值多项式的钨极x射线发射能谱模型，用于模拟钨极x射线能谱。bowtie滤线器：位于ct的x射线源(球管)之后的金属块，一般用铝制成，用于限制x射线的角分布；hvl：halfvaluelayer，半值层厚，是指将x射线的比释动能降低至初始值的一半所需的铝片等金属片厚度；qvl：quartervaluelayer，四分之一值层厚，是指将x射线的比释动能降低至初始值的四分之一所需的铝片等金属片厚度。管电压：ct的x射线球管电压，决定了出射x射线的能量。管电流：ct的x射线球管电流，决定了出射x射线的强度。dicom：digitalimagingandcommunicationsinmedicine，即医学数字成像和通信，是医学图像和相关信息的国际标准(iso12052)，它定义了质量能满足临床需要的可用于数据交换的医学图像格式。sfov：扫描视野。束宽：出射x射线束在z方向的宽度。扫描长度：扫描床相对机架移动的距离。现有的能谱ct扫描剂量计算方法主要包括ctdi体模剂量估算法和仿真人体体模直接测量法两种，前者采用了ctdi体模，通过测量体模中的ctdivol，最终计算得到dlp，由于ctdi体模与人体差别较大，剂量计算存在很大误差；后者则在较复杂的仿真人体体模内放置tld进行直接测量，其需要放置多个tld或进行多次测量，工作量大，耗时长，不能迅速给出给定扫描方案的剂量，尤其是无法给出具体受扫描者的剂量(不可能在人体内放置tld)。因此，业内亟需一种同时满足高计算精度和高计算效率需求的能谱ct辐射剂量计算方案。技术实现要素：为解决上述技术问题，本发明的目的在于：提供一种高精度和高效率的能谱计算机断层扫描剂量的蒙特卡罗模拟方法及系统。本发明所采取的第一技术方案是：一种能谱计算机断层扫描剂量的蒙特卡罗模拟方法，包括以下步骤：根据实测ct机的衰减层厚数据结合虚拟软射线源和虚拟硬化材料对ct机的x射线能谱进行模拟，得到模拟能谱；根据实测ct机的角分布数据和模拟能谱对虚拟滤线器厚度进行模拟，得到虚拟滤线器的形状；根据ct机的仿真体模扫描图像进行三维建模，得到三维仿真体模几何模型；根据模拟能谱、虚拟滤线器的形状和三维仿真体模几何模型进行蒙特卡罗模拟，得到仿真体模各器官的吸收剂量；根据仿真体模各器官的吸收剂量、给定的各管电压和给定的时间占比，计算出能谱ct扫描的吸收剂量。进一步，所述虚拟软射线源为利用tasmip能谱模型生成的虚拟钨极x射线源，所述虚拟硬化材料为铜片、铝片或其它具有x射线能谱硬化效果的材料，所述虚拟滤线器为bowtie滤线器。进一步，所述根据实测ct机的衰减层厚数据结合虚拟软射线源和虚拟硬化材料对ct机的x射线能谱进行模拟，得到模拟能谱这一步骤，具体包括：采用实验测量的方法，获得ct机扫描腔中心线处的半值层厚和/或四分之一值层厚；通过虚拟钨极x射线源生成初始能谱；将初始能谱穿过虚拟硬化材料后透过bowtie滤线器；计算透过bowtie滤线器中心后的半值层厚和/或四分之一值层厚；调整虚拟硬化材料的厚度，使得透过bowtie滤线器后的半值层厚和/或四分之一值层厚与ct机扫描腔中心线处的半值层厚和/或四分之一值层厚相匹配；选定透过厚度调整后的虚拟硬化材料后的能谱作为模拟能谱。进一步，所述根据实测ct机的角分布数据和模拟能谱对虚拟滤线器厚度进行模拟，得到虚拟滤线器的形状这一步骤，具体包括：采用实验测量的方法，得到实测ct机扫描腔垂向中心面上各个高度位置的空气吸收剂量；根据给定的bowtie滤线器中心的厚度计算透过bowtie滤线器中心后的空气吸收剂量，进而得到计算的空气吸收剂量值与实测的空气吸收剂量值的比值，并将得到的比值记为r；选定x射线的出射角序列；对出射角序列中的每个出射角度根据模拟能谱和r计算bowtie滤线器在每个出射角度对应的厚度，最终得到与出射角度序列对应的bowtie滤线器厚度序列，其中，bowtie滤线器在每个出射角度对应的厚度使得透过bowtie滤线器后计算得到的空气吸收剂量与该出射角度实测的空气吸收剂量比值等于r。进一步，所述根据ct机的仿真体模扫描图像进行三维建模，得到三维仿真体模几何模型这一步骤，具体包括：通过ct机扫描仿真体模，获得仿真体模的dicom图像；对仿真体模的dicom图像进行分割、像素合并以及器官标记；根据分割、像素合并以及器官标记的结果建立三维矩阵形式的仿真体模几何模型。进一步，所述根据模拟能谱、虚拟滤线器的形状和三维仿真体模几何模型进行蒙特卡罗模拟，得到仿真体模各器官的吸收剂量这一步骤，具体包括：根据三维仿真体模几何模型生成蒙特卡罗输运程序输入文件的几何描述；根据ct机的扫描参数、模拟能谱和虚拟滤线器的形状进行源定义，得到源粒子的起始位置、源粒子的出射角度、源粒子的能量以及出射粒子的权重；根据源定义的结果，采用蒙特卡罗输运程序模拟计算单圈轴向扫描条件下扫描腔中心线上的空气吸收剂量；根据源定义的结果，采用蒙特卡罗输运程序模拟计算螺旋扫描条件下三维仿真体模中各器官的吸收剂量。进一步，所述根据仿真体模各器官的吸收剂量、给定的各管电压和给定的时间占比，计算出能谱ct扫描的吸收剂量这一步骤，具体包括：采用与蒙特卡罗输运程序模拟相同的参数实测ct机单圈轴向扫描条件下扫描中心线上的空气吸收剂量；根据单圈轴向扫描条件下实测的空气吸收剂量和蒙特卡罗输运程序模拟计算的空气吸收剂量计算得到归一化因子；根据归一化因子和实际的扫描参数计算器官的实际吸收剂量；根据给定的时间占比，对给定的各管电压下器官的实际吸收剂量进行加权平均，最终得到能谱ct扫描的吸收剂量。本发明所采取的第二技术方案是：一种能谱计算机断层扫描剂量的蒙特卡罗模拟系统，包括：模拟能谱生成模块，用于根据实测ct机的衰减层厚数据结合虚拟软射线源和虚拟硬化材料对ct机的x射线能谱进行模拟，得到模拟能谱；虚拟滤线器形状确定模块，用于根据实测ct机的角分布数据和模拟能谱对虚拟滤线器厚度进行模拟，得到虚拟滤线器的形状；三维几何模型建模模块，用于根据ct机的仿真体模扫描图像进行三维建模，得到三维仿真体模几何模型；蒙特卡罗模拟模块，用于根据模拟能谱、虚拟滤线器的形状和三维仿真体模几何模型进行蒙特卡罗模拟，得到仿真体模各器官的吸收剂量；能谱ct扫描吸收剂量计算模块，用于根据仿真体模各器官的吸收剂量、给定的各管电压和给定的时间占比，计算出能谱ct扫描的吸收剂量。进一步，所述蒙特卡罗模拟模块包括：几何描述生成单元，用于根据三维仿真体模几何模型生成蒙特卡罗输运程序输入文件的几何描述；源定义单元，用于根据ct机的扫描参数、模拟能谱和虚拟滤线器的形状进行源定义，得到源粒子的起始位置、源粒子的出射角度、源粒子的能量以及出射粒子的权重；单圈轴向模拟单元，用于根据源定义的结果，采用蒙特卡罗输运程序模拟计算单圈轴向扫描条件下扫描腔中心线上的空气吸收剂量；器官吸收剂量计算单元，用于根据源定义的结果，采用蒙特卡罗输运程序模拟计算螺旋扫描条件下三维仿真体模中各器官的吸收剂量。本发明所采取的第三技术方案是：一种能谱计算机断层扫描剂量的蒙特卡罗模拟系统，包括：存储器，用于存储程序；处理器，用于加载所述程序以执行如第一技术方案所述的能谱计算机断层扫描剂量的蒙特卡罗模拟方法。本发明的有益效果是：本发明一种能谱计算机断层扫描剂量的蒙特卡罗模拟方法，根据ct机实测的衰减层厚数据和扫描腔的角分布数据，采用软件模拟的方法得到模拟能谱和虚拟滤线器的形状，再根据ct机的仿真体模扫描图像进行三维建模，然后通过蒙特卡罗模拟得到仿真体模各器官的吸收剂量，最后根据蒙特卡罗模拟的吸收剂量计算出能谱ct扫描的吸收剂量，通过软件模拟和蒙特卡罗模拟能在不进行仿真体模内实际剂量测量的情况下，以较高的精度和较快的速度实现能谱ct扫描剂量的模拟计算。附图说明图1为本发明一种能谱计算机断层扫描剂量的蒙特卡罗模拟方法的整体流程图；图2为本发明实施例一的具体流程图；图3为x射线能谱模拟过程示意图；图4为bowtie滤线器的形状测量过程示意图；图5为腹部某一层的dicom图像；图6为图5的mcnp几何模型图。具体实施方式参照图1，本发明一种能谱计算机断层扫描剂量的蒙特卡罗模拟方法，包括以下步骤：根据实测ct机的衰减层厚数据结合虚拟软射线源和虚拟硬化材料对ct机的x射线能谱进行模拟，得到模拟能谱；根据实测ct机的角分布数据和模拟能谱对虚拟滤线器厚度进行模拟，得到虚拟滤线器的形状；根据ct机的仿真体模扫描图像进行三维建模，得到三维仿真体模几何模型；根据模拟能谱、虚拟滤线器的形状和三维仿真体模几何模型进行蒙特卡罗模拟，得到仿真体模各器官的吸收剂量；根据仿真体模各器官的吸收剂量、给定的各管电压和给定的时间占比，计算出能谱ct扫描的吸收剂量。其中，虚拟软射线源、虚拟硬化材料和虚拟滤线器均属于计算机软件模拟的(即虚拟的)，用于对ct机的实际扫描进行模拟。虚拟软射线源发射的能谱相对于真实人体ct扫描的能谱来说较软，所以本发明加入了虚拟硬化材料来调整x射线光子的能量分布，减少低能x光子的占比，使得透过虚拟硬化材料后的能谱最大限度地与真实人体ct扫描的能谱相同或接近，提升了能谱ct扫描剂量计算的精度。ct机的衰减层厚数据主要包括半值层厚hvl和/或四分之一值层厚qvl。进一步作为优选的实施方式，所述虚拟软射线源为利用tasmip能谱模型生成的虚拟钨极x射线源，所述虚拟硬化材料为铜片、铝片或其它具有x射线能谱硬化效果的材料，所述虚拟滤线器为bowtie滤线器。进一步作为优选的实施方式，所述根据实测ct机的衰减层厚数据结合虚拟软射线源和虚拟硬化材料对ct机的x射线能谱进行模拟，得到模拟能谱这一步骤，具体包括：采用实验测量的方法，获得ct机扫描腔中心线处的半值层厚和/或四分之一值层厚；通过虚拟钨极x射线源生成初始能谱；将初始能谱穿过虚拟硬化材料后透过bowtie滤线器；计算透过bowtie滤线器中心后的半值层厚和/或四分之一值层厚；调整虚拟硬化材料的厚度，使得透过bowtie滤线器后的半值层厚和/或四分之一值层厚与ct机扫描腔中心线处的半值层厚和/或四分之一值层厚相匹配；选定透过厚度调整后的虚拟硬化材料后的能谱作为模拟能谱。进一步作为优选的实施方式，所述根据实测ct机的角分布数据和模拟能谱对虚拟滤线器厚度进行模拟，得到虚拟滤线器的形状这一步骤，具体包括：采用实验测量的方法，得到实测ct机扫描腔垂向中心面上各个高度位置的空气吸收剂量；根据给定的bowtie滤线器中心的厚度计算透过bowtie滤线器中心后的空气吸收剂量，进而得到计算的空气吸收剂量值与实测的空气吸收剂量值的比值，并将得到的比值记为r；选定x射线的出射角序列；对出射角序列中的每个出射角度根据模拟能谱和r计算bowtie滤线器在每个出射角度对应的厚度，最终得到与出射角度序列对应的bowtie滤线器厚度序列，其中，bowtie滤线器在每个出射角度对应的厚度使得透过bowtie滤线器后计算得到的空气吸收剂量与该出射角度实测的空气吸收剂量比值等于r。根据医学先验知识，ct机x射线出射角度的不同(即扫描腔的角分布不同)，会导致受扫描者最终的吸收剂量不同，反映在模拟过程中就是虚拟软射线源的出射角度的不同，透过bowtie滤线器的厚度不同。进一步作为优选的实施方式，所述根据ct机的仿真体模扫描图像进行三维建模，得到三维仿真体模几何模型这一步骤，具体包括：通过ct机扫描仿真体模，获得仿真体模的dicom图像；对仿真体模的dicom图像进行分割、像素合并以及器官标记；根据分割、像素合并以及器官标记的结果建立三维矩阵形式的仿真体模几何模型。进一步作为优选的实施方式，所述根据模拟能谱、虚拟滤线器的形状和三维仿真体模几何模型进行蒙特卡罗模拟，得到仿真体模各器官的吸收剂量这一步骤，具体包括：根据三维仿真体模几何模型生成蒙特卡罗输运程序输入文件的几何描述；根据ct机的扫描参数、模拟能谱和虚拟滤线器的形状进行源定义，得到源粒子的起始位置、源粒子的出射角度、源粒子的能量以及出射粒子的权重；根据源定义的结果，采用蒙特卡罗输运程序模拟计算单圈轴向扫描条件下扫描腔中心线上的空气吸收剂量；根据源定义的结果，采用蒙特卡罗输运程序模拟计算螺旋扫描条件下三维仿真体模中各器官的吸收剂量。进一步作为优选的实施方式，所述根据仿真体模各器官的吸收剂量、给定的各管电压和给定的时间占比，计算出能谱ct扫描的吸收剂量这一步骤，具体包括：采用与蒙特卡罗输运程序模拟相同的参数实测ct机单圈轴向扫描条件下扫描中心线上的空气吸收剂量；根据单圈轴向扫描条件下实测的空气吸收剂量和蒙特卡罗输运程序模拟计算的空气吸收剂量计算得到归一化因子；根据归一化因子和实际的扫描参数计算器官的实际吸收剂量；根据给定的时间占比，对给定的各管电压下器官的实际吸收剂量进行加权平均，最终得到能谱ct扫描的吸收剂量。其中，归一化因子等于单圈轴向扫描条件下实测的空气吸收剂量除以蒙特卡罗输运程序模拟计算的空气吸收剂量。与图1的方法相对应，本发明一种能谱计算机断层扫描剂量的蒙特卡罗模拟系统，包括：模拟能谱生成模块，用于根据实测ct机的衰减层厚数据结合虚拟软射线源和虚拟硬化材料对ct机的x射线能谱进行模拟，得到模拟能谱；虚拟滤线器形状确定模块，用于根据实测ct机的角分布数据和模拟能谱对虚拟滤线器厚度进行模拟，得到虚拟滤线器的形状；三维几何模型建模模块，用于根据ct机的仿真体模扫描图像进行三维建模，得到三维仿真体模几何模型；蒙特卡罗模拟模块，用于根据模拟能谱、虚拟滤线器的形状和三维仿真体模几何模型进行蒙特卡罗模拟，得到仿真体模各器官的吸收剂量；能谱ct扫描吸收剂量计算模块，用于根据仿真体模各器官的吸收剂量、给定的各管电压和给定的时间占比，计算出能谱ct扫描的吸收剂量。进一步作为优选的实施方式，所述蒙特卡罗模拟模块包括：几何描述生成单元，用于根据三维仿真体模几何模型生成蒙特卡罗输运程序输入文件的几何描述；源定义单元，用于根据ct机的扫描参数、模拟能谱和虚拟滤线器的形状进行源定义，得到源粒子的起始位置、源粒子的出射角度、源粒子的能量以及出射粒子的权重；单圈轴向模拟单元，用于根据源定义的结果，采用蒙特卡罗输运程序模拟计算单圈轴向扫描条件下扫描腔中心线上的空气吸收剂量；器官吸收剂量计算单元，用于根据源定义的结果，采用蒙特卡罗输运程序模拟计算螺旋扫描条件下三维仿真体模中各器官的吸收剂量。与图1的方法相对应，本发明一种能谱计算机断层扫描剂量的蒙特卡罗模拟系统，包括：存储器，用于存储程序；处理器，用于加载所述程序以执行如本发明所述的能谱计算机断层扫描剂量的蒙特卡罗模拟方法。下面结合说明书附图和具体实施例对本发明作进一步解释和说明。实施例一以软射线源为利用tasmip能谱模型生成的虚拟钨极x射线源(即初始能谱为tasmip能谱)，虚拟滤线器为bowtie滤线器为例，如图2所示，本实施例能谱ct扫描剂量的蒙特卡罗模拟技术方案主要实现步骤包括：一、能谱ct扫描x射线能谱和角分布的模拟为了使用计算机来重现ct扫描过程，必须获得穿过人体的x射线束的能谱和角分布。ct的能谱由管电压决定，对于特定的ct机型和给定的管电压来说，其扫描x射线能谱的模拟过程如下：(1)采用实验测量的方法，获得扫描腔中心线处的hvl和qvl；(2)以较软的tasmip能谱作为初始能谱，并设定bowtie滤线器的中心厚度为0.5mm，在虚拟钨极x射线源和bowtie滤线器之间加入虚拟硬化材料(铜片、铝片或其它具有x射线能谱硬化效果的材料)，以及将透过虚拟硬化材料后的能谱作为模拟ct能谱的候选能谱。(3)以较小的步长缓慢增加虚拟硬化材料的厚度，同时使用给定的质量衰减系数计算候选能谱透过bowtie滤线器后的hvl和qvl。当计算的hvl和qvl与实验测得的hvl和qvl匹配(可采取qvl相等，hvl最接近的原则，或仅使hvl相等)时，透过虚拟硬化材料和bowtie滤线器后的x射线能谱即可用作该管电压下的模拟能谱。上述x射线能谱模拟过程如图3所示，图3中，1为tasmip能谱，2为虚拟硬化材料，3为候选能谱，4为bowtie滤线器(中心的金属块为铝，中心厚度0.5mm)，5为用于计算hvl和qvl的铝片。6代表调整虚拟硬化材料的厚度或种类，使得hvl和/或qvl与实验测得的相近或相同。虚拟硬化材料2和铝片5的厚度可通过分立的虚拟硬化材料或铝片一层一层地叠加来增加。而x射线角分布的模拟则通过调节虚拟bowtie滤线器的形状的方式来实现，具体过程如下：(1)采用实验测量扫描腔垂向中心面上各个高度(或角度)位置的空气吸收剂量。(2)将bowtie滤线器中心的厚度设为0.5mm，使用预设的质量衰减系数和空气的质能衰减系数计算透过bowtie滤线器中心后的空气吸收剂量，将该剂量值与实测扫描腔中心的空气吸收剂量值的比值设为r。(3)选取一个x射线的出射角度序列，针对每个出射角度，采用前面得到的模拟能谱，以较小的步长缓慢增加bowtie滤线器在该角度上的厚度，然后使用预设的质量衰减系数和空气的质能衰减系数计算透过bowtie滤线器后的空气吸收剂量，当该剂量与该出射角度下实测扫描腔中心的空气吸收剂量的比值等于r时，记录此时的bowtie厚度并将其作为该出射角度对应的bowtie厚度。最终可以得到与出射角度序列对应的bowtie厚度序列，即虚拟bowtie滤线器的形状。bowtie滤线器的形状的测量过程如图4所示，图4中，4为bowtie滤线器，7为虚拟钨极x射线源，8为测量点，9为机架。二、仿真体模几何模型的建立仿真体模几何模型的建立具体过程如下：(1)、通过ct机扫描高仿真体模，获得较高分辨率的仿真体模全身的dicom图像；(2)、根据ct值的大小和解剖学判定理论，对dicom图像进行分割，并进行像素的合并，标记像素所属的器官或组织；(3)、根据分割、像素合并和标记结果，形成扫描区域的三维矩阵，即仿真体模几何模型。三维矩阵中每个元素代表该立方体元的器官或组织类型，每个立方体元内的材料成分、密度相同。该仿真体模几何模型用于后续生成蒙特卡罗模拟输入文件中的几何描述。三、ct扫描过程的蒙特卡罗模拟ct扫描过程的蒙特卡罗模拟具体过程如下：(1)、根据仿真体模几何模型生成蒙特卡罗输运程序输入文件中的几何描述。根据蒙特卡罗程序的不同，可采取不同的描述方式，以mcnp程序为例，其采用lattice重复结构描述方法。(2)、进行源定义：根据ct扫描参数，如扫描长度、扫描时间、每圈扫描时间等，随机抽样获得x射线源粒子的起始位置；根据有效视野、束宽，随机抽样获得源粒子的出射角度；根据出射角度及与之对应的bowtie厚度确定出射粒子权重。(3)、根据源定义的结果，使用蒙特卡罗输运程序模拟计算单圈轴向扫描条件下，扫描腔中心线上的空气吸收剂量ctdiair，记为cm。(4)、根据源定义的结果，使用蒙特卡罗输运程序模拟计算螺旋扫描条件下仿真体模几何模型中各器官的吸收剂量。将得到某器官的吸收剂量记为dm。四、能谱ct的吸收剂量计算(1)通过查阅ct机的说明书或直接试验测量得到ct机以与三的(3)相同的参数(管电压、扫描电流为ma0、旋转一圈时间为t0，机架只旋转一圈)进行扫描时ct扫描中心线上的空气吸收剂量ctdiair，记为c0。记nf为归一化因子，则nf＝c0/cm；(2)根据ct机的实际扫描圈数n，实际扫描电流ma，实际旋转一圈时间t，计算得到某器官实际吸收剂量d：(3)假设能谱ct使用的管电压分别为v1和v2，时间占比为t1：t2，通过前一步骤求得的某器官吸收剂量分别为d1和d2，则某器官最终的吸收剂量(即能谱ct扫描的吸收剂量)为(d1*t1+d2*t2)/(t1+t2)。实施例二本实施例对应用本发明的方法进行实际的能谱ct扫描的吸收剂量的原理、过程及结果进行说明。以交替使用80kvp和140kvp管电压的能谱ct扫描为例，对仿真体模进行由头顶至大腿部的螺旋扫描，扫描参数为管电流260ma，sfov＝500mm，束宽40mm，螺矩3.936cm，扫描长度93.5cm，单圈扫描时间0.6s，则其剂量模拟计算过程如下：1、查阅trm手册(ct机使用说明手册)得到某型ct机在80kvp和140kvp管电压，扫描视野为“largebody”时，hvl分别为6.1mmal和8.4mmal；分别生成80kvp和140kvp的tasmip能谱作为初始能谱，假设bowtie滤线器中心厚度为0.5mm(al)，以al为虚拟硬化材料，缓慢增加虚拟硬化材料的厚度，同时不断计算x射线穿过虚拟硬化材料和bowtie滤线器后的hvl，发现当虚拟硬化材料厚度分别为1cm和0.75cm时，80kvp和140kvp的hvl恰好分别为6.1mm和8.4mm，选择此时穿过虚拟硬化材料后的x射线能谱为ct扫描模拟能谱。2、由于尚无80kvp和140kvp管电压下的扫描腔垂向中心面上的空气吸收剂量测量结果，在本实施例中以120kvp管电压下的测量结果为基础对bowtie滤线器的形状进行计算(考虑到实际扫描时，管电压变化较为频繁，但bowtie滤线器一般不发生变化，因此这一替代方法是合理的)。测量的结果如表1所示，而bowtie滤线器的形状计算结果如表2所示。表1表23、通过预先进行的ct扫描获得仿真体模的全身图像，该全身图像涵盖了体模从头顶到大腿部位的区域。图5展示了腹部某一层的dicom图像。4、ct扫描原始图像每层的分辨率为512*512，通过对该图像进行合并(分辨率变为128*128)，并参照ct值和器官位置进行分割后，根据分割结果，利用蒙特卡罗输运程序mcnp5的lattice卡对体模进行了几何建模。建模完成后，体模的每一层将由128*128个尺寸相同的体元构成，每个体元具有单一的材料成分。与图5的dicom图像对应的mcnp计算所用的几何模型如图6所示。5、查阅trm手册得到80kvp管电压下的扫描参数为管电流260ma，sfov＝500mm，束宽40mm，单圈扫描时间为1s时，ctdiair为21.74mgy±40％；140kvp管电压，在同样的扫描参数下，ctdiair为82.71mgy±40％。6、采用mcnp5作为蒙特卡罗光子输运程序，修改其源定义代码，则源位置抽样方法(即源定义的方法)为：在直径为539mm，螺矩为3.936cm(即束宽*pitch)，总长度93.5cm的螺旋线上均匀抽样；出射角度在扇束角范围内均匀抽样，出射光子(即粒子)能量分别根据80kvp和140kvp的模拟能谱抽样；根据出射角度查表得到所穿过的bowtie厚度，根据预设的质量衰减系数计算粒子权重。蒙特卡罗输运程序首先对ctdiair模型进行模拟(即在扫描腔内仅放置一个电离室)，得到计数值分别为1.00923e-05(80kvp)和1.00481e-05(140kvp)；之后对仿真体模模型按给定扫描参数(管电流260ma，sfov＝500mm，束宽40mm，螺矩3.936cm，扫描长度93.5cm)进行光子输运模拟，并对特定器官进行计数，得到各器官的计数值为分别如下表3所示：表3mcnp5计数7、根据时间占比等比例关系，计算得到器官实际吸收剂量值如表4所示。设能谱ct扫描时，80kvp和140kvp时间比例为1：1(如果机型和扫描方式不同，该比例有所不同)，则可求得能谱ct扫描的器官剂量值，如表4所示。表4吸收剂量器官80kvp剂量(mgy)140kvp剂量(mgy)1:1能谱扫描剂量(mgy)肺9.11e+004.00e+0124.55495心脏9.07e+004.02e+0124.65735脑9.14e+004.01e+0124.64221肾上腺8.72e+003.95e+0124.11806膀胱8.76e+003.94e+0124.06839胸7.54e+003.39e+0120.72765结肠1.01e+014.37e+0126.85452食管8.76e+003.83e+0123.52733睾丸/卵巢1.06e+014.29e+0126.72926肾9.04e+004.07e+0124.8619眼晶状体1.17e+014.55e+0128.56769肝8.79e+003.95e+0124.13486胰腺9.21e+004.03e+0124.74995前列腺8.93e+003.83e+0123.63748唾液腺/腮腺1.28e+014.98e+0131.27026脾9.08e+004.10e+0125.04579胃8.93e+003.95e+0124.20083甲状腺1.50e+016.13e+0138.13466综上所述，本发明一种能谱计算机断层扫描剂量的蒙特卡罗模拟方法及系统，首先通过tasmip能谱和虚拟硬化材料相结合来根据实测的hvl和/或qvl得到x射线模拟能谱；然后通过ct机扫描仿真体模得到扫描图像，并进行合并、分割和标记形成三维矩阵形式的用于蒙特卡罗输运模拟的人体几何模型；再根据扫描参数抽样x射线源位置和光子出射角度，根据模拟能谱抽样光子能量，进行蒙特卡罗输运模拟，得到各器官吸收剂量；最后根据时间占比，通过对各管电压下吸收剂量的加权平均，得到能谱ct扫描的吸收剂量。本发明通过软件模拟和蒙特卡罗模拟，无需在较复杂的仿真人体体模内放置tld进行测量，能在不进行仿真体模内实际剂量测量的情况下，以较高的精度和较快的速度实现能谱ct扫描剂量的模拟计算，能应用于能谱ct进行局部或全身扫描时各个器官接收的x射线辐射剂量的估算。以上是对本发明的较佳实施进行了具体说明，但本发明并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这些等同的变形或替换均包含在本申请权利要求所限定的范围内。当前第1页12