static void
convertToPamPnm(struct pam * const outpamP,
jas_image_t * const jasperP,
int const jasperCmptNo[]) {
jas_matrix_t ** matrix; /* malloc'ed */
/* matrix[X] is the data for Plane X of the current row */
sample * jasperMaxval;
unsigned int row;
tuple * tuplerow;
jas_seqent_t ** jasperRow; /* malloc'ed */
/* A row of a plane of the raster from the Jasper library
This is an array of pointers into the 'matrix' data structures.
*/
bool singleMaxval;
createMatrices(outpamP, &matrix);
computeComponentMaxval(outpamP, jasperP, jasperCmptNo,
&jasperMaxval, &singleMaxval);
MALLOCARRAY(jasperRow, outpamP->depth);
if (jasperRow == NULL)
pm_error("Out of memory");
tuplerow = pnm_allocpamrow(outpamP);
for (row = 0; row < outpamP->height; ++row) {
unsigned int plane;
for (plane = 0; plane < outpamP->depth; ++plane) {
int rc;
rc = jas_image_readcmpt(jasperP, jasperCmptNo[plane], 0, row,
outpamP->width, 1,
matrix[plane]);
if (rc != 0)
pm_error("jas_image_readcmpt() of row %u plane %u "
"failed.",
row, plane);
jasperRow[plane] = jas_matrix_getref(matrix[plane], 0, 0);
}
if (singleMaxval)
copyRowSingleMaxval(jasperRow, tuplerow, outpamP);
else
copyRowAnyMaxval(jasperRow, tuplerow, outpamP, jasperMaxval);
pnm_writepamrow(outpamP, tuplerow);
}
pnm_freepamrow(tuplerow);
destroyMatrices(outpamP, matrix);
free(jasperRow);
free(jasperMaxval);
}
static int jpc_encclnpass(jpc_mqenc_t *mqenc, int bitpos, int orient, int vcausalflag, int segsymflag, jas_matrix_t *flags,
jas_matrix_t *data, int term, long *nmsedec)
{
int i;
int j;
int k;
int vscanlen;
int v;
int runlen;
jpc_fix_t *fp;
int width;
int height;
jpc_fix_t *dp;
int one;
int frowstep;
int drowstep;
int fstripestep;
int dstripestep;
jpc_fix_t *fstripestart;
jpc_fix_t *dstripestart;
jpc_fix_t *fvscanstart;
jpc_fix_t *dvscanstart;
*nmsedec = 0;
width = jas_matrix_numcols(data);
height = jas_matrix_numrows(data);
frowstep = jas_matrix_rowstep(flags);
drowstep = jas_matrix_rowstep(data);
fstripestep = frowstep << 2;
dstripestep = drowstep << 2;
one = 1 << (bitpos + JPC_NUMEXTRABITS);
fstripestart = jas_matrix_getref(flags, 1, 1);
dstripestart = jas_matrix_getref(data, 0, 0);
for (i = height; i > 0; i -= 4, fstripestart += fstripestep,
dstripestart += dstripestep) {
fvscanstart = fstripestart;
dvscanstart = dstripestart;
vscanlen = JAS_MIN(i, 4);
for (j = width; j > 0; --j, ++fvscanstart, ++dvscanstart) {
fp = fvscanstart;
if (vscanlen >= 4 && !((*fp) & (JPC_SIG | JPC_VISIT |
JPC_OTHSIGMSK)) && (fp += frowstep, !((*fp) & (JPC_SIG |
JPC_VISIT | JPC_OTHSIGMSK))) && (fp += frowstep, !((*fp) &
(JPC_SIG | JPC_VISIT | JPC_OTHSIGMSK))) && (fp += frowstep,
!((*fp) & (JPC_SIG | JPC_VISIT | JPC_OTHSIGMSK)))) {
dp = dvscanstart;
for (k = 0; k < vscanlen; ++k) {
v = (abs(*dp) & one) ? 1 : 0;
if (v) {
break;
}
dp += drowstep;
}
runlen = k;
if (runlen >= 4) {
jpc_mqenc_setcurctx(mqenc, JPC_AGGCTXNO);
jpc_mqenc_putbit(mqenc, 0);
continue;
}
jpc_mqenc_setcurctx(mqenc, JPC_AGGCTXNO);
jpc_mqenc_putbit(mqenc, 1);
jpc_mqenc_setcurctx(mqenc, JPC_UCTXNO);
jpc_mqenc_putbit(mqenc, runlen >> 1);
jpc_mqenc_putbit(mqenc, runlen & 1);
fp = fvscanstart + frowstep * runlen;
dp = dvscanstart + drowstep * runlen;
k = vscanlen - runlen;
switch (runlen) {
case 0:
goto clnpass_partial0;
break;
case 1:
goto clnpass_partial1;
break;
case 2:
goto clnpass_partial2;
break;
case 3:
goto clnpass_partial3;
break;
}
} else {
runlen = 0;
fp = fvscanstart;
dp = dvscanstart;
k = vscanlen;
goto clnpass_full0;
}
clnpass_step(fp, frowstep, dp, bitpos, one,
orient, nmsedec, mqenc, clnpass_full0:, clnpass_partial0:, vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
clnpass_step(fp, frowstep, dp, bitpos, one,
orient, nmsedec, mqenc, ;, clnpass_partial1:, 0);
static int jpc_encrawrefpass(jpc_bitstream_t *out, int bitpos, int vcausalflag, jas_matrix_t *flags,
jas_matrix_t *data, int term, long *nmsedec)
{
int i;
int j;
int k;
int one;
int vscanlen;
int width;
int height;
int frowstep;
int drowstep;
int fstripestep;
int dstripestep;
jpc_fix_t *fstripestart;
jpc_fix_t *dstripestart;
jpc_fix_t *fvscanstart;
jpc_fix_t *dvscanstart;
jpc_fix_t *dp;
jpc_fix_t *fp;
*nmsedec = 0;
width = jas_matrix_numcols(data);
height = jas_matrix_numrows(data);
frowstep = jas_matrix_rowstep(flags);
drowstep = jas_matrix_rowstep(data);
fstripestep = frowstep << 2;
dstripestep = drowstep << 2;
one = 1 << (bitpos + JPC_NUMEXTRABITS);
fstripestart = jas_matrix_getref(flags, 1, 1);
dstripestart = jas_matrix_getref(data, 0, 0);
for (i = height; i > 0; i -= 4, fstripestart += fstripestep,
dstripestart += dstripestep) {
fvscanstart = fstripestart;
dvscanstart = dstripestart;
vscanlen = JAS_MIN(i, 4);
for (j = width; j > 0; --j, ++fvscanstart, ++dvscanstart) {
fp = fvscanstart;
dp = dvscanstart;
k = vscanlen;
rawrefpass_step(fp, dp, bitpos, one, nmsedec,
out, vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
rawrefpass_step(fp, dp, bitpos, one, nmsedec,
out, vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
rawrefpass_step(fp, dp, bitpos, one, nmsedec,
out, vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
rawrefpass_step(fp, dp, bitpos, one, nmsedec,
out, vcausalflag);
}
}
if (term) {
jpc_bitstream_outalign(out, 0x2a);
}
return 0;
}
static int dec_clnpass(jpc_dec_t *dec, register jpc_mqdec_t *mqdec, int bitpos, int orient,
int vcausalflag, int segsymflag, jas_matrix_t *flags, jas_matrix_t *data)
{
int i;
int j;
int k;
int vscanlen;
int v;
int half;
int runlen;
int f;
int width;
int height;
int one;
int oneplushalf;
jpc_fix_t *fp;
int frowstep;
int fstripestep;
jpc_fix_t *fstripestart;
jpc_fix_t *fvscanstart;
jpc_fix_t *dp;
int drowstep;
int dstripestep;
jpc_fix_t *dstripestart;
jpc_fix_t *dvscanstart;
/* Avoid compiler warning about unused parameters. */
dec = 0;
one = 1 << bitpos;
half = one >> 1;
oneplushalf = one | half;
width = jas_matrix_numcols(data);
height = jas_matrix_numrows(data);
frowstep = jas_matrix_rowstep(flags);
drowstep = jas_matrix_rowstep(data);
fstripestep = frowstep << 2;
dstripestep = drowstep << 2;
fstripestart = jas_matrix_getref(flags, 1, 1);
dstripestart = jas_matrix_getref(data, 0, 0);
for (i = 0; i < height; i += 4, fstripestart += fstripestep,
dstripestart += dstripestep) {
fvscanstart = fstripestart;
dvscanstart = dstripestart;
vscanlen = JAS_MIN(4, height - i);
for (j = width; j > 0; --j, ++fvscanstart, ++dvscanstart) {
fp = fvscanstart;
if (vscanlen >= 4 && (!((*fp) & (JPC_SIG | JPC_VISIT |
JPC_OTHSIGMSK))) && (fp += frowstep, !((*fp) & (JPC_SIG |
JPC_VISIT | JPC_OTHSIGMSK))) && (fp += frowstep, !((*fp) &
(JPC_SIG | JPC_VISIT | JPC_OTHSIGMSK))) && (fp += frowstep,
!((*fp) & (JPC_SIG | JPC_VISIT | JPC_OTHSIGMSK)))) {
jpc_mqdec_setcurctx(mqdec, JPC_AGGCTXNO);
JPC_T1D_GETBIT(mqdec, v, "CLN", "AGG");
if (!v) {
continue;
}
jpc_mqdec_setcurctx(mqdec, JPC_UCTXNO);
JPC_T1D_GETBITNOSKEW(mqdec, v, "CLN", "RL");
runlen = v;
JPC_T1D_GETBITNOSKEW(mqdec, v, "CLN", "RL");
runlen = (runlen << 1) | v;
f = *(fp = fvscanstart + frowstep * runlen);
dp = dvscanstart + drowstep * runlen;
k = vscanlen - runlen;
switch (runlen) {
case 0:
goto clnpass_partial0;
break;
case 1:
goto clnpass_partial1;
break;
case 2:
goto clnpass_partial2;
break;
case 3:
goto clnpass_partial3;
break;
}
} else {
f = *(fp = fvscanstart);
dp = dvscanstart;
k = vscanlen;
goto clnpass_full0;
}
/* Process first sample in vertical scan. */
jpc_clnpass_step(f, fp, frowstep, dp, oneplushalf, orient,
mqdec, clnpass_full0:, clnpass_partial0:,
vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process second sample in vertical scan. */
f = *fp;
jpc_clnpass_step(f, fp, frowstep, dp, oneplushalf, orient,
mqdec, ;, clnpass_partial1:, 0);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process third sample in vertical scan. */
f = *fp;
jpc_clnpass_step(f, fp, frowstep, dp, oneplushalf, orient,
mqdec, ;, clnpass_partial2:, 0);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process fourth sample in vertical scan. */
f = *fp;
jpc_clnpass_step(f, fp, frowstep, dp, oneplushalf, orient,
mqdec, ;, clnpass_partial3:, 0);
static int dec_rawrefpass(jpc_dec_t *dec, jpc_bitstream_t *in, int bitpos, int vcausalflag,
jas_matrix_t *flags, jas_matrix_t *data)
{
int i;
int j;
int k;
int vscanlen;
int width;
int height;
int one;
int poshalf;
int neghalf;
jpc_fix_t *fp;
int frowstep;
int fstripestep;
jpc_fix_t *fstripestart;
jpc_fix_t *fvscanstart;
jpc_fix_t *dp;
int drowstep;
int dstripestep;
jpc_fix_t *dstripestart;
jpc_fix_t *dvscanstart;
/* Avoid compiler warning about unused parameters. */
dec = 0;
vcausalflag = 0;
width = jas_matrix_numcols(data);
height = jas_matrix_numrows(data);
frowstep = jas_matrix_rowstep(flags);
drowstep = jas_matrix_rowstep(data);
fstripestep = frowstep << 2;
dstripestep = drowstep << 2;
one = 1 << bitpos;
poshalf = one >> 1;
neghalf = (bitpos > 0) ? (-poshalf) : (-1);
fstripestart = jas_matrix_getref(flags, 1, 1);
dstripestart = jas_matrix_getref(data, 0, 0);
for (i = height; i > 0; i -= 4, fstripestart += fstripestep,
dstripestart += dstripestep) {
fvscanstart = fstripestart;
dvscanstart = dstripestart;
vscanlen = JAS_MIN(i, 4);
for (j = width; j > 0; --j, ++fvscanstart, ++dvscanstart) {
fp = fvscanstart;
dp = dvscanstart;
k = vscanlen;
/* Process first sample in vertical scan. */
jpc_rawrefpass_step(fp, dp, poshalf, neghalf, in,
vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process second sample in vertical scan. */
jpc_rawrefpass_step(fp, dp, poshalf, neghalf, in, 0);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process third sample in vertical scan. */
jpc_rawrefpass_step(fp, dp, poshalf, neghalf, in, 0);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process fourth sample in vertical scan. */
jpc_rawrefpass_step(fp, dp, poshalf, neghalf, in, 0);
}
}
return 0;
}
static gboolean
process (GeglOperation *operation,
GeglBuffer *output,
const GeglRectangle *result)
{
GeglChantO *o = GEGL_CHANT_PROPERTIES (operation);
GeglRectangle rect = {0,0,0,0};
jas_image_t *image;
gint width, height, depth;
gsize bpc;
guchar *data = NULL;
gboolean ret;
int components[3];
jas_matrix_t *matrices[3] = {NULL, NULL, NULL};
gint i;
gint row;
gboolean b;
gushort *ptr_s;
guchar *ptr_b;
image = NULL;
width = height = depth = 0;
if (!query_jp2 (o->path, &width, &height, &depth, &image))
return FALSE;
rect.height = height;
rect.width = width;
switch (depth)
{
case 8:
bpc = sizeof (guchar);
break;
case 16:
bpc = sizeof (gushort);
break;
default:
g_warning ("%s: Programmer stupidity error", G_STRLOC);
return FALSE;
}
data = (guchar *) g_malloc (width * height * 3 * bpc);
ptr_s = (gushort *) data;
ptr_b = data;
switch (depth)
{
case 16:
gegl_buffer_get (output, 1.0, &rect, babl_format ("R'G'B' u16"), data,
GEGL_AUTO_ROWSTRIDE);
break;
case 8:
default:
gegl_buffer_get (output, 1.0, &rect, babl_format ("R'G'B' u8"), data,
GEGL_AUTO_ROWSTRIDE);
}
ret = FALSE;
b = FALSE;
do
{
components[0] = jas_image_getcmptbytype
(image, JAS_IMAGE_CT_COLOR(JAS_IMAGE_CT_RGB_R));
components[1] = jas_image_getcmptbytype
(image, JAS_IMAGE_CT_COLOR(JAS_IMAGE_CT_RGB_G));
components[2] = jas_image_getcmptbytype
(image, JAS_IMAGE_CT_COLOR(JAS_IMAGE_CT_RGB_B));
if ((components[0] < 0) || (components[1] < 0) || (components[2] < 0))
{
g_warning (_("One or more of R, G, B components are missing "
"from '%s'"), o->path);
break;
}
if (jas_image_cmptsgnd (image, components[0]) ||
jas_image_cmptsgnd (image, components[1]) ||
jas_image_cmptsgnd (image, components[2]))
{
g_warning (_("One or more of R, G, B components have signed "
"data in '%s'"), o->path);
break;
}
for (i = 0; i < 3; i++)
matrices[i] = jas_matrix_create(1, width);
for (row = 0; row < height; row++)
{
gint plane, col;
jas_seqent_t *jrow[3] = {NULL, NULL, NULL};
for (plane = 0; plane < 3; plane++)
{
int r = jas_image_readcmpt (image, components[plane], 0, row,
width, 1, matrices[plane]);
if (r)
{
g_warning (_("Error reading row %d component %d from '%s'"),
row, plane, o->path);
b = TRUE;
break;
}
}
if (b)
break;
for (plane = 0; plane < 3; plane++)
jrow[plane] = jas_matrix_getref (matrices[plane], 0, 0);
for (col = 0; col < width; col++)
{
switch (depth)
{
case 16:
*ptr_s++ = (gushort) jrow[0][col];
*ptr_s++ = (gushort) jrow[1][col];
*ptr_s++ = (gushort) jrow[2][col];
break;
case 8:
default:
*ptr_b++ = (guchar) jrow[0][col];
*ptr_b++ = (guchar) jrow[1][col];
*ptr_b++ = (guchar) jrow[2][col];
}
}
}
if (b)
break;
switch (depth)
{
case 16:
gegl_buffer_set (output, &rect, babl_format ("R'G'B' u16"), data,
GEGL_AUTO_ROWSTRIDE);
break;
case 8:
default:
gegl_buffer_set (output, &rect, babl_format ("R'G'B' u8"), data,
GEGL_AUTO_ROWSTRIDE);
}
ret = TRUE;
}
while (FALSE); /* structured goto */
for (i = 0; i < 3; i++)
if (matrices[i])
jas_matrix_destroy (matrices[i]);
if (data)
g_free (data);
if (image)
jas_image_destroy (image);
return ret;
}
static int dec_rawsigpass(jpc_dec_t *dec, jpc_bitstream_t *in, int bitpos, int vcausalflag,
jas_matrix_t *flags, jas_matrix_t *data)
{
int i;
int j;
int k;
int one;
int half;
int oneplushalf;
int vscanlen;
int width;
int height;
jpc_fix_t *fp;
int frowstep;
int fstripestep;
jpc_fix_t *fstripestart;
jpc_fix_t *fvscanstart;
jpc_fix_t *dp;
int drowstep;
int dstripestep;
jpc_fix_t *dstripestart;
jpc_fix_t *dvscanstart;
width = jas_matrix_numcols(data);
height = jas_matrix_numrows(data);
frowstep = jas_matrix_rowstep(flags);
drowstep = jas_matrix_rowstep(data);
fstripestep = frowstep << 2;
dstripestep = drowstep << 2;
one = 1 << bitpos;
half = one >> 1;
oneplushalf = one | half;
fstripestart = jas_matrix_getref(flags, 1, 1);
dstripestart = jas_matrix_getref(data, 0, 0);
for (i = height; i > 0; i -= 4, fstripestart += fstripestep,
dstripestart += dstripestep) {
fvscanstart = fstripestart;
dvscanstart = dstripestart;
vscanlen = JAS_MIN(i, 4);
for (j = width; j > 0; --j, ++fvscanstart, ++dvscanstart) {
fp = fvscanstart;
dp = dvscanstart;
k = vscanlen;
/* Process first sample in vertical scan. */
jpc_rawsigpass_step(fp, frowstep, dp, oneplushalf,
in, vcausalflag);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process second sample in vertical scan. */
jpc_rawsigpass_step(fp, frowstep, dp, oneplushalf,
in, 0);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process third sample in vertical scan. */
jpc_rawsigpass_step(fp, frowstep, dp, oneplushalf,
in, 0);
if (--k <= 0) {
continue;
}
fp += frowstep;
dp += drowstep;
/* Process fourth sample in vertical scan. */
jpc_rawsigpass_step(fp, frowstep, dp, oneplushalf,
in, 0);
}
}
return 0;
}