1 /*=========================================================================
4 Module: $RCSfile: gdcmJpeg2000.cxx,v $
6 Date: $Date: 2005/07/05 20:58:27 $
7 Version: $Revision: 1.25 $
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)
39 std::ofstream out("/tmp/jpeg2000.jpc");
40 out.write((char*)inputdata,inputlength);
43 jas_init(); //important...
44 jas_stream_t *jasStream =
45 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);
70 // The following should serioulsy be rewritten I cannot belive we need to
71 // do a per pixel decompression, there should be a way to read a full
75 uint8_t *data8 = (uint8_t*)raw;
76 for ( i = 0; i < height; i++)
77 for ( j = 0; j < width; j++)
78 for ( k= 0; k < numcmpts; k++)
79 *data8++ = (uint8_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
83 uint16_t *data16 = (uint16_t*)raw;
84 for ( i = 0; i < height; i++)
85 for ( j = 0; j < width; j++)
86 for ( k= 0; k < numcmpts; k++)
87 *data16++ = (uint16_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
91 uint32_t *data32 = (uint32_t*)raw;
92 for ( i = 0; i < height; i++)
93 for ( j = 0; j < width; j++)
94 for ( k= 0; k < numcmpts; k++)
95 *data32++ = (uint32_t)(jas_image_readcmptsample(jasImage, k, j ,i ));
98 jas_image_destroy(jasImage);
99 jas_image_clearfmts();
102 //delete the jpeg temp buffer
108 //-----------------------------------------------------------------------------
109 } // end namespace gdcm