0

我在 Matlab 中有一个关于图像压缩的项目。到目前为止,我已经成功地实现了对图像的霍夫曼编码,这给了我一个二进制代码的向量。之后我运行霍夫曼解码,我得到一个包含压缩图像元素的向量。我的问题是我可以找到如何从这个向量中重建图像并创建图像文件。

任何帮助将不胜感激

更新

基于 Ben A. 的帮助,我取得了进展,但我仍然有一些问题。更准确地说。我有一个图像矩阵。在此图像矩阵上找到唯一符号(元素)后,我计算概率,然后使用此函数:

function [h,L,H]=Huffman_code(p,opt)
% Huffman code generator gives a Huffman code matrix h, 
%  average codeword length L & entropy H
% for a source with probability vector p given as argin(1) 
zero_one=['0'; '1']; 
if nargin>1&&opt>0, zero_one=['1'; '0']; end
if abs(sum(p)-1)>1e-6
  fprintf('\n The probabilities in p does not add up to 1!');
end  
M=length(p);  N=M-1; p=p(:); % Make p a column vector
h={zero_one(1),zero_one(2)};
if M>2
  pp(:,1)=p;
  for n=1:N
     % To sort in descending order
     [pp(1:M-n+1,n),o(1:M-n+1,n)]=sort(pp(1:M-n+1,n),1,'descend'); 
     if n==1, ord0=o; end  % Original descending order
     if M-n>1, pp(1:M-n,n+1)=[pp(1:M-1-n,n); sum(pp(M-n:M-n+1,n))]; end
  end
  for n=N:-1:2
     tmp=N-n+2; oi=o(1:tmp,n);
     for i=1:tmp, h1{oi(i)}=h{i}; end
     h=h1;   h{tmp+1}=h{tmp};
     h{tmp}=[h{tmp} zero_one(1)]; 
     h{tmp+1}=[h{tmp+1} zero_one(2)];
  end
  for i=1:length(ord0), h1{ord0(i)}=h{i}; end
  h=h1;
end
L=0; 
for n=1:M, L=L+p(n)*length(h{n}); end  % Average codeword length
H=-sum(p.*log2(p)); % Entropy by Eq.(9.1.4)

我计算图像的霍夫曼代码。

现在我使用这个功能:

function coded_seq=source_coding(src,symbols,codewords)
% Encode a data sequence src based on the given (symbols,codewords).
no_of_symbols=length(symbols);  coded_seq=[];
if length(codewords)<no_of_symbols
  error('The number of codewords must equal that of symbols');
end
for n=1:length(src)
   found=0;
   for i=1:no_of_symbols
      if src(n)==symbols(i), tmp=codewords{i}; found=1; break; end
   end
   if found==0, tmp='?'; end
   coded_seq=[coded_seq tmp];
end

在 src 我把我的图像(矩阵)放在哪里,我得到了我的图像的编码序列。

最后是这个功能:

function decoded_seq=source_decoding(coded_seq,h,symbols)
% Decode a coded_seq  based on the given (codewords,symbols).
M=length(h);  decoded_seq=[];
while ~isempty(coded_seq)
  lcs= length(coded_seq); found=0;
  for m=1:M
    codeword= h{m};
    lc= length(codeword);
    if lcs>=lc&codeword==coded_seq(1:lc)
      symbol=symbols(m); found=1;  break;
    end
    if found==0, symbol='?'; end
  end 
  decoded_seq=[decoded_seq symbol]; 
  coded_seq=coded_seq(lc+1:end); 
end

用于解码编码序列。问题是,最后作为编码序列,我得到一个 1x400 矩阵,我应该得到一个 225x400 的矩阵,这是我的图像尺寸。我错过了什么吗?也许我应该替换一些东西,因为我有一个矩阵而不是一个数字序列(为其编写代码)?

4

1 回答 1

1

你可能想看看这个:

http://www.mathworks.com/matlabcentral/answers/2158-huffman-coding-and-decoding-for-image-jpeg-bmp

这似乎就在你的小巷里。它最终会引导你到这里:

http://www.mathworks.com/matlabcentral/fileexchange/26384-ppt-for-chapter-9-of-matlabsimulink-for-digital-communication

于 2012-06-28T13:41:49.087 回答