Skip to content

Instantly share code, notes, and snippets.

@dzonesasaki
Created June 24, 2014 03:26
Show Gist options
  • Save dzonesasaki/c6370a933b9503944f52 to your computer and use it in GitHub Desktop.
Save dzonesasaki/c6370a933b9503944f52 to your computer and use it in GitHub Desktop.
irr_dat read
function [power,wlength,N_lensp,N_data]=loadIrrDat(fnameL)
fid=fopen(fnameL,'rt');
Nskp=5; % for iir file is 60
%dummyStr = fscanf(fid,'%s',(Nskphdr-1)*2);
for uilp=1:Nskp
[dummyStr] = fgetl(fid);
end
dummyStr = fgetl(fid);
%[dummyStr,dummyVal] = sscanf(dummyStr,['保存開始波長[nm]' '\t' '%f'])
lambdaStart = sscanf(dummyStr(18:end),'%f');
dummyStr = fgetl(fid);
lambdaEnd = sscanf(dummyStr(18:end),'%f');
dummyStr = fgetl(fid);
lambdaStep = sscanf(dummyStr(18:end),'%f');
N_lensp = (lambdaEnd-lambdaStart)/lambdaStep+1;
dummyStr = fgetl(fid);
dummyStr = fgetl(fid);
dummyStr = fgetl(fid);% 日付け(10byte)_(1byte)時間(8byte)\t(1byte)=20byte
N_data = (length(dummyStr)-1)/20;
if (N_data - round(N_data))~=0, error('undefined file');end
Nskp=48;
for uilp=1:Nskp
[dummyStr] = fgetl(fid);
end
Ndat=1000-350+1;
dat = fscanf(fid,'%f',[ (N_data+1) ,N_lensp]);
FCL=fclose(fid);
dat=dat';
power=dat(:,(2:end));
wlength=dat(:,1);
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment