[发明专利]基于交互式区域生长的胃部CT序列图像分割方法

权利要求书

1.一种基于交互式区域生长的胃部CT序列图像分割方法，包括如下步骤：

(1)输入胃部CT序列图像，对序列中每幅图像外围的脂肪组织区域用给定的掩膜图像减掉，去掉边缘的冗余部分，避免产生误分割；

(2)利用交互式区域生长对CT序列中的第一幅图像I进行分割；

(3)计算当前已分割得到的图像的目标区域的灰度均值avegray和灰度方差square。将已分割好的图像目标区域的中心及此中心的8邻域点投影到序列中与已分割好的图像相邻的下一幅图像II的相应位置，并将投影过来的点依次与avegray作差，将差值小于灰度方差square的点作为该幅图像的种子点，设投影完经过筛选后这些种子点的个数为kk；

(4)基于已选取的种子点再次利用交互式区域生长对序列中的当前图像II进行分割，得到该幅图像的分割结果；

(5)判断CT序列中的图像是否已全部分割完成，否则转步骤(3)，直到所有图像分割完成。

2.根据权利要求1所述的方法，其中步骤(2)所述的利用交互式区域生长对CT序列中的第一幅图像I进行分割，按如下步骤进行进行：

(2a)先在CT序列的第一幅图像I的待分割区域内人工选择k个种子点，设一个与图像I大小相同的存放分割结果的矩阵Y，再给定一个相似性阈值T，k取大于0的整数；

(2b)将图像I中的每个像素值依次与k个种子的灰度值作差，则Y中每个位置处将得到k个差值：Y(i，j，k)＝abs(I(i，j)-seed(k))，其中，I(i，j)表示图像I中第i行第j列处的像素值，seed(k)是第k个种子点的像素值，Y(i，j，k)表示图像I中第i行第j列处的第k个差值，abs表示求绝对值运算；

(2c)将k个差值Y(i，j，k)中小于或等于阈值T的值都置为1，大于阈值T的置为0，再继续对得到的k个结果进行或运算，得到矩阵Y中第i行第j列处的值；

Y(i，j)＝Y(i，j，1)|Y(i，j，2)|…|Y(i，j，k)； 

(2d)对矩阵Y进行形态学重构：将矩阵Y中值为1的各个点集合起来形成分割结果区域，最后进行区域标记得到分割结果。

3.根据权利要求1所述的方法，其中所述步骤(3)中计算当前已分割得到的图像的目标区域的灰度均值avegray是通过公式：avegray＝1/n∑I(p，q)计算，其中，n表示分割得到的目标区域内像素点的个数，p表示目标区域内每个点所在的行，q表示目标区域内每个点所在的列，I(p，q)表示当前已分割好的图像的目标区域内位于第p行第q列处的像素值。

4.根据权利要求1所述的方法，其中所述步骤(3)中计算当前已分割得到的图像的目标区域的灰度方差square是通过公式计算：

。

5.根据权利要求1所述的方法，其中步骤(4)所述的基于已选取的种子点再次利用交互式区域生长对序列中的当前图像II进行分割，按如下步骤进行：

(4a)将图像II中的每个像素值依次与kk个种子的灰度值作差，则矩阵Y中每个位置处将得到kk个差值：Y(i，j，kk)＝abs(II(i，j)-seed(kk))，

其中，II(i，j)表示图像II中第i行第j列处的像素值，seed(kk)是第kk个种子点的像素值，Y(i，j，kk)表示图像II中第i行第j列处的第kk个差值，abs表示求绝对值运算；

(4b)将kk个差值Y(i，j，kk)中小于阈值T的值均置为1，大于阈值T的值置为0，再继续对得到的kk个结果进行或运算，得到矩阵Y中第i行第j列处的值；

Y(i，j)＝Y(i，j，1)|Y(i，j，2)|…|Y(i，j，kk)

(4c)对矩阵Y进行形态学重构：将矩阵Y中值为1的各个点集合起来形成分割结果区域，最后进行区域标记得到分割结果。