[发明专利]自适应窗的胃部CT图像淋巴结跟踪检测系统及方法

权利要求书

1.自适应窗的胃部CT图像淋巴结跟踪检测系统，其特征是：依次包括用于对图像进行分类的基于模糊C均值的分类模块、用于对感兴趣区域进行外边界检测的基于边界先验的疑似淋巴结提取模块和用于对疑似淋巴结的窗口图像中的目标进行跟踪并检测出目标的自适应窗的淋巴结跟踪检测模块；

所述基于模糊C均值的分类模块依次包括用于求解模糊C均值目标函数的模糊C均值参数计算子模块和对计算结果样本分类的模糊C均值分类子模块；

所述基于边界先验的疑似淋巴结提取模块依次包括对图像各区域进行边界提取的边界提取子模块和对疑似淋巴结边界进行标记的疑似淋巴结提取子模块；

所述自适应窗的淋巴结跟踪检测模块依次包括疑似淋巴结面积排序子模块、疑似淋巴结窗口大小及中心位置计算子模块、血管跟踪子模块、淋巴结跟踪子模块、血管标记子模块和淋巴结标记子模块。

2.一种适用于权利要求1所述的跟踪检测系统的方法，其特征是：包括以下几个步骤：

第一步，基于模糊C均值的分类

用于对图像进行分类，分为3类区域，第一类为背景区域，第二类是脂肪区域，第三类为感兴趣区域，主要步骤为：

1.1、模糊C均值参数计算子步骤，用于求解模糊C均值目标函数的参数，其中v₁,…,v_c是c个聚类中心，m>1是个常量，μ_ik是隶属度，是样本x_k属于第i类的概率，并且0≤μ_ik≤1，U是c·n的隶属度矩阵；

1.2、模糊C均值分类子步骤，通过模糊C均值参数计算子模块，可以得到最终的参数，通过公式得到样本xj被分到第k类。其中c=3，可以把样本分为3类，第一类为背景区域中的样本，第二类是脂肪区域中的样本，第三类为感兴趣区域中的样本；

第二步，基于边界先验的疑似淋巴结提取

用于对感兴趣区域进行外边界检测，利用先验疑似淋巴结的外边界在脂肪区域中，能够提取出疑似淋巴结，主要步骤为：

2.1、边界提取子步骤，用于对感兴趣区域图中各个连通区域进行外边界提取，对于感兴趣区域图中每个像素的8领域中是脂肪区域和背景区域的像素标记为边界；

2.2、疑似淋巴结提取子步骤，利用疑似淋巴结的外边界在脂肪区域中这个先验，若感兴趣区域图的一个连通区域的外边界全都在脂肪区域中，则标记该区域为疑似淋巴结；

第三步，自适应窗的淋巴结跟踪检测

用于对疑似淋巴结的窗口图像中的目标进行跟踪并检测出目标，该模块按每个疑似淋巴结最大切片面积自适应定义每个疑似淋巴结的窗口并把最大切片面积所在CT图像作为当前帧；从当前帧开始向前和向后进行跟踪操作，由于淋巴结和血管在跟踪过程中其面积变化和中心偏移是有区别的，从而检测出淋巴结，主要步骤为：

3.1疑似淋巴结面积排序子步骤，通过基于边界先验的疑似淋巴结提取步骤，得到所有图像中的疑似淋巴结，对所有图像中的疑似淋巴结按照面积从大到小进行排序，得到A＝{a₁,a₂,…,a_n}，其中Area(a_k)≥Area(a_k+1)≥5,1≤k≤n-1.Area(a_k)是疑似淋巴结a_k的面积，即疑似淋巴结a_k所在区域的像素个数，n为所有图像的疑似淋巴结的总个数；

3.2、疑似淋巴结窗口大小及中心位置计算子步骤，对于疑似淋巴结a_k∈A的窗口大小为(2*l_k+1)·(2*l_k+1)，其中l_k的定义为：找到疑似淋巴结a_k中所有像素点到中心点(x,y)之间具有最远距离的像素点，该像素点与中心点(x,y)的距离即为l_k，中心点是疑似淋巴结a_k中所有像素的均值位置，即其中(x_j,y_j)∈K(a_k)，K(a_k)是疑似淋巴结a_k中的所有像素的位置集合，m为集合K(a_k)中包含像素总数；以疑似淋巴结a_k的中心像素点(x,y)为中心，大小为(2*l_k+1)·(2*l_k+1)的窗口作为a_k的窗W_f，其中f是疑似淋巴结a_k的窗口所在的帧；

