1 /*=========================================================================
4 Module: $RCSfile: gdcmJpeg2000.cxx,v $
6 Date: $Date: 2005/06/07 20:23:13 $
7 Version: $Revision: 1.24 $
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
31 * @param inputdata inputdata
32 * @param inputlength inputlength
33 * @return 1 on success, 0 on error
36 bool gdcm_read_JPEG2000_file (void* raw, char *inputdata, size_t inputlength)
38 jas_init(); //important...
39 jas_stream_t *jasStream =
40 jas_stream_memopen((char *)inputdata, inputlength);
43 if ((fmtid = jas_image_getfmt(jasStream)) < 0)
45 gdcmErrorMacro("unknown image format");
50 jas_image_t *jasImage = NULL;
51 if (!(jasImage = jas_image_decode(jasStream, fmtid, 0)))
53 gdcmErrorMacro("cannot decode image");
58 jas_stream_close(jasStream);
59 int numcmpts = jas_image_numcmpts(jasImage);
60 int width = jas_image_cmptwidth(jasImage, 0);
61 int height = jas_image_cmptheight(jasImage, 0);
62 int prec = jas_image_cmptprec(jasImage, 0);
65 // The following should serioulsy be rewritten I cannot belive we need to
66 // do a per pixel decompression, there should be a way to read a full
70 uint8_t *data8 = (uint8_t*)raw;
71 for ( i = 0; i < height; i++)
72 for ( j = 0; j < width; j++)
73 for ( k= 0; k < numcmpts; k++)
74 *data8++ = (uint8_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
78 uint16_t *data16 = (uint16_t*)raw;
79 for ( i = 0; i < height; i++)
80 for ( j = 0; j < width; j++)
81 for ( k= 0; k < numcmpts; k++)
82 *data16++ = (uint16_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
86 uint32_t *data32 = (uint32_t*)raw;
87 for ( i = 0; i < height; i++)
88 for ( j = 0; j < width; j++)
89 for ( k= 0; k < numcmpts; k++)
90 *data32++ = (uint32_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
93 jas_image_destroy(jasImage);
94 jas_image_clearfmts();
97 //delete the jpeg temp buffer
103 //-----------------------------------------------------------------------------
104 } // end namespace gdcm