一种图像纹理特征提取分类方法

1.本发明涉及图像纹理提取方法领域，具体为一种图像纹理特征提取分类方法。


背景技术：

2.纹理特征提取是图像处理、图像分析、计算机视觉等诸多领域的基础研究之一。纹理特征的类内差异性及类间相似性决定了仅依靠单一的提取算法无法准确真实的表现图像的纹理特征，因此，通过多种纹理特征提取算法的融合来提高提取的效率和准确率是当前研究的热点。这些年随着研究的不断深入，图像纹理的分类在诸多领域都得到了很重要的应用，例如在虚拟现实触觉再现领域，通过再现设备可以识别不同的材质，得到更加真实的触觉体验，增加使用者的沉浸感和真实感。
3.现有的纹理特征提取方法采用单一的纹理提取算法，此举不能准确的描述图像纹理的问题，而且冗余信息较多，现有的纹理特征提取方法的计算效率低、鲁棒性差、提取精度低，较难对错误的图像进行有效的区分，为此，提出一种图像纹理特征提取分类方法。


技术实现要素：

4.本发明的目的在于提供一种图像纹理特征提取分类方法，以解决上述背景技术中提出的现有的纹理特征提取方法采用单一的纹理提取算法，此举不能准确的描述图像纹理的问题，而且冗余信息较多，现有的纹理特征提取方法的计算效率低、鲁棒性差、提取精度低，较难对错误的图像进行有效的区分的问题。
5.为实现上述目的，本发明提供如下技术方案：一种图像纹理特征提取分类方法，包括以下步骤：
6.s1、对输入的图像进行gamma校正，将部分灰度区域映射到更宽或更窄的区域，对输入图像的灰度值进行指数变换；
7.s2、将校正后的图像进行自适应对比度增强算法(ace)处理，采用反锐化掩膜技术，把图像分为低频的反锐化掩膜以及高频部分，用原图减去反锐化掩膜获取高频部分，最后将高频部分加入反锐化掩膜并加上增益系数g(i,j)，得到增强的图像；
8.s3、使用基于t
‑
glcm和tamura融合的特征提取方法，对经过gamma校正和自适应对比度增强的图像使用基于t
‑
glcm和tamura融合的特征提取方法对图像纹理特征进行提取并生成特征矩阵；
9.s4、采用支持向量机(svm)对纹理材质进行分类，通过带高斯径向基函数的支持向量机对不同材质的纹理图像进行分类。
10.优选的，在s1中，gamma校正的表达式为:
[0011][0012]
其中，a为常数，v
in
为原图像，v
out
为校正后的图像，γ为校正值。
[0013]
优选的，在s2中，设x(i,j)是原像素，通过自适应对比度增强后的像素值为：
[0014]
f(i,j)＝m
x
(i,j)+g(i,j)[x(i,j)
‑
m
x
(i,j)]
[0015]
获取图像的反锐化掩膜即低频部分通过求取图像中心各像素点为中心的局部区域的像素平均值求得，中心像素点的低频部分的像素值通过下列公式计算得到：
[0016][0017]
其中，(i,j)、(k,l)分别为像素点和像素点周围邻域的坐标；
[0018]
局部增益g(i,j)的值是恒大于1的，所以[x(i,j)
‑
mx(i,j)]就会变大，局部增益的计算公式如下：
[0019][0020][0021]
其中，d为一个常数，d取值为图像的全局均方差；σ
x
(i,j)为局部均方差。
[0022]
优选的，在s3中，对图像纹理特征进行提取并生成特征矩阵的此阶段分为如下两个步骤：
[0023]
s3.1、基于改进的灰度共生矩阵(t
‑
glcm)特征提取：
[0024]
首先对传统的灰度共生矩阵(glcm)进行转置操作后与glcm取平均值得到glcm1；
[0025]
然后对glcm1从0
°
、45
°
、90
°
和135
°
四个方向上得到的各个矩阵进行相加取平均值得到最终的共生矩阵t
‑
glcm；采用的四种纹理特征统计量分别为能量(asm)、熵(ent)、对比度(con)和逆差分矩(idm)；
[0026]
s3.2、基于tamura特征提取，通过计算对比度、方向度和粗糙度三个特征对图像纹理进行分析。
[0027]
优选的，在s3.1中，所述能量(asm)、熵(ent)、对比度(con)和逆差分矩(idm)，四个特征统计量的计算方法公式分别为：
[0028]
能量(asm)：
[0029][0030]
熵(ent)：
[0031][0032]
对比度(con)：
[0033][0034]
逆差分矩(idm)：
[0035][0036]
优选的，在s3.2中，对比度是图像中最亮与最暗两部分的灰度之差，对比度的计算步骤如下：
[0037]
a、像素(x,y)的3*3邻域像素的平均值为：
[0038][0039]
其中，f(i,j)为像素点，(i,j)为灰度值，
[0040]
b、标准差为：
[0041][0042]
c、四阶矩为：
[0043][0044]
d、对比度为：
[0045][0046]
优选的，在s3.2中，方向度是图像中的纹理沿着某个方向呈现出某种形态，具有一定的方向性，通过计算梯度向量来计算出方向度，方向度的计算步骤如下：
[0047]
a、像素点(x,y)的梯度向量的模以及方向如下：
[0048]
|δg(x,y)|＝(|δ
h
|+|δ
v
|)/2
[0049][0050]
其中，δ
h
和δ
v
分别表示对图像进行卷积后的结果，矩阵分别为：
[0051]
和
[0052]
b、设定阈值t＝12，像素点(x,y)的方向角d(x,y)计算公式如下：
[0053][0054]
c、像素点(x,y)的3*3邻域所有像素的方向角均值：
[0055][0056]
d、像素点(x,y)的方向角为：
[0057]
d'(x,y)＝|d(x,y)
‑
μ(x,y)|。
[0058]
优选的，在s3.2中，粗糙度的计算方法是用线性量化代替指数量化又可以降低算法的复杂度，粗糙度的计算步骤如下：
[0059]
a、以像素(i,j)为中心，水平方向平均灰度方差为：
[0060]
e
n,u
(i,j)＝|a
n
(i
‑
n,j)
‑
a
n
(i+n,j)|
[0061]
b、垂直方向平均灰度方差为：
[0062]
e
n,v
(i,j)＝|a
n
(i,j
‑
n)
‑
a
n
(i,j+n)|
[0063]
其中，a(i,j)为平均灰度值，n为1,2,3,4,5，
[0064]
c、不分方向得到最大值：
[0065]
s
best
(i,j)＝max{e
n,m
(i,j)|n∈[1,5],m＝u,v}
[0066]
d、以像素(i,j)为中心的3*3邻域平均粗糙度为：
[0067][0068]
e、像素粗糙度与平均粗糙度差的绝对值就是该像素的粗糙度度量。
[0069]
优选的，所述glcm1计算公式为
[0070][0071]
式中得到的glcm1具有旋转不变性。
[0072]
优选的，所述高斯径向基函数的公式为:
[0073][0074]
其中，x'是内核函数中心，σ是函数的宽度参数，控制函数的径向范围。
[0075]
与现有技术相比，本发明的有益效果是：本发明提出基于t
‑
glcm和tamura融合特征的纹理材质分类方法，对传统灰度共生矩阵(glcm)进行改进，提出了具有改善旋转不变性和减少冗余信息的改进的灰度共生矩阵t
‑
glcm算法，通过gamma校正与自适应对比度增强(ace)处理再建立样本集与测试集，使用支持向量机(svm)进行分类，由实验结果可以得出，融合的特征提取方法较单一的纹理提取方法更加高效、分类精度更高且鲁棒性也更好，而且可以对错误的图像进行有效的区分。
附图说明
[0076]
图1为本发明的算法流程图。
具体实施方式
[0077]
下面将结合本发明实施例中的内容，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0078]
实施例
[0079]
本发明提供一种技术方案：一种图像纹理特征提取分类方法，包括以下步骤：
[0080]
s1、对输入的图像进行gamma校正，将部分灰度区域映射到更宽或更窄的区域，以达到增强图像的效果，目的是校正图像的亮度偏差,从而减少因图像亮度偏差对纹理提取精度的影响，方法是对输入图像的灰度值进行指数变换；gamma校正的表达式为:
[0081][0082]
其中，a为常数，v
in
为原图像，v
out
为校正后的图像，γ为校正值；
[0083]
s2、将校正后的图像进行自适应对比度增强算法(ace)处理，采用反锐化掩膜技
术，把图像分为低频的反锐化掩膜以及高频部分，用原图减去反锐化掩膜获取高频部分，最后将高频部分加入反锐化掩膜并加上增益系数g(i,j)，得到增强的图像，设x(i,j)是原像素，通过自适应对比度增强后的像素值为：
[0084]
f(i,j)＝m
x
(i,j)+g(i,j)[x(i,j)
‑
m
x
(i,j)]
[0085]
获取图像的反锐化掩膜即低频部分通过求取图像中心各像素点为中心的局部区域的像素平均值求得，中心像素点的低频部分的像素值可以通过下列公式计算得到：
[0086][0087]
其中，(i,j)、(k,l)分别为像素点和像素点周围邻域的坐标，局部增益g(i,j)的值是恒大于1的，所以[x(i,j)
‑
mx(i,j)]就会变大，局部增益的计算公式如下：
[0088][0089][0090]
其中，d为一个常数，它的值可以选择全局平均值或者是全局均方差，本文的d取值为图像的全局均方差；σ
x
(i,j)为局部均方差；
[0091]
s3、使用基于t
‑
glcm和tamura融合的特征提取方法，为了提高纹理特征分类的精准度，对经过gamma校正和自适应对比度增强的图像使用基于t
‑
glcm和tamura融合的特征提取方法对图像纹理特征进行提取并生成特征矩阵；本阶段具体分为两个步骤；
[0092]
步骤一：基于改进的灰度共生矩阵(t
‑
glcm)特征提取；直接采用传统的灰度共生矩阵(glcm)存在一个问题，其特征值随旋转角度的改变而变化且存在大量冗余信息；为了解决传统glcm存在的不足，提高特征提取的精度，本文提出了一种改进的灰度共生矩阵(t
‑
glcm)；改进的灰度共生矩阵的具体方法如下：首先对传统的灰度共生矩阵(glcm)进行转置操作后与glcm取平均值得到glcm1；
[0093]
然后对glcm1从0
°
、45
°
、90
°
和135
°
四个方向上得到的各个矩阵进行相加取平均值，得到最终的共生矩阵t
‑
glcm；
[0094][0095]
式中得到的glcm1具有旋转不变性；进一步得到共生矩阵t
‑
glcm，常用的纹理特征统计量有14个，经常采用的四种特征统计量分别为能量(asm)、熵(ent)、对比度(con)和逆差分矩(idm)，四个特征统计量的计算方法公式分别为：
[0096]
能量asm：
[0097][0098]
熵ent：
[0099][0100]
对比度con：
[0101][0102]
逆差分矩idm：
[0103][0104]
步骤二是基于tamura特征提取，方法为对该图像进行基于tamura特征提取，本发明应用对比度、方向度和粗糙度三个特征对图像纹理进行分析，图像中最亮与最暗两部分的灰度之差就是对比度；对比度越大，图像更清晰，反之，对比度越小，效果越模糊；对比度的计算方法如下：
[0105]
a、像素(x,y)的3*3邻域像素的平均值为：
[0106][0107]
其中，f(i,j)为像素点，(i,j)为灰度值，
[0108]
b、标准差为：
[0109][0110]
c、四阶矩为：
[0111][0112]
d、对比度为：
[0113][0114]
方向度是图像中的纹理沿着某个方向呈现出某种形态，具有一定的方向性，其可以通过计算梯度向量来计算出方向度；计算步骤如下：
[0115]
a、像素点(x,y)的梯度向量的模以及方向如下：
[0116]
|δg(x,y)|＝(|δ
h
|+|δ
v
|)/2
[0117][0118]
其中，δ
h
和δ
v
分别表示对图像进行卷积后的结果，矩阵分别为：
[0119]
和
[0120]
b、设定阈值t＝12，像素点(x,y)的方向角d(x,y)计算公式如下：
[0121][0122]
c、像素点(x,y)的3*3邻域所有像素的方向角均值：
[0123][0124]
d、像素点(x,y)的方向角为：
[0125]
d'(x,y)＝|d(x,y)
‑
μ(x,y)|；
[0126]
粗糙度的计算方法是用线性量化代替指数量化又可以降低算法的复杂度，粗糙度的计算步骤如下：
[0127]
a、以像素(i,j)为中心，水平方向平均灰度方差为：
[0128]
e
n,u
(i,j)＝|a
n
(i
‑
n,j)
‑
a
n
(i+n,j)|
[0129]
b、垂直方向平均灰度方差为：
[0130]
e
n,v
(i,j)＝|a
n
(i,j
‑
n)
‑
a
n
(i,j+n)|
[0131]
其中，a(i,j)为平均灰度值，n为1,2,3,4,5，
[0132]
c、不分方向得到最大值：
[0133]
s
best
(i,j)＝max{e
n,m
(i,j)|n∈[1,5],m＝u,v}
[0134]
d、以像素(i,j)为中心的3*3邻域平均粗糙度为：
[0135][0136]
e、像素粗糙度与平均粗糙度差的绝对值就是该像素的粗糙度度量；
[0137]
s4、利用支持向量机(svm)对不同材质的纹理信息进行分类识别，支持向量(svm)是一种在分类和回归的情况下进行预测的技术，在模式识别、分类和非线性回归方面得到了广泛的应用；支持向量机是一种与学习算法相关的监督模型，对于给定的一组训练集，每个不同的示例被标记为两个类别中的一个，svm将新的示例分配给其中一个类别，使得其成为非概率二进制线性分类器；依赖于被称为结构风险最小化(srm)的机制，利用支持向量机将输入向量映射到高维空间，计算两组样本之间的最佳分离超平面；svm需要一组参数，包括支持向量机的类型、核函数的类型、罚因子和训练停止的标准；其中，svm的核心是核函数；核函数的类型有很多，经过不断的研究表明，高斯径向基函数在纹理图像分类方面取得了较高的分类精度，并且高斯径向基函数计算复杂度低；其公式如下：
[0138][0139]
其中，x'是内核函数中心，σ是函数的宽度参数，控制函数的径向范围。
[0140]
所以本文选择带高斯径向基函数的支持向量机对不同材质的纹理图像进行分类，而且系统将采用网格搜索方法获得最优参数，以此提升训练精度，采用交叉验证方法来提高预测精度。
[0141]
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。