一种CT影像中头部倾斜的自动检测及校正方法

摘要

本发明属于放疗技术领域，涉及一种CT影像中头部倾斜的自动检测及校正方法、设备和存储介质。该方法包括如下步骤：对全身CT影像进行二值化操作；进行最大连通区域查找，删除所有非最大连通区域，并去除每层CT影像中的非主要身体部位，只保留主要身体部位；判断肩部与颈部的交界层；确定最佳颌骨层；对最佳颌骨层中的颌骨区域进行椭圆拟合；记录颌骨的质心和拟合椭圆的质心；确定体位，实施头部校正。本发明可以快速有效地检测CT影像中头部的倾角及朝向，进而为对头部影像的倾斜校正提供准确的参数信息；具有计算复杂度低、耗时短、抗噪性高、数据适应好的优点；本方法还可拓展应用于对不同场景下的人体倾斜校正，具有较好的可拓展性。

权利要求

1.一种全身CT影像中头部倾斜的自动检测及校正方法，适于在计算设备中执行，其特征在于：该方法包括如下步骤：(1)对全身CT影像进行二值化操作；

(2)把步骤(1)中得到的所有二值化影像逐层进行最大连通区域查找，删除所有非最大连通区域，并去除每层CT影像中的非主要身体部位，只保留主要身体部位；

(3)判断肩部与颈部的交界层；

(4)确定最佳颌骨层；其中，确定最佳颌骨层包括：以肩部与颈部的交界层为起始层，向头部方向取得所有包含下颌骨的CT图；对所有包含下颌骨的CT图以颌骨HU阈值做二值化得到颌骨区域，对所有颌骨区域按照最大连通面积进行排序，面积最大者则为最佳颌骨层；

(5)对最佳颌骨层中的颌骨区域进行椭圆拟合；

(6)根据对最佳颌骨层中的颌骨区域的椭圆拟合，判定椭圆长轴及其相对于水平轴在逆时针方向的水平倾角α,记录颌骨的质心和拟合椭圆的质心；

(7)确定体位，实施头部校正。

2.根据权利要求1所述的全身CT影像中头部倾斜的自动检测及校正方法，其特征在于：所述的CT影像为2DCT影像。

3.根据权利要求1所述的全身CT影像中头部倾斜的自动检测及校正方法，其特征在于：步骤(2)中，所述的主要身体部位为每层CT影像中的头颈、躯干和腿脚的区域。

4.根据权利要求1所述的全身CT影像中头部倾斜的自动检测及校正方法，其特征在于：步骤(3)中，判断肩部与颈部的交界层进一步包括：检测经步骤(2)处理过的各层CT影像数据中的身体宽度，按从上到下或从下到上的顺序求每连续两层的宽度差值绝对值，以差值的绝对值最大的那层为肩部与颈部分界层。

5.根据权利要求1所述的全身CT影像中头部倾斜的自动检测及校正方法，其特征在于：步骤(7)中，所述的确定体位进一步包括：根据椭圆质心与颌骨质心的相对位置以确定体位是仰卧还是俯卧；若椭圆质心在颌骨质心下方，则判定为仰卧；若椭圆质心在颌骨质心上方，则判定为俯卧。

6.根据权利要求1或5所述的全身CT影像中头部倾斜的自动检测及校正方法，其特征在于：步骤(7)中，所述的实施头部校正进一步包括：对于仰卧状态的头部倾斜校正方法为:若水平倾角α小于90度,则逆时针旋转(90-α)度，若水平倾角α大于90度,则顺时针旋转(α-90)度,最终得到头部垂直向上的校正图像；对于俯卧状态的头部倾斜校正方法为:若水平倾角α小于90度,则逆时针旋转(90-α)度；若水平倾角α大于90度,则顺时针旋转(α-90)度,此时的CT影像头部朝向为垂直向下；再顺时针或逆时针旋转180度，得到头部垂直向上的校正图像。

7.一种计算设备，包括：

一个或多个处理器；

存储器；以及一个或多个程序，其中所述一个或多个程序存储在所述存储器中并被配置为由一个或多个处理器执行，所述的一个或多个程序包括用于执行上述权利要求1-6中任一所述的全身CT影像中头部倾斜的自动检测及校正方法的指令。

8.一种存储一个或多个程序的计算机可读存储介质，所述的一个或多个程序包括指令，所述指令适于由存储器加载并执行上述权利要求1-6中任一所述的全身CT影像中头部倾斜的自动检测及校正方法。

