DWT - Implementation of a discrete wavelet transform

Data Structures
struct	dwt_local
union	opj_v4_t
struct	v4dwt_local

Macros
#define	OPJ_WS(i)   v->mem[(i)*2]
#define	OPJ_WD(i)   v->mem[(1+(i)*2)]
#define	VREG_INT_COUNT   4
	Number of int32 values in a SSE2 register. More...
#define	PARALLEL_COLS_53   (2*VREG_INT_COUNT)
	Number of columns that we can process in parallel in the vertical pass. More...

Typedefs
typedef void(*	DWT1DFN) (const opj_dwt_t *v)
	Virtual function type for wavelet transform in 1-D. More...

Local data structures
typedef struct dwt_local	opj_dwt_t
typedef struct v4dwt_local	opj_v4dwt_t
static const OPJ_FLOAT32	opj_dwt_alpha = 1.586134342f
static const OPJ_FLOAT32	opj_dwt_beta = 0.052980118f
static const OPJ_FLOAT32	opj_dwt_gamma = -0.882911075f
static const OPJ_FLOAT32	opj_dwt_delta = -0.443506852f
static const OPJ_FLOAT32	opj_K = 1.230174105f
static const OPJ_FLOAT32	opj_c13318 = 1.625732422f

Local static functions
static void	opj_dwt_deinterleave_h (OPJ_INT32 *a, OPJ_INT32 *b, OPJ_INT32 dn, OPJ_INT32 sn, OPJ_INT32 cas)
	Forward lazy transform (horizontal) More...
static void	opj_dwt_deinterleave_v (OPJ_INT32 *a, OPJ_INT32 *b, OPJ_INT32 dn, OPJ_INT32 sn, OPJ_INT32 x, OPJ_INT32 cas)
	Forward lazy transform (vertical) More...
static void	opj_dwt_encode_1 (OPJ_INT32 *a, OPJ_INT32 dn, OPJ_INT32 sn, OPJ_INT32 cas)
	Forward 5-3 wavelet transform in 1-D. More...
static void	opj_dwt_encode_1_real (OPJ_INT32 *a, OPJ_INT32 dn, OPJ_INT32 sn, OPJ_INT32 cas)
	Forward 9-7 wavelet transform in 1-D. More...
static void	opj_dwt_encode_stepsize (OPJ_INT32 stepsize, OPJ_INT32 numbps, opj_stepsize_t *bandno_stepsize)
	Explicit calculation of the Quantization Stepsizes. More...
static OPJ_BOOL	opj_dwt_decode_tile (opj_thread_pool_t *tp, opj_tcd_tilecomp_t *tilec, OPJ_UINT32 i)
	Inverse wavelet transform in 2-D. More...
static OPJ_BOOL	opj_dwt_decode_partial_tile (opj_tcd_tilecomp_t *tilec, OPJ_UINT32 numres)
static OPJ_BOOL	opj_dwt_encode_procedure (opj_tcd_tilecomp_t *tilec, void(*p_function)(OPJ_INT32 *, OPJ_INT32, OPJ_INT32, OPJ_INT32))
static OPJ_UINT32	opj_dwt_max_resolution (opj_tcd_resolution_t *OPJ_RESTRICT r, OPJ_UINT32 i)
static void	opj_v4dwt_decode (opj_v4dwt_t *OPJ_RESTRICT dwt)
static void	opj_v4dwt_interleave_h (opj_v4dwt_t *OPJ_RESTRICT dwt, OPJ_FLOAT32 *OPJ_RESTRICT a, OPJ_UINT32 width, OPJ_UINT32 remaining_height)
static void	opj_v4dwt_interleave_v (opj_v4dwt_t *OPJ_RESTRICT dwt, OPJ_FLOAT32 *OPJ_RESTRICT a, OPJ_UINT32 width, OPJ_UINT32 nb_elts_read)
static void	opj_v4dwt_decode_step1 (opj_v4_t *w, OPJ_UINT32 start, OPJ_UINT32 end, const OPJ_FLOAT32 c)
static void	opj_v4dwt_decode_step2 (opj_v4_t *l, opj_v4_t *w, OPJ_UINT32 start, OPJ_UINT32 end, OPJ_UINT32 m, OPJ_FLOAT32 c)

