sub_bmode.m

医学超声图像获取后data数据转换为二维图像的源代码

function BforDisplay = sub_BMode(Signal, SecondSample)

Dimension = size(Signal);
LineNumb = Dimension(1);
DotNumb =  Dimension(2);
TrueNumb = fix(DotNumb/SecondSample);
BforDisplay = zeros(LineNumb, TrueNumb);

for line = 1:LineNumb
    for dot = 1:TrueNumb
        BforDisplay(line,dot) = log(1+mean(Signal(line,fix((dot-1)*SecondSample)+1: fix(dot*SecondSample))));
    end
end


% 
% 
% for Index = 0 :TIME-1                                % TIME次独立采样    
%     PreFile = [DataFileDir, '\pre\sample_FB_0000',num2str(Index),'.dat'];
%     fid1=fopen(PreFile);
%     PostFile = [DataFileDir, '\post3\sample_FB_0000',num2str(Index),'.dat'];
% 
%     fid2=fopen(PostFile);
%     f1=fread(fid1,'int16');
%     f2=fread(fid2,'int16');
%     if Index ==0
%         N = length(f1);
%         p1 = zeros(size(f1));
%         p2 = zeros(size(f2));
%         ComputeLen = N/NWid;
%     end
%     p1 = p1+f1;
%     p2 = p2+f2;
%     fclose(fid1);
%     fclose(fid2); 
% end
% clear f1; clear f2;
% 
% p1 = p1./TIME;
% p2 = p2./TIME;
% 
% PreSignal=reshape(p1,ComputeLen,NWid);                % The precompression signal matrix
% PostSignal=reshape(p2,ComputeLen,NWid);               % The postcompression signal matrix
% clear p1;
% clear p2;
% 
% 
% 
% B = sub_BMode((PreSignal(1: 3500, :))', 10); 
% 
% figure(5)
% imagesc(B), colormap(gray);