[height,width]=size(data);
period = [];
for i = 1:width
    signal = data(:,i)-mean(data(:,i));
    [corr,shift] = xcorr(signal,8192*60,'coeff');
    figure
    p = plot(shift/8192,corr);
    [pksh,lcsh] = findpeaks(corr,'MinPeakDistance',3000);
    short = mean(diff(lcsh));
    [pklg,lclg] =findpeaks(corr,'MinPeakDistance',ceil(short),'MinPeakHeight',0.07);
    hold on
    peaks = plot(shift(lclg)/8192,pklg,'or');
    hold off
    period(i) = mean(diff(lclg))/8196;
end