Exported functions
OPJ_BOOL	opj_dwt_encode (opj_tcd_tilecomp_t *tilec)
	Forward 5-3 wavelet transform in 2-D. More...
OPJ_BOOL	opj_dwt_decode (opj_tcd_t *p_tcd, opj_tcd_tilecomp_t *tilec, OPJ_UINT32 numres)
	Inverse 5-3 wavelet transform in 2-D. More...
OPJ_UINT32	opj_dwt_getgain (OPJ_UINT32 orient)
	Get the gain of a subband for the reversible 5-3 DWT. More...
OPJ_FLOAT64	opj_dwt_getnorm (OPJ_UINT32 level, OPJ_UINT32 orient)
	Get the norm of a wavelet function of a subband at a specified level for the reversible 5-3 DWT. More...
OPJ_BOOL	opj_dwt_encode_real (opj_tcd_tilecomp_t *tilec)
	Forward 9-7 wavelet transform in 2-D. More...
OPJ_BOOL	opj_dwt_decode_real (opj_tcd_t *p_tcd, opj_tcd_tilecomp_t *OPJ_RESTRICT tilec, OPJ_UINT32 numres)
	Inverse 9-7 wavelet transform in 2-D. More...
OPJ_UINT32	opj_dwt_getgain_real (OPJ_UINT32 orient)
	Get the gain of a subband for the irreversible 9-7 DWT. More...
OPJ_FLOAT64	opj_dwt_getnorm_real (OPJ_UINT32 level, OPJ_UINT32 orient)
	Get the norm of a wavelet function of a subband at a specified level for the irreversible 9-7 DWT. More...
void	opj_dwt_calc_explicit_stepsizes (opj_tccp_t *tccp, OPJ_UINT32 prec)
	Explicit calculation of the Quantization Stepsizes. More...

Detailed Description

Macro Definition Documentation

#define OPJ_WD

(

i

)

v->mem[(1+(i)*2)]

#define OPJ_WS

(

i

)

v->mem[(i)*2]

#define PARALLEL_COLS_53   (2*VREG_INT_COUNT)

Number of columns that we can process in parallel in the vertical pass.

Referenced by opj_dwt_decode_tile(), opj_dwt_decode_v_func(), and opj_idwt53_v().

#define VREG_INT_COUNT   4

Number of int32 values in a SSE2 register.

Typedef Documentation

typedef void(* DWT1DFN) (const opj_dwt_t *v)

Virtual function type for wavelet transform in 1-D.

typedef struct dwt_local opj_dwt_t

typedef struct v4dwt_local opj_v4dwt_t

Function Documentation

void opj_dwt_calc_explicit_stepsizes	(	opj_tccp_t *	tccp,
		OPJ_UINT32	prec
	)

Explicit calculation of the Quantization Stepsizes.

Parameters

tccp	Tile-component coding parameters
prec	Precint analyzed

References J2K_CCP_QNTSTY_NOQNT, opj_tccp::numresolutions, opj_dwt_encode_stepsize(), opj_tccp::qmfbid, opj_tccp::qntsty, and opj_tccp::stepsizes.

Referenced by opj_j2k_setup_encoder().

OPJ_BOOL opj_dwt_decode	(	opj_tcd_t *	p_tcd,
		opj_tcd_tilecomp_t *	tilec,
		OPJ_UINT32	numres
	)

Inverse 5-3 wavelet transform in 2-D.

Apply a reversible inverse DWT transform to a component of an image.

Parameters

p_tcd	TCD handle
tilec	Tile component information (current tile)
numres	Number of resolution levels to decode

References opj_dwt_decode_partial_tile(), opj_dwt_decode_tile(), opj_tcd::thread_pool, and opj_tcd::whole_tile_decoding.

Referenced by main(), and opj_tcd_dwt_decode().

