-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathpsptf_ch1a_alok_pool_comp.m
123 lines (112 loc) · 3.75 KB
/
psptf_ch1a_alok_pool_comp.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
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
function psptf_ch1a_alok_pool_comp(f,cl,freq,n_contour,colorbar_flag,line_flag,mask_flag,coh_contour)
% function psptf_ch1a(f,cl,freq,n_contour,colorbar_flag,line_flag,mask_flag,coh_contour)
% Function to display time dependent coherence estimate in current subplot window.
%
% Copyright (C) 2008, 2016, David M. Halliday.
% This file is part of NeuroSpec.
%
% NeuroSpec is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2 of the License, or
% (at your option) any later version.
%
% NeuroSpec is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with NeuroSpec; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
%
% NeuroSpec is available at: http://www.neurospec.org/
% Contact: contact@neurospec.org
%
% Inputs
% f Frequency matrix from type 2 analysis, over range of offset values.
% cl cl structure from type 2 analysis.
% freq Frequency limit for plotting.
% n_contour Optional number of contours to use (default 10).
% colorbar_flag Optional flag to control drawing of colorbar - 0:No(default), 1:Yes.
% line_flag Optional flag to control drawing of contour lines - 0:No(default), 1:Yes.
% mask_flag Optional flag to mask coherence values below confidence limit - 0:No(default), 1:Yes.
% coh_contour Optional vector of contour heights for coherence.
%
% function psptf_ch1a(f,cl,freq,n_contour,colorbar_flag,line_flag,mask_flag,coh_contour)
% Revised plotting compatible with MATLAB R2014b graphics
% Check numbers of arguments.
if (nargin<3)
error('Not enough input arguments');
end
if (cl(1).type~=2)
error('Not type 2 analysis')
end
% Defaults
n_seg=length(f(1,1,:)); % No of time slice segments.
if (n_seg<3)
error('Less than 3 time offsets')
end
if (nargin<4)
n_contour=10;
end
if (nargin<5)
colorbar_flag=0;
end
if (nargin<6)
line_flag=0;
end
if nargin<7
mask_flag=0;
end
% Extract frequency range
plot_ind=find(f(:,1,1)<=freq);
% Frequency axis in Hz
freq_axis=f(plot_ind,1,1);
% Extract time axis values
dt=cl(1).dt; % Sampling interval (ms).
time_axis=[];
for ind=1:n_seg
time_axis=[time_axis; cl(ind).offset*dt];
end
% Contour values
if (nargin<8)
ch_max=max(max(f(plot_ind,4,:)));
ch_min=cl(1).ch_c95;
ch_int=(ch_max-ch_min)/(n_contour-1);
coh_contour=[0 ch_min:ch_int:ch_max];
else
coh_contour=[0 cl(1).ch_c95 coh_contour];
end
% Extract coherence
ch=squeeze(f(plot_ind,4,:));
% Mask values below confidence limit
if mask_flag
ch(find(ch<cl(1).ch_c95))=0.0; % Mask at zero
end
% Plot Time dependent Coherence.
[C,h]=contourf(time_axis,freq_axis,ch,coh_contour);
if (line_flag==0)
for ind=1:length(h)
set(h(ind),'LineStyle','none')
end
end
xlabel ('Offset (ms)')
ylabel ('Frequency (Hz)')
title('')
% Use jet colormap, no longer default from R2014b onwards
colormap(jet)
if (colorbar_flag)
H1=gca;
H=colorbar;
caxis([0 0.004]); %%%%%%%%%%%%%%% changed by alok %%%%%%%%%%%
% For older graphics can add confidence limit as line to colorbar
if verLessThan('matlab','8.4')
v=get(H);
axes(H)
for sect_ind=1:n_seg
ch_c95_vector=[v.XLim;cl(sect_ind).ch_c95,cl(sect_ind).ch_c95];
line(ch_c95_vector(1,:),ch_c95_vector(2,:),'Color','k')
end
axes(H1)
end
end