1 /*=========================================================================
4 Module: $RCSfile: gdcmJpeg2000.cxx,v $
6 Date: $Date: 2005/06/02 12:37:21 $
7 Version: $Revision: 1.22 $
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);
64 char *fmtname = jas_image_fmttostr(fmtid);
65 printf("%s %d %d %d %d %ld\n", fmtname, numcmpts, width, height, prec, (long) jas_image_rawsize(jasImage));
67 // The following should serioulsy be rewritten I cannot belive we need to
68 // do a per pixel decompression, there should be a way to read a full
72 uint8_t *data8 = (uint8_t*)raw;
73 for ( i = 0; i < height; i++)
74 for ( j = 0; j < width; j++)
75 for ( k= 0; k < numcmpts; k++)
76 *data8++ = (uint8_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
80 uint16_t *data16 = (uint16_t*)raw;
81 for ( i = 0; i < height; i++)
82 for ( j = 0; j < width; j++)
83 for ( k= 0; k < numcmpts; k++)
84 *data16++ = (uint16_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
88 uint32_t *data32 = (uint32_t*)raw;
89 for ( i = 0; i < height; i++)
90 for ( j = 0; j < width; j++)
91 for ( k= 0; k < numcmpts; k++)
92 *data32++ = (uint32_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
95 jas_image_destroy(jasImage);
96 jas_image_clearfmts();
99 //delete the jpeg temp buffer
105 //-----------------------------------------------------------------------------
106 } // end namespace gdcm