第3章MATLAB图形图像处理 2.pptVIP

? 【例 3.7.3 】 对 liftingbody.png 图像进行旋转，其效果如图 3.36 所示。 ? p = imread(liftingbody.png); ? angle = input(please input angle ： ); % 输入旋转角度为 45 度 ? pr = imrotate(p,angle,bilinear); % 使用双线性插值方法实 现旋转变换 ? subplot(1,2,1); ? imshow(p); ? subplot(1,2,2); ? imshow(pr); 3.8 图像边缘检测与分割 ? 图像分割 是将数字图像分割成互不相交 ( 不 重叠 ) 的有意义的子区域的过程。其目的是 使各个区域与景物中以某种方式描述的物 体相对应，以便进行高层次的 图像解释、 图像识别 等处理。 3.8.1 边缘检测概述 ? 图像分割 一般是基于图像像素灰度值的 基本特性— —不连续性与相似性中的一个或两个来进行的 。 区 域内部的像素一般具有灰度相似性，而在区域之间 的边界上一般具有灰度不连续性。 ? 根据图像灰度级的突变进行图像分割的方法中， 边 缘检测方法 通过检测包含不同区域的边缘来解决图 像分割问题。边缘检测最通用的方法是 检测灰度值 的不连续性 ，灰度值的不连续是指在不同区域之间 的边缘上像素灰度值的变化往往比较剧烈，一般利 用图像一阶导数的极大值或者二阶导数的过零点信 息提供判断边缘点的依据。 3.8.2 梯度算子 ? 对于图像中任一像素点，可以用二维函数 f(x,y) 来表示，而 梯度定义为向量 ? f ： ? 该向量的幅值为： ? 为简化计算，通常省略开方或通过取绝对值来近似计算， 即 ? 或 ? ? T T y f x f Gy Gx y x f ? ? ? ? ? ? ? ? ? ? ? ? ? ) , ( ? ? 2 1 2 2 2 1 2 2 ) ( ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y f x f G G f mag f y x 2 2 ) ( y x G G f mag f ? ? ? ? ? y x G G f ? ? ? ? 梯度向量的含义在于，它总是指向 f(x,y) 在点 (x,y) 处的最大 变化率方向。最大变化率是用方向 α 角来衡量的，即 ? 为了能估计出 Gx 及 Gy 值，通常使用一些经典的模板来做数 字化近似。这些模板有 Soble 模板、 Prewitt 模板、 Roberts 模 板 等。 ? MATLAB 图像处理工具箱中提供了 边缘函数（ edge ） 实现 对图像边缘的检测，其基本语法格式为： ? [e , s] = edge(p, method, para) ? 其中 p 为输入图像， method 是边缘检测的类型， para 为与 edge 对应的参数， e 为与 p 同样大小的逻辑矩阵，在检测到 边缘的位置时值为 1 ，其它位置则为 0 ， s 为一可选参数。 ) / arctan( ) , ( x y G G y x ? ? 1 ． Roberts 模板 ( 算子 ) ? Roberts 算子又称为 Roberts 交叉算子，它是由 Roberts 提出 的一种利用局部差分算子寻找边缘的模板，它在 2 × 2 邻域 上计算对角导数。在实际应用中，为了简化计算，常用梯 度函数的 Roberts 绝对值来近似。另外还可以用 Roberts 最 大值算子来计算。 Roberts 边缘检测算子如图所示： ? (a)Roberts-x (b)Roberts-y -1 0 0 1 0 -1 1 0 ? Roberts 边缘检测器的语法格式如下： ? [e , s] = edge(p, roberts, thresh, direct) ? 其中 p 为输入图像， thresh 为指定的阈值 T ，若 T 值 未指定，则函数 edge 会自动选择一个值。 direct 为 检测边缘的首选方向，通常可选作 horizontal( 水平 ) 、 vertical （垂直）或 both （默认值）。 2 ． Sobel 模板 ( 算子 ) ? Roberts 算子的一个主要问题是计算方向差分时对噪声敏 感。 Sobel 提出一种将方向差分运算与局部平均相结合的 方法，即 Sobel 算子。该算子是在以 f(x,y) 为中心的 3 × 3 邻域上计算 x 和 y 方向的偏导数，其模板如图所示： ? (c) Sobel-x (d) S