说明书

一种CT影像中头部倾斜的自动检测及校正方法

技术领域

[0001] 本发明属于放疗技术领域，涉及一种CT影像中头部倾斜的自动检测及校正方法、设备和存储介质。

背景技术

[0002] 自动化的计算机断层扫描(CT)医疗影像分析诊断技术通常对于要分析的影像形态有一定的要求。以自动检测头部CT影像中的眼睛算法为例，为降低算法的空间搜索时间和出错概率，通常都是假定被扫描对象是在某个特定的扫描体位的(如仰卧体位)。但在临床实践中被扫描对象经常因为各种原因(如昏迷、医嘱要求等)都是头部歪斜着参加扫描。由此导致的由各个不同的扫描体位得到的CT影像会大大增加影像分析算法的出错概率。为了能够增强影像分析算法的鲁棒性，通常需要有一套能够自适应地检测头部倾斜角度并进行倾斜校正的预处理算法。

[0003] 目前对医学影像做倾斜校正的方法主要服务于医学影像配准，在应用于针对CT影像中头部的倾斜检测及校正时，现有的基于图像几何矩、基于图像惯量矩阵、基于奇异值分解等的医学影像倾斜校正方法均存在计算复杂度高、耗时长、抗噪性较差、数据适应性较差等问题，在实际应用中的性能和鲁棒性不高。

发明内容

[0004] 本发明的目的在于为克服现有技术的缺陷，提供一种全身CT影像中头部倾斜角度自动检测及校正方法。

[0005] 为实现上述目的，本发明采用以下技术方案：

[0006] 本发明根据全身CT影像中头部下颌骨区域与其周边区域的影像的特殊对比及下颌骨区域天然表达的头部朝向信息等特征，提供一种基于CT影像中头部下颌骨朝向的头部倾斜角度自动检测及校正方法。

[0007] 一种全身CT影像中头部倾斜的自动检测及校正方法，适于在计算设备中执行，该方法包括如下步骤：

[0008] (1)对全身CT影像进行二值化操作；

[0009] (2)把步骤(1)中得到的所有二值化影像逐层进行最大连通区域查找，删除所有非最大连通区域，并去除每层CT影像中的非主要身体部位，只保留主要身体部位；

[0010] (3)判断肩部与颈部的交界层；

[0011] (4)确定最佳颌骨层；

[0012] (5)对最佳颌骨层中的颌骨区域进行椭圆拟合；

[0013] (6)根据对最佳颌骨层中的颌骨区域的椭圆拟合，判定椭圆长轴及其相对于水平轴在逆时针方向的水平倾角α,记录颌骨的质心和拟合椭圆的质心；

[0014] (7)确定体位，实施头部校正。

[0015] 本实施例进一步优选地，

[0016] 所述的CT影像为2D CT影像。

[0017] 所述的主要身体部位为每层CT影像中的头颈、躯干和腿脚的区域。

[0018] 步骤(3)中，判断肩部与颈部的交界层进一步包括：检测经步骤(2)处理过的各层CT影像数据中的身体宽度，按从上到下或从下到上的顺序求每连续两层的宽度差值绝对值，以差值的绝对值最大的那层为肩部与颈部分界层。步骤(4)中，确定最佳颌骨层进一步包括：以肩部与颈部的交界层为起始层，向头部方向取得所有包含下颌骨的CT图；对所有包含下颌骨的CT图以颌骨HU阈值做二值化得到颌骨区域，对所有颌骨区域按照最大连通面积进行排序，面积最大者则为最佳颌骨层。

[0019] 步骤(4)中，根据人体肩部与颈部的交界层的位置，向头部方向取N层图像，其中N＝“L(mm)/CT扫描层厚(mm)”的向下取整值，其中L的经验值为200mm足以包含所有下颌骨的CT图

[0020] 步骤(7)中，所述的确定体位进一步包括：根据椭圆质心与颌骨质心的相对位置以确定体位是仰卧还是俯卧；若椭圆质心在颌骨质心下方，则判定为仰卧；若椭圆质心在颌骨质心上方，则判定为俯卧。

[0021] 步骤(7)中，所述的实施头部校正进一步包括：对于仰卧状态的头部倾斜校正方法为:若水平倾角α小于90度,则逆时针旋转(90-α)度，若水平倾角α大于90度,则顺时针旋转(α-90)度,最终得到头部垂直向上的校正图像；对于俯卧状态的头部倾斜校正方法为:若水平倾角α小于90度,则逆时针旋转(90-α)度；若水平倾角α大于90度,则顺时针旋转(α-90)度,此时的CT图像头部朝向为垂直向下；再顺时针或逆时针旋转180度，得到头部垂直向上的校正图像。

[0022] 本发明还提供一种计算设备，包括：

[0023] 一个或多个处理器；

[0024] 存储器；以及一个或多个程序，其中所述一个或多个程序存储在所述存储器中并被配置为由一个或多个处理器执行，所述的一个或多个程序包括用于执行上述全身CT影像中头部倾斜的自动检测及校正方法的指令。

[0025] 本发明还提供一种存储一个或多个程序的计算机可读存储介质，所述的一个或多个程序包括指令，所述指令适于由存储器加载并执行上述全身CT影像中头部倾斜的自动检测及校正方法。

[0026] 本发明具有以下有益效果：

[0027] 1、本发明可以快速有效地检测CT影像中头部的倾角及朝向(朝上还是朝下、朝左还是朝右)，进而为对头部影像的倾斜校正提供准确的参数信息；

[0028] 2、本方法应用于CT头部影像具有计算复杂度低、耗时短、抗噪性高、数据适应好的优点；

[0029] 3、本方法对不同年龄、不同头部形状的人群的CT头部影像都具有良好的适应性；

[0030] 4、本方法可进一步拓展应用到其他自带朝向信息、具有明显组织边界的人体部位用于对不同场景下的人体倾斜校正，具有较好的可拓展性。

附图说明

[0031] 图1为本发明一个实施例中的CT影像中头部倾斜的自动检测及校正方法流程图。

[0032] 图2为本发明一个实施例中对仰卧体位的头部倾斜校正影像处理示意图。

[0033] 图3为本发明一个示例实施例中确定的最佳颌骨层示意图。

[0034] 图4为本发明一个示例实施例中对仰卧体位的最佳颌骨层中的颌骨进行椭圆拟合后的结果示意图。

[0035] 图5为本发明一个示例实施例中对仰卧体位的头部倾斜影像处理校正前、后对比示意图，其中

[0036] (a)为校正前的CT影像；

[0037] (b)为校正后的CT影像。

具体实施方式

[0038] 以下结合实施例和附图进一步说明本发明。

[0039] 一种全身CT影像中头部倾斜的自动检测及校正方法，适于在计算设备中执行，该方法包括如下步骤(流程图如图1所示)：

[0040] 步骤210，对全身CT影像进行二值化操作；

[0041] 本实施例进一步优选地，上述CT影像为2D CT影像；

[0042] 步骤220，把步骤210中得到的所有二值化影像逐层进行最大连通区域查找，删除所有非最大连通区域，并去除每层CT影像中的非主要身体部位，只保留主要身体部位；

[0043] 本实施例进一步优选地，上述主要身体部位为每层CT影像中的头颈、躯干和腿脚的区域，图2为包含下颌骨的CT影像示意图；。

[0044] 步骤230，判断肩部与颈部的交界层；

[0045] 本实施例进一步优选地，判断肩部与颈部的交界层进一步包括：检测经步骤步骤220，处理过的各层CT影像数据中的身体宽度，按从上到下或从下到上的顺序求每连续两层的宽度差值绝对值，以差值的绝对值最大的那层为肩部与颈部分界层；

[0046] 步骤240，确定最佳颌骨层；

[0047] 本实施例进一步优选地，确定最佳颌骨层包括：以肩部与颈部的交界层为起始层，向头部方向取得所有包含下颌骨的CT图；对所有包含下颌骨的CT图以颌骨HU阈值做二值化得到颌骨区域，对所有颌骨区域按照最大连通面积进行排序，面积最大者则为最佳颌骨层；如图3所示为一个示例实施例中确定的最佳颌骨层示意图。

[0048] 步骤250，对最佳颌骨层中的颌骨区域进行椭圆拟合；其中图4为对仰卧体位的颌骨进行椭圆拟合后的结果示意图。

[0049] 步骤260，根据对最佳颌骨层中的颌骨区域的椭圆拟合，判定椭圆长轴及其相对于水平轴在逆时针方向的水平倾角α,记录颌骨的质心和拟合椭圆的质心；

[0050] 步骤270，确定体位，实施头部校正：

[0051] 本实施例进一步优选地，确定体位进一步包括：根据椭圆质心与颌骨质心的相对位置以确定体位是仰卧还是俯卧；若椭圆质心在颌骨质心下方，则判定为仰卧；若椭圆质心在颌骨质心上方，则判定为俯卧；

