plot_intensity

 plot_intensity.m

 Plot timecourse of intensity

0001 % plot_intensity.m
0002 %
0003 % Plot timecourse of intensity
0004 
0005 dims = size(V);
0006 
0007 sum_int_V = sum(reshape(V,prod(dims(1:3)),dims(4)));
0008 
0009 pct_diff_max_min = (max(sum_intensity) - min(sum_intensity))/max(sum_intensity)*100;
0010 disp(sprintf('Pct difference between max and min: %1.2f', pct_diff_max_min))
0011 
0012 nplot = 5;
0013 
0014 subplot(nplot,1,1)
0015 plot(sum_intensity/mean(sum_intensity))
0016 
0017 subplot(nplot,1,2)
0018 plot(mean_intensity)
0019 
0020 subplot(nplot,1,3)
0021 %intens_rel_mean = (mean_intensity./sum_intensity);
0022 intens_rel_mean = (sum_intensity./mean_intensity);
0023 plot(intens_rel_mean)
0024 
0025 subplot(nplot,1,4)
0026 plot(intens_rel_mean/mean(intens_rel_mean))
0027 
0028 subplot(nplot,1,5)
0029 vox(1:dims(4)) = V(30,30,20,:);
0030 plot(vox/mean(vox))
0031 
0032 
0033