基于MATLAB的BP神经网络反演水质参数

点滴学习，贵在记录

基于MATLAB的BP神经网络反演水质参数

在做一个关于水质反演的项目时，常规模型精度较低，因此考虑采用bp神经网络来做反演。我要实现的是基于样本数据训练BP神经网络，然后基于此导入tif影像数据，然后得到水质参数的空间分布影像（带坐标系）。

读取样本点。将k个样本点分成训练集和测试集，我的样本点数据为k*4的矩阵，k为样本点个数，4为波段数。数据需要转置后并进行归一化处理才可导入net模型。建立神经网络采用newff()函数，训练后的数据性需要进行转成反演值。

% 神经网络建立
net=feedforwardnet(hiddenSizes,trainFcn)

hiddenSizes：神经元个数，默认单层10个。
trainFcn：训练函数，默认’trainlm’

%  数据归一化
[INPUT_train, ps_input] = mapminmax(input_train,0,1);

% 反归一化
T_sim = mapminmax('reverse',t_sim,ps_output);

模型训练完成便可做预测，对水质参数进行反演

首先读取影像，
影像为m×n×l的多光谱影像。需要转成m×n行，l 列的矩阵，然后转置成l行，m×n列的矩阵，代码实现如下

[A,R] = geotiffread('F:\DEMO\DEMO1.tif')
info = geotiffinfo('F:\DEMO\DEMO1.tif')
%转为m*n,l矩阵
a = double(reshape(A,[info.Height*info.Width,4]))
%影像转置
inputimg = a(:,:)';

归一化后做仿真测试，得到反演结果并输出tif格式的影像。
outimg = sim(net,Inputimg)

更新
新版本的MATLAB里面，前馈神经网络已经函数已经改为 feedforwardnet()

2022.1.17更新
注意！输出的结果需要reshape成输入影像的长宽。

outimg = reshape(outimg, [info.Height, info.Width])