CN104200230B - 一种基于小波变换和svm分类的土壤土质识别方法 - Google Patents
一种基于小波变换和svm分类的土壤土质识别方法 Download PDFInfo
- 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
- classification
- 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
Links
- 239000002689 soil Substances 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000009466 transformation Effects 0.000 title claims abstract description 9
- 238000012706 support-vector machine Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims description 29
- 238000012549 training Methods 0.000 claims description 17
- 239000004927 clay Substances 0.000 claims description 8
- -1 loam Substances 0.000 claims 1
- 239000004576 sand Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000007635 classification algorithm Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009313 farming Methods 0.000 description 2
- 230000035558 fertility Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 238000005273 aeration Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000004720 fertilization Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003971 tillage Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于小波变换和SVM分类的土壤土质识别方法。本发明通过对土壤可见光图像进行小波变换得到分层子图像，然后，计算分层子图像的能量值作为特性向量，最后通过支持向量机(SVM)分类器进行土壤可见光图像土质识别。该方法有效的降低了设备复杂度和设备成本，提高了土壤土质识别的实时性。
Description
技术领域
本发明涉及图像处理和模式识别领域，具体涉及一种基于小波变换和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>&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>&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>&rsqb;</mo>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</munder>
<msup>
<mrow>
<mo>&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>&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个测试样本都分类完成。
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)
| 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)
| 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 | 北京航空航天大学 | 一种复杂背景下小目标快速识别方法 |
-
2014
- 2014-09-11 CN CN201410460934.3A patent/CN104200230B/zh active Active
Patent Citations (4)
| 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)
| 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 |
|---|---|---|
| CN102609681B (zh) | 基于字典学习模型的人脸识别方法 | |
| CN105261000B (zh) | 一种基于端元提取与光谱解混的高光谱图像融合方法 | |
| CN102169545B (zh) | 一种高分辨率遥感图像的变化检测方法 | |
| CN101329736B (zh) | 基于特征选择和隐马尔科夫模型的图像分割方法 | |
| CN108830243A (zh) | 基于胶囊网络的高光谱图像分类方法 | |
| CN104732546B (zh) | 区域相似性和局部空间约束的非刚性sar图像配准方法 | |
| CN104282008B (zh) | 对图像进行纹理分割的方法和装置 | |
| CN102096819A (zh) | 利用稀疏表示与字典学习进行图像分割的方法 | |
| CN103366371B (zh) | 基于k分布和纹理特征的sar图像分割方法 | |
| CN103345760B (zh) | 一种医学图像对象形状模板标记点的自动生成方法 | |
| CN102074013B (zh) | 基于小波域多尺度Markov网模型的图像分割方法 | |
| CN103049760B (zh) | 基于图像分块和位置加权的稀疏表示目标识别方法 | |
| CN106960415A (zh) | 一种基于像素递归超分辨率模型恢复图像的方法 | |
| CN102122353A (zh) | 利用增量字典学习与稀疏表示进行图像分割的方法 | |
| CN105069774A (zh) | 基于多示例学习与图割优化的目标分割方法 | |
| CN102054170A (zh) | 基于极小化上界误差的视觉跟踪方法 | |
| CN108764330A (zh) | 基于超像素分割和卷积反卷积网络的sar图像分类方法 | |
| CN103745233A (zh) | 基于空间信息迁移的高光谱图像分类方法 | |
| CN108764357A (zh) | 基于压缩-激发的聚合残差网络高光谱图像分类方法 | |
| CN104021567B (zh) | 基于首数字定律的图像高斯模糊篡改检测方法 | |
| CN106097290A (zh) | 基于nmf图像融合的sar图像变化检测方法 | |
| CN104200230B (zh) | 一种基于小波变换和svm分类的土壤土质识别方法 | |
| CN104268557B (zh) | 基于协同训练和深度svm的极化sar分类方法 | |
| CN102722717B (zh) | 一种细胞分裂识别方法 | |
| CN104268555B (zh) | 基于模糊稀疏lssvm的极化sar图像分类方法 |
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 |