注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

matlab fread 读取二进制文件完整demo  

2013-05-23 13:09:26|  分类: 编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
function [data,header,kstnm] = sac2mat(sacfile);
%sac2mat  loads a SAC2000 file into Matlab, along with header information.
%  [data,header,kstnm] = sac2mat(sacfile);
%
%  sacfile must be a STRING, eg.  [fff,hhh,name]=sac2mat('dir1/file1.sac')
%
%Output:
%data= the data from the file
%header=the first 105 header variables, as described on the SAC Home Page
%kstnm= KSTNM, the station name (string)
% (header and kstnm are optional)
%
%Some common/useful headers:
%header(1:3) DELTA,DEPMIN,DEPMAX
%header(57) DEPMEN
%header(6:7) B,E
%header(8) Event origin time (seconds relative to reference time)
%header(32:39) STLA,STLO,STEL,STDP,EVLA,EVLO,EVEL,EVDP
%header(51:54) DIST,AZ,BAZ,GCARC        
%header(71:76)  NZYEAR,NZDAY,NZHOUR,NZMIN,NZSEC,NZMSEC
%header(80) NPTS

% Example: [data,header,station_name]=sac2mat('1409094992.sac'); 

sacfid = fopen(sacfile,'r');
if sacfid==-1;error='Not a valid path.'
return;end

%The first 70 header variables are floating point
header1 = fread(sacfid,70,'float32');

%The next 35 are integers (after that they're a mix of strings)
header2 = fread(sacfid,35,'int32');

%Outputs the number of points being read to the screen
header=[header1;header2];
npts = header(80)

%kstnm, the station name, is a string stored in the 110th header variable
%kstnm = KSTNM (First 3 letters of the station name)
fseek(sacfid,4*110,-1);
kstnm = char(fread(sacfid,3)');

%Now read the data...
fseek(sacfid,158*4,-1);
data(1:npts) = fread(sacfid,npts,'float32');
status = fclose(sacfid);
return
  评论这张
 
阅读(849)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016