一种确定门控信号的方法和装置的制造方法

文档序号:10726471阅读:470来源:国知局
一种确定门控信号的方法和装置的制造方法
【专利摘要】本发明公开了一种确定门控信号的方法和装置,属于医学领域。所述方法包括:获取正电子发射断层显像PET设备在检测时段内的每个预设时长的子时段中检测到的各湮灭点位置信息;根据所述每个子时段对应的各湮灭点位置信息,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息;根据在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,确定门控信号,基于所述门控信号进行PET图像重建。采用本发明,可以降低诊断系统成本。
【专利说明】
_种确定门控信号的方法和装置
技术领域
[0001] 本发明涉及医学领域,特别涉及一种确定门控信号的方法和装置。
【背景技术】
[0002] 在医学领域中,可以通过PET/CT(positron emission tomography/computed tomography,正电子发射断层显像/计算机断层扫描)设备对人体或动物体肿瘤、心脏系统 疾病和神经系统疾病进行早期诊断,即可以通过PET/CT设备对被扫描对象扫描得到的图像 确定发生病灶的位置。然而,在PET成像时,由于成像过程具有一定的时长,被扫描对象的呼 吸运动或心脏搏动使被发生病灶的组织(比如肺)在空间上周期性移位,导致所得的PET图 像模糊,使得难以或不可能识别PET图像中的小病变。
[0003] 为消除呼吸运动导致的PET图像模糊,往往在PET中采用门控技术,即通过特定的 设备采集呼吸运动对应的门控信号(其中,门控信号可以是用于反映被扫描对象在扫描时 长(比如12分钟)内的呼吸运动过程或心脏搏动过程的信号),目前,常用的采集门控信号的 方法往往需要依靠外部设备,比如,通过红外线采集装置或摄像机采集被扫描对象体表的 运动信号,其中,运动信号可以是被扫描对象处于平躺状态下,体表某部位的高度随着时间 进行变化的高度值,采集到的该运动信号即为门控信号。然后,可以根据得到的门控信号, 进行PET图像重建,即可以根据门控信号反映出的呼吸周期,将很多次重复的相同运动阶段 中的PET数据划分在一起,将划分在一起的PET数据重建为一个运动阶段的图像。
[0004] 在实现本发明的过程中,发明人发现现有技术至少存在以下问题:
[0005] 基于上述处理方式采集门控信号时,需要增加额外的外部设备,而外部设备的价 格往往比较昂贵,从而,导致诊断系统成本增加。

【发明内容】

