imgdouble.c

   1 /*
   2 =head1 NAME
   3
   4 imgdouble.c - implements double per sample images
   5
   6 =head1 SYNOPSIS
   7
   8   i_img *im = i_img_double_new(int x, int y, int channels);
   9   # use like a normal image
  10
  11 =head1 DESCRIPTION
  12
  13 Implements double/sample images.
  14
  15 This basic implementation is required so that we have some larger
  16 sample image type to work with.
  17
  18 =over
  19
  20 =cut
  21 */
  22
  23 #include "imager.h"
  24 #include "imageri.h"
  25
  26 static int i_ppix_ddoub(i_img *im, int x, int y, const i_color *val);
  27 static int i_gpix_ddoub(i_img *im, int x, int y, i_color *val);
  28 static int i_glin_ddoub(i_img *im, int l, int r, int y, i_color *vals);
  29 static int i_plin_ddoub(i_img *im, int l, int r, int y, const i_color *vals);
  30 static int i_ppixf_ddoub(i_img *im, int x, int y, const i_fcolor *val);
  31 static int i_gpixf_ddoub(i_img *im, int x, int y, i_fcolor *val);
  32 static int i_glinf_ddoub(i_img *im, int l, int r, int y, i_fcolor *vals);
  33 static int i_plinf_ddoub(i_img *im, int l, int r, int y, const i_fcolor *vals);
  34 static int i_gsamp_ddoub(i_img *im, int l, int r, int y, i_sample_t *samps,
  35                        int const *chans, int chan_count);
  36 static int i_gsampf_ddoub(i_img *im, int l, int r, int y, i_fsample_t *samps,
  37                         int const *chans, int chan_count);
  38
  39 /*
  40 =item IIM_base_16bit_direct
  41
  42 Base structure used to initialize a 16-bit/sample image.
  43
  44 Internal.
  45
  46 =cut
  47 */
  48 static i_img IIM_base_double_direct =
  49 {
  50   0, /* channels set */
  51   0, 0, 0, /* xsize, ysize, bytes */
  52   ~0U, /* ch_mask */
  53   i_double_bits, /* bits */
  54   i_direct_type, /* type */
  55   0, /* virtual */
  56   NULL, /* idata */
  57   { 0, 0, NULL }, /* tags */
  58   NULL, /* ext_data */
  59
  60   i_ppix_ddoub, /* i_f_ppix */
  61   i_ppixf_ddoub, /* i_f_ppixf */
  62   i_plin_ddoub, /* i_f_plin */
  63   i_plinf_ddoub, /* i_f_plinf */
  64   i_gpix_ddoub, /* i_f_gpix */
  65   i_gpixf_ddoub, /* i_f_gpixf */
  66   i_glin_ddoub, /* i_f_glin */
  67   i_glinf_ddoub, /* i_f_glinf */
  68   i_gsamp_ddoub, /* i_f_gsamp */
  69   i_gsampf_ddoub, /* i_f_gsampf */
  70
  71   NULL, /* i_f_gpal */
  72   NULL, /* i_f_ppal */
  73   NULL, /* i_f_addcolor */
  74   NULL, /* i_f_getcolor */
  75   NULL, /* i_f_colorcount */
  76   NULL, /* i_f_findcolor */
  77
  78   NULL, /* i_f_destroy */
  79 };
  80
  81 /*
  82 =item i_img_double_new(int x, int y, int ch)
  83
  84 =category Image creation
  85
  86 Creates a new double per sample image.
  87
  88 =cut
  89 */
  90 i_img *i_img_double_new_low(i_img *im, int x, int y, int ch) {
  91   int bytes;
  92
  93   mm_log((1,"i_img_double_new(x %d, y %d, ch %d)\n", x, y, ch));
  94
  95   if (x < 1 || y < 1) {
  96     i_push_error(0, "Image sizes must be positive");
  97     return NULL;
  98   }
  99   if (ch < 1 || ch > MAXCHANNELS) {
 100     i_push_errorf(0, "channels must be between 1 and %d", MAXCHANNELS);
 101     return NULL;
 102   }
 103   bytes = x * y * ch * sizeof(double);
 104   if (bytes / y / ch / sizeof(double) != x) {
 105     i_push_errorf(0, "integer overflow calculating image allocation");
 106     return NULL;
 107   }
 108
 109   *im = IIM_base_double_direct;
 110   i_tags_new(&im->tags);
 111   im->xsize = x;
 112   im->ysize = y;
 113   im->channels = ch;
 114   im->bytes = bytes;
 115   im->ext_data = NULL;
 116   im->idata = mymalloc(im->bytes);
 117   if (im->idata) {
 118     memset(im->idata, 0, im->bytes);
 119   }
 120   else {
 121     i_tags_destroy(&im->tags);
 122     im = NULL;
 123   }
 124
 125   return im;
 126 }
 127
 128 i_img *i_img_double_new(int x, int y, int ch) {
 129   i_img *im;
 130
 131   i_clear_error();
 132
 133   im = mymalloc(sizeof(i_img));
 134   if (im) {
 135     if (!i_img_double_new_low(im, x, y, ch)) {
 136       myfree(im);
 137       im = NULL;
 138     }
 139   }
 140
 141   mm_log((1, "(%p) <- i_img_double_new\n", im));
 142
 143   return im;
 144 }
 145
 146 static int i_ppix_ddoub(i_img *im, int x, int y, const i_color *val) {
 147   int off, ch;
 148
 149   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 150     return -1;
 151
 152   off = (x + y * im->xsize) * im->channels;
 153   if (I_ALL_CHANNELS_WRITABLE(im)) {
 154     for (ch = 0; ch < im->channels; ++ch)
 155       ((double*)im->idata)[off+ch] = Sample8ToF(val->channel[ch]);
 156   }
 157   else {
 158     for (ch = 0; ch < im->channels; ++ch)
 159       if (im->ch_mask & (1<<ch))
 160         ((double*)im->idata)[off+ch] = Sample8ToF(val->channel[ch]);
 161   }
 162
 163   return 0;
 164 }
 165
 166 static int i_gpix_ddoub(i_img *im, int x, int y, i_color *val) {
 167   int off, ch;
 168
 169   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 170     return -1;
 171
 172   off = (x + y * im->xsize) * im->channels;
 173   for (ch = 0; ch < im->channels; ++ch)
 174     val->channel[ch] = SampleFTo8(((double *)im->idata)[off+ch]);
 175
 176   return 0;
 177 }
 178
 179 static int i_ppixf_ddoub(i_img *im, int x, int y, const i_fcolor *val) {
 180   int off, ch;
 181
 182   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 183     return -1;
 184
 185   off = (x + y * im->xsize) * im->channels;
 186   if (I_ALL_CHANNELS_WRITABLE(im)) {
 187     for (ch = 0; ch < im->channels; ++ch)
 188       ((double *)im->idata)[off+ch] = val->channel[ch];
 189   }
 190   else {
 191     for (ch = 0; ch < im->channels; ++ch)
 192       if (im->ch_mask & (1 << ch))
 193         ((double *)im->idata)[off+ch] = val->channel[ch];
 194   }
 195
 196   return 0;
 197 }
 198
 199 static int i_gpixf_ddoub(i_img *im, int x, int y, i_fcolor *val) {
 200   int off, ch;
 201
 202   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 203     return -1;
 204
 205   off = (x + y * im->xsize) * im->channels;
 206   for (ch = 0; ch < im->channels; ++ch)
 207     val->channel[ch] = ((double *)im->idata)[off+ch];
 208
 209   return 0;
 210 }
 211
 212 static int i_glin_ddoub(i_img *im, int l, int r, int y, i_color *vals) {
 213   int ch, count, i;
 214   int off;
 215   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 216     if (r > im->xsize)
 217       r = im->xsize;
 218     off = (l+y*im->xsize) * im->channels;
 219     count = r - l;
 220     for (i = 0; i < count; ++i) {
 221       for (ch = 0; ch < im->channels; ++ch) {
 222         vals[i].channel[ch] = SampleFTo8(((double *)im->idata)[off]);
 223         ++off;
 224       }
 225     }
 226     return count;
 227   }
 228   else {
 229     return 0;
 230   }
 231 }
 232
 233 static int i_plin_ddoub(i_img *im, int l, int r, int y, const i_color *vals) {
 234   int ch, count, i;
 235   int off;
 236   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 237     if (r > im->xsize)
 238       r = im->xsize;
 239     off = (l+y*im->xsize) * im->channels;
 240     count = r - l;
 241     if (I_ALL_CHANNELS_WRITABLE(im)) {
 242       for (i = 0; i < count; ++i) {
 243         for (ch = 0; ch < im->channels; ++ch) {
 244           ((double *)im->idata)[off] = Sample8ToF(vals[i].channel[ch]);
 245           ++off;
 246         }
 247       }
 248     }
 249     else {
 250       for (i = 0; i < count; ++i) {
 251         for (ch = 0; ch < im->channels; ++ch) {
 252           if (im->ch_mask & (1 << ch))
 253             ((double *)im->idata)[off] = Sample8ToF(vals[i].channel[ch]);
 254           ++off;
 255         }
 256       }
 257     }
 258     return count;
 259   }
 260   else {
 261     return 0;
 262   }
 263 }
 264
 265 static int i_glinf_ddoub(i_img *im, int l, int r, int y, i_fcolor *vals) {
 266   int ch, count, i;
 267   int off;
 268   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 269     if (r > im->xsize)
 270       r = im->xsize;
 271     off = (l+y*im->xsize) * im->channels;
 272     count = r - l;
 273     for (i = 0; i < count; ++i) {
 274       for (ch = 0; ch < im->channels; ++ch) {
 275         vals[i].channel[ch] = ((double *)im->idata)[off];
 276         ++off;
 277       }
 278     }
 279     return count;
 280   }
 281   else {
 282     return 0;
 283   }
 284 }
 285
 286 static int i_plinf_ddoub(i_img *im, int l, int r, int y, const i_fcolor *vals) {
 287   int ch, count, i;
 288   int off;
 289   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 290     if (r > im->xsize)
 291       r = im->xsize;
 292     off = (l+y*im->xsize) * im->channels;
 293     count = r - l;
 294     if (I_ALL_CHANNELS_WRITABLE(im)) {
 295       for (i = 0; i < count; ++i) {
 296         for (ch = 0; ch < im->channels; ++ch) {
 297           ((double *)im->idata)[off] = vals[i].channel[ch];
 298           ++off;
 299         }
 300       }
 301     }
 302     else {
 303       for (i = 0; i < count; ++i) {
 304         for (ch = 0; ch < im->channels; ++ch) {
 305           if (im->ch_mask & (1 << ch))
 306             ((double *)im->idata)[off] = vals[i].channel[ch];
 307           ++off;
 308         }
 309       }
 310     }
 311     return count;
 312   }
 313   else {
 314     return 0;
 315   }
 316 }
 317
 318 static int i_gsamp_ddoub(i_img *im, int l, int r, int y, i_sample_t *samps,
 319                        int const *chans, int chan_count) {
 320   int ch, count, i, w;
 321   int off;
 322
 323   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 324     if (r > im->xsize)
 325       r = im->xsize;
 326     off = (l+y*im->xsize) * im->channels;
 327     w = r - l;
 328     count = 0;
 329
 330     if (chans) {
 331       /* make sure we have good channel numbers */
 332       for (ch = 0; ch < chan_count; ++ch) {
 333         if (chans[ch] < 0 || chans[ch] >= im->channels) {
 334           i_push_errorf(0, "No channel %d in this image", chans[ch]);
 335           return 0;
 336         }
 337       }
 338       for (i = 0; i < w; ++i) {
 339         for (ch = 0; ch < chan_count; ++ch) {
 340           *samps++ = SampleFTo8(((double *)im->idata)[off+chans[ch]]);
 341           ++count;
 342         }
 343         off += im->channels;
 344       }
 345     }
 346     else {
 347       for (i = 0; i < w; ++i) {
 348         for (ch = 0; ch < chan_count; ++ch) {
 349           *samps++ = SampleFTo8(((double *)im->idata)[off+ch]);
 350           ++count;
 351         }
 352         off += im->channels;
 353       }
 354     }
 355
 356     return count;
 357   }
 358   else {
 359     return 0;
 360   }
 361 }
 362
 363 static int i_gsampf_ddoub(i_img *im, int l, int r, int y, i_fsample_t *samps,
 364                         int const *chans, int chan_count) {
 365   int ch, count, i, w;
 366   int off;
 367
 368   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 369     if (r > im->xsize)
 370       r = im->xsize;
 371     off = (l+y*im->xsize) * im->channels;
 372     w = r - l;
 373     count = 0;
 374
 375     if (chans) {
 376       /* make sure we have good channel numbers */
 377       for (ch = 0; ch < chan_count; ++ch) {
 378         if (chans[ch] < 0 || chans[ch] >= im->channels) {
 379           i_push_errorf(0, "No channel %d in this image", chans[ch]);
 380           return 0;
 381         }
 382       }
 383       for (i = 0; i < w; ++i) {
 384         for (ch = 0; ch < chan_count; ++ch) {
 385           *samps++ = ((double *)im->idata)[off+chans[ch]];
 386           ++count;
 387         }
 388         off += im->channels;
 389       }
 390     }
 391     else {
 392       for (i = 0; i < w; ++i) {
 393         for (ch = 0; ch < chan_count; ++ch) {
 394           *samps++ = ((double *)im->idata)[off+ch];
 395           ++count;
 396         }
 397         off += im->channels;
 398       }
 399     }
 400
 401     return count;
 402   }
 403   else {
 404     return 0;
 405   }
 406 }
 407
 408
 409 /*
 410 =back
 411
 412 =head1 AUTHOR
 413
 414 Tony Cook <tony@develop-help.com>
 415
 416 =head1 SEE ALSO
 417
 418 Imager(3)
 419
 420 =cut
 421 */