=cut
*/
+#define IMAGER_NO_CONTEXT
#include "imager.h"
#include "imageri.h"
int const *chans, int chan_count);
static i_img_dim i_gsampf_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y, i_fsample_t *samps,
int const *chans, int chan_count);
+static i_img_dim
+i_psamp_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y, const i_sample_t *samps, const int *chans, int chan_count);
+static i_img_dim
+i_psampf_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y, const i_fsample_t *samps, const int *chans, int chan_count);
/*
=item IIM_base_16bit_direct
i_gsamp_bits_fb,
NULL, /* i_f_psamp_bits */
+
+ i_psamp_ddoub, /* i_f_psamp */
+ i_psampf_ddoub /* i_f_psampf */
};
/*
-=item i_img_double_new(i_img_dim x, i_img_dim y, int ch)
+=item im_img_double_new(ctx, x, y, ch)
+X<im_img_double_new API>X<i_img_double_new API>
=category Image creation/destruction
+=synopsis i_img *img = im_img_double_new(aIMCTX, width, height, channels);
=synopsis i_img *img = i_img_double_new(width, height, channels);
Creates a new double per sample image.
+Also callable as C<i_img_double_new(width, height, channels)>.
+
=cut
*/
-i_img *i_img_double_new(i_img_dim x, i_img_dim y, int ch) {
+i_img *
+im_img_double_new(pIMCTX, i_img_dim x, i_img_dim y, int ch) {
size_t bytes;
i_img *im;
- mm_log((1,"i_img_double_new(x %" i_DF ", y %" i_DF ", ch %d)\n",
+ im_log((aIMCTX, 1,"i_img_double_new(x %" i_DF ", y %" i_DF ", ch %d)\n",
i_DFc(x), i_DFc(y), ch));
if (x < 1 || y < 1) {
- i_push_error(0, "Image sizes must be positive");
+ im_push_error(aIMCTX, 0, "Image sizes must be positive");
return NULL;
}
if (ch < 1 || ch > MAXCHANNELS) {
- i_push_errorf(0, "channels must be between 1 and %d", MAXCHANNELS);
+ im_push_errorf(aIMCTX, 0, "channels must be between 1 and %d", MAXCHANNELS);
return NULL;
}
bytes = x * y * ch * sizeof(double);
if (bytes / y / ch / sizeof(double) != x) {
- i_push_errorf(0, "integer overflow calculating image allocation");
+ im_push_errorf(aIMCTX, 0, "integer overflow calculating image allocation");
return NULL;
}
- im = i_img_alloc();
+ im = im_img_alloc(aIMCTX);
*im = IIM_base_double_direct;
i_tags_new(&im->tags);
im->xsize = x;
im->ext_data = NULL;
im->idata = mymalloc(im->bytes);
memset(im->idata, 0, im->bytes);
- i_img_init(im);
+ im_img_init(aIMCTX, im);
return im;
}
/* make sure we have good channel numbers */
for (ch = 0; ch < chan_count; ++ch) {
if (chans[ch] < 0 || chans[ch] >= im->channels) {
- i_push_errorf(0, "No channel %d in this image", chans[ch]);
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "No channel %d in this image", chans[ch]);
return 0;
}
}
}
else {
if (chan_count <= 0 || chan_count > im->channels) {
- i_push_errorf(0, "chan_count %d out of range, must be >0, <= channels",
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "chan_count %d out of range, must be >0, <= channels",
chan_count);
return 0;
}
/* make sure we have good channel numbers */
for (ch = 0; ch < chan_count; ++ch) {
if (chans[ch] < 0 || chans[ch] >= im->channels) {
- i_push_errorf(0, "No channel %d in this image", chans[ch]);
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "No channel %d in this image", chans[ch]);
return 0;
}
}
}
else {
if (chan_count <= 0 || chan_count > im->channels) {
- i_push_errorf(0, "chan_count %d out of range, must be >0, <= channels",
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "chan_count %d out of range, must be >0, <= channels",
chan_count);
return 0;
}
}
}
+/*
+=item i_psamp_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y, const i_sample_t *samps, int *chans, int chan_count)
+
+Writes sample values to im for the horizontal line (l, y) to (r-1,y)
+for the channels specified by chans, an array of int with chan_count
+elements.
+
+Returns the number of samples written (which should be (r-l) *
+bits_set(chan_mask)
+
+=cut
+*/
+
+static
+i_img_dim
+i_psamp_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y,
+ const i_sample_t *samps, const int *chans, int chan_count) {
+ int ch;
+ i_img_dim count, i, w;
+
+ if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
+ i_img_dim offset;
+ if (r > im->xsize)
+ r = im->xsize;
+ offset = (l+y*im->xsize) * im->channels;
+ w = r - l;
+ count = 0;
+
+ if (chans) {
+ /* make sure we have good channel numbers */
+ /* and test if all channels specified are in the mask */
+ int all_in_mask = 1;
+ for (ch = 0; ch < chan_count; ++ch) {
+ if (chans[ch] < 0 || chans[ch] >= im->channels) {
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "No channel %d in this image", chans[ch]);
+ return -1;
+ }
+ if (!((1 << chans[ch]) & im->ch_mask))
+ all_in_mask = 0;
+ }
+ if (all_in_mask) {
+ for (i = 0; i < w; ++i) {
+ for (ch = 0; ch < chan_count; ++ch) {
+ ((double*)im->idata)[offset + chans[ch]] = Sample8ToF(*samps);
+ ++samps;
+ ++count;
+ }
+ offset += im->channels;
+ }
+ }
+ else {
+ for (i = 0; i < w; ++i) {
+ for (ch = 0; ch < chan_count; ++ch) {
+ if (im->ch_mask & (1 << (chans[ch])))
+ ((double*)im->idata)[offset + chans[ch]] = Sample8ToF(*samps);
+
+ ++samps;
+ ++count;
+ }
+ offset += im->channels;
+ }
+ }
+ }
+ else {
+ if (chan_count <= 0 || chan_count > im->channels) {
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "chan_count %d out of range, must be >0, <= channels",
+ chan_count);
+ return -1;
+ }
+ for (i = 0; i < w; ++i) {
+ unsigned mask = 1;
+ for (ch = 0; ch < chan_count; ++ch) {
+ if (im->ch_mask & mask)
+ ((double*)im->idata)[offset + ch] = Sample8ToF(*samps);
+
+ ++samps;
+ ++count;
+ mask <<= 1;
+ }
+ offset += im->channels;
+ }
+ }
+
+ return count;
+ }
+ else {
+ dIMCTXim(im);
+ i_push_error(0, "Image position outside of image");
+ return -1;
+ }
+}
+
+/*
+=item i_psampf_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y, const i_fsample_t *samps, int *chans, int chan_count)
+
+Writes sample values to im for the horizontal line (l, y) to (r-1,y)
+for the channels specified by chans, an array of int with chan_count
+elements.
+
+Returns the number of samples written (which should be (r-l) *
+bits_set(chan_mask)
+
+=cut
+*/
+
+static
+i_img_dim
+i_psampf_ddoub(i_img *im, i_img_dim l, i_img_dim r, i_img_dim y,
+ const i_fsample_t *samps, const int *chans, int chan_count) {
+ int ch;
+ i_img_dim count, i, w;
+
+ if (y >=0 && y < im->ysize && l < im->xsize && l >= 0) {
+ i_img_dim offset;
+ if (r > im->xsize)
+ r = im->xsize;
+ offset = (l+y*im->xsize) * im->channels;
+ w = r - l;
+ count = 0;
+
+ if (chans) {
+ /* make sure we have good channel numbers */
+ /* and test if all channels specified are in the mask */
+ int all_in_mask = 1;
+ for (ch = 0; ch < chan_count; ++ch) {
+ if (chans[ch] < 0 || chans[ch] >= im->channels) {
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "No channel %d in this image", chans[ch]);
+ return -1;
+ }
+ if (!((1 << chans[ch]) & im->ch_mask))
+ all_in_mask = 0;
+ }
+ if (all_in_mask) {
+ for (i = 0; i < w; ++i) {
+ for (ch = 0; ch < chan_count; ++ch) {
+ ((double*)im->idata)[offset + chans[ch]] = *samps;
+ ++samps;
+ ++count;
+ }
+ offset += im->channels;
+ }
+ }
+ else {
+ for (i = 0; i < w; ++i) {
+ for (ch = 0; ch < chan_count; ++ch) {
+ if (im->ch_mask & (1 << (chans[ch])))
+ ((double*)im->idata)[offset + chans[ch]] = *samps;
+
+ ++samps;
+ ++count;
+ }
+ offset += im->channels;
+ }
+ }
+ }
+ else {
+ if (chan_count <= 0 || chan_count > im->channels) {
+ dIMCTXim(im);
+ im_push_errorf(aIMCTX, 0, "chan_count %d out of range, must be >0, <= channels",
+ chan_count);
+ return -1;
+ }
+ for (i = 0; i < w; ++i) {
+ unsigned mask = 1;
+ for (ch = 0; ch < chan_count; ++ch) {
+ if (im->ch_mask & mask)
+ ((double*)im->idata)[offset + ch] = *samps;
+
+ ++samps;
+ ++count;
+ mask <<= 1;
+ }
+ offset += im->channels;
+ }
+ }
+
+ return count;
+ }
+ else {
+ dIMCTXim(im);
+ i_push_error(0, "Image position outside of image");
+ return -1;
+ }
+}
+
/*
=item i_img_to_drgb(im)
i_img *targ;
i_fcolor *line;
i_img_dim y;
+ dIMCTXim(im);
- targ = i_img_double_new(im->xsize, im->ysize, im->channels);
+ targ = im_img_double_new(aIMCTX, im->xsize, im->ysize, im->channels);
if (!targ)
return NULL;
line = mymalloc(sizeof(i_fcolor) * im->xsize);