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

function [data] = PlotCrossCorr6( wave, S );

waveFilt = BudapestbandPass(wave,1000,5,200);
[waveCorr,lags] = xcov(waveFilt(4201:4700,2:65),200,'coeff');
%[waveCorr,lags] = xcorr(waveFilt(1:100,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(lags,waveCorr(:,64.*62+20)); axis([-10 10 -1 1]); grid on;
%figure
for n=1:64
subplot(8,8,n); plot(lags, waveCorr(:,64.*(n-1)+28)); axis([-10 10 -1 1]); grid on; %axis off;
end
%figure
%plot(lags,waveCorr(:,64.*62+20)); axis([-100 100 -1 1]); grid on;
%figure
%for n=1:64
%subplot(8,8,n); plot(lags, waveCorr(:,64.*(n-1)+20)); axis([-100 100 -1 1]); grid on; %axis off;
%end

for n=1:64
[C(n),I(n)]=max(waveCorr(:,64*(n-1)+28))
end
J=I-I(28);
L=reshape(J,8,8)'
figure
interp_L=interp2(L,4,'linear');
imagesc(interp_L); colorbar('horiz'); %axes off;
axis('image')
axis('off')



   
 