| 3 |
michaesp |
1 |
% ----------------------------------------------------------------------
|
|
|
2 |
% Plot a set of trajectories
|
|
|
3 |
% ----------------------------------------------------------------------
|
|
|
4 |
|
|
|
5 |
% First image (otherwise first image is not correctly written)
|
|
|
6 |
filename = 'test.eps';
|
|
|
7 |
close;
|
|
|
8 |
fh=figure('Units','pixels','Position',[100 100 900 900])
|
|
|
9 |
set(gcf, 'PaperPosition', [2 1 15 10]);
|
|
|
10 |
print('-depsc2','-r0',filename);
|
|
|
11 |
|
|
|
12 |
% Read the trajectory file
|
|
|
13 |
tra = read_trajectory('trajectory.lsl');
|
|
|
14 |
|
|
|
15 |
% Open a new figure and set the geographical projection and region
|
|
|
16 |
figure(1);
|
|
|
17 |
clf;
|
|
|
18 |
load coast
|
|
|
19 |
h=axesm('MapProjection','stereo','origin',[ 90 70 ]);
|
|
|
20 |
gridm;
|
|
|
21 |
h=plotm(lat,long,'Color','k','LineWidth',1.5)
|
|
|
22 |
%axis([-1.5 1.5 -2 0]);
|
|
|
23 |
|
|
|
24 |
% Select a set of trajectories
|
|
|
25 |
tra.select = (tra.lon > 6) & (tra.lon < 10) & (tra.lat > 44) & (tra.lat < 49) & (tra.time == 0);
|
|
|
26 |
tra.select = tra.label(tra.select);
|
|
|
27 |
|
|
|
28 |
|
|
|
29 |
% Plot the trajectories
|
|
|
30 |
col.shift = 47;
|
|
|
31 |
col.name = [ pwd '/col_table.txt' ];
|
|
|
32 |
col.spacing = 900:-50:200;
|
|
|
33 |
col.field = 'p';
|
|
|
34 |
col.orientation = 'vert';
|
|
|
35 |
plot_trajectory(tra,col,[0 -48 -96]);
|
|
|
36 |
|
|
|
37 |
% Save figure
|
|
|
38 |
filename = 'trajectory.eps';
|
|
|
39 |
set(gcf, 'PaperPosition', [2 1 15 10]);
|
|
|
40 |
print('-depsc2','-r0',filename);
|