关于USGS 共享光谱库读取问题

欢迎学习交流！
邮箱： z…@1…6.com
网站： https://zephyrhours.github.io/

一、引言

最近有网友发邮件咨询： 如何读取USGS光谱库的不同物质的光谱？。再在这之前，其实之前我也没有接触过这些相关的光谱共享数据库，然后闲暇之余，大概看了一下，其实很简单，只要我们知道其中的数据构成和光谱文件的数据个数，其实就完全科研将光谱数据读取出来，下面将读取光谱的代码共享给各位科研朋友！

二、文件构成

下面是USGS光谱库的官方网址，如果感兴趣的朋友可以仔细看一下相关说明，这里不在过多赘述。

下面是笔者下载下来的光谱数据解压后的文件构成，读取的时候，直接以ASCIIdata数据为准，下面我们直接进入文件夹，以下面数据为例子进行展示：

三、光谱读取和显示

下面是笔者利用MATLAB读取和显示的代码，需要的朋友自取，具体代码和显示效果如下：

% Author: Zephyr Hou
% Time: 2022-04-20
% MATLAB code

filename='C:\Users\dream\Desktop\splib07a_Alizarin_crimson_(dk)_GDS780_ASDFRa_AREF.txt';
spectralfile = readtable(filename, 'HeaderLines',1);
spectral=table2array(spectralfile);

spectral = (spectral-min(spectral(:)))/(max(spectral(:))-min(spectral(:)))

figure;
plot(spectral);
xlabel('Bands') 
ylabel('Relative spectral response') 
title('splib07a')

显示效果如下：