Irreversible compression/decompression DWT: use 1/K constant as per standard
[openjpeg.git] / src / lib / openjp2 / dwt.h
1 /*
2  * The copyright in this software is being made available under the 2-clauses
3  * BSD License, included below. This software may be subject to other third
4  * party and contributor rights, including patent rights, and no such rights
5  * are granted under this license.
6  *
7  * Copyright (c) 2002-2014, Universite catholique de Louvain (UCL), Belgium
8  * Copyright (c) 2002-2014, Professor Benoit Macq
9  * Copyright (c) 2001-2003, David Janssens
10  * Copyright (c) 2002-2003, Yannick Verschueren
11  * Copyright (c) 2003-2007, Francois-Olivier Devaux
12  * Copyright (c) 2003-2014, Antonin Descampe
13  * Copyright (c) 2005, Herve Drolon, FreeImage Team
14  * All rights reserved.
15  *
16  * Redistribution and use in source and binary forms, with or without
17  * modification, are permitted provided that the following conditions
18  * are met:
19  * 1. Redistributions of source code must retain the above copyright
20  *    notice, this list of conditions and the following disclaimer.
21  * 2. Redistributions in binary form must reproduce the above copyright
22  *    notice, this list of conditions and the following disclaimer in the
23  *    documentation and/or other materials provided with the distribution.
24  *
25  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
26  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
28  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
29  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
30  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
31  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
32  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
33  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
34  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
35  * POSSIBILITY OF SUCH DAMAGE.
36  */
37
38 #ifndef OPJ_DWT_H
39 #define OPJ_DWT_H
40 /**
41 @file dwt.h
42 @brief Implementation of a discrete wavelet transform (DWT)
43
44 The functions in DWT.C have for goal to realize forward and inverse discret wavelet
45 transform with filter 5-3 (reversible) and filter 9-7 (irreversible). The functions in
46 DWT.C are used by some function in TCD.C.
47 */
48
49 /** @defgroup DWT DWT - Implementation of a discrete wavelet transform */
50 /*@{*/
51
52
53 /** @name Exported functions */
54 /*@{*/
55 /* ----------------------------------------------------------------------- */
56 /**
57 Forward 5-3 wavelet transform in 2-D.
58 Apply a reversible DWT transform to a component of an image.
59 @param p_tcd TCD handle
60 @param tilec Tile component information (current tile)
61 */
62 OPJ_BOOL opj_dwt_encode(opj_tcd_t *p_tcd,
63                         opj_tcd_tilecomp_t * tilec);
64
65 /**
66 Inverse 5-3 wavelet transform in 2-D.
67 Apply a reversible inverse DWT transform to a component of an image.
68 @param p_tcd TCD handle
69 @param tilec Tile component information (current tile)
70 @param numres Number of resolution levels to decode
71 */
72 OPJ_BOOL opj_dwt_decode(opj_tcd_t *p_tcd,
73                         opj_tcd_tilecomp_t* tilec,
74                         OPJ_UINT32 numres);
75
76 /**
77 Get the norm of a wavelet function of a subband at a specified level for the reversible 5-3 DWT.
78 @param level Level of the wavelet function
79 @param orient Band of the wavelet function
80 @return Returns the norm of the wavelet function
81 */
82 OPJ_FLOAT64 opj_dwt_getnorm(OPJ_UINT32 level, OPJ_UINT32 orient);
83 /**
84 Forward 9-7 wavelet transform in 2-D.
85 Apply an irreversible DWT transform to a component of an image.
86 @param p_tcd TCD handle
87 @param tilec Tile component information (current tile)
88 */
89 OPJ_BOOL opj_dwt_encode_real(opj_tcd_t *p_tcd,
90                              opj_tcd_tilecomp_t * tilec);
91 /**
92 Inverse 9-7 wavelet transform in 2-D.
93 Apply an irreversible inverse DWT transform to a component of an image.
94 @param p_tcd TCD handle
95 @param tilec Tile component information (current tile)
96 @param numres Number of resolution levels to decode
97 */
98 OPJ_BOOL opj_dwt_decode_real(opj_tcd_t *p_tcd,
99                              opj_tcd_tilecomp_t* OPJ_RESTRICT tilec,
100                              OPJ_UINT32 numres);
101
102 /**
103 Get the norm of a wavelet function of a subband at a specified level for the irreversible 9-7 DWT
104 @param level Level of the wavelet function
105 @param orient Band of the wavelet function
106 @return Returns the norm of the 9-7 wavelet
107 */
108 OPJ_FLOAT64 opj_dwt_getnorm_real(OPJ_UINT32 level, OPJ_UINT32 orient);
109 /**
110 Explicit calculation of the Quantization Stepsizes
111 @param tccp Tile-component coding parameters
112 @param prec Precint analyzed
113 */
114 void opj_dwt_calc_explicit_stepsizes(opj_tccp_t * tccp, OPJ_UINT32 prec);
115 /* ----------------------------------------------------------------------- */
116 /*@}*/
117
118 /*@}*/
119
120 #endif /* OPJ_DWT_H */