CN104200230B - 一种基于小波变换和svm分类的土壤土质识别方法 - Google Patents

一种基于小波变换和svm分类的土壤土质识别方法 Download PDF

Info

Publication number
CN104200230B
CN104200230B CN201410460934.3A CN201410460934A CN104200230B CN 104200230 B CN104200230 B CN 104200230B CN 201410460934 A CN201410460934 A CN 201410460934A CN 104200230 B CN104200230 B CN 104200230B
Authority
CN
China
Prior art keywords
mrow
soil
subgraph
layering
wavelet transformation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410460934.3A
Other languages
English (en)
Other versions
CN104200230A (zh
Inventor
庄树峰
尹振东
吴芝路
尹亮
蒋旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201410460934.3A priority Critical patent/CN104200230B/zh
Publication of CN104200230A publication Critical patent/CN104200230A/zh
Application granted granted Critical
Publication of CN104200230B publication Critical patent/CN104200230B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了一种基于小波变换和SVM分类的土壤土质识别方法。本发明通过对土壤可见光图像进行小波变换得到分层子图像，然后，计算分层子图像的能量值作为特性向量，最后通过支持向量机(SVM)分类器进行土壤可见光图像土质识别。该方法有效的降低了设备复杂度和设备成本，提高了土壤土质识别的实时性。

Description

一种基于小波变换和SVM分类的土壤土质识别方法
技术领域
本发明涉及图像处理和模式识别领域，具体涉及一种基于小波变换和SVM分类的土壤识别方法。
背景技术
土壤土质是土壤重要的物理特性之一，与土壤通气性，肥力、保水和耕作的难易程度有着密不可分的关系。土壤土质一般分为砂土、壤土和粘土三类。砂土保肥性弱，缺少粘粒和有机质。粘土养分丰富，有较高的有机质含量。壤土介于砂土和粘土之间，是较理想的土壤，适合耕作。因此，通过对土壤土质进行识别，能够对耕作和施肥提供有效的依据，为精准农业提供有力的技术支撑。土壤土质的确定，如果采取测定土壤粒级百分数或物理粘粒含量的方法，则成本高昂，设备复杂，实时性差。到目前为止，常采用观察与简易的测定方法，包括手指测定法，松紧度测定法，观察法。这些方法只能定性地分析土壤土质，误差大，而且需要对三种土质的物理特性有较深的辨识经验，需要人工对土壤进行操作。为此，提供一种成本低，实时性强，识别精度高的土壤土质识别方法尤为重要。
发明内容
本发明的目的在于克服现有土壤土质识别方法成本高昂，设备复杂，实时性差的不足，提供一种基于小波变换和SVM分类的土壤土质识别方法。
为了解决背景技术锁存在的问题，本发明采用以下技术方案：
一种基于小波变换和SVM分类的土壤土质识别方法，其特征在于，它包括以下步骤：
步骤一：对每幅包含512*512个像素的土壤可见光图像，取40个256*256个像素的子图像，其中20个子图像作为训练样本，剩下的20个图像作为测试样本；
步骤二：对每幅子图像做2层正交小波变换，就可以得到2*3+1＝7个分层子图像；
步骤三：对每个分层子图像做能量计算，取每个分层子图像的能量值作为该分层子图像的特征值，从而每幅子图像可以得到一个7维的特征向量；
步骤四：基于砂土、壤土、粘土三类土壤土质，对支持向量机SVM进行训练，输入训练样本，构造3个2类SVM分类器，每个类对应其中的一个，使得第k个2类分类器所构造的分类超平面能够把第k类与其他的k-1类分割开；
步骤五：输入测试样本，对测试数据计算对应于各个分类器的决策函数值并选择最大值所对应的类别作为测试样本所属的类别。
进一步的，在步骤三中，对每个经过小波变换得到的分层子图像作能量计算，用来作为特征值，取信息熵作为能量：
P(i,j)由下式可得
其中C(i,j)表示图像在(i,j)处的灰度值。
进一步的，在步骤四中，将步骤一得到的20个训练样本子图像输入支持向量机SVM，根据如下公式求解最优分类超平面：
式中，x为训练样本；i，j＝1,2，…，n，n为训练样本的个数，T表示矩阵求转置，ai是函数优化时的拉格朗日系数，即为支持向量，由此得到的判别函数为：
其中sgn(c)表示返回c的正负号，为非零的支持向量，b*为分类阈值。
进一步的，在步骤五中，将步骤一得到的20个测试样本子图像输入支持向量机SVM，根据步骤四中的判别函数计算出测试样本在各个类别中对应的判别函数值，选择其中的最大值所对应的类别作为该测试样本所述的类别，直到所有的20个测试样本都分类完成。
本发明对比现有技术，有如下的有益效果：本发明提供了一种基于小波变换特征提取和SVM分类算法的土壤土质识别技术，本技术结合图像处理和模式识别，通过对土壤可见光图像进行小波变换提取能量特征，组成特征向量，采用支持向量机训练分类器进行土质识别。该方法有效的降低了设备复杂度和设备成本，提高了土壤土质识别的实时性。
附图说明
图1是本发明的方法流程图；
图2是2层小波分解过程示意图；
图3是基于SVM的土壤土质识别过程图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步描述：
图1是本发明的方法流程图。一种基于小波变换特征提取和SVM分类算法的土壤土质识别方法，它包括以下步骤：
步骤一：对每幅包含512*512个像素的土壤可见光图像，取40个256*256个像素的子图像，其中20个子图像作为训练样本，剩下的20个图像作为测试样本；
步骤二：对每幅子图像做2层正交小波变换，就可以得到2*3+1＝7个分层子图像；
步骤三：对每个分层子图像做能量计算，取每个分层子图像的能量值作为该分层子图像的特征值，从而每幅子图像可以得到一个7维的特征向量；
步骤四：基于砂土、壤土、粘土三类土壤土质，对支持向量机SVM进行训练，输入训练样本，构造3个2类SVM分类器，每个类对应其中的一个，使得第k个2类分类器所构造的分类超平面能够把第k类与其他的k-1类分割开；
步骤五：输入测试样本，对测试数据计算对应于各个分类器的决策函数值并选择最大值所对应的类别作为测试样本所属的类别。
进一步的，在步骤三中，对每个对每个经过小波变换得到的分层子图像作能量计算，用来作为特征值，取信息熵作为能量：
P(i,j)由下式可得
其中C(i,j)表示图像在(i,j)处的灰度值。
进一步的，在步骤四中，将步骤一得到的20个训练样本子图像输入支持向量机SVM，根据如下公式求解最优分类超平面：
式中，x为训练样本；i，j＝1,2，…，n，n为训练样本的个数，T表示矩阵求转置，ai是函数优化时的拉格朗日系数，即为支持向量，由此得到的判别函数为：
其中sgn(c)表示返回c的正负号，为非零的支持向量，b*为分类阈值。
进一步的，在步骤五中，将步骤一得到的20个测试样本子图像输入支持向量机SVM，根据步骤四中的判别函数计算出测试样本在各个类别中对应的判别函数值，选择其中的最大值所对应的类别作为该测试样本所述的类别，直到所有的20个测试样本都分类完成。
图2是2层小波分解过程示意图。说明了一幅原图像进行2层正交小波变换的过程，将图像看成二维的离散函数，在小波基中分别做投影，得到各个方向的小波系数，也就是分尺度、分方向的分层子图像。每幅原图像通过2层正交小波变换就可以分解为7个分层子图像，其中分层子图像1为原始图像的低分辨率图像，分层子图像3和5为水平方向的高频成分，分层子图像2和6为垂直方向的高频成分，分层子图像4和7为对角方向的高频成分。
图3是基于SVM的土壤土质识别过程，包括训练阶段和测试阶段，训练阶段通过输入训练样本对SVM进行训练。测试阶段即为在训练阶段结束后输入测试样本到SVM，对测试样本进行分类。图3中特征提取对应于本发明方法步骤二和步骤三，训练阶段对应于步骤四，测试阶段对应于步骤五。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于小波变换和SVM分类的土壤土质识别方法，其特征在于，它包括以下步骤：
步骤一：对每幅包含512*512个像素的土壤可见光图像，取40个256*256个像素的子图像，其中20个子图像作为训练样本，剩下的20个图像作为测试样本；
步骤二：对每幅子图像做2层正交小波变换，就可以得到2*3+1＝7个分层子图像；
步骤三：对每个分层子图像做能量计算，取每个分层子图像的能量值作为该分层子图像的特征值，从而每幅子图像可以得到一个7维的特征向量；
具体为：对每个经过小波变换得到的分层子图像作能量计算，用来作为特征值，取信息熵作为能量：
<mrow> <mi>e</mi> <mo>=</mo> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </munder> <mo>|</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>|</mo> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mo>|</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow>
P(i,j)由下式可得
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>&amp;lsqb;</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </munder> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>/</mo> <mn>4</mn> </mrow> </mfrac> </mrow>
其中C(i,j)表示图像在(i,j)处的灰度值；
步骤四：基于砂土、壤土、粘土三类土壤土质，对支持向量机SVM进行训练，输入训练样本，构造3个2类SVM分类器，每个类对应其中的一个，使得第k个2类分类器所构造的分类超平面能够把第k类与其他的k-1类分割开；
步骤五：输入测试样本，对测试数据计算对应于各个分类器的决策函数值并选择最大值所对应的类别作为测试样本所属的类别。
2.根据权利要求1所述的一种基于小波变换和SVM分类的土壤土质识别方法，其特征在于，在步骤五中，将步骤一得到的20个测试样本子图像输入支持向量机SVM，根据步骤四中的判别函数计算出测试样本在各个类别中对应的判别函数值，选择其中的最大值所对应的类别作为该测试样本所述的类别，直到所有的20个测试样本都分类完成。
CN201410460934.3A 2014-09-11 2014-09-11 一种基于小波变换和svm分类的土壤土质识别方法 Active CN104200230B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410460934.3A CN104200230B (zh) 2014-09-11 2014-09-11 一种基于小波变换和svm分类的土壤土质识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410460934.3A CN104200230B (zh) 2014-09-11 2014-09-11 一种基于小波变换和svm分类的土壤土质识别方法

Publications (2)

Publication Number Publication Date
CN104200230A CN104200230A (zh) 2014-12-10
CN104200230B true CN104200230B (zh) 2018-04-27

Family

ID=52085519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410460934.3A Active CN104200230B (zh) 2014-09-11 2014-09-11 一种基于小波变换和svm分类的土壤土质识别方法

Country Status (1)

Country Link
CN (1) CN104200230B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105181933B (zh) * 2015-09-11 2017-04-05 北华航天工业学院 预测土壤压缩系数的方法
SG10201902144RA (en) * 2019-03-11 2020-10-29 Housing And Dev Board Apparatus, system and method for classification of soil and soil types

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101488188A (zh) * 2008-11-10 2009-07-22 西安电子科技大学 基于混合核函数的svm分类器的sar图像分类方法
CN101630369A (zh) * 2009-07-30 2010-01-20 南京航空航天大学 一种基于小波分形特征的行人检测方法
CN102063627A (zh) * 2010-12-31 2011-05-18 宁波大学 基于多小波变换的自然图像和计算机生成图像的识别方法
CN102521591A (zh) * 2011-11-29 2012-06-27 北京航空航天大学 一种复杂背景下小目标快速识别方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101488188A (zh) * 2008-11-10 2009-07-22 西安电子科技大学 基于混合核函数的svm分类器的sar图像分类方法
CN101630369A (zh) * 2009-07-30 2010-01-20 南京航空航天大学 一种基于小波分形特征的行人检测方法
CN102063627A (zh) * 2010-12-31 2011-05-18 宁波大学 基于多小波变换的自然图像和计算机生成图像的识别方法
CN102521591A (zh) * 2011-11-29 2012-06-27 北京航空航天大学 一种复杂背景下小目标快速识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
小波分析理论与影像分析;朱长青;《测绘出版社》;19980331;第90-92页 *
数字图像模式识别工程项目研究;蒋先刚;《西南交通大学出版社》;20140331;第303-304页 *
水稻卫星遥感不确定性研究;黄敬峰等;《浙江大学出版社》;20130831;第36-37页 *

Also Published As

Publication number Publication date
CN104200230A (zh) 2014-12-10

Similar Documents

Publication Publication Date Title
Zhuang et al. Strategies combining spectral angle mapper and change vector analysis to unsupervised change detection in multispectral images
CN103984959B (zh) 一种基于数据与任务驱动的图像分类方法
CN104021394B (zh) 基于AdaBoost算法的绝缘子图像识别方法
CN103810704B (zh) 基于支持向量机和判别随机场的sar图像变化检测方法
CN105931253A (zh) 一种基于半监督学习相结合的图像分割方法
CN107229917A (zh) 一种基于迭代聚类的多幅遥感影像共性显著目标检测方法
CN103745233B (zh) 基于空间信息迁移的高光谱图像分类方法
CN102867183B (zh) 一种车辆遗撒物检测方法、装置及智能交通监控系统
CN106778687A (zh) 基于局部评估和全局优化的注视点检测方法
CN105261004A (zh) 基于均值漂移和邻域信息的模糊c均值图像分割方法
CN105930846A (zh) 基于邻域信息和svgdl的极化sar图像分类方法
CN104282008B (zh) 对图像进行纹理分割的方法和装置
CN105335975B (zh) 基于低秩分解和直方图统计的极化sar图像分割方法
CN105913073B (zh) 基于深度增量支持向量机的sar图像目标识别方法
CN104102928B (zh) 一种基于纹理基元的遥感图像分类方法
CN106960415A (zh) 一种基于像素递归超分辨率模型恢复图像的方法
CN104732552B (zh) 基于非平稳条件场的sar图像分割方法
CN106127221A (zh) 基于极化‑纹理特征与dpl的极化sar图像分类方法
CN105069774A (zh) 基于多示例学习与图割优化的目标分割方法
CN105046268A (zh) 基于Wishart深度网络的极化SAR图像分类方法
CN102982544A (zh) 多前景目标图像交互式分割方法
CN104200230B (zh) 一种基于小波变换和svm分类的土壤土质识别方法
CN104268557B (zh) 基于协同训练和深度svm的极化sar分类方法
CN104463130B (zh) 一种基于评估反馈系统的车牌图像光照处理方法
CN108090913A (zh) 一种基于对象级Gauss-Markov随机场的图像语义分割方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant