一种改进的IS图像分割算法.docVIP

一种改进的IS图像分割算法 　　摘要：由于交互式分割中Live-wire算法运算速度较慢和操作复杂的问题，该文提出了改进的IS分割算法，在其代价函数中，使用Canny算子代替Laplace算子，并加入了像素边界函数、内边界函数和外边界函数等约束条件，较大的提高了算法运算速度的同时也提高了精度。 　　关键词：Live-wire交互式分割；canny算子；IS交互式分割；Dijkstra算法 　　中图分类号：TP18 文献标识码：A文章编号：1009-3044(2011)16-3914-03 　　An Improved Image Segmentation Algorithm 　　WANG Wei-nan1, SONG Jin2 　　(1.Electronic and Information Engineering College, Lanzhou Jiaotong University, Lanzhou 730070; 2.China Electronics Technology Group Corporation Article 13 Institute, Taiyuan 030006, China) 　　Abstract：In view of the problems of the precision of segmentation and that the speed of interactive segmentation algorithm is slow, we have made the improvement to original Live-wire algorithm. The improved algorithm uses the canny-operator to replace the laplace-operator, and uses the automatic segmentation algorithm in ITK to segment excess. It can enhance operating speed, at the same time, can enhance the precision greatly. 　　Key words: Live-wire interactive segmentation algorithms;Canny operator;IS Interactive segmentation algorithms;Dijkstra algorithms 　　1 Live-wire交互式分割的不足 　　Live-wire算法是一种典型半自动交互式分割算法，该交互分割中将人工干预减少到最小程度，且分割精度高，被广泛应用于医学图像处理。这种方法可通过参考图像集的训练得到描述目标边界的数学函数分割，对于与参考图像集类似的医学图像（如相同部位的CT图或MRI图），当人工选定起始点与目标点后，算法会自动定位光标在该区间中的边界。Live-wire算法中的边界检测实际是一个图的搜索问题，即在人工选择的点之间寻找全局的最优路径。 　　本算法有两个不足：第一，用户指定控制点后，动态规划会自动计算这点到其他点的最短路径，运算速度不理想。第二，是分割结果对c2(f)的mean和δ值的依赖性很强，每次分割之前都要进行训练来确定上述值的大小，操作复杂。 　　2 Intelligent Scissors（IS）交互式分割 　　将全相图用f(x,y)来表示，其中x, y为图中各像素点的坐标值，f(x,y)为该点的灰度值，于此同时，x, y和f(x,y)均为有限多个离散值，这样，原图就被转换为矩阵式的数字图像。 　　假设p，q为数字图像中相邻的一对像素，用l(p,q)表示像素p到像素q的局部代价函数，则: 　　式中l(p,q)代表p到其邻接像素q的局部代价；fz,fG,fD,fp,fl,fo分别表示拉普拉斯零点穿越函数、梯度幅度函数、梯度方向函数、像素边界值函数、内临界函数和外临界函数，各组分函数前面的系数w为各组分函数在局部代价函数中的权重。通常情况下wz=0.3, wG==0.3, wD=0.1, wp=0.1, wl=0.1, wo=0.1 　　1) 拉普拉斯零点穿越函数 　　其中，Il(q)为像素点q的拉普拉斯函数值。在实际的计算过程中，将两个相邻像素点的正负号相反拉普拉斯函数值中绝对值较小的那个像素点认为是零点穿越点。 　　需要注意的是，该方法中采用的拉式零点穿越函数将会增强图像中的高频部分，同时可能增加噪点，引入冗余信息，影响后期的处理和判断，本文中将其改进为canny算子。 　　3 Canny算子边缘检测原理 　　在图像与一对称的二维高斯函数做卷