matlab读取spa格式光谱.docxVIP

matlab读取spa格式光谱

clc;clear;closeall

%%

path=E:\读研资料熟期红枣光谱\10团5连10排胶板土;

A=dir(fullfile(path,*.spa));

A=struct2cell(A);

B=A(1,:);

[m,n]=size(A);

fori=1:n

filename=strcat(path,\,B(i));

address=getStart(char(filename));%获取吸光度值的指针位置

fid=fopen(char(filename),r);

%Findthepointsnumber

fseek(fid,hex2dec(234),bof);%fseek:更改当前文件（fid）指针

%fid:表示当前文件

%hex2dec:将16进制字符转换成10进制数字

%将指针移到hex2dec(234)之后

%bof:表示从文件开始偏移

Number_of_DataPoints=fread(fid,1,int32);%fread以二进制形式，从文件读出数据

%1:读出1个数据，构成列向量

%int32:精度为32位整型数

%FindthemaximumandminimumofWavenumber(cm-1)range

fseek(fid,576,bof);

Maximum_Wavenumber=fread(fid,1,single);

Minimum_Wavenumber=fread(fid,1,single);

Interval=(Maximum_Wavenumber-Minimum_Wavenumber)/(Number_of_DataPoints-1);

Wavenumber=linspace(Minimum_Wavenumber,Maximum_Wavenumber,Number_of_DataPoints).;

Wavenumber=flipud(Wavenumber);

%FindtheY-Axisdatatype:%TransmittanceorAbsorbance

fseek(fid,hex2dec(316),bof);Y_Label=char(fread(fid,10,uchar));

%Howtodefinetheoffsetforspectraldatastillremainsunresolved.

fseek(fid,address,bof);

spectrum(:,i)=fread(fid,Number_of_DataPoints,single);%double);%float64,%real*8

end

figure(1),plot(Wavenumber,spectrum);set(gcf,color,w);

set(gca,xdir,rev,xcolor,b,ycolor,b,xlim,[round(Minimum_Wavenumber),round(Maximum_Wavenumber)]);

xlabel(Wavenumber/cm^{-1});ylabel(Y_Label);

%%保存CSV格式

%fori=1:n

%a=[wavenumber,spectrum(:,i)];

%index=B(i);

%index=char(index);

%index(:,end-3:end)=[];

%dlmwrite([E:\读研资料熟期红枣光谱\csv\,index,.csv],a,precision,19)

%end