git.maemo.org Git - opencv/blob - 3rdparty/libjasper/jpc_tsfb.c

   1 /*
   2  * Copyright (c) 1999-2000 Image Power, Inc. and the University of
   3  *   British Columbia.
   4  * Copyright (c) 2001-2004 Michael David Adams.
   5  * All rights reserved.
   6  */
   7
   8 /* __START_OF_JASPER_LICENSE__
   9  *
  10  * JasPer License Version 2.0
  11  *
  12  * Copyright (c) 2001-2006 Michael David Adams
  13  * Copyright (c) 1999-2000 Image Power, Inc.
  14  * Copyright (c) 1999-2000 The University of British Columbia
  15  *
  16  * All rights reserved.
  17  *
  18  * Permission is hereby granted, free of charge, to any person (the
  19  * "User") obtaining a copy of this software and associated documentation
  20  * files (the "Software"), to deal in the Software without restriction,
  21  * including without limitation the rights to use, copy, modify, merge,
  22  * publish, distribute, and/or sell copies of the Software, and to permit
  23  * persons to whom the Software is furnished to do so, subject to the
  24  * following conditions:
  25  *
  26  * 1.  The above copyright notices and this permission notice (which
  27  * includes the disclaimer below) shall be included in all copies or
  28  * substantial portions of the Software.
  29  *
  30  * 2.  The name of a copyright holder shall not be used to endorse or
  31  * promote products derived from the Software without specific prior
  32  * written permission.
  33  *
  34  * THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS
  35  * LICENSE.  NO USE OF THE SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER
  36  * THIS DISCLAIMER.  THE SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS
  37  * "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING
  38  * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
  39  * PARTICULAR PURPOSE AND NONINFRINGEMENT OF THIRD PARTY RIGHTS.  IN NO
  40  * EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, OR ANY SPECIAL
  41  * INDIRECT OR CONSEQUENTIAL DAMAGES, OR ANY DAMAGES WHATSOEVER RESULTING
  42  * FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT,
  43  * NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
  44  * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.  NO ASSURANCES ARE
  45  * PROVIDED BY THE COPYRIGHT HOLDERS THAT THE SOFTWARE DOES NOT INFRINGE
  46  * THE PATENT OR OTHER INTELLECTUAL PROPERTY RIGHTS OF ANY OTHER ENTITY.
  47  * EACH COPYRIGHT HOLDER DISCLAIMS ANY LIABILITY TO THE USER FOR CLAIMS
  48  * BROUGHT BY ANY OTHER ENTITY BASED ON INFRINGEMENT OF INTELLECTUAL
  49  * PROPERTY RIGHTS OR OTHERWISE.  AS A CONDITION TO EXERCISING THE RIGHTS
  50  * GRANTED HEREUNDER, EACH USER HEREBY ASSUMES SOLE RESPONSIBILITY TO SECURE
  51  * ANY OTHER INTELLECTUAL PROPERTY RIGHTS NEEDED, IF ANY.  THE SOFTWARE
  52  * IS NOT FAULT-TOLERANT AND IS NOT INTENDED FOR USE IN MISSION-CRITICAL
  53  * SYSTEMS, SUCH AS THOSE USED IN THE OPERATION OF NUCLEAR FACILITIES,
  54  * AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL
  55  * SYSTEMS, DIRECT LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH
  56  * THE FAILURE OF THE SOFTWARE OR SYSTEM COULD LEAD DIRECTLY TO DEATH,
  57  * PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE ("HIGH
  58  * RISK ACTIVITIES").  THE COPYRIGHT HOLDERS SPECIFICALLY DISCLAIM ANY
  59  * EXPRESS OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES.
  60  *
  61  * __END_OF_JASPER_LICENSE__
  62  */
  63
  64 /*
  65  * Tree-Structured Filter Bank (TSFB) Library
  66  *
  67  * $Id: jpc_tsfb.c,v 1.2 2008-05-26 09:40:52 vp153 Exp $
  68  */
  69
  70 /******************************************************************************\
  71 * Includes.
  72 \******************************************************************************/
  73
  74 #include <assert.h>
  75
  76 #include "jasper/jas_malloc.h"
  77 #include "jasper/jas_seq.h"
  78
  79 #include "jpc_tsfb.h"
  80 #include "jpc_cod.h"
  81 #include "jpc_cs.h"
  82 #include "jpc_util.h"
  83 #include "jpc_math.h"
  84
  85 void jpc_tsfb_getbands2(jpc_tsfb_t *tsfb, int locxstart, int locystart,
  86   int xstart, int ystart, int xend, int yend, jpc_tsfb_band_t **bands,
  87   int numlvls);
  88
  89 /******************************************************************************\
  90 *
  91 \******************************************************************************/
  92
  93 jpc_tsfb_t *jpc_cod_gettsfb(int qmfbid, int numlvls)
  94 {
  95         jpc_tsfb_t *tsfb;
  96
  97         if (!(tsfb = malloc(sizeof(jpc_tsfb_t))))
  98                 return 0;
  99
 100         if (numlvls > 0) {
 101                 switch (qmfbid) {
 102                 case JPC_COX_INS:
 103                         tsfb->qmfb = &jpc_ns_qmfb2d;
 104                         break;
 105                 default:
 106                 case JPC_COX_RFT:
 107                         tsfb->qmfb = &jpc_ft_qmfb2d;
 108                         break;
 109                 }
 110         } else {
 111                 tsfb->qmfb = 0;
 112         }
 113         tsfb->numlvls = numlvls;
 114         return tsfb;
 115 }
 116
 117 void jpc_tsfb_destroy(jpc_tsfb_t *tsfb)
 118 {
 119         free(tsfb);
 120 }
 121
 122 int jpc_tsfb_analyze(jpc_tsfb_t *tsfb, jas_seq2d_t *a)
 123 {
 124         return (tsfb->numlvls > 0) ? jpc_tsfb_analyze2(tsfb, jas_seq2d_getref(a,
 125           jas_seq2d_xstart(a), jas_seq2d_ystart(a)), jas_seq2d_xstart(a),
 126           jas_seq2d_ystart(a), jas_seq2d_width(a),
 127           jas_seq2d_height(a), jas_seq2d_rowstep(a), tsfb->numlvls - 1) : 0;
 128 }
 129
 130 int jpc_tsfb_analyze2(jpc_tsfb_t *tsfb, int *a, int xstart, int ystart,
 131   int width, int height, int stride, int numlvls)
 132 {
 133         if (width > 0 && height > 0) {
 134                 if ((*tsfb->qmfb->analyze)(a, xstart, ystart, width, height, stride))
 135                         return -1;
 136                 if (numlvls > 0) {
 137                         if (jpc_tsfb_analyze2(tsfb, a, JPC_CEILDIVPOW2(xstart,
 138                           1), JPC_CEILDIVPOW2(ystart, 1), JPC_CEILDIVPOW2(
 139                           xstart + width, 1) - JPC_CEILDIVPOW2(xstart, 1),
 140                           JPC_CEILDIVPOW2(ystart + height, 1) -
 141                           JPC_CEILDIVPOW2(ystart, 1), stride, numlvls - 1)) {
 142                                 return -1;
 143                         }
 144                 }
 145         }
 146         return 0;
 147 }
 148
 149 int jpc_tsfb_synthesize(jpc_tsfb_t *tsfb, jas_seq2d_t *a)
 150 {
 151         return (tsfb->numlvls > 0) ? jpc_tsfb_synthesize2(tsfb,
 152           jas_seq2d_getref(a, jas_seq2d_xstart(a), jas_seq2d_ystart(a)),
 153           jas_seq2d_xstart(a), jas_seq2d_ystart(a), jas_seq2d_width(a),
 154           jas_seq2d_height(a), jas_seq2d_rowstep(a), tsfb->numlvls - 1) : 0;
 155 }
 156
 157 int jpc_tsfb_synthesize2(jpc_tsfb_t *tsfb, int *a, int xstart, int ystart,
 158   int width, int height, int stride, int numlvls)
 159 {
 160         if (numlvls > 0) {
 161                 if (jpc_tsfb_synthesize2(tsfb, a, JPC_CEILDIVPOW2(xstart, 1),
 162                   JPC_CEILDIVPOW2(ystart, 1), JPC_CEILDIVPOW2(xstart + width,
 163                   1) - JPC_CEILDIVPOW2(xstart, 1), JPC_CEILDIVPOW2(ystart +
 164                   height, 1) - JPC_CEILDIVPOW2(ystart, 1), stride, numlvls -
 165                   1)) {
 166                         return -1;
 167                 }
 168         }
 169         if (width > 0 && height > 0) {
 170                 if ((*tsfb->qmfb->synthesize)(a, xstart, ystart, width, height, stride)) {
 171                         return -1;
 172                 }
 173         }
 174         return 0;
 175 }
 176
 177 int jpc_tsfb_getbands(jpc_tsfb_t *tsfb, uint_fast32_t xstart,
 178   uint_fast32_t ystart, uint_fast32_t xend, uint_fast32_t yend,
 179   jpc_tsfb_band_t *bands)
 180 {
 181         jpc_tsfb_band_t *band;
 182
 183         band = bands;
 184         if (tsfb->numlvls > 0) {
 185                 jpc_tsfb_getbands2(tsfb, xstart, ystart, xstart, ystart, xend, yend,
 186                   &band, tsfb->numlvls);
 187         } else {
 188
 189                 band->xstart = xstart;
 190                 band->ystart = ystart;
 191                 band->xend = xend;
 192                 band->yend = yend;
 193                 band->locxstart = xstart;
 194                 band->locystart = ystart;
 195                 band->locxend = band->locxstart + band->xend - band->xstart;
 196                 band->locyend = band->locystart + band->yend - band->ystart;
 197                 band->orient = JPC_TSFB_LL;
 198                 band->synenergywt = JPC_FIX_ONE;
 199                 ++band;
 200         }
 201         return band - bands;
 202 }
 203
 204 void jpc_tsfb_getbands2(jpc_tsfb_t *tsfb, int locxstart, int locystart,
 205   int xstart, int ystart, int xend, int yend, jpc_tsfb_band_t **bands,
 206   int numlvls)
 207 {
 208         int newxstart;
 209         int newystart;
 210         int newxend;
 211         int newyend;
 212         jpc_tsfb_band_t *band;
 213
 214         newxstart = JPC_CEILDIVPOW2(xstart, 1);
 215         newystart = JPC_CEILDIVPOW2(ystart, 1);
 216         newxend = JPC_CEILDIVPOW2(xend, 1);
 217         newyend = JPC_CEILDIVPOW2(yend, 1);
 218
 219         if (numlvls > 0) {
 220
 221                 jpc_tsfb_getbands2(tsfb, locxstart, locystart, newxstart, newystart,
 222                   newxend, newyend, bands, numlvls - 1);
 223
 224                 band = *bands;
 225                 band->xstart = JPC_FLOORDIVPOW2(xstart, 1);
 226                 band->ystart = newystart;
 227                 band->xend = JPC_FLOORDIVPOW2(xend, 1);
 228                 band->yend = newyend;
 229                 band->locxstart = locxstart + newxend - newxstart;
 230                 band->locystart = locystart;
 231                 band->locxend = band->locxstart + band->xend - band->xstart;
 232                 band->locyend = band->locystart + band->yend - band->ystart;
 233                 band->orient = JPC_TSFB_HL;
 234                 band->synenergywt = jpc_dbltofix(tsfb->qmfb->hpenergywts[
 235                   tsfb->numlvls - numlvls] * tsfb->qmfb->lpenergywts[
 236                   tsfb->numlvls - numlvls]);
 237                 ++(*bands);
 238
 239                 band = *bands;
 240                 band->xstart = newxstart;
 241                 band->ystart = JPC_FLOORDIVPOW2(ystart, 1);
 242                 band->xend = newxend;
 243                 band->yend = JPC_FLOORDIVPOW2(yend, 1);
 244                 band->locxstart = locxstart;
 245                 band->locystart = locystart + newyend - newystart;
 246                 band->locxend = band->locxstart + band->xend - band->xstart;
 247                 band->locyend = band->locystart + band->yend - band->ystart;
 248                 band->orient = JPC_TSFB_LH;
 249                 band->synenergywt = jpc_dbltofix(tsfb->qmfb->lpenergywts[
 250                   tsfb->numlvls - numlvls] * tsfb->qmfb->hpenergywts[
 251                   tsfb->numlvls - numlvls]);
 252                 ++(*bands);
 253
 254                 band = *bands;
 255                 band->xstart = JPC_FLOORDIVPOW2(xstart, 1);
 256                 band->ystart = JPC_FLOORDIVPOW2(ystart, 1);
 257                 band->xend = JPC_FLOORDIVPOW2(xend, 1);
 258                 band->yend = JPC_FLOORDIVPOW2(yend, 1);
 259                 band->locxstart = locxstart + newxend - newxstart;
 260                 band->locystart = locystart + newyend - newystart;
 261                 band->locxend = band->locxstart + band->xend - band->xstart;
 262                 band->locyend = band->locystart + band->yend - band->ystart;
 263                 band->orient = JPC_TSFB_HH;
 264                 band->synenergywt = jpc_dbltofix(tsfb->qmfb->hpenergywts[
 265                   tsfb->numlvls - numlvls] * tsfb->qmfb->hpenergywts[
 266                   tsfb->numlvls - numlvls]);
 267                 ++(*bands);
 268
 269         } else {
 270
 271                 band = *bands;
 272                 band->xstart = xstart;
 273                 band->ystart = ystart;
 274                 band->xend = xend;
 275                 band->yend = yend;
 276                 band->locxstart = locxstart;
 277                 band->locystart = locystart;
 278                 band->locxend = band->locxstart + band->xend - band->xstart;
 279                 band->locyend = band->locystart + band->yend - band->ystart;
 280                 band->orient = JPC_TSFB_LL;
 281                 band->synenergywt = jpc_dbltofix(tsfb->qmfb->lpenergywts[
 282                   tsfb->numlvls - numlvls - 1] * tsfb->qmfb->lpenergywts[
 283                   tsfb->numlvls - numlvls - 1]);
 284                 ++(*bands);
 285
 286         }
 287
 288 }