]> Creatis software - gdcm.git/blob - src/gdcmopenjpeg/libopenjpeg/jp2.c
ENH: A user does not know (and should not) what we internally use for representing...
[gdcm.git] / src / gdcmopenjpeg / libopenjpeg / jp2.c
1 /*
2 * Copyright (c) 2003-2004, Yannick Verschueren
3 * Copyright (c) 2003-2004,  Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
4 * All rights reserved.
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 * 1. Redistributions of source code must retain the above copyright
10 *    notice, this list of conditions and the following disclaimer.
11 * 2. Redistributions in binary form must reproduce the above copyright
12 *    notice, this list of conditions and the following disclaimer in the
13 *    documentation and/or other materials provided with the distribution.
14 *
15 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
16 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18 * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25 * POSSIBILITY OF SUCH DAMAGE.
26 */
27
28 #include <stdio.h>
29 #include <stdlib.h>
30 #include <string.h>
31
32 #include "j2k.h"
33 #include "jp2.h"
34 #include "cio.h"
35
36 #define JPIP_JPIP 0x6a706970
37
38 #define JP2_JP   0x6a502020
39 #define JP2_FTYP 0x66747970
40 #define JP2_JP2H 0x6a703268
41 #define JP2_IHDR 0x69686472
42 #define JP2_COLR 0x636f6c72
43 #define JP2_JP2C 0x6a703263
44 #define JP2_URL  0x75726c20
45 #define JP2_DBTL 0x6474626c
46 #define JP2_BPCC 0x62706363
47 #define JP2_JP2  0x6a703220
48
49 /*
50
51 * Read box headers
52 *
53 */
54
55 int jp2_read_boxhdr(jp2_box_t * box)
56 {
57   box->init_pos = cio_tell();
58   box->length = cio_read(4);
59   box->type = cio_read(4);
60   if (box->length == 1) {
61     if (cio_read(4) != 0) {
62       fprintf(stderr, "Error: Cannot handle box sizes higher than 2^32\n");
63       return 1;
64     }
65     box->length = cio_read(4);
66     if (box->length == 0)
67       box->length = cio_numbytesleft() + 12;
68   } else if (box->length == 0) {
69     box->length = cio_numbytesleft() + 8;
70   }
71   return 0;
72 }
73
74 /*
75
76 * Initialisation of a Standard JP2 structure
77 */
78
79 int jp2_init_stdjp2(jp2_struct_t * jp2_struct)
80 {
81
82   jp2_struct->comps =
83     (jp2_comps_t *) malloc(jp2_struct->numcomps * sizeof(jp2_comps_t));
84
85   jp2_struct->precedence = 0;   /* PRECEDENCE*/
86   jp2_struct->approx = 0;   /* APPROX*/
87
88   jp2_struct->brand = JP2_JP2;   /* BR         */
89   jp2_struct->minversion = 0;   /* MinV       */
90   jp2_struct->numcl = 1;
91   jp2_struct->cl = (unsigned int *) malloc(jp2_struct->numcl * sizeof(int));
92   jp2_struct->cl[0] = JP2_JP2;   /* CL0 : JP2  */
93
94   jp2_struct->C = 7;      /* C : Always 7*/
95   jp2_struct->UnkC = 0;      /* UnkC, colorspace specified in colr box*/
96   jp2_struct->IPR = 0;      /* IPR, no intellectual property*/
97
98   return 0;
99 }
100
101
102 void jp2_write_url(char *Idx_file)
103 {
104   unsigned int i;
105   char str[256];
106   jp2_box_t box;
107
108   sprintf(str, "%s", Idx_file);
109
110
111   box.init_pos = cio_tell();
112   cio_skip(4);
113   cio_write(JP2_URL, 4);   /* DBTL*/
114   cio_write(0, 1);      /* VERS*/
115   cio_write(0, 3);      /* FLAG*/
116
117   for (i = 0; i < strlen(str); i++) {
118     cio_write(str[i], 1);
119   }
120
121   box.length = cio_tell() - box.init_pos;
122   cio_seek(box.init_pos);
123   cio_write(box.length, 4);   /*    L       */
124   cio_seek(box.init_pos + box.length);
125 }
126
127 /*
128 * Read the IHDR box
129 *
130 * Image Header box
131 *
132 */
133 int jp2_read_ihdr(jp2_struct_t * jp2_struct)
134 {
135   jp2_box_t box;
136
137   jp2_read_boxhdr(&box);
138   if (JP2_IHDR != box.type) {
139     fprintf(stderr, "Error: Expected IHDR Marker\n");
140     return 1;
141   }
142
143   jp2_struct->h = cio_read(4);   /* HEIGHT*/
144   jp2_struct->w = cio_read(4);   /* WIDTH*/
145   jp2_struct->numcomps = cio_read(2);   /* NC*/
146
147   jp2_struct->bpc = cio_read(1);   /* BPC*/
148
149   jp2_struct->C = cio_read(1);   /* C */
150   jp2_struct->UnkC = cio_read(1);   /* UnkC*/
151   jp2_struct->IPR = cio_read(1);   /* IPR*/
152
153   if (cio_tell() - box.init_pos != box.length) {
154     fprintf(stderr, "Error with IHDR Box\n");
155     return 1;
156   }
157   return 0;
158 }
159
160 void jp2_write_ihdr(jp2_struct_t * jp2_struct)
161 {
162   jp2_box_t box;
163
164   box.init_pos = cio_tell();
165   cio_skip(4);
166   cio_write(JP2_IHDR, 4);   /* IHDR*/
167
168   cio_write(jp2_struct->h, 4);   /* HEIGHT*/
169   cio_write(jp2_struct->w, 4);   /* WIDTH*/
170   cio_write(jp2_struct->numcomps, 2);   /* NC*/
171
172   cio_write(jp2_struct->bpc, 1);   /* BPC  */
173
174   cio_write(jp2_struct->C, 1);   /* C : Always 7*/
175   cio_write(jp2_struct->UnkC, 1);   /* UnkC, colorspace unknow*/
176   cio_write(jp2_struct->IPR, 1);   /* IPR, no intellectual property*/
177
178   box.length = cio_tell() - box.init_pos;
179   cio_seek(box.init_pos);
180   cio_write(box.length, 4);   /*    L       */
181   cio_seek(box.init_pos + box.length);
182 }
183
184
185 void jp2_write_bpcc(jp2_struct_t * jp2_struct)
186 {
187   unsigned int i;
188   jp2_box_t box;
189
190   box.init_pos = cio_tell();
191   cio_skip(4);
192   cio_write(JP2_BPCC, 4);   /* BPCC*/
193
194   for (i = 0; i < jp2_struct->numcomps; i++)
195     cio_write(jp2_struct->comps[i].bpcc, 1);
196
197   box.length = cio_tell() - box.init_pos;
198   cio_seek(box.init_pos);
199   cio_write(box.length, 4);   /*    L       */
200   cio_seek(box.init_pos + box.length);
201 }
202
203
204 int jp2_read_bpcc(jp2_struct_t * jp2_struct)
205 {
206   unsigned int i;
207   jp2_box_t box;
208
209   jp2_read_boxhdr(&box);
210   if (JP2_BPCC != box.type) {
211     fprintf(stderr, "Error: Expected BPCC Marker\n");
212     return 1;
213   }
214
215   for (i = 0; i < jp2_struct->numcomps; i++)
216     jp2_struct->comps[i].bpcc = cio_read(1);
217
218   if (cio_tell() - box.init_pos != box.length) {
219     fprintf(stderr, "Error with BPCC Box\n");
220     return 1;
221   }
222   return 0;
223 }
224
225 void jp2_write_colr(jp2_struct_t * jp2_struct)
226 {
227   jp2_box_t box;
228
229   box.init_pos = cio_tell();
230   cio_skip(4);
231   cio_write(JP2_COLR, 4);   /* COLR*/
232
233   cio_write(jp2_struct->meth, 1);   /* METH*/
234   cio_write(jp2_struct->precedence, 1);   /* PRECEDENCE*/
235   cio_write(jp2_struct->approx, 1);   /* APPROX*/
236
237   if (jp2_struct->meth == 1)
238     cio_write(jp2_struct->enumcs, 4);   /* EnumCS*/
239   else
240     cio_write(0, 1);      /* PROFILE (??)*/
241
242   box.length = cio_tell() - box.init_pos;
243   cio_seek(box.init_pos);
244   cio_write(box.length, 4);   /*    L       */
245   cio_seek(box.init_pos + box.length);
246 }
247
248 int jp2_read_colr(jp2_struct_t * jp2_struct)
249 {
250   jp2_box_t box;
251   int skip_len;
252
253   jp2_read_boxhdr(&box);
254   do {
255     if (JP2_COLR != box.type) {
256       cio_skip(box.length - 8);
257       jp2_read_boxhdr(&box);
258     }
259   } while (JP2_COLR != box.type);
260
261   jp2_struct->meth = cio_read(1);   /* METH*/
262   jp2_struct->precedence = cio_read(1);   /* PRECEDENCE*/
263   jp2_struct->approx = cio_read(1);   /* APPROX*/
264
265   if (jp2_struct->meth == 1)
266     jp2_struct->enumcs = cio_read(4);   /* EnumCS*/
267   else {
268     /* SKIP PROFILE     */
269     skip_len = box.init_pos + box.length - cio_tell();
270     if (skip_len < 0) {
271       fprintf(stderr, "Error with JP2H box size\n");
272       return 1;
273     }
274     cio_skip(box.init_pos + box.length - cio_tell());
275   }
276
277   if (cio_tell() - box.init_pos != box.length) {
278     fprintf(stderr, "Error with BPCC Box\n");
279     return 1;
280   }
281   return 0;
282 }
283
284 /*
285 * Write the JP2H box
286 *
287 * JP2 Header box
288 *
289 */
290 void jp2_write_jp2h(jp2_struct_t * jp2_struct)
291 {
292   jp2_box_t box;
293
294   box.init_pos = cio_tell();
295   cio_skip(4);;
296   cio_write(JP2_JP2H, 4);   /* JP2H */
297
298   jp2_write_ihdr(jp2_struct);
299
300   if (jp2_struct->bpc == 255)
301     jp2_write_bpcc(jp2_struct);
302   jp2_write_colr(jp2_struct);
303
304   box.length = cio_tell() - box.init_pos;
305   cio_seek(box.init_pos);
306   cio_write(box.length, 4);   /*    L       */
307   cio_seek(box.init_pos + box.length);
308 }
309
310
311 /*
312 * Read the JP2H box
313 *
314 * JP2 Header box
315 *
316 */
317 int jp2_read_jp2h(jp2_struct_t * jp2_struct)
318 {
319   jp2_box_t box;
320   int skip_len;
321
322   jp2_read_boxhdr(&box);
323   do {
324     if (JP2_JP2H != box.type) {
325       if (box.type == JP2_JP2C) {
326    fprintf(stderr, "Error: Expected JP2H Marker\n");
327    return 1;
328       }
329       cio_skip(box.length - 8);
330       jp2_read_boxhdr(&box);
331     }
332   } while (JP2_JP2H != box.type);
333
334   if (jp2_read_ihdr(jp2_struct))
335     return 1;
336
337   if (jp2_struct->bpc == 255) {
338     if (jp2_read_bpcc(jp2_struct))
339       return 1;
340   }
341
342   if (jp2_read_colr(jp2_struct))
343     return 1;
344
345   skip_len = box.init_pos + box.length - cio_tell();
346   if (skip_len < 0) {
347     fprintf(stderr, "Error with JP2H box size\n");
348     return 1;
349   }
350   cio_skip(box.init_pos + box.length - cio_tell());
351
352   return 0;
353 }
354
355 /*
356 * Write the FTYP box
357 *
358 * File type box
359 *
360 */
361 void jp2_write_ftyp(jp2_struct_t * jp2_struct)
362 {
363   unsigned int i;
364   jp2_box_t box;
365
366   box.init_pos = cio_tell();
367   cio_skip(4);
368   cio_write(JP2_FTYP, 4);   /* FTYP       */
369
370   cio_write(jp2_struct->brand, 4);   /* BR         */
371   cio_write(jp2_struct->minversion, 4);   /* MinV       */
372
373   for (i = 0; i < jp2_struct->numcl; i++)
374     cio_write(jp2_struct->cl[i], 4);   /* CL           */
375
376   box.length = cio_tell() - box.init_pos;
377   cio_seek(box.init_pos);
378   cio_write(box.length, 4);   /*    L       */
379   cio_seek(box.init_pos + box.length);
380 }
381
382 /*
383 * Read the FTYP box
384 *
385 * File type box
386 *
387 */
388 int jp2_read_ftyp(jp2_struct_t * jp2_struct)
389 {
390   int i;
391   jp2_box_t box;
392
393   jp2_read_boxhdr(&box);
394
395   if (JP2_FTYP != box.type) {
396     fprintf(stderr, "Error: Excpected FTYP Marker\n");
397     return 1;
398   }
399
400   jp2_struct->brand = cio_read(4);   /* BR              */
401   jp2_struct->minversion = cio_read(4);   /* MinV            */
402   jp2_struct->numcl = (box.length - 16) / 4;
403   jp2_struct->cl =
404     (unsigned int *) malloc(jp2_struct->numcl * sizeof(unsigned int));
405
406   for (i = 0; i < (int) jp2_struct->numcl; i++)
407     jp2_struct->cl[i] = cio_read(4);   /* CLi */
408
409   if (cio_tell() - box.init_pos != box.length) {
410     fprintf(stderr, "Error with FTYP Box\n");
411     return 1;
412   }
413   return 0;
414 }
415
416 int jp2_write_jp2c(int j2k_codestream_len, int *j2k_codestream_offset,
417          char *j2k_codestream)
418 {
419   jp2_box_t box;
420
421   box.init_pos = cio_tell();
422   cio_skip(4);
423   cio_write(JP2_JP2C, 4);   /* JP2C*/
424
425   *j2k_codestream_offset = cio_tell();
426   memcpy(cio_getbp(), j2k_codestream, j2k_codestream_len);
427
428   box.length = 8 + j2k_codestream_len;
429   cio_seek(box.init_pos);
430   cio_write(box.length, 4);   /*    L       */
431   cio_seek(box.init_pos + box.length);
432
433   return box.length;
434 }
435
436
437 int jp2_read_jp2c(unsigned int *j2k_codestream_len,
438         unsigned int *j2k_codestream_offset)
439 {
440   jp2_box_t box;
441
442   jp2_read_boxhdr(&box);
443   do {
444     if (JP2_JP2C != box.type) {
445       cio_skip(box.length - 8);
446       jp2_read_boxhdr(&box);
447     }
448   } while (JP2_JP2C != box.type);
449
450   *j2k_codestream_offset = cio_tell();
451   *j2k_codestream_len = box.length - 8;
452
453   return 0;
454 }
455
456 void jp2_write_jp()
457 {
458   jp2_box_t box;
459
460   box.init_pos = cio_tell();
461   cio_skip(4);
462   cio_write(JP2_JP, 4);      /* JP*/
463   cio_write(0x0d0a870a, 4);
464
465   box.length = cio_tell() - box.init_pos;
466   cio_seek(box.init_pos);
467   cio_write(box.length, 4);   /*    L       */
468   cio_seek(box.init_pos + box.length);
469 }
470
471 /*
472 * Read the JP box
473 *
474 * JPEG 2000 signature
475 *
476 * return 1 if error else 0
477 */
478 int jp2_read_jp()
479 {
480   jp2_box_t box;
481
482   jp2_read_boxhdr(&box);
483   if (JP2_JP != box.type) {
484     fprintf(stderr, "Error: Expected JP Marker\n");
485     return 1;
486   }
487   if (0x0d0a870a != cio_read(4)) {
488     fprintf(stderr, "Error with JP Marker\n");
489     return 1;
490   }
491   if (cio_tell() - box.init_pos != box.length) {
492     fprintf(stderr, "Error with JP Box size\n");
493     return 1;
494   }
495   return 0;
496
497 }
498
499
500 int jp2_read_struct(unsigned char *src, jp2_struct_t * jp2_struct, int len)
501 {
502   cio_init(src, len);
503
504   if (jp2_read_jp())
505     return 1;
506   if (jp2_read_ftyp(jp2_struct))
507     return 1;
508   if (jp2_read_jp2h(jp2_struct))
509     return 1;
510   if (jp2_read_jp2c
511       (&jp2_struct->j2k_codestream_len,
512        &jp2_struct->j2k_codestream_offset))
513     return 1;
514   return 0;
515 }
516
517 int jp2_wrap_j2k(jp2_struct_t * jp2_struct, char *j2k_codestream,
518        char *output)
519 {
520   (void)output;
521   jp2_write_jp();
522   jp2_write_ftyp(jp2_struct);
523   jp2_write_jp2h(jp2_struct);
524
525   jp2_write_jp2c(jp2_struct->j2k_codestream_len,
526        &jp2_struct->j2k_codestream_offset, j2k_codestream);
527
528   return cio_tell();
529 }