/
sc_plotallclusters.m
95 lines (63 loc) · 3.34 KB
/
sc_plotallclusters.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
function sc_plotallclusters(features,mua,s_opt)
psize=0.65;
%plot(mua.ts_spike+1.5,mua.waveforms./1000);
xpos=[0 0 0 1 1 1 2 2 2];
ypos=[1 2 3 1 2 3 1 2 3];
for i=1:features.Nclusters
xo=(xpos(i)*(psize+.01))+.05;
yo=-(ypos(i)*(psize+.01))+1;
%im=-sqrt(min(features.clusterimages(:,:,i),70));
im=-((features.clusterimages(:,:,i)./max(max(features.clusterimages(:,:,i))) ).^(.3));
imagesc( linspace(1,1+psize,features.imagesize)+xo , linspace(0,psize,features.imagesize)+yo , im );
% plot ISI distribution
try numel(features.isiplots); % update the isi plots in case this is an older file
catch
features=sc_updateclusterimages(features,mua,s_opt);
end;
if numel(features.isiplots{i})>1
stairs(linspace(0,psize,features.isioptions(1).nbins)+xo+1,(features.isiplots{i}.*0)+yo-psize+1.02-(psize/2),'k','color',[1 1 1].*.8,'LineWidth',2);
stairs(linspace(0,psize,features.isioptions(1).nbins)+xo+1,(features.isiplots{i}/4)+yo-psize+1.02-(psize/2),'k','color',[.9 .9 1].*.6,'LineWidth',2);
end;
% text(xpos,ypos-psize-.1,'0');
% text(xpos+psize/2,ypos-psize-.1,'ms');
% text(xpos+psize-0.1,ypos-psize-.1,num2str(features.isioptions(1).tmax));
text(xo+1.01,yo+0.02,num2str(i),'color',[0 0 0]);
plot(xo+1.06,yo+0.03,features.clusterfstrs{i},'MarkerSize',22,'color',features.colors(i,:));
if i==1
text(xo+1.1 ,yo+0.02,['N: ',num2str(sum(features.clusters==i)),' (MUA/null cluster)'],'color',[0 0 0]);
else
if features.clusterlabels(i)>1 % if there is a label
text(xo+1.1 ,yo+0.02,['N: ',num2str(sum(features.clusters==i)),' ',features.labelcategories{features.clusterlabels(i)}],'color',[0 0 0]);
else
text(xo+1.1 ,yo+0.02,['N: ',num2str(sum(features.clusters==i)),' ',features.labelcategories{features.clusterlabels(i)}],'color',[0 0 0]);
text(xo+1.35 ,yo+0.02,'(no tag)', 'color',[0 0 0],'BackgroundColor',[1 .9 .7]);
end;
end;
if i>0 % plot options symbol
plot( [1.0 0.9]+xo+psize , [psize-0.1 psize]+yo,'k'); % diagonal line
text(xo+0.96+psize,yo+psize-0.02,'+');
end;
if i>1 % plot tag symbol
plot( [1 1.1]+xo , [psize-0.1 psize]+yo,'k'); % diagonal line
text(xo+1.01,yo+psize-0.02,'tag');
end;
c=features.colors(i,:);
plot( [1 1]+xo , [0 psize]+yo,'k','color',c);
plot( [1+psize 1+psize]+xo , [0 psize]+yo,'k','color',c);
plot( [1 1+psize]+xo , [0 0]+yo,'k','color',c);
plot( [1 1+psize]+xo , [psize psize]+yo,'k','color',c);
if features.highlight>0
%npoints=numel(mua.ts_spike); % crashes with tetrode data
npoints=size(mua.waveforms,2);
xa= (linspace(0,psize,npoints));
ya= ( mua.waveforms(features.highlight,:)) *features.waveformscale .*psize;
if features.plotgroup==1 % display number of wveforms
for j=1:10
yb=( mua.waveforms(features.highlight_multiple(j),:)) *features.waveformscale .*psize;
plot(xa+xo+1,yb+yo+(psize/2),'k--','color',[.1 .1 .1]);
end;
end;
plot(xa+xo+1,ya+yo+(psize/2),'color',features.colors( features.clusters(features.highlight) ,:),'LineWidth',2);
end;
end;
colormap gray;