Commit c9e22a3e authored by tilman.metz's avatar tilman.metz

Modified script for toy-example seismograms

Updated the seismograms in documentation
with seismograms from the updated script
and cleanup of old figures
parent f5db1d79
......@@ -81,7 +81,7 @@ For a comparison of observed and inverted data, it is necessary to apply a lowpa
\begin{figure}[h!]
\begin{center}
\includegraphics[width=\textwidth]{fig_toy/seismo_trace_toy}
\caption[Toy example - observed, initial and inverted seismograms]{Multi-component seismograms exemplarily for source at ($x_s$, $z_s$, $y_s$) = (32\,m, 96\,m, 92\,m) and receiver at ($x_r$, $z_r$, $y_r$) = (56\,m, 40\,m, 24\,m), normalised to one trace a) initial vs. observed data filtered below 200\,Hz and b) observed vs. inverted data filtered below 320\,Hz} \label{fig:toy_seismo1}
\caption[Toy example - observed, initial and inverted seismograms]{Multi-component seismograms exemplarily for source at ($x_s$, $z_s$, $y_s$) = (32\,m, 96\,m, 92\,m) and receiver at ($x_r$, $z_r$, $y_r$) = (48\,m, 40\,m, 24\,m), normalised to one trace in each subfigure. left) initial vs. observed data filtered below 200\,Hz and right) observed vs. inverted data filtered below 320\,Hz} \label{fig:toy_seismo1}
\end{center}
\end{figure}
\subsubsection*{Discussion}
......
......@@ -2,20 +2,23 @@
%receiver and x-, y- and z-component
%Input: binary format
clear all; clc;
close all
tracenum=50;
file_inp1='../par/su/cal_toy_vx_it1.su.shot4'; % initial data
file_inp4='../par/su_obs/obs_toy_vx_it1.su.shot4'; %observed data
file_inp4='../par/su_obs/obs_toy_vx.su.shot4.filt_200'; %filtered with lpb 200 Hz
file_inp10='../par/su_obs/obs_toy_vx.su.shot4.filt_320'; %filtered with lpb 320 Hz
file_inp7='../par/su/cal_toy_vx_it60.su.shot4'; % inverted data
file_inp2='../par/su/cal_toy_vy_it1.su.shot4';
file_inp5='../par/su_obs/obs_toy_vy_it1.su.shot4';
file_inp5='../par/su_obs/obs_toy_vy.su.shot4.filt_200';
file_inp11='../par/su_obs/obs_toy_vy.su.shot4.filt_320'; %filtered with lpb 320 Hz
file_inp8='../par/su/cal_toy_vy_it60.su.shot4';
file_inp3='../par/su/cal_toy_vz_it1.su.shot4';
file_inp6='../par/su_obs/obs_toy_vz_it1.su.shot4';
file_inp6='../par/su_obs/obs_toy_vz.su.shot4.filt_200';
file_inp12='../par/su_obs/obs_toy_vz.su.shot4.filt_320'; %filtered with lpb 320 Hz
file_inp9='../par/su/cal_toy_vz_it60.su.shot4';
......@@ -28,7 +31,9 @@ tr6 = su2matlab(file_inp6);
tr7 = su2matlab(file_inp7);
tr8 = su2matlab(file_inp8);
tr9 = su2matlab(file_inp9);
tr10 = su2matlab(file_inp10);
tr11 = su2matlab(file_inp11);
tr12 = su2matlab(file_inp12);
nt = tr1.ns; % number of time samples
......@@ -39,36 +44,48 @@ ns = tr1.ns; % number of samples per trace
fig=55;
%--------------------------------------------------------------------------
trace1=tr1(1).trace;
trace2=tr2(1).trace;
trace3=tr3(1).trace;
trace4=tr4(1).trace;
trace5=tr5(1).trace;
trace6=tr6(1).trace;
trace7=tr7(1).trace;
trace8=tr8(1).trace;
trace9=tr9(1).trace;
recx=tr2(tracenum).gx
recy=tr2(tracenum).gy
recx=tr5(tracenum).gx
recy=tr5(tracenum).gy
soux=tr2(tracenum).sx
souy=tr2(tracenum).sy
soux=tr5(tracenum).sx
souy=tr5(tracenum).sy
trace1=tr1(tracenum).trace;
trace2=tr2(tracenum).trace;
trace3=tr3(tracenum).trace;
trace4=tr4(tracenum).trace;
trace5=tr5(tracenum).trace;
trace6=tr6(tracenum).trace;
trace7=tr7(tracenum).trace;
trace8=tr8(tracenum).trace;
trace9=tr9(tracenum).trace;
trace10=tr10(tracenum).trace;
trace11=tr11(tracenum).trace;
trace12=tr12(tracenum).trace;
t=(dt:dt:ns*dt)./10^6; % time in seconds
linewidth=1;
figure(fig)
subplot(1,2,1)
plot(t,trace1/max(trace6)+4.8,'b-','LineWidth',1.5);
plot(t,trace1/max(trace5)+4.8,'b-','LineWidth',linewidth);
hold on
plot(t,trace4/max(trace6)+4.8,'k-','LineWidth',1.5);
plot(t,trace4/max(trace5)+4.8,'k-','LineWidth',linewidth);
hold on
plot(t,(trace2/max(trace6))+2.4,'b-','LineWidth',1.5);
plot(t,(trace2/max(trace5))+2.4,'b-','LineWidth',linewidth);
hold on
plot(t,(trace5/max(trace6))+2.4,'k-','LineWidth',1.5);
plot(t,(trace5/max(trace5))+2.4,'k-','LineWidth',linewidth);
hold on
plot(t,trace3/max(trace6),'b-','LineWidth',1.5);
plot(t,trace3/max(trace5),'b-','LineWidth',linewidth);
hold on
plot(t,trace6/max(trace6),'k-','LineWidth',1.5);
plot(t,trace6/max(trace5),'k-','LineWidth',linewidth);
xlabel('time in s');
xlim([0.015 0.05]);
......@@ -76,68 +93,59 @@ ylabel('normalized amplitude');
ylim([-2 7]);
set(gca,'ytick',[])
set(get(gca,'Ylabel'),'FontSize',12);
set(get(gca,'Ylabel'),'FontSize',10);
set(get(gca,'Ylabel'),'FontWeight','normal');
set(get(gca,'Xlabel'),'FontSize',12);
set(get(gca,'Xlabel'),'FontSize',10);
set(get(gca,'Xlabel'),'FontWeight','normal');
set(gca,'FontSize',12);
set(gca,'FontSize',10);
set(gca,'FontWeight','normal');
set(gca,'Linewidth',1.0);
legend('starting','observed','Location','NorthWest');
txt1 = 'x-component';
text(0.04,5.2,txt1);
text(0.038,5.3,txt1,'FontSize',8);
txt2 = 'y-component';
text(0.04,2.7,txt2);
text(0.038,2.9,txt2,'FontSize',8);
txt3 = 'z-component';
text(0.04,0.3,txt3);
text(0.038,0.5,txt3,'FontSize',8);
subplot(1,2,2)
plot(t,trace4/max(trace6)+4.8,'k-','LineWidth',1.5);
plot(t,trace10/max(trace8)+4.8,'k-','LineWidth',linewidth);
hold on
plot(t,trace7/max(trace6)+4.8,'r-','LineWidth',1.5);
plot(t,trace7/max(trace8)+4.8,'r-','LineWidth',linewidth);
hold on
plot(t,(trace5/max(trace6))+2.4,'k-','LineWidth',1.5);
plot(t,(trace11/max(trace8))+2.4,'k-','LineWidth',linewidth);
hold on
plot(t,(trace8/max(trace6))+2.4,'r-','LineWidth',1.5);
plot(t,(trace8/max(trace8))+2.4,'r-','LineWidth',linewidth);
hold on
plot(t,trace6/max(trace6),'k-','LineWidth',1.5);
plot(t,trace12/max(trace8),'k-','LineWidth',linewidth);
hold on
plot(t,trace9/max(trace6),'r-','LineWidth',1.5);
plot(t,trace9/max(trace8),'r-','LineWidth',linewidth);
xlabel('time in s');
xlim([0.015 0.05]);
ylabel('normalized amplitude');
ylim([-2 7]);
set(get(gca,'Ylabel'),'FontSize',12);
set(get(gca,'Ylabel'),'FontSize',10);
set(get(gca,'Ylabel'),'FontWeight','normal');
set(get(gca,'Xlabel'),'FontSize',12);
set(get(gca,'Xlabel'),'FontSize',10);
set(get(gca,'Xlabel'),'FontWeight','normal');
set(gca,'FontSize',12);
set(gca,'FontSize',10);
set(gca,'FontWeight','normal');
set(gca,'Linewidth',1.0);
set(gca,'ytick',[]);
hleg1 = legend('observed','inverted','Location', 'NorthWest');
txt1 = 'x-component';
text(0.04,5.2,txt1);
text(0.038,5.3,txt1,'FontSize',8);
txt2 = 'y-component';
text(0.04,2.7,txt2);
text(0.038,2.9,txt2,'FontSize',8);
txt3 = 'z-component';
text(0.04,0.3,txt3);
text(0.038,0.5,txt3,'FontSize',8);
% find rec, source position:
% SU=su2matlab('/data14/sdunkl/toy/results1/su/cal_toynew_vy_it43.su.shot3')
% recx=SU(tracenum).gx
% recy=SU(tracenum).gx
% soux=SU(tracenum).sx
% souy=SU(tracenum).sy
% exportfig(fig, ['seismo_trace_toy.eps'],'bounds','tight', 'color','rgb', ...
% 'preview','none', 'resolution',200, 'lockaxes',1);
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment