dispcvpavg.m

基于matlab的反演程序,用于地球物理勘探中射线追踪及偏移成像程序.

function dispcvpavg
% Display of all the Cross Over point averages (left and rigth) for each shot
% with the arrival time curves 
f=gcf;
fbtime=refdata('get','fbtime');
fbcoord=refdata('get','fbcoord');
nshots = refdata('get','nshots');
cvpi = refdata('get', 'cvpi');
cvpj = refdata('get', 'cvpj');
shotcoord=refdata('get','shotcoord');

cvpavg=refdata('get','cvpavg');
% Call the Cross Over point averaging function
if( length(cvpavg) == 0 )
[cvpavg, cvpstd, cvpfold] = avgcvp(cvpi, cvpj, nshots);

figure('menubar','none');
hold on;
for n=1:nshots
	plot(fbcoord(n,:),fbtime(n,:))
	valid=find(~isnan(fbtime(n,:)));
    if (isnan(cvpavg(n,1)) ~=1)
	timei=interp1(fbcoord(n,valid),fbtime(n,valid),cvpavg(n,1));
	plot(cvpavg(n,1),timei,'color','c','linestyle','o')
    end
   if (isnan(cvpavg(n,2)) ~=1)
 	timej=interp1(fbcoord(n,valid),fbtime(n,valid),cvpavg(n,2));
	plot(cvpavg(n,2),timej,'color','r','linestyle','*')
    end
end

% Addition of the shot number to the display
xy=axis;
t=xy(4)-xy(3);
d=t/40;

for n=10:10:nshots
  str=sprintf('%d',n); 
  text(shotcoord(n),xy(3)+d,str)
end
text(xy(1)+100,xy(3)+2.5*d,'shot number')
xlabel('Coordinate (m)');
ylabel('Traveltime (ms)');
title('Traveltime for all shots with their corresponding cross over point averages (left: blue, rigth: red)');
set(gcf,'units','pixels','position',[0 0 864 576]);
figure(f); set(gcf,'menubar','none');