第十五节、韦伯局部描述符(WLD，附源码)

一 WLD原理
二代码实现
三代码下载

纹理作为一种重要的视觉线索，是图像中普遍存在而又难以描述的特征，图像的纹理特征一般是指图像上地物重复排列造成的灰度值有规则的分布。纹理特征的关键在于纹理特征的提取方法。目前，用于纹理特征提取的方法有很多，最具有代表性的是有基于二阶概率密度的灰度共生矩阵、符合人眼视觉特性的小波变换、纹理谱法以及基于图像结构基元的纹理元方法等。
为了更有效地描述图像局部纹理特征，又先后提取了局部二值模式(Local Binary Pattern,LBP)和韦伯局部描述符(Weber Local Descriptor,WLD),这两种方法都是基于邻域像素点间的灰度变化特征来描述图像纹理的。因两者易于理解、便于计算且具有较好的局部特征描述能力而被广泛地应用于纹理分类、目标检测、人脸识别、合成孔径雷达(Synthetic Aperture Radar,SAR) 图像索引、指纹活性检测、图像伪造检测等众多领域中。

回到顶部

一 WLD原理

韦伯定律是反映心理量和物理量之间关系的定律，它表明能够引起感觉差异的差别阈值与原始刺激的强度之比是一个常量，即：

$\frac{ΔI}{I}=k$

其中 $k$ 是一个常量； $ΔI$ 表示差别阈值； $I$ 表示原始刺激的强度。由此可以推知，刺激的变化所引起的感觉差异不仅与刺激变化的大小有关，还与原始刺激的强度有关。局部图像描述符WLB就是根据该定律提出的，它包含两个算子：差分激励算子和方向算子。WLD计算除边缘像素点外的每个像素点的差分激励和方向，并以其二维分布直方图来联合表征图像的纹理特征。

1、差分激励算子

差分激励反映局部窗内灰度变化的强度信息。通过计算局部窗内邻域像素点与中心像素点间的灰度差值和与中心像素点灰度值的比值 $G_{ratio}(x_c)$ ，再利用反正切变换将分布在 $[-P,+∞]$ 范围内的 $G_{ratio}(x_c)$ 映射到区间 $(-\frac{\pi}{2},\frac{\pi}{2})$ 内得到，差分激励 $ξ(x_c)$ 的计算式为：

$ξ(x_c)=arctan(G_{ratio}(x_c))=arctan(\sum\limits_{i=0}^{P-1}\frac{(x_i-x_c)}{x_c})$

其中, $x_c$ 和 $x_i,i=0,1,...,P-1$ 分别表示中心像素点和邻域像素点的灰度值； $P$ 表示邻域像素点个数。

2、方向算子

方向反映局部窗内灰度变化的空间分布信息。通过局部窗内水平方向与垂直方向上邻域像素点的灰度差值比值的反正切变换来描述。其计算式为:

$Φ(x_c)=arctan(\frac{D_V}{D_H})$

其中, $D_H$ 和 $D_V$ 分别表示水平方向上和垂直方向上中心像素点两侧的邻域像素点间的灰度差异，若对于 $3\times3$ 像素的局部窗口如下图所示，则 $D_H=x_7-x_3$ ， $D_V=x_5-x_1$ 。

为了能够更加有效的区分局部窗口的灰度分布变换，进一步将方向由 $Φ(x_c)\in(-\frac{\pi}{2},\frac{\pi}{2})$ 变换到了 $Φ'(x_c)\in[0,2\pi]$ ，其变换公式为：

$Φ'(x_c)=\begin{cases} Φ(x_c) \qquad D_H>0,D_V>0 \\ \pi+Φ(x_c) \qquad D_H<0,D_V>0 \\\pi+Φ(x_c) \qquad D_H<0,D_V<0 \\ 2\pi+Φ(x_c) \qquad D_H>0,D_V<0 \end{cases}$

3、WLD直方图

