%% backscatter correction

%

b_AB=c_AB-a_AB;

N=size(b_AB,1);

for k=1:N
    b_730=b_AB(k,73);
    a_730=a_AB(k,73);
    rr=a_730./b_730;
    a_AB_corr(k,:)=a_AB(k,:)-b_AB(k,:)*rr;
end




b_AC=c_AC-a_AC;

N=size(b_AC,1);

for k=1:N
    b_730=b_AC(k,73);
    a_730=a_AC(k,73);
    rr=a_730./b_730;
    a_AC_corr(k,:)=a_AC(k,:)-b_AC(k,:)*rr;
end


b_AD=c_AD-a_AD;

N=size(b_AD,1);

for k=1:N
    b_730=b_AD(k,73);
    a_730=a_AD(k,73);
    rr=a_730./b_730;
    a_AD_corr(k,:)=a_AD(k,:)-b_AD(k,:)*rr;
end



b_AE=c_AE-a_AE;

N=size(b_AE,1);

for k=1:N
    b_730=b_AE(k,73);
    a_730=a_AE(k,73);
    rr=a_730./b_730;
    a_AE_corr(k,:)=a_AE(k,:)-b_AE(k,:)*rr;
end



