?數(shù)據(jù)的排列： 橫： 樣本個(gè)數(shù)? 縱：光譜數(shù)據(jù)（350：2500）
一 、數(shù)據(jù)導(dǎo)入Matlab
數(shù)據(jù)導(dǎo)入的兩種方法：
1） 1----2
2）3 （即在空白工作區(qū)內(nèi) 點(diǎn)擊??右鍵---新的空白文件夾---雙擊點(diǎn)開---直接將數(shù)據(jù)粘貼進(jìn)去）
注意：數(shù)據(jù)導(dǎo)入過程中不要導(dǎo)入編號(hào)，只導(dǎo)入需要的數(shù)據(jù),例如: 42個(gè)樣本數(shù)據(jù)，350-2500的光譜，導(dǎo)入后數(shù)據(jù)為2151*42
工作表的重新命名和電腦桌面對(duì)文件夾的命名一樣袒啼，可直接雙擊文件，進(jìn)行名稱的重新編輯
數(shù)據(jù)的分析
1 查找代碼所在位置
2? 相關(guān)性分析 （農(nóng)學(xué)數(shù)據(jù)與光譜數(shù)據(jù)的相關(guān)性）
1）定義函數(shù)function ,該函數(shù)具體用法可百度查閱
2 代碼的運(yùn)行
代碼運(yùn)行時(shí)，第一次建議不要使用快捷鍵Enter , 可使用 編輯器---運(yùn)行（防止代碼文件位置不統(tǒng)一造成無法運(yùn)行或運(yùn)行出錯(cuò)）
3 運(yùn)行結(jié)果分析
左邊?------當(dāng)前文件夾------處理后文件直接帶點(diǎn)擊開拱礁，可查看生理指標(biāo)與反射率的相關(guān)性，第一例為原始光譜辕漂，第二列為一階導(dǎo)數(shù)
右邊---工作區(qū)---cor代表相關(guān)系數(shù)呢灶，若不明白，可以根據(jù)自己個(gè)人數(shù)據(jù)雙擊點(diǎn)開查看钉嘹。
相關(guān)性圖片的查看
代碼運(yùn)行后會(huì)出現(xiàn)一張?jiān)紨?shù)據(jù)和光譜數(shù)據(jù)的 相關(guān)性圖鸯乃，可根據(jù)圖查看數(shù)據(jù)相關(guān)性強(qiáng)弱
光譜指數(shù)的計(jì)算以及等勢(shì)圖的調(diào)整和查看
一、 以NDVI為例
和相關(guān)系數(shù)代碼運(yùn)行一樣跋涣，
第一步：粘貼第一行數(shù)據(jù)缨睡。
第二步:將數(shù)據(jù)粘貼在命令窗口?
?第三步：加; 代碼的運(yùn)行
代碼具體解析如下：
function[ndvicor] = ndvicorr1(parameter,refle)#定義函數(shù)
ndvicor = zeros(2151,2151);# 光譜個(gè)數(shù)鸟悴，如果你為350：1000的波段，則需更改為zeros(651:651)
for i = 1:2151# 循環(huán)
ndvi = (refle - refle(:,i)) ./ (refle + refle(:,i));#NDVI 的計(jì)算公式宏蛉，你可根據(jù)自己需求更改為其他植被指數(shù)的算法
ndvicor(i,:) = corr(ndvi,parameter).^2;#相關(guān)系數(shù)
end
ndvicor((1350-349):(1400-349),:) = NaN;#刪除不必要的范圍
ndvicor(:,(1350-349):(1400-349)) = NaN;#刪除不必要的范圍
ndvicor((1800-349):(1950-349),:) = NaN;#刪除不必要的范圍
ndvicor(:,(1800-349):(1950-349)) = NaN;#刪除不必要的范圍
figure?
contourf(350:2500,350:2500,ndvicor)#畫出等試圖 （橫坐標(biāo)遣臼，縱坐標(biāo)，據(jù)NDVI與農(nóng)學(xué)參數(shù)的決定系數(shù)）
colormap(jet);#為等試圖顏色組合拾并，具體顏色組合可查閱百度也可自行進(jìn)行具體的操作揍堰，見后面
#輸出了決定系數(shù)最大的前20個(gè)波段組合。
ndvicor2 = ndvicor;
ndvicor_max = zeros(20,3);
for i = 1:20
max1 = max(max(ndvicor2));
[x,y] = find(ndvicor2 == max1);
ndvicor2(ndvicor2 == max1) = 0;
a1 = [x,y];
ndvicor_max(i,:) = [max1,(a1(:,1)'+ 349)];
end
ndvicor_max
等試圖的具體調(diào)整
可代碼（我不會(huì)嗅义，Matlab有手動(dòng)調(diào)節(jié)屏歹，可不許代碼）
手動(dòng)調(diào)整具體如下:
