PGC_pgcMATLAB_pgcseismic_pgc_seismicdata

clc;clear %% PGC Test shot=altreadsegy('shot20.sgy'); %reading the sht gather r=zeros(1000,2400); %zeros matrix for replacing m=shot(1:50,1:1800); %m matrix (part of shot as ref Window m2=shot(1:100,1800:2400); %% Matrices for scaning the shot matrix f=shot(50:100,1:1800); n=shot(100:300,1:1800); p=shot(300:500,1:1800); q=shot(500:800,1:1800); l=shot(800:1000,1:1800); k=shot(1:500,1800:2400); h=shot(500:1000,1800:2400); %% calculating the RMSs and using some constant numbers for efficient gain z1=mean(rms(m))/mean(rms(n)); n=n*z1*5; z2=mean(rms(m))/mean(rms(p)); p=p*z2; z3=mean(rms(m))/mean(rms(q)); q=q*z3; z4=mean(rms(m))/mean(rms(l)); l=l*z4; z5=mean(rms(m2))/mean(rms(k)); k=k*z5*3; z6=mean(rms(m))/mean(rms(h)); h=h*z6; %% Replacing the results r(1:50,1:1800)=m*5; r(50:100,1:1800)=f; r(100:300,1:1800)=n; r(300:500,1:1800)=p; r(500:800,1:1800)=q; r(800:1000,1:1800)=l; r(1:500,1800:2400)=k; r(500:1000,1800:2400)=h; %% Printing the result and comparing them figure(1) subplot(2,1,1); imagesc(abs(shot)); title('Before PGC'); subplot(2,1,2); imagesc(abs(r)); title('After PGC');