脑瘫了,下了这么一个玩意,代码都是屁来的,还想参考下
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));