本申请涉及医学图像处理
技术领域：
：，特别涉及pet/ct的图像融合绘制方法、装置、设备及存储介质。
背景技术：
：：正电子发射断层显像/x线计算机体层成像仪(positronemissiontomography/computedtomography，pet/ct)是正电子发射断层显像(positronemissiontomography，pet)和x线计算机体层成像仪(computedtomography，ct)组合体，将pet和ct设计为一体，由一个工作站控制。单pet进行核医学显像时，有其它诊断设备无法比拟的早期发现及灵敏性高等优越特性，但因药物及其原理所限，其定位精度不够，将pet和ct设计为一体，扫描时根据需求同时进行pet显像和ct显像，并将两种图像融合在一起，以达到更好的鉴别和定位。但是，目前的pet/ct的图像融合显示的图像是将根据最大密度投影(maximalintensityprojection，mip)算法获取的pet图像与根据体绘制(volumerendering，vr)算法绘制的ct图像直接加权平均融合成pet/ct图像，由于上述pet/ct图像只是简单的图像融合，不具有深度信息，从图像中看不出pet的最大密度位置在图像的具体深度。技术实现要素：有鉴于此，本申请提出pet/ct的图像融合绘制方法、装置、设备及存储介质以解决上述技术问题。为了达到上述目的，本申请所采用的技术方案为：本申请提供pet/ct的图像融合绘制方法、装置、设备及机器可读存储介质。根据本申请实施例的第一方面，提供一种pet/ct的图像绘制方法，所述方法包括：在采样路径上对pet体数据进行采样，获取所述采样路径上的最大密度点坐标及所述最大密度点处的pet颜色值；在所述采样路径上，根据预设步长对ct体数据进行采样，采样至最大密度点坐标之前，根据各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值；在所述最大密度点对应的位置，将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值；继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值；根据所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像。在一些例子中，所述根据所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像，包括：若在所述采样路径上，所述最大密度点在采样结束点之后，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第一目标颜色值确定；或若在所述采样路径上，所述最大密度点为采样结束点，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第二目标颜色值确定；或若在所述采样路径上，所述最大密度点在采样结束点之前，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第三目标颜色值确定。在一些例子中，所述各采样点的ct颜色值及各采样点的阻光度根据所述被采样点的ct体数据通过传递函数转换得到。在一些例子中，所述方法还包括：根据最大密度投影算法或光线投影算法对pet体数据进行采样；根据体绘制算法对ct体数据进行采样。在一些例子中，所述采样路径为光线投射路径，所述光线投射路径穿越pet体数据或ct体数据进行采样。在一些例子中，所述将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值，包括：所述第二目标颜色值为所述pet颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加而得到；其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述pet颜色值确定；当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度根据本申请实施例的第二方面，提供一种pet/ct的图像绘制装置，所述装置包括：pet采集模块，用于在采样路径上对pet体数据进行采样，获取所述采样路径上的最大密度点坐标及所述最大密度点处的pet颜色值；ct采集模块，用于在所述采样路径上，根据预设步长对ct体数据进行采样，采样至最大密度点坐标之前，根据各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值；在所述最大密度点对应的位置，将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值；继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值；图像绘制模块，用于根据所述第一目标颜色值，所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像。在一些例子中，所述ct采样模块具体用于：将所述pet颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加而得到，得到所述最大密度点对应的第二目标颜色值；其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述pet颜色值确定；当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度。根据本申请实施例的第三方面，提出了一种pet/ct的图像绘制设备，包括：处理器；用于存储所述处理器可执行指令的存储器；其中，所述处理器被配置为执行如所述第一方面所述任一一项方法所述的操作。根据本申请实施例的第四方面，提出一种存储介质，所述机器可读存储介质上存储有若干计算机指令，所述计算机指令被执行时进行如所述第一方面所述任一一项方法所述的操作。本申请提出一种新的pet/ct的图像融合绘制方案，在采样路径上对pet体数据进行采样时，获取所述采样路径上的最大密度点；在所述采样路径上对ct体数据进行采样时，在所述最大密点对应的位置处，将该位置处的pet颜色值融合进来，使绘制的pet/ct图像保持了正确的叠加顺序，确保绘制的pet/ct图像能正确的看出pet的最大密度位置在ct图像中的空间关系。附图说明图1为本申请示例性示出的现有技术中的ct图像和pet图像的示意图；图2为本申请示例性示出的现有技术中pet/ct图像的示意图；图3为本申请实施例示意性示出的一种pet/ct的图像融合绘制方法的部分流程图；图4为本申请实施例示意性示出的一种pet/ct图像的示意图；图5为本申请实施例示意性示出的一种pet/ct的图像融合绘制装置的框图；图6为本申请实施例示意性示出的一种pet/ct的图像融合绘制设备的的逻辑框图。具体实施方式以下将结合附图所示的具体实施方式对本申请进行详细描述。但这些实施方式并不限制本申请，本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本申请的保护范围内。在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。x线计算机体层成像仪(computedtomography，ct)是利用x射线束对人体横断面进行扫描，由探测器接收穿过人体的x射线。x射线在穿过人体时，不同的组织对x射线具有不同的衰减系数或吸收系数，经过计算机的重建就得到人体横断面的断层图像。正电子发射断层显像(positronemissiontomography，pet)的原理是：放射性元素在衰变的过程中会释放出正电子，正电子在遇到电子后发生湮灭，从而在相反的方向上产生一对能量为511kev的光子，由高灵敏度的探测器捕捉到这对光子后经过计算机处理后可以得到正电子与电子湮灭发生的位置。在人体代谢中必需的物质中用半衰期较短的同位素进行标记，注入人体后，通过检测该物质在人体内的聚集情况，来反映代谢活动的情况，从而达到诊断的目的。通过观察代谢物的聚集情况能够发现早期癌症等病变。正电子发射断层显像/x线计算机体层成像仪(positronemissiontomography/computedtomography，pet/ct)是pet和ct的组合体，将pet和ct设计为一体，由一个工作站控制。单pet进行核医学显像时，有其它诊断设备无法比拟的早期发现及灵敏性高等优越特性，但因药物及其原理所限，其定位精度不够，将pet和ct设计为一体，扫描时根据需求同时进行pet显像和ct显像，并将两种图像融合在一起，以达到更好的鉴别和定位。但是，目前的pet/ct的图像融合显示的图像是将根据最大密度投影(maximalintensityprojection，mip)算法获取的pet图像与根据体绘制(volumerendering，vr)算法绘制的ct图像直接加权平均融合成pet/ct图像，由于现有技术中pet/ct图像只是简单的图像融合，不具有深度信息，所以从图像中看不出pet的最大密度位置在图像的具体深度。参照图1，为本申请示例性示出的现有技术中ct图像和pet的图像，图中ct图像110利用体绘制算法获得，pet图像120利用最大密度投影算法获得，将图1中ct图像110与pet图像120融合得到图2所示的pet/ct图像。参照图2，为本申请示例示意性示出的现有技术中pet/ct图像的示意图，从图中可以看出，膀胱121与尾骨111相互重叠，无法得到膀胱121与尾骨111的深度关系，实际上，我们知道，膀胱121应该位于尾骨111之前，但是通过现有的pet/ct图像融合方案，无法得到正确的深度信息。为了解决上述技术问题，本申请提出一种pet/ct的图像融合绘制方法、装置、设备及存储介质，参照图3，为本申请实施例示意性示出的一种pet/ct的图像融合绘制方法的部分流程图，部分步骤如下：s310：在采样路径上对pet体数据进行采样时，获取所述采样路径上的最大密度点坐标及所述最大密度点处的pet颜色值；本申请实施例提出的pet体数据为pet/ct设备中的pet或单独的pet设备扫描出的断层图像序列，它是由多幅连续的断层图像组成。在一些例子中，可以根据光线投影算法、最大密度投影算法等对pet体数据进行采样，本申请并不限制采样的算法。以最大密度投影算法为例，投射出若干条平行光线穿越pet体数据，此时采样路径为光线投射路径，所述光线投射路径垂直与图像投影面，获取出每条光线投射路径上的最大密度点的位置信息。在一些例子中，在采样路径上对pet体数据进行采样时，还可以获取所述采样路径上的最大密度点的点值，所述最大密度点的点值为pet体数据中最大密度点的pet值，通过颜色查找表(lookuptable,lut)将点值映射为pet颜色值。s320：在所述采样路径上，根据预设步长对ct体数据进行采样。本申请实施例提出的ct体数据为pet/ct设备中的ct或单独的ct设备扫描出的断层图像序列，它是由多幅连续的断层图像组成。在一些例子中，可以根据体绘制算法对ct体数据进行采样，体绘制算法可以包括光线投射算法、抛雪球算法及错切-变形法等，本申请并不限制对ct体数据进行采样的算法。所述采样路径与pet/ct图像的投影面垂直，pet/ct图像的每个像素点对应一条采样路径。若采用体绘制算法中的光线投影算法进行采样，那所述采样路径即为光线投影路径，具体采样时，由多束平行光线穿越ct体数据。s330：采样至最大密度点坐标之前，根据各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值。本申请实施例提出的ct颜色值、pet颜色值、第一、第二及第三目标颜色值可以指rgb颜色值、yuv颜色值或灰度值。本申请中，所述各采样点的ct颜色值及各采样点的阻光度根据所述被采样点的ct体数据通过传递函数得到。具体可以是：根据ct设备扫描得到的ct体数据，将所述ct体数据通过传递函数转换为颜色值和阻光度。具体的，从采样开始点为第一个采样点开始采样，第i个采样点的第一目标颜色值可以根据以下公式(1)计算得到：其中，colour1i为第i个采样点的第一目标颜色值；colour0为初始颜色值；colourcti为第i个采样点对应的ct颜色值；αi为i个采样点的阻光度；ai为采样至第i个采样点光线剩余的阻光度。s340：在所述最大密度点对应的位置，将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值。具体的，所述第二目标颜色值为所述pet颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加而得到；其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述pet颜色值确定；当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度。具体的，所述第二目标颜色值可以通过公式(2)计算得到：其中，colour2为第二目标颜色值；colourpet为最大密度点对应的pet颜色值；所述采样路径上最大密度点为第n个采样点；colour1n-1为最大密度点上一个采样点的第一目标颜色值；an为在最大密度点的位置处光线剩余的阻光度；δ为动态融合系数，所述动态融合系数的计算公式可以为：κ*colourpet/255，其中，κ为预设融合系数，所述κ可以由用户设定以调整最大密度点处的pet颜色值的透明度。s350：继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值。具体的，对最大密度点后的ct体数据继续采样，直至采样终点，以最大密度点下一个采样点，即第n+1个采样点继续采样，第m个采样点的第三目标颜色值可以根据如下公式(4)得到：其中，所述colour3m为第m个采样点的第三目标颜色值；colour2为第二目标颜色值；所述采样路径上最大密度点为第n个采样点；colourctm为第m个采样点对应的ct颜色值；αm为m个采样点的阻光度；am为采样至第m个采样点光线剩余的阻光度。s260：根据所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像。在绘制pet/ct图像时，根据每条采样路径上采样开始点，及采样结束点，确定最终显示的pet/ct图像的各像素点的颜色值。所述采样开始点可以为所述采样路径首次接触所述体数据的点，或由用户设定的采样开始点。所述采样结束点可以为采样路径离开所述ct体数据的点。在一些例子中，若所述采样路径为光线投影路径，而在光线采样路径上，光线的剩余阻光度随着采样的进行越来越高，当光线的剩余阻光度达到预设的值时，光线采样路径上剩余的采样点对最终绘制的图像不会做出太大贡献，所以，一些例子中，所述采样结束点为光线剩余阻光度达到预设值的采样点。步骤s260中，绘制pet/ct图像时，若在某条采样路径上，所述最大密度点在采样结束点之后，即未采样至最大密度点处，采样即结束，那么绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第一目标颜色值确定，即所述采样路径对应像素点的颜色值为采样结束点的第一目标颜色值。在一些例子中，若在某条采样路径上，所述最大密度点为采样结束点，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第二目标颜色值确，即所述采样路径对应像素点的颜色值为最大密度点对应的第二目标颜色值。在一些例子中，若在某条采样路径上，所述最大密度点在采样结束点之前，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第三目标颜色值确定，即所述采样路径对应像素点的颜色值为采样结束点的第三目标颜色值。参照图4，为本申请实施例示意性示出的一张pet/ct图像，图中膀胱410遮挡部分尾骨420，可以很清楚地看出膀胱410位于尾骨420的前方。所以可以明显看出，本申请提出的pet/ct的图像融合绘制方案，在采样路径上对pet体数据进行采样时，获取所述采样路径上的最大密度点；在所述采样路径上对ct体数据进行采样时，在所述最大密度点对应的位置处，将所述最大密度点对应的pet颜色值融合进来。使融合绘制的pet/ct图像保持了正确的叠加顺序，能确保pet/ct图像中能正确的看出pet的最大密度位置在ct图像中的空间关系。此外，各个步骤的描述，可以实现为软件、硬件或者其结合的形式，例如，本领域技术人员可以将其实现为软件代码的形式，可以为能够实现所述步骤对应的逻辑功能的计算机可执行指令。当其以软件的方式实现时，所述的可执行指令可以存储在存储器中，并被设备中的处理器执行。与前述pet/ct的图像融合绘制方法的实施例相对应，本申请还提供了pet/ct的图像融合显示装置及pet/ct的图像融合显示设备的实施例。参见图5，为本申请一种pet/ct的图像融合显示装置500的一个实施例框图，该装置包括：pet采样模块510，用于在采样路径上对pet体数据进行采样，获取所述采样路径上的最大密度点坐标及所述最大密度点处的pet颜色值；ct采样模块520，用于在所述采样路径上，根据预设步长对ct体数据进行采样，采样至最大密度点坐标之前，根据各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值；在所述最大密度点对应的位置，将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值；继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值；图像绘制模块530，用于根据所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像。在一些例子中，所述将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值，ct采样模块520具体用于：将所述pet颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加，得到所述最大密度点对应的第二目标颜色值；其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述pet颜色值确定；当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度。在一些例子中，所述根据所述第一目标颜色值，所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像，具体为：若在所述采样路径上，所述最大密度点在采样结束点之前，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第一目标颜色值确定；或若在所述采样路径上，所述最大密度点为采样结束点，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第二目标颜色值确定；或若在所述采样路径上，所述最大密度点在采样结束点之后，绘制的pet/ct图像中，所述采样路径对应的像素点根据所述第三目标颜色值确定。在一些例子中，所述各采样点的ct颜色值及各采样点的阻光度根据所述被采样点的ct体数据通过传递函数转换得到。在一些例子中，根据最大密度投影算法或光线投影算法对pet体数据进行采样；根据体绘制算法对ct体数据进行采样。在一些例子中，所述采样路径为光线投射路径，所述光线投射路径穿越pet体数据或ct体数据进行采样。参见图6，为本申请提出一种pet/ct的图像绘制设备的一个实施例示意图，该pet/ct的图像绘制设备600可以包括：存储处理器可执行指令的存储器602；其中，所述处理器601耦合于所述存储器602，用于读取所述存储器602存储的程序指令，并作为响应，执行如下操作：在采样路径上对pet体数据进行采样，获取所述采样路径上的最大密度点坐标及所述最大密度点处的pet颜色值；在所述采样路径上，根据预设步长对ct体数据进行采样，采样至最大密度点坐标之前，根据各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值；在所述最大密度点对应的位置，将所述pet颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值；继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的ct颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值；根据所述第一目标颜色值，所述第二目标颜色值及所述第三目标颜色值，绘制pet/ct图像。需要说明的是，图6所述的pet/ct的图像绘制设备可以是pet/ct设备、与pet/ct设备连接的控制台设备、pet设备、ct设备、与pet设备和/或ct设备连接的控制台设备、或任意具有处理能力的电子设备，例如工作站或计算机等。在本申请实施例中，计算机可读存储介质可以是多种形式，比如，在不同的例子中，所述机器可读存储介质可以是：ram(radomaccessmemory，随机存取存储器)、易失存储器、非易失性存储器、闪存、存储驱动器(如硬盘驱动器)、固态硬盘、任何类型的存储盘(如光盘、dvd等)，或者类似的存储介质，或者它们的组合。特殊的，所述的计算机可读介质还可以是纸张或者其他合适的能够打印程序的介质。使用这些介质，这些程序可以被通过电学的方式获取到(例如，光学扫描)、可以被以合适的方式编译、解释和处理，然后可以被存储到计算机介质中。本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本
技术领域：
：中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本申请的真正范围和精神由本申请的权利要求指出。应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本申请的范围仅由所附的权利要求来限制。当前第1页12当前第1页12