《一种判断待检对象发生运动的方法和装置.pdf》由会员分享,可在线阅读,更多相关《一种判断待检对象发生运动的方法和装置.pdf(13页完整版)》请在专利查询网上搜索。
1、(10)授权公告号 CN 102805634 B (45)授权公告日 2015.01.21 CN 102805634 B (21)申请号 201110145206.X (22)申请日 2011.05.31 A61B 6/03(2006.01) (73)专利权人 上海西门子医疗器械有限公司 地址 201318 上海市南汇区周祝公路 278 号 (72)发明人 杜春华 王薇 US 2004/0136490 A1,2004.06.15, 全文 . CN 1895174 A,2007.01.17, US 2005/0175140 A1,2005.08.11, CN 102755170 A,2012.1。
2、0.31, CN 1787077 A,2006.06.14, 全文 . JP P2006-376 A,2006.01.05, 全文 . CN 1918605 A,2007.02.21, 全文 . US 2005/0238135 A1,2005.10.27, 全文 . CN 101222703 A,2008.07.16, 全文 . Oliver Morin, et al.MEGVOLTAGE CONE-BEAM CT:SYSTEM DESCRIPTION AND CLINICAL APPLICATIONS.Medical Dosimetry .2006, 第 31 卷 ( 第 1 期 ), 第 。
3、51-61 页 . Eberhard Hansis, et al.Motion-Compensated Iterative Sparse Data Reconstruction for Interventional 3-D Coronary Artery Imaging. 2008 IEEE Nuclear Science Symposium Conference Record .2008, 第 4280-4284 页 . (54) 发明名称 一种判断待检对象发生运动的方法和装置 (57) 摘要 本发明涉及医疗设备领域, 特别是一种判断 待检对象发生运动的方法, 所述方法包括 : 获取 所述原。
4、始投影数据组的补偿数据组 ; 计算所述原 始投影数据组与所述补偿数据组的距离矩阵 ; 取 所述距离矩阵的某一列, 找出该列中最小值所对 应的行值, 若该行值与该列值的差值的绝对值大 于一预设阈值, 则所述原始投影数据组中存在运 动伪影。本发明还提供一种判断待检对象发生运 动的装置。采用本发明能直接从 0, 180 +2 的原始投影数据组中就判断出待检对象在本次扫 描中是否发生运动, 并进一步得到发生运动的具 体位置, 同时减少了待检对象接收的 X 射线辐射 剂量。 (51)Int.Cl. (56)对比文件 审查员 戚永娟 权利要求书 2 页 说明书 7 页 附图 3 页 (19)中华人民共和国。
5、国家知识产权局 (12)发明专利 权利要求书2页 说明书7页 附图3页 (10)授权公告号 CN 102805634 B CN 102805634 B 1/2 页 2 1. 一种判断待检对象发生运动的方法, 其特征在于, 所述方法包括 : 获取原始投影数据组的补偿数据组 ; 计算所述原始投影数据组与所述补偿数据组的距离矩阵 ; 取所述距离矩阵的某一列, 找出该列中最小值所对应的行值, 若该行值与该列的列值 的差值的绝对值大于一预设阈值, 则判断待检对象发生运动。 2. 根据权利要求 1 所述的方法, 其特征在于, 所述距离矩阵为欧几里得距离矩阵。 3. 根据权利要求 1 所述的方法, 其特征在。
6、于, 所述获取补偿数据组包括 : 读取所述原始投影数据组中与一投影数据对应的第一投影角度和第一探测器通道 ; 根据所述第一投影角度和所述第一探测器通道来计算一第二投影角度和第二探测器 通道, 其中所述第二投影角度下的第二探测器通道与所述第一投影角度下的第一探测器通 道位于同一 X 射线传播路径上 ; 将所述第二投影角度下第二探测器通道所采集到的投影数据作为所述补偿数据组中 的数据。 4. 根据权利要求 3 所述的方法, 其特征在于, 所述计算第二投影角度和第二探测器通 道包括 : 根据所述第一探测器通道按如下公式来计算所述X射线传播路径偏离X射线束中心的 角度 : -n1*, 其中, n1 为。
7、所述第一探测器通道, 为所述探测器通道的角度步距, 为所述 X 射 线传播路径与 X 射线束中心的夹角, 为所述 X 射线束的一半张角 ; 根据所述夹角按如下公式来计算所述第二探测器通道 : n2 (+)/, 其中, n2 为所述第二探测器通道, 且所述第二探测器通道与所述第一探测器通道关于 所述探测器的中心呈镜像对称 ; 根据所述夹角按如下公式来计算所述第二投影角度 : 2 1pi-2, 其中, 2 为所述第二投影角度, 1 为所述第一投影角度。 5. 根据权利要求 4 所述的方法, 其特征在于, 所述第一投影角度的范围为 0, 180 +2。 6. 根据权利要求 1 所述的方法, 其特征在。
8、于, 所述预设阈值为 20。 7. 根据权利要求 1 所述的方法, 其特征在于, 所述方法进一步包括 : 若该行值与该列值的差值的绝对值大于所述预设阈值, 则根据该行值得到所述待检对 象发生运动的具体位置。 8. 一种判断待检对象发生运动的装置, 其特征在于, 所述装置 (10) 包括 : 一补偿数据获取组件 (11), 用于获取原始投影数据组的补偿数据组 ; 一距离矩阵计算组件 (12), 用于计算所述原始投影数据组与所述补偿数据组的距离矩 阵 ; 一运动判断组件 (13), 取所述距离矩阵的某一列, 找出该列中最小值所对应的行值, 若 该行值与该列的列值的差值的绝对值大于一预设阈值, 则判。
9、断待检对象发生运动。 权 利 要 求 书 CN 102805634 B 2 2/2 页 3 9. 根据权利要求 8 所述的装置, 其特征在于, 所述补偿数据获取组件 (11) 包括 : 一读取模块 (111), 用于读取所述原始投影数据组中与一投影数据对应的第一投影角 度和第一探测器通道 ; 一计算模块 (112), 用于根据所述第一投影角度和所述第一探测器通道来计算一第二 投影角度和一第二探测器通道, 其中所述第二投影角度下的第二探测器通道与所述第一投 影角度下的第一探测器通道位于同一 X 射线传播路径上 ; 一补偿数据获取模块 (113), 其将所述第二投影角度下第二探测器通道所采集到的投。
10、 影数据作为所述补偿数据组中的数据。 10. 根据权利要求 9 所述的装置, 其特征在于, 所述计算模块 (112) 包括 : 一角度计算单元, 用于根据所述第一探测器通道按如下公式来计算所述 X 射线传播路 径偏离 X 射线束中心的角度 : -n1*, 其中, n1 为所述第一探测器通道, 为所述探测器通道的角度步距, 为所述 X 射 线传播路径与 X 射线束中心的夹角, 为所述 X 射线束的一半张角 ; 一第二通道计算单元, 用于根据所述夹角按如下公式来计算所述第二探测器通道 : n2 (+)/, 其中, n2 为所述第二探测器通道, 且所述第二探测器通道与所述第一探测器通道关于 所述探测。
11、器的中心呈镜像对称 ; 一第二角度计算单元, 用于根据所述夹角按如下公式来计算所述第二投影角度 : 2 1pi-2, 其中, 2 为所述第二投影角度, 1 为所述第一投影角度。 11. 根据权利要求 8 所述的装置, 其特征在于, 所述装置 (10) 进一步包括 : 一运动位置检测组件 (14), 用于在该行值与该列值的差值的绝对值大于所述预设阈值 时, 根据该行值得到所述待检对象发生运动的具体位置。 权 利 要 求 书 CN 102805634 B 3 1/7 页 4 一种判断待检对象发生运动的方法和装置 技术领域 0001 本发明涉及医疗设备领域, 特别是一种判断待检对象发生运动的方法和装。
12、置。 背景技术 0002 随着医疗水平的进步, X 射线计算机断层成像 (Computed Tomography, CT) 设备得 到了广泛的应用。当对待检对象进行 CT 扫描时, 旋转机架绕 Z 轴旋转 360 度, 并同时带动 包括 X 射线源和探测器的数据采集系统同步围绕待检对象旋转 360, 以从 0, 360 的 投影角度来获取 X 射线投影数据, 其中 Z 轴为 CT 检查床进出旋转机架的水平方向。在每个 投影角度处, 探测器的 M 个通道同时接收 X 射线投影数据, 并将其传送至 CT 的图像重建系 统用于生成重建图像。然而由于人体不可避免会发生运动, 所以在对人体进行 CT 扫。
13、描时, 重建图像会因人体的运动而存在运动伪影, 从而影响重建图像的质量。 0003 公开号为 102018524A 的发明专利申请公开了一种伪影检测方法和装置, 其通过 计算 0, 360 范围内每一投影角度所对应的组中每个投影数据分别与前一投影角度或 后一投影角度所对应的组中对应相同通道的投影数据的差值来判断重构图像中有无伪影, 若所得差值中的任一差值大于等于预先设置的阈值上限值或小于等于预先设置的阈值下 限值, 则判断重构图像中存在伪影。然而该发明专利申请适用于检查 CT 扫描仪是否有故 障, 以避免人工检查的误差。 发明内容 0004 本发明的目的在于提供一种判断待检对象是否发生运动的方。
14、法和装置, 其只需数 据采集系统围绕待检对象进行 0, 180 +2 投影即可判断出待检对象是否发生运动, 以 减少待检对象接收的辐射剂量。 0005 有鉴于此, 本发明提出一种判断待检对象发生运动的方法, 所述方法包括 : 获取所 述原始投影数据组的补偿数据组 ; 计算所述原始投影数据组与所述补偿数据组的距离矩 阵 ; 取所述距离矩阵的某一列, 找出该列中最小值所对应的行值, 若该行值与该列值的差值 的绝对值大于一预设阈值, 则判断待检对象发生运动。可见本发明能直接从原始投影数据 组中就判断出待检对象在本次扫描中是否发生运动, 准确率和可靠性更高。 0006 根据本发明的一个实施例, 所述距。
15、离矩阵为欧几里得距离矩阵。 0007 根据本发明的再一个实施例, 所述获取补偿数据组包括 : 读取所述原始投影数据 组中与一投影数据对应的第一投影角度和第一探测器通道 ; 根据所述第一投影角度和所述 第一探测器通道来计算一第二投影角度和第二探测器通道, 其中所述第二投影角度下的第 二探测器通道与所述第一投影角度下的第一探测器通道位于同一 X 射线传播路径上 ; 将所 述第二投影角度下第二探测器通道所采集到的投影数据作为所述补偿数据组中的数据。 0008 根据本发明的又一个实施例, 所述计算第二投影角度和第二探测器通道包括 : 根 据所述第一探测器通道按如下公式来计算所述 X 射线传播路径偏离 。
16、X 射线束中心的角度 : -n1*, 其中, n1为所述第一探测器通道, 为所述探测器通道的角度步距, 说 明 书 CN 102805634 B 4 2/7 页 5 为所述 X 射线传播路径与 X 射线束中心的夹角, 为所述 X 射线束的一半张角 ; 根据所述 夹角按如下公式来计算所述第二探测器通道 : n2 (+)/, 其中, n2 为所述第二探 测器通道, 且所述第二探测器通道与所述第一探测器通道关于所述探测器的中心呈镜像对 称 ; 根据所述夹角按如下公式来计算所述第二投影角度 : 2 1pi-2, 其中, 2 为 所述第二投影角度, 1 为所述第一投影角度。 0009 根据本发明的又一个。
17、实施例, 所述第一投影角度的范围为 0, 180 +2。这样 数据采集系统只需围绕待检对象进行 0, 180 +2 的投影即可得到所述原始投影数据, 相比 360的投影方式减少了待检对象接收的辐射剂量, 同时减少了原始投影数据组的数 据量, 减少了数据存储空间。 0010 根据本发明的又一个实施例, 所述预设阈值为 20。 0011 进一步, 所述方法包括 : 若该行值与该列值的差值的绝对值大于所述预设阈值, 则 根据该行值得到所述待检对象发生运动的具体位置。 0012 本发明还提供一种计算机程序, 使计算机执行所述判断待检对象发生运动的方 法。 0013 本发明还提供一种计算机可读的记录介质。
18、, 存储了所述的计算机程序。 0014 本发明还提供一种判断待检对象发生运动的装置, 所述装置包括 : 一补偿数据获 取组件, 用于获取所述原始投影数据组的补偿数据组 ; 一距离矩阵计算组件, 用于计算所述 原始投影数据组与所述补偿数据组的距离矩阵 ; 一运动判断组件, 取所述距离矩阵的某一 列, 找出该列中最小值所对应的行值, 若该行值与该列值的差值的绝对值大于一预设阈值, 则判断待检对象发生运动。 可见本发明能直接从原始投影数据组中就判断出待检对象在本 次扫描中是否发生运动, 所以准确率和可靠性更高。 0015 根据本发明的再一个实施例, 所述补偿数据获取组件包括 : 一读取模块, 用于读。
19、取 所述原始投影数据组中与一投影数据对应的第一投影角度和第一探测器通道 ; 一计算模 块, 用于根据所述第一投影角度和所述第一探测器通道来计算一第二投影角度和一第二探 测器通道, 其中所述第二投影角度下的第二探测器通道与所述第一投影角度下的第一探测 器通道位于同一 X 射线传播路径上 ; 一补偿数据获取模块, 其将所述第二投影角度下第二 探测器通道所采集到的投影数据作为所述补偿数据组中的数据。 0016 根据本发明的又一个实施例, 所述计算模块包括 : 一角度计算单元, 用于根据所述 第一探测器通道按如下公式来计算所述 X 射线传播路径偏离 X 射线束中心的角度 : -n1*, 其中, n1 。
20、为所述第一探测器通道, 为所述探测器通道的角度步距, 为所 述 X 射线传播路径与 X 射线束中心的夹角, 为所述 X 射线束的一半张角 ; 一第二通道计 算单元, 用于根据所述夹角按如下公式来计算所述第二探测器通道 : n2 (+)/, 其中, n2 为所述第二探测器通道, 且所述第二探测器通道与所述第一探测器通道关于所述 探测器的中心呈镜像对称 ; 一第二角度计算单元, 用于根据所述夹角按如下公式来计算所 述第二投影角度 : 2 1pi-2, 其中, 2 为所述第二投影角度, 1 为所述第一投 影角度。 0017 进一步, 所述装置包括 : 一运动位置检测组件, 用于在该行值与该列值的差值。
21、的绝 对值大于所述预设阈值时, 根据该行值得到所述待检对象发生运动的具体位置。 0018 从上述方案中可以看出, 由于本发明能直接从原始投影数据组中就判断出待检对 说 明 书 CN 102805634 B 5 3/7 页 6 象在本次扫描中是否发生运动, 所以准确率和可靠性更高 ; 而且所述原始投影数据为数据 采集系统围绕待检对象进行 0, 180 +2 的投影, 相比 360的投影方式减少了待检对 象接收的辐射剂量, 同时减少了原始投影数据组的数据量, 减少了数据存储空间。另外, 本 发明还能得到待检对象发生运动的具体位置。 附图说明 0019 下面将通过参照附图详细描述本发明的实施例, 使。
22、本领域技术人员更清楚本发明 的上述及其它特征和优点, 附图中 : 0020 图 1 为本发明判断待检对象发生运动的方法示意图。 0021 图 2 为本发明获取补偿数据的原理示意图。 0022 图 3a 为待检对象未发生运动时本发明距离矩阵中某列数据的示意图。 0023 图 3b 为待检对象发生运动时本发明距离矩阵中某列数据的示意图。 0024 图 4 为本发明判断待检对象发生运动的装置示意图。 0025 图 5 为本发明补偿数据获取组件的示意图。 0026 附图标记 0027 步骤 101 108 判断待检对象发生运动的步骤 0028 1X 射线源 0029 3 探测器 0030 10 本发明。
23、判断待检对象发生运动的装置 0031 11 补偿数据获取组件 0032 12 距离矩阵计算组件 0033 13 运动判断组件 0034 14 运动位置检测组件 0035 111 读取模块 0036 112 计算模块 0037 113 补偿数据获取模块 具体实施方式 0038 为使本发明的目的、 技术方案和优点更加清楚, 以下举实施例对本发明进一步详 细说明。 0039 以下以 CT 设备为例, 具体说明本发明判断待检对象发生运动的方法和装置。对 目前的双排 CT 而言, 通常数据采集系统围绕待检对象同步旋转一周的时间少于 1s, 如为 0.6s, 而人体发生运动频率通常为一次 1 秒, 所以在。
24、扫描过程中, 并非每次扫描得到的原始 投影数据组中都含有待检对象的运动数据。除非待检对象患有帕金森综合症, 这种症状的 病人的运动频率非常快, 需要扫描时间很短的高端 CT 来对其进行扫描。若原始投影数据组 中都含有待检对象的运动数据, 则重建图像会因这些运动数据而存在运动伪影, 从而影响 重建图像的质量, 所以需要检测原始数据组中是否含有待检对象的运动数据。 0040 有鉴于此, 本发明提出一种判断待检对象是否发生运动的方法, 如图 1 所示。所述 方法包括 : 说 明 书 CN 102805634 B 6 4/7 页 7 0041 步骤 101, 读取原始投影数据组 R 中与一投影数据 v。
25、1 对应的第一投影角度 1 和 第一探测器通道 n1。 0042 本发明是基于如图 2 所示的数据补偿的原理来进行的。图 2 所示数据采集系统包 括一种 X 射线源 1(F1 或 F2) 和探测器 3, 探测器 3 的通道数为 M, 其起始探测器通道为 C-1, 结束探测器通道为 C-M。其中, X 射线源 1 发射一 X 射线束, X 射线束张角为 2, 探测器 3 与 X 射线源 1 一起围绕一待检对象顺时针同步旋转, 即 X 射线源从 F1 处顺时针旋转到 F2 处, 以从复数个视点来采集 X 射线投影数据。理论上, 在同一 X 射线传播路经 F1F2 上, 由于 X射线从F1到F2方向。
26、经过待检对象的衰减与从F2到F1方向经过待检对象的衰减相同, 所 以在投影角度为 1( 从 X 轴正半轴旋转到 F1 点的旋转角度 ) 时, 探测器 3 的一个通道 n1 上接收到的 X 射线投影数据 v1 应该等于投影角度为 2( 从 X 轴正半轴旋转到 F2 点的旋 转角度 ) 时, 探测器 3 的另一个通道 n2 接收的 X 射线投影数据 v2, 即 v1 v2。本发明据 此来获取补充数据组 R_C 的数据。 0043 步骤102, 根据第一探测器通道n1按如下公式(1)来计算X射线传播路径F1F2偏 离 X 射线束中心的角度 : 0044 -n1* (1) 0045 其中, 为探测器 。
27、3 中相邻两个通道间的角度步距, 为所述 X 射线束的一半 张角。 0046 步骤 103, 根据夹角 按如下公式 (2) 来计算第二探测器通道 n2 : 0047 n2 (+)/ (2) 0048 这里, 可以将公式 (1) 变形为如下公式 (3) : 0049 n1 (-)/ (3) 0050 比较公式 (2) 和 (3) 可以得出 : 第一探测器通道 n1 与第二探测器通道 n2 关于探 测器 3 的中心是镜像对称的。 0051 步骤 104, 根据夹角 按如下公式 (4) 来计算第二投影角度 2 : 0052 2 1pi- (4) 0053 由于 2, 所以 0054 2 1pi-2 。
28、(5) 0055 图 2 的实施例中 2 1-pi-2。 0056 步骤105, 将在第二投影角度2下第二探测器通道n2所采集到的投影数据v2作 为原始投影数据组 R 的补偿数据组 R_C 中的数据。 0057 采用上述步骤 101 105, 无需包括数据采集系统围绕待检对象旋转一周 (360 ), 而只需进行 0, 180 +2 的投影即可获得补偿数据组 R_C 中的全部数据。于 是, 本发明中第一投影角度的范围为0, 180+2, 根据公式(5)可知第二投影角度的范 围也为 0, 180 +2。这样就减少了待检对象接收的辐射剂量, 同时减少了原始投影数 据的数据量, 减少了数据存储空间。 。
29、0058 步骤 106, 计算原始投影数据组 R 与补偿数据组 R_C 的距离矩阵 S。 0059 根据本发明的一个实施例, 所述距离矩阵为欧几里得距离矩阵。 0060 在 CT 扫描中, 数据采集系统是根据一定的角度间隔来采集 X 射线投影数据的, 即 读片。如在一次扫描过程中本发明将 0, 180 +2 分为 N 个角度间隔, 如 N 2000, 则 说 明 书 CN 102805634 B 7 5/7 页 8 数据采集系统围绕待检对象每旋转 (180 +2)/2000), 探测器 3 就读取一次 X 射线投 影数据, 总共需要读片 2000 次。 0061 假设距离矩阵S为一N行N列的矩。
30、阵, 其中第i行第j列的数据可以是原始投影数 据组R中第i次读片的数据与补偿数据组R_C中第j次读片的数据的欧几里得(Euclidean) 距离。也就是说, 距离矩阵 S 中第 i 行的数据是原始投影数据组 R 中第 i 次读片的数据分 别与补偿数据组 R_C 中第 1 N 次读片的数据的欧几里得距离, 距离矩阵 S 中第 j 列的数 据是补偿数据组 R_C 中第 j 次读片的数据分别与原始投影数据组 R 中第 1 N 次读片的数 据的欧几里得距离。 0062 当然距离矩阵 S 第 i 行第 j 列的数据也可以是补偿数据组 R_C 中第 i 次读片的数 据与原始投影数据组 R 中第 j 次读片。
31、的数据的欧几里得 (Euclidean) 距离。 0063 步骤 107, 取距离矩阵 S 的某一列 j, 找出该列中最小值所对应的行值 i, 若该行值 i 与该列值 j 的差值的绝对值大于一预设阈值, 则所述待检对象发生运动。 0064 假设 J 1000, 则距离矩阵 S 中第 1000 列数据为补偿数据组 R_C 中第 1000 次读 片的数据分别与原始投影数据组 R 中第 1 N 次读片的数据的欧几里得距离。理论上, 距 离矩阵 S 中第 j 列的最小数据为 Sj, j, 即 i j。然而实际上, 由于各种原因如机械误差 ( 例 如探测器通道的大小、 X 射线源焦点到扫描中心的距离、 。
32、X 射线源焦点到探测器的距离、 X 射 线束张角等的误差 )、 条状伪影检测的准确率、 飞焦点、 旋转机架的速度等, 使得 S 中第 j 列 的最小数据并不一定为 Sj, j, 而是在 Sj, j附近的一个预设阈值范围内。根据本发明的一个实 施例, 所述预设阈值为 20。如图 3a 所示为待检对象未发生运动时 S 中第 1000 列数据的示 意图, 其中横坐标为 S 的行值 i, 纵坐标为 Si, j的相对大小, 图中在 i 1000 附近有一个圆 圈, 该圆圈示意了一预设阈值范围, 即当S中第1000列的最小值在该圆圈内, 则认为在该次 扫描过程中待检对象没有发生运动。反之, 若第 1000。
33、 列的最小值不在该圆圈内, 如图 3b 所 示, 则认为该次扫描过程中待检对象发生了运动。可见本发明能直接从原始投影数据组 R 中就判断出待检对象在本次扫描中是否发生运动, 准确率和可靠性更高。 0065 进一步, 本发明的方法还包括 : 0066 步骤108, , 若该行值i与该列值j的差值的绝对值大于所述预设阈值, 则根据该行 值 i 得到所述待检对象发生运动的具体位置。 0067 如图 3b 所示 S 中第 1000 列最小值为 i 1250 左右, 其相比 i 1000 附近的圆 圈有很大的偏移, 可知在本次扫描过程中约第 1250 次读片时待检对象发生了运动, 且纵坐 标有一个明显的。
34、跳跃 (sharp break)。 0068 本发明还提供一种计算机程序, 使计算机执行所述判断待检对象发生运动的方 法。 0069 本发明还提供一种计算机可读的记录介质, 存储用于使一机器执行如本文所述的 判断待检对象发生运动的方法指令。 具体地, 可以提供配有存储介质的系统或者装置, 在该 存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码, 且使该系统或者 装置的计算机 ( 或 CPU 或 MPU) 读出并执行存储在存储介质中的程序代码。这种情况下, 从 存储介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能, 因此程序代码 和存储程序代码的存储介质构成了本发明的一。
35、部分。 0070 用于提供程序代码的存储介质实施例包括软盘、 硬盘、 磁光盘、 光盘 ( 如 CD-ROM、 说 明 书 CN 102805634 B 8 6/7 页 9 CD-R、 CD-RW、 DVD-ROM、 DVD-RAM、 DVD-RW、 DVD+RW)、 磁带、 非易失性存储卡和 ROM。可选择地, 可以由通信网络从服务器计算机上下载程序代码。 0071 此外, 应该清楚的是, 不仅可以通过执行计算机所读出的程序代码, 而且可以通过 基于程序代码的指令使计算机上操作的操作系统等来完成部分或者全部的实际操作, 从而 实现上述实施例中任意一项实施例的功能。 0072 此外, 可以理解的。
36、是, 将由存储介质读出的程序代码写到插入计算机内的扩展板 中所设置的存储器中或者写到与计算机相连接的扩展单元中设置的存储器中, 随后基于程 序代码的指令使安装在扩展板或者扩展单元上的 CPU 等来执行部分和全部实际操作, 从而 实现上述实施例中任一实施例的功能。 0073 本发明还提供一种判断待检对象发生运动的装置, 如图 4 所示, 所述装置 10 包 括 : 0074 一补偿数据获取组件 11, 用于获取原始投影数据组 R 的补偿数据组 R_C。 0075 具体地, 如图 5 所示补偿数据获取组 11 件包括 : 0076 一读取模块 111, 用于读取原始投影数据组 R 中与一投影数据 。
37、v1 对应的第一投影 角度 1 和第一探测器通道 n1。 0077 一计算模块 112, 用于根据第一投影角度 1 和第一探测器通道 n1 来计算一第二 投影角度 2 和第二探测器通道 n2, 其中在第二投影角度 2 下的第二探测器通道 n2 与 在第一投影角度 1 下的第一探测器通道 n1 位于同一 X 射线传播路径 F1F2 上。具体地包 括 : 0078 一角度计算单元, 用于根据第一探测器通道 n1 按公式 (1) 来计算 X 射线传播路径 F1F2 偏离 X 射线束中心的角度 ; 0079 一第二通道计算单元, 用于根据夹角 按公式 (2) 来计算第二探测器通道 n2 ; 0080 。
38、一第二角度计算单元, 用于根据夹角 按公式 (5) 来计算第二投影角度 2。 0081 一补偿数据获取模块113, 其将计算模块112得到的第二投影角度2下第二探测 器通道 n2 所采集到的投影数据 v2 作为补偿数据组 R_C 中的数据。 0082 一距离矩阵计算组件 12, 用于计算原始投影数据组 R 与补偿数据组 R_C 的距离矩 阵 S。 0083 一运动判断组件 13, 取所述距离矩阵中的第 j 列, 找出该列中最小值所对应的行 值 i, 若该行值 i 与该列值 j 的差值的绝对值大于一预设阈值, 则认为待检对象在该次扫描 过程中发生运动。根据本发明的一个实施例, 预设阈值为 20。。
39、具体判断方式已在本发明的 方法中做了详细说明, 这里不再赘述。 0084 进一步, 所述装置 10 还包括 : 0085 一运动位置检测组件 14, 用于在该行值 i 与该列值 j 的差值的绝对值大于所述预 设阈值时, 根据该行值得到所述待检对象发生运动的具体位置。 0086 从上述方案中可以看出, 由于本发明能直接从原始投影数据组中就判断出待检对 象在本次扫描中是否发生运动, 所以准确率和可靠性更高 ; 而且所述原始投影数据组为数 据采集系统围绕待检对象进行 0, 180 +2 的投影, 相比 360的投影方式减少了待检 对象接收的辐射剂量, 同时减少了原始投影数据的数据量, 减少了数据存储。
40、空间。 0087 本发明涉及医疗设备领域, 特别是一种判断待检对象发生运动的方法, 所述方法 说 明 书 CN 102805634 B 9 7/7 页 10 包括 : 获取所述原始投影数据组的补偿数据组 ; 计算所述原始投影数据组与所述补偿数据 组的距离矩阵 ; 取所述距离矩阵的某一列, 找出该列中最小值所对应的行值, 若该行值与该 列值的差值的绝对值大于一预设阈值, 则所述原始投影数据组中存在运动伪影。本发明还 提供一种判断待检对象发生运动的装置。 采用本发明能直接从0, 180+2的原始投影 数据组中就判断出待检对象在本次扫描中是否发生运动, 并进一步得到发生运动的具体位 置, 同时减少了待检对象接收的 X 射线辐射剂量。 0088 以上所述仅为本发明的较佳实施例而已, 并不用以限制本发明, 凡在本发明的精 神和原则之内, 所作的任何修改、 等同替换、 改进等, 均应包含在本发明的保护范围之内。 说 明 书 CN 102805634 B 10 1/3 页 11 图 1 说 明 书 附 图 CN 102805634 B 11 2/3 页 12 图 2 说 明 书 附 图 CN 102805634 B 12 3/3 页 13 图 3 图 4 图 5 说 明 书 附 图 CN 102805634 B 13 。