-
Notifications
You must be signed in to change notification settings - Fork 5
/
test_fig5.m
166 lines (145 loc) · 7.48 KB
/
test_fig5.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
% Script to plot Figure 6
% BY Yangkang Chen
% Dec, 22, 2021
clc;clear;close all;
addpath(genpath('subroutines/'));
NOs=[1,20,10,25,11,2];
labels={... %P-arrival sample NO from the SEGY file
'FORGE\_78-32\_iDASv3-P11\_UTC190423150554.sgy',... %24169
'FORGE\_78-32\_iDASv3-P11\_UTC190426070723.sgy',... %24811
'FORGE\_78-32\_iDASv3-P11\_UTC190426062208.sgy',... %26090
'FORGE\_78-32\_iDASv3-P11\_UTC190426110008.sgy',... %4921
'FORGE\_78-32\_iDASv3-P11\_UTC190426062553.sgy',... %8934
'FORGE\_78-32\_iDASv3-P11\_UTC190423182409.sgy'}; %4210
eq=zeros(2000,960);
[n1,n2]=size(eq);
t=[0:n1]*0.0005;
ngap=50;
x=1:n2*3+2*ngap;
%% first one
ii=1;
if ~ismember(NOs(ii),[14,16,17,27,47,52])
load(strcat('mat_raw/eq-',num2str(NOs(ii)),'.mat'));
end
eq=d1;
load(strcat('mat_bpsomffk/eq-',num2str(NOs(ii)),'.mat'));
comp1=[eq,zeros(n1,ngap),d1,zeros(n1,ngap),eq-d1];
%% second one
ii=2;
if ~ismember(NOs(ii),[14,16,17,27,47,52])
load(strcat('mat_raw/eq-',num2str(NOs(ii)),'.mat'));
end
eq=d1;
load(strcat('mat_bpsomffk/eq-',num2str(NOs(ii)),'.mat'));
comp2=[eq,zeros(n1,ngap),d1,zeros(n1,ngap),eq-d1];
%% third one
ii=3;
if ~ismember(NOs(ii),[14,16,17,27,47,52])
load(strcat('mat_raw/eq-',num2str(NOs(ii)),'.mat'));
end
eq=d1;
load(strcat('mat_bpsomffk/eq-',num2str(NOs(ii)),'.mat'));
comp3=[eq,zeros(n1,ngap),d1,zeros(n1,ngap),eq-d1];
%% fourth one
ii=4;
if ~ismember(NOs(ii),[14,16,17,27,47,52])
load(strcat('mat_raw/eq-',num2str(NOs(ii)),'.mat'));
end
eq=d1;
load(strcat('mat_bpsomffk/eq-',num2str(NOs(ii)),'.mat'));
comp4=[eq,zeros(n1,ngap),d1,zeros(n1,ngap),eq-d1];
%% fifth one
ii=5;
if ~ismember(NOs(ii),[14,16,17,27,47,52])
load(strcat('mat_raw/eq-',num2str(NOs(ii)),'.mat'));
end
eq=d1;
load(strcat('mat_bpsomffk/eq-',num2str(NOs(ii)),'.mat'));
comp5=[eq,zeros(n1,ngap),d1,zeros(n1,ngap),eq-d1];
%% sixth one
ii=6;
if ~ismember(NOs(ii),[14,16,17,27,47,52])
load(strcat('mat_raw/eq-',num2str(NOs(ii)),'.mat'));
end
eq=d1;
load(strcat('mat_bpsomffk/eq-',num2str(NOs(ii)),'.mat'));
comp6=[eq,zeros(n1,ngap),d1,zeros(n1,ngap),eq-d1];
%% set xticks
xts1=[300,600,960];
xts2=xts1+ngap+n2;
xts3=xts1+ngap*2+n2*2;
xts=[xts1,xts2,xts3];
%% combined figure
figure('units','normalized','Position',[0.0 0.0 0.5, 1],'color','w');
subplot(3,2,1);das_imagesc(comp1,95,1,x,t);
ylabel('Time (s)','Fontsize',10,'fontweight','bold');
xlabel('Channel','Fontsize',10,'fontweight','bold');
set(gca,'Linewidth',2,'Fontsize',10,'Fontweight','bold');
text(n2/2,-0.05,'Raw data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap+n2,-0.05,'Denoised data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap*2+n2*2,-0.05,'Noise','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(-200,-0.1,'(a)','color','k','Fontsize',15,'fontweight','bold','HorizontalAlignment','center');
text(0.1,0.95,labels{1},'color','b','Fontsize',10,'fontweight','bold','HorizontalAlignment','left');
xticks(xts);set(gca,'xticklabel',{'300','600','960'});
subplot(3,2,2);das_imagesc(comp2,95,1,x,t);
ylabel('Time (s)','Fontsize',10,'fontweight','bold');
xlabel('Channel','Fontsize',10,'fontweight','bold');
set(gca,'Linewidth',2,'Fontsize',10,'Fontweight','bold');
text(n2/2,-0.05,'Raw data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap+n2,-0.05,'Denoised data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap*2+n2*2,-0.05,'Noise','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(-200,-0.1,'(b)','color','k','Fontsize',15,'fontweight','bold','HorizontalAlignment','center');
text(0.1,0.95,labels{2},'color','b','Fontsize',10,'fontweight','bold','HorizontalAlignment','left');
xticks(xts);set(gca,'xticklabel',{'300','600','960'});
subplot(3,2,3);das_imagesc(comp3,95,1,x,t);
ylabel('Time (s)','Fontsize',10,'fontweight','bold');
xlabel('Channel','Fontsize',10,'fontweight','bold');
set(gca,'Linewidth',2,'Fontsize',10,'Fontweight','bold');
text(n2/2,-0.05,'Raw data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap+n2,-0.05,'Denoised data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap*2+n2*2,-0.05,'Noise','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(-200,-0.1,'(c)','color','k','Fontsize',15,'fontweight','bold','HorizontalAlignment','center');
text(0.1,0.95,labels{3},'color','b','Fontsize',10,'fontweight','bold','HorizontalAlignment','left');
xticks(xts);set(gca,'xticklabel',{'300','600','960'});
subplot(3,2,4);das_imagesc(comp4,95,1,x,t);
ylabel('Time (s)','Fontsize',10,'fontweight','bold');
xlabel('Channel','Fontsize',10,'fontweight','bold');
set(gca,'Linewidth',2,'Fontsize',10,'Fontweight','bold');
text(n2/2,-0.05,'Raw data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap+n2,-0.05,'Denoised data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap*2+n2*2,-0.05,'Noise','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(-200,-0.1,'(d)','color','k','Fontsize',15,'fontweight','bold','HorizontalAlignment','center');
text(0.1,0.95,labels{4},'color','b','Fontsize',10,'fontweight','bold','HorizontalAlignment','left');
xticks(xts);set(gca,'xticklabel',{'300','600','960'});
inds1=1:400;
subplot(3,2,5);das_imagesc(comp5,95,1,x,t);
ylabel('Time (s)','Fontsize',10,'fontweight','bold');
xlabel('Channel','Fontsize',10,'fontweight','bold');
set(gca,'Linewidth',2,'Fontsize',10,'Fontweight','bold');
text(n2/2,-0.05,'Raw data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap+n2,-0.05,'Denoised data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap*2+n2*2,-0.05,'Noise','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(-200,-0.1,'(e)','color','k','Fontsize',15,'fontweight','bold','HorizontalAlignment','center');
text(0.1,0.95,labels{5},'color','b','Fontsize',10,'fontweight','bold','HorizontalAlignment','left');
annotation(gcf,'textarrow',[0.310 0.289],...
[0.252 0.296],'Color','r','TextColor','r','HorizontalAlignment','center',...
'String',{'Visible signal'},...
'LineWidth',2,...
'FontSize',10,'fontweight','bold');
xticks(xts);set(gca,'xticklabel',{'300','600','960'});
subplot(3,2,6);das_imagesc(comp6,95,1,x,t);
ylabel('Time (s)','Fontsize',10,'fontweight','bold');
xlabel('Channel','Fontsize',10,'fontweight','bold');
set(gca,'Linewidth',2,'Fontsize',10,'Fontweight','bold');
text(n2/2,-0.05,'Raw data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap+n2,-0.05,'Denoised data','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(n2/2+ngap*2+n2*2,-0.05,'Noise','color','k','Fontsize',10,'fontweight','bold','HorizontalAlignment','center');
text(-200,-0.1,'(f)','color','k','Fontsize',15,'fontweight','bold','HorizontalAlignment','center');
text(0.1,0.95,labels{6},'color','b','Fontsize',10,'fontweight','bold','HorizontalAlignment','left');
annotation(gcf,'textarrow',[0.740 0.719],...
[0.252 0.296],'Color','r','TextColor','r','HorizontalAlignment','center',...
'String',{'Visible signal'},...
'LineWidth',2,...
'FontSize',10,'fontweight','bold');
xticks(xts);set(gca,'xticklabel',{'300','600','960'});
print(gcf,'-depsc','-r300','fig5.eps');