[0052] 实施头部校正进一步包括：对于仰卧状态的头部倾斜校正方法为:若水平倾角α小于90度,则逆时针旋转(90-α)度，若水平倾角α大于90度,则顺时针旋转(α-90)度,最终得到头部垂直向上的校正图像；对于俯卧状态的头部倾斜校正方法为:若水平倾角α小于90度,则逆时针旋转(90-α)度；若水平倾角α大于90度,则顺时针旋转(α-90)度,此时的CT图像头部朝向为垂直向下；再顺时针或逆时针旋转180度，得到头部垂直向上的校正图像。如图5所示为对仰卧体位的头部倾斜影像处理校正前、后对比示意图。

[0053] 本发明还提供一种计算设备，包括：

[0054] 一个或多个处理器；

[0055] 存储器；以及一个或多个程序，其中所述一个或多个程序存储在所述存储器中并被配置为由一个或多个处理器执行，上述一个或多个程序包括用于执行上述全身CT影像中头部倾斜的自动检测及校正方法的指令，该方法包括如下步骤：

[0056] (1)对全身CT影像进行二值化操作；

[0057] (2)把步骤(1)中得到的所有二值化影像逐层进行最大连通区域查找，删除所有非最大连通区域，并去除每层CT影像中的非主要身体部位，只保留主要身体部位；

[0058] (3)判断肩部与颈部的交界层；

[0059] (4)确定最佳颌骨层；

[0060] (5)对最佳颌骨层中的颌骨区域进行椭圆拟合；

[0061] (6)根据对最佳颌骨层中的颌骨区域的椭圆拟合，判定椭圆长轴及其相对于水平轴在逆时针方向的水平倾角α,记录颌骨的质心和拟合椭圆的质心；

[0062] (7)确定体位，实施头部校正。

[0063] 本发明还提供一种存储一个或多个程序的计算机可读存储介质，所述的一个或多个程序包括指令，上述指令适于由存储器加载并执行上述全身CT影像中头部倾斜的自动检测及校正方法，该方法包括如下步骤：

[0064] (1)对全身CT影像进行二值化操作；

[0065] (2)把步骤(1)中得到的所有二值化影像逐层进行最大连通区域查找，删除所有非最大连通区域，并去除每层CT影像中的非主要身体部位，只保留主要身体部位；

[0066] (3)判断肩部与颈部的交界层；

[0067] (4)确定最佳颌骨层；

[0068] (5)对最佳颌骨层中的颌骨区域进行椭圆拟合；

[0069] (6)根据对最佳颌骨层中的颌骨区域的椭圆拟合，判定椭圆长轴及其相对于水平轴在逆时针方向的水平倾角α,记录颌骨的质心和拟合椭圆的质心；

[0070] (7)确定体位，实施头部校正。

[0071] 应当理解，这里描述的各种技术可结合硬件或软件，或者它们的组合一起实现。从而，本发明的方法和设备，或者本发明的方法和设备的某些方面或部分可采取嵌入有形媒介，例如软盘、CD-ROM、硬盘驱动器或者其它任意机器可读的存储介质中的程序代码(即指令)的形式，其中当程序被载入诸如计算机之类的机器，并被该机器执行时，该机器变成实践本发明的设备。

[0072] 以示例而非限制的方式，计算机可读介质包括计算机存储介质和通信介质。计算机存储介质存储诸如计算机可读指令、数据结构、程序模块或其它数据等信息。通信介质一般以诸如载波或其它传输机制等已调制数据信号来体现计算机可读指令、数据结构、程序模块或其它数据，并且包括任何信息传递介质。以上的任一种的组合也包括在计算机可读介质的范围之内。

[0073] 本领域的技术人员可以理解，可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件，以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。

[0074] 如在此所使用的那样，除非另行规定，使用序数词“第一”、“第二”、“第三”等等来描述普通对象仅仅表示涉及类似对象的不同实例，并且并不意图暗示这样被描述的对象必须具有时间上、空间上、排序方面或者以任意其它方式的给定顺序。

[0075] 此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在下面的权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。

[0076] 上述的对实施例的描述是为便于该技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对这些实施例做出各种修改，并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此，本发明不限于这里的实施例，本领域技术人员根据本发明的揭示，不脱离本发明范畴所做出的改进和修改都应该在本发明的保护范围之内。