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 "image.h"
  24 #include "imagei.h"
  25
  26 static int i_ppix_ddoub(i_img *im, int x, int y, 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, i_color *vals);
  30 static int i_ppixf_ddoub(i_img *im, int x, int y, 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, 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 Creates a new double per sample image.
  85
  86 =cut
  87 */
  88 i_img *i_img_double_new_low(i_img *im, int x, int y, int ch) {
  89   mm_log((1,"i_img_double_new(x %d, y %d, ch %d)\n", x, y, ch));
  90
  91   if (x < 1 || y < 1) {
  92     i_push_error(0, "Image sizes must be positive");
  93     return NULL;
  94   }
  95   if (ch < 1 || ch > MAXCHANNELS) {
  96     i_push_errorf(0, "channels must be between 1 and %d", MAXCHANNELS);
  97     return NULL;
  98   }
  99
 100   *im = IIM_base_double_direct;
 101   i_tags_new(&im->tags);
 102   im->xsize = x;
 103   im->ysize = y;
 104   im->channels = ch;
 105   im->bytes = x * y * ch * sizeof(double);
 106   im->ext_data = NULL;
 107   im->idata = mymalloc(im->bytes);
 108   if (im->idata) {
 109     memset(im->idata, 0, im->bytes);
 110   }
 111   else {
 112     i_tags_destroy(&im->tags);
 113     im = NULL;
 114   }
 115
 116   return im;
 117 }
 118
 119 i_img *i_img_double_new(int x, int y, int ch) {
 120   i_img *im;
 121
 122   i_clear_error();
 123
 124   im = mymalloc(sizeof(i_img));
 125   if (im) {
 126     if (!i_img_double_new_low(im, x, y, ch)) {
 127       myfree(im);
 128       im = NULL;
 129     }
 130   }
 131
 132   mm_log((1, "(%p) <- i_img_double_new\n", im));
 133
 134   return im;
 135 }
 136
 137 static int i_ppix_ddoub(i_img *im, int x, int y, i_color *val) {
 138   int off, ch;
 139
 140   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 141     return -1;
 142
 143   off = (x + y * im->xsize) * im->channels;
 144   for (ch = 0; ch < im->channels; ++ch)
 145     ((double*)im->idata)[off+ch] = Sample8ToF(val->channel[ch]);
 146
 147   return 0;
 148 }
 149
 150 static int i_gpix_ddoub(i_img *im, int x, int y, i_color *val) {
 151   int off, ch;
 152
 153   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 154     return -1;
 155
 156   off = (x + y * im->xsize) * im->channels;
 157   for (ch = 0; ch < im->channels; ++ch)
 158     val->channel[ch] = SampleFTo8(((double *)im->idata)[off+ch]);
 159
 160   return 0;
 161 }
 162
 163 static int i_ppixf_ddoub(i_img *im, int x, int y, i_fcolor *val) {
 164   int off, ch;
 165
 166   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 167     return -1;
 168
 169   off = (x + y * im->xsize) * im->channels;
 170   for (ch = 0; ch < im->channels; ++ch)
 171     ((double *)im->idata)[off+ch] = val->channel[ch];;
 172
 173   return 0;
 174 }
 175
 176 static int i_gpixf_ddoub(i_img *im, int x, int y, i_fcolor *val) {
 177   int off, ch;
 178
 179   if (x < 0 || x >= im->xsize || y < 0 || y > im->ysize)
 180     return -1;
 181
 182   off = (x + y * im->xsize) * im->channels;
 183   for (ch = 0; ch < im->channels; ++ch)
 184     val->channel[ch] = ((double *)im->idata)[off+ch];
 185
 186   return 0;
 187 }
 188
 189 static int i_glin_ddoub(i_img *im, int l, int r, int y, i_color *vals) {
 190   int ch, count, i;
 191   int off;
 192   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 193     if (r > im->xsize)
 194       r = im->xsize;
 195     off = (l+y*im->xsize) * im->channels;
 196     count = r - l;
 197     for (i = 0; i < count; ++i) {
 198       for (ch = 0; ch < im->channels; ++ch) {
 199         vals[i].channel[ch] = SampleFTo8(((double *)im->idata)[off]);
 200         ++off;
 201       }
 202     }
 203     return count;
 204   }
 205   else {
 206     return 0;
 207   }
 208 }
 209
 210 static int i_plin_ddoub(i_img *im, int l, int r, int y, i_color *vals) {
 211   int ch, count, i;
 212   int off;
 213   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 214     if (r > im->xsize)
 215       r = im->xsize;
 216     off = (l+y*im->xsize) * im->channels;
 217     count = r - l;
 218     for (i = 0; i < count; ++i) {
 219       for (ch = 0; ch < im->channels; ++ch) {
 220         ((double *)im->idata)[off] = Sample8ToF(vals[i].channel[ch]);
 221         ++off;
 222       }
 223     }
 224     return count;
 225   }
 226   else {
 227     return 0;
 228   }
 229 }
 230
 231 static int i_glinf_ddoub(i_img *im, int l, int r, int y, i_fcolor *vals) {
 232   int ch, count, i;
 233   int off;
 234   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 235     if (r > im->xsize)
 236       r = im->xsize;
 237     off = (l+y*im->xsize) * im->channels;
 238     count = r - l;
 239     for (i = 0; i < count; ++i) {
 240       for (ch = 0; ch < im->channels; ++ch) {
 241         vals[i].channel[ch] = ((double *)im->idata)[off];
 242         ++off;
 243       }
 244     }
 245     return count;
 246   }
 247   else {
 248     return 0;
 249   }
 250 }
 251
 252 static int i_plinf_ddoub(i_img *im, int l, int r, int y, i_fcolor *vals) {
 253   int ch, count, i;
 254   int off;
 255   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 256     if (r > im->xsize)
 257       r = im->xsize;
 258     off = (l+y*im->xsize) * im->channels;
 259     count = r - l;
 260     for (i = 0; i < count; ++i) {
 261       for (ch = 0; ch < im->channels; ++ch) {
 262         ((double *)im->idata)[off] = vals[i].channel[ch];
 263         ++off;
 264       }
 265     }
 266     return count;
 267   }
 268   else {
 269     return 0;
 270   }
 271 }
 272
 273 static int i_gsamp_ddoub(i_img *im, int l, int r, int y, i_sample_t *samps,
 274                        int const *chans, int chan_count) {
 275   int ch, count, i, w;
 276   int off;
 277
 278   if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
 279     if (r > im->xsize)
 280       r = im->xsize;
 281     off = (l+y*im->xsize) * im->channels;
 282     w = r - l;
 283     count = 0;
 284
 285     if (chans) {
 286       /* make sure we have good channel numbers */
 287       for (ch = 0; ch < chan_count; ++ch) {
 288         if (chans[ch] < 0 || chans[ch] >= im->channels) {
 289           i_push_errorf(0, "No channel %d in this image", chans[ch]);
 290           return 0;
 291         }
 292       }
 293       for (i = 0; i < w; ++i) {
 294         for (ch = 0; ch < chan_count; ++ch) {
 295           *samps++ = SampleFTo8(((double *)im->idata)[off+chans[ch]]);
 296           ++count;
 297         }
 298         off += im->channels;
 299       }
 300     }
 301     else {
 302       for (i = 0; i < w; ++i) {
 303         for (ch = 0; ch < chan_count; ++ch) {
 304           *samps++ = SampleFTo8(((double *)im->idata)[off+ch]);
 305           ++count;
 306         }
 307         off += im->channels;
 308       }
 309     }
 310
 311     return count;
 312   }
 313   else {
 314     return 0;
 315   }
 316 }
 317
 318 static int i_gsampf_ddoub(i_img *im, int l, int r, int y, i_fsample_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++ = ((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++ = ((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
 364 /*
 365 =back
 366
 367 =head1 AUTHOR
 368
 369 Tony Cook <tony@develop-help.com>
 370
 371 =head1 SEE ALSO
 372
 373 Imager(3)
 374
 375 =cut
 376 */