【土壤分类】基于matlab GUI多类SVM土壤分类【含Matlab源码 1398期】

海神之光

已于 2023-10-28 15:49:04 修改

阅读量996

点赞数 2

分类专栏： Matlab图像处理（进阶版）文章标签： matlab

于 2021-10-16 16:32:54 首次发布

代码事宜私信博主

本文链接：https://blog.csdn.net/TIQCmatlab/article/details/120800517

版权

Matlab图像处理（进阶版）专栏收录该内容

797 篇文章 958 订阅

订阅专栏

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

⛄一、SVM简介

支持向量机(Support Vector Machine)是Cortes和Vapnik于1995年首先提出的，它在解决小样本、非线性及高维模式识别中表现出许多特有的优势，并能够推广应用到函数拟合等其他机器学习问题中。
1 数学部分
1.1 二维空间
在这里插入图片描述

2 算法部分

⛄二、部分源代码

% Project Title: Soil Detection & Classification

function varargout = SoilDetect_GUI(varargin)
% SOILDETECT_GUI MATLAB code for SoilDetect_GUI.fig
% SOILDETECT_GUI, by itself, creates a new SOILDETECT_GUI or raises the existing
% singleton*.
%
% H = SOILDETECT_GUI returns the handle to a new SOILDETECT_GUI or the handle to
% the existing singleton*.
%
% SOILDETECT_GUI(‘CALLBACK’,hObject,eventData,handles,…) calls the local
% function named CALLBACK in SOILDETECT_GUI.M with the given input arguments.
%
% SOILDETECT_GUI(‘Property’,‘Value’,…) creates a new SOILDETECT_GUI or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before SoilDetect_GUI_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to SoilDetect_GUI_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE’s Tools menu. Choose “GUI allows only one
% instance to run (singleton)”.
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help SoilDetect_GUI

% Last Modified by GUIDE v2.5 28-Aug-2021 14:31:16

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct(‘gui_Name’, mfilename, …
‘gui_Singleton’, gui_Singleton, …
‘gui_OpeningFcn’, @SoilDetect_GUI_OpeningFcn, …
‘gui_OutputFcn’, @SoilDetect_GUI_OutputFcn, …
‘gui_LayoutFcn’, [] , …
‘gui_Callback’, []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT

% — Executes just before SoilDetect_GUI is made visible.
function SoilDetect_GUI_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to SoilDetect_GUI (see VARARGIN)

% Choose default command line output for SoilDetect_GUI
handles.output = hObject;
handles.output = hObject;
ss = ones(300,400);
axes(handles.axes1);
imshow(ss);
axes(handles.axes2);
imshow(ss);
% Update handles structure
guidata(hObject, handles);

% UIWAIT makes SoilDetect_GUI wait for user response (see UIRESUME)
% uiwait(handles.figure1);

% — Outputs from this function are returned to the command line.
function varargout = SoilDetect_GUI_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;

% — Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
clc
[filename, pathname] = uigetfile({‘.’;‘.bmp’;'.jpg’;‘*.gif’}, ‘Pick a Soil Image’);
I = imread([pathname,filename]);
I = imresize(I,[256,256]);
I2 = imresize(I,[300,400]);
axes(handles.axes1);
imshow(I2);title(‘Query Image’);
ss = ones(300,400);
axes(handles.axes2);
imshow(ss);
handles.ImgData1 = I;
guidata(hObject,handles);

% — Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
I3 = handles.ImgData1;
I4 = imadjust(I3,stretchlim(I3));
I5 = imresize(I4,[300,400]);
axes(handles.axes2);
imshow(I5);title(’ Contrast Enhanced ');
handles.ImgData2 = I4;
%% Feature Extraction Part
% queryimage = I4;
%queryImage = imresize(queryImage, [256 256]);
hsvHist = hsvHistogram(I4);
autoCorrelogram = colorAutoCorrelogram(I4);
color_moments = colorMoments(I4);
% for gabor filters we need gary scale image
img = double(rgb2gray(I4))/255;
[meanAmplitude, msEnergy] = gaborWavelet(img, 4, 6); % 4 = number of scales, 6 = number of orientations
wavelet_moments = waveletTransform(I4);
% construct the queryImage feature vector
Feature_Vector = [hsvHist autoCorrelogram color_moments meanAmplitude msEnergy wavelet_moments];
whos Feature_Vector
F1 = mean2(hsvHist(😃);
F2 = mean2(autoCorrelogram(😃);
F3 = mean2(color_moments(😃);
F4 = mean2(meanAmplitude(😃);
F5 = mean2(msEnergy(😃);
F6 = mean2(wavelet_moments(😃);

set(handles.edit3,'string',F1);
set(handles.edit4,'string',F2);
set(handles.edit5,'string',F3);
set(handles.edit6,'string',F4);
set(handles.edit7,'string',F5);
set(handles.edit8,'string',F6);

handles.ImgData3 = Feature_Vector;

guidata(hObject,handles);
%%
% — Executes on button press in pushbutton3.
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
load(‘TrainFeat_Soil.mat’)

test = handles.ImgData3;
result = multisvm(TrainFeat,Train_Label,test);
disp(result);

if result == 1
A1 = ‘Clay’;
set(handles.edit1,‘string’,A1);
helpdlg(’ Clay ‘);
disp(’ Clay ‘);
elseif result == 2
A2 = ‘Clayey Peat’;
set(handles.edit1,‘string’,A2);
helpdlg(’ Clayey Peat ‘);
disp(‘Clayey Peat’);
elseif result == 3
A3 = ‘Clayey Sand’;
set(handles.edit1,‘string’,A3);
helpdlg(’ Clayey Sand ‘);
disp(’ Clayey Sand ‘);
elseif result == 4
A4 = ‘Humus Clay’;
set(handles.edit1,‘string’,A4);
helpdlg(’ Humus Clay ‘);
disp(’ Humus Clay ‘);
elseif result == 5
A5 = ‘Peat’;
set(handles.edit1,‘string’,A5);
helpdlg(’ Peat ‘);
disp(’ Peat ‘);
elseif result == 6
A6 = ‘Sandy Clay’;
set(handles.edit1,‘string’,A6);
helpdlg(’ Sandy Clay ‘);
disp(‘Sandy Clay’);
elseif result == 7
A7 = ‘Silty Sand’;
set(handles.edit1,‘string’,A7);
helpdlg(’ Silty Sand ‘);
disp(’ Silty Sand ');
end
guidata(hObject,handles);
% — Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton4 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
load(‘Accuracy_Data.mat’)
Accuracy_Percent= zeros(200,1);
itr = 500;
hWaitBar = waitbar(0,‘Evaluating Maximum Accuracy with 500 iterations’);
for i = 1:itr
data = Train_Feat;
%groups = ismember(Train_Label,1);
groups = ismember(Train_Label,0);
[train,test] = crossvalind(‘HoldOut’,groups);
cp = classperf(groups);
svmStruct = svmtrain(data(train,:),groups(train),‘showplot’,false,‘kernel_function’,‘linear’);
classes = svmclassify(svmStruct,data(test,:),‘showplot’,false);
classperf(cp,classes,test);
Accuracy = cp.CorrectRate;
Accuracy_Percent(i) = Accuracy.*100;
sprintf(‘Accuracy of Linear Kernel is: %g%%’,Accuracy_Percent(i))
waitbar(i/itr);
end

⛄三、运行结果

在这里插入图片描述

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]胡晓艳,宋海燕.基于支持向量机和近红外光谱特性的土壤质地分类[J].山西农业科学. 2017,45(10)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流人群疏散病毒扩散晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合