libopenjpeg/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_header_t* opj_image_header_create0(void) {
  35         opj_image_header_t *image_header = (opj_image_header_t*)opj_calloc(1, sizeof(opj_image_header_t));
  36         return image_header;
  37 }
  38
  39 opj_image_t* OPJ_CALLCONV opj_image_create(int numcmpts, opj_image_cmptparm_t *cmptparms, OPJ_COLOR_SPACE clrspc) {
  40         int compno;
  41         opj_image_t *image = NULL;
  42
  43         image = (opj_image_t*) opj_calloc(1, sizeof(opj_image_t));
  44         if(image) {
  45                 image->color_space = clrspc;
  46                 image->numcomps = numcmpts;
  47                 /* allocate memory for the per-component information */
  48                 image->comps = (opj_image_comp_t*)opj_malloc(image->numcomps * sizeof(opj_image_comp_t));
  49                 if(!image->comps) {
  50                         fprintf(stderr,"Unable to allocate memory for image.\n");
  51                         opj_image_destroy(image);
  52                         return NULL;
  53                 }
  54                 /* create the individual image components */
  55                 for(compno = 0; compno < numcmpts; compno++) {
  56                         opj_image_comp_t *comp = &image->comps[compno];
  57                         comp->dx = cmptparms[compno].dx;
  58                         comp->dy = cmptparms[compno].dy;
  59                         comp->w = cmptparms[compno].w;
  60                         comp->h = cmptparms[compno].h;
  61                         comp->x0 = cmptparms[compno].x0;
  62                         comp->y0 = cmptparms[compno].y0;
  63                         comp->prec = cmptparms[compno].prec;
  64                         comp->bpp = cmptparms[compno].bpp;
  65                         comp->sgnd = cmptparms[compno].sgnd;
  66                         comp->data = (int*) opj_calloc(comp->w * comp->h, sizeof(int));
  67                         if(!comp->data) {
  68                                 fprintf(stderr,"Unable to allocate memory for image.\n");
  69                                 opj_image_destroy(image);
  70                                 return NULL;
  71                         }
  72                 }
  73         }
  74
  75         return image;
  76 }
  77
  78 // TODO remove this one
  79 void OPJ_CALLCONV opj_image_destroy(opj_image_t *image) {
  80         int i;
  81         if(image) {
  82                 if(image->comps) {
  83                         /* image components */
  84                         for(i = 0; i < image->numcomps; i++) {
  85                                 opj_image_comp_t *image_comp = &image->comps[i];
  86                                 if(image_comp->data) {
  87                                         opj_free(image_comp->data);
  88                                 }
  89                         }
  90                         opj_free(image->comps);
  91                 }
  92                 opj_free(image);
  93         }
  94 }
  95
  96
  97 void OPJ_CALLCONV opj_image_header_destroy(opj_image_header_t *image) {
  98         if(image) {
  99                 if(image->comps) {
 100                         /* image components */
 101                         opj_free(image->comps);
 102                 }
 103                 opj_free(image);
 104         }
 105 }
 106
 107 /**
 108  * Updates the components characteristics of the image from the coding parameters.
 109  *
 110  * @param p_image_header        the image header to update.
 111  * @param p_cp                          the coding parameters from which to update the image.
 112  */
 113 void opj_image_comp_header_update(opj_image_header_t * p_image_header, const opj_cp_v2_t * p_cp)
 114 {
 115         OPJ_UINT32 i, l_width, l_height;
 116         OPJ_INT32 l_x0, l_y0, l_x1, l_y1;
 117         OPJ_INT32 l_comp_x0, l_comp_y0, l_comp_x1, l_comp_y1;
 118         opj_image_comp_header_t* l_img_comp = NULL;
 119
 120         l_x0 = int_max(p_cp->tx0 , p_image_header->x0);
 121         l_y0 = int_max(p_cp->ty0 , p_image_header->y0);
 122         l_x1 = int_min(p_cp->tx0 + p_cp->tw * p_cp->tdx, p_image_header->x1);
 123         l_y1 = int_min(p_cp->ty0 + p_cp->th * p_cp->tdy, p_image_header->y1);
 124
 125         l_img_comp = p_image_header->comps;
 126         for     (i = 0; i < p_image_header->numcomps; ++i) {
 127                 l_comp_x0 = int_ceildiv(l_x0, l_img_comp->dx);
 128                 l_comp_y0 = int_ceildiv(l_y0, l_img_comp->dy);
 129                 l_comp_x1 = int_ceildiv(l_x1, l_img_comp->dx);
 130                 l_comp_y1 = int_ceildiv(l_y1, l_img_comp->dy);
 131                 l_width = int_ceildivpow2(l_comp_x1 - l_comp_x0, l_img_comp->factor);
 132                 l_height = int_ceildivpow2(l_comp_y1 - l_comp_y0, l_img_comp->factor);
 133                 l_img_comp->w = l_width;
 134                 l_img_comp->h = l_height;
 135                 l_img_comp->x0 = l_x0;
 136                 l_img_comp->y0 = l_y0;
 137                 ++l_img_comp;
 138         }
 139 }
 140