[发明专利]PET-CT扫描设备的运动信号分析方法、系统及计算机可读存储介质

说明书

本发明提供了一种PET‑CT扫描设备的运动信号分析方法、系统及计算机可读存储介质，通过识别区域划分模块对原始数据进行区域分割标定目标识别区域；通过识别区域评估模块获取符合事件的排数信息，计算获取位于目标识别区域的符合事件的时间谱和频率谱，再由识别区域划分模块对在频率谱的位置信息Lsubgt;V/subgt;进行带通滤波后获得运动信号，由此实现高排数的计算机断层扫描设备的运动信号的准确提取。

技术领域

本发明涉及电子计算机断层扫描技术领域，尤其涉及一种PET-CT扫描设备的运动信号分析方法、系统及计算机可读存储介质。

背景技术

呼吸或心脏运动伪影是PET CT扫描中常见的伪影形式。它的出现直接降低了病人肺部、肝部等多个重要部位的特性分析。为了消除这种伪影，一般都需要提取病人扫描时的呼吸状态，再进行门控、配准等处理，进而消除呼吸运动的影响。

由于PET系统的排数的增长，如目前商用的72排或96排机型，机器一次扫描的轴向长度可以到达300mm以上，在扫描范围内，运动部分可能只是轴向的部分区域而非运动区域的信号则会形成噪声。因此，如果使用常规的数据驱动方法进行运动分析，往往造成运动信号淹没在噪声之中，从而无法获取到合适的运动信息。

发明内容

为了克服上述技术缺陷，本发明的目的在于提供一种可以在高排数探测器的数据中准确的提取出运动信号的PET-CT扫描设备的运动信号分析方法、系统及计算机可读存储介质。

本发明公开了一种PET-CT扫描设备的运动信号分析方法，包括如下步骤：获取符合事件的探测器排数；对于一待评估扫描过程，预设若干个连续的时间范围阈值，并以探测器排号范围标定的轴向范围为目标识别区域，获取每个所述时间范围阈值的所有目标识别区域；统计在每个所述时间范围阈值内落入所述目标识别区域的符合事件，记为识别符合事件；获取所述识别符合事件的特征参数；对所有所述时间范围阈值内的所述特征参数进行频域变换；设定一运动目标频率范围，获取所述特征参数在所述运动目标频率范围内的频谱功率的最大值V_max和平均值V_AVE，最大值V_max和平均值V_AVE的比值作为运动信号的评价参数S；遍历所有目标识别区域内的评价参数S，选取最大值S_max获取其位置信息L_V、时间谱和频率谱；对最大值S_max对应的频率谱的位置信息L_V进行带通滤波后进行时域变换获得运动信号。

优选地，所述以探测器排号范围标定的轴向范围为目标识别区域包括：遍历所述待评估扫描过程中所有可能的排数识别长度L_d，0＜L_d＜L；遍历所述待评估扫描过程中所有可能的偏移排数L_off，0＜L_off＜L-L_d；则所有L_off和L_d的组合范围记为所述目标识别区域。

优选地，所述获取符合事件的探测器排数包括：从原始扫描数据中提取符合事件，获取每一个符合事件的探测器晶体序号，分析所述探测器晶体序号所位于的探测器排号；若一符合事件的起始排号为W₀，结束排号为W₁，则所述符合事件的排数为L＝W₁-W₀+1。

优选地，所述特征参数包括每个所述时间范围阈值内所有所述识别符合事件的轴向平均位置或者所有所述识别符合事件的总数。