src/lib/openjp2/image.c

   1 /*
   2  * Copyright (c) 2005, Herve Drolon, FreeImage Team
   3  * All rights reserved.
   4  *
   5  * Redistribution and use in source and binary forms, with or without
   6  * modification, are permitted provided that the following conditions
   7  * are met:
   8  * 1. Redistributions of source code must retain the above copyright
   9  *    notice, this list of conditions and the following disclaimer.
  10  * 2. Redistributions in binary form must reproduce the above copyright
  11  *    notice, this list of conditions and the following disclaimer in the
  12  *    documentation and/or other materials provided with the distribution.
  13  *
  14  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
  15  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  16  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  17  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  18  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  19  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  20  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  21  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  22  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  23  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  24  * POSSIBILITY OF SUCH DAMAGE.
  25  */
  26
  27 #include "opj_includes.h"
  28
  29 opj_image_t* opj_image_create0(void) {
  30         opj_image_t *image = (opj_image_t*)opj_calloc(1, sizeof(opj_image_t));
  31         return image;
  32 }
  33
  34 opj_image_t* OPJ_CALLCONV opj_image_create(OPJ_UINT32 numcmpts, opj_image_cmptparm_t *cmptparms, OPJ_COLOR_SPACE clrspc) {
  35         OPJ_UINT32 compno;
  36         opj_image_t *image = NULL;
  37
  38         image = (opj_image_t*) opj_calloc(1, sizeof(opj_image_t));
  39         if(image) {
  40                 image->color_space = clrspc;
  41                 image->numcomps = numcmpts;
  42                 /* allocate memory for the per-component information */
  43                 image->comps = (opj_image_comp_t*)opj_malloc(image->numcomps * sizeof(opj_image_comp_t));
  44                 if(!image->comps) {
  45                         fprintf(stderr,"Unable to allocate memory for image.\n");
  46                         opj_image_destroy(image);
  47                         return NULL;
  48                 }
  49                 /* create the individual image components */
  50                 for(compno = 0; compno < numcmpts; compno++) {
  51                         opj_image_comp_t *comp = &image->comps[compno];
  52                         comp->dx = cmptparms[compno].dx;
  53                         comp->dy = cmptparms[compno].dy;
  54                         comp->w = cmptparms[compno].w;
  55                         comp->h = cmptparms[compno].h;
  56                         comp->x0 = cmptparms[compno].x0;
  57                         comp->y0 = cmptparms[compno].y0;
  58                         comp->prec = cmptparms[compno].prec;
  59                         comp->bpp = cmptparms[compno].bpp;
  60                         comp->sgnd = cmptparms[compno].sgnd;
  61                         comp->data = (int*) opj_calloc(comp->w * comp->h, sizeof(int));
  62                         if(!comp->data) {
  63                                 fprintf(stderr,"Unable to allocate memory for image.\n");
  64                                 opj_image_destroy(image);
  65                                 return NULL;
  66                         }
  67                 }
  68         }
  69
  70         return image;
  71 }
  72
  73 void OPJ_CALLCONV opj_image_destroy(opj_image_t *image) {
  74         if(image) {
  75                 if(image->comps) {
  76                         OPJ_UINT32 compno;
  77
  78                         /* image components */
  79                         for(compno = 0; compno < image->numcomps; compno++) {
  80                                 opj_image_comp_t *image_comp = &(image->comps[compno]);
  81                                 if(image_comp->data) {
  82                                         opj_free(image_comp->data);
  83                                 }
  84                         }
  85                         opj_free(image->comps);
  86                 }
  87
  88                 if(image->icc_profile_buf) {
  89                         opj_free(image->icc_profile_buf);
  90                 }
  91
  92                 opj_free(image);
  93         }
  94 }
  95
  96 /**
  97  * Updates the components characteristics of the image from the coding parameters.
  98  *
  99  * @param p_image_header        the image header to update.
 100  * @param p_cp                          the coding parameters from which to update the image.
 101  */
 102 void opj_image_comp_header_update(opj_image_t * p_image_header, const struct opj_cp * p_cp)
 103 {
 104         OPJ_UINT32 i, l_width, l_height;
 105         OPJ_INT32 l_x0, l_y0, l_x1, l_y1;
 106         OPJ_INT32 l_comp_x0, l_comp_y0, l_comp_x1, l_comp_y1;
 107         opj_image_comp_t* l_img_comp = NULL;
 108
 109         l_x0 = opj_int_max(p_cp->tx0 , p_image_header->x0);
 110         l_y0 = opj_int_max(p_cp->ty0 , p_image_header->y0);
 111         l_x1 = opj_int_min(p_cp->tx0 + p_cp->tw * p_cp->tdx, p_image_header->x1);
 112         l_y1 = opj_int_min(p_cp->ty0 + p_cp->th * p_cp->tdy, p_image_header->y1);
 113
 114         l_img_comp = p_image_header->comps;
 115         for     (i = 0; i < p_image_header->numcomps; ++i) {
 116                 l_comp_x0 = opj_int_ceildiv(l_x0, l_img_comp->dx);
 117                 l_comp_y0 = opj_int_ceildiv(l_y0, l_img_comp->dy);
 118                 l_comp_x1 = opj_int_ceildiv(l_x1, l_img_comp->dx);
 119                 l_comp_y1 = opj_int_ceildiv(l_y1, l_img_comp->dy);
 120                 l_width = opj_int_ceildivpow2(l_comp_x1 - l_comp_x0, l_img_comp->factor);
 121                 l_height = opj_int_ceildivpow2(l_comp_y1 - l_comp_y0, l_img_comp->factor);
 122                 l_img_comp->w = l_width;
 123                 l_img_comp->h = l_height;
 124                 l_img_comp->x0 = l_comp_x0/*l_x0*/;
 125                 l_img_comp->y0 = l_comp_y0/*l_y0*/;
 126                 ++l_img_comp;
 127         }
 128 }
 129
 130
 131 /**
 132  * Copy only header of image and its component header (no data are copied)
 133  * if dest image have data, they will be freed
 134  *
 135  * @param       p_image_src             the src image
 136  * @param       p_image_dest    the dest image
 137  *
 138  */
 139 void opj_copy_image_header(const opj_image_t* p_image_src, opj_image_t* p_image_dest)
 140 {
 141         OPJ_UINT32 compno;
 142
 143         /* preconditions */
 144         assert(p_image_src != 00);
 145         assert(p_image_dest != 00);
 146
 147         p_image_dest->x0 = p_image_src->x0;
 148         p_image_dest->y0 = p_image_src->y0;
 149         p_image_dest->x1 = p_image_src->x1;
 150         p_image_dest->y1 = p_image_src->y1;
 151
 152         if (p_image_dest->comps){
 153                 for(compno = 0; compno < p_image_dest->numcomps; compno++) {
 154                         opj_image_comp_t *image_comp = &(p_image_dest->comps[compno]);
 155                         if(image_comp->data) {
 156                                 opj_free(image_comp->data);
 157                         }
 158                 }
 159                 opj_free(p_image_dest->comps);
 160                 p_image_dest->comps = NULL;
 161         }
 162
 163         p_image_dest->numcomps = p_image_src->numcomps;
 164
 165         p_image_dest->comps = (opj_image_comp_t*) opj_malloc(p_image_dest->numcomps * sizeof(opj_image_comp_t));
 166         if (!p_image_dest->comps){
 167                 p_image_dest->comps = NULL;
 168                 p_image_dest->numcomps = 0;
 169                 return;
 170         }
 171
 172         for (compno=0; compno < p_image_dest->numcomps; compno++){
 173                 memcpy( &(p_image_dest->comps[compno]),
 174                                 &(p_image_src->comps[compno]),
 175                                 sizeof(opj_image_comp_t));
 176                 p_image_dest->comps[compno].data = NULL;
 177         }
 178
 179         p_image_dest->color_space = p_image_src->color_space;
 180         p_image_dest->icc_profile_len = p_image_src->icc_profile_len;
 181
 182         if (p_image_dest->icc_profile_len) {
 183                 p_image_dest->icc_profile_buf = (OPJ_BYTE*)opj_malloc(p_image_dest->icc_profile_len);
 184                 if (!p_image_dest->icc_profile_buf){
 185                         p_image_dest->icc_profile_buf = NULL;
 186                         p_image_dest->icc_profile_len = 0;
 187                         return;
 188                 }
 189                 memcpy( p_image_dest->icc_profile_buf,
 190                                 p_image_src->icc_profile_buf,
 191                                 p_image_src->icc_profile_len);
 192                 }
 193                 else
 194                         p_image_dest->icc_profile_buf = NULL;
 195
 196         return;
 197 }
 198
 199 opj_image_t* OPJ_CALLCONV opj_image_tile_create(OPJ_UINT32 numcmpts, opj_image_cmptparm_t *cmptparms, OPJ_COLOR_SPACE clrspc) {
 200         OPJ_UINT32 compno;
 201         opj_image_t *image = 00;
 202
 203         image = (opj_image_t*) opj_malloc(sizeof(opj_image_t));
 204         if (image)
 205         {
 206                 memset(image,0,sizeof(opj_image_t));
 207
 208                 image->color_space = clrspc;
 209                 image->numcomps = numcmpts;
 210
 211                 /* allocate memory for the per-component information */
 212                 image->comps = (opj_image_comp_t*)opj_malloc(image->numcomps * sizeof(opj_image_comp_t));
 213                 if (!image->comps) {
 214                         opj_image_destroy(image);
 215                         return 00;
 216                 }
 217                 memset(image->comps,0,image->numcomps * sizeof(opj_image_comp_t));
 218
 219                 /* create the individual image components */
 220                 for(compno = 0; compno < numcmpts; compno++) {
 221                         opj_image_comp_t *comp = &image->comps[compno];
 222                         comp->dx = cmptparms[compno].dx;
 223                         comp->dy = cmptparms[compno].dy;
 224                         comp->w = cmptparms[compno].w;
 225                         comp->h = cmptparms[compno].h;
 226                         comp->x0 = cmptparms[compno].x0;
 227                         comp->y0 = cmptparms[compno].y0;
 228                         comp->prec = cmptparms[compno].prec;
 229                         comp->sgnd = cmptparms[compno].sgnd;
 230                         comp->data = 0;
 231                 }
 232         }
 233
 234         return image;
 235 }