CN106548464B - 一种构建相位图像的方法及装置 - Google Patents
一种构建相位图像的方法及装置 Download PDFInfo
- Publication number
- CN106548464B CN106548464B CN201610976174.0A CN201610976174A CN106548464B CN 106548464 B CN106548464 B CN 106548464B CN 201610976174 A CN201610976174 A CN 201610976174A CN 106548464 B CN106548464 B CN 106548464B
- Authority
- CN
- China
- Prior art keywords
- phase
- image
- target
- phase image
- transition matrix
- 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 22
- 239000011159 matrix material Substances 0.000 claims abstract description 55
- 230000007704 transition Effects 0.000 claims abstract description 47
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 18
- 230000004927 fusion Effects 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000002591 computed tomography Methods 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 9
- 230000001419 dependent effect Effects 0.000 abstract description 6
- 230000033001 locomotion Effects 0.000 description 28
- 230000003902 lesion Effects 0.000 description 22
- 238000012545 processing Methods 0.000 description 20
- 230000006870 function Effects 0.000 description 14
- 210000000056 organ Anatomy 0.000 description 14
- 230000015654 memory Effects 0.000 description 13
- 239000002207 metabolite Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 201000010099 disease Diseases 0.000 description 6
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 5
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- 238000012937 correction Methods 0.000 description 4
- VYFYYTLLBUKUHU-UHFFFAOYSA-N dopamine Chemical compound NCCC1=CC=C(O)C(O)=C1 VYFYYTLLBUKUHU-UHFFFAOYSA-N 0.000 description 4
- 210000002216 heart Anatomy 0.000 description 4
- 210000004072 lung Anatomy 0.000 description 4
- 230000018199 S phase Effects 0.000 description 3
- 238000002372 labelling Methods 0.000 description 3
- 210000004185 liver Anatomy 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 229960003638 dopamine Drugs 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000009206 nuclear medicine Methods 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 239000012217 radiopharmaceutical Substances 0.000 description 2
- 229940121896 radiopharmaceutical Drugs 0.000 description 2
- 230000002799 radiopharmaceutical effect Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- WQZGKKKJIJFFOK-VFUOTHLCSA-N beta-D-glucose Chemical compound OC[C@H]1O[C@@H](O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-VFUOTHLCSA-N 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000000069 prophylactic effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000000941 radioactive substance Substances 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Nuclear Medicine (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种构建相位图像的方法及装置,属于医疗器械领域。所述方法包括:根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T‑1,T为相位总数目,其值为大于1的整数;构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵将所述T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。所述装置包括:第一构建模块、第二构建模块和合成模块。本发明能够降低构建的相位图像的噪声。
Description
技术领域
本发明涉及医疗器械领域,特别涉及一种构建相位图像的方法及装置。
背景技术
PET/CT(Positron emission tomography/Computedtomography,正电子发射断层显像/X线计算机体层成像仪)设备是一种高端核医学影像诊断设备,其包括PET设备和CT设备。PET设备用于构建人体器官内的病灶图像,CT设备用于拍摄人体的CT图像,这样医生根据该CT图像并结合该病灶图像,能及早发现病灶、准确定位病灶,以及帮助医生诊断疾病。
目前在构造病灶图像时,PET设备先对病人进行扫描多个运动周期,每个运动周期划分成T个相位,然后根据每个运动周期的第t个相位内PET设备采集的采集数据,构建第t个相位的相位图像,t为大于或等于0且小于T的整数;其中,第t个相位的相位图像就是包括人体器官的病灶的病灶图像。
在实现本发明的过程中,发明人发现现有技术至少存在以下问题:
目前在重建某个相位的相位图像时,利用该相位内采集的采集数据来构建该相位的相位图像,使得构建出的相位图像的噪声较高。
发明内容
为了降低构建的相位图像的噪声,本发明提供了一种构建相位图像的方法及装置。所述技术方案如下:
第一方面,提供了一种构建相位图像的方法,所述方法包括:
根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数;
构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;
根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵,将所述T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。
可选的,所述构建所述第t个相位图像的转换矩阵,包括:
获取电子计算机断层扫描4D-CT设备构建的所述第t个相位的CT图像和所述目标相位的CT图像;
根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵。
可选的,所述根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵,包括:
根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
可选的,所述根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵,将所述T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数,包括:
按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
可选的,所述将所述T个相位图像合成目标相位图像之后,还包括:
将所述目标相位图像和所述目标相位的CT图像融合成一幅图像,显示所述融合的一幅图像。
第二方面,提供了一种构建相位图像的装置,所述装置包括:
第一构建模块,用于根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数;
第二构建模块,用于构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;
合成模块,用于根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵将所述T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。
可选的,所述第二构建模块包括:
获取单元,用于获取电子计算机断层扫描4D-CT设备构建的所述第t个相位的CT图像和所述目标相位的CT图像;
构建单元,用于根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵。
可选的,所述构建单元,用于根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
可选的,所述合成模块,用于
按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
可选的,所述装置还包括:
融合模块,用于将所述目标相位图像和所述目标相位的CT图像融合成一幅图像,显示所述融合的一幅图像。
本发明提供的技术方案的有益效果是:
通过T个相位图像中的每个相位图像的转换矩阵,将该T个相位的相位图像合成一个目标相位图像,这样利用了全部相位内各探测器对的采集数据来构建一幅目标相位图像,从而降低了目标相位图像的噪声。
附图说明
图1-1是本发明实施例提供的一种PET设备结构示意图;
图1-2是本发明实施例提供的一种γ光子对产生示意图;
图1-3是本发明实施例1提供的一种构建相位图像的方法流程图;
图2-1是本发明实施例2提供的一种构建相位图像的方法流程图;
图2-2是本发明实施例2提供的一种PET/CT设备结构示意图;
图3是本发明实施例3提供的一种构建相位图像的装置结构示意图;
图4是本发明实施例4提供的一种处理终端结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
在PET/CT系统中,CT扫描通常可以在很短时间内完成,获得的图像几乎是某时刻的快照。但是PET扫描速度较慢,每个体位通常会花费几分钟时间,因此不可能在病人屏气状态下完成数据采集。受心脏搏动和呼吸运动的影响,PET图像和CT图像对相同病灶在位置及相位上存在一定程度的不匹配。一方面,PET采集会将运动的病灶在不同位置时采集的数据叠加,对应的成像反映的是病灶位置的一个平均效应,这不可避免的带来了分辨率的下降,导致图像模糊和定量分析SUV值(标准化摄取值)准确性降低。另一方面,CT图像与PET图像配准融合产生偏差(特别在运动幅度最大的膈膜附近),用瞬间的CT图像去对平均的PET图像做衰减校正必然出现误差,会在PET图像上产生局部伪影,可能影响胸腹部肿瘤的准确诊断与治疗计划的制定。
为了去除心脏搏动或校正呼吸运动所造成的影响,目前在PET/CT上通常使用门控技术,即通过特定的设备采集心脏或呼吸的运动信号,并以此为基础产生门控信号并发送给PET/CT采集控制系统。PET/CT采集控制系统接收到门控信号后,根据一定的阈值触发PET,与运动信号同步采集病灶的图像信息。采集到的数据按周期时相归类分析,保证在运动周期内的同一相位区间内病灶近似保持静止。这样将相同时相内的PET图像进行重建,减小了运动伪影,改善图像质量,提高了病灶定位的准确度和SUV值测量准确性。所有周期相位的图像综合分析,可以得到病灶三维空间影像随时间变化的运动信息,对临床诊断、分期和病灶放疗靶区的勾画有着重要的意义。另外,选择与CT位置最相匹配的相位,其衰减校正受运动影响最小,这提高了图像融合的精度和衰减校正的准确性,最终进一步提高了病灶SUV值测量准确性。
在PET/CT系统中,PET模块是一种高端核医学影像诊断设备,其可以对病人进行扫描,然后构造病人体内的病灶图像,帮助医生及早发现病灶。参见图1-1,PET设备包括处理终端1、一个或多个探测器环2。每个探测器环2由多个探测器3组成,在实际实现时每个探测器环2可以由若干探测器3组成,这些探测器3可以为闪烁探测器,终端设备1可以为计算机等。
在实际扫描过程中,可先利用放射性核素对代谢物质进行标记,然后可以采用注射方式向病人注射利用放射性核素标记的代谢物质。该代谢物质进入病人体内,会被病人体内的器官所吸收,以及用于标记该代谢物质的放射性核素也连同该代谢物质一起被吸收。其中,需要说明的是:如果病人体内的某个器官有病灶,该病灶对放射性核素标记的代谢物质摄取高。例如,假设病人的肝内生长了肿瘤,由于肿瘤细胞的生长速度远大于病人体内的正常细胞,使该肿瘤对代谢物质的摄取要远多于病人体内的其他正常部位。
仍参见图1-1,病人在注射放射性药物之后,将病人放置在探测器环2中,让探测器环2对病人进行扫描。由于放射性核素的衰变,放射性示踪剂在被病人体内的器官和病灶部位吸收后,用于标记的放射性核素很快就会湮灭,并产生一对方向完全相反的一对γ光子,两者的方向互成180度。例如,参见图1-2,放射性核素在湮灭时产生了一对γ光子,分别为γ光子1和γ光子2,且γ光子1的方向和γ光子2的方向互成180度。在探测器环2扫描病人的过程中,每个放射性核素湮灭时产生的一对γ光子会同时被探测器环2上的两个探测器所接收。处理终端1将该两个探测器组成一对探测器对,并记录该对探测器对的标识信息和接收该γ光子对的时间等数据,该探测器对的标识信息可以由该两个探测器的位置信息组成。
其中,放射性核素可以为18F、11C等,放射性药物可以为放射性核素标记的葡萄糖、多巴胺等。让病人注射的代谢物质与病人需要检测的器官有关。如果病人需要检测的器官为心脏、肺、肝脏等器官,让病人注射的放射性药物可以为FDG。如果病人需要检测的器官为大脑,让病人注射的代谢物质可以为多巴胺。
其中,需要说明的是:PET设备可以对病人扫描一段时间,该段时间为N个运动周期,N为大于1的整数,例如N可以为10、15或20等数值。由于受病人体内的心脏搏动或肺部的呼吸运动的影响,病人体内的器官和病灶的位置是按照运动周期周期性地变动。如果该器官是心脏,则该器官的运动周期为心脏的跳动周期,如果该器官为肺或肝脏等器官,则该器官的运动周期为肺的呼吸周期。在利用PET设备扫描病人的实际操作中,会使用探测器环2对病人扫描N个运动周期,处理终端1会记录下该N个运动周期的每个运动周期内接收γ光子对的探测器对的标识信息和该探测器对接收的γ光子对的时间等数据。处理终端1还将每个运动周期划分成T个时间段,每个时间段称为一个相位,所以每个运动周期包括T个相位,T为大于1的整数,例如T可以为4、5或6等数值。
在扫描结束后,对于第t个相位,处理终端1记录了N个第t个相位内接收γ光子对的各探测器对和各探测器对接收γ光子对的时间等数据。然后,处理终端1根据该第t个相位内接收γ光子对的各探测器对接收γ光子对的时间,统计该各探测器对接收γ光子对的数目,得到各探测器对的采集数据。最后,处理终端1根据该N个运动周期内的第t个相位内各探测器对的采集数据,构建出一个第t个相位的相位图像,该第t个相位的相位图像即为包括病人的器官和病灶的图像,可以呈现给医生,帮助医生进行疾病的诊断。在实现本发明技术方向时,还发现目前处理终端1只是利用了T分之一的采集数据来构造第t个相位的相位图像,剩下的T-1份采集数据都未用于构造第t个相位的相位图像,导致构造的相位图像的噪声较大。在本发明中可以通过如下任一实施例来构造出某个相位的相位图像,且构造的相位图像的噪声较低。
实施例1
参见图1-3,本发明实施例提供了一种构建相位图像的方法,所述方法包括:
步骤101:根据第t个相位内接收γ光子对的各探测器对的采集数据,构建第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数。
步骤102:构建第t个相位图像的转换矩阵,该转换矩阵用于将第t个相位图像转换为目标相位的相位图像,目标相位是所述T个相位中的一个相位。
步骤103:根据该T个相位图像、该T个相位图像中的每个相位图像的转换矩阵,将该T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。
可选的,上述步骤102,可以包括:
1021:获取4D-CT设备构建的第t个相位的CT图像和目标相位的CT图像。
1022:根据第t个相位的CT图像和目标相位的CT图像,构建第t个相位的转换矩阵。
可选的,上述步骤1022,可以为:
根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
可选的,上述步骤103,可以为:
按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
可选的,在上述步骤103之后,该方法还可以包括:
将目标相位图像和目标相位的CT图像融合成一幅图像,显示融合的一幅图像。
在本发明实施例中,通过T个相位图像中的每个相位图像的转换矩阵,将该T个相位的相位图像合成一个目标相位图像,这样利用了全部相位内各探测器对的采集数据来构建一幅目标相位图像,从而降低了目标相位图像的噪声。
实施例2
本发明实施例提供了一种构建相位图像的方法。本实施例的执行主体可以为上述实施例介绍的处理终端,且该构建方法是在利用PET设备对病人进行扫描后执行的。参见图2-1,该方法包括:
步骤201:对于N个运动周期中的每个运动周期,确定在该运动周期的第t相位内接收γ光子对的探测器对以及获取该各探测器对的采集数据。
其中,在对病人扫描的过程中由于处理终端记录了接收γ光子对的探测器对的标识信息和各探测器对接收γ光子对的时间等信息,因此根据记录的数据可以确定出在该运动周期的第t个相位内接收γ光子对的各探测器对,根据记录该各探测器对接收γ光子对的时间,统计出各探测器对接收γ光子对的数目,即得到在该运动周期的第t相位内该各探测器对的采集数据。运动周期包括N个,都按上述方式得到N个第t相位内各探测器对的采集数据。在本实施例中,t=0、1、2……T-1。
步骤202:根据N个运动周期的第t个相位内的各探测器对的采集数据,构建第t个相位的相位图像。
在PET设备中事先定义了如下公开(1)所示的相位图像构建模型:
在上述公式中y=[y1,y2,…,yM]T,表示各探测器对的采集数据,M为探测器对的总数目,A为系统矩阵,用数学的形式表达了PET设备中每个空间位置的点源被探测器探测到的概率,反映了系统的物理特性,r表示噪声的平均值,x表示待构建的相位图像,E[.]表示期望值算子。
根据公式(1)的构建模型可以得出在第t个相位的相位图像构建模型,如下公式(2)所示。
E[yt]=A·xt+rt……(2)
在上述公式(2)中,t=0,1,…,T-1,T为一个运动周期包括的相位总数,yt为第t个相位内各探测器对的采集数据,xt为第t相位的待构建相位图像,rt为第t相位的噪声的平均值,yt=[y1t,y2t,…,yMt]T。
其中,E[yt]是一个求平均值运算算子,表示去除掉统计噪声后的期望值。所对应的其中,表示探测器对1在每个运动周期的第t相位内的采集数据的平均值,的含义与的含义相同,就不一一说明。
根据上述公式(2),可以得出第t相位的相位图像xt的估计值如下公式(3)所示:
其中,t=0,1,2……,T-1,m=1,2……,M,L(xt,yt)是针对相位t的log-likelihood函数,如下公式(4)所示:
Pr为系统矩阵表示的概率函数,针对的是物体在第t相位时采集到的数据。
由于每个运动周期具有T个相位,因此对每个相位都按上述方式得到T个相位的相位图像,求解公式(4)可以通过传统的OSEM或MLEM算法。
步骤203:构建第t个相位图像的转换矩阵,该转换矩阵用于将第t个相位图像转换为目标相位的相位图像,目标相位是T个相位中的一个相位。
本步骤可以通过如下第一和第二两个步骤来实现,包括:
第一步:获取4D-CT设备构建的第t个相位的CT图像和目标相位的CT图像。
参见图2-2,通常PET设备和4D-CT设备都是结合使用,两者并排设置,在对病人扫描时,首先病人先经过4D-CT设备再经过PET设备,PET设备和4D-CT设备依次对病人扫描。其中,4D-CT设备在对病人扫描的过程中立即就可以构建出病人的CT图像。在本步骤中,可以从4D-CT设备中获取第t个相位的CT图像和目标相位的CT图像。
第二步:根据第t个相位的CT图像和目标相位的CT图像,构建第t个相位的转换矩阵。
本步骤可以为:根据第t个相位的CT图像和目标相位的CT图像,通过B样条函数构建第t个相位的转换矩阵。
在实际实现时,对第t个相位的CT图像中的任一个像素点,为了便于说明称该像素点为第一像素点,从目标相位的CT图像中获取该第一像素点对应的目标像素点,第一像素点在第t个相位的CT图像中的位置与目标像素点在目标相位的CT图像中的位置不同,通过B样条函数计算出第一像素点和目标像素点的差异信息;按上述方法计算出第t个相位的CT图像中的其他每个像素点的差异信息,将第t个相位的CT图像中的每个像素点的差异信息组成第t个相位的转换换矩阵。
步骤204:根据该T个相位图像、该T个相位图像中的每个相位图像的转换矩阵,将该T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。
具体的,按如下公式将该T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中各探测器对的采集数据。
其中,在合成得到目标相位图像后,还可以将目标相位图像与目标相位的CT图像融合成一幅图像,并显示该一幅图像给医生。
在本发明实施例中,将T个相位的相位图像合成一个目标相位图像,这样利用了全部相位内各探测器对的采集数据来构建一幅目标相位图像,从而降低了目标相位图像的噪声。将该目标相位图像呈现给医生,能够更好地帮助医生诊断疾病。
实施例3
参见图3,本发明实施例提供了一种构建相位图像的装置300,所述装置300包括:
第一构建模块301,用于根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数;
第二构建模块302,用于构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;
合成模块303,用于根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵将所述T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。
可选的,所述第二构建模块302包括:
获取单元,用于获取电子计算机断层扫描4D-CT设备构建的所述第t个相位的CT图像和所述目标相位的CT图像;
构建单元,用于根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵。
可选的,所述构建单元,用于根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
可选的,所述合成模块303,用于
按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
进一步地,所述装置300还包括:
融合模块,用于将所述目标相位图像和所述目标相位的CT图像融合成一幅图像,显示所述融合的一幅图像。
在本发明实施例中,将T个相位的相位图像合成一个目标相位图像,这样利用了全部相位内各探测器对的采集数据来构建一幅目标相位图像,从而降低了目标相位图像的噪声。将该目标相位图像呈现给医生,能够更好地帮助医生诊断疾病。
实施例4
参见图4,本发明实施例提供了一种处理终端400。该处理终端400用于实施上述实施例中提供的构建相位图像的方法。具体来讲:
处理终端400可以包括处理器410、收发器420、存储器430、输入单元440、显示单元450、音频电路460以及电源470等部件,如图4所示,本领域技术人员可以理解,图4中示出的终端结构并不构成对终端的限定,可以包括比图示更多或更少的部件或者组合某些部件,或者不同的部件布置。其中:
处理器410可以是终端的控制中心,利用各种接口和线路连接整个终端设备的各个部分,如收发器420和存储器430等,通过运行或执行存储在存储器430内的软件程序和/或模块,以及调用存储在存储器430内的数据,执行处理终端400的各种功能和处理数据,从而对处理终端400进行整体监控。可选的,处理器410可包括一个或多个处理核心。在本发明中,处理器410可以用于确定门控信号的相关处理。收发器420可以用于接收和发送数据,终端可以通过收发器420接收和发送数据,终端可以通过因特网收发数据,收发器可以是网卡。
存储器430可用于存储软件程序以及模块,处理器410通过运行存储在存储器430的软件程序以及模块,从而执行各种功能应用以及数据处理。存储器430可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序(比如确定门控信号功能等)等;存储数据区可存储根据终端的使用所创建的数据(比如湮灭点位置信息等)等。此外,存储器430可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。输入单元440可以用于接收输入的数字或字符信息,以及产生与用户设置以及功能控制有关的键盘、鼠标、操作杆、光学或者轨迹球信号输入。显示单元450可以用于显示由用户输入的信息或提供给用户的信息以及终端的各种图形用户接口,这些图形用户接口可以由图形、文本、图标、视频和其任意组合来构成。显示单元450可包括显示面板451,可选的,可以采用LCD(LiquidCrystal Display,液晶显示器)、OLED(Organic Light-Emitting Diode,有机发光二极管)等形式来配置显示面板451。音频电路460、扬声器461,传声器462可提供用户与终端之间的音频接口,音频电路460可将接收到的音频数据转换为电信号。电源470可以通过电源管理系统与处理器410逻辑相连,从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。电源470还可以包括一个或一个以上的直流或交流电源、再充电系统、电源故障检测电路、电源转换器或者逆变器、电源状态指示器等任意组件。
具体在本实施例中,处理终端400还包括有存储器,以及一个或者一个以上的程序,其中一个或者一个以上程序存储于存储器中,且经配置以由一个或者一个以上处理器执行。上述一个或者一个以上程序包含用于进行以下操作的指令:
根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数;
构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;
根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵,将所述T个相位图像合成目标相位图像,不同相位图像合成的权重可以是经验数值或是相关相位内采集数据的总数。
可选的,所述构建所述第t个相位图像的转换矩阵,包括:
获取电子计算机断层扫描4D-CT设备构建的所述第t个相位的CT图像和所述目标相位的CT图像;
根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵。
可选的,所述根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵,包括:
根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
可选的,所述根据所述T个相位图像、所述T个相位图像中的每个相位图像的转换矩阵以及在所述T个相位中的每个相位内所述各探测器对的采集数据,将所述T个相位图像合成目标相位图像,包括:
按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
可选的,所述将所述T个相位图像合成目标相位图像之后,还包括:
将所述目标相位图像和所述目标相位的CT图像融合成一幅图像,显示所述融合的一幅图像。
在本发明实施例中,将T个相位的相位图像合成一个目标相位图像,这样利用了全部相位内各探测器对的采集数据来构建一幅目标相位图像,从而降低了目标相位图像的噪声。将该目标相位图像呈现给医生,能够更好地帮助医生诊断疾病。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种构建相位图像的方法,其特征在于,所述方法包括:
根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数;
构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;
按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
2.如权利要求1所述的方法,其特征在于,所述构建所述第t个相位图像的转换矩阵,包括:
获取电子计算机断层扫描4D-CT设备构建的所述第t个相位的CT图像和所述目标相位的CT图像;
根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵。
3.如权利要求2所述的方法,其特征在于,所述根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵,包括:
根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
4.如权利要求1所述的方法,其特征在于,所述将所述T个相位图像合成目标相位图像之后,还包括:
将所述目标相位图像和所述目标相位的CT图像融合成一幅图像,显示所述融合的一幅图像。
5.一种构建相位图像的装置,其特征在于,所述装置包括:
第一构建模块,用于根据第t个相位内接收γ光子对的各探测器对的采集数据,构建所述第t个相位图像,t=0,1,2…T-1,T为相位总数目,其值为大于1的整数;
第二构建模块,用于构建所述第t个相位图像的转换矩阵,所述转换矩阵用于将所述第t个相位图像转换为目标相位的相位图像,所述目标相位是所述T个相位中的一个相位;
合成模块,用于按如下公式将所述T个相位图像合成目标相位图像
在上述公式中,为所述第t个相位图像,Mt->e为所述第t个相位图像的转换矩阵,ymt为在第t相位内所述各探测器对的采集数据,ym为在所有相位中所述各探测器对的采集数据。
6.如权利要求5所述的装置,其特征在于,所述第二构建模块包括:
获取单元,用于获取电子计算机断层扫描4D-CT设备构建的所述第t个相位的CT图像和所述目标相位的CT图像;
构建单元,用于根据所述第t个相位的CT图像和所述目标相位的CT图像,构建所述第t个相位的转换矩阵。
7.如权利要求6所述的装置,其特征在于,所述构建单元,用于根据所述第t个相位的CT图像和所述目标相位的CT图像,通过B样条函数构建所述第t个相位的转换矩阵。
8.如权利要求5所述的装置,其特征在于,所述装置还包括:
融合模块,用于将所述目标相位图像和所述目标相位的CT图像融合成一幅图像,显示所述融合的一幅图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610976174.0A CN106548464B (zh) | 2016-11-07 | 2016-11-07 | 一种构建相位图像的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610976174.0A CN106548464B (zh) | 2016-11-07 | 2016-11-07 | 一种构建相位图像的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106548464A CN106548464A (zh) | 2017-03-29 |
CN106548464B true CN106548464B (zh) | 2019-04-02 |
Family
ID=58394223
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610976174.0A Active CN106548464B (zh) | 2016-11-07 | 2016-11-07 | 一种构建相位图像的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106548464B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107638188B (zh) * | 2017-09-28 | 2021-03-19 | 江苏赛诺格兰医疗科技有限公司 | 图像衰减校正方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101454804A (zh) * | 2006-05-26 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | 动态计算机断层摄影成像 |
CN101856234A (zh) * | 2009-04-08 | 2010-10-13 | 株式会社东芝 | X射线诊断装置以及图像重构处理装置 |
CN102144927A (zh) * | 2010-02-10 | 2011-08-10 | 清华大学 | 基于运动补偿的ct设备和方法 |
CN104331914A (zh) * | 2014-09-18 | 2015-02-04 | 东莞南方医大松山湖科技园有限公司 | 重建心脏相位图像的方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102008058488B4 (de) * | 2008-11-21 | 2018-09-20 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Aufbereitung von kombinierten MR-Emissionstomographieaufnahmen |
-
2016
- 2016-11-07 CN CN201610976174.0A patent/CN106548464B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101454804A (zh) * | 2006-05-26 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | 动态计算机断层摄影成像 |
CN101856234A (zh) * | 2009-04-08 | 2010-10-13 | 株式会社东芝 | X射线诊断装置以及图像重构处理装置 |
CN102144927A (zh) * | 2010-02-10 | 2011-08-10 | 清华大学 | 基于运动补偿的ct设备和方法 |
CN104331914A (zh) * | 2014-09-18 | 2015-02-04 | 东莞南方医大松山湖科技园有限公司 | 重建心脏相位图像的方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106548464A (zh) | 2017-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107638188B (zh) | 图像衰减校正方法及装置 | |
US8897518B2 (en) | Functional imaging | |
US7737406B2 (en) | Compensating for truncated CT images for use as attenuation maps in emission tomography | |
US10311604B2 (en) | System and method for generating attenuation map | |
JP3800101B2 (ja) | 断層像作成装置及び断層像作成方法並びに放射線検査装置 | |
US7729467B2 (en) | Methods and systems for attentuation correction in medical imaging | |
US8658979B2 (en) | Nuclear image reconstruction | |
CN104093360A (zh) | 核成像系统 | |
JP7254656B2 (ja) | 医用画像処理装置、医用画像診断装置及び核医学診断装置 | |
CN107348969A (zh) | 一种pet数据处理方法、系统及pet成像设备 | |
CN107610198A (zh) | Pet图像衰减校正方法及装置 | |
JP2020076584A (ja) | 医用画像処理装置 | |
US20140133707A1 (en) | Motion information estimation method and image generation apparatus using the same | |
JP2009540935A (ja) | 対象物の位置を決定するためのシステムおよび方法 | |
CN106548473B (zh) | 一种构建相位图像的方法及装置 | |
JP2004313785A (ja) | 断層撮影システム及びx線投影システムの組み合わせ装置 | |
CN106548464B (zh) | 一种构建相位图像的方法及装置 | |
CN110215226B (zh) | 图像衰减校正方法、装置、计算机设备和存储介质 | |
JP2007512513A (ja) | 血液中のトレーサの濃度を決定する装置及び方法 | |
CN111316327B (zh) | 对移动对象的pet数据的衰减校正 | |
TWI430777B (zh) | 雙光子發射性斷層掃描系統及方法 | |
Nekolla et al. | Clinical impact of cardiac-gated PET imaging | |
Ay et al. | Pet/Ct | |
De Vasconcelos Caribé | Evaluation of image quality and reconstruction parameters in recent PET-CT and PET-MR systems | |
Manwell | Data-driven patient motion compensation in cardiac positron emission tomography |
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 |