1 /*=========================================================================
4 Module: $RCSfile: gdcmJpeg2000.cxx,v $
6 Date: $Date: 2005/05/22 18:38:52 $
7 Version: $Revision: 1.20 $
9 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
10 l'Image). All rights reserved. See Doc/License.txt or
11 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
13 This software is distributed WITHOUT ANY WARRANTY; without even
14 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 PURPOSE. See the above copyright notices for more information.
17 =========================================================================*/
18 #include "gdcmFileHelper.h"
19 #include "gdcmDebug.h"
23 #include <jasper/jasper.h>
27 //-----------------------------------------------------------------------------
29 * \brief routine for JPEG decompression
30 * @param fp pointer to an already open file descriptor
31 * JPEG2000 encoded image
32 * @param image_buffer to receive uncompressed pixels
33 * @return 1 on success, 0 on error
34 * @warning : not yet made
37 bool gdcm_read_JPEG2000_file (std::ifstream* fp, void* raw, size_t inputlength)
39 jas_init(); //important...
40 // FIXME this is really ugly but it seems I have to load the complete
41 // jpeg2000 stream to use jasper:
42 uint8_t *inputdata = new uint8_t[inputlength];
43 //fp is already 'seek' to proper pos
44 fp->read((char*)inputdata, inputlength);
45 jas_stream_t *jasStream = jas_stream_memopen((char *)inputdata, inputlength);
48 if ((fmtid = jas_image_getfmt(jasStream)) < 0)
50 gdcmErrorMacro("unknown image format");
55 jas_image_t *jasImage = NULL;
56 if (!(jasImage = jas_image_decode(jasStream, fmtid, 0)))
58 gdcmErrorMacro("cannot decode image");
63 jas_stream_close(jasStream);
64 int numcmpts = jas_image_numcmpts(jasImage);
65 int width = jas_image_cmptwidth(jasImage, 0);
66 int height = jas_image_cmptheight(jasImage, 0);
67 int prec = jas_image_cmptprec(jasImage, 0);
69 char *fmtname = jas_image_fmttostr(fmtid);
70 printf("%s %d %d %d %d %ld\n", fmtname, numcmpts, width, height, prec, (long) jas_image_rawsize(jasImage));
72 // The following should serioulsy be rewritten I cannot belive we need to
73 // do a per pixel decompression, there should be a way to read a full
77 uint8_t *data8 = (uint8_t*)raw;
78 for ( i = 0; i < height; i++)
79 for ( j = 0; j < width; j++)
80 for ( k= 0; k < numcmpts; k++)
82 (uint8_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
86 uint16_t *data16 = (uint16_t*)raw;
87 for ( i = 0; i < height; i++)
88 for ( j = 0; j < width; j++)
89 for ( k= 0; k < numcmpts; k++)
91 (uint16_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
95 uint32_t *data32 = (uint32_t*)raw;
96 for ( i = 0; i < height; i++)
97 for ( j = 0; j < width; j++)
98 for ( k= 0; k < numcmpts; k++)
100 (uint32_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
103 jas_image_destroy(jasImage);
104 jas_image_clearfmts();
107 //delete the jpeg temp buffer
113 //-----------------------------------------------------------------------------
114 } // end namespace gdcm