dt = 4e-8;
xzero = 0;
t = [0:2502]*dt;

fid=fopen('rat_liver.bin','rb','b');
raw=zeros(64,64,2502);

load bpf2;
m = (length(bpf)-1)/2;
for k = 1:64,
k
    for l = 1:64,
        q = 2.4414e-6*fread(fid,2502,'short');
        if (k == l)
            q = -0.65*[q(3:2502)-mean(q(1000:2502)) ; zeros(2,1)]; 
        else
            q = q - mean(q(1000:2502));
        end
        q = conv(q,bpf);
        raw(k,l,:) = q(1+m:2502+m);
    end
end

fclose(fid);


% Dip filter
pitch = 315e-6;
cassumed = 1500;
rho1 = 30*pi/180;
rho2 = -30*pi/180;
[f1,f2]=freqspace(41,'meshgrid');
Hd = zeros(41,41);
Hd(f1 > 0 & f2 < f1*pitch*sin(rho1)/(cassumed*dt) & f2 > f1*pitch*sin(rho2)/(cassumed*dt)) = 1;
Hd(f1 == 0 & f2 == 0) = 1;
Hd(f1 < 0 & f2 > f1*pitch*sin(rho1)/(cassumed*dt) & f2 < f1*pitch*sin(rho2)/(cassumed*dt)) = 1;
dipfilt = fwind1(Hd,hamming(41));

for k = 1:64,
    k
    raw(k,:,:) = filter2(dipfilt,squeeze(raw(k,:,:)));
end