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
| clear
figure (1)
clf
load bathy_ship
imagesc (x,y,z)
axis xy
axis equal
hold on
load NavTot
plot (NavTot(:,9), NavTot(:,8), 'k.')
date0=datenum(2013,11,1,0,0,0);
datev=datenum(NavTot(:,2),NavTot(:,3),NavTot(:,4), NavTot(:,5), NavTot(:,6), NavTot(:,7))-date0;
p=[21 11 2013 2 36 43 1100
21 11 2013 3 35 25 1000
21 11 2013 5 21 20 1200
21 11 2013 5 31 40 1200
21 11 2013 5 33 19 1200
21 11 2013 5 35 55 1200
21 11 2013 5 41 13 1200
21 11 2013 6 51 06 1200
21 11 2013 6 54 31 1200
21 11 2013 7 50 59 1200
21 11 2013 8 20 50 1200
21 11 2013 2 39 09 1200
21 11 2013 4 41 30 1200
21 11 2013 5 18 40 1200
21 11 2013 7 54 21 1200
21 11 2013 7 55 17 1200
21 11 2013 8 55 17 1200
22 11 2013 4 58 19 1000
22 11 2013 5 14 30 1300
22 11 2013 5 28 25 1200
22 11 2013 6 37 37 1200
22 11 2013 6 58 25 1200
22 11 2013 7 58 25 1200
22 11 2013 8 28 25 1200
22 11 2013 9 23 40 1200
22 11 2013 4 49 27 1200
22 11 2013 5 37 22 1200
22 11 2013 5 49 27 1200
22 11 2013 6 19 27 1200
22 11 2013 8 19 27 1200
22 11 2013 9 49 27 1200
23 11 2013 1 16 35 1100
23 11 2013 1 23 28 1100
23 11 2013 4 17 56 1200
23 11 2013 5 22 46 1200
23 11 2013 5 27 52 1200
23 11 2013 5 01 06 1100
24 11 2013 15 35 30 1100
24 11 2013 16 12 26 1000
24 11 2013 17 22 52 1100
24 11 2013 17 38 32 1100
24 11 2013 18 08 37 1100
24 11 2013 18 25 59 1100
24 11 2013 20 43 08 1100
25 11 2013 7 16 36 1100
24 11 2013 16 38 03 1100
24 11 2013 17 21 43 1000
24 11 2013 17 26 25 1200
24 11 2013 0 01 20 1200]
datep=datenum(p(:,3), p(:,2), p(:,1), p(:,4), p(:,5), p(:,6))-date0;
lonn=interp1(datev,NavTot(:,9),datep);
latn=interp1(datev,NavTot(:,8),datep);
plot(lonn,latn,'g*') |
Partager