第七章_医学图像分割.pptVIP

第七章　医学图像分割 韩贵来 象素分类 象素分类 医学图像分割 阈值分割法——基于区域的并行分割算法 区域分割法——基于区域的串行分割算法 边缘分割法——基于边界的并行分割算法 边缘跟踪与边界拟合——基于边界的串行分割算法 阈值分割 全阈值分割结果 全阈值分割的MATLAB实现 I=imread(source.bmp); figure(1) imhist(I); T＝110;%=80; p=find(I=T); I(p)=0; p=find(IT); I(p)=255; figure(2) imshow(I); 全阈值分割结果 灰度级门限化 灰度级门限化 基本公式 条件(1)表明分割区域要覆盖整个图像且各区域互不重叠； 条件(2)表明每个区域具有相同的性质； 条件(3)表明相邻的两个区域性质相异不能合并为一个区域。 区域生长----分割区域 区域生长----分割区域 区域生长----分割区域 区域生长 区域生长 7.3.2 区域分裂合并 7.3.2 区域分裂合并 7.3.2 区域分裂合并 7.3.2 区域分裂合并 理想的边缘 斜坡数字边缘 边缘及检测原理 方向模板 Kirsch边缘算子 Kirsch边缘算子 Kirsch边缘检测算子需要求出f(x，y)8个方向的平均差分的最大值，计算8个方向模板，该算子可以检测出边缘的方向性信息，并能较好地抑制边缘检测的噪声。 Kirsch算子边缘检测结果 拉普拉斯算子 拉普拉斯算子 从以上几个图片可见： 1、Laplacian算子虽对边缘有响应，但对拐角、线条、线端点和孤立点响应更强。即对noise 增强作用更大 2、8的比4的增强作用大 Cany算子 梯度算子边缘检测的MATLAB实现 I = imread(‘cameraman.tif); BW1 = edge(I,sobel); %应用Sobel算子进行滤波 BW2 = edge(I,roberts); %应用Roberts算子进行滤波 BW3 = edge(I,prewitt); %应用Prewitt算子进行滤波 BW4 = edge(I,log); %应用LOG算子进行滤波 subplot(2,3,1),imshow(BW1) subplot(2,3,2),imshow(BW2) subplot(2,3,3),imshow(BW3) subplot(2,3,4),imshow(BW4) 7.5 边界跟踪 7.5 边界跟踪 7.5 边界跟踪 图像分割:边缘连接法 Hough变换 问题的提出 Hough变换的基本思想 算法实现 Hough变换的扩展 图像分割:边缘连接法 Hough变换问题的提出 在找出边界点集之后，需要连接，形成完整的边界图形描述 图像分割:边缘连接法 Hough变换的基本思想 对于边界上的n个点的点集，找出共线的点集和直线方程。 对于任意两点的直线方程：y = ax + b，构造一个参数a，b的平面，从而有如下结论： 图像分割:边缘连接法 Hough变换的基本思想 xy平面上的任意一条直线y = ax + b ，对应在参数ab平面上都有一个点 过xy平面一个点(x,y)的所有直线，构成参数ab平面上的一条直线。 图像分割:边缘连接法 Hough变换的基本思想 如果点(x1,y1)与点(x2,y2)共线，那么这两点在参数ab平面上的直线将有一个交点 在参数ab平面上相交直线最多的点，对应的xy平面上的直线就是我们的解 图像分割:边缘连接法 Hough变换的基本思想 图像分割:边缘连接法 Hough变换算法实现 由于垂直直线a,为无穷大，我们改用极坐标形式： xcos ? + ysin? = ? 参数平面为?,? ，对应不是直线而是正弦曲线 使用交点累加器，或交点统计直方图，找出相交线段最多的参数空间的点 然后找出该点对应的xy平面的直线线段 图像分割:边缘连接法 Hough变换的扩展 Hough变换不只对直线，也可以用于圆： （x – c1）2 + (y - c2)2 = c32 这时需要三个参数的参数空间。 曲面拟合法 微分算子进行边缘提取的共同缺陷：提升噪声。 曲面拟合法可以较好地抑制噪声。 基本思路：用一个平面或曲面去逼近一个图像面积元，然后利用这个平面或曲面的梯度代替点的梯度，从而实现边缘检测。 曲面拟合法 一次平面拟合 　如下图构成图像面积元： 　 　 　现 用 一 次 平 面 ax + by + c 去拟合该面积元 ?s 上的 4 个相邻像素，即用 p(x, y) = ax + by + c 去逼近f (x,y) 。 　　 已知 f (x, y) 与 p (x, y) 之间的均方差为 ： 　按均方差最