-
Notifications
You must be signed in to change notification settings - Fork 1
/
load_counts.m
69 lines (66 loc) · 1.93 KB
/
load_counts.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
data=dlmread('count.txt');
% [time, block]
maxval = max(max(max(data)))
% 82 time steps, 36 blocks, 180 degrees
%figure('Color', [0.2 0.2 0.2])
%
% figure
% colormap(hot);
% j=0
% for i=1:4:36
% j=j+1;
% subplot(3,3,j)
% %HeatMap(squeeze(data1(:,i,:))', 'LabelsWithMarkers', true, 'DisplayRange', 40, 'Colormap', redgreencmap(100, 'Interpolation', 'Sigmoid'))
% plot(data(:,i)')
% title(sprintf('block %d', i))
% xlabel('Time Step')
% ylabel('Count')
% set(gca,'YDir','normal', ...
% 'xgrid', 'on', 'xcolor',[.2 .2 .2], ...
% 'ygrid', 'on', 'ycolor',[.2 .2 .2], ...
% 'ytick', 0:10:100, ...
% 'color', 'w')
% % ylim([15, 95])
% end
% saveas(gca, 'degrees.png', 'png');
j=0
for i=1:1:36
j=j+1;
subplot(6,6,j)
%HeatMap(squeeze(data1(:,i,:))', 'LabelsWithMarkers', true, 'DisplayRange', 40, 'Colormap', redgreencmap(100, 'Interpolation', 'Sigmoid'))
plot(data(:,i)')
title(sprintf('block %d', i))
%xlabel('Time Step')
ylabel('Count')
% set(gca,'YDir','normal', ...
% 'xgrid', 'on', 'xcolor',[.2 .2 .2], ...
% 'ygrid', 'on', 'ycolor',[.2 .2 .2], ...
% 'ytick', 70:10:100, ...
% 'color', 'w')
ylim([0, maxval])
%xlim([1, 576])
end
saveas(gca, 'count1.png', 'png');
%
% %colorbar
%
% figure
% colormap(hot);
% for i=1:36
% %HeatMap(squeeze(data1(:,i,:))', 'LabelsWithMarkers', true, 'DisplayRange', 40, 'Colormap', redgreencmap(100, 'Interpolation', 'Sigmoid'))
% imagesc(squeeze(data1(i,:,:))', [0 maxval])
% title(sprintf('block %d', i))
% xlabel('Time Step')
% ylabel('Degree')
% set(gca,'YDir','normal', ...
% 'xgrid', 'on', 'xcolor',[.2 .2 .2], ...
% 'ygrid', 'on', 'ycolor',[.2 .2 .2], ...
% 'ytick', 0:10:180, ...
% 'color', 'w')
% ylim([0, maxval])
% saveas(gca, sprintf('block %02d.png', i), 'png');
% % pause
% end
%
% [m,I] = max(data1, [], 3);
% I