dt = 2.5e-8;
xzero = 2.495e-5;
t = [0:2007]*dt + xzero;

fid=fopen('mouseab.bin','rb','b');
raw=zeros(64,64,2008);

for k = 1:64,
    for l = 1:64,
        q = 6.1035e-6*interp(fread(fid,502,'short'),4);
        if (k == l)
            q = -0.2*[q(4:2008)-mean(q(1000:2008)) ; zeros(3,1)]; 
        else
            q = q - mean(q(1000:2008));
        end
        raw(k,l,:) = q;
    end
end

fclose(fid);