static OPJ_BOOL opj_dwt_decode_partial_tile ( opj_tcd_tilecomp_t *  tilec, OPJ_UINT32  numres  )

static

References opj_tcd_resolution::bands, dwt_local::cas, opj_tcd_tilecomp::data_win, dwt_local::dn, dwt_local::mem, opj_aligned_32_malloc(), opj_aligned_free(), opj_dwt_decode_partial_1(), opj_dwt_decode_partial_1_parallel(), opj_dwt_get_band_coordinates(), opj_dwt_init_sparse_array(), opj_dwt_interleave_partial_h(), opj_dwt_interleave_partial_v(), opj_dwt_max_resolution(), opj_dwt_segment_grow(), OPJ_FALSE, opj_sparse_array_int32_free(), opj_sparse_array_int32_read(), opj_sparse_array_int32_write(), OPJ_TRUE, opj_uint_max(), opj_uint_min(), opj_uint_subs(), OPJ_UNUSED, opj_tcd_tilecomp::resolutions, SIZE_MAX, dwt_local::sn, opj_tcd_resolution::win_x0, opj_tcd_tilecomp::win_x0, opj_tcd_resolution::win_x1, opj_tcd_tilecomp::win_x1, opj_tcd_resolution::win_y0, opj_tcd_tilecomp::win_y0, opj_tcd_resolution::win_y1, opj_tcd_tilecomp::win_y1, opj_tcd_band::x0, opj_tcd_resolution::x0, opj_tcd_resolution::x1, opj_tcd_band::y0, opj_tcd_resolution::y0, and opj_tcd_resolution::y1.

Referenced by opj_dwt_decode().

OPJ_BOOL opj_dwt_decode_real	(	opj_tcd_t *	p_tcd,
		opj_tcd_tilecomp_t *OPJ_RESTRICT	tilec,
		OPJ_UINT32	numres
	)

Inverse 9-7 wavelet transform in 2-D.

Apply an irreversible inverse DWT transform to a component of an image.

Parameters

p_tcd	TCD handle
tilec	Tile component information (current tile)
numres	Number of resolution levels to decode

References opj_dwt_decode_partial_97(), opj_dwt_decode_tile_97(), and opj_tcd::whole_tile_decoding.

Referenced by opj_tcd_dwt_decode().

static OPJ_BOOL opj_dwt_decode_tile ( opj_thread_pool_t *  tp, opj_tcd_tilecomp_t *  tilec, OPJ_UINT32  i  )

static

Inverse wavelet transform in 2-D.

References dwt_local::cas, opj_tcd_tilecomp::data, dwt_local::dn, opj_dwd_decode_h_job_t::h, opj_dwd_decode_h_job_t::max_j, opj_dwd_decode_v_job_t::max_j, dwt_local::mem, opj_dwd_decode_h_job_t::min_j, opj_dwd_decode_v_job_t::min_j, opj_tcd_tilecomp::minimum_num_resolutions, opj_aligned_32_malloc(), opj_aligned_free(), opj_dwt_decode_h_func(), opj_dwt_decode_v_func(), opj_dwt_max_resolution(), OPJ_FALSE, opj_free(), opj_idwt53_h(), opj_idwt53_v(), opj_malloc(), OPJ_RESTRICT, opj_thread_pool_get_thread_count(), opj_thread_pool_submit_job(), opj_thread_pool_wait_completion(), OPJ_TRUE, PARALLEL_COLS_53, opj_tcd_tilecomp::resolutions, opj_dwd_decode_v_job_t::rh, opj_dwd_decode_h_job_t::rw, SIZE_MAX, dwt_local::sn, opj_dwd_decode_h_job_t::tiledp, opj_dwd_decode_v_job_t::tiledp, opj_dwd_decode_v_job_t::v, opj_dwd_decode_h_job_t::w, opj_dwd_decode_v_job_t::w, opj_tcd_resolution::x0, opj_tcd_resolution::x1, opj_tcd_resolution::y0, and opj_tcd_resolution::y1.

Referenced by opj_dwt_decode().

