-
Notifications
You must be signed in to change notification settings - Fork 3
/
example_3IncAngles.m
51 lines (47 loc) · 1.96 KB
/
example_3IncAngles.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
% Script to create the plots TB H-pol for three incidence angles e.g. 20, 30, 40°
% Region min to max coordinates [47.5 7.5 50 10]: Case study Neckar Catchment.
% This script has been used to produce Figure 2 for JORS paper.
% More information: https://github.com/pablosaa/Process_SMOS_L1C
%
% Copiright 2017, Pablo Saavedra G.
%
% Geographical limits:
lat_lon_box = [47.5 7.5 50 10]; % Limits used in Figure 2 JORS paper
latlim = lat_lon_box([1,3]);
lonlim = lat_lon_box([2,4]);
% selecting DBL file to process and plot:
file_example = ['SM_OPER_MIR_SCLF1C_20150702T042618_20150702T051937_620_001_1.DBL'];
if exist(file_example,'file'),
TSF = Process_SMOSxL1C(file_example,lat_lon_box); % no MAT-file output!
else
% Open file browser GUI:
TSF = Process_SMOSxL1C;
end
idx_pol =1; % 1 for H, 2 for V
inc_ang = [35 40 45]; % Incidence Angles indexes (index=1 is 0 deg.)
N_ang = length(inc_ang);
figure;
set(gcf,'PaperPositionMode','auto','Position',[117 500 1326 395]);
for idx=1:N_ang,
ax(idx) = subplot(1,N_ang,idx);
scatter(TSF.GridPoint_Longitude,TSF.GridPoint_Latitude,250,...
TSF.TB_Fixed_IncAngle(:,inc_ang(idx)+1,idx_pol),'o','filled');
tl=title(sprintf('TB_H [K] @ \\theta_{inc} = %2d^\\circ',...
inc_ang(idx)),...
'FontName','utopia','FontSize',15);
tmp=get(ax(idx),'XTickLabel');
set(ax(idx),'XTickLabel',[tmp repmat(' °E',size(tmp,1),1)]);
tmp=get(ax(idx),'YTickLabel');
set(ax(idx),'YTickLabel',[tmp repmat(' °N',size(tmp,1),1)]);
if idx==2, xlabel('LONGITUDE','FontSize',15); end
if idx==1, ylabel('LATITUDE','FontSize',15); end
end
set(ax,'FontName','utopia','FontSize',15,'LineWidth',1,'Color','none',...
'Box','on','GridLineStyle','--','XLim',lonlim,'YLim',latlim,...
'TickDir','out','XMinorTick','on','YMinorTick','on','Layer','top',...
'LineWidth',1.5,'TickLength',[0.02 0.05],'CLim',[205 285]);
colormap(jet);
hb = colorbar('FontSize',15,'FontName','utopia');
set(hb,'Position',[0.9178 0.1164 0.0105 0.7570]);
return;
% end of script