高光谱

position = zeros(16,2);
position(1,1:2) = [887,718];
position(2,1:2) = [1175,718];
position(3,1:2) = [1419,702];
position(4,1:2) = [711,914];
position(5,1:2) = [915,910];
position(6,1:2) = [1107,906];
position(7,1:2) = [1279,906];
position(8,1:2) = [1435,902];

position(9,1:2) = [1095,1118];
position(10,1:2) = [1283,1122];
position(11,1:2) = [1475,1110];
position(12,1:2) = [915,1298];
position(13,1:2) = [1095,1298];
position(14,1:2) = [1267,1294];
position(15,1:2) = [1475,1274];
position(16,1:2) = [1211,1038];


x = 650:10:1100;
point = zeros(1,46);

%测试窗函数
pic = I{65};
figure,imagesc(pic),colormap gray,axis tight equal;
for k = 1:16
pic(position(k,2):position(k,2)+29,position(k,1):position(k,1)+29) = max(pic(:));
end
figure,imagesc(pic),colormap gray,axis tight equal;

figure;
for k = 1:16
for i = 65:110
pic = I{i};
point(i-64) = mean(mean(pic(position(k,2):position(k,2)+29 , position(k,1):position(k,1)+29)));
% point(i-64) = pic(position(k,2),position(k,1));
end
subplot(4,4,k);
plot(x,point);
hold on;
end

figure;
%str = {'材料1','材料2','材料3','材料4','材料5','材料6','材料7','材料8','材料9','材料10','材料11','材料12','材料13','材料14','材料15'};
str = {'材料1','材料2','材料3','材料4','材料5','材料6','材料7','材料8','材料9','聚酰亚胺渗碳膜（防杂光）','双面镀铝薄膜','低反射率涂层绸','聚酰亚胺镀铝膜','防原子氧布（低轨卫星用）','氟46膜（低吸收高反射）'};
type = {'b-','g-','r-','m-','k-','b:','g:','r:','m:','k:','b-.','g-.','r-.','m-.','k-.'};
for k = 1:15
for i = 65:110
pic = I{i};
point(i-64) = mean(mean(pic(position(k,2):position(k,2)+29 , position(k,1):position(k,1)+29)));
% point(i-64) = pic(position(k,2),position(k,1));
end
% subplot(4,4,k);

%%%%%%g归一化处理,归一到10000%%%
max_number = max(point);
point = point/max_number;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
plot(x,point,type{k});
hold on;
end
legend(str);
%suptitle('光谱图');
title('不同材料归一化光谱图','fontsize',12,'fontweight','b','color','k');
xlabel('波长/nm','fontsize',12,'fontweight','b','color','k');
ylabel('灰度归一化','fontsize',12,'fontweight','b','color','k');


figure;
str = {'普通布料','双面镀铝薄膜','防原子氧布（低轨卫星用）','氟46膜（低吸收高反射）'};
for k = 1:4

switch k
case 1
t = 2;
case 2
t = 11;
case 3
t = 14;
case 4
t = 15;
end
for i = 65:110
pic = I{i};
point(i-64) = mean(mean(pic(position(t,2):position(t,2)+29 , position(t,1):position(t,1)+29)));
% point(i-64) = pic(position(k,2),position(k,1));
end
if k == 1
suptitle('光谱图');
end
subplot(2,2,k);
plot(x,point,'r.-');
xlabel('波长/nm','fontsize',8,'fontweight','b','color','k');
ylabel('灰度值','fontsize',8,'fontweight','b','color','k');
legend(str{k});
end


figure;
%str = {'材料1','材料2','材料3','材料4','材料5','材料6','材料7','材料8','材料9','材料10','材料11','材料12','材料13','材料14','材料15'};
str = {'材料1','材料2','材料3','材料4','材料5','材料6','材料7','材料8','材料9','聚酰亚胺渗碳膜（防杂光）','双面镀铝薄膜','低反射率涂层绸','聚酰亚胺镀铝膜','防原子氧布（低轨卫星用）','氟46膜（低吸收高反射）'};
type = {'b-','g-','r-','m-','k-','b:','g:','r:','m:','k:','b-.','g-.','r-.','m-.','k-.'};
for k = 1:15
for i = 65:110
pic = I{i};
point(i-64) = mean(mean(pic(position(k,2):position(k,2)+29 , position(k,1):position(k,1)+29))) / mean(mean(pic(position(16,2):position(16,2)+29 , position(16,1):position(16,1)+29)));
% point(i-64) = pic(position(k,2),position(k,1));
end
if k == 1
suptitle('光谱图(与背景的比值)');
end
subplot(4,4,k);
plot(x,point,type{3});
legend(str{k});
xlabel('波长/nm','fontweight','b','fontsize',8,'color','k');
ylabel('灰度与背景的比值','fontweight','b','fontsize',8,'color','k');
end
%suptitle('光谱图');

标签：pic,point,光谱,材料,29,-.,position
来源： https://www.cnblogs.com/wolf-code/p/15385584.html