版本16和19间的区别 (跳过第3版)
于2008-06-10 15:40:00修订的的版本16
大小: 2533
编辑: czk
备注:
于2008-06-10 16:08:48修订的的版本19
大小: 2937
编辑: czk
备注:
删除的内容标记成这样。 加入的内容标记成这样。
行号 99: 行号 99:
例子:

{{attachment:figure1009.png}}

峰值检测:
 1. 找到最大值,记下它的位置
 1. 将上一位置周围邻域中的值设为0
 1. 重复前两个步骤,直到检测到所需数量的直线或者低于某一特定的阈值

线检测和链接
 1. 找到小的裂口,将其合并
 1. 找出比最小阈值长的直线段

例子:

{{attachment:figure1010.png}}

点、线、边缘的检测

点、线、边缘的检测常用空域滤波的方法进行。

$$R = w_1z_1+w_2z_2+\dots +w_9z_9=\sum_{i=1}^{9}w_iz_i$$

1. 检测孤立点

检测孤立点采用的滤波模板是

-1 -1 -1
-1  8 -1
-1 -1 -1

检测的标准是判断滤波后的绝对值是否大于某一阈值

$$|R| \ge T$$

在matlab中可以用如下方法实现:

g = abs(imfilter(f, w)) >= T;

另一种检测方法是看一个邻域内的最大值和最小值的差是否超过某个阈值。实现方法是:

g = imsubtract(ordfilt2(f, m*n, ones(m,n)), ordfilt2(f, 1, ones(m,n))) >= T;

figure1002.png

2. 线检测

线检测的模板可以是:

figure1001.png

如果对某一方向的线感兴趣,可以使用相应的模板做滤波后进行阈值处理。

figure1003.png

3. 边缘检测

边缘所满足的条件:

  • 像素值的一阶导数的绝对值大于指定阈值
  • 像素值的二阶导数有零交叉。

figure1004.png

Sobel:

$$g = \sqrt{ [(z_7+2z_8+z_9) - (z_1 + 2z_2 + z_3)]^2 + [(z_3+2z_6+z_9)-(z_1+2z_4+z_7))]^2 }$$

g = edge(f, 'sobel', T, dir)

f是输入图像,T是一个阈值,dir是方向('horizontol'或者'vertical'或者'both')

Prewitt:

g = edge(f, 'prewitt', T, dir) 

Roberts:

g = edge(f, 'roberts', T, dir)

Laplacian of Gaussian:

g = edge(f, 'log', T, sigma)

sigma为标准差,默认是2。

检测零交叉:

g = edge(f, 'zerocross', T, H)

H指定卷积的模板。

Canny:

g = edge(f, 'canny', T, sigma)

figure1005.png

figure1006.png

Hough变换

给定(xi,yi)点,经过该点的直线参数a和b满足:yi=axi+b,a为斜率,b为截距。将经过该点的所有可能的a和b画在ab构成的参数坐标系内,可以得到一条直线。

figure1007.png

原图像上的直线y = ax+b上所有点在参数空间内对应的直线会交在同一点,这点的坐标就是所求的参数值。

由于a可能为无穷大,所以在计算上不方便。我们将直线方程写成参数形式:

$$x\cos\theta + y\sin\theta = \rho $$

原始图像上的一个点,对应于参数空间上的一条曲线:

figure1008.png

例子:

figure1009.png

峰值检测:

  1. 找到最大值,记下它的位置
  2. 将上一位置周围邻域中的值设为0
  3. 重复前两个步骤,直到检测到所需数量的直线或者低于某一特定的阈值

线检测和链接

  1. 找到小的裂口,将其合并
  2. 找出比最小阈值长的直线段

例子:

figure1010.png

阈值处理

区域分割

分水岭算法

图像如何分割 (2008-06-22 15:33:18由czk编辑)

ch3n2k.com | Copyright (c) 2004-2020 czk.