MATLAB SCRIPT 1.3
load SeismicData_C.mat
whos
%% Extract geometrical acquisition parameter
load ('SeismicData_C.mat','H')
[sx,sy,gx,gy,shot_gathers,num_trace_per_sg,sz,gz] = extracting_geometry(H);
%% Wiggle display shot number 4 to 6
load SeismicData_C.mat
shot_num=4:6;
p=0;
[Dshot,dt,dx,t,offset]=extracting_shots(D,H,shot_num ,p);
scale=1;
whos
%% Extract geometrical acquisition parameter
load ('SeismicData_C.mat','H')
[sx,sy,gx,gy,shot_gathers,num_trace_per_sg,sz,gz] = extracting_geometry(H);
%% Wiggle display shot number 4 to 6
load SeismicData_C.mat
shot_num=4:6;
p=0;
[Dshot,dt,dx,t,offset]=extracting_shots(D,H,shot_num ,p);
scale=1;

Figure 1: variable area display shot number 4 to 6
%% Display grey scaled variable density display
mwigb(Dshot,scale,offset,t)
xlabel('Offset(ft)','FontSize',14)
ylabel('Time(s)','FontSize',14)

Figure 2: grey scaled density display shot no 4 to 6
%% Display color variable density display
figure,simage_display(Dshot,offset,t,0)
xlabel('Offset(ft)','FontSize',14)
ylabel('Time(s)','FontSize',14)

Figure 3: colored variable density display shot no 4 to 6

Figure 4: wiggle on top of colored plotting area display shot gather no 4 to 6
No comments:
Post a Comment