static void opj_dwt_deinterleave_h ( OPJ_INT32 *  a, OPJ_INT32 *  b, OPJ_INT32  dn, OPJ_INT32  sn, OPJ_INT32  cas  )

static

Forward lazy transform (horizontal)

Referenced by opj_dwt_encode_procedure().

static void opj_dwt_deinterleave_v ( OPJ_INT32 *  a, OPJ_INT32 *  b, OPJ_INT32  dn, OPJ_INT32  sn, OPJ_INT32  x, OPJ_INT32  cas  )

static

Forward lazy transform (vertical)

Referenced by opj_dwt_encode_procedure().

OPJ_BOOL opj_dwt_encode

(

opj_tcd_tilecomp_t *

tilec

)

Forward 5-3 wavelet transform in 2-D.

Apply a reversible DWT transform to a component of an image.

Parameters

tilec

Tile component information (current tile)

References opj_dwt_encode_1(), and opj_dwt_encode_procedure().

Referenced by main(), and opj_tcd_dwt_encode().

static void opj_dwt_encode_1 ( OPJ_INT32 *  a, OPJ_INT32  dn, OPJ_INT32  sn, OPJ_INT32  cas  )

static

Forward 5-3 wavelet transform in 1-D.

References OPJ_D, OPJ_D_, OPJ_DD_, OPJ_S, OPJ_S_, and OPJ_SS_.

Referenced by opj_dwt_encode().

static void opj_dwt_encode_1_real ( OPJ_INT32 *  a, OPJ_INT32  dn, OPJ_INT32  sn, OPJ_INT32  cas  )

static

Forward 9-7 wavelet transform in 1-D.

References OPJ_D, OPJ_D_, OPJ_DD_, opj_int_fix_mul(), OPJ_S, OPJ_S_, and OPJ_SS_.

Referenced by opj_dwt_encode_real().

static INLINE OPJ_BOOL opj_dwt_encode_procedure ( opj_tcd_tilecomp_t *  tilec, void(*)(OPJ_INT32 *, OPJ_INT32, OPJ_INT32, OPJ_INT32)  p_function  )

static

References opj_tcd_tilecomp::data, opj_tcd_tilecomp::numresolutions, opj_dwt_deinterleave_h(), opj_dwt_deinterleave_v(), opj_dwt_max_resolution(), OPJ_FALSE, opj_free(), opj_malloc(), OPJ_TRUE, opj_tcd_tilecomp::resolutions, SIZE_MAX, opj_tcd_resolution::x0, opj_tcd_tilecomp::x0, opj_tcd_resolution::x1, opj_tcd_tilecomp::x1, opj_tcd_resolution::y0, and opj_tcd_resolution::y1.

Referenced by opj_dwt_encode(), and opj_dwt_encode_real().

OPJ_BOOL opj_dwt_encode_real

(

opj_tcd_tilecomp_t *

tilec

)

Forward 9-7 wavelet transform in 2-D.

Apply an irreversible DWT transform to a component of an image.

Parameters

tilec

Tile component information (current tile)

References opj_dwt_encode_1_real(), and opj_dwt_encode_procedure().

Referenced by opj_tcd_dwt_encode().

static void opj_dwt_encode_stepsize ( OPJ_INT32  stepsize, OPJ_INT32  numbps, opj_stepsize_t *  bandno_stepsize  )

static

Explicit calculation of the Quantization Stepsizes.

References opj_stepsize::expn, opj_stepsize::mant, and opj_int_floorlog2().

Referenced by opj_dwt_calc_explicit_stepsizes().

OPJ_UINT32 opj_dwt_getgain

(

OPJ_UINT32

orient

)

Get the gain of a subband for the reversible 5-3 DWT.

Parameters

orient

Number that identifies the subband (0->LL, 1->HL, 2->LH, 3->HH)

Returns

Returns 0 if orient = 0, returns 1 if orient = 1 or 2, returns 2 otherwise

Referenced by opj_tcd_init_tile().

OPJ_UINT32 opj_dwt_getgain_real

