2 * Copyright (c) 2002-2007, Communications and Remote Sensing Laboratory, Universite catholique de Louvain (UCL), Belgium
3 * Copyright (c) 2002-2007, Professor Benoit Macq
4 * Copyright (c) 2001-2003, David Janssens
5 * Copyright (c) 2002-2003, Yannick Verschueren
6 * Copyright (c) 2003-2007, Francois-Olivier Devaux and Antonin Descampe
7 * Copyright (c) 2005, Herve Drolon, FreeImage Team
8 * Copyright (c) 2010-2011, Kaori Hagihara
11 * Redistribution and use in source and binary forms, with or without
12 * modification, are permitted provided that the following conditions
14 * 1. Redistributions of source code must retain the above copyright
15 * notice, this list of conditions and the following disclaimer.
16 * 2. Redistributions in binary form must reproduce the above copyright
17 * notice, this list of conditions and the following disclaimer in the
18 * documentation and/or other materials provided with the distribution.
20 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
21 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
22 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
23 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
24 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
25 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
26 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
27 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
28 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
29 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
30 * POSSIBILITY OF SUCH DAMAGE.
32 #include "opj_includes.h"
34 /** @defgroup JP2 JP2 - JPEG-2000 file format reader/writer */
37 /** @name Local static functions */
42 @param cinfo Codec context info
43 @param cio Input stream
45 @return Returns true if successful, returns false otherwise
47 static opj_bool jp2_read_boxhdr(opj_common_ptr cinfo, opj_cio_t *cio, opj_jp2_box_t *box);
48 /*static void jp2_write_url(opj_cio_t *cio, char *Idx_file);*/
50 Read the IHDR box - Image Header box
52 @param cio Input buffer stream
53 @return Returns true if successful, returns false otherwise
55 static opj_bool jp2_read_ihdr(opj_jp2_t *jp2, opj_cio_t *cio);
56 static void jp2_write_ihdr(opj_jp2_t *jp2, opj_cio_t *cio);
57 static void jp2_write_bpcc(opj_jp2_t *jp2, opj_cio_t *cio);
58 static opj_bool jp2_read_bpcc(opj_jp2_t *jp2, opj_cio_t *cio);
59 static void jp2_write_colr(opj_jp2_t *jp2, opj_cio_t *cio);
61 Write the FTYP box - File type box
63 @param cio Output buffer stream
65 static void jp2_write_ftyp(opj_jp2_t *jp2, opj_cio_t *cio);
67 Read the FTYP box - File type box
69 @param cio Input buffer stream
70 @return Returns true if successful, returns false otherwise
72 static opj_bool jp2_read_ftyp(opj_jp2_t *jp2, opj_cio_t *cio);
73 static int jp2_write_jp2c(opj_jp2_t *jp2, opj_cio_t *cio, opj_image_t *image, opj_codestream_info_t *cstr_info);
74 static opj_bool jp2_read_jp2c(opj_jp2_t *jp2, opj_cio_t *cio, unsigned int *j2k_codestream_length, unsigned int *j2k_codestream_offset);
75 static void jp2_write_jp(opj_cio_t *cio);
77 Read the JP box - JPEG 2000 signature
79 @param cio Input buffer stream
80 @return Returns true if successful, returns false otherwise
82 static opj_bool jp2_read_jp(opj_jp2_t *jp2, opj_cio_t *cio);
84 Decode the structure of a JP2 file
86 @param cio Input buffer stream
87 @param color Collector for profile, cdef and pclr data
88 @return Returns true if successful, returns false otherwise
90 static opj_bool jp2_read_struct(opj_jp2_t *jp2, opj_cio_t *cio,
91 opj_jp2_color_t *color);
93 Apply collected palette data
94 @param color Collector for profile, cdef and pclr data
97 static void jp2_apply_pclr(opj_jp2_color_t *color, opj_image_t *image);
100 @param jp2 JP2 handle
101 @param cio Input buffer stream
103 @param color Collector for profile, cdef and pclr data
104 @return Returns true if successful, returns false otherwise
106 static opj_bool jp2_read_pclr(opj_jp2_t *jp2, opj_cio_t *cio,
107 opj_jp2_box_t *box, opj_jp2_color_t *color);
109 Collect component mapping data
110 @param jp2 JP2 handle
111 @param cio Input buffer stream
113 @param color Collector for profile, cdef and pclr data
114 @return Returns true if successful, returns false otherwise
116 static opj_bool jp2_read_cmap(opj_jp2_t *jp2, opj_cio_t *cio,
117 opj_jp2_box_t *box, opj_jp2_color_t *color);
119 Collect colour specification data
120 @param jp2 JP2 handle
121 @param cio Input buffer stream
123 @param color Collector for profile, cdef and pclr data
124 @return Returns true if successful, returns false otherwise
126 static opj_bool jp2_read_colr(opj_jp2_t *jp2, opj_cio_t *cio,
127 opj_jp2_box_t *box, opj_jp2_color_t *color);
129 Write file Index (superbox)
130 @param[in] offset_jp2c offset of jp2c box
131 @param[in] length_jp2c length of jp2c box
132 @param[in] offset_idx offset of cidx box
133 @param[in] length_idx length of cidx box
134 @param[in] cio file output handle
135 @return length of fidx box
137 static int write_fidx( int offset_jp2c, int length_jp2c, int offset_idx, int length_idx, opj_cio_t *cio);
139 Write index Finder box
140 @param[in] offset offset of fidx box
141 @param[in] length length of fidx box
142 @param[in] cio file output handle
144 static void write_iptr( int offset, int length, opj_cio_t *cio);
147 @param[in] offset_jp2c offset of jp2c box
148 @param[in] length_jp2c length of jp2c box
149 @param[in] offset_idx offset of cidx box
150 @param[in] length_idx length of cidx box
151 @param[in] cio file output handle
153 static void write_prxy( int offset_jp2c, int length_jp2c, int offset_idx, int length_idx, opj_cio_t *cio);
158 /* ----------------------------------------------------------------------- */
160 static opj_bool jp2_read_boxhdr(opj_common_ptr cinfo, opj_cio_t *cio, opj_jp2_box_t *box) {
161 box->init_pos = cio_tell(cio);
162 box->length = cio_read(cio, 4);
163 box->type = cio_read(cio, 4);
164 if (box->length == 1) {
165 if (cio_read(cio, 4) != 0) {
166 opj_event_msg(cinfo, EVT_ERROR, "Cannot handle box sizes higher than 2^32\n");
169 box->length = cio_read(cio, 4);
170 if (box->length == 0)
171 box->length = cio_numbytesleft(cio) + 12;
173 else if (box->length == 0) {
174 box->length = cio_numbytesleft(cio) + 8;
181 static void jp2_write_url(opj_cio_t *cio, char *Idx_file) {
185 box.init_pos = cio_tell(cio);
187 cio_write(cio, JP2_URL, 4); /* DBTL */
188 cio_write(cio, 0, 1); /* VERS */
189 cio_write(cio, 0, 3); /* FLAG */
192 for (i = 0; i < strlen(Idx_file); i++) {
193 cio_write(cio, Idx_file[i], 1);
197 box.length = cio_tell(cio) - box.init_pos;
198 cio_seek(cio, box.init_pos);
199 cio_write(cio, box.length, 4); /* L */
200 cio_seek(cio, box.init_pos + box.length);
204 static opj_bool jp2_read_ihdr(opj_jp2_t *jp2, opj_cio_t *cio) {
207 opj_common_ptr cinfo = jp2->cinfo;
209 jp2_read_boxhdr(cinfo, cio, &box);
210 if (JP2_IHDR != box.type) {
211 opj_event_msg(cinfo, EVT_ERROR, "Expected IHDR Marker\n");
215 jp2->h = cio_read(cio, 4); /* HEIGHT */
216 jp2->w = cio_read(cio, 4); /* WIDTH */
217 jp2->numcomps = cio_read(cio, 2); /* NC */
218 jp2->comps = (opj_jp2_comps_t*) opj_malloc(jp2->numcomps * sizeof(opj_jp2_comps_t));
220 jp2->bpc = cio_read(cio, 1); /* BPC */
222 jp2->C = cio_read(cio, 1); /* C */
223 jp2->UnkC = cio_read(cio, 1); /* UnkC */
224 jp2->IPR = cio_read(cio, 1); /* IPR */
226 if (cio_tell(cio) - box.init_pos != box.length) {
227 opj_event_msg(cinfo, EVT_ERROR, "Error with IHDR Box\n");
234 static void jp2_write_ihdr(opj_jp2_t *jp2, opj_cio_t *cio) {
237 box.init_pos = cio_tell(cio);
239 cio_write(cio, JP2_IHDR, 4); /* IHDR */
241 cio_write(cio, jp2->h, 4); /* HEIGHT */
242 cio_write(cio, jp2->w, 4); /* WIDTH */
243 cio_write(cio, jp2->numcomps, 2); /* NC */
245 cio_write(cio, jp2->bpc, 1); /* BPC */
247 cio_write(cio, jp2->C, 1); /* C : Always 7 */
248 cio_write(cio, jp2->UnkC, 1); /* UnkC, colorspace unknown */
249 cio_write(cio, jp2->IPR, 1); /* IPR, no intellectual property */
251 box.length = cio_tell(cio) - box.init_pos;
252 cio_seek(cio, box.init_pos);
253 cio_write(cio, box.length, 4); /* L */
254 cio_seek(cio, box.init_pos + box.length);
257 static void jp2_write_bpcc(opj_jp2_t *jp2, opj_cio_t *cio) {
261 box.init_pos = cio_tell(cio);
263 cio_write(cio, JP2_BPCC, 4); /* BPCC */
265 for (i = 0; i < jp2->numcomps; i++) {
266 cio_write(cio, jp2->comps[i].bpcc, 1);
269 box.length = cio_tell(cio) - box.init_pos;
270 cio_seek(cio, box.init_pos);
271 cio_write(cio, box.length, 4); /* L */
272 cio_seek(cio, box.init_pos + box.length);
276 static opj_bool jp2_read_bpcc(opj_jp2_t *jp2, opj_cio_t *cio) {
280 opj_common_ptr cinfo = jp2->cinfo;
282 jp2_read_boxhdr(cinfo, cio, &box);
283 if (JP2_BPCC != box.type) {
284 opj_event_msg(cinfo, EVT_ERROR, "Expected BPCC Marker\n");
288 for (i = 0; i < jp2->numcomps; i++) {
289 jp2->comps[i].bpcc = cio_read(cio, 1);
292 if (cio_tell(cio) - box.init_pos != box.length) {
293 opj_event_msg(cinfo, EVT_ERROR, "Error with BPCC Box\n");
300 static void jp2_write_colr(opj_jp2_t *jp2, opj_cio_t *cio) {
303 box.init_pos = cio_tell(cio);
305 cio_write(cio, JP2_COLR, 4); /* COLR */
307 cio_write(cio, jp2->meth, 1); /* METH */
308 cio_write(cio, jp2->precedence, 1); /* PRECEDENCE */
309 cio_write(cio, jp2->approx, 1); /* APPROX */
314 cio_write(cio, jp2->enumcs, 4); /* EnumCS */
316 box.length = cio_tell(cio) - box.init_pos;
317 cio_seek(cio, box.init_pos);
318 cio_write(cio, box.length, 4); /* L */
319 cio_seek(cio, box.init_pos + box.length);
322 static void jp2_free_pclr(opj_jp2_color_t *color)
324 opj_free(color->jp2_pclr->channel_sign);
325 opj_free(color->jp2_pclr->channel_size);
326 opj_free(color->jp2_pclr->entries);
328 if(color->jp2_pclr->cmap) opj_free(color->jp2_pclr->cmap);
330 opj_free(color->jp2_pclr); color->jp2_pclr = NULL;
333 static void free_color_data(opj_jp2_color_t *color)
337 jp2_free_pclr(color);
341 if(color->jp2_cdef->info) opj_free(color->jp2_cdef->info);
342 opj_free(color->jp2_cdef);
344 if(color->icc_profile_buf) opj_free(color->icc_profile_buf);
347 static void jp2_apply_pclr(opj_jp2_color_t *color, opj_image_t *image)
349 opj_image_comp_t *old_comps, *new_comps;
350 unsigned char *channel_size, *channel_sign;
351 unsigned int *entries;
352 opj_jp2_cmap_comp_t *cmap;
355 unsigned short i, nr_channels, cmp, pcol;
358 channel_size = color->jp2_pclr->channel_size;
359 channel_sign = color->jp2_pclr->channel_sign;
360 entries = color->jp2_pclr->entries;
361 cmap = color->jp2_pclr->cmap;
362 nr_channels = color->jp2_pclr->nr_channels;
364 old_comps = image->comps;
365 new_comps = (opj_image_comp_t*)
366 opj_malloc(nr_channels * sizeof(opj_image_comp_t));
368 for(i = 0; i < nr_channels; ++i)
370 pcol = cmap[i].pcol; cmp = cmap[i].cmp;
372 new_comps[pcol] = old_comps[cmp];
374 if(cmap[i].mtyp == 0) /* Direct use */
376 old_comps[cmp].data = NULL; continue;
378 /* Palette mapping: */
379 new_comps[pcol].data = (int*)
380 opj_malloc(old_comps[cmp].w * old_comps[cmp].h * sizeof(int));
381 new_comps[pcol].prec = channel_size[i];
382 new_comps[pcol].sgnd = channel_sign[i];
384 top_k = color->jp2_pclr->nr_entries - 1;
386 for(i = 0; i < nr_channels; ++i)
389 if(cmap[i].mtyp == 0) continue;
391 /* Palette mapping: */
392 cmp = cmap[i].cmp; pcol = cmap[i].pcol;
393 src = old_comps[cmp].data;
394 dst = new_comps[pcol].data;
395 max = new_comps[pcol].w * new_comps[pcol].h;
397 for(j = 0; j < max; ++j)
400 if((k = src[j]) < 0) k = 0; else if(k > top_k) k = top_k;
402 dst[j] = entries[k * nr_channels + pcol];
405 max = image->numcomps;
406 for(i = 0; i < max; ++i)
408 if(old_comps[i].data) opj_free(old_comps[i].data);
411 image->comps = new_comps;
412 image->numcomps = nr_channels;
414 jp2_free_pclr(color);
418 static opj_bool jp2_read_pclr(opj_jp2_t *jp2, opj_cio_t *cio,
419 opj_jp2_box_t *box, opj_jp2_color_t *color)
421 opj_jp2_pclr_t *jp2_pclr;
422 unsigned char *channel_size, *channel_sign;
423 unsigned int *entries;
424 unsigned short nr_entries, nr_channels;
428 OPJ_ARG_NOT_USED(box);
429 OPJ_ARG_NOT_USED(jp2);
431 /* Part 1, I.5.3.4: 'There shall be at most one Palette box inside
432 * a JP2 Header box' :
434 if(color->jp2_pclr) return OPJ_FALSE;
436 nr_entries = (unsigned short)cio_read(cio, 2); /* NE */
437 nr_channels = (unsigned short)cio_read(cio, 1);/* NPC */
439 entries = (unsigned int*)
440 opj_malloc(nr_channels * nr_entries * sizeof(unsigned int));
441 channel_size = (unsigned char*)opj_malloc(nr_channels);
442 channel_sign = (unsigned char*)opj_malloc(nr_channels);
444 jp2_pclr = (opj_jp2_pclr_t*)opj_malloc(sizeof(opj_jp2_pclr_t));
445 jp2_pclr->channel_sign = channel_sign;
446 jp2_pclr->channel_size = channel_size;
447 jp2_pclr->entries = entries;
448 jp2_pclr->nr_entries = nr_entries;
449 jp2_pclr->nr_channels = nr_channels;
450 jp2_pclr->cmap = NULL;
452 color->jp2_pclr = jp2_pclr;
454 for(i = 0; i < nr_channels; ++i)
456 uc = cio_read(cio, 1); /* Bi */
457 channel_size[i] = (uc & 0x7f) + 1;
458 channel_sign[i] = (uc & 0x80)?1:0;
461 for(j = 0; j < nr_entries; ++j)
463 for(i = 0; i < nr_channels; ++i)
466 *entries++ = cio_read(cio, channel_size[i]>>3);
471 }/* jp2_read_pclr() */
473 static opj_bool jp2_read_cmap(opj_jp2_t *jp2, opj_cio_t *cio,
474 opj_jp2_box_t *box, opj_jp2_color_t *color)
476 opj_jp2_cmap_comp_t *cmap;
477 unsigned short i, nr_channels;
479 OPJ_ARG_NOT_USED(box);
480 OPJ_ARG_NOT_USED(jp2);
482 /* Need nr_channels: */
483 if(color->jp2_pclr == NULL) return OPJ_FALSE;
485 /* Part 1, I.5.3.5: 'There shall be at most one Component Mapping box
486 * inside a JP2 Header box' :
488 if(color->jp2_pclr->cmap) return OPJ_FALSE;
490 nr_channels = color->jp2_pclr->nr_channels;
491 cmap = (opj_jp2_cmap_comp_t*)
492 opj_malloc(nr_channels * sizeof(opj_jp2_cmap_comp_t));
494 for(i = 0; i < nr_channels; ++i)
496 cmap[i].cmp = (unsigned short)cio_read(cio, 2);
497 cmap[i].mtyp = cio_read(cio, 1);
498 cmap[i].pcol = cio_read(cio, 1);
501 color->jp2_pclr->cmap = cmap;
504 }/* jp2_read_cmap() */
506 static void jp2_apply_cdef(opj_image_t *image, opj_jp2_color_t *color)
508 opj_jp2_cdef_info_t *info;
510 unsigned short i, n, cn, typ, asoc, acn;
512 color_space = image->color_space;
513 info = color->jp2_cdef->info;
514 n = color->jp2_cdef->n;
516 for(i = 0; i < n; ++i)
518 /* WATCH: acn = asoc - 1 ! */
519 if((asoc = info[i].asoc) == 0) continue;
521 cn = info[i].cn; typ = info[i].typ; acn = asoc - 1;
525 opj_image_comp_t saved;
527 memcpy(&saved, &image->comps[cn], sizeof(opj_image_comp_t));
528 memcpy(&image->comps[cn], &image->comps[acn], sizeof(opj_image_comp_t));
529 memcpy(&image->comps[acn], &saved, sizeof(opj_image_comp_t));
531 info[i].asoc = cn + 1;
532 info[acn].asoc = info[acn].cn + 1;
535 if(color->jp2_cdef->info) opj_free(color->jp2_cdef->info);
537 opj_free(color->jp2_cdef); color->jp2_cdef = NULL;
539 }/* jp2_apply_cdef() */
541 static opj_bool jp2_read_cdef(opj_jp2_t *jp2, opj_cio_t *cio,
542 opj_jp2_box_t *box, opj_jp2_color_t *color)
544 opj_jp2_cdef_info_t *info;
547 OPJ_ARG_NOT_USED(box);
548 OPJ_ARG_NOT_USED(jp2);
550 /* Part 1, I.5.3.6: 'The shall be at most one Channel Definition box
551 * inside a JP2 Header box.'
553 if(color->jp2_cdef) return OPJ_FALSE;
555 if((n = (unsigned short)cio_read(cio, 2)) == 0) return OPJ_FALSE; /* szukw000: FIXME */
557 info = (opj_jp2_cdef_info_t*)
558 opj_malloc(n * sizeof(opj_jp2_cdef_info_t));
560 color->jp2_cdef = (opj_jp2_cdef_t*)opj_malloc(sizeof(opj_jp2_cdef_t));
561 color->jp2_cdef->info = info;
562 color->jp2_cdef->n = n;
564 for(i = 0; i < n; ++i)
566 info[i].cn = (unsigned short)cio_read(cio, 2);
567 info[i].typ = (unsigned short)cio_read(cio, 2);
568 info[i].asoc = (unsigned short)cio_read(cio, 2);
572 }/* jp2_read_cdef() */
574 static opj_bool jp2_read_colr(opj_jp2_t *jp2, opj_cio_t *cio,
575 opj_jp2_box_t *box, opj_jp2_color_t *color)
578 opj_common_ptr cinfo;
580 /* Part 1, I.5.3.3 : 'A conforming JP2 reader shall ignore all Colour
581 * Specification boxes after the first.'
583 if(color->jp2_has_colr) return OPJ_FALSE;
587 jp2->meth = cio_read(cio, 1); /* METH */
588 jp2->precedence = cio_read(cio, 1); /* PRECEDENCE */
589 jp2->approx = cio_read(cio, 1); /* APPROX */
593 jp2->enumcs = cio_read(cio, 4); /* EnumCS */
599 skip_len = box->init_pos + box->length - cio_tell(cio);
602 opj_event_msg(cinfo, EVT_ERROR, "Error with COLR box size\n");
607 unsigned char *start;
609 start = cio_getbp(cio);
610 color->icc_profile_buf = (unsigned char*)opj_malloc(skip_len);
611 color->icc_profile_len = skip_len;
613 cio_skip(cio, box->init_pos + box->length - cio_tell(cio));
615 memcpy(color->icc_profile_buf, start, skip_len);
619 if (cio_tell(cio) - box->init_pos != box->length)
621 opj_event_msg(cinfo, EVT_ERROR, "Error with COLR Box\n");
624 color->jp2_has_colr = 1;
627 }/* jp2_read_colr() */
629 opj_bool jp2_read_jp2h(opj_jp2_t *jp2, opj_cio_t *cio, opj_jp2_color_t *color)
634 opj_common_ptr cinfo = jp2->cinfo;
636 jp2_read_boxhdr(cinfo, cio, &box);
639 if (JP2_JP2H != box.type)
641 if (box.type == JP2_JP2C)
643 opj_event_msg(cinfo, EVT_ERROR, "Expected JP2H Marker\n");
646 cio_skip(cio, box.length - 8);
648 if(cio->bp >= cio->end) return OPJ_FALSE;
650 jp2_read_boxhdr(cinfo, cio, &box);
652 } while(JP2_JP2H != box.type);
654 if (!jp2_read_ihdr(jp2, cio))
656 jp2h_end = box.init_pos + box.length;
660 if (!jp2_read_bpcc(jp2, cio))
663 jp2_read_boxhdr(cinfo, cio, &box);
665 while(cio_tell(cio) < jp2h_end)
667 if(box.type == JP2_COLR)
669 if( !jp2_read_colr(jp2, cio, &box, color))
671 cio_seek(cio, box.init_pos + 8);
672 cio_skip(cio, box.length - 8);
674 jp2_read_boxhdr(cinfo, cio, &box);
677 if(box.type == JP2_CDEF)
679 if( !jp2_read_cdef(jp2, cio, &box, color))
681 cio_seek(cio, box.init_pos + 8);
682 cio_skip(cio, box.length - 8);
684 jp2_read_boxhdr(cinfo, cio, &box);
687 if(box.type == JP2_PCLR)
689 if( !jp2_read_pclr(jp2, cio, &box, color))
691 cio_seek(cio, box.init_pos + 8);
692 cio_skip(cio, box.length - 8);
694 jp2_read_boxhdr(cinfo, cio, &box);
697 if(box.type == JP2_CMAP)
699 if( !jp2_read_cmap(jp2, cio, &box, color))
701 cio_seek(cio, box.init_pos + 8);
702 cio_skip(cio, box.length - 8);
704 jp2_read_boxhdr(cinfo, cio, &box);
707 cio_seek(cio, box.init_pos + 8);
708 cio_skip(cio, box.length - 8);
709 jp2_read_boxhdr(cinfo, cio, &box);
711 }/* while(cio_tell(cio) < box_end) */
713 cio_seek(cio, jp2h_end);
715 /* Part 1, I.5.3.3 : 'must contain at least one' */
716 return (color->jp2_has_colr == 1);
718 }/* jp2_read_jp2h() */
720 opj_image_t* opj_jp2_decode(opj_jp2_t *jp2, opj_cio_t *cio,
721 opj_codestream_info_t *cstr_info)
723 opj_common_ptr cinfo;
724 opj_image_t *image = NULL;
725 opj_jp2_color_t color;
731 memset(&color, 0, sizeof(opj_jp2_color_t));
735 if(!jp2_read_struct(jp2, cio, &color))
737 free_color_data(&color);
738 opj_event_msg(cinfo, EVT_ERROR, "Failed to decode jp2 structure\n");
743 image = j2k_decode(jp2->j2k, cio, cstr_info);
747 free_color_data(&color);
748 opj_event_msg(cinfo, EVT_ERROR, "Failed to decode J2K image\n");
752 /* Set Image Color Space */
753 if (jp2->enumcs == 16)
754 image->color_space = CLRSPC_SRGB;
755 else if (jp2->enumcs == 17)
756 image->color_space = CLRSPC_GRAY;
757 else if (jp2->enumcs == 18)
758 image->color_space = CLRSPC_SYCC;
760 image->color_space = CLRSPC_UNKNOWN;
764 jp2_apply_cdef(image, &color);
768 /* Part 1, I.5.3.4: Either both or none : */
769 if( !color.jp2_pclr->cmap)
770 jp2_free_pclr(&color);
772 jp2_apply_pclr(&color, image);
774 if(color.icc_profile_buf)
776 image->icc_profile_buf = color.icc_profile_buf;
777 color.icc_profile_buf = NULL;
778 image->icc_profile_len = color.icc_profile_len;
782 }/* opj_jp2_decode() */
785 void jp2_write_jp2h(opj_jp2_t *jp2, opj_cio_t *cio) {
788 box.init_pos = cio_tell(cio);
790 cio_write(cio, JP2_JP2H, 4); /* JP2H */
792 jp2_write_ihdr(jp2, cio);
794 if (jp2->bpc == 255) {
795 jp2_write_bpcc(jp2, cio);
797 jp2_write_colr(jp2, cio);
799 box.length = cio_tell(cio) - box.init_pos;
800 cio_seek(cio, box.init_pos);
801 cio_write(cio, box.length, 4); /* L */
802 cio_seek(cio, box.init_pos + box.length);
805 static void jp2_write_ftyp(opj_jp2_t *jp2, opj_cio_t *cio) {
809 box.init_pos = cio_tell(cio);
811 cio_write(cio, JP2_FTYP, 4); /* FTYP */
813 cio_write(cio, jp2->brand, 4); /* BR */
814 cio_write(cio, jp2->minversion, 4); /* MinV */
816 for (i = 0; i < jp2->numcl; i++) {
817 cio_write(cio, jp2->cl[i], 4); /* CL */
820 box.length = cio_tell(cio) - box.init_pos;
821 cio_seek(cio, box.init_pos);
822 cio_write(cio, box.length, 4); /* L */
823 cio_seek(cio, box.init_pos + box.length);
826 static opj_bool jp2_read_ftyp(opj_jp2_t *jp2, opj_cio_t *cio) {
830 opj_common_ptr cinfo = jp2->cinfo;
832 jp2_read_boxhdr(cinfo, cio, &box);
834 if (JP2_FTYP != box.type) {
835 opj_event_msg(cinfo, EVT_ERROR, "Expected FTYP Marker\n");
839 jp2->brand = cio_read(cio, 4); /* BR */
840 jp2->minversion = cio_read(cio, 4); /* MinV */
841 jp2->numcl = (box.length - 16) / 4;
842 jp2->cl = (unsigned int *) opj_malloc(jp2->numcl * sizeof(unsigned int));
844 for (i = 0; i < (int)jp2->numcl; i++) {
845 jp2->cl[i] = cio_read(cio, 4); /* CLi */
848 if (cio_tell(cio) - box.init_pos != box.length) {
849 opj_event_msg(cinfo, EVT_ERROR, "Error with FTYP Box\n");
856 static int jp2_write_jp2c(opj_jp2_t *jp2, opj_cio_t *cio, opj_image_t *image, opj_codestream_info_t *cstr_info) {
857 unsigned int j2k_codestream_offset, j2k_codestream_length;
860 opj_j2k_t *j2k = jp2->j2k;
862 box.init_pos = cio_tell(cio);
864 cio_write(cio, JP2_JP2C, 4); /* JP2C */
867 j2k_codestream_offset = cio_tell(cio);
868 if(!j2k_encode(j2k, cio, image, cstr_info)) {
869 opj_event_msg(j2k->cinfo, EVT_ERROR, "Failed to encode image\n");
872 j2k_codestream_length = cio_tell(cio) - j2k_codestream_offset;
874 jp2->j2k_codestream_offset = j2k_codestream_offset;
875 jp2->j2k_codestream_length = j2k_codestream_length;
877 box.length = 8 + jp2->j2k_codestream_length;
878 cio_seek(cio, box.init_pos);
879 cio_write(cio, box.length, 4); /* L */
880 cio_seek(cio, box.init_pos + box.length);
885 static opj_bool jp2_read_jp2c(opj_jp2_t *jp2, opj_cio_t *cio, unsigned int *j2k_codestream_length, unsigned int *j2k_codestream_offset) {
888 opj_common_ptr cinfo = jp2->cinfo;
890 jp2_read_boxhdr(cinfo, cio, &box);
892 if(JP2_JP2C != box.type) {
893 cio_skip(cio, box.length - 8);
894 jp2_read_boxhdr(cinfo, cio, &box);
896 } while(JP2_JP2C != box.type);
898 *j2k_codestream_offset = cio_tell(cio);
899 *j2k_codestream_length = box.length - 8;
904 static void jp2_write_jp(opj_cio_t *cio) {
907 box.init_pos = cio_tell(cio);
909 cio_write(cio, JP2_JP, 4); /* JP2 signature */
910 cio_write(cio, 0x0d0a870a, 4);
912 box.length = cio_tell(cio) - box.init_pos;
913 cio_seek(cio, box.init_pos);
914 cio_write(cio, box.length, 4); /* L */
915 cio_seek(cio, box.init_pos + box.length);
918 static opj_bool jp2_read_jp(opj_jp2_t *jp2, opj_cio_t *cio) {
921 opj_common_ptr cinfo = jp2->cinfo;
923 jp2_read_boxhdr(cinfo, cio, &box);
924 if (JP2_JP != box.type) {
925 opj_event_msg(cinfo, EVT_ERROR, "Expected JP Marker\n");
928 if (0x0d0a870a != cio_read(cio, 4)) {
929 opj_event_msg(cinfo, EVT_ERROR, "Error with JP Marker\n");
932 if (cio_tell(cio) - box.init_pos != box.length) {
933 opj_event_msg(cinfo, EVT_ERROR, "Error with JP Box size\n");
941 static opj_bool jp2_read_struct(opj_jp2_t *jp2, opj_cio_t *cio,
942 opj_jp2_color_t *color) {
943 if (!jp2_read_jp(jp2, cio))
945 if (!jp2_read_ftyp(jp2, cio))
947 if (!jp2_read_jp2h(jp2, cio, color))
949 if (!jp2_read_jp2c(jp2, cio, &jp2->j2k_codestream_length, &jp2->j2k_codestream_offset))
956 static int write_fidx( int offset_jp2c, int length_jp2c, int offset_idx, int length_idx, opj_cio_t *cio)
960 lenp = cio_tell( cio);
961 cio_skip( cio, 4); /* L [at the end] */
962 cio_write( cio, JPIP_FIDX, 4); /* IPTR */
964 write_prxy( offset_jp2c, length_jp2c, offset_idx, offset_jp2c, cio);
966 len = cio_tell( cio)-lenp;
967 cio_seek( cio, lenp);
968 cio_write( cio, len, 4); /* L */
969 cio_seek( cio, lenp+len);
974 static void write_prxy( int offset_jp2c, int length_jp2c, int offset_idx, int length_idx, opj_cio_t *cio)
978 lenp = cio_tell( cio);
979 cio_skip( cio, 4); /* L [at the end] */
980 cio_write( cio, JPIP_PRXY, 4); /* IPTR */
982 cio_write( cio, offset_jp2c, 8); /* OOFF */
983 cio_write( cio, length_jp2c, 4); /* OBH part 1 */
984 cio_write( cio, JP2_JP2C, 4); /* OBH part 2 */
986 cio_write( cio, 1,1); /* NI */
988 cio_write( cio, offset_idx, 8); /* IOFF */
989 cio_write( cio, length_idx, 4); /* IBH part 1 */
990 cio_write( cio, JPIP_CIDX, 4); /* IBH part 2 */
992 len = cio_tell( cio)-lenp;
993 cio_seek( cio, lenp);
994 cio_write( cio, len, 4); /* L */
995 cio_seek( cio, lenp+len);
998 static void write_iptr( int offset, int length, opj_cio_t *cio)
1002 lenp = cio_tell( cio);
1003 cio_skip( cio, 4); /* L [at the end] */
1004 cio_write( cio, JPIP_IPTR, 4); /* IPTR */
1006 cio_write( cio, offset, 8);
1007 cio_write( cio, length, 8);
1009 len = cio_tell( cio)-lenp;
1010 cio_seek( cio, lenp);
1011 cio_write( cio, len, 4); /* L */
1012 cio_seek( cio, lenp+len);
1016 /* ----------------------------------------------------------------------- */
1017 /* JP2 decoder interface */
1018 /* ----------------------------------------------------------------------- */
1020 opj_jp2_t* jp2_create_decompress(opj_common_ptr cinfo) {
1021 opj_jp2_t *jp2 = (opj_jp2_t*) opj_calloc(1, sizeof(opj_jp2_t));
1024 /* create the J2K codec */
1025 jp2->j2k = j2k_create_decompress(cinfo);
1026 if(jp2->j2k == NULL) {
1027 jp2_destroy_decompress(jp2);
1034 void jp2_destroy_decompress(opj_jp2_t *jp2) {
1036 /* destroy the J2K codec */
1037 j2k_destroy_decompress(jp2->j2k);
1040 opj_free(jp2->comps);
1049 void jp2_setup_decoder(opj_jp2_t *jp2, opj_dparameters_t *parameters) {
1050 /* setup the J2K codec */
1051 j2k_setup_decoder(jp2->j2k, parameters);
1052 /* further JP2 initializations go here */
1055 /* ----------------------------------------------------------------------- */
1056 /* JP2 encoder interface */
1057 /* ----------------------------------------------------------------------- */
1059 opj_jp2_t* jp2_create_compress(opj_common_ptr cinfo) {
1060 opj_jp2_t *jp2 = (opj_jp2_t*)opj_malloc(sizeof(opj_jp2_t));
1063 /* create the J2K codec */
1064 jp2->j2k = j2k_create_compress(cinfo);
1065 if(jp2->j2k == NULL) {
1066 jp2_destroy_compress(jp2);
1073 void jp2_destroy_compress(opj_jp2_t *jp2) {
1075 /* destroy the J2K codec */
1076 j2k_destroy_compress(jp2->j2k);
1079 opj_free(jp2->comps);
1088 void jp2_setup_encoder(opj_jp2_t *jp2, opj_cparameters_t *parameters, opj_image_t *image) {
1092 if(!jp2 || !parameters || !image)
1095 /* setup the J2K codec */
1096 /* ------------------- */
1098 /* Check if number of components respects standard */
1099 if (image->numcomps < 1 || image->numcomps > 16384) {
1100 opj_event_msg(jp2->cinfo, EVT_ERROR, "Invalid number of components specified while setting up JP2 encoder\n");
1104 j2k_setup_encoder(jp2->j2k, parameters, image);
1106 /* setup the JP2 codec */
1107 /* ------------------- */
1111 jp2->brand = JP2_JP2; /* BR */
1112 jp2->minversion = 0; /* MinV */
1114 jp2->cl = (unsigned int*) opj_malloc(jp2->numcl * sizeof(unsigned int));
1115 jp2->cl[0] = JP2_JP2; /* CL0 : JP2 */
1117 /* Image Header box */
1119 jp2->numcomps = image->numcomps; /* NC */
1120 jp2->comps = (opj_jp2_comps_t*) opj_malloc(jp2->numcomps * sizeof(opj_jp2_comps_t));
1121 jp2->h = image->y1 - image->y0; /* HEIGHT */
1122 jp2->w = image->x1 - image->x0; /* WIDTH */
1124 depth_0 = image->comps[0].prec - 1;
1125 sign = image->comps[0].sgnd;
1126 jp2->bpc = depth_0 + (sign << 7);
1127 for (i = 1; i < image->numcomps; i++) {
1128 int depth = image->comps[i].prec - 1;
1129 sign = image->comps[i].sgnd;
1130 if (depth_0 != depth)
1133 jp2->C = 7; /* C : Always 7 */
1134 jp2->UnkC = 0; /* UnkC, colorspace specified in colr box */
1135 jp2->IPR = 0; /* IPR, no intellectual property */
1137 /* BitsPerComponent box */
1139 for (i = 0; i < image->numcomps; i++) {
1140 jp2->comps[i].bpcc = image->comps[i].prec - 1 + (image->comps[i].sgnd << 7);
1143 if (image->color_space == 1)
1144 jp2->enumcs = 16; /* sRGB as defined by IEC 61966-2.1 */
1145 else if (image->color_space == 2)
1146 jp2->enumcs = 17; /* greyscale */
1147 else if (image->color_space == 3)
1148 jp2->enumcs = 18; /* YUV */
1149 jp2->precedence = 0; /* PRECEDENCE */
1150 jp2->approx = 0; /* APPROX */
1152 jp2->jpip_on = parameters->jpip_on;
1155 opj_bool opj_jp2_encode(opj_jp2_t *jp2, opj_cio_t *cio, opj_image_t *image, opj_codestream_info_t *cstr_info) {
1157 int pos_iptr, pos_cidx, pos_jp2c, len_jp2c, end_pos, pos_fidx, len_fidx;
1161 /* JPEG 2000 Signature box */
1164 jp2_write_ftyp(jp2, cio);
1165 /* JP2 Header box */
1166 jp2_write_jp2h(jp2, cio);
1169 pos_iptr = cio_tell( cio);
1170 cio_skip( cio, 24); /* IPTR further ! */
1172 pos_jp2c = cio_tell( cio);
1176 if(!(len_jp2c = jp2_write_jp2c( jp2, cio, image, cstr_info))){
1177 opj_event_msg(jp2->cinfo, EVT_ERROR, "Failed to encode image\n");
1182 pos_cidx = cio_tell( cio);
1184 write_cidx( pos_jp2c+8, cio, image, *cstr_info, len_jp2c-8);
1186 pos_fidx = cio_tell( cio);
1187 len_fidx = write_fidx( pos_jp2c, len_jp2c, pos_cidx, cio_tell(cio), cio);
1189 end_pos = cio_tell( cio);
1191 cio_seek( cio, pos_iptr);
1192 write_iptr( pos_fidx, len_fidx, cio);
1194 cio_seek( cio, end_pos);