%filters data and plots power spectrum for 64 electrode array
%requires raw data 

function [waveCorr, time] = PlotAutoCorr4(wave, time);
%extra=[0];
rev=fliplr(time);
t=cat(2,-rev,time);
t=t';
waveFilt = BudapestbandPass(wave(:,2:65),1000,10,40,1000);
[waveCorr,lags] = xcov(waveFilt(200:700,1:64),200,'coeff');
%[waveCorr,lags] = xcorr(waveFilt(200:700,2:65),200,'coeff');

%Lrotate=ones(8,8);
%for m=1:8
 %   for n=1:8
  %      Lrotate(m,n)=n*8-m+1;
  % end
  %end
%Lrotate=reshape(Lrotate,1,64);

%Rrotate=ones(8,8);
%for n=1:8
 %   for m=1:8
  %      Rrotate(m,n)=8*(8-m)+n;
  % end
  %end
%Rrotate=reshape(Rrotate,1,64);

%figure
%plot(time,waveCorr(:,64.*12+12)); axis([-0.1 0.1 -1 1]);time(800:1200)
figure
for n=1:64
%subplot(8,8,n); plot(lags, waveCorr(:,64*(n-1)+n)); axis([-50 50 -1 1]); axis off;
subplot(8,8,n); plot(lags, waveCorr(:,64*(n-1)+n),('r')); axis([-50 50 -1 1]); axis off;

end