(

OPJ_UINT32

orient

)

Get the gain of a subband for the irreversible 9-7 DWT.

Parameters

orient

Number that identifies the subband (0->LL, 1->HL, 2->LH, 3->HH)

Returns

Returns the gain of the 9-7 wavelet transform

Referenced by opj_tcd_init_tile().

OPJ_FLOAT64 opj_dwt_getnorm	(	OPJ_UINT32	level,
		OPJ_UINT32	orient
	)

Get the norm of a wavelet function of a subband at a specified level for the reversible 5-3 DWT.

Parameters

level	Level of the wavelet function
orient	Band of the wavelet function

Returns

Returns the norm of the wavelet function

Referenced by opj_t1_getwmsedec().

OPJ_FLOAT64 opj_dwt_getnorm_real	(	OPJ_UINT32	level,
		OPJ_UINT32	orient
	)

Get the norm of a wavelet function of a subband at a specified level for the irreversible 9-7 DWT.

Parameters

level	Level of the wavelet function
orient	Band of the wavelet function

Returns

Returns the norm of the 9-7 wavelet

Referenced by opj_t1_getwmsedec().

static OPJ_UINT32 opj_dwt_max_resolution ( opj_tcd_resolution_t *OPJ_RESTRICT  r, OPJ_UINT32  i  )

static

Referenced by opj_dwt_decode_partial_97(), opj_dwt_decode_partial_tile(), opj_dwt_decode_tile(), opj_dwt_decode_tile_97(), and opj_dwt_encode_procedure().

static void opj_v4dwt_decode ( opj_v4dwt_t *OPJ_RESTRICT  dwt )

static

References opj_c13318, opj_dwt_alpha, opj_dwt_beta, opj_dwt_delta, opj_dwt_gamma, opj_int_min(), opj_K, opj_v4dwt_decode_step1(), and opj_v4dwt_decode_step2().

Referenced by opj_dwt_decode_partial_97(), and opj_dwt_decode_tile_97().

static void opj_v4dwt_decode_step1 ( opj_v4_t *  w, OPJ_UINT32  start, OPJ_UINT32  end, const OPJ_FLOAT32  c  )

static

References OPJ_RESTRICT.

Referenced by opj_v4dwt_decode().

static void opj_v4dwt_decode_step2 ( opj_v4_t *  l, opj_v4_t *  w, OPJ_UINT32  start, OPJ_UINT32  end, OPJ_UINT32  m, OPJ_FLOAT32  c  )

static

References opj_uint_min().

Referenced by opj_v4dwt_decode().

static void opj_v4dwt_interleave_h ( opj_v4dwt_t *OPJ_RESTRICT  dwt, OPJ_FLOAT32 *OPJ_RESTRICT  a, OPJ_UINT32  width, OPJ_UINT32  remaining_height  )

static

References OPJ_RESTRICT.

Referenced by opj_dwt_decode_tile_97().

static void opj_v4dwt_interleave_v ( opj_v4dwt_t *OPJ_RESTRICT  dwt, OPJ_FLOAT32 *OPJ_RESTRICT  a, OPJ_UINT32  width, OPJ_UINT32  nb_elts_read  )

static

References OPJ_RESTRICT.

Referenced by opj_dwt_decode_tile_97().

Variable Documentation

const OPJ_FLOAT32 opj_c13318 = 1.625732422f

static

Referenced by opj_v4dwt_decode().

const OPJ_FLOAT32 opj_dwt_alpha = 1.586134342f

static

Referenced by opj_v4dwt_decode().

const OPJ_FLOAT32 opj_dwt_beta = 0.052980118f

static

Referenced by opj_v4dwt_decode().

const OPJ_FLOAT32 opj_dwt_delta = -0.443506852f

static

Referenced by opj_v4dwt_decode().

const OPJ_FLOAT32 opj_dwt_gamma = -0.882911075f

static

Referenced by opj_v4dwt_decode().

const OPJ_FLOAT32 opj_K = 1.230174105f

static

Referenced by opj_v4dwt_decode().