CN109272476B - Pet/ct的图像融合绘制方法、装置、设备及存储介质 - Google Patents
Pet/ct的图像融合绘制方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN109272476B CN109272476B CN201810879996.6A CN201810879996A CN109272476B CN 109272476 B CN109272476 B CN 109272476B CN 201810879996 A CN201810879996 A CN 201810879996A CN 109272476 B CN109272476 B CN 109272476B
- Authority
- CN
- China
- Prior art keywords
- color value
- sampling
- point
- pet
- target color
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000003286 fusion draw glass process Methods 0.000 title abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 251
- 230000004927 fusion Effects 0.000 claims description 71
- 238000009877 rendering Methods 0.000 claims description 31
- 238000004422 calculation algorithm Methods 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000000903 blocking effect Effects 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 6
- 238000005266 casting Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims 1
- 238000002600 positron emission tomography Methods 0.000 description 138
- 238000002591 computed tomography Methods 0.000 description 131
- 238000010586 diagram Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 238000009206 nuclear medicine Methods 0.000 description 4
- 238000012879 PET imaging Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 239000002207 metabolite Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/50—Lighting effects
-
- 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/10104—Positron emission tomography [PET]
-
- 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/20—Special algorithmic details
- G06T2207/20068—Projection on vertical or horizontal image axis
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请公开了一种PET/CT的图像融合绘制方法、装置、设备及存储介质，该方法包括：在采样路径上对PET体数据进行采样时，获取所述采样路径上的最大密度点及所述最大密度点处的PET颜色值；在所述采样路径上对CT体数据进行采样时，在所述最大密点对应的位置处，将该位置处的PET颜色值融合进来，使绘制的PET/CT图像保持了正确的叠加顺序，确保绘制的PET/CT图像能正确的看出PET的最大密度位置在CT图像中的空间关系。
Description
技术领域
本申请涉及医学图像处理技术领域，特别涉及PET/CT的图像融合绘制方法、装置、设备及存储介质。
背景技术
正电子发射断层显像/X线计算机体层成像仪(Positron Emission Tomography/Computedtomography，PET/CT)是正电子发射断层显像(Positron Emission Tomography，PET)和X线计算机体层成像仪(Computed Tomography，CT)组合体，将PET和CT设计为一体，由一个工作站控制。单PET进行核医学显像时，有其它诊断设备无法比拟的早期发现及灵敏性高等优越特性，但因药物及其原理所限，其定位精度不够，将PET和CT设计为一体，扫描时根据需求同时进行PET显像和CT显像，并将两种图像融合在一起，以达到更好的鉴别和定位。
但是，目前的PET/CT的图像融合显示的图像是将根据最大密度投影(MaximalIntensity Projection，MIP)算法获取的PET图像与根据体绘制(Volume Rendering，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线计算机体层成像仪(Computed Tomography，CT)是利用X射线束对人体横断面进行扫描，由探测器接收穿过人体的X射线。X射线在穿过人体时，不同的组织对X射线具有不同的衰减系数或吸收系数，经过计算机的重建就得到人体横断面的断层图像。
正电子发射断层显像(Positron Emission Tomography，PET)的原理是：放射性元素在衰变的过程中会释放出正电子，正电子在遇到电子后发生湮灭，从而在相反的方向上产生一对能量为511KeV的光子，由高灵敏度的探测器捕捉到这对光子后经过计算机处理后可以得到正电子与电子湮灭发生的位置。在人体代谢中必需的物质中用半衰期较短的同位素进行标记，注入人体后，通过检测该物质在人体内的聚集情况，来反映代谢活动的情况，从而达到诊断的目的。通过观察代谢物的聚集情况能够发现早期癌症等病变。
正电子发射断层显像/X线计算机体层成像仪(Positron Emission Tomography/Computedtomography，PET/CT)是PET和CT的组合体，将PET和CT设计为一体，由一个工作站控制。单PET进行核医学显像时，有其它诊断设备无法比拟的早期发现及灵敏性高等优越特性，但因药物及其原理所限，其定位精度不够，将PET和CT设计为一体，扫描时根据需求同时进行PET显像和CT显像，并将两种图像融合在一起，以达到更好的鉴别和定位。
但是，目前的PET/CT的图像融合显示的图像是将根据最大密度投影(MaximalIntensity Projection，MIP)算法获取的PET图像与根据体绘制(Volume Rendering，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值，通过颜色查找表(Look Up Table,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个采样点的第一目标颜色值；
Colour 0为初始颜色值；
ColourCTi为第i个采样点对应的CT颜色值；
αi为i个采样点的阻光度；
Ai为采样至第i个采样点光线剩余的阻光度。
S340：在所述最大密度点对应的位置，将所述PET颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值。
具体的，所述第二目标颜色值为所述PET颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加而得到；其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述PET颜色值确定；当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度。
具体的，所述第二目标颜色值可以通过公式(2)计算得到：
其中，Colour 2为第二目标颜色值；
Colour PET为最大密度点对应的PET颜色值；
所述采样路径上最大密度点为第n个采样点；
Colour 1n-1为最大密度点上一个采样点的第一目标颜色值；
An为在最大密度点的位置处光线剩余的阻光度；
δ为动态融合系数，所述动态融合系数的计算公式可以为：
κ*Colour PET/255，其中，κ为预设融合系数，所述κ可以由用户设定以调整最大密度点处的PET颜色值的透明度。
S350：继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的CT颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值。
具体的，对最大密度点后的CT体数据继续采样，直至采样终点，以最大密度点下一个采样点，即第n+1个采样点继续采样，第m个采样点的第三目标颜色值可以根据如下公式(4)得到：
其中，所述Colour 3m为第m个采样点的第三目标颜色值；
Colour 2为第二目标颜色值；
所述采样路径上最大密度点为第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(Radom Access Memory，随机存取存储器)、易失存储器、非易失性存储器、闪存、存储驱动器(如硬盘驱动器)、固态硬盘、任何类型的存储盘(如光盘、dvd等)，或者类似的存储介质，或者它们的组合。特殊的，所述的计算机可读介质还可以是纸张或者其他合适的能够打印程序的介质。使用这些介质，这些程序可以被通过电学的方式获取到(例如，光学扫描)、可以被以合适的方式编译、解释和处理，然后可以被存储到计算机介质中。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本申请的真正范围和精神由本申请的权利要求指出。
应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本申请的范围仅由所附的权利要求来限制。
Claims (10)
1.一种PET/CT的图像绘制方法，其特征在于，所述方法包括：
在采样路径上对PET体数据进行采样，获取所述采样路径上的最大密度点坐标及所述最大密度点处的PET颜色值；
在所述采样路径上，根据预设步长对CT体数据进行采样，采样至最大密度点坐标之前，根据各采样点的CT颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值；
所述计算所述各采样点的第一目标颜色值，包括：对于第i个采样点，将第1个采样点到所述第i个 采样点之间的每个采样点的CT颜色值、阻光度和光线剩余的阻光度相乘，得到每个采样点对应的CT颜色值、阻光度和光线剩余的阻光度的乘积，将乘积求和并与初始颜色值相加，得到所述第i个采样点的第一目标颜色值，所述第i个采样点为所述各采样点中任一个采样点；
在所述最大密度点对应的位置，将所述PET颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值；
继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的CT颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值；
根据所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值，绘制PET/CT图像，包括：依据最大密度点与采样结束点之间的关系，从三种目标颜色值：所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值中，选取一种目标颜色值用于确定绘制的PET/CT图像中采样路径对应的像素点。
2.根据权利要求1所述的方法，其特征在于，所述根据所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值，绘制PET/CT图像，包括：
若在所述采样路径上，所述最大密度点在采样结束点之后，绘制的PET/CT图像中，所述采样路径对应的像素点根据所述第一目标颜色值确定；或
若在所述采样路径上，所述最大密度点为采样结束点，绘制的PET/CT图像中，所述采样路径对应的像素点根据所述第二目标颜色值确定；或
若在所述采样路径上，所述最大密度点在采样结束点之前，绘制的PET/CT图像中，所述采样路径对应的像素点根据所述第三目标颜色值确定。
3.根据权利要求1所述的方法，其特征在于，所述各采样点的CT颜色值及各采样点的阻光度根据所述采样点的CT体数据通过传递函数转换得到。
4.根据权利要求1所述的方法，其特征在于，所述方法还包括：根据最大密度投影算法或光线投影算法对PET体数据进行采样；或/和
根据体绘制算法对CT体数据进行采样。
5.根据权利要求1所述的方法，其特征在于，所述采样路径为光线投射路径，所述光线投射路径穿越PET体数据或CT体数据进行采样。
6.根据权利要求5所述的方法，其特征在于，所述将所述PET颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值，包括：
所述第二目标颜色值为所述PET颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加而得到；
其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述PET颜色值确定；
当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度。
7.一种PET/CT的图像绘制装置，其特征在于，所述装置包括：
PET采集模块，用于在采样路径上对PET体数据进行采样，获取所述采样路径上的最大密度点坐标及所述最大密度点处的PET颜色值；
CT采集模块，用于在所述采样路径上，根据预设步长对CT体数据进行采样，采样至最大密度点坐标之前，根据各采样点的CT颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述各采样点的第一目标颜色值；所述计算所述各采样点的第一目标颜色值，包括：对于第i个采样点，将第1个采样点到所述第i个 采样点之间的每个采样点的CT颜色值、阻光度和光线剩余的阻光度相乘，得到每个采样点对应的CT颜色值、阻光度和光线剩余的阻光度的乘积，将乘积求和并与初始颜色值相加，得到所述第i个采样点的第一目标颜色值，所述第i个采样点为所述各采样点中任一个采样点；在所述最大密度点对应的位置，将所述PET颜色值与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行融合计算，得到所述最大密度点对应的第二目标颜色值；继续在所述采样路径上进行采样，根据所述第二目标颜色值、各采样点的CT颜色值、各采样点的阻光度及各采样点的光线剩余阻光度，计算所述采样点的第三目标颜色值；
图像绘制模块，用于根据所述第一目标颜色值，所述第二目标颜色值及所述第三目标颜色值，绘制PET/CT图像，包括：依据最大密度点与采样结束点之间的关系，从三种目标颜色值：所述第一目标颜色值、所述第二目标颜色值及所述第三目标颜色值中，选取一种目标颜色值用于确定绘制的PET/CT图像中采样路径对应的像素点。
8.根据权利要求7所述的装置，其特征在于，所述CT采集模块具体用于：
将所述PET颜色值与融合系数的乘积，与距离所述最大密度点最近的上一个采样点的所述第一目标颜色值进行叠加而得到，得到所述最大密度点对应的第二目标颜色值；
其中，当所述最大密度点的坐标处光线的剩余阻光度大于动态融合系数时，所述融合系数为动态融合系数，其中，所述动态融合系数根据所述PET颜色值确定；
当所述最大密度点的坐标处光线的剩余阻光度不大于动态融合系数时，所述融合系数为所述光线的剩余阻光度。
9.一种PET/CT的图像绘制设备，其特征在于，包括：
存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其中，所述处理器执行所述程序时实现如所述权利要求1至6任意一项所述的方法。
10.一种机器可读存储介质，其特征在于，所述机器可读存储介质上存储有若干计算机指令，所述计算机指令被执行时实现如所述权利要求1至6任意一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810879996.6A CN109272476B (zh) | 2018-08-03 | 2018-08-03 | Pet/ct的图像融合绘制方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810879996.6A CN109272476B (zh) | 2018-08-03 | 2018-08-03 | Pet/ct的图像融合绘制方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109272476A CN109272476A (zh) | 2019-01-25 |
| CN109272476B true CN109272476B (zh) | 2021-12-31 |
Family
ID=65153084
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810879996.6A Active CN109272476B (zh) | 2018-08-03 | 2018-08-03 | Pet/ct的图像融合绘制方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109272476B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112397198B (zh) * | 2020-11-17 | 2024-08-16 | 上海商汤善萃医疗科技有限公司 | 图像处理方法及装置、电子设备和存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6807247B2 (en) * | 2002-03-06 | 2004-10-19 | Siemens Corporate Research, Inc. | Visualization of volume—volume fusion |
| CN101004838A (zh) * | 2007-01-19 | 2007-07-25 | 哈尔滨工程大学 | 阻光度传递函数的构造方法 |
| US7711162B2 (en) * | 2005-01-27 | 2010-05-04 | Siemens Medical Solutions Usa, Inc. | Accelerated texture-based fusion renderer |
| CN104143190A (zh) * | 2014-07-24 | 2014-11-12 | 东软集团股份有限公司 | Ct图像中组织的分割方法及系统 |
-
2018
- 2018-08-03 CN CN201810879996.6A patent/CN109272476B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6807247B2 (en) * | 2002-03-06 | 2004-10-19 | Siemens Corporate Research, Inc. | Visualization of volume—volume fusion |
| US7711162B2 (en) * | 2005-01-27 | 2010-05-04 | Siemens Medical Solutions Usa, Inc. | Accelerated texture-based fusion renderer |
| CN101004838A (zh) * | 2007-01-19 | 2007-07-25 | 哈尔滨工程大学 | 阻光度传递函数的构造方法 |
| CN104143190A (zh) * | 2014-07-24 | 2014-11-12 | 东软集团股份有限公司 | Ct图像中组织的分割方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109272476A (zh) | 2019-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5172347B2 (ja) | 反復的制約デコンボリューションによる核医学２次元平面画像の復元 | |
| EP1522045B1 (en) | Motion artifact correction of tomographical images | |
| US7348564B2 (en) | Multi modality imaging methods and apparatus | |
| EP2668639B1 (en) | Truncation compensation for iterative cone-beam ct reconstruction for spect/ct systems | |
| CN108242068B (zh) | 用于成像系统的方法和ct/pet系统 | |
| EP2174294B1 (en) | Motion correction in nuclear imaging | |
| US8421021B2 (en) | Motion correction of SPECT images | |
| EP1934942B1 (en) | Iterative reconstruction with enhanced noise control filtering | |
| JP4414410B2 (ja) | 画像再構成方法 | |
| Martinez-Möller et al. | Attenuation correction for PET/MR: problems, novel approaches and practical solutions | |
| US20180061031A1 (en) | Model-Based Scatter in Multi-Modality Multi-Energy SPECT Reconstruction | |
| Funama et al. | A newly-developed metal artifact reduction algorithm improves the visibility of oral cavity lesions on 320-MDCT volume scans | |
| JP6662880B2 (ja) | 放射線放出撮像システム、記憶媒体及び撮像方法 | |
| US11816764B2 (en) | Partial volume correction in multi-modality emission tomography | |
| JP2008537113A (ja) | 粗い角度及びスライスリビニングによる三次元飛行時間ｐｅｔ | |
| WO2015198189A1 (en) | Reconstruction with multiple photopeaks in quantitative single photon emission computed tomography | |
| JP4347651B2 (ja) | マルチ・モダリティ・イメージング方法及び装置 | |
| JP6526428B2 (ja) | 医用画像処理装置、医用画像処理方法および医用画像診断装置 | |
| CN109272476B (zh) | Pet/ct的图像融合绘制方法、装置、设备及存储介质 | |
| KR102026735B1 (ko) | 영상 촬영 장치의 검출기의 시스템 응답 및 시스템 응답을 이용하여 의료 영상을 생성하는 방법 및 장치 | |
| EP4148680A1 (en) | Attenuation correction-based weighting for tomographic inconsistency detection | |
| JP2008185335A (ja) | 医用画像アニメーション表示装置 | |
| KR101331504B1 (ko) | 핵의학 방출 영상 시스템 및 단일 투과 영상을 이용한 감쇠보정 방법 | |
| WO2021112821A1 (en) | Network determination of limited-angle reconstruction | |
| Us | Reduction of Limited Angle Artifacts in Medical Tomography via Image Reconstruction |
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 | ||
| CB02 | Change of applicant information |
Address after: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant after: DongSoft Medical System Co., Ltd. Address before: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant before: Dongruan Medical Systems Co., Ltd., Shenyang |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |