bwtdec.asv

来自「it describes about the Bwt transform for」· ASV 代码 · 共 56 行

function bwtdec
%bwtdec is a BWT d
% type bwtenc at the command prompt and 
% a user interface is dispalyed.
% currently,its just text file.u can select just text files 
% for now.
%the BWT transformed file is written as bwt.cmp in the 
%current directory
%MA Imran [e-mail]
%please send me ur feedback and suggestion 

dec_bwt_file=fopen('bwt.cmp','r');
dec_bwt_read=fread(dec_bwt_file,'uint8');
fclose(dec_bwt_file);
disp('BWT DECODING Started');

encoded_data=dec_bwt_read(1:length(dec_bwt_read)-1);
primary_index=dec_bwt_read(length(dec_bwt_read));


sorted_data=sort(encoded_data);
vector_flag=ones(1,length(encoded_data))';
vector=zeros(1,length(encoded_data))';
%%%%%%%preparing vector table
for i=1:length(sorted_data)
    for j=1:length(sorted_data)
        if(encoded_data(j)==sorted_data(i) && vector_flag(j))
%             clc;
%             encoded_data(j);
%             sorted_data(i);
            vector_flag(j);
            vector(i)=j;
            vector_flag(j)=0;
            break
            
        end
    end
end

index=primary_index;
reconst_data=zeros(1,length(encoded_data));
%getting original data back
for i=1:length(encoded_data)
       reconst_data(i)=encoded_data(index);
       index=vector(index);
   end
  
   
%    char(lexi_sorted_data(:,length(a)));
%    primary_index=find(ind==2);
%    char(reconst_data);
   rec_file=fopen('orig_file.txt','w');
   fwrite(rec_file,reconst_data,'uint8');
   fclose(rec_file);
   disp('BWT Decosing over..');
   disp('file writeen back as orig_file.txt')