Gregor Weiss
008ccfa4a9
Cleaning up by removal of redundant information like parameter linked list, file extension, and file pointer including associated functionality. Start constructing an API taking void pointers to the original data for compression.
169 lines
No EOL
14 KiB
C
Executable file
169 lines
No EOL
14 KiB
C
Executable file
/*================================================================================================*\
|
|
|| ||
|
|
|| /$$$$$$$ /$$ /$$ /$$ /$$ ||
|
|
|| | $$__ $$|__/ | $$ /$ | $$| $$ ||
|
|
|| | $$ \ $$ /$$ /$$$$$$ | $$ /$$$| $$| $$$$$$$ /$$$$$$ /$$$$$$ /$$$$$$ ||
|
|
|| | $$$$$$$ | $$ /$$__ $$ | $$/$$ $$ $$| $$__ $$ /$$__ $$ /$$__ $$ /$$__ $$ ||
|
|
|| | $$__ $$| $$| $$ \ $$ | $$$$_ $$$$| $$ \ $$| $$ \ $$| $$ \ $$| $$ \ $$ ||
|
|
|| | $$ \ $$| $$| $$ | $$ | $$$/ \ $$$| $$ | $$| $$ | $$| $$ | $$| $$ | $$ ||
|
|
|| | $$$$$$$/| $$| $$$$$$$ | $$/ \ $$| $$ | $$| $$$$$$/| $$$$$$/| $$$$$$$/ ||
|
|
|| |_______/ |__/ \____ $$ |__/ \__/|__/ |__/ \______/ \______/ | $$____/ ||
|
|
|| /$$ \ $$ | $$ ||
|
|
|| | $$$$$$/ | $$ ||
|
|
|| \______/ |__/ ||
|
|
|| ||
|
|
|| DESCRIPTION: ||
|
|
|| ------------ ||
|
|
|| ||
|
|
|| This file describes a set of functions that can be used to performe the forward/ ||
|
|
|| inverse discrete wavelet transform on 1- to 4-dimensional IEEE 754 data-sets. ||
|
|
|| For more information please refere to JPEG2000 by D. S. Taubman and M. W. ||
|
|
|| Marcellin. ||
|
|
|| ||
|
|
|| -------------------------------------------------------------------------------------------- ||
|
|
|| Copyright (c) 2023, High Performance Computing Center - University of Stuttgart ||
|
|
|| ||
|
|
|| Redistribution and use in source and binary forms, with or without modification, are ||
|
|
|| permitted provided that the following conditions are met: ||
|
|
|| ||
|
|
|| (1) Redistributions of source code must retain the above copyright notice, this list of ||
|
|
|| conditions and the following disclaimer. ||
|
|
|| ||
|
|
|| (2) Redistributions in binary form must reproduce the above copyright notice, this list ||
|
|
|| of conditions and the following disclaimer in the documentation and/or other ||
|
|
|| materials provided with the distribution. ||
|
|
|| ||
|
|
|| THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS “AS IS” AND ANY EXPRESS ||
|
|
|| OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF ||
|
|
|| MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE ||
|
|
|| COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, ||
|
|
|| EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF ||
|
|
|| SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) ||
|
|
|| HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR ||
|
|
|| TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, ||
|
|
|| EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. ||
|
|
|| ||
|
|
\*================================================================================================*/
|
|
#ifndef DWT_H
|
|
#define DWT_H
|
|
/************************************************************************************************\
|
|
|| _ _ _ ____ _ _ _ ___ ____ ||
|
|
|| | |\ | | | | | | \ |___ ||
|
|
|| | | \| |___ |___ |__| |__/ |___ ||
|
|
|| ||
|
|
\************************************************************************************************/
|
|
#include "types.h"
|
|
|
|
/************************************************************************************************\
|
|
|| _ _ ____ ____ ____ ____ ____ ||
|
|
|| |\/| |__| | |__/ | | [__ ||
|
|
|| | | | | |___ | \ |__| ___] ||
|
|
|| ||
|
|
\************************************************************************************************/
|
|
/*----------------------------------------------------------------------------------------------*\
|
|
! !
|
|
! DESCRIPTION: !
|
|
! ------------ !
|
|
! !
|
|
! Maximum number of wavelet layers for which the energy gain is calculated. !
|
|
! !
|
|
\*----------------------------------------------------------------------------------------------*/
|
|
#define MAX_DECOMP_LEVELS 4
|
|
|
|
/*----------------------------------------------------------------------------------------------*\
|
|
! !
|
|
! DESCRIPTION: !
|
|
! ------------ !
|
|
! !
|
|
! These macros define the irrational coefficients for the high and low pass !
|
|
! synthesis filters associated with the (9-7) Cohen-Daubechies-Feauveau-Wavelet !
|
|
! as well as the coefficients for its coressponding lifting scheme. !
|
|
! !
|
|
\*----------------------------------------------------------------------------------------------*/
|
|
#define DWT_9X7_H0 1.115087052457000f // Low pass synthesis filter ...
|
|
#define DWT_9X7_H1 0.591271763114250f // ... coefficients
|
|
#define DWT_9X7_H2 -0.057543526228500f
|
|
#define DWT_9X7_H3 -0.091271763114250f
|
|
|
|
#define DWT_9X7_G0 1.205898036472721f // High pass synthesis filter ...
|
|
#define DWT_9X7_G1 -0.533728236885750f // ... coefficients
|
|
#define DWT_9X7_G2 -0.156446533057980f
|
|
#define DWT_9X7_G3 0.033728236885750f
|
|
#define DWT_9X7_G4 0.053497514821620f
|
|
|
|
#define ALPHA -1.586134342059924f // Lifting coefficients.
|
|
#define BETA -0.052980118572961f
|
|
#define GAMMA 0.882911075530934f
|
|
#define DELTA 0.360523644801462f
|
|
#define KAPPA_H 1.230174104914001f
|
|
#define KAPPA_L 0.812893066115961f
|
|
|
|
/*----------------------------------------------------------------------------------------------*\
|
|
! !
|
|
! DESCRIPTION: !
|
|
! ------------ !
|
|
! !
|
|
! These macros define the irrational coefficients for the high and low pass !
|
|
! synthesis filters associated with the (5-3) LeGall-Wavelet. !
|
|
! !
|
|
\*----------------------------------------------------------------------------------------------*/
|
|
#define DWT_5X3_H0 1.0f // Low pass synthesis filter ...
|
|
#define DWT_5X3_H1 0.5f // ... coefficients
|
|
|
|
#define DWT_5X3_G0 0.75f // High pass synthesis filter
|
|
#define DWT_5X3_G1 -0.25f // ... coefficients
|
|
#define DWT_5X3_G2 -0.125f
|
|
|
|
/*----------------------------------------------------------------------------------------------*\
|
|
! !
|
|
! DESCRIPTION: !
|
|
! ------------ !
|
|
! !
|
|
! These macros define the irrational coefficients for the high and low pass !
|
|
! synthesis filters associated with the Haar-Wavelet. !
|
|
! !
|
|
\*----------------------------------------------------------------------------------------------*/
|
|
#define DWT_HAAR_H0 1 // Low pass synthesis filter ...
|
|
#define DWT_HAAR_H1 1 // ... coefficients
|
|
|
|
#define DWT_HAAR_G0 0.5 // High pass synthesis filter ...
|
|
#define DWT_HAAR_G1 -0.5 // ... coefficients
|
|
|
|
/************************************************************************************************\
|
|
|| ____ _ _ ___ ____ ____ _ _ ____ _ _ _ ____ ____ _ ____ ___ _ ____ ____ ||
|
|
|| |___ \/ | |___ |__/ |\ | |__| | | | |__| |__/ | |__| |__] | |___ [__ ||
|
|
|| |___ _/\_ | |___ | \ | \| | | |___ \/ | | | \ | | | |__] |___ |___ ___] ||
|
|
|| ||
|
|
\************************************************************************************************/
|
|
/*----------------------------------------------------------------------------------------------*\
|
|
! !
|
|
! DESCRIPTION: !
|
|
! ------------ !
|
|
! !
|
|
! This array defines a lookup table used to store the energy gain factors for the !
|
|
! one dimensional CDF(9/7), LeGall-(5/3) and Haar wavelet transform. Each LUT !
|
|
! contains energy gain factors for the 6 low-pass and 6 high-pass decomposition !
|
|
! including level zero. !
|
|
! !
|
|
\*----------------------------------------------------------------------------------------------*/
|
|
extern double DWT_ENERGY_GAIN_LUT[3][2 * MAX_DECOMP_LEVELS + 2];
|
|
|
|
|
|
/************************************************************************************************\
|
|
|| ___ _ _ ___ _ _ ____ ____ _ _ _ _ ____ ___ _ ____ _ _ ____ ||
|
|
|| |__] | | |__] | | | |___ | | |\ | | | | | | |\ | [__ ||
|
|
|| | |__| |__] |___ | |___ | |__| | \| |___ | | |__| | \| ___] ||
|
|
|| ||
|
|
\************************************************************************************************/
|
|
uchar initialize_gain_lut ();
|
|
//==========|==========================|======================|======|=======|====================
|
|
bwc_float get_dwt_energy_gain (bwc_codec *const field,
|
|
uchar const highband_flag,
|
|
uint16 const level);
|
|
//==========|==========================|======================|======|=======|====================
|
|
uchar forward_wavelet_transform (bwc_codec *const field,
|
|
bwc_parameter *const parameter);
|
|
//==========|==========================|======================|======|=======|====================
|
|
uchar inverse_wavelet_transform (bwc_codec *const field,
|
|
bwc_parameter *const parameter);
|
|
#endif |