小波变换C代码(wavelet transform code)

WaveCode copyright (c) 1998 by Charles Bloom v1.0 released May 23, 1998 ----------------------------------------------------------------------- Introduction version 1.0 does nothing revolutionary ; it's our initial implementation. The advantage of WaveCode over existing transform coders is its modularity. Choose from many transforms : Wavelet [Daubechies 7/9] (no flag) DCT (tree structures) (-d) S + Prediction (-s) And many coders : (use -c#) : Bitplane (simplified EZW) Significance Map Significance Map 2 Order 1 Order 0 Vector Quantizing Fractal -style VQ The Bitplane Coder with Wavelet Transform performance roughly identically to SPIHT (Said & Pearlman's Improved EZW). We also allow any choice of # of wavelet levels (-l# option), and with more levels (SPIHT uses l=3) we get even better performance. Our coders and transforms are quite competitive. Our quantizers are very poor in v1.0, offering as a best option uniform quantization. We also have no bit allocation handling. Perhaps the nicest thing in WaveCode are the (-w) write options, with which you can output a raw image of the wavelet (or DCT, or S+P) data, a map of the errors, as well as the standard reconstructed image. So far as I know this is the first time anyone has really *looked* at DCT data (in tree structured form where you can make sense of it visually) and let me tell you, it is a beautiful sight. When comparing DCT and Wavelet be sure to use (-l3) on the wavelets, because using more levels (4 is default) gives the wavelets an unfair advantage. ----------------------------------------------------------------------- Some Details All images input and output are raw byteplanes. Adobe Photoshop, for example, is convenient for viewing and converting raw images. The "coder_param" is used in the following way by the coder modules : VQ : sets maximum allowed MSE in a fuzzy match. other : ignored The "tune_param" is for my use only. :^) If anyone has efficient code for a 23/25 wavelet transform, I'd love to get my hands on it. ----------------------------------------------------------------------- Acknowledgements: The DCT code is stolen from the IJG (Indepedent JPEG Group) distribution. The Floating (DWT) Wavelet Transform is stolen from Robert Estes. All other code is original. Thanks to the many people on comp.compression who have helped me get into image coding. Among them are Tom Lane, Michael Schindler, Matt Timmermans, Jim Andrew, Peter Lindstrom, Andrew Watson, and many I'm sure I'm forgetting. ----------------------------------------------------------------------- Apologies: I've committed the cardinal sin of introducing yet another imperfect Wavelet coder into the mix. I know, I'm sorry, there are *no* practical wavelet coders (the way IJG JPEG is practical) and this code is not it either. Other nice code available on the net is Geoff Davis's package : we implement here some of the todos on GD's list, but our quantizers and (lack of) bit allocation makes it a trade-off ; perhaps my time would have been better spent adding his code... Robert Estes' code is also nice because of its simplicity and readability. All other wavelet code I've found is either baroque (LiftPack) or full of Fortran (EPIC/EPWIC) or not ANSI (wavelet_tutorial). Anyway, I apologize for contributing to this mix - if the trend holds up, soon everyone in America will be writing their own wavelet compressor, but they'll all be imperfect, and we'll still be using JPEG to send images !