3.3、血管跟踪子步骤，从当前帧开始向前向后进行血管的跟踪，下一帧所能正确跟踪出的血管是由上一帧的血管在上一帧窗口的相应位置所决定的，由于血管的面积变化不大，且两帧之间的位移不大，所以在下一帧的相应位置可以跟踪到血管，我们分为三种情况，如果在相应位置上有多个不同的连通区域，并且在相应位置所含像素最多的那个连通区域的面积和当前帧血管的面积变化不大，那么我们选择在相应位置所含像素最多的那个连通区域即为跟踪到的目标；如果在相应位置上只有一个连通区域，并且该连通区域在该窗口中的面积与当前帧血管的面积变化不大，那么该连通区域即为跟踪到的目标；如果在相应位置上没有连通区域存在，这种情况是由于血管突然转走，以后不再出现该血管区域，那么血管跟踪算法结束；如果在下一帧中可以跟踪到目标，那么我们首先取一个比原来窗口更大的新的窗口W，该新窗口W的边界线上不包含该跟踪到的目标在W中对应的目标中的像素，如果在下一帧跟踪到的目标区域和新的窗口W对应的目标区域的面积不相同，那么我们需要做扩窗操作，重新得到的窗口大小是由新的窗口W对应的目标区域中各个像素与中心点之间的最大距离得到的；

3.4、淋巴结跟踪子步骤，从当前帧开始向前向后在分割后的感兴趣区域图上的相应位置各取相应大小的15个窗，该相应位置和相应大小与当前帧相同，得到31个窗口图像，去除这31个窗口图像的多余冗余信息，既每一帧最多保留一个跟踪区域，我们所关注的区域是疑似淋巴结a_k所在窗口中的区域位置O，对于某一帧，若有多个疑似淋巴结都在区域位置O中有像素，那么跟踪结果是在区域位置O中有像素的疑似淋巴结的中心距离窗口中心最近的那个疑似淋巴结；若只有一个疑似淋巴结在区域位置O中有像素，则跟踪结果就是该疑似淋巴结；若没有疑似淋巴结在区域位置O中有像素，则跟踪结果为空；其它窗口图像也是按照此种方法得到跟踪结果，从而可以去除多余信息；接下来，计算所有窗口图像到当前帧窗口图像的相似度，对于一个窗口图像的特征为窗口中所有灰度值非零像素在直方图的十个区间中的像素个数所组成的十维向量，即对于窗口图像W_f，设该窗口图像中最大的像素灰度值为d＝255，最小的像素灰度值为c＝0，十个区间为则窗口图像f的每个灰度值非零像素一定落在这十个区间中的一个，可以得到一个十维向量既可以得到每个窗口图像的特征向量分别为计算每个窗口的特征向量与疑似淋巴结a_k所在窗口图像f的特征的相似度为得到一个长度为31维的相似性向量s=(||b→f-15-b→f||2||b→f||2,||b→f-14-b→f||2||b→f||2,...,||b→f+15-b→f||2||b→f||2);]]>现在开始找该疑似淋巴结的开始帧和结束帧，从s_f开始，向前找到第一个大于t₁的相似度所在的帧为开始帧low_f，我们找到中的最大值，如果最大值大于则更改开始帧low_f为具有该最大值的帧中离f最近的那个最大值所在的帧；找结束帧的方法和开始帧相似，从s_f开始，向后找到第一个大于t₁的相似度所在的帧为结束帧high_f，我们找到中的最大值，如果最大值大于则更改结束帧high_f为具有该最大值的帧中离f最近的那个最大值所在的帧；

3.5、血管标记子步骤，通过血管跟踪子步骤，得到该疑似淋巴结a_k的跟踪窗口序列图，若窗口的目标面积变化不大，且最终得到的窗口比原来的窗口大2倍，则判断该疑似淋巴结a_k为血管，否则不是血管；如果疑似淋巴结a_k是血管，则a_k的窗口序列中的血管有同一个标记，该标记不同于其它序列，从A中删除血管a_k的窗口序列图中对应的血管；

3.6、淋巴结标记子步骤，通过淋巴结跟踪子步骤可以得到该疑似淋巴结a_k的窗口序列图，若疑似淋巴结a_k窗口序列图中的目标的面积变化程度大于阈值t₂，则该疑似淋巴结a_k是淋巴结，否则不是淋巴结；若疑似淋巴结a_k是淋巴结，则疑似淋巴结a_k的窗口序列图中的淋巴结有一个相同的标记，且该标记不同于其它序列，并从A中删除淋巴结a_k的窗口序列图中对应的淋巴结。