JPEG2000MATLAB代码

function y = im2jpeg2k（x n q）  
%IM2JPEG2K Compresses an image using a JPEG 2000 approximation. 
%   Y = IM2JPEG2K（X N Q） compresses image X using an N-scale JPEG 
%   2K wavelet transform implicit or explicit coefficient 
%   quantization and Huffman symbol coding augmented by zero 
%   run-length coding. If quantization vector Q contains two 
%   elements they are assumed to be implicit quantization 
%   parameters; else it is assumed to contain explicit subband step 
%   sizes.  Y is an encoding structure containing Huffman-encoded 
%   data and additional parameters needed by JPEG2K2IM for decoding. 
% 
%   See also JPEG2K2IM. 
 
%   Copyright 2002-2004 R. C. Gonzalez R. E. Woods & S. L. Eddins 
%   Digital Image Processing Using MATLAB Prentice-Hall 2004 
%   $Revision: 1.5 $  $Date: 2003/10/26 18:38:13 $ 
 
global RUNS 
 
error（nargchk（3 3 nargin））;         % Check input arguments 
 
if ndims（x） ~= 2 || ~isreal（x） || ~isnumeric（x） || ~isa（x ‘uint8‘） 
%   error（‘The input must be a UINT8 image.‘）;    
end 
 
if length（q） ~= 2 & length（q） ~= 3 * n + 1 
   error（‘The quantization step size vector is bad.‘）;   
end 
 
% Level shift the input and compute its wavelet transform. 
x = double（x） - 128; 
[c s] = wavefast（x n ‘jpeg9.7‘）; 
 
% Quantize the wavelet coefficients. 
q = stepsize（n q）; 
sgn = sign（c）;     sgn（find（sgn == 0）） = 1;     c = abs（c）; 
for k = 1:n 
   qi = 3 * k - 2; 
   c = wavepaste（‘h‘ c s k wavecopy（‘h‘ c s k） / q（qi））; 
   c = wavepaste（‘v‘ c s k wavecopy（‘v‘ c s k） / q（qi + 1））; 
   c = wavepaste（‘d‘ c s k wavecopy（‘d‘ c s k） / q（qi + 2））; 
end 
c = wavepaste（‘a‘ c s k wavecopy（‘a‘ c s k） / q（qi + 3））; 
c = floor（c）;       c = c .* sgn; 
 
% Run-length code zero runs of more than 10. Begin by creating 
% a special code for 0 runs （‘zrc‘） and end-of-code （‘eoc‘） and 
% making a run-length table.  
zrc = min（c（:）） - 1;      eoc = zrc - 1;      RUNS = [65535]; 
 
% Find the run transition points: ‘plus‘ contains the index of the 
% start of a