OPJViewer/source/imagj2k.cpp

   1 /*\r
   2  * Copyright (c) 2007, Digital Signal Processing Laboratory, Universit� degli studi di Perugia (UPG), Italy\r
   3  * All rights reserved.\r
   4  *\r
   5  * Redistribution and use in source and binary forms, with or without\r
   6  * modification, are permitted provided that the following conditions\r
   7  * are met:\r
   8  * 1. Redistributions of source code must retain the above copyright\r
   9  *    notice, this list of conditions and the following disclaimer.\r
  10  * 2. Redistributions in binary form must reproduce the above copyright\r
  11  *    notice, this list of conditions and the following disclaimer in the\r
  12  *    documentation and/or other materials provided with the distribution.\r
  13  *\r
  14  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'\r
  15  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE\r
  16  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE\r
  17  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE\r
  18  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR\r
  19  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF\r
  20  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS\r
  21  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN\r
  22  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)\r
  23  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE\r
  24  * POSSIBILITY OF SUCH DAMAGE.\r
  25  */\r
  26 /////////////////////////////////////////////////////////////////////////////\r
  27 // Name:        imagj2k.cpp\r
  28 // Purpose:     wxImage JPEG 2000 codestream handler\r
  29 // Author:      Giuseppe Baruffa - based on imagjpeg.cpp, Vaclav Slavik\r
  30 // RCS-ID:      $Id: imagj2k.cpp,v 0.00 2007/02/08 23:59:00 MW Exp $\r
  31 // Copyright:   (c) Giuseppe Baruffa\r
  32 // Licence:     wxWindows licence\r
  33 /////////////////////////////////////////////////////////////////////////////\r
  34 \r
  35 // For compilers that support precompilation, includes "wx.h".\r
  36 #include "wx/wxprec.h"\r
  37 \r
  38 #ifdef __BORLANDC__\r
  39     #pragma hdrstop\r
  40 #endif\r
  41 \r
  42 #if wxUSE_IMAGE && wxUSE_LIBOPENJPEG\r
  43 \r
  44 #include "imagj2k.h"\r
  45 \r
  46 #ifndef WX_PRECOMP\r
  47     #include "wx/log.h"\r
  48     #include "wx/app.h"\r
  49     #include "wx/intl.h"\r
  50     #include "wx/bitmap.h"\r
  51     #include "wx/module.h"\r
  52 #endif\r
  53 \r
  54 #include "wx/filefn.h"\r
  55 #include "wx/wfstream.h"\r
  56 \r
  57 // ----------------------------------------------------------------------------\r
  58 // types\r
  59 // ----------------------------------------------------------------------------\r
  60 \r
  61 \r
  62 //-----------------------------------------------------------------------------\r
  63 // wxJ2KHandler\r
  64 //-----------------------------------------------------------------------------\r
  65 \r
  66 IMPLEMENT_DYNAMIC_CLASS(wxJ2KHandler,wxImageHandler)\r
  67 \r
  68 #if wxUSE_STREAMS\r
  69 \r
  70 //------------- JPEG 2000 Data Source Manager\r
  71 \r
  72 #define J2K_CFMT 0\r
  73 #define JP2_CFMT 1\r
  74 #define JPT_CFMT 2\r
  75 #define MJ2_CFMT 3\r
  76 #define PXM_DFMT 0\r
  77 #define PGX_DFMT 1\r
  78 #define BMP_DFMT 2\r
  79 #define YUV_DFMT 3\r
  80 \r
  81 #define MAX_MESSAGE_LEN 200\r
  82 \r
  83 /* sample error callback expecting a FILE* client object */\r
  84 void j2k_error_callback(const char *msg, void *client_data) {\r
  85         int message_len = strlen(msg) - 1;\r
  86         if (msg[message_len] != '\n')\r
  87                 message_len = MAX_MESSAGE_LEN;\r
  88     /*wxMutexGuiEnter();\r
  89         wxLogMessage(wxT("[ERROR] %.*s"), message_len, msg);\r
  90     wxMutexGuiLeave();*/\r
  91 }\r
  92 /* sample warning callback expecting a FILE* client object */\r
  93 void j2k_warning_callback(const char *msg, void *client_data) {\r
  94         int message_len = strlen(msg) - 1;\r
  95         if (msg[message_len] != '\n')\r
  96                 message_len = MAX_MESSAGE_LEN;\r
  97     /*wxMutexGuiEnter();\r
  98         wxLogMessage(wxT("[WARNING] %.*s"), message_len, msg);\r
  99     wxMutexGuiLeave();*/\r
 100 }\r
 101 /* sample debug callback expecting no client object */\r
 102 void j2k_info_callback(const char *msg, void *client_data) {\r
 103         int message_len = strlen(msg) - 1;\r
 104         if (msg[message_len] != '\n')\r
 105                 message_len = MAX_MESSAGE_LEN;\r
 106     /*wxMutexGuiEnter();\r
 107         wxLogMessage(wxT("[INFO] %.*s"), message_len, msg);\r
 108     wxMutexGuiLeave();*/\r
 109 }\r
 110 \r
 111 // load the j2k codestream\r
 112 bool wxJ2KHandler::LoadFile(wxImage *image, wxInputStream& stream, bool verbose, int index)\r
 113 {\r
 114         opj_dparameters_t parameters;   /* decompression parameters */\r
 115         opj_event_mgr_t event_mgr;              /* event manager */\r
 116         opj_image_t *opjimage = NULL;\r
 117         unsigned char *src = NULL;\r
 118     unsigned char *ptr;\r
 119         int file_length;\r
 120         int shiftbpp;\r
 121 \r
 122         // destroy the image\r
 123     image->Destroy();\r
 124 \r
 125         /* handle to a decompressor */\r
 126         opj_dinfo_t* dinfo = NULL;      \r
 127         opj_cio_t *cio = NULL;\r
 128 \r
 129 \r
 130         /* configure the event callbacks (not required) */\r
 131         memset(&event_mgr, 0, sizeof(opj_event_mgr_t));\r
 132         event_mgr.error_handler = j2k_error_callback;\r
 133         event_mgr.warning_handler = j2k_warning_callback;\r
 134         event_mgr.info_handler = j2k_info_callback;\r
 135 \r
 136         /* set decoding parameters to default values */\r
 137         opj_set_default_decoder_parameters(&parameters);\r
 138 \r
 139         /* prepare parameters */\r
 140         strncpy(parameters.infile, "", sizeof(parameters.infile)-1);\r
 141         strncpy(parameters.outfile, "", sizeof(parameters.outfile)-1);\r
 142         parameters.decod_format = J2K_CFMT;\r
 143         parameters.cod_format = BMP_DFMT;\r
 144         if (m_reducefactor)\r
 145                 parameters.cp_reduce = m_reducefactor;\r
 146         if (m_qualitylayers)\r
 147                 parameters.cp_layer = m_qualitylayers;\r
 148         /*if (n_components)\r
 149                 parameters. = n_components;*/\r
 150 \r
 151         /* JPWL only */\r
 152 #ifdef USE_JPWL\r
 153         parameters.jpwl_exp_comps = m_expcomps;\r
 154         parameters.jpwl_max_tiles = m_maxtiles;\r
 155         parameters.jpwl_correct = m_enablejpwl;\r
 156 #endif /* USE_JPWL */\r
 157 \r
 158         /* get a decoder handle */\r
 159         dinfo = opj_create_decompress(CODEC_J2K);\r
 160 \r
 161         /* find length of the stream */\r
 162         stream.SeekI(0, wxFromEnd);\r
 163         file_length = (int) stream.TellI();\r
 164 \r
 165         /* get data */\r
 166         stream.SeekI(0, wxFromStart);\r
 167     src = (unsigned char *) malloc(file_length);\r
 168         stream.Read(src, file_length);\r
 169 \r
 170         /* catch events using our callbacks and give a local context */\r
 171         opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);\r
 172 \r
 173         /* setup the decoder decoding parameters using user parameters */\r
 174         opj_setup_decoder(dinfo, &parameters);\r
 175 \r
 176         /* open a byte stream */\r
 177         cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);\r
 178 \r
 179         /* decode the stream and fill the image structure */\r
 180         opjimage = opj_decode(dinfo, cio);\r
 181         if (!opjimage) {\r
 182                 wxMutexGuiEnter();\r
 183                 wxLogError(wxT("J2K: failed to decode image!"));\r
 184                 wxMutexGuiLeave();\r
 185                 opj_destroy_decompress(dinfo);\r
 186                 opj_cio_close(cio);\r
 187                 opj_image_destroy(opjimage);\r
 188                 free(src);\r
 189                 return false;\r
 190         }\r
 191 \r
 192         /* close the byte stream */\r
 193         opj_cio_close(cio);\r
 194 \r
 195         // check image components\r
 196         if ((opjimage->numcomps != 1) && (opjimage->numcomps != 3)) {\r
 197                 wxMutexGuiEnter();\r
 198                 wxLogError(wxT("J2K: weird number of components"));\r
 199                 wxMutexGuiLeave();\r
 200                 opj_destroy_decompress(dinfo);\r
 201                 free(src);\r
 202                 return false;\r
 203         }\r
 204 \r
 205         // check image depth (only on the first one, for now)\r
 206         shiftbpp = opjimage->comps[0].prec - 8;\r
 207 \r
 208         // prepare image size\r
 209     image->Create(opjimage->comps[0].w, opjimage->comps[0].h, true );\r
 210 \r
 211         // access image raw data\r
 212     image->SetMask( false );\r
 213     ptr = image->GetData();\r
 214 \r
 215         // RGB color picture\r
 216         if (opjimage->numcomps == 3) {\r
 217                 int row, col;\r
 218                 int *r = opjimage->comps[0].data;\r
 219                 int *g = opjimage->comps[1].data;\r
 220                 int *b = opjimage->comps[2].data;\r
 221                 if (shiftbpp > 0) {\r
 222                         for (row = 0; row < opjimage->comps[0].h; row++) {\r
 223                                 for (col = 0; col < opjimage->comps[0].w; col++) {\r
 224                                         \r
 225                                         *(ptr++) = (*(r++)) >> shiftbpp;\r
 226                                         *(ptr++) = (*(g++)) >> shiftbpp;\r
 227                                         *(ptr++) = (*(b++)) >> shiftbpp;\r
 228 \r
 229                                 }\r
 230                         }\r
 231 \r
 232                 } else if (shiftbpp < 0) {\r
 233                         for (row = 0; row < opjimage->comps[0].h; row++) {\r
 234                                 for (col = 0; col < opjimage->comps[0].w; col++) {\r
 235                                         \r
 236                                         *(ptr++) = (*(r++)) << -shiftbpp;\r
 237                                         *(ptr++) = (*(g++)) << -shiftbpp;\r
 238                                         *(ptr++) = (*(b++)) << -shiftbpp;\r
 239 \r
 240                                 }\r
 241                         }\r
 242                         \r
 243                 } else {\r
 244                         for (row = 0; row < opjimage->comps[0].h; row++) {\r
 245                                 for (col = 0; col < opjimage->comps[0].w; col++) {\r
 246 \r
 247                                         *(ptr++) = *(r++);\r
 248                                         *(ptr++) = *(g++);\r
 249                                         *(ptr++) = *(b++);\r
 250                                 \r
 251                                 }\r
 252                         }\r
 253                 }\r
 254         }\r
 255 \r
 256         // B/W picture\r
 257         if (opjimage->numcomps == 1) {\r
 258                 int row, col;\r
 259                 int *y = opjimage->comps[0].data;\r
 260                 if (shiftbpp > 0) {\r
 261                         for (row = 0; row < opjimage->comps[0].h; row++) {\r
 262                                 for (col = 0; col < opjimage->comps[0].w; col++) {\r
 263                                         \r
 264                                         *(ptr++) = (*(y)) >> shiftbpp;\r
 265                                         *(ptr++) = (*(y)) >> shiftbpp;\r
 266                                         *(ptr++) = (*(y++)) >> shiftbpp;\r
 267 \r
 268                                 }\r
 269                         }\r
 270                 } else if (shiftbpp < 0) {\r
 271                         for (row = 0; row < opjimage->comps[0].h; row++) {\r
 272                                 for (col = 0; col < opjimage->comps[0].w; col++) {\r
 273                                         \r
 274                                         *(ptr++) = (*(y)) << -shiftbpp;\r
 275                                         *(ptr++) = (*(y)) << -shiftbpp;\r
 276                                         *(ptr++) = (*(y++)) << -shiftbpp;\r
 277 \r
 278                                 }\r
 279                         }\r
 280                 } else {\r
 281                         for (row = 0; row < opjimage->comps[0].h; row++) {\r
 282                                 for (col = 0; col < opjimage->comps[0].w; col++) {\r
 283                                         \r
 284                                         *(ptr++) = *(y);\r
 285                                         *(ptr++) = *(y);\r
 286                                         *(ptr++) = *(y++);\r
 287 \r
 288                                 }\r
 289                         }\r
 290                 }\r
 291         }\r
 292 \r
 293     wxMutexGuiEnter();\r
 294     wxLogMessage(wxT("J2K: image loaded."));\r
 295     wxMutexGuiLeave();\r
 296 \r
 297         /* close openjpeg structs */\r
 298         opj_destroy_decompress(dinfo);\r
 299         opj_image_destroy(opjimage);\r
 300         free(src);\r
 301 \r
 302         if (!image->Ok())\r
 303                 return false;\r
 304         else\r
 305                 return true;\r
 306 \r
 307 }\r
 308 \r
 309 // save the j2k codestream\r
 310 bool wxJ2KHandler::SaveFile( wxImage *image, wxOutputStream& stream, bool verbose )\r
 311 {\r
 312     wxLogError(wxT("J2K: Couldn't save image -> not implemented."));\r
 313     return false;\r
 314 }\r
 315 \r
 316 #ifdef __VISUALC__\r
 317     #pragma warning(default:4611)\r
 318 #endif /* VC++ */\r
 319 \r
 320 // recognize the 0xFF4F JPEG 2000 SOC marker\r
 321 bool wxJ2KHandler::DoCanRead( wxInputStream& stream )\r
 322 {\r
 323     unsigned char hdr[2];\r
 324 \r
 325     if ( !stream.Read(hdr, WXSIZEOF(hdr)) )\r
 326         return false;\r
 327 \r
 328     return hdr[0] == 0xFF && hdr[1] == 0x4F;\r
 329 }\r
 330 \r
 331 #endif   // wxUSE_STREAMS\r
 332 \r
 333 #endif   // wxUSE_LIBOPENJPEG\r