[0006] 为了解决现有技术的问题,本发明实施例提供了一种确定门控信号的方法和装 置。所述技术方案如下:
[0007] 第一方面,提供了一种确定门控信号的方法,所述方法包括:
[0008] 获取正电子发射断层显像PET设备在检测时段内的每个预设时长的子时段中检测 到的各湮灭点位置信息;
[0009] 根据所述每个子时段对应的各湮灭点位置信息,确定病灶组织的运动方向,确定 在每个子时段所述病灶组织在所述运动方向上的运动幅度信息;
[0010] 根据在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,确定门控信 号,基于所述门控信号进行PET图像重建。
[0011] 可选的,所述根据所述每个子时段对应的各湮灭点位置信息,确定病灶组织的运 动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括:
[0012] 对于每个子时段,根据所述子时段对应的各湮灭点位置信息,确定所述子时段中 不同湮灭点位置信息对应的湮灭点数目;
[0013] 根据所述每个子时段中不同湮灭点位置信息对应的湮灭点数目,确定病灶组织的 运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息。
[0014] 可选的,所述对于每个子时段,根据所述子时段对应的各湮灭点位置信息,确定所 述子时段中不同湮灭点位置信息对应的湮灭点数目,包括:
[0015] 对于每个子时段,根据所述子时段对应的各湮灭点位置信息中的径向位置、角度 位置和轴向位置标识,确定每个径向位置-角度位置-轴向位置标识组合对应的湮灭点数 目,其中,所述轴向位置标识是用于表示湮灭点的轴向位置的信息;
[0016] 所述根据所述每个子时段中不同湮灭点位置信息对应的湮灭点数目,确定病灶组 织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括:
[0017] 对于每个子时段,将所述子时段中所述每个径向位置-角度位置-轴向位置标识组 合中满足预设接近度条件的多个径向位置-角度位置组合合并为一个新的径向位置-角度 位置组合,并将多个径向位置-角度位置组合对应的湮灭点数目之和,作为所述新的径向位 置-角度位置组合对应的湮灭点数目,得到所述子时段中新的径向位置-角度位置-轴向位 置标识组合对应的湮灭点数目;
[0018] 根据所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合对应的 湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向 上的运动幅度信息。
[0019] 可选的,所述根据所述每个子时段中所述新的径向位置-角度位置-轴向位置标识 组合对应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所 述运动方向上的运动幅度信息,包括:
[0020] 对所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合进行数据 重组处理,得到数据重组后的径向位置-角度位置-轴向位置标识组合,根据所述每个子时 段中所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子 时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目;
[0021] 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组 合对应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述 运动方向上的运动幅度信息。
[0022] 可选的,所述根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向 位置标识组合对应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶 组织在所述运动方向上的运动幅度信息,包括:
[0023] 根据核素衰变公式和所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子时段中所述数据重组后的径向位 置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目;
[0024] 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病 灶组织在所述运动方向上的运动幅度信息。
[0025] 可选的,所述根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向 位置标识组合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子 时段所述病灶组织在所述运动方向上的运动幅度信息,包括:
[0026] 对于所述每个子时段中每个所述数据重组后的径向位置-角度位置-轴向位置标 识组合中的径向位置-角度位置组合,根据所述径向位置-角度位置组合所属的数据重组后 的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定所述径向 位置-角度位置组合对应的衰变修正后的湮灭点数目;
[0027] 将每个子时段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度 位置组合,所属的数据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮 灭点数目设置为零,得到所述每个子时段中消除噪声后的径向位置-角度位置-轴向位置标 识组合对应的衰变修正后的湮灭点数目;
[0028] 根据所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病 灶组织在所述运动方向上的运动幅度信息。
[0029] 可选的,所述根据所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向 位置标识组合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子 时段所述病灶组织在所述运动方向上的运动幅度信息,包括:
[0030] 根据预设的主成分分析PCA算法,对由所述每个子时段中所述消除噪声后的径向 位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目构成的矩阵进行主成 分分析,得到至少一个特征值,以及所述至少一个特征值对应的特征向量,其中,所述特征 值越大表示湮灭点数目在所述特征值的特征向量对应的方向上变化程度越大;
[0031] 选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将所述 第一特征向量或所述第二特征向量对应的方向作为病灶组织的运动方向;
[0032] 对于每个子时段,根据所述子时段中所述消除噪声后的径向位置-角度位置-轴向 位置标识组合对应的衰变修正后的湮灭点数目,以及所述第一特征向量或所述第二特征向 量,确定在所述子时段所述病灶组织在所述运动方向上的运动幅度信息。
[0033] 可选的,所述根据在每个子时段所述病灶组织在所述运动方向上的运动幅度信 息,确定门控信号,包括:
[0034] 对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度 信息构成的运动幅度信息-时间函数进行插值处理,得到插值处理后的运动幅度信息-时间 函数;
[0035]将所述插值处理后的运动幅度信息-时间函数作为门控信号。
[0036] 可选的,所述根据在每个子时段所述病灶组织在所述运动方向上的运动幅度信 息,确定门控信号,包括:
[0037] 对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度 信息构成的运动幅度信息-时间函数进行频域变换处理,得到所述运动幅度信息-时间函数 对应的频谱;
[0038] 基于具有预设中心频率和预设带宽的带通滤波器,对所述频谱进行滤波,得到滤 波后的频谱;
[0039]对所述滤波后的频谱进行所述频域变换的逆变换,得到门控信号。
[0040] 第二方面,提供了一种确定门控信号的装置,所述装置包括:
[0041] 获取模块,用于获取正电子发射断层显像PET设备在检测时段内的每个预设时长 的子时段中检测到的各湮灭点位置信息;
[0042] 第一确定模块,用于根据所述每个子时段对应的各湮灭点位置信息,确定病灶组 织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息;
[0043] 第二确定模块,用于根据在每个子时段所述病灶组织在所述运动方向上的运动幅 度信息,确定门控信号,基于所述门控信号进行PET图像重建。
[0044] 可选的,所述第一确定模块,包括:
[0045] 第一确定子模块,用于对于每个子时段,根据所述子时段对应的各湮灭点位置信 息,确定所述子时段中不同湮灭点位置信息对应的湮灭点数目;
[0046] 第二确定子模块,用于根据所述每个子时段中不同湮灭点位置信息对应的湮灭点 数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运 动幅度信息。
[0047]可选的,所述第一确定子模块,用于:
[0048] 对于每个子时段,根据所述子时段对应的各湮灭点位置信息中的径向位置、角度 位置和轴向位置标识,确定每个径向位置-角度位置-轴向位置标识组合对应的湮灭点数 目,其中,所述轴向位置标识是用于表示湮灭点的轴向位置的信息;
[0049] 所述第二确定子模块,用于:
[0050] 对于每个子时段,将所述子时段中所述每个径向位置-角度位置-轴向位置标识组 合中满足预设接近度条件的多个径向位置-角度位置组合合并为一个新的径向位置-角度 位置组合,并将多个径向位置-角度位置组合对应的湮灭点数目之和,作为所述新的径向位 置-角度位置组合对应的湮灭点数目,得到所述子时段中新的径向位置-角度位置-轴向位 置标识组合对应的湮灭点数目;
[0051] 根据所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合对应的 湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向 上的运动幅度信息。
[0052] 可选的,所述第二确定子模块,用于:
[0053]对所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合进行数据 重组处理,得到数据重组后的径向位置-角度位置-轴向位置标识组合,根据所述每个子时 段中所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子 时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目;
[0054] 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组 合对应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述 运动方向上的运动幅度信息。
[0055] 可选的,所述第二确定子模块,用于:
[0056] 根据核素衰变公式和所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子时段中所述数据重组后的径向位 置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目;
[0057] 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病 灶组织在所述运动方向上的运动幅度信息。
[0058] 可选的,所述第二确定子模块,用于:
[0059] 对于所述每个子时段中每个所述数据重组后的径向位置-角度位置-轴向位置标 识组合中的径向位置-角度位置组合,根据所述径向位置-角度位置组合所属的数据重组后 的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定所述径向 位置-角度位置组合对应的衰变修正后的湮灭点数目;
[0060] 将每个子时段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度 位置组合,所属的数据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮 灭点数目设置为零,得到所述每个子时段中消除噪声后的径向位置-角度位置-轴向位置标 识组合对应的衰变修正后的湮灭点数目;
[0061] 根据所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病 灶组织在所述运动方向上的运动幅度信息。
[0062] 可选的,所述第二确定子模块,用于:
[0063] 根据预设的主成分分析PCA算法,对由所述每个子时段中所述消除噪声后的径向 位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目构成的矩阵进行主成 分分析,得到至少一个特征值,以及所述至少一个特征值对应的特征向量,其中,所述特征 值越大表示湮灭点数目在所述特征值的特征向量对应的方向上变化程度越大;
[0064] 选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将所述 第一特征向量或所述第二特征向量对应的方向作为病灶组织的运动方向;
[0065] 对于每个子时段,根据所述子时段中所述消除噪声后的径向位置-角度位置-轴向 位置标识组合对应的衰变修正后的湮灭点数目,以及所述第一特征向量或所述第二特征向 量,确定在所述子时段所述病灶组织在所述运动方向上的运动幅度信息。
[0066]可选的,所述第二确定模块,用于:
[0067]对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度 信息构成的运动幅度信息-时间函数进行插值处理,得到插值处理后的运动幅度信息-时间 函数;
[0068]将所述插值处理后的运动幅度信息-时间函数作为门控信号。
[0069]可选的,所述第二确定模块,包括:
[0070] 频域变换子模块,用于对所述每个子时段和在每个子时段所述病灶组织在所述运 动方向上的运动幅度信息构成的运动幅度信息-时间函数进行频域变换处理,得到所述运 动幅度信息-时间函数对应的频谱;
[0071] 滤波子模块,用于基于具有预设中心频率和预设带宽的带通滤波器,对所述频谱 进行滤波,得到滤波后的频谱;
[0072]频域逆变换子模块,用于对所述滤波后的频谱进行所述频域变换的逆变换,得到 门控信号。
[0073] 本发明实施例提供的技术方案带来的有益效果是:
[0074] 本发明实施例中,可以根据PET设备检测到的各湮灭点位置信息确定门控信号,无 需在PET设备检测各湮灭点位置信息的时候,增加额外的外部设备采集门控信号,从而,可 以降低诊断系统成本。
【附图说明】
[0075] 为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使 用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于 本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他 的附图。
[0076] 图1是本发明实施例提供的一种确定门控信号的方法流程图;
[0077] 图2是本发明实施例提供的一种终端的结构示意图;
[0078] 图3是本发明实施例提供的一种探测晶体对的示意图;
[0079] 图4是本发明实施例提供的一种确定门控信号的装置结构示意图;
[0080] 图5是本发明实施例提供的一种确定门控信号的装置结构示意图;
[0081] 图6是本发明实施例提供的一种确定门控信号的装置结构示意图。
【具体实施方式】
[0082]为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方 式作进一步地详细描述。
[0083]本发明实施例提供了一种确定门控信号的方法,如图1所示,该方法的执行主体为 终端。其中,该终端可以是具有确定门控信号功能的终端,比如,可以是与PET (posi tron emi ss ion tomography,正电子发射断层显像)设备相连,可以根据PET设备采集的数据确定 门控信号的计算机设备。
[0084] 终端可以包括处理器210、收发器220、存储器230、输入单元240、显示单元250、音 频电路260以及电源270等部件,如图2所示,本领域技术人员可以理解,图2中示出的终端结 构并不构成对终端的限定,可以包括比图示更多或更少的部件或者组合某些部件,或者不 同的部件布置。其中:
[0085]处理器210可以是终端的控制中心,利用各种接口和线路连接整个终端设备的各 个部分,如收发器220和存储器230等,通过运行或执行存储在存储器230内的软件程序和/ 或模块,以及调用存储在存储器230内的数据,执行终端的各种功能和处理数据,从而对终 端进行整体监控。可选的,处理器210可包括一个或多个处理核心。在本发明中,处理器210 可以用于确定门控信号的相关处理。收发器220可以用于接收和发送数据,终端可以通过收 发器220接收和发送数据,终端可以通过因特网收发数据,收发器可以是网卡。
[0086] 存储器230可用于存储软件程序以及模块,处理器210通过运行存储在存储器230 的软件程序以及模块,从而执行各种功能应用以及数据处理。存储器230可主要包括存储程 序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序(比 如确定门控信号功能等)等;存储数据区可存储根据终端的使用所创建的数据(比如湮灭点 位置信息等)等。此外,存储器230可以包括高速随机存取存储器,还可以包括非易失性存储 器,例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。输入单元240可以 用于接收输入的数字或字符信息,以及产生与用户设置以及功能控制有关的键盘、鼠标、操 作杆、光学或者轨迹球信号输入。显示单元250可以用于显示由用户输入的信息或提供给用 户的信息以及终端的各种图形用户接口,这些图形用户接口可以由图形、文本、图标、视频 和其任意组合来构成。显示单元250可包括显示面板251,可选的,可以采用LCD(Liquid Crystal Display,液晶显示器)、OLED(Organic Light-Emitting Diode,有机发光二极管) 等形式来配置显示面板251。音频电路260、扬声器261,传声器262可提供用户与终端之间的 音频接口,音频电路260可将接收到的音频数据转换为电信号。电源270可以通过电源管理 系统与处理器210逻辑相连,从而通过电源管理系统实现管理充电、放电、以及功耗管理等 功能。电源270还可以包括一个或一个以上的直流或交流电源、再充电系统、电源故障检测 电路、电源转换器或者逆变器、电源状态指示器等任意组件。终端还包括有存储器,以及一 个或者一个以上的程序,其中一个或者一个以上程序存储于存储器中,且经配置以由一个 或者一个以上处理器执行述一个或者一个以上程序来执行下述各个实施例所述的确定门 控信号的方法。
[0087]下面将结合【具体实施方式】,对图1所示的处理流程进行详细的说明,内容可以如 下:
[0088]步骤101,获取正电子发射断层显像PET设备在检测时段内的每个预设时长的子时 段中检测到的各湮灭点位置信息。
[0089] 其中,湮灭点可以是正电子与体内的负电子发生湮灭产生一对光子的位置点,湮 灭点位置信息可以是用于表示湮灭点位置的信息,比如,可以是几近同时接收到一对光子 的两个探测晶体的编号和轴向位置标识,其中,轴向位置标识可以是用于表示湮灭点的轴 向位置的信息,轴向位置标识可以用PET设备中每个探测器环的编号表示,例如,每个探测 器环包括60个探测晶体,编号依次为1,2,…60,PET设备共包括15个探测器环,编号依次是 1,2,···15,某个湮灭点位置信息可以是(40,50,2),即表示在该湮灭点产生的一对光子,被 第2个探测器环的40号和50号探测晶体几近同时检测到。
[0090] 湮灭点位置信息还可以是径向位置、角度位置和轴向位置标识,即由径向位置、角 度位置和轴向位置标识三个参数组合作为一个湮灭点位置信息,其中,由探测器环编号和 轴向位置标识组成的湮灭点位置信息与由径向位置、角度位置和轴向位置标识组成的湮灭 点位置信息具有一一对应关系,探测器环编号与径向位置-角度位置组合之间可以互相转 换,角度位置可以是两个探测器环编号对应的探测晶体之间的连线与竖直线之间的夹角, 取值在0-180之间,径向位置可以是探测器环的圆心到两个探测晶体之间的连线的距离。
[0091] 需要说明的一点是,PET设备检测到的湮灭点位置信息中的轴向位置标识可以用 探测器环的编号表示。具体来讲,当几近同时检测到一对光子的探测晶体位于同一个探测 器环时,该湮灭点对应的轴向位置标识即是该探测器环的编号。但是,几近同时检测到一对 光子的探测晶体还可能位于不同的探测器环,如图3所示,此时,该湮灭点对应的轴向位置 标识可以用探测晶体位于的两个探测器环的编号组合表示,也可以用预先设置的这两个探 测器环组合对应的标识表示,其中,本地中可以预先设置有每两个探测器环组合对应的轴 向位置标识,可以是大于PET设备所包含的探测器环数目的数字,例如,PET设备共有15个探 测器环,编号依次是1-15,可以预先设置16表示两个探测晶体位于第1、2探测器环、17表示 两个探测晶体位于第1、3探测器环、18表示两个探测晶体位于第1、4探测器环,以此类推,每 两个探测器环组合都预先对应有一个编号。这样,PET设备检测到湮灭点时,可以根据检测 到该湮灭点的探测晶体所在的探测器环,为该湮灭点分配相应的轴向位置标识。
[0092]在实施中,在医学领域,可以通过PET设备对人体或动物体肿瘤、心脏系统疾病和 神经系统疾病进行早期诊断,即可以通过PET设备对被扫描对象扫描得到的图像确定发生 病灶的位置。具体的,首先可以对被扫描对象注射放射性核素(如18F、nC),其中,放射性核 素是正电子核素,放射性核素会随着代谢进入体内的组织(其中,发生病灶的组织会吸收大 量的放射性核素),然后,正电子将会与体内的负电子发生湮灭,产生一对反向运动的光子, 这一对反向运动的光子将会被位于对称位置的两个探测晶体检测到,其中,PET设备每个探 测器环上具有相同个数的探测晶体,且每个探测器环上的每个探测晶体具有唯一的编号, PET设备每检测到一对光子时,即可对表示该湮灭点位置的信息(即湮灭点位置信息)进行 存储,即可以对几近同时检测到一对光子的两个探测晶体的编号和轴向位置标识进行存 储,最后,终端可以根据PET设备检测到的湮灭点位置信息,形成PET图像。进而,医生可以根 据PET图像确定该扫描对象是否有发生病灶的组织,其中,PET设备在检测各湮灭点位置信 息时,需要一定的检测时长(一般检测时长大约在12分钟左右),在PET设备对被扫描对象扫 描期间,被扫描对象会进行多次呼吸运动或心脏搏动,此种情况下,若被扫描对象有发生病 灶的组织时,由于病灶组织的运动,将会导致得到的PET图像发生模糊。
[0093]目前,可以通过门控技术解决由于病灶组织运动导致得到的PET图像模糊的问题, 下面将具体介绍终端确定门控信号的详细过程。
[0094] 终端可以将PET设备的整个检测时段划分为多个具有相同预设时长的子时段,其 中,预设时长可以是0.5秒,也可以是1秒。终端可以获取PET设备在检测时段内检测到的多 个湮灭点位置信息,其中,PET设备检测到的各湮灭点位置信息可以是按照各湮灭点位置信 息对应的检测时间依次排列的列表式数据,进而,终端可以基于各湮灭点位置信息对应的 检测时间,确定在检测时段内的每个预设时长的子时段中检测到的各湮灭点位置信息。
[0095] 步骤102,根据每个子时段对应的各湮灭点位置信息,确定病灶组织的运动方向, 确定在每个子时段病灶组织在运动方向上的运动幅度信息。
[0096] 其中,运动幅度信息可以是反映病灶组织在运动方向上的运动幅度的信息。
[0097] 在实施中,对于每个子时段,终端可以根据每个子时段对应的各湮灭点位置信息, 确定病灶组织的运动方向,以及根据每个子时段对应的各湮灭点位置信息,确定病灶组织 在运动方向上的运动幅度信息随时间的变化趋势,即确定病灶组织在每个子时段的运动幅 度信息。也就是说,当湮灭点位置信息包含探测晶体编号和轴向位置标识时,可以根据各湮 灭点对应的探测晶体编号和轴向位置标识,确定病灶组织的运动方向,以及确定在每个子 时段病灶组织在运动方向上的运动幅度信息。当湮灭点位置信息包含径向位置、角度位置 和轴向位置标识时,可以根据各湮灭点对应的径向位置、角度位置和轴向位置标识,确定病 灶组织的运动方向,以及确定在每个子时段病灶组织在运动方向上的运动幅度信息。
[0098] 可选的,可以根据每个子时段中的不同的湮灭点位置信息的出现次数(湮灭点数 目),确定病灶组织的运动方向,以及确定在每个子时段所述病灶组织在所述运动方向上的 运动幅度信息,相应的,步骤102的处理过程可以如下:对于每个子时段,根据子时段对应的 各湮灭点位置信息,确定子时段中不同湮灭点位置信息对应的湮灭点数目;根据每个子时 段中不同湮灭点位置信息对应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时 段所述病灶组织在运动方向上的运动幅度信息。
[0099] 在实施中,终端获取到每个子时段对应的各湮灭点位置信息(其中,每个子时段对 应的各湮灭点位置信息中包含重复的湮灭点位置信息)后,可以对每个不同的湮灭点位置 信息的出现次数进行统计,获取每个不同的湮灭点位置信息的出现次数(即每个子时段中 不同湮灭点位置信息对应的湮灭点数目),即可得到每个子时段对应的正弦图Sinogram数 据,其中,Sinogram数据可以包括每个不同的湮灭点位置信息及其对应的出现次数(湮灭点 数目),Sinogram数据可以用来反映湮灭点的分布情况。例如,第一子时段的湮灭点位置信 息包括(1,10,1)、(5,10,1)、(1,10,3)、(10,20,1)、(1,10,3)、(10,20,1)、(5,10,1)、(10, 20,1)、(1,10,3)、(1,10,3)、(1,10,2)、(1,10,1)、(1,10,2)、(5,10,1)、(1,10,1)、(5,10, 1),对该子时段的湮灭点位置信息的出现次数进行统计,得到的Sinogram数据包含的信息 为:3个(1,10,1)、4个(5,10,1)、4个(1,10,3)、3个(10,20,1)、2个(1,10,2),其中^11(^瓜111 数据可以是以表格的形式记录每个湮灭点位置信息及其对应的湮灭点数目,需要说明的 是,在实际操作中,每个子时段会包括大量的湮灭点位置信息,上述示例只是为了说明由列 表式数据到Sinogram数据的处理过程。终端得到每个子时段的Sinogram数据后,可以根据 每个子时段的Sinogram数据,确定病灶组织的运动方向,确定在每个子时段所述病灶组织 在运动方向上的运动幅度信息。
[0100] 可选的,终端可以根据每个子时段对应的各湮灭点位置信息中的径向位置、角度 位置和轴向位置标识,确定不同湮灭点位置信息对应的湮灭点数目,相应的,处理过程可以 如下:对于每个子时段,根据子时段对应的各湮灭点位置信息中的径向位置、角度位置和轴 向位置标识,确定每个径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,其中,轴 向位置标识是用于表示湮灭点的轴向位置的信息。相应的,处理过程可以如下:对于每个子 时段,将子时段中每个径向位置-角度位置-轴向位置标识组合中满足预设接近度条件的多 个径向位置-角度位置组合合并为一个新的径向位置-角度位置组合,并将多个径向位置-角度位置组合对应的湮灭点数目之和,作为新的径向位置-角度位置组合对应的湮灭点数 目,得到子时段中新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目;根据每 个子时段中新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定病灶组织 的运动方向,确定在每个子时段病灶组织在所述运动方向上的运动幅度信息。
[0101] 在实施中,终端得到每个子时段对应的湮灭点位置信息后,其中,湮灭点位置信息 可以是径向位置、角度位置和轴向位置标识,可以统计每个径向位置-角度位置-轴向位置 标识组合的出现次数,即可以统计每个径向位置-角度位置-轴向位置标识组合对应的湮灭 点数目,也即可以根据每个子时段对应的各湮灭点位置信息中的径向位置、角度位置和轴 向位置标识,确定每个子时段对应的Sinogram数据,其中,Sinogram数据的格式可以是预先 存储的,每确定一个径向位置-角度位置-轴向位置标识组合对应的湮灭点数目后,可以将 该湮灭点数目添加到Sinogram数据中的相应位置,即每个子时段的Sinogram数据中包含的 各个径向位置-角度位置-轴向位置标识组合是固定的,每个径向位置-角度位置-轴向位置 标识组合对应的湮灭点数目是变化的,是根据每个子时段的各湮灭点位置信息确定出的, Sinogram数据可以是三维表格的格式。例如,对于某轴向位置标识对应的子Sinogram数据, 如表1所示,其中,按照从左到右,从上到下数值依次增大的顺序排列,每个径向位置和角度 位置组合对应的数目即为该湮灭点位置信息对应的湮灭点数目。例如,对于轴向位置标识 18,位置(0,0)对应的5即表示径向位置为0、角度位置为0、轴向位置标识为18的湮灭点位置 信息对应的湮灭点数目为5。
[0102] 表1
[0105] 终端得到每个子时段对应的Sinogram数据后,对于每个轴向位置标识对应的子 Sinogram数据,可以使子Sinogram数据按照预设的接近度条件,沿着径向位置方向和角度 位置方向进行压缩。具体的,终端可以将每相邻(即满足预设接近条件)的预设数目个径向 位置-角度位置组合合并为一个新的径向位置-角度位置组合,其中,新的径向位置-角度位 置组合中的径向位置可以是多个径向位置-角度位置组合中的径向位置的均值,新的径向 位置-角度位置组合中的角度位置可以是多个径向位置-角度位置组合中的角度位置的均 值,并且可以将预设数目个径向位置-角度位置组合对应的湮灭点数目进行相加,得到的湮 灭点数目之和作为新的径向位置-角度位置组合对应的湮灭点数目,对于每个轴向位置标 识,终端都可按照上述数据压缩方法沿着径向位置方向和角度位置方向进行数据压缩,得 到子时段中新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目。例如,如表1所 示的子Sinogram数据,预设数目可以是4,终端可以将(0,0)、(1,0)、(0,5)、(1,5)合并为一 个新的径向位置-角度位置组合_(0.5,2.5),将(2,0)、(3,0)、(2,5)、(3,5)合并为一个新的 径向位置-角度位置组合-(2.5,2.5),将(0,10)、(1,10)、(0,12)、(1,12)合并为一个新的径 向位置-角度位置组合_(0.5,11),将(2,10)、(3,10)、(2,12)、(3,12)合并为一个新的径向 位置-角度位置组合-(2.5,11)。这样,可以使数据压缩后的Sinogram数据中的每个径向位 置-角度位置-轴向位置标识组合对应的湮灭点数目增大,即可以提高每个径向位置-角度 位置-轴向位置标识组合的数据统计量,进而,可以降低噪声的干扰。
[0106] 终端对每个子时段的Sinogram数据进行数据压缩后,可以根据每个子时段中新的 径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定病灶组织的运动方向,确 定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息。具体的,可以根据预设 的主成分分析PCA算法,对由每个子时段中新的径向位置-角度位置-轴向位置标识组合对 应的湮灭点数目构成的矩阵进行主成分分析,得到至少一个特征值,以及至少一个特征值 对应的特征向量,其中,特征值越大表示湮灭点数目在特征值的特征向量对应的方向上变 化程度越大;选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将 第一特征向量或第二特征向量对应的方向作为病灶组织的运动方向;对于每个子时段,根 据子时段中新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目以及第一特征 向量或第二特征向量,确定在子时段病灶组织在运动方向上的运动幅度信息。
[0107] 可选的,终端还可以对每个子时段中新的径向位置-角度位置-轴向位置标识组合 进行数据重组处理,相应的,处理过程可以如下:对每个子时段中新的径向位置-角度位置-轴向位置标识组合进行数据重组处理,得到数据重组后的径向位置-角度位置-轴向位置标 识组合,根据每个子时段中新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数 目,确定每个子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点 数目;根据每个子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭 点数目,确定病灶组织的运动方向,确定在每个子时段病灶组织在所述运动方向上的运动 幅度信息。
[0108] 在实施中,对于每个子时段,终端对子时段的Sinogram数据进行数据压缩后,首先 可以根据预设的数据重组算法(比如SSRB(Single Slice Rebin,单层重组)算法、FORE (Fourier Rebin,傅里叶重组)算法),为湮灭点位置信息表示几近同时接收到光子对的探 测器晶体位于不同探测器环的湮灭点(即轴向位置标识大于探测器环数目的湮灭点位置信 息对应的湮灭点,或者轴向位置标识用两个探测器环编号组合表示的湮灭点位置信息对应 的湮灭点),确定数据重组后的径向位置-角度位置-轴向位置标识组合,也就是说,可以根 据预设的数据重组算法,确定这些湮灭点对应的数据重组后的湮灭点位置信息,例如,PET 设备共包括15个探测器环,某湮灭点对应的湮灭点位置信息为(1,10,17),其中,预先设置 的轴向位置标识17表示两个不同的探测器环为1、3,可以根据预设的数据重组算法,确定该 湮灭点对应的新的轴向位置标识为2,即经过数据重组处理后,该湮灭点的湮灭点位置信息 为(10,20,2),即将大于探测器环数的轴向位置标识转换为小于探测器环数的轴向位置标 识。然后,终端可以对数据重组的各径向位置-角度位置-轴向位置标识组合和未经数据重 组的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目进行重新整合,确定每个子 时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,即得到数 据重组后的Sinogram数据。
[0109] 终端得到数据重组后的Sinogram数据后,可以根据每个子时段中数据重组后的径 向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定病灶组织的运动方向,确定 在每个子时段病灶组织在运动方向上的运动幅度信息。具体的,可以根据预设的主成分分 析PCA算法,对由每个子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应 的湮灭点数目构成的矩阵进行主成分分析,得到至少一个特征值,以及至少一个特征值对 应的特征向量,其中,特征值越大表示湮灭点数目在特征值的特征向量对应的方向上变化 程度越大;选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将第 一特征向量或第二特征向量对应的方向作为病灶组织的运动方向;对于每个子时段,根据 子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目以及第 一特征向量或第二特征向量,确定在子时段病灶组织在运动方向上的运动幅度信息。
[0110] 可选的,考虑到核素衰变对统计数据的影响,得到每个子时段中数据重组后的径 向位置-角度位置-轴向位置标识组合对应的湮灭点数目后,可以对每个子时段对应的各湮 灭点数目进行衰变修正,相应的,处理过程可以如下:根据核素衰变公式和每个子时段中数 据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定每个子时段中 数据重组后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目;根 据每个子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后 的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段病灶组织在所述运动方向上 的运动幅度信息。
[0111] 在实施中,对于每个子时段,得到该子时段下数据重组后的径向位置-角度位置- 轴向位置标识组合对应的湮灭点数目后,可以根据核素衰变公式,对该子时段下数据重组 后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目进行衰变修正,得到该子时 段下数据重组后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数 目,具体的,可以按照公式(1)得到每个子时段中数据重组后的径向位置-角度位置-轴向位 置标识组合对应的衰变修正后的湮灭点数目,
[0112] 5?/)-5(〇·^*'"(2)/?1/2 (1)
[0113] 其中,Sit)为子时段t下数据重组后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目,S(t)为子时段t下数据重组后的径向位置-角度位置-轴向 位置标识组合对应的衰变修正前的湮灭点数目,Tl/2为核素的半衰期,为核素 衰变公式。
[0114] 终端得到每个子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目后,可以根据每个子时段中数据重组后的径向位置-角度位 置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在 每个子时段病灶组织在所述运动方向上的运动幅度信息。具体的,可以根据预设的主成分 分析PCA算法,对由每个子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目构成的矩阵进行主成分分析,得到至少一个特征值,以及至 少一个特征值对应的特征向量,其中,特征值越大表示湮灭点数目在特征值的特征向量对 应的方向上变化程度越大;选取最大特征值对应的第一特征向量或次大特征值对应的第二 特征向量,将第一特征向量或第二特征向量对应的方向作为病灶组织的运动方向;对于每 个子时段,根据子时段中数据重组后的径向位置-角度位置-轴向位置标识组合对应的衰变 修正后的湮灭点数目以及第一特征向量或第二特征向量,确定在子时段病灶组织在运动方 向上的运动幅度信息。
[0115] 可选的,考虑到噪声对确定运动方向的影响,得到每个子时段中数据重组后的径 向位置-角度位置-轴向位置标识组合对应的湮灭点数目后,可以对每个子时段对应的各湮 灭点数目进行消噪处理,相应的,处理过程可以如下:对于每个子时段中每个数据重组后的 径向位置-角度位置-轴向位置标识组合中的径向位置-角度位置组合,根据径向位置-角度 位置组合所属的数据重组后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后 的湮灭点数目,确定径向位置-角度位置组合对应的衰变修正后的湮灭点数目;将每个子时 段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度位置组合,所属的数 据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮灭点数目设置为零, 得到每个子时段中消除噪声后的径向位置-角度位置-轴向位置标识组合对应的衰变修正 后的湮灭点数目;根据每个子时段中消除噪声后的径向位置-角度位置-轴向位置标识组合 对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段病灶组织 在运动方向上的运动幅度信息。
[0116] 在实施中,终端中可以预先存储有数目阈值,其中,每个子时间段可以对应有相应 的数目阈值,具体的,终端中可以预先存储有百分比,每个子时段对应的数目阈值即是该子 时段中的最大湮灭点数目与预设百分比的乘积。PET设备采集的湮灭点位置信息可能包含 噪声信息,考虑到噪声信息对确定运动方向的影响,可以对每个子时段中每个数据重组后 的径向位置-角度位置-轴向位置标识组合对应的衰减修正后的湮灭点数目进行消噪处理, 具体的,首先对于每个子时段中每个数据重组后的径向位置-角度位置-轴向位置标识组合 中的径向位置-角度位置组合,可以对径向位置-角度位置组合所属的数据重组后的径向位 置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目进行相加,将其作为径向 位置-角度位置组合对应的衰变修正后的湮灭点数目,然后,对于每个子时段,判断径向位 置-角度位置组合对应的衰变修正后的湮灭点数目与预设数目阈值的大小关系,确定该子 时段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度位置组合,进而,可 以确定衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度位置组合所属的数 据重组后的径向位置-角度位置-轴向位置标识组合,可以将其所属的数据重组后的径向位 置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目设置为零,衰变修正后的 湮灭点数目超过预设数目阈值的径向位置-角度位置组合,所属的数据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮灭点数目保持不变,得到每个子时段中消 除噪声后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目。最 后,终端可以根据每个子时段中消除噪声后的径向位置-角度位置-轴向位置标识组合对应 的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段病灶组织在运 动方向上的运动幅度信息。
[0117] 可选的,可以利用主成分分析算法,确定病灶组织的运动方向以及在该运动方向 上的运动幅度信息,相应的,处理过程可以如下:根据预设的主成分分析PCA算法,对由每个 子时段中消除噪声后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭 点数目构成的矩阵进行主成分分析,得到至少一个特征值,以及至少一个特征值对应的特 征向量,其中,特征值越大表示湮灭点数目在特征值的特征向量对应的方向上变化程度越 大;选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将第一特征 向量或第二特征向量对应的方向作为病灶组织的运动方向;对于每个子时段,根据子时段 中所述消除噪声后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点 数目,以及第一特征向量或第二特征向量,确定在子时段病灶组织在运动方向上的运动幅 度信息。
[0118] 在实施中,对应每个子时段,终端得到该子时段中消除噪声后的径向位置-角度位 置-轴向位置标识组合对应的衰变修正后的湮灭点数目后,可以将其按照一定的排列顺序, 组成η XI维的向量,可以用X(ti)表示,其中,η为正整数,η可以是Sinogram数据格式中的径 向位置的数量、角度位置的数量和轴向位置标识的数量的乘积,表示第i个子时段,i = 1, 2,…M,其中,Μ是子时段的数目,每一维的数值即是该维对应的消除噪声后的径向位置-角 度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,进而,可以由每个子时段的X (t〇组成矩阵,可以用矩阵S表示,即5=[乂^ 1)4&2)4&3)-7(切)]。进而,终端可以利用 PCA(Principal Component Analysis,主成分分析)算法,对矩阵S进行主成分分析。具体 的,终端获得矩阵S后,可以计算矩阵S中每一行数值的均值,将矩阵S每行的每个数据减去 该行对应的均值,得到矩阵f,进而,可以根据公式(2)计算矩阵f的协方差矩阵Σ,
(2 ) L0120」其中,S' 1表示矩阵S'的转置矩阵,进而,可以计算协方差矩阵Σ的特征值和特征 向量,其中,特征值越大表示湮灭点数目在特征值的特征向量对应的方向上变化程度越大, 在本方案中,最大特征值对应的特征向量对应的方向可以是运动强度最大的运动对应的运 动方向(由于特征值的大小反应数据的变化程度,由于呼吸运动或心脏搏动等运动导致病 灶组织运动,进而,导致湮灭点的分布发生变化,即矩阵S中的每一位置的数据发生变化,因 此,可以将根据矩阵S得到的最大特征值对应的特征向量的方向作为运动强度最大的运动 对应的运动方向),以此类推,可以根据获取的特征向量确定某运动的运动方向,比如,一般 情况下,呼吸运动的运动强度大于心脏搏动的运动强度,可以将最大特征值对应的特征向 量对应的方向作为呼吸运动的运动方向,将次大特征值对应的特征向量对应的方向作为心 脏搏动的运动方向,特殊情况下,如果心脏搏动的运动强度大于呼吸运动的运动强度,则可 以将最大特征值对应的特征向量对应的方向作为心脏搏动的运动方向,将次大特征值对应 的特征向量对应的方向作为呼吸运动的运动方向。进而,终端可以获取最大特征值对应的 第一特征向量或次大特征值对应的第二特征向量,将第一特征向量或第二特征向量对应的 方向作为病灶组织的运动方向,即可以根据需要,选取其中的一个特征向量对应的方向作 为病灶组织的运动方向,优选的,可以将第一特征向量对应的方向作为病灶组织的运动方 向。
[0121] 终端确定病灶组织的运动方向后,可以计算上述矩阵f的每一列在第一特征向量 或第二特征向量的投影,即可以计算第一特征向量或者第二特征向量与减去均值后的X (t〇的点积,其中,可以将得到的每个子时段对应的投影,作为在子时段病灶组织在运动方 向上的运动幅度信息,可以记为W(t),即W(t)可以用于表示病灶组织在t子时段的运动幅度 信息。
[0122] 步骤103,根据在每个子时段病灶组织在运动方向上的运动幅度信息,确定门控信 号,基于门控信号进行PET图像重建。
[0123] 其中,门控信号可以是由每个子时段对应的运动幅度信息确定出的运动信号。
[0124] 在实施中,终端确定出病灶组织在每个子时段的在运动方向上的运动幅度信息 后,可以将每个子时段对应的运动幅度信息,作为门控信号。可以从门控信号中,提取被扫 描对象呼气最大时(或者吸气最大时)对应的各个时刻,将每相邻两个时刻之间的时长等分 为预设数目个子时长,将每个时长内相对应的子时长中检测到的湮灭点位置信息划分在一 起,其中,可以认为划分在一起的湮灭点位置信息是被扫描对象处于相同运动阶段时检测 到的,进而,可以根据划分在一起的湮灭点位置信息进行PET图像重建,得到相应运动阶段 的PET图像。例如,被扫描对象呼气最大时对应的时刻分别是tl、t2、t3、t4、t5,分别将t2-tl、t3-t2、t4_t3、t5_t4的时长划分为3等分,编号分别为1-3,可以将在所有编号为1的子时 长内检测到的各湮灭点位置信息划分在一起,根据这些湮灭点位置进行PET图像重建,得到 对应该运动阶段的PET图像,将在所有编号为2的子时长内检测到的各湮灭点位置信息划分 在一起,根据这些湮灭点位置进行PET图像重建,得到对应该运动阶段的PET图像,在所有编 号为3的子时长内检测到的各湮灭点位置信息划分在一起,根据这些湮灭点位置进行PET图 像重建,得到对应该运动阶段的PET图像,进而,可以得到多个运动阶段的PET图像。
[0125] 可选的,终端可以将经过插值处理的运动幅度信息,作为门控信号,相应的,步骤 103的处理过程可以如下:对每个子时段和在每个子时段病灶组织在运动方向上的运动幅 度信息构成的运动幅度信息-时间函数进行插值处理,得到插值处理后的运动幅度信息-时 间函数;将插值处理后的运动幅度信息-时间函数作为门控信号。
[0126] 在实施中,终端得到终端确定出病灶组织在每个子时段的在运动方向上的运动幅 度信息W(t)后,可以基于预设的插值算法对W(t)进行插值处理,得到插值处理后的运动幅 度信息-时间函数^(t),具体的,可以对W(t)进行线性插值,通过相邻时刻对应的运动幅度 信息,确定该相邻时刻的中间时刻对应的运动幅度信息,本方案中也可以采用非线性插值 算法,不对其进行限定。得到插值处理后的运动幅度信息-时间函数矿(t),可以将,(t)作 为门控信号。这样,可以更准确的获取最大运动幅度信息对应的时刻,以便进行PET图像重 建。
[0127] 可选的,还可以将经过滤波后的每个子时段对应的运动幅度信息作为门控信号, 相应的,步骤103的处理过程可以如下:对每个子时段和在每个子时段病灶组织在运动方向 上的运动幅度信息构成的运动幅度信息-时间函数进行频域变换处理,得到运动幅度信息-时间函数对应的频谱;基于具有预设中心频率和预设带宽的带通滤波器,对频谱进行滤波, 得到滤波后的频谱;对滤波后的频谱进行频域变换的逆变换,得到门控信号。
[0128] 在实施中,终端得到每个子时段t对应的运动幅度信息W(t)后,可以对W(t)进行快 速傅里叶变换(FFT),进而,得到W(t)的频谱,其中,W(t)的频谱可以按照公式(3)计算得到,
[0129] H(f)= |FFT(ff(t)) (3)
[0130] 其中,H( f)即是W( t)的频谱,f表示频率,进而,终端可以对运动幅度信息W( t)的频 谱进行滤波,其中,可以通过具有预设中心频率和预设带宽的带通滤波器(比如高斯滤波 器)对其进行滤波,得到滤波后的频谱,之后,对滤波后的频谱进行反傅里叶变换处理,得到 门控信号。
[0131] 另外,当运动幅度信息W(t)是呼吸运动对应的运动信号时,预设中心频率可以为 预设呼吸频率值(比如,预设呼吸频率值可以是〇. 5Hz)时,得到的门控信号为呼吸运动对应 的门控信号,当运动幅度信息W(t)是心脏搏动对应的运动信号时,预设中心频率为预设心 跳频率值(比如,预设心跳频率值可以是1Hz)时,得到的门控信号为心脏搏动对应的门控信 号。
[0132] 本发明实施例中,可以根据PET设备检测到的各湮灭点位置信息确定门控信号,无 需在PET设备检测各湮灭点位置信息的时候,增加额外的外部设备采集门控信号,从而,可 以降低诊断系统成本。
[0133] 基于相同的技术构思,本发明实施例还提供了一种确定门控信号的装置,如图4所 示,该装置包括:
[0134] 获取模块410,用于获取正电子发射断层显像PET设备在检测时段内的每个预设时 长的子时段中检测到的各湮灭点位置信息;
[0135] 第一确定模块420,用于根据所述每个子时段对应的各湮灭点位置信息,确定病灶 组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息;
[0136] 第二确定模块430,用于根据在每个子时段所述病灶组织在所述运动方向上的运 动幅度信息,确定门控信号,基于所述门控信号进行PET图像重建。
[0137] 可选的,如图5所示,所述第一确定模块420,包括:
[0138] 第一确定子模块4201,用于对于每个子时段,根据所述子时段对应的各湮灭点位 置信息,确定所述子时段中不同湮灭点位置信息对应的湮灭点数目;
[0139]第二确定子模块4202,用于根据所述每个子时段中不同湮灭点位置信息对应的湮 灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上 的运动幅度信息。
[0140]可选的,所述第一确定子模块4201,用于:
[0141] 对于每个子时段,根据所述子时段对应的各湮灭点位置信息中的径向位置、角度 位置和轴向位置标识,确定每个径向位置-角度位置-轴向位置标识组合对应的湮灭点数 目,其中,所述轴向位置标识是用于表示湮灭点的轴向位置的信息;
[0142] 所述第二确定子模块4202,用于:
[0143] 对于每个子时段,将所述子时段中所述每个径向位置-角度位置-轴向位置标识组 合中满足预设接近度条件的多个径向位置-角度位置组合合并为一个新的径向位置-角度 位置组合,并将多个径向位置-角度位置组合对应的湮灭点数目之和,作为所述新的径向位 置-角度位置组合对应的湮灭点数目,得到所述子时段中新的径向位置-角度位置-轴向位 置标识组合对应的湮灭点数目;
[0144] 根据所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合对应的 湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向 上的运动幅度信息。
[0145] 可选的,所述第二确定子模块4202,用于:
[0146] 对所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合进行数据 重组处理,得到数据重组后的径向位置-角度位置-轴向位置标识组合,根据所述每个子时 段中所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子 时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目;
[0147] 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组 合对应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述 运动方向上的运动幅度信息。
[0148] 可选的,所述第二确定子模块4202,用于:
[0149] 根据核素衰变公式和所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子时段中所述数据重组后的径向位 置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目;
[0150] 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病 灶组织在所述运动方向上的运动幅度信息。
[0151]可选的,所述第二确定子模块4202,用于:
[0152] 对于所述每个子时段中每个所述数据重组后的径向位置-角度位置-轴向位置标 识组合中的径向位置-角度位置组合,根据所述径向位置-角度位置组合所属的数据重组后 的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定所述径向 位置-角度位置组合对应的衰变修正后的湮灭点数目;
[0153] 将每个子时段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度 位置组合,所属的数据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮 灭点数目设置为零,得到所述每个子时段中消除噪声后的径向位置-角度位置-轴向位置标 识组合对应的衰变修正后的湮灭点数目;
[0154] 根据所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病 灶组织在所述运动方向上的运动幅度信息。
[0155] 可选的,所述第二确定子模块4202,用于:
[0156] 根据预设的主成分分析PCA算法,对由所述每个子时段中所述消除噪声后的径向 位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目构成的矩阵进行主成 分分析,得到至少一个特征值,以及所述至少一个特征值对应的特征向量,其中,所述特征 值越大表示湮灭点数目在所述特征值的特征向量对应的方向上变化程度越大;
[0157] 选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将所述 第一特征向量或所述第二特征向量对应的方向作为病灶组织的运动方向;
[0158] 对于每个子时段,根据所述子时段中所述消除噪声后的径向位置-角度位置-轴向 位置标识组合对应的衰变修正后的湮灭点数目,以及所述第一特征向量或所述第二特征向 量,确定在所述子时段所述病灶组织在所述运动方向上的运动幅度信息。
[0159]可选的,所述第二确定模块430,用于:
[0160] 对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度 信息构成的运动幅度信息-时间函数进行插值处理,得到插值处理后的运动幅度信息-时间 函数;
[0161] 将所述插值处理后的运动幅度信息-时间函数作为门控信号。
[0162] 可选的,如图6所示,所述第二确定模块430,包括:
[0163] 频域变换子模块4301,用于对所述每个子时段和在每个子时段所述病灶组织在所 述运动方向上的运动幅度信息构成的运动幅度信息-时间函数进行频域变换处理,得到所 述运动幅度信息-时间函数对应的频谱;
[0164] 滤波子模块4302,用于基于具有预设中心频率和预设带宽的带通滤波器,对所述 频谱进行滤波,得到滤波后的频谱;
[0165] 频域逆变换子模块4303,用于对所述滤波后的频谱进行所述频域变换的逆变换, 得到门控信号。
[0166] 本发明实施例中,可以根据PET设备检测到的各湮灭点位置信息确定门控信号,无 需在PET设备检测各湮灭点位置信息的时候,增加额外的外部设备采集门控信号,从而,可 以降低诊断系统成本。
[0167] 需要说明的是:上述实施例提供的确定门控信号的装置在确定门控信号时,仅以 上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不 同的功能模块完成,即将终端的内部结构划分成不同的功能模块,以完成以上描述的全部 或者部分功能。另外,上述实施例提供的确定门控信号的装置与确定门控信号的方法实施 例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
[0168] 本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件 来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读 存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
[0169] 以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和 原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
【主权项】
1. 一种确定门控信号的方法,其特征在于,所述方法包括: 获取正电子发射断层显像PET设备在检测时段内的每个预设时长的子时段中检测到的 各湮灭点位置信息; 根据所述每个子时段对应的各湮灭点位置信息,确定病灶组织的运动方向,确定在每 个子时段所述病灶组织在所述运动方向上的运动幅度信息; 根据在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,确定门控信号, 基于所述门控信号进行PET图像重建。2. 根据权利要求1所述的方法,其特征在于,所述根据所述每个子时段对应的各湮灭点 位置信息,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上 的运动幅度信息,包括: 对于每个子时段,根据所述子时段对应的各湮灭点位置信息,确定所述子时段中不同 湮灭点位置信息对应的湮灭点数目; 根据所述每个子时段中不同湮灭点位置信息对应的湮灭点数目,确定病灶组织的运动 方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息。3. 根据权利要求2所述的方法,其特征在于,所述对于每个子时段,根据所述子时段对 应的各湮灭点位置信息,确定所述子时段中不同湮灭点位置信息对应的湮灭点数目,包括: 对于每个子时段,根据所述子时段对应的各湮灭点位置信息中的径向位置、角度位置 和轴向位置标识,确定每个径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,其 中,所述轴向位置标识是用于表示湮灭点的轴向位置的信息; 所述根据所述每个子时段中不同湮灭点位置信息对应的湮灭点数目,确定病灶组织的 运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括: 对于每个子时段,将所述子时段中所述每个径向位置-角度位置-轴向位置标识组合中 满足预设接近度条件的多个径向位置-角度位置组合合并为一个新的径向位置-角度位置 组合,并将多个径向位置-角度位置组合对应的湮灭点数目之和,作为所述新的径向位置-角度位置组合对应的湮灭点数目,得到所述子时段中新的径向位置-角度位置-轴向位置标 识组合对应的湮灭点数目; 根据所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭 点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的 运动幅度信息。4. 根据权利要求3所述的方法,其特征在于,所述根据所述每个子时段中所述新的径向 位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定病灶组织的运动方向,确定在 每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括: 对所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合进行数据重组 处理,得到数据重组后的径向位置-角度位置-轴向位置标识组合,根据所述每个子时段中 所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子时段 中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目; 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对 应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动 方向上的运动幅度信息。5. 根据权利要求4所述的方法,其特征在于,所述根据所述每个子时段中所述数据重组 后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定病灶组织的运动方 向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括: 根据核素衰变公式和所述每个子时段中所述数据重组后的径向位置-角度位置-轴向 位置标识组合对应的湮灭点数目,确定所述每个子时段中所述数据重组后的径向位置-角 度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目; 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组 织在所述运动方向上的运动幅度信息。6. 根据权利要求5所述的方法,其特征在于,所述根据所述每个子时段中所述数据重组 后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定病灶组 织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括: 对于所述每个子时段中每个所述数据重组后的径向位置-角度位置-轴向位置标识组 合中的径向位置-角度位置组合,根据所述径向位置-角度位置组合所属的数据重组后的径 向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定所述径向位 置-角度位置组合对应的衰变修正后的湮灭点数目; 将每个子时段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度位置 组合,所属的数据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮灭点 数目设置为零,得到所述每个子时段中消除噪声后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目; 根据所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组 织在所述运动方向上的运动幅度信息。7. 根据权利要求6所述的方法,其特征在于,所述根据所述每个子时段中所述消除噪声 后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定病灶组 织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息,包括: 根据预设的主成分分析PCA算法,对由所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目构成的矩阵进行主成分分 析,得到至少一个特征值,以及所述至少一个特征值对应的特征向量,其中,所述特征值越 大表示湮灭点数目在所述特征值的特征向量对应的方向上变化程度越大; 选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将所述第一 特征向量或所述第二特征向量对应的方向作为病灶组织的运动方向; 对于每个子时段,根据所述子时段中所述消除噪声后的径向位置-角度位置-轴向位置 标识组合对应的衰变修正后的湮灭点数目,以及所述第一特征向量或所述第二特征向量, 确定在所述子时段所述病灶组织在所述运动方向上的运动幅度信息。8. 根据权利要求1所述的方法,其特征在于,所述根据在每个子时段所述病灶组织在所 述运动方向上的运动幅度信息,确定门控信号,包括: 对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度信息 构成的运动幅度信息-时间函数进行插值处理,得到插值处理后的运动幅度信息-时间函 数; 将所述插值处理后的运动幅度信息-时间函数作为门控信号。9. 根据权利要求1所述的方法,其特征在于,所述根据在每个子时段所述病灶组织在所 述运动方向上的运动幅度信息,确定门控信号,包括: 对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度信息 构成的运动幅度信息-时间函数进行频域变换处理,得到所述运动幅度信息-时间函数对应 的频谱; 基于具有预设中心频率和预设带宽的带通滤波器,对所述频谱进行滤波,得到滤波后 的频谱; 对所述滤波后的频谱进行所述频域变换的逆变换,得到门控信号。10. -种确定门控信号的装置,其特征在于,所述装置包括: 获取模块,用于获取正电子发射断层显像PET设备在检测时段内的每个预设时长的子 时段中检测到的各湮灭点位置信息; 第一确定模块,用于根据所述每个子时段对应的各湮灭点位置信息,确定病灶组织的 运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动幅度信息; 第二确定模块,用于根据在每个子时段所述病灶组织在所述运动方向上的运动幅度信 息,确定门控信号,基于所述门控信号进行PET图像重建。11. 根据权利要求10所述的装置,其特征在于,所述第一确定模块,包括: 第一确定子模块,用于对于每个子时段,根据所述子时段对应的各湮灭点位置信息,确 定所述子时段中不同湮灭点位置信息对应的湮灭点数目; 第二确定子模块,用于根据所述每个子时段中不同湮灭点位置信息对应的湮灭点数 目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的运动 幅度信息。12. 根据权利要求11所述的装置,其特征在于,所述第一确定子模块,用于: 对于每个子时段,根据所述子时段对应的各湮灭点位置信息中的径向位置、角度位置 和轴向位置标识,确定每个径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,其 中,所述轴向位置标识是用于表示湮灭点的轴向位置的信息; 所述第二确定子模块,用于: 对于每个子时段,将所述子时段中所述每个径向位置-角度位置-轴向位置标识组合中 满足预设接近度条件的多个径向位置-角度位置组合合并为一个新的径向位置-角度位置 组合,并将多个径向位置-角度位置组合对应的湮灭点数目之和,作为所述新的径向位置-角度位置组合对应的湮灭点数目,得到所述子时段中新的径向位置-角度位置-轴向位置标 识组合对应的湮灭点数目; 根据所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭 点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动方向上的 运动幅度信息。13. 根据权利要求12所述的装置,其特征在于,所述第二确定子模块,用于: 对所述每个子时段中所述新的径向位置-角度位置-轴向位置标识组合进行数据重组 处理,得到数据重组后的径向位置-角度位置-轴向位置标识组合,根据所述每个子时段中 所述新的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目,确定所述每个子时段 中所述数据重组后的径向位置-角度位置-轴向位置标识组合对应的湮灭点数目; 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对 应的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组织在所述运动 方向上的运动幅度信息。14. 根据权利要求13所述的装置,其特征在于,所述第二确定子模块,用于: 根据核素衰变公式和所述每个子时段中所述数据重组后的径向位置-角度位置-轴向 位置标识组合对应的湮灭点数目,确定所述每个子时段中所述数据重组后的径向位置-角 度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目; 根据所述每个子时段中所述数据重组后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组 织在所述运动方向上的运动幅度信息。15. 根据权利要求14所述的装置,其特征在于,所述第二确定子模块,用于: 对于所述每个子时段中每个所述数据重组后的径向位置-角度位置-轴向位置标识组 合中的径向位置-角度位置组合,根据所述径向位置-角度位置组合所属的数据重组后的径 向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目,确定所述径向位 置-角度位置组合对应的衰变修正后的湮灭点数目; 将每个子时段中衰变修正后的湮灭点数目未超过预设数目阈值的径向位置-角度位置 组合,所属的数据重组后的径向位置-角度位置-轴向位置标识组合的衰变修正后的湮灭点 数目设置为零,得到所述每个子时段中消除噪声后的径向位置-角度位置-轴向位置标识组 合对应的衰变修正后的湮灭点数目; 根据所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组合对 应的衰变修正后的湮灭点数目,确定病灶组织的运动方向,确定在每个子时段所述病灶组 织在所述运动方向上的运动幅度信息。16. 根据权利要求15所述的装置,其特征在于,所述第二确定子模块,用于: 根据预设的主成分分析PCA算法,对由所述每个子时段中所述消除噪声后的径向位置-角度位置-轴向位置标识组合对应的衰变修正后的湮灭点数目构成的矩阵进行主成分分 析,得到至少一个特征值,以及所述至少一个特征值对应的特征向量,其中,所述特征值越 大表示湮灭点数目在所述特征值的特征向量对应的方向上变化程度越大; 选取最大特征值对应的第一特征向量或次大特征值对应的第二特征向量,将所述第一 特征向量或所述第二特征向量对应的方向作为病灶组织的运动方向; 对于每个子时段,根据所述子时段中所述消除噪声后的径向位置-角度位置-轴向位置 标识组合对应的衰变修正后的湮灭点数目,以及所述第一特征向量或所述第二特征向量, 确定在所述子时段所述病灶组织在所述运动方向上的运动幅度信息。17. 根据权利要求10所述的装置,其特征在于,所述第二确定模块,用于: 对所述每个子时段和在每个子时段所述病灶组织在所述运动方向上的运动幅度信息 构成的运动幅度信息-时间函数进行插值处理,得到插值处理后的运动幅度信息-时间函 数; 将所述插值处理后的运动幅度信息-时间函数作为门控信号。18.根据权利要求10所述的装置,其特征在于,所述第二确定模块,包括: 频域变换子模块,用于对所述每个子时段和在每个子时段所述病灶组织在所述运动方 向上的运动幅度信息构成的运动幅度信息-时间函数进行频域变换处理,得到所述运动幅 度信息-时间函数对应的频谱; 滤波子模块,用于基于具有预设中心频率和预设带宽的带通滤波器,对所述频谱进行 滤波,得到滤波后的频谱; 频域逆变换子模块,用于对所述滤波后的频谱进行所述频域变换的逆变换,得到门控 信号。
【文档编号】G06T5/00GK106097384SQ201610374708
【公开日】2016年11月9日
【申请日】2016年5月30日
【发明人】朱传伟, 李楠
【申请人】江苏赛诺格兰医疗科技有限公司
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1