基于matlab图像处理教程系列之图像分割（一）

dapmood

5 O. l/ @# M: a8 d6 u# ?6 ~; r8 W" M
一、简介
1. 点检测
  s! G4 m, I0 j" s1 }g = abs(imfilter(tofloat(f), w)) >= T;

; v* E, j# d% d2 A+ T0 h
& t! r- d" v/ \/ r5 O$ `$ `
f 是输入图像，w 是适合点检测的模板:
0 A2 a8 b- M3 w* k1 Y, D
，g 是包含检测点的图像。imfilter 把输出转换为输入所属的类，如果输入是整数类，并且 abs 操作不接受整数数据，那么在滤波操作中用 tofloat(f)来防止对数值的过早截取。输出图像 g 是 logical 类；值是 0 和 1。如果 T 值没有给出，那么通常基于滤波结果来选取。在那种情况下，先前的一串指令分成三个基本步骤：1) 计算滤波后的图像 abs(imfilter(tofloat(f),w))；2) 从滤波后的图像的数据中找出 T 的值；3) 把滤波后的图像与T做比较。

2.线检测
: ^5 \! C4 n# C% Eg = imfilter(tofloat(f), w);
其中w为线检测模板：
3 D; V+ [5 x( [8 P4 l5 q


( l! h& e( e4 s  {1 X如果图中的第一个模板在图像上移动，就会对水平线(一个像素宽)的响应更强烈。对于恒定的背景，当线通过模板的中间一行时可能产生更大的响应。同样，图 中的第 2 个模板对+45°线响应最好，第 3 个模板对垂直线响应最好，第 4 个模板对 –45°线响应最好。注意，每个模板的优先方向都用比其他可能方向要大的系数加权。每个模板的系数之和为 0，这表明在恒定亮度区域中，模板的响应为 0。若每方向模板都应用于同一图像，其中一方向模板在图像中心点响应比其他模板在图像中心点响应都大，我们就说那个点与模板 i 方向的线更相关。

• =imread('gantrycrane.png');
• I=rgb2gray(I);
• h1=[-1, -1. -1; 2, 2, 2; -1, -1, -1];
• h2=[-1, -1, 2; -1, 2, -1; 2, -1, -1];
• h3=[-1, 2, -1; -1, 2, -1; -1, 2, -1];
• h4=[2, -1, -1; -1, 2, -1; -1, -1, 2];
• J1=imfilter(I, h1);
• J2=imfilter(I, h2);
• J3=imfilter(I, h3);
• J4=imfilter(I, h4);
• J=J1+J2+J3+J4;
• h=figure;
• set(h,'color','white')
• subplot(121); imshow(I);
• subplot(122); imshow(J);
  / D! }- z% @" I0 G6 _% y3 ]

1 M1 r% D- D9 o

/ D2 [3 Z# X+ N! l  q
3.边缘检测
  ~4 d2 \7 D/ V6 v7 |0 K
! R! r: O7 B+ f) E# T4 F1 i[g，t]=edge(f, ‘method’, parameters)
2 Z# M& ?4 Y( U
f 是输入图像，method 是以下列出方法中的一种，parameters 是下边说明
( B+ [, M& [' V9 V5 T
8 N4 b& _% o% r的附加参数。输出 g 是在 f 中被检测到的边缘点的位置为 1，而在其他地方为 0 的逻辑数组。

参数 t 是可选择的；由 edge 给出阈值，以决定哪个梯度值足够大到可以被称作边缘点。

（1）Sobel 边缘检测算子
* h% A" ?9 b8 {: f+ F$ {) P
Sobel 检测算子的一般调用语法是：[g,t] = edge(f, ‘sobel’, T, dir);
% k0 w% i/ c; A- c; n' Q7 A
; G- s9 M# b. P8 |其中 f 是输入图像，T 是指定的阈值，dir 是指定的检测边缘的首选方向：’horizontal’、 ‘vertical’或’both’(默认值)。

) |: D; u$ j6 e5 J（2）Prewitt 边缘检测算子

一般调用语法是：[g,t] = edge(f, ‘prewitt’, T, dir);

# A$ d0 \6 a- w% L7 C- N( s0 C4 G这个函数的参数和 Sobel 参数相同。Prewitt 检测算子相比 Sobel 检测算子在计算上要简单一点，但是比较容易产生噪声。
2 N# A1 m* k/ a- N5 O/ b
（3）Roberts 边缘检测算子

一般调用语法是：[g,t] = edge(f, ‘roberts’, T, dir);

; Z- c: ^+ B" @: e这个函数的参数和 Sobel 参数相同。Roberts 检测算子是数字图像处理中最古老的边缘检测算子中的一种，并且也是最简单的一种。因为在部分功能上有限制，所以这种检测算子的使用明显少于其他几种算子(比如，Roberts 检测算子是非对称的，而且不能检测诸如 45°倍数的边缘)。然而，在简单和速度为主导因素的情况下，Roberts 检测算子还是经常用在硬件实现方面。

4 Z4 T5 P. F% M& f) O) x（4） LoG 检测算子

, n/ H" y" d6 S% RLoG 检测算子的一般调用语法是：[g,t ]= edge(f, ‘log’, T, sigma);

9 G# A6 W! l( R( W其中的 sigma 是标准差，其他参数和前边解释过的一样。sigma 的默认值是 2。和以前一样，edge 忽略一切不比 T 强的边缘。如果 T 值没有给出或为空[ ]，edge 会自动选值。将 T置为 0，将产生封闭的轮廓，这是我们熟知的 LoG 方法的典型特征。
3 K; C8 S# `) P$ z! |- E7 B8 q
（5）零交叉检测算子

& @! J. K' D+ x5 O# @2 O这种检测算子基于的概念与 LoG 方法相同，但是卷积使用特殊的滤波函数 H 来完成，调用语法为：[g,t] = edge(f, ‘zerocross’, T, H);
4 Z  T! a! n$ n) z2 d
其他参数和 LoG 解释的一样。
- x6 l  I6 b" R( E/ X
) I4 ]( {* f1 w& j（6）Canny 检测算子
8 A& i2 y2 o; K/ S( D( Y
5 g. H& F/ f5 [6 r% ICanny检测算子(Canny[1986])是 edge 函数中最强的边缘检测算子。Canny检测算子的语法是：[g,t] = edge(f, ‘canny’, T, sigma);

在这里，T 是向量。T=[T 1 ，T 2 ]，包含在前边步骤 3)的两个阈值，sigma 是平滑滤波器的标准差。如果 t 包括在输出参量中，t 就是二元矢量，其中包含该算法用到的两个阈值。语法中的其余参数和其他方法中解释的一样，包括：如果 T 没有指定，就自动计算阈值。sigma 的默认值是 1。
! S, U* t/ N: Z$ m$ b7 Z" a
( k8 |. }) M  @( h0 {* Q+ E
二、源代码

• I=imread('cameraman.tif');
• I=im2double(I);
• [J, thresh]=edge(I, 'prewitt', [], 'both');
• h=figure;
• set(h,'color','w')
• subplot(121); imshow(I);
• subplot(122); imshow(J);
  5 b+ Q0 g- m' t- \% y; w1 j' k% r

* P: \. U6 ]. h! U


7 u/ X5 ]/ K' v7 j" k9 I

youOK

基于matlab图像处理教程系列之图像分割（一）

Jame33

g 是包含检测点的图像。imfilter 把输出转换为输入所属的类，如果输入是整数类，并且 abs 操作不接受整数数据，那么在滤波操作中用 tofloat(f)来防止对数值的过早截取