地震数据剖面图-matlab

这里主要有四个参数需要设置
wigb.m地震数据画图函数下载地址

a: seismic data 地震数据
%        scale: multiple data by scale  扫描频率【1-3】范围内任选一个数字
%        x: x-axis;  所有道的位置信息
%        z: vertical axis (time or depth)  时间
%        amx：这个参数可以不用管它

这里有一个时序信号；数据一行为所有道的采集数据，列为道数。

%加载地震数据的，也可以使用fopen()函数加载数据
load earthquak_data.mat     %加载地震数据
rec=earthquak_data(1:1:1000,13:36);


[M,N]=size(rec);%地震数据是按照时间序列存放的，比如第一次采样时，所有道的数据是一行，列为道数
dt=0.002;  %这里为采样时间间隔
t=0:dt:(M-1)*dt; %这里需要计算正确数据的时间

%下面这句是用来设置每道的道间距
for i=1:N
    often_offset(i)=2*(i-1);
end

figure
xlabel('Traces');
ylabel('Time (s)');
wigb(rec,1,often_offset,t);