WLD采用均匀量化技术，将方向 $Φ'(x_c)$ 均匀地量化为 $T$ 个方向，将差分激励均匀地划分为 $M$ 个频段，分别对应于图像中的高频、中频和低频变化，再将划分的每个频段上将差分激励均匀地量化为 $S$ 格，形成一个 $T\times{C}=T\times(M\times{S})$ 的二维直方图，并通过编码将其转化为一维向量用于表示图像的纹理特征。

二维直方图如上图所示，横轴表示方向，纵轴表示差分激励。每个小矩形表示在该方向下所在差分激励区间像素的数量，数量不同，颜色不同。

下图是WLD直方图得到的过程：

1、在每个主方向上得到差分激励子直方图，得到 $H(0)$ 至 $H(T-1)$ ;

2、将 $H(k)$ 分成 $M$ 个子区间，即 $l_m,m=0,...,M-1$ ,将 $H(k)$ 中的 ${l_m}$ 对应放置在 $t=k$ 和 $m=i$ 处；

3、将 $m=i$ 的一行子直方图拼接成一个直方图，即 $H_i$ ；

4、将 $H_i$ 组合成一个直方图，即WLD直方图。

然后可以把WLD直方图用于后续分类。

4、缺点

我们来看下图所示的3个局部灰度分布示例。从3个局部窗内的灰度分布看，(a)~(c)分别表示了高频、中频和低频3种变换。按照WLD的计算方法，他们的 $ΔI$ 都等于0，即差分激励等于0，而且，在垂直方向上的灰度值差值也为0，表明方向也等于0.这意味着，WLD将无法区分这3个纹理模式。

存在上述问题的主要原因有:

WLD的差分激励算子是利用各向同性的边缘检测滤波器———拉普拉斯算子统计局部窗内 $P$ 个邻域像素点与中心像素点间的灰度差值之和 $ΔI$ ，导致了灰度变化的正负差值相互抵消,换言之,局部窗内的灰度变化信息没有充分体现;
WLD 的方向算子仅表达了水平方向和垂直方向上邻域像素点间灰度变化梯度的空间分布方位,不能充分反映局部窗内灰度变化的空间分布结构信息,难以体现纹理的内在变化特征;

5、改进

亲爱的读者和支持者们，自动博客加入了打赏功能，陆陆续续收到了各位老铁的打赏。在此，我想由衷地感谢每一位对我们博客的支持和打赏。你们的慷慨与支持，是我们前行的动力与源泉。

日期	姓名	金额
2023-09-06	*源	19
2023-09-11	*朝科	88
2023-09-21	*号	5
2023-09-16	*真	60
2023-10-26	*通	9.9
2023-11-04	*慎	0.66
2023-11-24	*恩	0.01
2023-12-30	I*B	1
2024-01-28	*兴	20
2024-02-01	QYing	20
2024-02-11	*督	6
2024-02-18	一*x	1
2024-02-20	c*l	18.88
2024-01-01	*I	5
2024-04-08	*程	150
2024-04-18	*超	20
2024-04-26	.*V	30
2024-05-08	D*W	5
2024-05-29	*辉	20
2024-05-30	*雄	10
2024-06-08	*:	10
2024-06-23	小狮子	666
2024-06-28	*s	6.66
2024-06-29	*炼	1
2024-06-30	*!	1
2024-07-08	*方	20
2024-07-18	A*1	6.66
2024-07-31	*北	12
2024-08-13	*基	1
2024-08-23	n*s	2
2024-09-02	*源	50
2024-09-04	*J	2
2024-09-06	*强	8.8
2024-09-09	*波	1
2024-09-10	*口	1
2024-09-10	*波	1
2024-09-12	*波	10
2024-09-18	*明	1.68
2024-09-26	B*h	10
2024-09-30	岁	10
2024-10-02	M*i	1
2024-10-14	*朋	10
2024-10-22	*海	10
2024-10-23	*南	10
2024-10-26	*节	6.66
2024-10-27	*o	5
2024-10-28	W*F	6.66
2024-10-29	R*n	6.66
2024-11-02	*球	6
2024-11-021	*鑫	6.66
2024-11-25	*沙	5
2024-11-29	C*n	2.88