对 segy地震数据做增益处理matlab代码.zip

close all clear clc data=readsegy('Untitled+GathEP-C.sgy'); % 构造时间增益函数 % 这里应用的是指数补偿增益 x = 1:length(data); % y = 1.1.^(x*1e-2); % y = 1.2.^(x*1e-2); y = 1.3.^(x*1e-2)-0.5; % y = 1.4.^(x*1e-2); % y = 1.5.^(x*1e-2); % 设置最大增益倍数 ymax = 25; % 限制增益倍数 y(y>ymax) = ymax; % 绘制时间增益函数图像 figure(1), plot(y,'k','lineWidth',2); axis([0 inf 0 30]); data_new = data.*y; % 绘制增益后的探地雷达图像 figure(2) imagesc(data_new); colorbar; xlabel('trace number'); ylabel('samples'); title('revised image'); % 绘制原始探地雷达图像 figure(3) imagesc(data); colorbar xlabel('trace number'); ylabel('samples'); title('original image');