PageRenderTime 460ms CodeModel.GetById 47ms RepoModel.GetById 10ms app.codeStats 1ms

/src/FreeImage/Source/LibOpenJPEG/tcd.c

https://bitbucket.org/cabalistic/ogredeps/
C | 1524 lines | 1151 code | 237 blank | 136 comment | 274 complexity | 9290d021d70d51991d16f4b8f83f4ab7 MD5 | raw file
Possible License(s): LGPL-3.0, BSD-3-Clause, CPL-1.0, Unlicense, GPL-2.0, GPL-3.0, LGPL-2.0, MPL-2.0-no-copyleft-exception, BSD-2-Clause, LGPL-2.1 
    

  1. /*
    2. 

  2.  * Copyright (c) 2002-2007, Communications and Remote Sensing Laboratory, Universite catholique de Louvain (UCL), Belgium
    3. 

  3.  * Copyright (c) 2002-2007, Professor Benoit Macq
    4. 

  4.  * Copyright (c) 2001-2003, David Janssens
    5. 

  5.  * Copyright (c) 2002-2003, Yannick Verschueren
    6. 

  6.  * Copyright (c) 2003-2007, Francois-Olivier Devaux and Antonin Descampe
    7. 

  7.  * Copyright (c) 2005, Herve Drolon, FreeImage Team
    8. 

  8.  * Copyright (c) 2006-2007, Parvatha Elangovan
    9. 

  9.  * All rights reserved.
    10. 

  10.  *
    11. 

  11.  * Redistribution and use in source and binary forms, with or without
    12. 

  12.  * modification, are permitted provided that the following conditions
    13. 

  13.  * are met:
    14. 

  14.  * 1. Redistributions of source code must retain the above copyright
    15. 

  15.  *    notice, this list of conditions and the following disclaimer.
    16. 

  16.  * 2. Redistributions in binary form must reproduce the above copyright
    17. 

  17.  *    notice, this list of conditions and the following disclaimer in the
    18. 

  18.  *    documentation and/or other materials provided with the distribution.
    19. 

  19.  *
    20. 

  20.  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
    21. 

  21.  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
    22. 

  22.  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
    23. 

  23.  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
    24. 

  24.  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
    25. 

  25.  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
    26. 

  26.  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
    27. 

  27.  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
    28. 

  28.  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
    29. 

  29.  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
    30. 

  30.  * POSSIBILITY OF SUCH DAMAGE.
    31. 

  31.  */
    32. 

  32. 

  33. #include "opj_includes.h"
    34. 

  34. 

  35. void tcd_dump(FILE *fd, opj_tcd_t *tcd, opj_tcd_image_t * img) {
    36. 

  36. 	int tileno, compno, resno, bandno, precno;/*, cblkno;*/
    37. 

  37. 

  38. 	fprintf(fd, "image {\n");
    39. 

  39. 	fprintf(fd, "  tw=%d, th=%d x0=%d x1=%d y0=%d y1=%d\n", 
    40. 

  40. 		img->tw, img->th, tcd->image->x0, tcd->image->x1, tcd->image->y0, tcd->image->y1);
    41. 

  41. 

  42. 	for (tileno = 0; tileno < img->th * img->tw; tileno++) {
    43. 

  43. 		opj_tcd_tile_t *tile = &tcd->tcd_image->tiles[tileno];
    44. 

  44. 		fprintf(fd, "  tile {\n");
    45. 

  45. 		fprintf(fd, "    x0=%d, y0=%d, x1=%d, y1=%d, numcomps=%d\n",
    46. 

  46. 			tile->x0, tile->y0, tile->x1, tile->y1, tile->numcomps);
    47. 

  47. 		for (compno = 0; compno < tile->numcomps; compno++) {
    48. 

  48. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    49. 

  49. 			fprintf(fd, "    tilec {\n");
    50. 

  50. 			fprintf(fd,
    51. 

  51. 				"      x0=%d, y0=%d, x1=%d, y1=%d, numresolutions=%d\n",
    52. 

  52. 				tilec->x0, tilec->y0, tilec->x1, tilec->y1, tilec->numresolutions);
    53. 

  53. 			for (resno = 0; resno < tilec->numresolutions; resno++) {
    54. 

  54. 				opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    55. 

  55. 				fprintf(fd, "\n   res {\n");
    56. 

  56. 				fprintf(fd,
    57. 

  57. 					"          x0=%d, y0=%d, x1=%d, y1=%d, pw=%d, ph=%d, numbands=%d\n",
    58. 

  58. 					res->x0, res->y0, res->x1, res->y1, res->pw, res->ph, res->numbands);
    59. 

  59. 				for (bandno = 0; bandno < res->numbands; bandno++) {
    60. 

  60. 					opj_tcd_band_t *band = &res->bands[bandno];
    61. 

  61. 					fprintf(fd, "        band {\n");
    62. 

  62. 					fprintf(fd,
    63. 

  63. 						"          x0=%d, y0=%d, x1=%d, y1=%d, stepsize=%f, numbps=%d\n",
    64. 

  64. 						band->x0, band->y0, band->x1, band->y1, band->stepsize, band->numbps);
    65. 

  65. 					for (precno = 0; precno < res->pw * res->ph; precno++) {
    66. 

  66. 						opj_tcd_precinct_t *prec = &band->precincts[precno];
    67. 

  67. 						fprintf(fd, "          prec {\n");
    68. 

  68. 						fprintf(fd,
    69. 

  69. 							"            x0=%d, y0=%d, x1=%d, y1=%d, cw=%d, ch=%d\n",
    70. 

  70. 							prec->x0, prec->y0, prec->x1, prec->y1, prec->cw, prec->ch);
    71. 

  71. 						/*
    72. 

  72. 						for (cblkno = 0; cblkno < prec->cw * prec->ch; cblkno++) {
    73. 

  73. 							opj_tcd_cblk_t *cblk = &prec->cblks[cblkno];
    74. 

  74. 							fprintf(fd, "            cblk {\n");
    75. 

  75. 							fprintf(fd,
    76. 

  76. 								"              x0=%d, y0=%d, x1=%d, y1=%d\n",
    77. 

  77. 								cblk->x0, cblk->y0, cblk->x1, cblk->y1);
    78. 

  78. 							fprintf(fd, "            }\n");
    79. 

  79. 						}
    80. 

  80. 						*/
    81. 

  81. 						fprintf(fd, "          }\n");
    82. 

  82. 					}
    83. 

  83. 					fprintf(fd, "        }\n");
    84. 

  84. 				}
    85. 

  85. 				fprintf(fd, "      }\n");
    86. 

  86. 			}
    87. 

  87. 			fprintf(fd, "    }\n");
    88. 

  88. 		}
    89. 

  89. 		fprintf(fd, "  }\n");
    90. 

  90. 	}
    91. 

  91. 	fprintf(fd, "}\n");
    92. 

  92. }
    93. 

  93. 

  94. /* ----------------------------------------------------------------------- */
    95. 

  95. 

  96. /**
    97. 

  97. Create a new TCD handle
    98. 

  98. */
    99. 

  99. opj_tcd_t* tcd_create(opj_common_ptr cinfo) {
    100. 

  100. 	/* create the tcd structure */
    101. 

  101. 	opj_tcd_t *tcd = (opj_tcd_t*)opj_malloc(sizeof(opj_tcd_t));
    102. 

  102. 	if(!tcd) return NULL;
    103. 

  103. 	tcd->cinfo = cinfo;
    104. 

  104. 	tcd->tcd_image = (opj_tcd_image_t*)opj_malloc(sizeof(opj_tcd_image_t));
    105. 

  105. 	if(!tcd->tcd_image) {
    106. 

  106. 		opj_free(tcd);
    107. 

  107. 		return NULL;
    108. 

  108. 	}
    109. 

  109. 

  110. 	return tcd;
    111. 

  111. }
    112. 

  112. 

  113. /**
    114. 

  114. Destroy a previously created TCD handle
    115. 

  115. */
    116. 

  116. void tcd_destroy(opj_tcd_t *tcd) {
    117. 

  117. 	if(tcd) {
    118. 

  118. 		opj_free(tcd->tcd_image);
    119. 

  119. 		opj_free(tcd);
    120. 

  120. 	}
    121. 

  121. }
    122. 

  122. 

  123. /* ----------------------------------------------------------------------- */
    124. 

  124. 

  125. void tcd_malloc_encode(opj_tcd_t *tcd, opj_image_t * image, opj_cp_t * cp, int curtileno) {
    126. 

  126. 	int tileno, compno, resno, bandno, precno, cblkno;
    127. 

  127. 

  128. 	tcd->image = image;
    129. 

  129. 	tcd->cp = cp;
    130. 

  130. 	tcd->tcd_image->tw = cp->tw;
    131. 

  131. 	tcd->tcd_image->th = cp->th;
    132. 

  132. 	tcd->tcd_image->tiles = (opj_tcd_tile_t *) opj_malloc(sizeof(opj_tcd_tile_t));
    133. 

  133. 	
    134. 

  134. 	for (tileno = 0; tileno < 1; tileno++) {
    135. 

  135. 		opj_tcp_t *tcp = &cp->tcps[curtileno];
    136. 

  136. 		int j;
    137. 

  137. 

  138. 		/* cfr p59 ISO/IEC FDIS15444-1 : 2000 (18 august 2000) */
    139. 

  139. 		int p = curtileno % cp->tw;	/* si numerotation matricielle .. */
    140. 

  140. 		int q = curtileno / cp->tw;	/* .. coordonnees de la tile (q,p) q pour ligne et p pour colonne */
    141. 

  141. 

  142. 		/* opj_tcd_tile_t *tile=&tcd->tcd_image->tiles[tileno]; */
    143. 

  143. 		opj_tcd_tile_t *tile = tcd->tcd_image->tiles;
    144. 

  144. 

  145. 		/* 4 borders of the tile rescale on the image if necessary */
    146. 

  146. 		tile->x0 = int_max(cp->tx0 + p * cp->tdx, image->x0);
    147. 

  147. 		tile->y0 = int_max(cp->ty0 + q * cp->tdy, image->y0);
    148. 

  148. 		tile->x1 = int_min(cp->tx0 + (p + 1) * cp->tdx, image->x1);
    149. 

  149. 		tile->y1 = int_min(cp->ty0 + (q + 1) * cp->tdy, image->y1);
    150. 

  150. 		tile->numcomps = image->numcomps;
    151. 

  151. 		/* tile->PPT=image->PPT;  */
    152. 

  152. 

  153. 		/* Modification of the RATE >> */
    154. 

  154. 		for (j = 0; j < tcp->numlayers; j++) {
    155. 

  155. 			tcp->rates[j] = tcp->rates[j] ? 
    156. 

  156. 				cp->tp_on ? 
    157. 

  157. 					(((float) (tile->numcomps 
    158. 

  158. 					* (tile->x1 - tile->x0) 
    159. 

  159. 					* (tile->y1 - tile->y0)
    160. 

  160. 					* image->comps[0].prec))
    161. 

  161. 					/(tcp->rates[j] * 8 * image->comps[0].dx * image->comps[0].dy)) - (((tcd->cur_totnum_tp - 1) * 14 )/ tcp->numlayers)
    162. 

  162. 					:
    163. 

  163. 				((float) (tile->numcomps 
    164. 

  164. 					* (tile->x1 - tile->x0) 
    165. 

  165. 					* (tile->y1 - tile->y0) 
    166. 

  166. 					* image->comps[0].prec))/ 
    167. 

  167. 					(tcp->rates[j] * 8 * image->comps[0].dx * image->comps[0].dy)
    168. 

  168. 					: 0;
    169. 

  169. 

  170. 			if (tcp->rates[j]) {
    171. 

  171. 				if (j && tcp->rates[j] < tcp->rates[j - 1] + 10) {
    172. 

  172. 					tcp->rates[j] = tcp->rates[j - 1] + 20;
    173. 

  173. 				} else {
    174. 

  174. 					if (!j && tcp->rates[j] < 30)
    175. 

  175. 						tcp->rates[j] = 30;
    176. 

  176. 				}
    177. 

  177. 				
    178. 

  178. 				if(j == (tcp->numlayers-1)){
    179. 

  179. 					tcp->rates[j] = tcp->rates[j]- 2;
    180. 

  180. 				}
    181. 

  181. 			}
    182. 

  182. 		}
    183. 

  183. 		/* << Modification of the RATE */
    184. 

  184. 		
    185. 

  185. 		tile->comps = (opj_tcd_tilecomp_t *) opj_malloc(image->numcomps * sizeof(opj_tcd_tilecomp_t));
    186. 

  186. 		for (compno = 0; compno < tile->numcomps; compno++) {
    187. 

  187. 			opj_tccp_t *tccp = &tcp->tccps[compno];
    188. 

  188. 

  189. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    190. 

  190. 

  191. 			/* border of each tile component (global) */
    192. 

  192. 			tilec->x0 = int_ceildiv(tile->x0, image->comps[compno].dx);
    193. 

  193. 			tilec->y0 = int_ceildiv(tile->y0, image->comps[compno].dy);
    194. 

  194. 			tilec->x1 = int_ceildiv(tile->x1, image->comps[compno].dx);
    195. 

  195. 			tilec->y1 = int_ceildiv(tile->y1, image->comps[compno].dy);
    196. 

  196. 			
    197. 

  197. 			tilec->data = (int *) opj_aligned_malloc((tilec->x1 - tilec->x0) * (tilec->y1 - tilec->y0) * sizeof(int));
    198. 

  198. 			tilec->numresolutions = tccp->numresolutions;
    199. 

  199. 

  200. 			tilec->resolutions = (opj_tcd_resolution_t *) opj_malloc(tilec->numresolutions * sizeof(opj_tcd_resolution_t));
    201. 

  201. 			
    202. 

  202. 			for (resno = 0; resno < tilec->numresolutions; resno++) {
    203. 

  203. 				int pdx, pdy;
    204. 

  204. 				int levelno = tilec->numresolutions - 1 - resno;
    205. 

  205. 				int tlprcxstart, tlprcystart, brprcxend, brprcyend;
    206. 

  206. 				int tlcbgxstart, tlcbgystart, brcbgxend, brcbgyend;
    207. 

  207. 				int cbgwidthexpn, cbgheightexpn;
    208. 

  208. 				int cblkwidthexpn, cblkheightexpn;
    209. 

  209. 

  210. 				opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    211. 

  211. 				
    212. 

  212. 				/* border for each resolution level (global) */
    213. 

  213. 				res->x0 = int_ceildivpow2(tilec->x0, levelno);
    214. 

  214. 				res->y0 = int_ceildivpow2(tilec->y0, levelno);
    215. 

  215. 				res->x1 = int_ceildivpow2(tilec->x1, levelno);
    216. 

  216. 				res->y1 = int_ceildivpow2(tilec->y1, levelno);
    217. 

  217. 				
    218. 

  218. 				res->numbands = resno == 0 ? 1 : 3;
    219. 

  219. 				/* p. 35, table A-23, ISO/IEC FDIS154444-1 : 2000 (18 august 2000) */
    220. 

  220. 				if (tccp->csty & J2K_CCP_CSTY_PRT) {
    221. 

  221. 					pdx = tccp->prcw[resno];
    222. 

  222. 					pdy = tccp->prch[resno];
    223. 

  223. 				} else {
    224. 

  224. 					pdx = 15;
    225. 

  225. 					pdy = 15;
    226. 

  226. 				}
    227. 

  227. 				/* p. 64, B.6, ISO/IEC FDIS15444-1 : 2000 (18 august 2000)  */
    228. 

  228. 				tlprcxstart = int_floordivpow2(res->x0, pdx) << pdx;
    229. 

  229. 				tlprcystart = int_floordivpow2(res->y0, pdy) << pdy;
    230. 

  230. 				
    231. 

  231. 				brprcxend = int_ceildivpow2(res->x1, pdx) << pdx;
    232. 

  232. 				brprcyend = int_ceildivpow2(res->y1, pdy) << pdy;
    233. 

  233. 				
    234. 

  234. 				res->pw = (brprcxend - tlprcxstart) >> pdx;
    235. 

  235. 				res->ph = (brprcyend - tlprcystart) >> pdy;
    236. 

  236. 				
    237. 

  237. 				if (resno == 0) {
    238. 

  238. 					tlcbgxstart = tlprcxstart;
    239. 

  239. 					tlcbgystart = tlprcystart;
    240. 

  240. 					brcbgxend = brprcxend;
    241. 

  241. 					brcbgyend = brprcyend;
    242. 

  242. 					cbgwidthexpn = pdx;
    243. 

  243. 					cbgheightexpn = pdy;
    244. 

  244. 				} else {
    245. 

  245. 					tlcbgxstart = int_ceildivpow2(tlprcxstart, 1);
    246. 

  246. 					tlcbgystart = int_ceildivpow2(tlprcystart, 1);
    247. 

  247. 					brcbgxend = int_ceildivpow2(brprcxend, 1);
    248. 

  248. 					brcbgyend = int_ceildivpow2(brprcyend, 1);
    249. 

  249. 					cbgwidthexpn = pdx - 1;
    250. 

  250. 					cbgheightexpn = pdy - 1;
    251. 

  251. 				}
    252. 

  252. 				
    253. 

  253. 				cblkwidthexpn = int_min(tccp->cblkw, cbgwidthexpn);
    254. 

  254. 				cblkheightexpn = int_min(tccp->cblkh, cbgheightexpn);
    255. 

  255. 				
    256. 

  256. 				for (bandno = 0; bandno < res->numbands; bandno++) {
    257. 

  257. 					int x0b, y0b, i;
    258. 

  258. 					int gain, numbps;
    259. 

  259. 					opj_stepsize_t *ss = NULL;
    260. 

  260. 

  261. 					opj_tcd_band_t *band = &res->bands[bandno];
    262. 

  262. 

  263. 					band->bandno = resno == 0 ? 0 : bandno + 1;
    264. 

  264. 					x0b = (band->bandno == 1) || (band->bandno == 3) ? 1 : 0;
    265. 

  265. 					y0b = (band->bandno == 2) || (band->bandno == 3) ? 1 : 0;
    266. 

  266. 					
    267. 

  267. 					if (band->bandno == 0) {
    268. 

  268. 						/* band border (global) */
    269. 

  269. 						band->x0 = int_ceildivpow2(tilec->x0, levelno);
    270. 

  270. 						band->y0 = int_ceildivpow2(tilec->y0, levelno);
    271. 

  271. 						band->x1 = int_ceildivpow2(tilec->x1, levelno);
    272. 

  272. 						band->y1 = int_ceildivpow2(tilec->y1, levelno);
    273. 

  273. 					} else {
    274. 

  274. 						/* band border (global) */
    275. 

  275. 						band->x0 = int_ceildivpow2(tilec->x0 - (1 << levelno) * x0b, levelno + 1);
    276. 

  276. 						band->y0 = int_ceildivpow2(tilec->y0 - (1 << levelno) * y0b, levelno + 1);
    277. 

  277. 						band->x1 = int_ceildivpow2(tilec->x1 - (1 << levelno) * x0b, levelno + 1);
    278. 

  278. 						band->y1 = int_ceildivpow2(tilec->y1 - (1 << levelno) * y0b, levelno + 1);
    279. 

  279. 					}
    280. 

  280. 					
    281. 

  281. 					ss = &tccp->stepsizes[resno == 0 ? 0 : 3 * (resno - 1) + bandno + 1];
    282. 

  282. 					gain = tccp->qmfbid == 0 ? dwt_getgain_real(band->bandno) : dwt_getgain(band->bandno);					
    283. 

  283. 					numbps = image->comps[compno].prec + gain;
    284. 

  284. 					
    285. 

  285. 					band->stepsize = (float)((1.0 + ss->mant / 2048.0) * pow(2.0, numbps - ss->expn));
    286. 

  286. 					band->numbps = ss->expn + tccp->numgbits - 1;	/* WHY -1 ? */
    287. 

  287. 					
    288. 

  288. 					band->precincts = (opj_tcd_precinct_t *) opj_malloc(3 * res->pw * res->ph * sizeof(opj_tcd_precinct_t));
    289. 

  289. 					
    290. 

  290. 					for (i = 0; i < res->pw * res->ph * 3; i++) {
    291. 

  291. 						band->precincts[i].imsbtree = NULL;
    292. 

  292. 						band->precincts[i].incltree = NULL;
    293. 

  293. 						band->precincts[i].cblks.enc = NULL;
    294. 

  294. 					}
    295. 

  295. 					
    296. 

  296. 					for (precno = 0; precno < res->pw * res->ph; precno++) {
    297. 

  297. 						int tlcblkxstart, tlcblkystart, brcblkxend, brcblkyend;
    298. 

  298. 

  299. 						int cbgxstart = tlcbgxstart + (precno % res->pw) * (1 << cbgwidthexpn);
    300. 

  300. 						int cbgystart = tlcbgystart + (precno / res->pw) * (1 << cbgheightexpn);
    301. 

  301. 						int cbgxend = cbgxstart + (1 << cbgwidthexpn);
    302. 

  302. 						int cbgyend = cbgystart + (1 << cbgheightexpn);
    303. 

  303. 

  304. 						opj_tcd_precinct_t *prc = &band->precincts[precno];
    305. 

  305. 

  306. 						/* precinct size (global) */
    307. 

  307. 						prc->x0 = int_max(cbgxstart, band->x0);
    308. 

  308. 						prc->y0 = int_max(cbgystart, band->y0);
    309. 

  309. 						prc->x1 = int_min(cbgxend, band->x1);
    310. 

  310. 						prc->y1 = int_min(cbgyend, band->y1);
    311. 

  311. 

  312. 						tlcblkxstart = int_floordivpow2(prc->x0, cblkwidthexpn) << cblkwidthexpn;
    313. 

  313. 						tlcblkystart = int_floordivpow2(prc->y0, cblkheightexpn) << cblkheightexpn;
    314. 

  314. 						brcblkxend = int_ceildivpow2(prc->x1, cblkwidthexpn) << cblkwidthexpn;
    315. 

  315. 						brcblkyend = int_ceildivpow2(prc->y1, cblkheightexpn) << cblkheightexpn;
    316. 

  316. 						prc->cw = (brcblkxend - tlcblkxstart) >> cblkwidthexpn;
    317. 

  317. 						prc->ch = (brcblkyend - tlcblkystart) >> cblkheightexpn;
    318. 

  318. 

  319. 						prc->cblks.enc = (opj_tcd_cblk_enc_t*) opj_calloc((prc->cw * prc->ch), sizeof(opj_tcd_cblk_enc_t));
    320. 

  320. 						prc->incltree = tgt_create(prc->cw, prc->ch);
    321. 

  321. 						prc->imsbtree = tgt_create(prc->cw, prc->ch);
    322. 

  322. 						
    323. 

  323. 						for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    324. 

  324. 							int cblkxstart = tlcblkxstart + (cblkno % prc->cw) * (1 << cblkwidthexpn);
    325. 

  325. 							int cblkystart = tlcblkystart + (cblkno / prc->cw) * (1 << cblkheightexpn);
    326. 

  326. 							int cblkxend = cblkxstart + (1 << cblkwidthexpn);
    327. 

  327. 							int cblkyend = cblkystart + (1 << cblkheightexpn);
    328. 

  328. 							
    329. 

  329. 							opj_tcd_cblk_enc_t* cblk = &prc->cblks.enc[cblkno];
    330. 

  330. 

  331. 							/* code-block size (global) */
    332. 

  332. 							cblk->x0 = int_max(cblkxstart, prc->x0);
    333. 

  333. 							cblk->y0 = int_max(cblkystart, prc->y0);
    334. 

  334. 							cblk->x1 = int_min(cblkxend, prc->x1);
    335. 

  335. 							cblk->y1 = int_min(cblkyend, prc->y1);
    336. 

  336. 							cblk->data = (unsigned char*) opj_calloc(8192+2, sizeof(unsigned char));
    337. 

  337. 							/* FIXME: mqc_init_enc and mqc_byteout underrun the buffer if we don't do this. Why? */
    338. 

  338. 							cblk->data += 2;
    339. 

  339. 							cblk->layers = (opj_tcd_layer_t*) opj_calloc(100, sizeof(opj_tcd_layer_t));
    340. 

  340. 							cblk->passes = (opj_tcd_pass_t*) opj_calloc(100, sizeof(opj_tcd_pass_t));
    341. 

  341. 						}
    342. 

  342. 					}
    343. 

  343. 				}
    344. 

  344. 			}
    345. 

  345. 		}
    346. 

  346. 	}
    347. 

  347. 	
    348. 

  348. 	/* tcd_dump(stdout, tcd, &tcd->tcd_image); */
    349. 

  349. }
    350. 

  350. 

  351. void tcd_free_encode(opj_tcd_t *tcd) {
    352. 

  352. 	int tileno, compno, resno, bandno, precno, cblkno;
    353. 

  353. 

  354. 	for (tileno = 0; tileno < 1; tileno++) {
    355. 

  355. 		opj_tcd_tile_t *tile = tcd->tcd_image->tiles;
    356. 

  356. 

  357. 		for (compno = 0; compno < tile->numcomps; compno++) {
    358. 

  358. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    359. 

  359. 

  360. 			for (resno = 0; resno < tilec->numresolutions; resno++) {
    361. 

  361. 				opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    362. 

  362. 

  363. 				for (bandno = 0; bandno < res->numbands; bandno++) {
    364. 

  364. 					opj_tcd_band_t *band = &res->bands[bandno];
    365. 

  365. 

  366. 					for (precno = 0; precno < res->pw * res->ph; precno++) {
    367. 

  367. 						opj_tcd_precinct_t *prc = &band->precincts[precno];
    368. 

  368. 

  369. 						if (prc->incltree != NULL) {
    370. 

  370. 							tgt_destroy(prc->incltree);
    371. 

  371. 							prc->incltree = NULL;
    372. 

  372. 						}
    373. 

  373. 						if (prc->imsbtree != NULL) {
    374. 

  374. 							tgt_destroy(prc->imsbtree);	
    375. 

  375. 							prc->imsbtree = NULL;
    376. 

  376. 						}
    377. 

  377. 						for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    378. 

  378. 							opj_free(prc->cblks.enc[cblkno].data - 2);
    379. 

  379. 							opj_free(prc->cblks.enc[cblkno].layers);
    380. 

  380. 							opj_free(prc->cblks.enc[cblkno].passes);
    381. 

  381. 						}
    382. 

  382. 						opj_free(prc->cblks.enc);
    383. 

  383. 					} /* for (precno */
    384. 

  384. 					opj_free(band->precincts);
    385. 

  385. 					band->precincts = NULL;
    386. 

  386. 				} /* for (bandno */
    387. 

  387. 			} /* for (resno */
    388. 

  388. 			opj_free(tilec->resolutions);
    389. 

  389. 			tilec->resolutions = NULL;
    390. 

  390. 		} /* for (compno */
    391. 

  391. 		opj_free(tile->comps);
    392. 

  392. 		tile->comps = NULL;
    393. 

  393. 	} /* for (tileno */
    394. 

  394. 	opj_free(tcd->tcd_image->tiles);
    395. 

  395. 	tcd->tcd_image->tiles = NULL;
    396. 

  396. }
    397. 

  397. 

  398. void tcd_init_encode(opj_tcd_t *tcd, opj_image_t * image, opj_cp_t * cp, int curtileno) {
    399. 

  399. 	int tileno, compno, resno, bandno, precno, cblkno;
    400. 

  400. 

  401. 	for (tileno = 0; tileno < 1; tileno++) {
    402. 

  402. 		opj_tcp_t *tcp = &cp->tcps[curtileno];
    403. 

  403. 		int j;
    404. 

  404. 		/* cfr p59 ISO/IEC FDIS15444-1 : 2000 (18 august 2000) */
    405. 

  405. 		int p = curtileno % cp->tw;
    406. 

  406. 		int q = curtileno / cp->tw;
    407. 

  407. 

  408. 		opj_tcd_tile_t *tile = tcd->tcd_image->tiles;
    409. 

  409. 		
    410. 

  410. 		/* 4 borders of the tile rescale on the image if necessary */
    411. 

  411. 		tile->x0 = int_max(cp->tx0 + p * cp->tdx, image->x0);
    412. 

  412. 		tile->y0 = int_max(cp->ty0 + q * cp->tdy, image->y0);
    413. 

  413. 		tile->x1 = int_min(cp->tx0 + (p + 1) * cp->tdx, image->x1);
    414. 

  414. 		tile->y1 = int_min(cp->ty0 + (q + 1) * cp->tdy, image->y1);
    415. 

  415. 		
    416. 

  416. 		tile->numcomps = image->numcomps;
    417. 

  417. 		/* tile->PPT=image->PPT; */
    418. 

  418. 

  419. 		/* Modification of the RATE >> */
    420. 

  420. 		for (j = 0; j < tcp->numlayers; j++) {
    421. 

  421. 			tcp->rates[j] = tcp->rates[j] ? 
    422. 

  422. 				cp->tp_on ? 
    423. 

  423. 					(((float) (tile->numcomps 
    424. 

  424. 					* (tile->x1 - tile->x0) 
    425. 

  425. 					* (tile->y1 - tile->y0)
    426. 

  426. 					* image->comps[0].prec))
    427. 

  427. 					/(tcp->rates[j] * 8 * image->comps[0].dx * image->comps[0].dy)) - (((tcd->cur_totnum_tp - 1) * 14 )/ tcp->numlayers)
    428. 

  428. 					:
    429. 

  429. 				((float) (tile->numcomps 
    430. 

  430. 					* (tile->x1 - tile->x0) 
    431. 

  431. 					* (tile->y1 - tile->y0) 
    432. 

  432. 					* image->comps[0].prec))/ 
    433. 

  433. 					(tcp->rates[j] * 8 * image->comps[0].dx * image->comps[0].dy)
    434. 

  434. 					: 0;
    435. 

  435. 

  436. 			if (tcp->rates[j]) {
    437. 

  437. 				if (j && tcp->rates[j] < tcp->rates[j - 1] + 10) {
    438. 

  438. 					tcp->rates[j] = tcp->rates[j - 1] + 20;
    439. 

  439. 				} else {
    440. 

  440. 					if (!j && tcp->rates[j] < 30)
    441. 

  441. 						tcp->rates[j] = 30;
    442. 

  442. 				}
    443. 

  443. 			}
    444. 

  444. 		}
    445. 

  445. 		/* << Modification of the RATE */
    446. 

  446. 

  447. 		/* tile->comps=(opj_tcd_tilecomp_t*)opj_realloc(tile->comps,image->numcomps*sizeof(opj_tcd_tilecomp_t)); */
    448. 

  448. 		for (compno = 0; compno < tile->numcomps; compno++) {
    449. 

  449. 			opj_tccp_t *tccp = &tcp->tccps[compno];
    450. 

  450. 			
    451. 

  451. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    452. 

  452. 

  453. 			/* border of each tile component (global) */
    454. 

  454. 			tilec->x0 = int_ceildiv(tile->x0, image->comps[compno].dx);
    455. 

  455. 			tilec->y0 = int_ceildiv(tile->y0, image->comps[compno].dy);
    456. 

  456. 			tilec->x1 = int_ceildiv(tile->x1, image->comps[compno].dx);
    457. 

  457. 			tilec->y1 = int_ceildiv(tile->y1, image->comps[compno].dy);
    458. 

  458. 			
    459. 

  459. 			tilec->data = (int *) opj_aligned_malloc((tilec->x1 - tilec->x0) * (tilec->y1 - tilec->y0) * sizeof(int));
    460. 

  460. 			tilec->numresolutions = tccp->numresolutions;
    461. 

  461. 			/* tilec->resolutions=(opj_tcd_resolution_t*)opj_realloc(tilec->resolutions,tilec->numresolutions*sizeof(opj_tcd_resolution_t)); */
    462. 

  462. 			for (resno = 0; resno < tilec->numresolutions; resno++) {
    463. 

  463. 				int pdx, pdy;
    464. 

  464. 

  465. 				int levelno = tilec->numresolutions - 1 - resno;
    466. 

  466. 				int tlprcxstart, tlprcystart, brprcxend, brprcyend;
    467. 

  467. 				int tlcbgxstart, tlcbgystart, brcbgxend, brcbgyend;
    468. 

  468. 				int cbgwidthexpn, cbgheightexpn;
    469. 

  469. 				int cblkwidthexpn, cblkheightexpn;
    470. 

  470. 				
    471. 

  471. 				opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    472. 

  472. 

  473. 				/* border for each resolution level (global) */
    474. 

  474. 				res->x0 = int_ceildivpow2(tilec->x0, levelno);
    475. 

  475. 				res->y0 = int_ceildivpow2(tilec->y0, levelno);
    476. 

  476. 				res->x1 = int_ceildivpow2(tilec->x1, levelno);
    477. 

  477. 				res->y1 = int_ceildivpow2(tilec->y1, levelno);	
    478. 

  478. 				res->numbands = resno == 0 ? 1 : 3;
    479. 

  479. 

  480. 				/* p. 35, table A-23, ISO/IEC FDIS154444-1 : 2000 (18 august 2000) */
    481. 

  481. 				if (tccp->csty & J2K_CCP_CSTY_PRT) {
    482. 

  482. 					pdx = tccp->prcw[resno];
    483. 

  483. 					pdy = tccp->prch[resno];
    484. 

  484. 				} else {
    485. 

  485. 					pdx = 15;
    486. 

  486. 					pdy = 15;
    487. 

  487. 				}
    488. 

  488. 				/* p. 64, B.6, ISO/IEC FDIS15444-1 : 2000 (18 august 2000)  */
    489. 

  489. 				tlprcxstart = int_floordivpow2(res->x0, pdx) << pdx;
    490. 

  490. 				tlprcystart = int_floordivpow2(res->y0, pdy) << pdy;
    491. 

  491. 				brprcxend = int_ceildivpow2(res->x1, pdx) << pdx;
    492. 

  492. 				brprcyend = int_ceildivpow2(res->y1, pdy) << pdy;
    493. 

  493. 				
    494. 

  494. 				res->pw = (brprcxend - tlprcxstart) >> pdx;
    495. 

  495. 				res->ph = (brprcyend - tlprcystart) >> pdy;
    496. 

  496. 				
    497. 

  497. 				if (resno == 0) {
    498. 

  498. 					tlcbgxstart = tlprcxstart;
    499. 

  499. 					tlcbgystart = tlprcystart;
    500. 

  500. 					brcbgxend = brprcxend;
    501. 

  501. 					brcbgyend = brprcyend;
    502. 

  502. 					cbgwidthexpn = pdx;
    503. 

  503. 					cbgheightexpn = pdy;
    504. 

  504. 				} else {
    505. 

  505. 					tlcbgxstart = int_ceildivpow2(tlprcxstart, 1);
    506. 

  506. 					tlcbgystart = int_ceildivpow2(tlprcystart, 1);
    507. 

  507. 					brcbgxend = int_ceildivpow2(brprcxend, 1);
    508. 

  508. 					brcbgyend = int_ceildivpow2(brprcyend, 1);
    509. 

  509. 					cbgwidthexpn = pdx - 1;
    510. 

  510. 					cbgheightexpn = pdy - 1;
    511. 

  511. 				}
    512. 

  512. 				
    513. 

  513. 				cblkwidthexpn = int_min(tccp->cblkw, cbgwidthexpn);
    514. 

  514. 				cblkheightexpn = int_min(tccp->cblkh, cbgheightexpn);
    515. 

  515. 				
    516. 

  516. 				for (bandno = 0; bandno < res->numbands; bandno++) {
    517. 

  517. 					int x0b, y0b;
    518. 

  518. 					int gain, numbps;
    519. 

  519. 					opj_stepsize_t *ss = NULL;
    520. 

  520. 

  521. 					opj_tcd_band_t *band = &res->bands[bandno];
    522. 

  522. 

  523. 					band->bandno = resno == 0 ? 0 : bandno + 1;
    524. 

  524. 					x0b = (band->bandno == 1) || (band->bandno == 3) ? 1 : 0;
    525. 

  525. 					y0b = (band->bandno == 2) || (band->bandno == 3) ? 1 : 0;
    526. 

  526. 					
    527. 

  527. 					if (band->bandno == 0) {
    528. 

  528. 						/* band border */
    529. 

  529. 						band->x0 = int_ceildivpow2(tilec->x0, levelno);
    530. 

  530. 						band->y0 = int_ceildivpow2(tilec->y0, levelno);
    531. 

  531. 						band->x1 = int_ceildivpow2(tilec->x1, levelno);
    532. 

  532. 						band->y1 = int_ceildivpow2(tilec->y1, levelno);
    533. 

  533. 					} else {
    534. 

  534. 						band->x0 = int_ceildivpow2(tilec->x0 - (1 << levelno) * x0b, levelno + 1);
    535. 

  535. 						band->y0 = int_ceildivpow2(tilec->y0 - (1 << levelno) * y0b, levelno + 1);
    536. 

  536. 						band->x1 = int_ceildivpow2(tilec->x1 - (1 << levelno) * x0b, levelno + 1);
    537. 

  537. 						band->y1 = int_ceildivpow2(tilec->y1 - (1 << levelno) * y0b, levelno + 1);
    538. 

  538. 					}
    539. 

  539. 					
    540. 

  540. 					ss = &tccp->stepsizes[resno == 0 ? 0 : 3 * (resno - 1) + bandno + 1];
    541. 

  541. 					gain = tccp->qmfbid == 0 ? dwt_getgain_real(band->bandno) : dwt_getgain(band->bandno);
    542. 

  542. 					numbps = image->comps[compno].prec + gain;
    543. 

  543. 					band->stepsize = (float)((1.0 + ss->mant / 2048.0) * pow(2.0, numbps - ss->expn));
    544. 

  544. 					band->numbps = ss->expn + tccp->numgbits - 1;	/* WHY -1 ? */
    545. 

  545. 					
    546. 

  546. 					for (precno = 0; precno < res->pw * res->ph; precno++) {
    547. 

  547. 						int tlcblkxstart, tlcblkystart, brcblkxend, brcblkyend;
    548. 

  548. 

  549. 						int cbgxstart = tlcbgxstart + (precno % res->pw) * (1 << cbgwidthexpn);
    550. 

  550. 						int cbgystart = tlcbgystart + (precno / res->pw) * (1 << cbgheightexpn);
    551. 

  551. 						int cbgxend = cbgxstart + (1 << cbgwidthexpn);
    552. 

  552. 						int cbgyend = cbgystart + (1 << cbgheightexpn);
    553. 

  553. 						
    554. 

  554. 						opj_tcd_precinct_t *prc = &band->precincts[precno];
    555. 

  555. 

  556. 						/* precinct size (global) */
    557. 

  557. 						prc->x0 = int_max(cbgxstart, band->x0);
    558. 

  558. 						prc->y0 = int_max(cbgystart, band->y0);
    559. 

  559. 						prc->x1 = int_min(cbgxend, band->x1);
    560. 

  560. 						prc->y1 = int_min(cbgyend, band->y1);
    561. 

  561. 

  562. 						tlcblkxstart = int_floordivpow2(prc->x0, cblkwidthexpn) << cblkwidthexpn;
    563. 

  563. 						tlcblkystart = int_floordivpow2(prc->y0, cblkheightexpn) << cblkheightexpn;
    564. 

  564. 						brcblkxend = int_ceildivpow2(prc->x1, cblkwidthexpn) << cblkwidthexpn;
    565. 

  565. 						brcblkyend = int_ceildivpow2(prc->y1, cblkheightexpn) << cblkheightexpn;
    566. 

  566. 						prc->cw = (brcblkxend - tlcblkxstart) >> cblkwidthexpn;
    567. 

  567. 						prc->ch = (brcblkyend - tlcblkystart) >> cblkheightexpn;
    568. 

  568. 

  569. 						opj_free(prc->cblks.enc);
    570. 

  570. 						prc->cblks.enc = (opj_tcd_cblk_enc_t*) opj_calloc(prc->cw * prc->ch, sizeof(opj_tcd_cblk_enc_t));
    571. 

  571. 

  572. 						if (prc->incltree != NULL) {
    573. 

  573. 							tgt_destroy(prc->incltree);
    574. 

  574. 						}
    575. 

  575. 						if (prc->imsbtree != NULL) {
    576. 

  576. 							tgt_destroy(prc->imsbtree);
    577. 

  577. 						}
    578. 

  578. 						
    579. 

  579. 						prc->incltree = tgt_create(prc->cw, prc->ch);
    580. 

  580. 						prc->imsbtree = tgt_create(prc->cw, prc->ch);
    581. 

  581. 

  582. 						for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    583. 

  583. 							int cblkxstart = tlcblkxstart + (cblkno % prc->cw) * (1 << cblkwidthexpn);
    584. 

  584. 							int cblkystart = tlcblkystart + (cblkno / prc->cw) * (1 << cblkheightexpn);
    585. 

  585. 							int cblkxend = cblkxstart + (1 << cblkwidthexpn);
    586. 

  586. 							int cblkyend = cblkystart + (1 << cblkheightexpn);
    587. 

  587. 

  588. 							opj_tcd_cblk_enc_t* cblk = &prc->cblks.enc[cblkno];
    589. 

  589. 

  590. 							/* code-block size (global) */
    591. 

  591. 							cblk->x0 = int_max(cblkxstart, prc->x0);
    592. 

  592. 							cblk->y0 = int_max(cblkystart, prc->y0);
    593. 

  593. 							cblk->x1 = int_min(cblkxend, prc->x1);
    594. 

  594. 							cblk->y1 = int_min(cblkyend, prc->y1);
    595. 

  595. 							cblk->data = (unsigned char*) opj_calloc(8192+2, sizeof(unsigned char));
    596. 

  596. 							/* FIXME: mqc_init_enc and mqc_byteout underrun the buffer if we don't do this. Why? */
    597. 

  597. 							cblk->data += 2;
    598. 

  598. 							cblk->layers = (opj_tcd_layer_t*) opj_calloc(100, sizeof(opj_tcd_layer_t));
    599. 

  599. 							cblk->passes = (opj_tcd_pass_t*) opj_calloc(100, sizeof(opj_tcd_pass_t));
    600. 

  600. 						}
    601. 

  601. 					} /* precno */
    602. 

  602. 				} /* bandno */
    603. 

  603. 			} /* resno */
    604. 

  604. 		} /* compno */
    605. 

  605. 	} /* tileno */
    606. 

  606. 

  607. 	/* tcd_dump(stdout, tcd, &tcd->tcd_image); */
    608. 

  608. }
    609. 

  609. 

  610. void tcd_malloc_decode(opj_tcd_t *tcd, opj_image_t * image, opj_cp_t * cp) {
    611. 

  611. 	int i, j, tileno, p, q;
    612. 

  612. 	unsigned int x0 = 0, y0 = 0, x1 = 0, y1 = 0, w, h;
    613. 

  613. 

  614. 	tcd->image = image;
    615. 

  615. 	tcd->tcd_image->tw = cp->tw;
    616. 

  616. 	tcd->tcd_image->th = cp->th;
    617. 

  617. 	tcd->tcd_image->tiles = (opj_tcd_tile_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tcd_tile_t));
    618. 

  618. 

  619. 	/* 
    620. 

  620. 	Allocate place to store the decoded data = final image
    621. 

  621. 	Place limited by the tile really present in the codestream 
    622. 

  622. 	*/
    623. 

  623. 

  624. 	for (j = 0; j < cp->tileno_size; j++) {
    625. 

  625. 		opj_tcd_tile_t *tile;
    626. 

  626. 		
    627. 

  627. 		tileno = cp->tileno[j];		
    628. 

  628. 		tile = &(tcd->tcd_image->tiles[cp->tileno[tileno]]);		
    629. 

  629. 		tile->numcomps = image->numcomps;
    630. 

  630. 		tile->comps = (opj_tcd_tilecomp_t*) opj_calloc(image->numcomps, sizeof(opj_tcd_tilecomp_t));
    631. 

  631. 	}
    632. 

  632. 

  633. 	for (i = 0; i < image->numcomps; i++) {
    634. 

  634. 		for (j = 0; j < cp->tileno_size; j++) {
    635. 

  635. 			opj_tcd_tile_t *tile;
    636. 

  636. 			opj_tcd_tilecomp_t *tilec;
    637. 

  637. 			
    638. 

  638. 			/* cfr p59 ISO/IEC FDIS15444-1 : 2000 (18 august 2000) */
    639. 

  639. 			
    640. 

  640. 			tileno = cp->tileno[j];
    641. 

  641. 			
    642. 

  642. 			tile = &(tcd->tcd_image->tiles[cp->tileno[tileno]]);
    643. 

  643. 			tilec = &tile->comps[i];
    644. 

  644. 			
    645. 

  645. 			p = tileno % cp->tw;	/* si numerotation matricielle .. */
    646. 

  646. 			q = tileno / cp->tw;	/* .. coordonnees de la tile (q,p) q pour ligne et p pour colonne */
    647. 

  647. 			
    648. 

  648. 			/* 4 borders of the tile rescale on the image if necessary */
    649. 

  649. 			tile->x0 = int_max(cp->tx0 + p * cp->tdx, image->x0);
    650. 

  650. 			tile->y0 = int_max(cp->ty0 + q * cp->tdy, image->y0);
    651. 

  651. 			tile->x1 = int_min(cp->tx0 + (p + 1) * cp->tdx, image->x1);
    652. 

  652. 			tile->y1 = int_min(cp->ty0 + (q + 1) * cp->tdy, image->y1);
    653. 

  653. 

  654. 			tilec->x0 = int_ceildiv(tile->x0, image->comps[i].dx);
    655. 

  655. 			tilec->y0 = int_ceildiv(tile->y0, image->comps[i].dy);
    656. 

  656. 			tilec->x1 = int_ceildiv(tile->x1, image->comps[i].dx);
    657. 

  657. 			tilec->y1 = int_ceildiv(tile->y1, image->comps[i].dy);
    658. 

  658. 

  659. 			x0 = j == 0 ? tilec->x0 : int_min(x0, (unsigned int) tilec->x0);
    660. 

  660. 			y0 = j == 0 ? tilec->y0 : int_min(y0,	(unsigned int) tilec->y0);
    661. 

  661. 			x1 = j == 0 ? tilec->x1 : int_max(x1,	(unsigned int) tilec->x1);
    662. 

  662. 			y1 = j == 0 ? tilec->y1 : int_max(y1,	(unsigned int) tilec->y1);
    663. 

  663. 		}
    664. 

  664. 

  665. 		w = int_ceildivpow2(x1 - x0, image->comps[i].factor);
    666. 

  666. 		h = int_ceildivpow2(y1 - y0, image->comps[i].factor);
    667. 

  667. 

  668. 		image->comps[i].w = w;
    669. 

  669. 		image->comps[i].h = h;
    670. 

  670. 		image->comps[i].x0 = x0;
    671. 

  671. 		image->comps[i].y0 = y0;
    672. 

  672. 	}
    673. 

  673. }
    674. 

  674. 

  675. void tcd_malloc_decode_tile(opj_tcd_t *tcd, opj_image_t * image, opj_cp_t * cp, int tileno, opj_codestream_info_t *cstr_info) {
    676. 

  676. 	int compno, resno, bandno, precno, cblkno;
    677. 

  677. 	opj_tcp_t *tcp;
    678. 

  678. 	opj_tcd_tile_t *tile;
    679. 

  679. 

  680. 	OPJ_ARG_NOT_USED(cstr_info);
    681. 

  681. 

  682. 	tcd->cp = cp;
    683. 

  683. 	
    684. 

  684. 	tcp = &(cp->tcps[cp->tileno[tileno]]);
    685. 

  685. 	tile = &(tcd->tcd_image->tiles[cp->tileno[tileno]]);
    686. 

  686. 	
    687. 

  687. 	tileno = cp->tileno[tileno];
    688. 

  688. 	
    689. 

  689. 	for (compno = 0; compno < tile->numcomps; compno++) {
    690. 

  690. 		opj_tccp_t *tccp = &tcp->tccps[compno];
    691. 

  691. 		opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    692. 

  692. 		
    693. 

  693. 		/* border of each tile component (global) */
    694. 

  694. 		tilec->x0 = int_ceildiv(tile->x0, image->comps[compno].dx);
    695. 

  695. 		tilec->y0 = int_ceildiv(tile->y0, image->comps[compno].dy);
    696. 

  696. 		tilec->x1 = int_ceildiv(tile->x1, image->comps[compno].dx);
    697. 

  697. 		tilec->y1 = int_ceildiv(tile->y1, image->comps[compno].dy);
    698. 

  698. 

  699. 		tilec->numresolutions = tccp->numresolutions;
    700. 

  700. 		tilec->resolutions = (opj_tcd_resolution_t *) opj_malloc(tilec->numresolutions * sizeof(opj_tcd_resolution_t));
    701. 

  701. 		
    702. 

  702. 		for (resno = 0; resno < tilec->numresolutions; resno++) {
    703. 

  703. 			int pdx, pdy;
    704. 

  704. 			int levelno = tilec->numresolutions - 1 - resno;
    705. 

  705. 			int tlprcxstart, tlprcystart, brprcxend, brprcyend;
    706. 

  706. 			int tlcbgxstart, tlcbgystart, brcbgxend, brcbgyend;
    707. 

  707. 			int cbgwidthexpn, cbgheightexpn;
    708. 

  708. 			int cblkwidthexpn, cblkheightexpn;
    709. 

  709. 			
    710. 

  710. 			opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    711. 

  711. 			
    712. 

  712. 			/* border for each resolution level (global) */
    713. 

  713. 			res->x0 = int_ceildivpow2(tilec->x0, levelno);
    714. 

  714. 			res->y0 = int_ceildivpow2(tilec->y0, levelno);
    715. 

  715. 			res->x1 = int_ceildivpow2(tilec->x1, levelno);
    716. 

  716. 			res->y1 = int_ceildivpow2(tilec->y1, levelno);
    717. 

  717. 			res->numbands = resno == 0 ? 1 : 3;
    718. 

  718. 			
    719. 

  719. 			/* p. 35, table A-23, ISO/IEC FDIS154444-1 : 2000 (18 august 2000) */
    720. 

  720. 			if (tccp->csty & J2K_CCP_CSTY_PRT) {
    721. 

  721. 				pdx = tccp->prcw[resno];
    722. 

  722. 				pdy = tccp->prch[resno];
    723. 

  723. 			} else {
    724. 

  724. 				pdx = 15;
    725. 

  725. 				pdy = 15;
    726. 

  726. 			}			
    727. 

  727. 			
    728. 

  728. 			/* p. 64, B.6, ISO/IEC FDIS15444-1 : 2000 (18 august 2000)  */
    729. 

  729. 			tlprcxstart = int_floordivpow2(res->x0, pdx) << pdx;
    730. 

  730. 			tlprcystart = int_floordivpow2(res->y0, pdy) << pdy;
    731. 

  731. 			brprcxend = int_ceildivpow2(res->x1, pdx) << pdx;
    732. 

  732. 			brprcyend = int_ceildivpow2(res->y1, pdy) << pdy;
    733. 

  733. 			
    734. 

  734. 			res->pw = (res->x0 == res->x1) ? 0 : ((brprcxend - tlprcxstart) >> pdx);
    735. 

  735. 			res->ph = (res->y0 == res->y1) ? 0 : ((brprcyend - tlprcystart) >> pdy);
    736. 

  736. 			
    737. 

  737. 			if (resno == 0) {
    738. 

  738. 				tlcbgxstart = tlprcxstart;
    739. 

  739. 				tlcbgystart = tlprcystart;
    740. 

  740. 				brcbgxend = brprcxend;
    741. 

  741. 				brcbgyend = brprcyend;
    742. 

  742. 				cbgwidthexpn = pdx;
    743. 

  743. 				cbgheightexpn = pdy;
    744. 

  744. 			} else {
    745. 

  745. 				tlcbgxstart = int_ceildivpow2(tlprcxstart, 1);
    746. 

  746. 				tlcbgystart = int_ceildivpow2(tlprcystart, 1);
    747. 

  747. 				brcbgxend = int_ceildivpow2(brprcxend, 1);
    748. 

  748. 				brcbgyend = int_ceildivpow2(brprcyend, 1);
    749. 

  749. 				cbgwidthexpn = pdx - 1;
    750. 

  750. 				cbgheightexpn = pdy - 1;
    751. 

  751. 			}
    752. 

  752. 			
    753. 

  753. 			cblkwidthexpn = int_min(tccp->cblkw, cbgwidthexpn);
    754. 

  754. 			cblkheightexpn = int_min(tccp->cblkh, cbgheightexpn);
    755. 

  755. 			
    756. 

  756. 			for (bandno = 0; bandno < res->numbands; bandno++) {
    757. 

  757. 				int x0b, y0b;
    758. 

  758. 				int gain, numbps;
    759. 

  759. 				opj_stepsize_t *ss = NULL;
    760. 

  760. 				
    761. 

  761. 				opj_tcd_band_t *band = &res->bands[bandno];
    762. 

  762. 				band->bandno = resno == 0 ? 0 : bandno + 1;
    763. 

  763. 				x0b = (band->bandno == 1) || (band->bandno == 3) ? 1 : 0;
    764. 

  764. 				y0b = (band->bandno == 2) || (band->bandno == 3) ? 1 : 0;
    765. 

  765. 				
    766. 

  766. 				if (band->bandno == 0) {
    767. 

  767. 					/* band border (global) */
    768. 

  768. 					band->x0 = int_ceildivpow2(tilec->x0, levelno);
    769. 

  769. 					band->y0 = int_ceildivpow2(tilec->y0, levelno);
    770. 

  770. 					band->x1 = int_ceildivpow2(tilec->x1, levelno);
    771. 

  771. 					band->y1 = int_ceildivpow2(tilec->y1, levelno);
    772. 

  772. 				} else {
    773. 

  773. 					/* band border (global) */
    774. 

  774. 					band->x0 = int_ceildivpow2(tilec->x0 - (1 << levelno) * x0b, levelno + 1);
    775. 

  775. 					band->y0 = int_ceildivpow2(tilec->y0 - (1 << levelno) * y0b, levelno + 1);
    776. 

  776. 					band->x1 = int_ceildivpow2(tilec->x1 - (1 << levelno) * x0b, levelno + 1);
    777. 

  777. 					band->y1 = int_ceildivpow2(tilec->y1 - (1 << levelno) * y0b, levelno + 1);
    778. 

  778. 				}
    779. 

  779. 				
    780. 

  780. 				ss = &tccp->stepsizes[resno == 0 ? 0 : 3 * (resno - 1) + bandno + 1];
    781. 

  781. 				gain = tccp->qmfbid == 0 ? dwt_getgain_real(band->bandno) : dwt_getgain(band->bandno);
    782. 

  782. 				numbps = image->comps[compno].prec + gain;
    783. 

  783. 				band->stepsize = (float)(((1.0 + ss->mant / 2048.0) * pow(2.0, numbps - ss->expn)) * 0.5);
    784. 

  784. 				band->numbps = ss->expn + tccp->numgbits - 1;	/* WHY -1 ? */
    785. 

  785. 				
    786. 

  786. 				band->precincts = (opj_tcd_precinct_t *) opj_malloc(res->pw * res->ph * sizeof(opj_tcd_precinct_t));
    787. 

  787. 				
    788. 

  788. 				for (precno = 0; precno < res->pw * res->ph; precno++) {
    789. 

  789. 					int tlcblkxstart, tlcblkystart, brcblkxend, brcblkyend;
    790. 

  790. 					int cbgxstart = tlcbgxstart + (precno % res->pw) * (1 << cbgwidthexpn);
    791. 

  791. 					int cbgystart = tlcbgystart + (precno / res->pw) * (1 << cbgheightexpn);
    792. 

  792. 					int cbgxend = cbgxstart + (1 << cbgwidthexpn);
    793. 

  793. 					int cbgyend = cbgystart + (1 << cbgheightexpn);
    794. 

  794. 					
    795. 

  795. 					opj_tcd_precinct_t *prc = &band->precincts[precno];
    796. 

  796. 					/* precinct size (global) */
    797. 

  797. 					prc->x0 = int_max(cbgxstart, band->x0);
    798. 

  798. 					prc->y0 = int_max(cbgystart, band->y0);
    799. 

  799. 					prc->x1 = int_min(cbgxend, band->x1);
    800. 

  800. 					prc->y1 = int_min(cbgyend, band->y1);
    801. 

  801. 					
    802. 

  802. 					tlcblkxstart = int_floordivpow2(prc->x0, cblkwidthexpn) << cblkwidthexpn;
    803. 

  803. 					tlcblkystart = int_floordivpow2(prc->y0, cblkheightexpn) << cblkheightexpn;
    804. 

  804. 					brcblkxend = int_ceildivpow2(prc->x1, cblkwidthexpn) << cblkwidthexpn;
    805. 

  805. 					brcblkyend = int_ceildivpow2(prc->y1, cblkheightexpn) << cblkheightexpn;
    806. 

  806. 					prc->cw = (brcblkxend - tlcblkxstart) >> cblkwidthexpn;
    807. 

  807. 					prc->ch = (brcblkyend - tlcblkystart) >> cblkheightexpn;
    808. 

  808. 

  809. 					prc->cblks.dec = (opj_tcd_cblk_dec_t*) opj_malloc(prc->cw * prc->ch * sizeof(opj_tcd_cblk_dec_t));
    810. 

  810. 

  811. 					prc->incltree = tgt_create(prc->cw, prc->ch);
    812. 

  812. 					prc->imsbtree = tgt_create(prc->cw, prc->ch);
    813. 

  813. 					
    814. 

  814. 					for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    815. 

  815. 						int cblkxstart = tlcblkxstart + (cblkno % prc->cw) * (1 << cblkwidthexpn);
    816. 

  816. 						int cblkystart = tlcblkystart + (cblkno / prc->cw) * (1 << cblkheightexpn);
    817. 

  817. 						int cblkxend = cblkxstart + (1 << cblkwidthexpn);
    818. 

  818. 						int cblkyend = cblkystart + (1 << cblkheightexpn);					
    819. 

  819. 

  820. 						opj_tcd_cblk_dec_t* cblk = &prc->cblks.dec[cblkno];
    821. 

  821. 						cblk->data = NULL;
    822. 

  822. 						cblk->segs = NULL;
    823. 

  823. 						/* code-block size (global) */
    824. 

  824. 						cblk->x0 = int_max(cblkxstart, prc->x0);
    825. 

  825. 						cblk->y0 = int_max(cblkystart, prc->y0);
    826. 

  826. 						cblk->x1 = int_min(cblkxend, prc->x1);
    827. 

  827. 						cblk->y1 = int_min(cblkyend, prc->y1);
    828. 

  828. 						cblk->numsegs = 0;
    829. 

  829. 					}
    830. 

  830. 				} /* precno */
    831. 

  831. 			} /* bandno */
    832. 

  832. 		} /* resno */
    833. 

  833. 	} /* compno */
    834. 

  834. 	/* tcd_dump(stdout, tcd, &tcd->tcd_image); */
    835. 

  835. }
    836. 

  836. 

  837. void tcd_makelayer_fixed(opj_tcd_t *tcd, int layno, int final) {
    838. 

  838. 	int compno, resno, bandno, precno, cblkno;
    839. 

  839. 	int value;			/*, matrice[tcd_tcp->numlayers][tcd_tile->comps[0].numresolutions][3]; */
    840. 

  840. 	int matrice[10][10][3];
    841. 

  841. 	int i, j, k;
    842. 

  842. 

  843. 	opj_cp_t *cp = tcd->cp;
    844. 

  844. 	opj_tcd_tile_t *tcd_tile = tcd->tcd_tile;
    845. 

  845. 	opj_tcp_t *tcd_tcp = tcd->tcp;
    846. 

  846. 

  847. 	/*matrice=(int*)opj_malloc(tcd_tcp->numlayers*tcd_tile->comps[0].numresolutions*3*sizeof(int)); */
    848. 

  848. 	
    849. 

  849. 	for (compno = 0; compno < tcd_tile->numcomps; compno++) {
    850. 

  850. 		opj_tcd_tilecomp_t *tilec = &tcd_tile->comps[compno];
    851. 

  851. 		for (i = 0; i < tcd_tcp->numlayers; i++) {
    852. 

  852. 			for (j = 0; j < tilec->numresolutions; j++) {
    853. 

  853. 				for (k = 0; k < 3; k++) {
    854. 

  854. 					matrice[i][j][k] =
    855. 

  855. 						(int) (cp->matrice[i * tilec->numresolutions * 3 + j * 3 + k] 
    856. 

  856. 						* (float) (tcd->image->comps[compno].prec / 16.0));
    857. 

  857. 				}
    858. 

  858. 			}
    859. 

  859. 		}
    860. 

  860.         
    861. 

  861. 		for (resno = 0; resno < tilec->numresolutions; resno++) {
    862. 

  862. 			opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    863. 

  863. 			for (bandno = 0; bandno < res->numbands; bandno++) {
    864. 

  864. 				opj_tcd_band_t *band = &res->bands[bandno];
    865. 

  865. 				for (precno = 0; precno < res->pw * res->ph; precno++) {
    866. 

  866. 					opj_tcd_precinct_t *prc = &band->precincts[precno];
    867. 

  867. 					for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    868. 

  868. 						opj_tcd_cblk_enc_t *cblk = &prc->cblks.enc[cblkno];
    869. 

  869. 						opj_tcd_layer_t *layer = &cblk->layers[layno];
    870. 

  870. 						int n;
    871. 

  871. 						int imsb = tcd->image->comps[compno].prec - cblk->numbps;	/* number of bit-plan equal to zero */
    872. 

  872. 						/* Correction of the matrix of coefficient to include the IMSB information */
    873. 

  873. 						if (layno == 0) {
    874. 

  874. 							value = matrice[layno][resno][bandno];
    875. 

  875. 							if (imsb >= value) {
    876. 

  876. 								value = 0;
    877. 

  877. 							} else {
    878. 

  878. 								value -= imsb;
    879. 

  879. 							}
    880. 

  880. 						} else {
    881. 

  881. 							value =	matrice[layno][resno][bandno] -	matrice[layno - 1][resno][bandno];
    882. 

  882. 							if (imsb >= matrice[layno - 1][resno][bandno]) {
    883. 

  883. 								value -= (imsb - matrice[layno - 1][resno][bandno]);
    884. 

  884. 								if (value < 0) {
    885. 

  885. 									value = 0;
    886. 

  886. 								}
    887. 

  887. 							}
    888. 

  888. 						}
    889. 

  889. 						
    890. 

  890. 						if (layno == 0) {
    891. 

  891. 							cblk->numpassesinlayers = 0;
    892. 

  892. 						}
    893. 

  893. 						
    894. 

  894. 						n = cblk->numpassesinlayers;
    895. 

  895. 						if (cblk->numpassesinlayers == 0) {
    896. 

  896. 							if (value != 0) {
    897. 

  897. 								n = 3 * value - 2 + cblk->numpassesinlayers;
    898. 

  898. 							} else {
    899. 

  899. 								n = cblk->numpassesinlayers;
    900. 

  900. 							}
    901. 

  901. 						} else {
    902. 

  902. 							n = 3 * value + cblk->numpassesinlayers;
    903. 

  903. 						}
    904. 

  904. 						
    905. 

  905. 						layer->numpasses = n - cblk->numpassesinlayers;
    906. 

  906. 						
    907. 

  907. 						if (!layer->numpasses)
    908. 

  908. 							continue;
    909. 

  909. 						
    910. 

  910. 						if (cblk->numpassesinlayers == 0) {
    911. 

  911. 							layer->len = cblk->passes[n - 1].rate;
    912. 

  912. 							layer->data = cblk->data;
    913. 

  913. 						} else {
    914. 

  914. 							layer->len = cblk->passes[n - 1].rate - cblk->passes[cblk->numpassesinlayers - 1].rate;
    915. 

  915. 							layer->data = cblk->data + cblk->passes[cblk->numpassesinlayers - 1].rate;
    916. 

  916. 						}
    917. 

  917. 						if (final)
    918. 

  918. 							cblk->numpassesinlayers = n;
    919. 

  919. 					}
    920. 

  920. 				}
    921. 

  921. 			}
    922. 

  922. 		}
    923. 

  923. 	}
    924. 

  924. }
    925. 

  925. 

  926. void tcd_rateallocate_fixed(opj_tcd_t *tcd) {
    927. 

  927. 	int layno;
    928. 

  928. 	for (layno = 0; layno < tcd->tcp->numlayers; layno++) {
    929. 

  929. 		tcd_makelayer_fixed(tcd, layno, 1);
    930. 

  930. 	}
    931. 

  931. }
    932. 

  932. 

  933. void tcd_makelayer(opj_tcd_t *tcd, int layno, double thresh, int final) {
    934. 

  934. 	int compno, resno, bandno, precno, cblkno, passno;
    935. 

  935. 	
    936. 

  936. 	opj_tcd_tile_t *tcd_tile = tcd->tcd_tile;
    937. 

  937. 

  938. 	tcd_tile->distolayer[layno] = 0;	/* fixed_quality */
    939. 

  939. 	
    940. 

  940. 	for (compno = 0; compno < tcd_tile->numcomps; compno++) {
    941. 

  941. 		opj_tcd_tilecomp_t *tilec = &tcd_tile->comps[compno];
    942. 

  942. 		for (resno = 0; resno < tilec->numresolutions; resno++) {
    943. 

  943. 			opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    944. 

  944. 			for (bandno = 0; bandno < res->numbands; bandno++) {
    945. 

  945. 				opj_tcd_band_t *band = &res->bands[bandno];
    946. 

  946. 				for (precno = 0; precno < res->pw * res->ph; precno++) {
    947. 

  947. 					opj_tcd_precinct_t *prc = &band->precincts[precno];
    948. 

  948. 					for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    949. 

  949. 						opj_tcd_cblk_enc_t *cblk = &prc->cblks.enc[cblkno];
    950. 

  950. 						opj_tcd_layer_t *layer = &cblk->layers[layno];
    951. 

  951. 						
    952. 

  952. 						int n;
    953. 

  953. 						if (layno == 0) {
    954. 

  954. 							cblk->numpassesinlayers = 0;
    955. 

  955. 						}
    956. 

  956. 						n = cblk->numpassesinlayers;
    957. 

  957. 						for (passno = cblk->numpassesinlayers; passno < cblk->totalpasses; passno++) {
    958. 

  958. 							int dr;
    959. 

  959. 							double dd;
    960. 

  960. 							opj_tcd_pass_t *pass = &cblk->passes[passno];
    961. 

  961. 							if (n == 0) {
    962. 

  962. 								dr = pass->rate;
    963. 

  963. 								dd = pass->distortiondec;
    964. 

  964. 							} else {
    965. 

  965. 								dr = pass->rate - cblk->passes[n - 1].rate;
    966. 

  966. 								dd = pass->distortiondec - cblk->passes[n - 1].distortiondec;
    967. 

  967. 							}
    968. 

  968. 							if (!dr) {
    969. 

  969. 								if (dd != 0)
    970. 

  970. 									n = passno + 1;
    971. 

  971. 								continue;
    972. 

  972. 							}
    973. 

  973. 							if (dd / dr >= thresh)
    974. 

  974. 								n = passno + 1;
    975. 

  975. 						}
    976. 

  976. 						layer->numpasses = n - cblk->numpassesinlayers;
    977. 

  977. 						
    978. 

  978. 						if (!layer->numpasses) {
    979. 

  979. 							layer->disto = 0;
    980. 

  980. 							continue;
    981. 

  981. 						}
    982. 

  982. 						if (cblk->numpassesinlayers == 0) {
    983. 

  983. 							layer->len = cblk->passes[n - 1].rate;
    984. 

  984. 							layer->data = cblk->data;
    985. 

  985. 							layer->disto = cblk->passes[n - 1].distortiondec;
    986. 

  986. 						} else {
    987. 

  987. 							layer->len = cblk->passes[n - 1].rate -	cblk->passes[cblk->numpassesinlayers - 1].rate;
    988. 

  988. 							layer->data = cblk->data + cblk->passes[cblk->numpassesinlayers - 1].rate;
    989. 

  989. 							layer->disto = cblk->passes[n - 1].distortiondec - cblk->passes[cblk->numpassesinlayers - 1].distortiondec;
    990. 

  990. 						}
    991. 

  991. 						
    992. 

  992. 						tcd_tile->distolayer[layno] += layer->disto;	/* fixed_quality */
    993. 

  993. 						
    994. 

  994. 						if (final)
    995. 

  995. 							cblk->numpassesinlayers = n;
    996. 

  996. 					}
    997. 

  997. 				}
    998. 

  998. 			}
    999. 

  999. 		}
    1000. 

  1000. 	}
    1001. 

  1001. }
    1002. 

  1002. 

  1003. opj_bool tcd_rateallocate(opj_tcd_t *tcd, unsigned char *dest, int len, opj_codestream_info_t *cstr_info) {
    1004. 

  1004. 	int compno, resno, bandno, precno, cblkno, passno, layno;
    1005. 

  1005. 	double min, max;
    1006. 

  1006. 	double cumdisto[100];	/* fixed_quality */
    1007. 

  1007. 	const double K = 1;		/* 1.1; fixed_quality */
    1008. 

  1008. 	double maxSE = 0;
    1009. 

  1009. 

  1010. 	opj_cp_t *cp = tcd->cp;
    1011. 

  1011. 	opj_tcd_tile_t *tcd_tile = tcd->tcd_tile;
    1012. 

  1012. 	opj_tcp_t *tcd_tcp = tcd->tcp;
    1013. 

  1013. 

  1014. 	min = DBL_MAX;
    1015. 

  1015. 	max = 0;
    1016. 

  1016. 	
    1017. 

  1017. 	tcd_tile->numpix = 0;		/* fixed_quality */
    1018. 

  1018. 	
    1019. 

  1019. 	for (compno = 0; compno < tcd_tile->numcomps; compno++) {
    1020. 

  1020. 		opj_tcd_tilecomp_t *tilec = &tcd_tile->comps[compno];
    1021. 

  1021. 		tilec->numpix = 0;
    1022. 

  1022. 

  1023. 		for (resno = 0; resno < tilec->numresolutions; resno++) {
    1024. 

  1024. 			opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    1025. 

  1025. 

  1026. 			for (bandno = 0; bandno < res->numbands; bandno++) {
    1027. 

  1027. 				opj_tcd_band_t *band = &res->bands[bandno];
    1028. 

  1028. 

  1029. 				for (precno = 0; precno < res->pw * res->ph; precno++) {
    1030. 

  1030. 					opj_tcd_precinct_t *prc = &band->precincts[precno];
    1031. 

  1031. 

  1032. 					for (cblkno = 0; cblkno < prc->cw * prc->ch; cblkno++) {
    1033. 

  1033. 						opj_tcd_cblk_enc_t *cblk = &prc->cblks.enc[cblkno];
    1034. 

  1034. 

  1035. 						for (passno = 0; passno < cblk->totalpasses; passno++) {
    1036. 

  1036. 							opj_tcd_pass_t *pass = &cblk->passes[passno];
    1037. 

  1037. 							int dr;
    1038. 

  1038. 							double dd, rdslope;
    1039. 

  1039. 							if (passno == 0) {
    1040. 

  1040. 								dr = pass->rate;
    1041. 

  1041. 								dd = pass->distortiondec;
    1042. 

  1042. 							} else {
    1043. 

  1043. 								dr = pass->rate - cblk->passes[passno - 1].rate;
    1044. 

  1044. 								dd = pass->distortiondec - cblk->passes[passno - 1].distortiondec;
    1045. 

  1045. 							}
    1046. 

  1046. 							if (dr == 0) {
    1047. 

  1047. 								continue;
    1048. 

  1048. 							}
    1049. 

  1049. 							rdslope = dd / dr;
    1050. 

  1050. 							if (rdslope < min) {
    1051. 

  1051. 								min = rdslope;
    1052. 

  1052. 							}
    1053. 

  1053. 							if (rdslope > max) {
    1054. 

  1054. 								max = rdslope;
    1055. 

  1055. 							}
    1056. 

  1056. 						} /* passno */
    1057. 

  1057. 						
    1058. 

  1058. 						/* fixed_quality */
    1059. 

  1059. 						tcd_tile->numpix += ((cblk->x1 - cblk->x0) * (cblk->y1 - cblk->y0));
    1060. 

  1060. 						tilec->numpix += ((cblk->x1 - cblk->x0) * (cblk->y1 - cblk->y0));
    1061. 

  1061. 					} /* cbklno */
    1062. 

  1062. 				} /* precno */
    1063. 

  1063. 			} /* bandno */
    1064. 

  1064. 		} /* resno */
    1065. 

  1065. 		
    1066. 

  1066. 		maxSE += (((double)(1 << tcd->image->comps[compno].prec) - 1.0) 
    1067. 

  1067. 			* ((double)(1 << tcd->image->comps[compno].prec) -1.0)) 
    1068. 

  1068. 			* ((double)(tilec->numpix));
    1069. 

  1069. 	} /* compno */
    1070. 

  1070. 	
    1071. 

  1071. 	/* index file */
    1072. 

  1072. 	if(cstr_info) {
    1073. 

  1073. 		opj_tile_info_t *tile_info = &cstr_info->tile[tcd->tcd_tileno];
    1074. 

  1074. 		tile_info->numpix = tcd_tile->numpix;
    1075. 

  1075. 		tile_info->distotile = tcd_tile->distotile;
    1076. 

  1076. 		tile_info->thresh = (double *) opj_malloc(tcd_tcp->numlayers * sizeof(double));
    1077. 

  1077. 	}
    1078. 

  1078. 	
    1079. 

  1079. 	for (layno = 0; layno < tcd_tcp->numlayers; layno++) {
    1080. 

  1080. 		double lo = min;
    1081. 

  1081. 		double hi = max;
    1082. 

  1082. 		int success = 0;
    1083. 

  1083. 		int maxlen = tcd_tcp->rates[layno] ? int_min(((int) ceil(tcd_tcp->rates[layno])), len) : len;
    1084. 

  1084. 		double goodthresh = 0;
    1085. 

  1085. 		double stable_thresh = 0;
    1086. 

  1086. 		int i;
    1087. 

  1087. 		double distotarget;		/* fixed_quality */
    1088. 

  1088. 		
    1089. 

  1089. 		/* fixed_quality */
    1090. 

  1090. 		distotarget = tcd_tile->distotile - ((K * maxSE) / pow((float)10, tcd_tcp->distoratio[layno] / 10));
    1091. 

  1091.         
    1092. 

  1092. 		/* Don't try to find an optimal threshold but rather take everything not included yet, if
    1093. 

  1093. 		  -r xx,yy,zz,0   (disto_alloc == 1 and rates == 0)
    1094. 

  1094. 		  -q xx,yy,zz,0	  (fixed_quality == 1 and distoratio == 0)
    1095. 

  1095. 		  ==> possible to have some lossy layers and the last layer for sure lossless */
    1096. 

  1096. 		if ( ((cp->disto_alloc==1) && (tcd_tcp->rates[layno]>0)) || ((cp->fixed_quality==1) && (tcd_tcp->distoratio[layno]>0))) {
    1097. 

  1097. 			opj_t2_t *t2 = t2_create(tcd->cinfo, tcd->image, cp);
    1098. 

  1098. 			double thresh = 0;
    1099. 

  1099. 

  1100. 			for (i = 0; i < 128; i++) {
    1101. 

  1101. 				int l = 0;
    1102. 

  1102. 				double distoachieved = 0;	/* fixed_quality */
    1103. 

  1103. 				thresh = (lo + hi) / 2;
    1104. 

  1104. 				
    1105. 

  1105. 				tcd_makelayer(tcd, layno, thresh, 0);
    1106. 

  1106. 				
    1107. 

  1107. 				if (cp->fixed_quality) {	/* fixed_quality */
    1108. 

  1108. 					if(cp->cinema){
    1109. 

  1109. 						l = t2_encode_packets(t2,tcd->tcd_tileno, tcd_tile, layno + 1, dest, maxlen, cstr_info,tcd->cur_tp_num,tcd->tp_pos,tcd->cur_pino,THRESH_CALC, tcd->cur_totnum_tp);
    1110. 

  1110. 						if (l == -999) {
    1111. 

  1111. 							lo = thresh;
    1112. 

  1112. 							continue;
    1113. 

  1113. 						}else{
    1114. 

  1114.            		distoachieved =	layno == 0 ? 
    1115. 

  1115. 							tcd_tile->distolayer[0]	: cumdisto[layno - 1] + tcd_tile->distolayer[layno];
    1116. 

  1116. 							if (distoachieved < distotarget) {
    1117. 

  1117. 								hi=thresh; 
    1118. 

  1118. 								stable_thresh = thresh;
    1119. 

  1119. 								continue;
    1120. 

  1120. 							}else{
    1121. 

  1121. 								lo=thresh;
    1122. 

  1122. 							}
    1123. 

  1123. 						}
    1124. 

  1124. 					}else{
    1125. 

  1125. 						distoachieved =	(layno == 0) ? 
    1126. 

  1126. 							tcd_tile->distolayer[0]	: (cumdisto[layno - 1] + tcd_tile->distolayer[layno]);
    1127. 

  1127. 						if (distoachieved < distotarget) {
    1128. 

  1128. 							hi = thresh;
    1129. 

  1129. 							stable_thresh = thresh;
    1130. 

  1130. 							continue;
    1131. 

  1131. 						}
    1132. 

  1132. 						lo = thresh;
    1133. 

  1133. 					}
    1134. 

  1134. 				} else {
    1135. 

  1135. 					l = t2_encode_packets(t2, tcd->tcd_tileno, tcd_tile, layno + 1, dest, maxlen, cstr_info,tcd->cur_tp_num,tcd->tp_pos,tcd->cur_pino,THRESH_CALC, tcd->cur_totnum_tp);
    1136. 

  1136. 					/* TODO: what to do with l ??? seek / tell ??? */
    1137. 

  1137. 					/* opj_event_msg(tcd->cinfo, EVT_INFO, "rate alloc: len=%d, max=%d\n", l, maxlen); */
    1138. 

  1138. 					if (l == -999) {
    1139. 

  1139. 						lo = thresh;
    1140. 

  1140. 						continue;
    1141. 

  1141. 					}
    1142. 

  1142. 					hi = thresh;
    1143. 

  1143. 					stable_thresh = thresh;
    1144. 

  1144. 				}
    1145. 

  1145. 			}
    1146. 

  1146. 			success = 1;
    1147. 

  1147. 			goodthresh = stable_thresh == 0? thresh : stable_thresh;
    1148. 

  1148. 			t2_destroy(t2);
    1149. 

  1149. 		} else {
    1150. 

  1150. 			success = 1;
    1151. 

  1151. 			goodthresh = min;
    1152. 

  1152. 		}
    1153. 

  1153. 		
    1154. 

  1154. 		if (!success) {
    1155. 

  1155. 			return OPJ_FALSE;
    1156. 

  1156. 		}
    1157. 

  1157. 		
    1158. 

  1158. 		if(cstr_info) {	/* Threshold for Marcela Index */
    1159. 

  1159. 			cstr_info->tile[tcd->tcd_tileno].thresh[layno] = goodthresh;
    1160. 

  1160. 		}
    1161. 

  1161. 		tcd_makelayer(tcd, layno, goodthresh, 1);
    1162. 

  1162.         
    1163. 

  1163. 		/* fixed_quality */
    1164. 

  1164. 		cumdisto[layno] = (layno == 0) ? tcd_tile->distolayer[0] : (cumdisto[layno - 1] + tcd_tile->distolayer[layno]);	
    1165. 

  1165. 	}
    1166. 

  1166. 

  1167. 	return OPJ_TRUE;
    1168. 

  1168. }
    1169. 

  1169. 

  1170. int tcd_encode_tile(opj_tcd_t *tcd, int tileno, unsigned char *dest, int len, opj_codestream_info_t *cstr_info) {
    1171. 

  1171. 	int compno;
    1172. 

  1172. 	int l, i, numpacks = 0;
    1173. 

  1173. 	opj_tcd_tile_t *tile = NULL;
    1174. 

  1174. 	opj_tcp_t *tcd_tcp = NULL;
    1175. 

  1175. 	opj_cp_t *cp = NULL;
    1176. 

  1176. 

  1177. 	opj_tcp_t *tcp = &tcd->cp->tcps[0];
    1178. 

  1178. 	opj_tccp_t *tccp = &tcp->tccps[0];
    1179. 

  1179. 	opj_image_t *image = tcd->image;
    1180. 

  1180. 	
    1181. 

  1181. 	opj_t1_t *t1 = NULL;		/* T1 component */
    1182. 

  1182. 	opj_t2_t *t2 = NULL;		/* T2 component */
    1183. 

  1183. 

  1184. 	tcd->tcd_tileno = tileno;
    1185. 

  1185. 	tcd->tcd_tile = tcd->tcd_image->tiles;
    1186. 

  1186. 	tcd->tcp = &tcd->cp->tcps[tileno];
    1187. 

  1187. 

  1188. 	tile = tcd->tcd_tile;
    1189. 

  1189. 	tcd_tcp = tcd->tcp;
    1190. 

  1190. 	cp = tcd->cp;
    1191. 

  1191. 

  1192. 	if(tcd->cur_tp_num == 0){
    1193. 

  1193. 		tcd->encoding_time = opj_clock();	/* time needed to encode a tile */
    1194. 

  1194. 		/* INDEX >> "Precinct_nb_X et Precinct_nb_Y" */
    1195. 

  1195. 		if(cstr_info) {
    1196. 

  1196. 			opj_tcd_tilecomp_t *tilec_idx = &tile->comps[0];	/* based on component 0 */
    1197. 

  1197. 			for (i = 0; i < tilec_idx->numresolutions; i++) {
    1198. 

  1198. 				opj_tcd_resolution_t *res_idx = &tilec_idx->resolutions[i];
    1199. 

  1199. 				
    1200. 

  1200. 				cstr_info->tile[tileno].pw[i] = res_idx->pw;
    1201. 

  1201. 				cstr_info->tile[tileno].ph[i] = res_idx->ph;
    1202. 

  1202. 				
    1203. 

  1203. 				numpacks += res_idx->pw * res_idx->ph;
    1204. 

  1204. 				
    1205. 

  1205. 				cstr_info->tile[tileno].pdx[i] = tccp->prcw[i];
    1206. 

  1206. 				cstr_info->tile[tileno].pdy[i] = tccp->prch[i];
    1207. 

  1207. 			}
    1208. 

  1208. 			cstr_info->tile[tileno].packet = (opj_packet_info_t*) opj_calloc(cstr_info->numcomps * cstr_info->numlayers * numpacks, sizeof(opj_packet_info_t));
    1209. 

  1209. 		}
    1210. 

  1210. 		/* << INDEX */
    1211. 

  1211. 		
    1212. 

  1212. 		/*---------------TILE-------------------*/
    1213. 

  1213. 		
    1214. 

  1214. 		for (compno = 0; compno < tile->numcomps; compno++) {
    1215. 

  1215. 			int x, y;
    1216. 

  1216. 			
    1217. 

  1217. 			int adjust = image->comps[compno].sgnd ? 0 : 1 << (image->comps[compno].prec - 1);
    1218. 

  1218. 			int offset_x = int_ceildiv(image->x0, image->comps[compno].dx);
    1219. 

  1219. 			int offset_y = int_ceildiv(image->y0, image->comps[compno].dy);
    1220. 

  1220. 			
    1221. 

  1221. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    1222. 

  1222. 			int tw = tilec->x1 - tilec->x0;
    1223. 

  1223. 			int w = int_ceildiv(image->x1 - image->x0, image->comps[compno].dx);
    1224. 

  1224. 			
    1225. 

  1225. 			/* extract tile data */
    1226. 

  1226. 			
    1227. 

  1227. 			if (tcd_tcp->tccps[compno].qmfbid == 1) {
    1228. 

  1228. 				for (y = tilec->y0; y < tilec->y1; y++) {
    1229. 

  1229. 					/* start of the src tile scanline */
    1230. 

  1230. 					int *data = &image->comps[compno].data[(tilec->x0 - offset_x) + (y - offset_y) * w];
    1231. 

  1231. 					/* start of the dst tile scanline */
    1232. 

  1232. 					int *tile_data = &tilec->data[(y - tilec->y0) * tw];
    1233. 

  1233. 					for (x = tilec->x0; x < tilec->x1; x++) {
    1234. 

  1234. 						*tile_data++ = *data++ - adjust;
    1235. 

  1235. 					}
    1236. 

  1236. 				}
    1237. 

  1237. 			} else if (tcd_tcp->tccps[compno].qmfbid == 0) {
    1238. 

  1238. 				for (y = tilec->y0; y < tilec->y1; y++) {
    1239. 

  1239. 					/* start of the src tile scanline */
    1240. 

  1240. 					int *data = &image->comps[compno].data[(tilec->x0 - offset_x) + (y - offset_y) * w];
    1241. 

  1241. 					/* start of the dst tile scanline */
    1242. 

  1242. 					int *tile_data = &tilec->data[(y - tilec->y0) * tw];
    1243. 

  1243. 					for (x = tilec->x0; x < tilec->x1; x++) {
    1244. 

  1244. 						*tile_data++ = (*data++ - adjust) << 11;
    1245. 

  1245. 					}
    1246. 

  1246. 					
    1247. 

  1247. 				}
    1248. 

  1248. 			}
    1249. 

  1249. 		}
    1250. 

  1250. 		
    1251. 

  1251. 		/*----------------MCT-------------------*/
    1252. 

  1252. 		if (tcd_tcp->mct) {
    1253. 

  1253. 			int samples = (tile->comps[0].x1 - tile->comps[0].x0) * (tile->comps[0].y1 - tile->comps[0].y0);
    1254. 

  1254. 			if (tcd_tcp->tccps[0].qmfbid == 0) {
    1255. 

  1255. 				mct_encode_real(tile->comps[0].data, tile->comps[1].data, tile->comps[2].data, samples);
    1256. 

  1256. 			} else {
    1257. 

  1257. 				mct_encode(tile->comps[0].data, tile->comps[1].data, tile->comps[2].data, samples);
    1258. 

  1258. 			}
    1259. 

  1259. 		}
    1260. 

  1260. 		
    1261. 

  1261. 		/*----------------DWT---------------------*/
    1262. 

  1262. 		
    1263. 

  1263. 		for (compno = 0; compno < tile->numcomps; compno++) {
    1264. 

  1264. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    1265. 

  1265. 			if (tcd_tcp->tccps[compno].qmfbid == 1) {
    1266. 

  1266. 				dwt_encode(tilec);
    1267. 

  1267. 			} else if (tcd_tcp->tccps[compno].qmfbid == 0) {
    1268. 

  1268. 				dwt_encode_real(tilec);
    1269. 

  1269. 			}
    1270. 

  1270. 		}
    1271. 

  1271. 		
    1272. 

  1272. 		/*------------------TIER1-----------------*/
    1273. 

  1273. 		t1 = t1_create(tcd->cinfo);
    1274. 

  1274. 		t1_encode_cblks(t1, tile, tcd_tcp);
    1275. 

  1275. 		t1_destroy(t1);
    1276. 

  1276. 		
    1277. 

  1277. 		/*-----------RATE-ALLOCATE------------------*/
    1278. 

  1278. 		
    1279. 

  1279. 		/* INDEX */
    1280. 

  1280. 		if(cstr_info) {
    1281. 

  1281. 			cstr_info->index_write = 0;
    1282. 

  1282. 		}
    1283. 

  1283. 		if (cp->disto_alloc || cp->fixed_quality) {	/* fixed_quality */
    1284. 

  1284. 			/* Normal Rate/distortion allocation */
    1285. 

  1285. 			tcd_rateallocate(tcd, dest, len, cstr_info);
    1286. 

  1286. 		} else {
    1287. 

  1287. 			/* Fixed layer allocation */
    1288. 

  1288. 			tcd_rateallocate_fixed(tcd);
    1289. 

  1289. 		}
    1290. 

  1290. 	}
    1291. 

  1291. 	/*--------------TIER2------------------*/
    1292. 

  1292. 

  1293. 	/* INDEX */
    1294. 

  1294. 	if(cstr_info) {
    1295. 

  1295. 		cstr_info->index_write = 1;
    1296. 

  1296. 	}
    1297. 

  1297. 

  1298. 	t2 = t2_create(tcd->cinfo, image, cp);
    1299. 

  1299. 	l = t2_encode_packets(t2,tileno, tile, tcd_tcp->numlayers, dest, len, cstr_info,tcd->tp_num,tcd->tp_pos,tcd->cur_pino,FINAL_PASS,tcd->cur_totnum_tp);
    1300. 

  1300. 	t2_destroy(t2);
    1301. 

  1301. 	
    1302. 

  1302. 	/*---------------CLEAN-------------------*/
    1303. 

  1303. 

  1304. 	
    1305. 

  1305. 	if(tcd->cur_tp_num == tcd->cur_totnum_tp - 1){
    1306. 

  1306. 		tcd->encoding_time = opj_clock() - tcd->encoding_time;
    1307. 

  1307. 		opj_event_msg(tcd->cinfo, EVT_INFO, "- tile encoded in %f s\n", tcd->encoding_time);
    1308. 

  1308. 

  1309. 		/* cleaning memory */
    1310. 

  1310. 		for (compno = 0; compno < tile->numcomps; compno++) {
    1311. 

  1311. 			opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    1312. 

  1312. 			opj_aligned_free(tilec->data);
    1313. 

  1313. 		}
    1314. 

  1314. 	}
    1315. 

  1315. 

  1316. 	return l;
    1317. 

  1317. }
    1318. 

  1318. 

  1319. opj_bool tcd_decode_tile(opj_tcd_t *tcd, unsigned char *src, int len, int tileno, opj_codestream_info_t *cstr_info) {
    1320. 

  1320. 	int l;
    1321. 

  1321. 	int compno;
    1322. 

  1322. 	int eof = 0;
    1323. 

  1323. 	double tile_time, t1_time, dwt_time;
    1324. 

  1324. 	opj_tcd_tile_t *tile = NULL;
    1325. 

  1325. 

  1326. 	opj_t1_t *t1 = NULL;		/* T1 component */
    1327. 

  1327. 	opj_t2_t *t2 = NULL;		/* T2 component */
    1328. 

  1328. 	
    1329. 

  1329. 	tcd->tcd_tileno = tileno;
    1330. 

  1330. 	tcd->tcd_tile = &(tcd->tcd_image->tiles[tileno]);
    1331. 

  1331. 	tcd->tcp = &(tcd->cp->tcps[tileno]);
    1332. 

  1332. 	tile = tcd->tcd_tile;
    1333. 

  1333. 	
    1334. 

  1334. 	tile_time = opj_clock();	/* time needed to decode a tile */
    1335. 

  1335. 	opj_event_msg(tcd->cinfo, EVT_INFO, "tile %d of %d\n", tileno + 1, tcd->cp->tw * tcd->cp->th);
    1336. 

  1336. 

  1337. 	/* INDEX >>  */
    1338. 

  1338. 	if(cstr_info) {
    1339. 

  1339. 		int resno, compno, numprec = 0;
    1340. 

  1340. 		for (compno = 0; compno < cstr_info->numcomps; compno++) {
    1341. 

  1341. 			opj_tcp_t *tcp = &tcd->cp->tcps[0];
    1342. 

  1342. 			opj_tccp_t *tccp = &tcp->tccps[compno];
    1343. 

  1343. 			opj_tcd_tilecomp_t *tilec_idx = &tile->comps[compno];	
    1344. 

  1344. 			for (resno = 0; resno < tilec_idx->numresolutions; resno++) {
    1345. 

  1345. 				opj_tcd_resolution_t *res_idx = &tilec_idx->resolutions[resno];
    1346. 

  1346. 				cstr_info->tile[tileno].pw[resno] = res_idx->pw;
    1347. 

  1347. 				cstr_info->tile[tileno].ph[resno] = res_idx->ph;
    1348. 

  1348. 				numprec += res_idx->pw * res_idx->ph;
    1349. 

  1349. 				if (tccp->csty & J2K_CP_CSTY_PRT) {
    1350. 

  1350. 					cstr_info->tile[tileno].pdx[resno] = tccp->prcw[resno];
    1351. 

  1351. 					cstr_info->tile[tileno].pdy[resno] = tccp->prch[resno];
    1352. 

  1352. 				}
    1353. 

  1353. 				else {
    1354. 

  1354. 					cstr_info->tile[tileno].pdx[resno] = 15;
    1355. 

  1355. 					cstr_info->tile[tileno].pdy[resno] = 15;
    1356. 

  1356. 				}
    1357. 

  1357. 			}
    1358. 

  1358. 		}
    1359. 

  1359. 		cstr_info->tile[tileno].packet = (opj_packet_info_t *) opj_malloc(cstr_info->numlayers * numprec * sizeof(opj_packet_info_t));
    1360. 

  1360. 		cstr_info->packno = 0;
    1361. 

  1361. 	}
    1362. 

  1362. 	/* << INDEX */
    1363. 

  1363. 	
    1364. 

  1364. 	/*--------------TIER2------------------*/
    1365. 

  1365. 	
    1366. 

  1366. 	t2 = t2_create(tcd->cinfo, tcd->image, tcd->cp);
    1367. 

  1367. 	l = t2_decode_packets(t2, src, len, tileno, tile, cstr_info);
    1368. 

  1368. 	t2_destroy(t2);
    1369. 

  1369. 

  1370. 	if (l == -999) {
    1371. 

  1371. 		eof = 1;
    1372. 

  1372. 		opj_event_msg(tcd->cinfo, EVT_ERROR, "tcd_decode: incomplete bistream\n");
    1373. 

  1373. 	}
    1374. 

  1374. 	
    1375. 

  1375. 	/*------------------TIER1-----------------*/
    1376. 

  1376. 	
    1377. 

  1377. 	t1_time = opj_clock();	/* time needed to decode a tile */
    1378. 

  1378. 	t1 = t1_create(tcd->cinfo);
    1379. 

  1379. 	for (compno = 0; compno < tile->numcomps; ++compno) {
    1380. 

  1380. 		opj_tcd_tilecomp_t* tilec = &tile->comps[compno];
    1381. 

  1381. 		/* The +3 is headroom required by the vectorized DWT */
    1382. 

  1382. 		tilec->data = (int*) opj_aligned_malloc((((tilec->x1 - tilec->x0) * (tilec->y1 - tilec->y0))+3) * sizeof(int));
    1383. 

  1383. 		t1_decode_cblks(t1, tilec, &tcd->tcp->tccps[compno]);
    1384. 

  1384. 	}
    1385. 

  1385. 	t1_destroy(t1);
    1386. 

  1386. 	t1_time = opj_clock() - t1_time;
    1387. 

  1387. 	opj_event_msg(tcd->cinfo, EVT_INFO, "- tiers-1 took %f s\n", t1_time);
    1388. 

  1388. 	
    1389. 

  1389. 	/*----------------DWT---------------------*/
    1390. 

  1390. 

  1391. 	dwt_time = opj_clock();	/* time needed to decode a tile */
    1392. 

  1392. 	for (compno = 0; compno < tile->numcomps; compno++) {
    1393. 

  1393. 		opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    1394. 

  1394. 		int numres2decode;
    1395. 

  1395. 

  1396. 		if (tcd->cp->reduce != 0) {
    1397. 

  1397. 			tcd->image->comps[compno].resno_decoded =
    1398. 

  1398. 				tile->comps[compno].numresolutions - tcd->cp->reduce - 1;
    1399. 

  1399. 			if (tcd->image->comps[compno].resno_decoded < 0) {				
    1400. 

  1400. 				opj_event_msg(tcd->cinfo, EVT_ERROR, "Error decoding tile. The number of resolutions to remove [%d+1] is higher than the number "
    1401. 

  1401. 					" of resolutions in the original codestream [%d]\nModify the cp_reduce parameter.\n", tcd->cp->reduce, tile->comps[compno].numresolutions);
    1402. 

  1402. 				return OPJ_FALSE;
    1403. 

  1403. 			}
    1404. 

  1404. 		}
    1405. 

  1405. 

  1406. 		numres2decode = tcd->image->comps[compno].resno_decoded + 1;
    1407. 

  1407. 		if(numres2decode > 0){
    1408. 

  1408. 			if (tcd->tcp->tccps[compno].qmfbid == 1) {
    1409. 

  1409. 				dwt_decode(tilec, numres2decode);
    1410. 

  1410. 			} else {
    1411. 

  1411. 				dwt_decode_real(tilec, numres2decode);
    1412. 

  1412. 			}
    1413. 

  1413. 		}
    1414. 

  1414. 	}
    1415. 

  1415. 	dwt_time = opj_clock() - dwt_time;
    1416. 

  1416. 	opj_event_msg(tcd->cinfo, EVT_INFO, "- dwt took %f s\n", dwt_time);
    1417. 

  1417. 

  1418. 	/*----------------MCT-------------------*/
    1419. 

  1419. 

  1420. 	if (tcd->tcp->mct) {
    1421. 

  1421. 		int n = (tile->comps[0].x1 - tile->comps[0].x0) * (tile->comps[0].y1 - tile->comps[0].y0);
    1422. 

  1422. 

  1423. 		if (tile->numcomps >= 3 ){
    1424. 

  1424. 			if (tcd->tcp->tccps[0].qmfbid == 1) {
    1425. 

  1425. 				mct_decode(
    1426. 

  1426. 						tile->comps[0].data,
    1427. 

  1427. 						tile->comps[1].data,
    1428. 

  1428. 						tile->comps[2].data,
    1429. 

  1429. 						n);
    1430. 

  1430. 			} else {
    1431. 

  1431. 				mct_decode_real(
    1432. 

  1432. 						(float*)tile->comps[0].data,
    1433. 

  1433. 						(float*)tile->comps[1].data,
    1434. 

  1434. 						(float*)tile->comps[2].data,
    1435. 

  1435. 						n);
    1436. 

  1436. 			}
    1437. 

  1437. 		} else{
    1438. 

  1438. 			opj_event_msg(tcd->cinfo, EVT_WARNING,"Number of components (%d) is inconsistent with a MCT. Skip the MCT step.\n",tile->numcomps);
    1439. 

  1439. 		}
    1440. 

  1440. 	}
    1441. 

  1441. 

  1442. 	/*---------------TILE-------------------*/
    1443. 

  1443. 

  1444. 	for (compno = 0; compno < tile->numcomps; ++compno) {
    1445. 

  1445. 		opj_tcd_tilecomp_t* tilec = &tile->comps[compno];
    1446. 

  1446. 		opj_image_comp_t* imagec = &tcd->image->comps[compno];
    1447. 

  1447. 		opj_tcd_resolution_t* res = &tilec->resolutions[imagec->resno_decoded];
    1448. 

  1448. 		int adjust = imagec->sgnd ? 0 : 1 << (imagec->prec - 1);
    1449. 

  1449. 		int min = imagec->sgnd ? -(1 << (imagec->prec - 1)) : 0;
    1450. 

  1450. 		int max = imagec->sgnd ?  (1 << (imagec->prec - 1)) - 1 : (1 << imagec->prec) - 1;
    1451. 

  1451. 

  1452. 		int tw = tilec->x1 - tilec->x0;
    1453. 

  1453. 		int w = imagec->w;
    1454. 

  1454. 

  1455. 		int offset_x = int_ceildivpow2(imagec->x0, imagec->factor);
    1456. 

  1456. 		int offset_y = int_ceildivpow2(imagec->y0, imagec->factor);
    1457. 

  1457. 

  1458. 		int i, j;
    1459. 

  1459. 		if(!imagec->data){
    1460. 

  1460. 			imagec->data = (int*) opj_malloc(imagec->w * imagec->h * sizeof(int));
    1461. 

  1461. 		}
    1462. 

  1462. 		if(tcd->tcp->tccps[compno].qmfbid == 1) {
    1463. 

  1463. 			for(j = res->y0; j < res->y1; ++j) {
    1464. 

  1464. 				for(i = res->x0; i < res->x1; ++i) {
    1465. 

  1465. 					int v = tilec->data[i - res->x0 + (j - res->y0) * tw];
    1466. 

  1466. 					v += adjust;
    1467. 

  1467. 					imagec->data[(i - offset_x) + (j - offset_y) * w] = int_clamp(v, min, max);
    1468. 

  1468. 				}
    1469. 

  1469. 			}
    1470. 

  1470. 		}else{
    1471. 

  1471. 			for(j = res->y0; j < res->y1; ++j) {
    1472. 

  1472. 				for(i = res->x0; i < res->x1; ++i) {
    1473. 

  1473. 					float tmp = ((float*)tilec->data)[i - res->x0 + (j - res->y0) * tw];
    1474. 

  1474. 					int v = lrintf(tmp);
    1475. 

  1475. 					v += adjust;
    1476. 

  1476. 					imagec->data[(i - offset_x) + (j - offset_y) * w] = int_clamp(v, min, max);
    1477. 

  1477. 				}
    1478. 

  1478. 			}
    1479. 

  1479. 		}
    1480. 

  1480. 		opj_aligned_free(tilec->data);
    1481. 

  1481. 	}
    1482. 

  1482. 

  1483. 	tile_time = opj_clock() - tile_time;	/* time needed to decode a tile */
    1484. 

  1484. 	opj_event_msg(tcd->cinfo, EVT_INFO, "- tile decoded in %f s\n", tile_time);
    1485. 

  1485. 

  1486. 	if (eof) {
    1487. 

  1487. 		return OPJ_FALSE;
    1488. 

  1488. 	}
    1489. 

  1489. 	
    1490. 

  1490. 	return OPJ_TRUE;
    1491. 

  1491. }
    1492. 

  1492. 

  1493. void tcd_free_decode(opj_tcd_t *tcd) {
    1494. 

  1494. 	opj_tcd_image_t *tcd_image = tcd->tcd_image;	
    1495. 

  1495. 	opj_free(tcd_image->tiles);
    1496. 

  1496. }
    1497. 

  1497. 

  1498. void tcd_free_decode_tile(opj_tcd_t *tcd, int tileno) {
    1499. 

  1499. 	int compno,resno,bandno,precno;
    1500. 

  1500. 

  1501. 	opj_tcd_image_t *tcd_image = tcd->tcd_image;
    1502. 

  1502. 

  1503. 	opj_tcd_tile_t *tile = &tcd_image->tiles[tileno];
    1504. 

  1504. 	for (compno = 0; compno < tile->numcomps; compno++) {
    1505. 

  1505. 		opj_tcd_tilecomp_t *tilec = &tile->comps[compno];
    1506. 

  1506. 		for (resno = 0; resno < tilec->numresolutions; resno++) {
    1507. 

  1507. 			opj_tcd_resolution_t *res = &tilec->resolutions[resno];
    1508. 

  1508. 			for (bandno = 0; bandno < res->numbands; bandno++) {
    1509. 

  1509. 				opj_tcd_band_t *band = &res->bands[bandno];
    1510. 

  1510. 				for (precno = 0; precno < res->ph * res->pw; precno++) {
    1511. 

  1511. 					opj_tcd_precinct_t *prec = &band->precincts[precno];
    1512. 

  1512. 					if (prec->imsbtree != NULL) tgt_destroy(prec->imsbtree);
    1513. 

  1513. 					if (prec->incltree != NULL) tgt_destroy(prec->incltree);
    1514. 

  1514. 				}
    1515. 

  1515. 				opj_free(band->precincts);
    1516. 

  1516. 			}
    1517. 

  1517. 		}
    1518. 

  1518. 		opj_free(tilec->resolutions);
    1519. 

  1519. 	}
    1520. 

  1520. 	opj_free(tile->comps);
    1521. 

  1521. }
    1522. 

  1522. 

  1523. 

  1524. 

  1525.