LAB 05 - CMD SORTING

MATLAB SCRIPT LAB 5.1

load SeismicData_gain_bpf_sdecon_gain_C
[Dsort,Hsort]=ssort(Ds_gain,H);
save SeismicData_gain_bpf_sdecon_gain_sorted_C.mat Dsort Hsort

%% CMD sorting
[cmps,fold_cmp]=extracting_cmp_fold_num(Dsort,Hsort);
figure,stem(cmps,fold_cmp,'-')
xlabel('CMP numbers','FontSize',14)
ylabel('Fold','FontSize',14)
set(gca,'YMinorGrid','on')

Figure 1: the fold (number of traces per CMP) versus CMP numbers 

%% Display CMP gather
cmp_num=208; %change according to desired CMP number
[Dcmp,dt,t,cdp,jj,cmp_offset]=extracting_cmp(Dsort,Hsort,cmp_num);
scale=2;
mwigb(Dcmp,scale,cdp,t)
xlabel(['CMP:',num2str(cmp_num),''],'FontSize',14)
ylabel('Time(s)','FontSize',14)​

Figure 2: CMP gather number 208