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 } 289
