| /* $Id: tif_pixarlog.c,v 1.53 2017-05-17 09:53:06 erouault Exp $ */ |
| |
| /* |
| * Copyright (c) 1996-1997 Sam Leffler |
| * Copyright (c) 1996 Pixar |
| * |
| * Permission to use, copy, modify, distribute, and sell this software and |
| * its documentation for any purpose is hereby granted without fee, provided |
| * that (i) the above copyright notices and this permission notice appear in |
| * all copies of the software and related documentation, and (ii) the names of |
| * Pixar, Sam Leffler and Silicon Graphics may not be used in any advertising or |
| * publicity relating to the software without the specific, prior written |
| * permission of Pixar, Sam Leffler and Silicon Graphics. |
| * |
| * THE SOFTWARE IS PROVIDED "AS-IS" AND WITHOUT WARRANTY OF ANY KIND, |
| * EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY |
| * WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. |
| * |
| * IN NO EVENT SHALL PIXAR, SAM LEFFLER OR SILICON GRAPHICS BE LIABLE FOR |
| * ANY SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY KIND, |
| * OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, |
| * WHETHER OR NOT ADVISED OF THE POSSIBILITY OF DAMAGE, AND ON ANY THEORY OF |
| * LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE |
| * OF THIS SOFTWARE. |
| */ |
| |
| #include "tiffiop.h" |
| #ifdef PIXARLOG_SUPPORT |
| |
| /* |
| * TIFF Library. |
| * PixarLog Compression Support |
| * |
| * Contributed by Dan McCoy. |
| * |
| * PixarLog film support uses the TIFF library to store companded |
| * 11 bit values into a tiff file, which are compressed using the |
| * zip compressor. |
| * |
| * The codec can take as input and produce as output 32-bit IEEE float values |
| * as well as 16-bit or 8-bit unsigned integer values. |
| * |
| * On writing any of the above are converted into the internal |
| * 11-bit log format. In the case of 8 and 16 bit values, the |
| * input is assumed to be unsigned linear color values that represent |
| * the range 0-1. In the case of IEEE values, the 0-1 range is assumed to |
| * be the normal linear color range, in addition over 1 values are |
| * accepted up to a value of about 25.0 to encode "hot" highlights and such. |
| * The encoding is lossless for 8-bit values, slightly lossy for the |
| * other bit depths. The actual color precision should be better |
| * than the human eye can perceive with extra room to allow for |
| * error introduced by further image computation. As with any quantized |
| * color format, it is possible to perform image calculations which |
| * expose the quantization error. This format should certainly be less |
| * susceptible to such errors than standard 8-bit encodings, but more |
| * susceptible than straight 16-bit or 32-bit encodings. |
| * |
| * On reading the internal format is converted to the desired output format. |
| * The program can request which format it desires by setting the internal |
| * pseudo tag TIFFTAG_PIXARLOGDATAFMT to one of these possible values: |
| * PIXARLOGDATAFMT_FLOAT = provide IEEE float values. |
| * PIXARLOGDATAFMT_16BIT = provide unsigned 16-bit integer values |
| * PIXARLOGDATAFMT_8BIT = provide unsigned 8-bit integer values |
| * |
| * alternately PIXARLOGDATAFMT_8BITABGR provides unsigned 8-bit integer |
| * values with the difference that if there are exactly three or four channels |
| * (rgb or rgba) it swaps the channel order (bgr or abgr). |
| * |
| * PIXARLOGDATAFMT_11BITLOG provides the internal encoding directly |
| * packed in 16-bit values. However no tools are supplied for interpreting |
| * these values. |
| * |
| * "hot" (over 1.0) areas written in floating point get clamped to |
| * 1.0 in the integer data types. |
| * |
| * When the file is closed after writing, the bit depth and sample format |
| * are set always to appear as if 8-bit data has been written into it. |
| * That way a naive program unaware of the particulars of the encoding |
| * gets the format it is most likely able to handle. |
| * |
| * The codec does it's own horizontal differencing step on the coded |
| * values so the libraries predictor stuff should be turned off. |
| * The codec also handle byte swapping the encoded values as necessary |
| * since the library does not have the information necessary |
| * to know the bit depth of the raw unencoded buffer. |
| * |
| * NOTE: This decoder does not appear to update tif_rawcp, and tif_rawcc. |
| * This can cause problems with the implementation of CHUNKY_STRIP_READ_SUPPORT |
| * as noted in http://trac.osgeo.org/gdal/ticket/3894. FrankW - Jan'11 |
| */ |
| |
| #include "tif_predict.h" |
| #include "zlib.h" |
| |
| #include <stdio.h> |
| #include <stdlib.h> |
| #include <math.h> |
| |
| /* Tables for converting to/from 11 bit coded values */ |
| |
| #define TSIZE 2048 /* decode table size (11-bit tokens) */ |
| #define TSIZEP1 2049 /* Plus one for slop */ |
| #define ONE 1250 /* token value of 1.0 exactly */ |
| #define RATIO 1.004 /* nominal ratio for log part */ |
| |
| #define CODE_MASK 0x7ff /* 11 bits. */ |
| |
| static float Fltsize; |
| static float LogK1, LogK2; |
| |
| #define REPEAT(n, op) { int i; i=n; do { i--; op; } while (i>0); } |
| |
| static void |
| horizontalAccumulateF(uint16 *wp, int n, int stride, float *op, |
| float *ToLinearF) |
| { |
| register unsigned int cr, cg, cb, ca, mask; |
| register float t0, t1, t2, t3; |
| |
| if (n >= stride) { |
| mask = CODE_MASK; |
| if (stride == 3) { |
| t0 = ToLinearF[cr = (wp[0] & mask)]; |
| t1 = ToLinearF[cg = (wp[1] & mask)]; |
| t2 = ToLinearF[cb = (wp[2] & mask)]; |
| op[0] = t0; |
| op[1] = t1; |
| op[2] = t2; |
| n -= 3; |
| while (n > 0) { |
| wp += 3; |
| op += 3; |
| n -= 3; |
| t0 = ToLinearF[(cr += wp[0]) & mask]; |
| t1 = ToLinearF[(cg += wp[1]) & mask]; |
| t2 = ToLinearF[(cb += wp[2]) & mask]; |
| op[0] = t0; |
| op[1] = t1; |
| op[2] = t2; |
| } |
| } else if (stride == 4) { |
| t0 = ToLinearF[cr = (wp[0] & mask)]; |
| t1 = ToLinearF[cg = (wp[1] & mask)]; |
| t2 = ToLinearF[cb = (wp[2] & mask)]; |
| t3 = ToLinearF[ca = (wp[3] & mask)]; |
| op[0] = t0; |
| op[1] = t1; |
| op[2] = t2; |
| op[3] = t3; |
| n -= 4; |
| while (n > 0) { |
| wp += 4; |
| op += 4; |
| n -= 4; |
| t0 = ToLinearF[(cr += wp[0]) & mask]; |
| t1 = ToLinearF[(cg += wp[1]) & mask]; |
| t2 = ToLinearF[(cb += wp[2]) & mask]; |
| t3 = ToLinearF[(ca += wp[3]) & mask]; |
| op[0] = t0; |
| op[1] = t1; |
| op[2] = t2; |
| op[3] = t3; |
| } |
| } else { |
| REPEAT(stride, *op = ToLinearF[*wp&mask]; wp++; op++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[stride] += *wp; *op = ToLinearF[*wp&mask]; wp++; op++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| static void |
| horizontalAccumulate12(uint16 *wp, int n, int stride, int16 *op, |
| float *ToLinearF) |
| { |
| register unsigned int cr, cg, cb, ca, mask; |
| register float t0, t1, t2, t3; |
| |
| #define SCALE12 2048.0F |
| #define CLAMP12(t) (((t) < 3071) ? (uint16) (t) : 3071) |
| |
| if (n >= stride) { |
| mask = CODE_MASK; |
| if (stride == 3) { |
| t0 = ToLinearF[cr = (wp[0] & mask)] * SCALE12; |
| t1 = ToLinearF[cg = (wp[1] & mask)] * SCALE12; |
| t2 = ToLinearF[cb = (wp[2] & mask)] * SCALE12; |
| op[0] = CLAMP12(t0); |
| op[1] = CLAMP12(t1); |
| op[2] = CLAMP12(t2); |
| n -= 3; |
| while (n > 0) { |
| wp += 3; |
| op += 3; |
| n -= 3; |
| t0 = ToLinearF[(cr += wp[0]) & mask] * SCALE12; |
| t1 = ToLinearF[(cg += wp[1]) & mask] * SCALE12; |
| t2 = ToLinearF[(cb += wp[2]) & mask] * SCALE12; |
| op[0] = CLAMP12(t0); |
| op[1] = CLAMP12(t1); |
| op[2] = CLAMP12(t2); |
| } |
| } else if (stride == 4) { |
| t0 = ToLinearF[cr = (wp[0] & mask)] * SCALE12; |
| t1 = ToLinearF[cg = (wp[1] & mask)] * SCALE12; |
| t2 = ToLinearF[cb = (wp[2] & mask)] * SCALE12; |
| t3 = ToLinearF[ca = (wp[3] & mask)] * SCALE12; |
| op[0] = CLAMP12(t0); |
| op[1] = CLAMP12(t1); |
| op[2] = CLAMP12(t2); |
| op[3] = CLAMP12(t3); |
| n -= 4; |
| while (n > 0) { |
| wp += 4; |
| op += 4; |
| n -= 4; |
| t0 = ToLinearF[(cr += wp[0]) & mask] * SCALE12; |
| t1 = ToLinearF[(cg += wp[1]) & mask] * SCALE12; |
| t2 = ToLinearF[(cb += wp[2]) & mask] * SCALE12; |
| t3 = ToLinearF[(ca += wp[3]) & mask] * SCALE12; |
| op[0] = CLAMP12(t0); |
| op[1] = CLAMP12(t1); |
| op[2] = CLAMP12(t2); |
| op[3] = CLAMP12(t3); |
| } |
| } else { |
| REPEAT(stride, t0 = ToLinearF[*wp&mask] * SCALE12; |
| *op = CLAMP12(t0); wp++; op++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[stride] += *wp; t0 = ToLinearF[wp[stride]&mask]*SCALE12; |
| *op = CLAMP12(t0); wp++; op++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| static void |
| horizontalAccumulate16(uint16 *wp, int n, int stride, uint16 *op, |
| uint16 *ToLinear16) |
| { |
| register unsigned int cr, cg, cb, ca, mask; |
| |
| if (n >= stride) { |
| mask = CODE_MASK; |
| if (stride == 3) { |
| op[0] = ToLinear16[cr = (wp[0] & mask)]; |
| op[1] = ToLinear16[cg = (wp[1] & mask)]; |
| op[2] = ToLinear16[cb = (wp[2] & mask)]; |
| n -= 3; |
| while (n > 0) { |
| wp += 3; |
| op += 3; |
| n -= 3; |
| op[0] = ToLinear16[(cr += wp[0]) & mask]; |
| op[1] = ToLinear16[(cg += wp[1]) & mask]; |
| op[2] = ToLinear16[(cb += wp[2]) & mask]; |
| } |
| } else if (stride == 4) { |
| op[0] = ToLinear16[cr = (wp[0] & mask)]; |
| op[1] = ToLinear16[cg = (wp[1] & mask)]; |
| op[2] = ToLinear16[cb = (wp[2] & mask)]; |
| op[3] = ToLinear16[ca = (wp[3] & mask)]; |
| n -= 4; |
| while (n > 0) { |
| wp += 4; |
| op += 4; |
| n -= 4; |
| op[0] = ToLinear16[(cr += wp[0]) & mask]; |
| op[1] = ToLinear16[(cg += wp[1]) & mask]; |
| op[2] = ToLinear16[(cb += wp[2]) & mask]; |
| op[3] = ToLinear16[(ca += wp[3]) & mask]; |
| } |
| } else { |
| REPEAT(stride, *op = ToLinear16[*wp&mask]; wp++; op++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[stride] += *wp; *op = ToLinear16[*wp&mask]; wp++; op++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| /* |
| * Returns the log encoded 11-bit values with the horizontal |
| * differencing undone. |
| */ |
| static void |
| horizontalAccumulate11(uint16 *wp, int n, int stride, uint16 *op) |
| { |
| register unsigned int cr, cg, cb, ca, mask; |
| |
| if (n >= stride) { |
| mask = CODE_MASK; |
| if (stride == 3) { |
| op[0] = wp[0]; op[1] = wp[1]; op[2] = wp[2]; |
| cr = wp[0]; cg = wp[1]; cb = wp[2]; |
| n -= 3; |
| while (n > 0) { |
| wp += 3; |
| op += 3; |
| n -= 3; |
| op[0] = (uint16)((cr += wp[0]) & mask); |
| op[1] = (uint16)((cg += wp[1]) & mask); |
| op[2] = (uint16)((cb += wp[2]) & mask); |
| } |
| } else if (stride == 4) { |
| op[0] = wp[0]; op[1] = wp[1]; |
| op[2] = wp[2]; op[3] = wp[3]; |
| cr = wp[0]; cg = wp[1]; cb = wp[2]; ca = wp[3]; |
| n -= 4; |
| while (n > 0) { |
| wp += 4; |
| op += 4; |
| n -= 4; |
| op[0] = (uint16)((cr += wp[0]) & mask); |
| op[1] = (uint16)((cg += wp[1]) & mask); |
| op[2] = (uint16)((cb += wp[2]) & mask); |
| op[3] = (uint16)((ca += wp[3]) & mask); |
| } |
| } else { |
| REPEAT(stride, *op = *wp&mask; wp++; op++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[stride] += *wp; *op = *wp&mask; wp++; op++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| static void |
| horizontalAccumulate8(uint16 *wp, int n, int stride, unsigned char *op, |
| unsigned char *ToLinear8) |
| { |
| register unsigned int cr, cg, cb, ca, mask; |
| |
| if (n >= stride) { |
| mask = CODE_MASK; |
| if (stride == 3) { |
| op[0] = ToLinear8[cr = (wp[0] & mask)]; |
| op[1] = ToLinear8[cg = (wp[1] & mask)]; |
| op[2] = ToLinear8[cb = (wp[2] & mask)]; |
| n -= 3; |
| while (n > 0) { |
| n -= 3; |
| wp += 3; |
| op += 3; |
| op[0] = ToLinear8[(cr += wp[0]) & mask]; |
| op[1] = ToLinear8[(cg += wp[1]) & mask]; |
| op[2] = ToLinear8[(cb += wp[2]) & mask]; |
| } |
| } else if (stride == 4) { |
| op[0] = ToLinear8[cr = (wp[0] & mask)]; |
| op[1] = ToLinear8[cg = (wp[1] & mask)]; |
| op[2] = ToLinear8[cb = (wp[2] & mask)]; |
| op[3] = ToLinear8[ca = (wp[3] & mask)]; |
| n -= 4; |
| while (n > 0) { |
| n -= 4; |
| wp += 4; |
| op += 4; |
| op[0] = ToLinear8[(cr += wp[0]) & mask]; |
| op[1] = ToLinear8[(cg += wp[1]) & mask]; |
| op[2] = ToLinear8[(cb += wp[2]) & mask]; |
| op[3] = ToLinear8[(ca += wp[3]) & mask]; |
| } |
| } else { |
| REPEAT(stride, *op = ToLinear8[*wp&mask]; wp++; op++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[stride] += *wp; *op = ToLinear8[*wp&mask]; wp++; op++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| |
| static void |
| horizontalAccumulate8abgr(uint16 *wp, int n, int stride, unsigned char *op, |
| unsigned char *ToLinear8) |
| { |
| register unsigned int cr, cg, cb, ca, mask; |
| register unsigned char t0, t1, t2, t3; |
| |
| if (n >= stride) { |
| mask = CODE_MASK; |
| if (stride == 3) { |
| op[0] = 0; |
| t1 = ToLinear8[cb = (wp[2] & mask)]; |
| t2 = ToLinear8[cg = (wp[1] & mask)]; |
| t3 = ToLinear8[cr = (wp[0] & mask)]; |
| op[1] = t1; |
| op[2] = t2; |
| op[3] = t3; |
| n -= 3; |
| while (n > 0) { |
| n -= 3; |
| wp += 3; |
| op += 4; |
| op[0] = 0; |
| t1 = ToLinear8[(cb += wp[2]) & mask]; |
| t2 = ToLinear8[(cg += wp[1]) & mask]; |
| t3 = ToLinear8[(cr += wp[0]) & mask]; |
| op[1] = t1; |
| op[2] = t2; |
| op[3] = t3; |
| } |
| } else if (stride == 4) { |
| t0 = ToLinear8[ca = (wp[3] & mask)]; |
| t1 = ToLinear8[cb = (wp[2] & mask)]; |
| t2 = ToLinear8[cg = (wp[1] & mask)]; |
| t3 = ToLinear8[cr = (wp[0] & mask)]; |
| op[0] = t0; |
| op[1] = t1; |
| op[2] = t2; |
| op[3] = t3; |
| n -= 4; |
| while (n > 0) { |
| n -= 4; |
| wp += 4; |
| op += 4; |
| t0 = ToLinear8[(ca += wp[3]) & mask]; |
| t1 = ToLinear8[(cb += wp[2]) & mask]; |
| t2 = ToLinear8[(cg += wp[1]) & mask]; |
| t3 = ToLinear8[(cr += wp[0]) & mask]; |
| op[0] = t0; |
| op[1] = t1; |
| op[2] = t2; |
| op[3] = t3; |
| } |
| } else { |
| REPEAT(stride, *op = ToLinear8[*wp&mask]; wp++; op++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[stride] += *wp; *op = ToLinear8[*wp&mask]; wp++; op++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| /* |
| * State block for each open TIFF |
| * file using PixarLog compression/decompression. |
| */ |
| typedef struct { |
| TIFFPredictorState predict; |
| z_stream stream; |
| tmsize_t tbuf_size; /* only set/used on reading for now */ |
| uint16 *tbuf; |
| uint16 stride; |
| int state; |
| int user_datafmt; |
| int quality; |
| #define PLSTATE_INIT 1 |
| |
| TIFFVSetMethod vgetparent; /* super-class method */ |
| TIFFVSetMethod vsetparent; /* super-class method */ |
| |
| float *ToLinearF; |
| uint16 *ToLinear16; |
| unsigned char *ToLinear8; |
| uint16 *FromLT2; |
| uint16 *From14; /* Really for 16-bit data, but we shift down 2 */ |
| uint16 *From8; |
| |
| } PixarLogState; |
| |
| static int |
| PixarLogMakeTables(PixarLogState *sp) |
| { |
| |
| /* |
| * We make several tables here to convert between various external |
| * representations (float, 16-bit, and 8-bit) and the internal |
| * 11-bit companded representation. The 11-bit representation has two |
| * distinct regions. A linear bottom end up through .018316 in steps |
| * of about .000073, and a region of constant ratio up to about 25. |
| * These floating point numbers are stored in the main table ToLinearF. |
| * All other tables are derived from this one. The tables (and the |
| * ratios) are continuous at the internal seam. |
| */ |
| |
| int nlin, lt2size; |
| int i, j; |
| double b, c, linstep, v; |
| float *ToLinearF; |
| uint16 *ToLinear16; |
| unsigned char *ToLinear8; |
| uint16 *FromLT2; |
| uint16 *From14; /* Really for 16-bit data, but we shift down 2 */ |
| uint16 *From8; |
| |
| c = log(RATIO); |
| nlin = (int)(1./c); /* nlin must be an integer */ |
| c = 1./nlin; |
| b = exp(-c*ONE); /* multiplicative scale factor [b*exp(c*ONE) = 1] */ |
| linstep = b*c*exp(1.); |
| |
| LogK1 = (float)(1./c); /* if (v >= 2) token = k1*log(v*k2) */ |
| LogK2 = (float)(1./b); |
| lt2size = (int)(2./linstep) + 1; |
| FromLT2 = (uint16 *)_TIFFmalloc(lt2size*sizeof(uint16)); |
| From14 = (uint16 *)_TIFFmalloc(16384*sizeof(uint16)); |
| From8 = (uint16 *)_TIFFmalloc(256*sizeof(uint16)); |
| ToLinearF = (float *)_TIFFmalloc(TSIZEP1 * sizeof(float)); |
| ToLinear16 = (uint16 *)_TIFFmalloc(TSIZEP1 * sizeof(uint16)); |
| ToLinear8 = (unsigned char *)_TIFFmalloc(TSIZEP1 * sizeof(unsigned char)); |
| if (FromLT2 == NULL || From14 == NULL || From8 == NULL || |
| ToLinearF == NULL || ToLinear16 == NULL || ToLinear8 == NULL) { |
| if (FromLT2) _TIFFfree(FromLT2); |
| if (From14) _TIFFfree(From14); |
| if (From8) _TIFFfree(From8); |
| if (ToLinearF) _TIFFfree(ToLinearF); |
| if (ToLinear16) _TIFFfree(ToLinear16); |
| if (ToLinear8) _TIFFfree(ToLinear8); |
| sp->FromLT2 = NULL; |
| sp->From14 = NULL; |
| sp->From8 = NULL; |
| sp->ToLinearF = NULL; |
| sp->ToLinear16 = NULL; |
| sp->ToLinear8 = NULL; |
| return 0; |
| } |
| |
| j = 0; |
| |
| for (i = 0; i < nlin; i++) { |
| v = i * linstep; |
| ToLinearF[j++] = (float)v; |
| } |
| |
| for (i = nlin; i < TSIZE; i++) |
| ToLinearF[j++] = (float)(b*exp(c*i)); |
| |
| ToLinearF[2048] = ToLinearF[2047]; |
| |
| for (i = 0; i < TSIZEP1; i++) { |
| v = ToLinearF[i]*65535.0 + 0.5; |
| ToLinear16[i] = (v > 65535.0) ? 65535 : (uint16)v; |
| v = ToLinearF[i]*255.0 + 0.5; |
| ToLinear8[i] = (v > 255.0) ? 255 : (unsigned char)v; |
| } |
| |
| j = 0; |
| for (i = 0; i < lt2size; i++) { |
| if ((i*linstep)*(i*linstep) > ToLinearF[j]*ToLinearF[j+1]) |
| j++; |
| FromLT2[i] = (uint16)j; |
| } |
| |
| /* |
| * Since we lose info anyway on 16-bit data, we set up a 14-bit |
| * table and shift 16-bit values down two bits on input. |
| * saves a little table space. |
| */ |
| j = 0; |
| for (i = 0; i < 16384; i++) { |
| while ((i/16383.)*(i/16383.) > ToLinearF[j]*ToLinearF[j+1]) |
| j++; |
| From14[i] = (uint16)j; |
| } |
| |
| j = 0; |
| for (i = 0; i < 256; i++) { |
| while ((i/255.)*(i/255.) > ToLinearF[j]*ToLinearF[j+1]) |
| j++; |
| From8[i] = (uint16)j; |
| } |
| |
| Fltsize = (float)(lt2size/2); |
| |
| sp->ToLinearF = ToLinearF; |
| sp->ToLinear16 = ToLinear16; |
| sp->ToLinear8 = ToLinear8; |
| sp->FromLT2 = FromLT2; |
| sp->From14 = From14; |
| sp->From8 = From8; |
| |
| return 1; |
| } |
| |
| #define DecoderState(tif) ((PixarLogState*) (tif)->tif_data) |
| #define EncoderState(tif) ((PixarLogState*) (tif)->tif_data) |
| |
| static int PixarLogEncode(TIFF* tif, uint8* bp, tmsize_t cc, uint16 s); |
| static int PixarLogDecode(TIFF* tif, uint8* op, tmsize_t occ, uint16 s); |
| |
| #define PIXARLOGDATAFMT_UNKNOWN -1 |
| |
| static int |
| PixarLogGuessDataFmt(TIFFDirectory *td) |
| { |
| int guess = PIXARLOGDATAFMT_UNKNOWN; |
| int format = td->td_sampleformat; |
| |
| /* If the user didn't tell us his datafmt, |
| * take our best guess from the bitspersample. |
| */ |
| switch (td->td_bitspersample) { |
| case 32: |
| if (format == SAMPLEFORMAT_IEEEFP) |
| guess = PIXARLOGDATAFMT_FLOAT; |
| break; |
| case 16: |
| if (format == SAMPLEFORMAT_VOID || format == SAMPLEFORMAT_UINT) |
| guess = PIXARLOGDATAFMT_16BIT; |
| break; |
| case 12: |
| if (format == SAMPLEFORMAT_VOID || format == SAMPLEFORMAT_INT) |
| guess = PIXARLOGDATAFMT_12BITPICIO; |
| break; |
| case 11: |
| if (format == SAMPLEFORMAT_VOID || format == SAMPLEFORMAT_UINT) |
| guess = PIXARLOGDATAFMT_11BITLOG; |
| break; |
| case 8: |
| if (format == SAMPLEFORMAT_VOID || format == SAMPLEFORMAT_UINT) |
| guess = PIXARLOGDATAFMT_8BIT; |
| break; |
| } |
| |
| return guess; |
| } |
| |
| #define TIFF_SIZE_T_MAX ((size_t) ~ ((size_t)0)) |
| #define TIFF_TMSIZE_T_MAX (tmsize_t)(TIFF_SIZE_T_MAX >> 1) |
| |
| static tmsize_t |
| multiply_ms(tmsize_t m1, tmsize_t m2) |
| { |
| if( m1 == 0 || m2 > TIFF_TMSIZE_T_MAX / m1 ) |
| return 0; |
| return m1 * m2; |
| } |
| |
| static tmsize_t |
| add_ms(tmsize_t m1, tmsize_t m2) |
| { |
| /* if either input is zero, assume overflow already occurred */ |
| if (m1 == 0 || m2 == 0) |
| return 0; |
| else if (m1 > TIFF_TMSIZE_T_MAX - m2) |
| return 0; |
| |
| return m1 + m2; |
| } |
| |
| static int |
| PixarLogFixupTags(TIFF* tif) |
| { |
| (void) tif; |
| return (1); |
| } |
| |
| static int |
| PixarLogSetupDecode(TIFF* tif) |
| { |
| static const char module[] = "PixarLogSetupDecode"; |
| TIFFDirectory *td = &tif->tif_dir; |
| PixarLogState* sp = DecoderState(tif); |
| tmsize_t tbuf_size; |
| |
| assert(sp != NULL); |
| |
| /* This function can possibly be called several times by */ |
| /* PredictorSetupDecode() if this function succeeds but */ |
| /* PredictorSetup() fails */ |
| if( (sp->state & PLSTATE_INIT) != 0 ) |
| return 1; |
| |
| /* Make sure no byte swapping happens on the data |
| * after decompression. */ |
| tif->tif_postdecode = _TIFFNoPostDecode; |
| |
| /* for some reason, we can't do this in TIFFInitPixarLog */ |
| |
| sp->stride = (td->td_planarconfig == PLANARCONFIG_CONTIG ? |
| td->td_samplesperpixel : 1); |
| tbuf_size = multiply_ms(multiply_ms(multiply_ms(sp->stride, td->td_imagewidth), |
| td->td_rowsperstrip), sizeof(uint16)); |
| /* add one more stride in case input ends mid-stride */ |
| tbuf_size = add_ms(tbuf_size, sizeof(uint16) * sp->stride); |
| if (tbuf_size == 0) |
| return (0); /* TODO: this is an error return without error report through TIFFErrorExt */ |
| sp->tbuf = (uint16 *) _TIFFmalloc(tbuf_size); |
| if (sp->tbuf == NULL) |
| return (0); |
| sp->tbuf_size = tbuf_size; |
| if (sp->user_datafmt == PIXARLOGDATAFMT_UNKNOWN) |
| sp->user_datafmt = PixarLogGuessDataFmt(td); |
| if (sp->user_datafmt == PIXARLOGDATAFMT_UNKNOWN) { |
| _TIFFfree(sp->tbuf); |
| sp->tbuf = NULL; |
| sp->tbuf_size = 0; |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "PixarLog compression can't handle bits depth/data format combination (depth: %d)", |
| td->td_bitspersample); |
| return (0); |
| } |
| |
| if (inflateInit(&sp->stream) != Z_OK) { |
| _TIFFfree(sp->tbuf); |
| sp->tbuf = NULL; |
| sp->tbuf_size = 0; |
| TIFFErrorExt(tif->tif_clientdata, module, "%s", sp->stream.msg ? sp->stream.msg : "(null)"); |
| return (0); |
| } else { |
| sp->state |= PLSTATE_INIT; |
| return (1); |
| } |
| } |
| |
| /* |
| * Setup state for decoding a strip. |
| */ |
| static int |
| PixarLogPreDecode(TIFF* tif, uint16 s) |
| { |
| static const char module[] = "PixarLogPreDecode"; |
| PixarLogState* sp = DecoderState(tif); |
| |
| (void) s; |
| assert(sp != NULL); |
| sp->stream.next_in = tif->tif_rawdata; |
| assert(sizeof(sp->stream.avail_in)==4); /* if this assert gets raised, |
| we need to simplify this code to reflect a ZLib that is likely updated |
| to deal with 8byte memory sizes, though this code will respond |
| appropriately even before we simplify it */ |
| sp->stream.avail_in = (uInt) tif->tif_rawcc; |
| if ((tmsize_t)sp->stream.avail_in != tif->tif_rawcc) |
| { |
| TIFFErrorExt(tif->tif_clientdata, module, "ZLib cannot deal with buffers this size"); |
| return (0); |
| } |
| return (inflateReset(&sp->stream) == Z_OK); |
| } |
| |
| static int |
| PixarLogDecode(TIFF* tif, uint8* op, tmsize_t occ, uint16 s) |
| { |
| static const char module[] = "PixarLogDecode"; |
| TIFFDirectory *td = &tif->tif_dir; |
| PixarLogState* sp = DecoderState(tif); |
| tmsize_t i; |
| tmsize_t nsamples; |
| int llen; |
| uint16 *up; |
| |
| switch (sp->user_datafmt) { |
| case PIXARLOGDATAFMT_FLOAT: |
| nsamples = occ / sizeof(float); /* XXX float == 32 bits */ |
| break; |
| case PIXARLOGDATAFMT_16BIT: |
| case PIXARLOGDATAFMT_12BITPICIO: |
| case PIXARLOGDATAFMT_11BITLOG: |
| nsamples = occ / sizeof(uint16); /* XXX uint16 == 16 bits */ |
| break; |
| case PIXARLOGDATAFMT_8BIT: |
| case PIXARLOGDATAFMT_8BITABGR: |
| nsamples = occ; |
| break; |
| default: |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "%d bit input not supported in PixarLog", |
| td->td_bitspersample); |
| return 0; |
| } |
| |
| llen = sp->stride * td->td_imagewidth; |
| |
| (void) s; |
| assert(sp != NULL); |
| |
| sp->stream.next_in = tif->tif_rawcp; |
| sp->stream.avail_in = (uInt) tif->tif_rawcc; |
| |
| sp->stream.next_out = (unsigned char *) sp->tbuf; |
| assert(sizeof(sp->stream.avail_out)==4); /* if this assert gets raised, |
| we need to simplify this code to reflect a ZLib that is likely updated |
| to deal with 8byte memory sizes, though this code will respond |
| appropriately even before we simplify it */ |
| sp->stream.avail_out = (uInt) (nsamples * sizeof(uint16)); |
| if (sp->stream.avail_out != nsamples * sizeof(uint16)) |
| { |
| TIFFErrorExt(tif->tif_clientdata, module, "ZLib cannot deal with buffers this size"); |
| return (0); |
| } |
| /* Check that we will not fill more than what was allocated */ |
| if ((tmsize_t)sp->stream.avail_out > sp->tbuf_size) |
| { |
| TIFFErrorExt(tif->tif_clientdata, module, "sp->stream.avail_out > sp->tbuf_size"); |
| return (0); |
| } |
| do { |
| int state = inflate(&sp->stream, Z_PARTIAL_FLUSH); |
| if (state == Z_STREAM_END) { |
| break; /* XXX */ |
| } |
| if (state == Z_DATA_ERROR) { |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "Decoding error at scanline %lu, %s", |
| (unsigned long) tif->tif_row, sp->stream.msg ? sp->stream.msg : "(null)"); |
| if (inflateSync(&sp->stream) != Z_OK) |
| return (0); |
| continue; |
| } |
| if (state != Z_OK) { |
| TIFFErrorExt(tif->tif_clientdata, module, "ZLib error: %s", |
| sp->stream.msg ? sp->stream.msg : "(null)"); |
| return (0); |
| } |
| } while (sp->stream.avail_out > 0); |
| |
| /* hopefully, we got all the bytes we needed */ |
| if (sp->stream.avail_out != 0) { |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "Not enough data at scanline %lu (short " TIFF_UINT64_FORMAT " bytes)", |
| (unsigned long) tif->tif_row, (TIFF_UINT64_T) sp->stream.avail_out); |
| return (0); |
| } |
| |
| tif->tif_rawcp = sp->stream.next_in; |
| tif->tif_rawcc = sp->stream.avail_in; |
| |
| up = sp->tbuf; |
| /* Swap bytes in the data if from a different endian machine. */ |
| if (tif->tif_flags & TIFF_SWAB) |
| TIFFSwabArrayOfShort(up, nsamples); |
| |
| /* |
| * if llen is not an exact multiple of nsamples, the decode operation |
| * may overflow the output buffer, so truncate it enough to prevent |
| * that but still salvage as much data as possible. |
| */ |
| if (nsamples % llen) { |
| TIFFWarningExt(tif->tif_clientdata, module, |
| "stride %lu is not a multiple of sample count, " |
| "%lu, data truncated.", (unsigned long) llen, (unsigned long) nsamples); |
| nsamples -= nsamples % llen; |
| } |
| |
| for (i = 0; i < nsamples; i += llen, up += llen) { |
| switch (sp->user_datafmt) { |
| case PIXARLOGDATAFMT_FLOAT: |
| horizontalAccumulateF(up, llen, sp->stride, |
| (float *)op, sp->ToLinearF); |
| op += llen * sizeof(float); |
| break; |
| case PIXARLOGDATAFMT_16BIT: |
| horizontalAccumulate16(up, llen, sp->stride, |
| (uint16 *)op, sp->ToLinear16); |
| op += llen * sizeof(uint16); |
| break; |
| case PIXARLOGDATAFMT_12BITPICIO: |
| horizontalAccumulate12(up, llen, sp->stride, |
| (int16 *)op, sp->ToLinearF); |
| op += llen * sizeof(int16); |
| break; |
| case PIXARLOGDATAFMT_11BITLOG: |
| horizontalAccumulate11(up, llen, sp->stride, |
| (uint16 *)op); |
| op += llen * sizeof(uint16); |
| break; |
| case PIXARLOGDATAFMT_8BIT: |
| horizontalAccumulate8(up, llen, sp->stride, |
| (unsigned char *)op, sp->ToLinear8); |
| op += llen * sizeof(unsigned char); |
| break; |
| case PIXARLOGDATAFMT_8BITABGR: |
| horizontalAccumulate8abgr(up, llen, sp->stride, |
| (unsigned char *)op, sp->ToLinear8); |
| op += llen * sizeof(unsigned char); |
| break; |
| default: |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "Unsupported bits/sample: %d", |
| td->td_bitspersample); |
| return (0); |
| } |
| } |
| |
| return (1); |
| } |
| |
| static int |
| PixarLogSetupEncode(TIFF* tif) |
| { |
| static const char module[] = "PixarLogSetupEncode"; |
| TIFFDirectory *td = &tif->tif_dir; |
| PixarLogState* sp = EncoderState(tif); |
| tmsize_t tbuf_size; |
| |
| assert(sp != NULL); |
| |
| /* for some reason, we can't do this in TIFFInitPixarLog */ |
| |
| sp->stride = (td->td_planarconfig == PLANARCONFIG_CONTIG ? |
| td->td_samplesperpixel : 1); |
| tbuf_size = multiply_ms(multiply_ms(multiply_ms(sp->stride, td->td_imagewidth), |
| td->td_rowsperstrip), sizeof(uint16)); |
| if (tbuf_size == 0) |
| return (0); /* TODO: this is an error return without error report through TIFFErrorExt */ |
| sp->tbuf = (uint16 *) _TIFFmalloc(tbuf_size); |
| if (sp->tbuf == NULL) |
| return (0); |
| if (sp->user_datafmt == PIXARLOGDATAFMT_UNKNOWN) |
| sp->user_datafmt = PixarLogGuessDataFmt(td); |
| if (sp->user_datafmt == PIXARLOGDATAFMT_UNKNOWN) { |
| TIFFErrorExt(tif->tif_clientdata, module, "PixarLog compression can't handle %d bit linear encodings", td->td_bitspersample); |
| return (0); |
| } |
| |
| if (deflateInit(&sp->stream, sp->quality) != Z_OK) { |
| TIFFErrorExt(tif->tif_clientdata, module, "%s", sp->stream.msg ? sp->stream.msg : "(null)"); |
| return (0); |
| } else { |
| sp->state |= PLSTATE_INIT; |
| return (1); |
| } |
| } |
| |
| /* |
| * Reset encoding state at the start of a strip. |
| */ |
| static int |
| PixarLogPreEncode(TIFF* tif, uint16 s) |
| { |
| static const char module[] = "PixarLogPreEncode"; |
| PixarLogState *sp = EncoderState(tif); |
| |
| (void) s; |
| assert(sp != NULL); |
| sp->stream.next_out = tif->tif_rawdata; |
| assert(sizeof(sp->stream.avail_out)==4); /* if this assert gets raised, |
| we need to simplify this code to reflect a ZLib that is likely updated |
| to deal with 8byte memory sizes, though this code will respond |
| appropriately even before we simplify it */ |
| sp->stream.avail_out = (uInt)tif->tif_rawdatasize; |
| if ((tmsize_t)sp->stream.avail_out != tif->tif_rawdatasize) |
| { |
| TIFFErrorExt(tif->tif_clientdata, module, "ZLib cannot deal with buffers this size"); |
| return (0); |
| } |
| return (deflateReset(&sp->stream) == Z_OK); |
| } |
| |
| static void |
| horizontalDifferenceF(float *ip, int n, int stride, uint16 *wp, uint16 *FromLT2) |
| { |
| int32 r1, g1, b1, a1, r2, g2, b2, a2, mask; |
| float fltsize = Fltsize; |
| |
| #define CLAMP(v) ( (v<(float)0.) ? 0 \ |
| : (v<(float)2.) ? FromLT2[(int)(v*fltsize)] \ |
| : (v>(float)24.2) ? 2047 \ |
| : LogK1*log(v*LogK2) + 0.5 ) |
| |
| mask = CODE_MASK; |
| if (n >= stride) { |
| if (stride == 3) { |
| r2 = wp[0] = (uint16) CLAMP(ip[0]); |
| g2 = wp[1] = (uint16) CLAMP(ip[1]); |
| b2 = wp[2] = (uint16) CLAMP(ip[2]); |
| n -= 3; |
| while (n > 0) { |
| n -= 3; |
| wp += 3; |
| ip += 3; |
| r1 = (int32) CLAMP(ip[0]); wp[0] = (uint16)((r1-r2) & mask); r2 = r1; |
| g1 = (int32) CLAMP(ip[1]); wp[1] = (uint16)((g1-g2) & mask); g2 = g1; |
| b1 = (int32) CLAMP(ip[2]); wp[2] = (uint16)((b1-b2) & mask); b2 = b1; |
| } |
| } else if (stride == 4) { |
| r2 = wp[0] = (uint16) CLAMP(ip[0]); |
| g2 = wp[1] = (uint16) CLAMP(ip[1]); |
| b2 = wp[2] = (uint16) CLAMP(ip[2]); |
| a2 = wp[3] = (uint16) CLAMP(ip[3]); |
| n -= 4; |
| while (n > 0) { |
| n -= 4; |
| wp += 4; |
| ip += 4; |
| r1 = (int32) CLAMP(ip[0]); wp[0] = (uint16)((r1-r2) & mask); r2 = r1; |
| g1 = (int32) CLAMP(ip[1]); wp[1] = (uint16)((g1-g2) & mask); g2 = g1; |
| b1 = (int32) CLAMP(ip[2]); wp[2] = (uint16)((b1-b2) & mask); b2 = b1; |
| a1 = (int32) CLAMP(ip[3]); wp[3] = (uint16)((a1-a2) & mask); a2 = a1; |
| } |
| } else { |
| REPEAT(stride, wp[0] = (uint16) CLAMP(ip[0]); wp++; ip++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[0] = (uint16)(((int32)CLAMP(ip[0])-(int32)CLAMP(ip[-stride])) & mask); |
| wp++; ip++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| static void |
| horizontalDifference16(unsigned short *ip, int n, int stride, |
| unsigned short *wp, uint16 *From14) |
| { |
| register int r1, g1, b1, a1, r2, g2, b2, a2, mask; |
| |
| /* assumption is unsigned pixel values */ |
| #undef CLAMP |
| #define CLAMP(v) From14[(v) >> 2] |
| |
| mask = CODE_MASK; |
| if (n >= stride) { |
| if (stride == 3) { |
| r2 = wp[0] = CLAMP(ip[0]); g2 = wp[1] = CLAMP(ip[1]); |
| b2 = wp[2] = CLAMP(ip[2]); |
| n -= 3; |
| while (n > 0) { |
| n -= 3; |
| wp += 3; |
| ip += 3; |
| r1 = CLAMP(ip[0]); wp[0] = (uint16)((r1-r2) & mask); r2 = r1; |
| g1 = CLAMP(ip[1]); wp[1] = (uint16)((g1-g2) & mask); g2 = g1; |
| b1 = CLAMP(ip[2]); wp[2] = (uint16)((b1-b2) & mask); b2 = b1; |
| } |
| } else if (stride == 4) { |
| r2 = wp[0] = CLAMP(ip[0]); g2 = wp[1] = CLAMP(ip[1]); |
| b2 = wp[2] = CLAMP(ip[2]); a2 = wp[3] = CLAMP(ip[3]); |
| n -= 4; |
| while (n > 0) { |
| n -= 4; |
| wp += 4; |
| ip += 4; |
| r1 = CLAMP(ip[0]); wp[0] = (uint16)((r1-r2) & mask); r2 = r1; |
| g1 = CLAMP(ip[1]); wp[1] = (uint16)((g1-g2) & mask); g2 = g1; |
| b1 = CLAMP(ip[2]); wp[2] = (uint16)((b1-b2) & mask); b2 = b1; |
| a1 = CLAMP(ip[3]); wp[3] = (uint16)((a1-a2) & mask); a2 = a1; |
| } |
| } else { |
| REPEAT(stride, wp[0] = CLAMP(ip[0]); wp++; ip++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[0] = (uint16)((CLAMP(ip[0])-CLAMP(ip[-stride])) & mask); |
| wp++; ip++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| |
| static void |
| horizontalDifference8(unsigned char *ip, int n, int stride, |
| unsigned short *wp, uint16 *From8) |
| { |
| register int r1, g1, b1, a1, r2, g2, b2, a2, mask; |
| |
| #undef CLAMP |
| #define CLAMP(v) (From8[(v)]) |
| |
| mask = CODE_MASK; |
| if (n >= stride) { |
| if (stride == 3) { |
| r2 = wp[0] = CLAMP(ip[0]); g2 = wp[1] = CLAMP(ip[1]); |
| b2 = wp[2] = CLAMP(ip[2]); |
| n -= 3; |
| while (n > 0) { |
| n -= 3; |
| r1 = CLAMP(ip[3]); wp[3] = (uint16)((r1-r2) & mask); r2 = r1; |
| g1 = CLAMP(ip[4]); wp[4] = (uint16)((g1-g2) & mask); g2 = g1; |
| b1 = CLAMP(ip[5]); wp[5] = (uint16)((b1-b2) & mask); b2 = b1; |
| wp += 3; |
| ip += 3; |
| } |
| } else if (stride == 4) { |
| r2 = wp[0] = CLAMP(ip[0]); g2 = wp[1] = CLAMP(ip[1]); |
| b2 = wp[2] = CLAMP(ip[2]); a2 = wp[3] = CLAMP(ip[3]); |
| n -= 4; |
| while (n > 0) { |
| n -= 4; |
| r1 = CLAMP(ip[4]); wp[4] = (uint16)((r1-r2) & mask); r2 = r1; |
| g1 = CLAMP(ip[5]); wp[5] = (uint16)((g1-g2) & mask); g2 = g1; |
| b1 = CLAMP(ip[6]); wp[6] = (uint16)((b1-b2) & mask); b2 = b1; |
| a1 = CLAMP(ip[7]); wp[7] = (uint16)((a1-a2) & mask); a2 = a1; |
| wp += 4; |
| ip += 4; |
| } |
| } else { |
| REPEAT(stride, wp[0] = CLAMP(ip[0]); wp++; ip++) |
| n -= stride; |
| while (n > 0) { |
| REPEAT(stride, |
| wp[0] = (uint16)((CLAMP(ip[0])-CLAMP(ip[-stride])) & mask); |
| wp++; ip++) |
| n -= stride; |
| } |
| } |
| } |
| } |
| |
| /* |
| * Encode a chunk of pixels. |
| */ |
| static int |
| PixarLogEncode(TIFF* tif, uint8* bp, tmsize_t cc, uint16 s) |
| { |
| static const char module[] = "PixarLogEncode"; |
| TIFFDirectory *td = &tif->tif_dir; |
| PixarLogState *sp = EncoderState(tif); |
| tmsize_t i; |
| tmsize_t n; |
| int llen; |
| unsigned short * up; |
| |
| (void) s; |
| |
| switch (sp->user_datafmt) { |
| case PIXARLOGDATAFMT_FLOAT: |
| n = cc / sizeof(float); /* XXX float == 32 bits */ |
| break; |
| case PIXARLOGDATAFMT_16BIT: |
| case PIXARLOGDATAFMT_12BITPICIO: |
| case PIXARLOGDATAFMT_11BITLOG: |
| n = cc / sizeof(uint16); /* XXX uint16 == 16 bits */ |
| break; |
| case PIXARLOGDATAFMT_8BIT: |
| case PIXARLOGDATAFMT_8BITABGR: |
| n = cc; |
| break; |
| default: |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "%d bit input not supported in PixarLog", |
| td->td_bitspersample); |
| return 0; |
| } |
| |
| llen = sp->stride * td->td_imagewidth; |
| /* Check against the number of elements (of size uint16) of sp->tbuf */ |
| if( n > (tmsize_t)(td->td_rowsperstrip * llen) ) |
| { |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "Too many input bytes provided"); |
| return 0; |
| } |
| |
| for (i = 0, up = sp->tbuf; i < n; i += llen, up += llen) { |
| switch (sp->user_datafmt) { |
| case PIXARLOGDATAFMT_FLOAT: |
| horizontalDifferenceF((float *)bp, llen, |
| sp->stride, up, sp->FromLT2); |
| bp += llen * sizeof(float); |
| break; |
| case PIXARLOGDATAFMT_16BIT: |
| horizontalDifference16((uint16 *)bp, llen, |
| sp->stride, up, sp->From14); |
| bp += llen * sizeof(uint16); |
| break; |
| case PIXARLOGDATAFMT_8BIT: |
| horizontalDifference8((unsigned char *)bp, llen, |
| sp->stride, up, sp->From8); |
| bp += llen * sizeof(unsigned char); |
| break; |
| default: |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "%d bit input not supported in PixarLog", |
| td->td_bitspersample); |
| return 0; |
| } |
| } |
| |
| sp->stream.next_in = (unsigned char *) sp->tbuf; |
| assert(sizeof(sp->stream.avail_in)==4); /* if this assert gets raised, |
| we need to simplify this code to reflect a ZLib that is likely updated |
| to deal with 8byte memory sizes, though this code will respond |
| appropriately even before we simplify it */ |
| sp->stream.avail_in = (uInt) (n * sizeof(uint16)); |
| if ((sp->stream.avail_in / sizeof(uint16)) != (uInt) n) |
| { |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "ZLib cannot deal with buffers this size"); |
| return (0); |
| } |
| |
| do { |
| if (deflate(&sp->stream, Z_NO_FLUSH) != Z_OK) { |
| TIFFErrorExt(tif->tif_clientdata, module, "Encoder error: %s", |
| sp->stream.msg ? sp->stream.msg : "(null)"); |
| return (0); |
| } |
| if (sp->stream.avail_out == 0) { |
| tif->tif_rawcc = tif->tif_rawdatasize; |
| TIFFFlushData1(tif); |
| sp->stream.next_out = tif->tif_rawdata; |
| sp->stream.avail_out = (uInt) tif->tif_rawdatasize; /* this is a safe typecast, as check is made already in PixarLogPreEncode */ |
| } |
| } while (sp->stream.avail_in > 0); |
| return (1); |
| } |
| |
| /* |
| * Finish off an encoded strip by flushing the last |
| * string and tacking on an End Of Information code. |
| */ |
| |
| static int |
| PixarLogPostEncode(TIFF* tif) |
| { |
| static const char module[] = "PixarLogPostEncode"; |
| PixarLogState *sp = EncoderState(tif); |
| int state; |
| |
| sp->stream.avail_in = 0; |
| |
| do { |
| state = deflate(&sp->stream, Z_FINISH); |
| switch (state) { |
| case Z_STREAM_END: |
| case Z_OK: |
| if ((tmsize_t)sp->stream.avail_out != tif->tif_rawdatasize) { |
| tif->tif_rawcc = |
| tif->tif_rawdatasize - sp->stream.avail_out; |
| TIFFFlushData1(tif); |
| sp->stream.next_out = tif->tif_rawdata; |
| sp->stream.avail_out = (uInt) tif->tif_rawdatasize; /* this is a safe typecast, as check is made already in PixarLogPreEncode */ |
| } |
| break; |
| default: |
| TIFFErrorExt(tif->tif_clientdata, module, "ZLib error: %s", |
| sp->stream.msg ? sp->stream.msg : "(null)"); |
| return (0); |
| } |
| } while (state != Z_STREAM_END); |
| return (1); |
| } |
| |
| static void |
| PixarLogClose(TIFF* tif) |
| { |
| PixarLogState* sp = (PixarLogState*) tif->tif_data; |
| TIFFDirectory *td = &tif->tif_dir; |
| |
| assert(sp != 0); |
| /* In a really sneaky (and really incorrect, and untruthful, and |
| * troublesome, and error-prone) maneuver that completely goes against |
| * the spirit of TIFF, and breaks TIFF, on close, we covertly |
| * modify both bitspersample and sampleformat in the directory to |
| * indicate 8-bit linear. This way, the decode "just works" even for |
| * readers that don't know about PixarLog, or how to set |
| * the PIXARLOGDATFMT pseudo-tag. |
| */ |
| |
| if (sp->state&PLSTATE_INIT) { |
| /* We test the state to avoid an issue such as in |
| * http://bugzilla.maptools.org/show_bug.cgi?id=2604 |
| * What appends in that case is that the bitspersample is 1 and |
| * a TransferFunction is set. The size of the TransferFunction |
| * depends on 1<<bitspersample. So if we increase it, an access |
| * out of the buffer will happen at directory flushing. |
| * Another option would be to clear those targs. |
| */ |
| td->td_bitspersample = 8; |
| td->td_sampleformat = SAMPLEFORMAT_UINT; |
| } |
| } |
| |
| static void |
| PixarLogCleanup(TIFF* tif) |
| { |
| PixarLogState* sp = (PixarLogState*) tif->tif_data; |
| |
| assert(sp != 0); |
| |
| (void)TIFFPredictorCleanup(tif); |
| |
| tif->tif_tagmethods.vgetfield = sp->vgetparent; |
| tif->tif_tagmethods.vsetfield = sp->vsetparent; |
| |
| if (sp->FromLT2) _TIFFfree(sp->FromLT2); |
| if (sp->From14) _TIFFfree(sp->From14); |
| if (sp->From8) _TIFFfree(sp->From8); |
| if (sp->ToLinearF) _TIFFfree(sp->ToLinearF); |
| if (sp->ToLinear16) _TIFFfree(sp->ToLinear16); |
| if (sp->ToLinear8) _TIFFfree(sp->ToLinear8); |
| if (sp->state&PLSTATE_INIT) { |
| if (tif->tif_mode == O_RDONLY) |
| inflateEnd(&sp->stream); |
| else |
| deflateEnd(&sp->stream); |
| } |
| if (sp->tbuf) |
| _TIFFfree(sp->tbuf); |
| _TIFFfree(sp); |
| tif->tif_data = NULL; |
| |
| _TIFFSetDefaultCompressionState(tif); |
| } |
| |
| static int |
| PixarLogVSetField(TIFF* tif, uint32 tag, va_list ap) |
| { |
| static const char module[] = "PixarLogVSetField"; |
| PixarLogState *sp = (PixarLogState *)tif->tif_data; |
| int result; |
| |
| switch (tag) { |
| case TIFFTAG_PIXARLOGQUALITY: |
| sp->quality = (int) va_arg(ap, int); |
| if (tif->tif_mode != O_RDONLY && (sp->state&PLSTATE_INIT)) { |
| if (deflateParams(&sp->stream, |
| sp->quality, Z_DEFAULT_STRATEGY) != Z_OK) { |
| TIFFErrorExt(tif->tif_clientdata, module, "ZLib error: %s", |
| sp->stream.msg ? sp->stream.msg : "(null)"); |
| return (0); |
| } |
| } |
| return (1); |
| case TIFFTAG_PIXARLOGDATAFMT: |
| sp->user_datafmt = (int) va_arg(ap, int); |
| /* Tweak the TIFF header so that the rest of libtiff knows what |
| * size of data will be passed between app and library, and |
| * assume that the app knows what it is doing and is not |
| * confused by these header manipulations... |
| */ |
| switch (sp->user_datafmt) { |
| case PIXARLOGDATAFMT_8BIT: |
| case PIXARLOGDATAFMT_8BITABGR: |
| TIFFSetField(tif, TIFFTAG_BITSPERSAMPLE, 8); |
| TIFFSetField(tif, TIFFTAG_SAMPLEFORMAT, SAMPLEFORMAT_UINT); |
| break; |
| case PIXARLOGDATAFMT_11BITLOG: |
| TIFFSetField(tif, TIFFTAG_BITSPERSAMPLE, 16); |
| TIFFSetField(tif, TIFFTAG_SAMPLEFORMAT, SAMPLEFORMAT_UINT); |
| break; |
| case PIXARLOGDATAFMT_12BITPICIO: |
| TIFFSetField(tif, TIFFTAG_BITSPERSAMPLE, 16); |
| TIFFSetField(tif, TIFFTAG_SAMPLEFORMAT, SAMPLEFORMAT_INT); |
| break; |
| case PIXARLOGDATAFMT_16BIT: |
| TIFFSetField(tif, TIFFTAG_BITSPERSAMPLE, 16); |
| TIFFSetField(tif, TIFFTAG_SAMPLEFORMAT, SAMPLEFORMAT_UINT); |
| break; |
| case PIXARLOGDATAFMT_FLOAT: |
| TIFFSetField(tif, TIFFTAG_BITSPERSAMPLE, 32); |
| TIFFSetField(tif, TIFFTAG_SAMPLEFORMAT, SAMPLEFORMAT_IEEEFP); |
| break; |
| } |
| /* |
| * Must recalculate sizes should bits/sample change. |
| */ |
| tif->tif_tilesize = isTiled(tif) ? TIFFTileSize(tif) : (tmsize_t)(-1); |
| tif->tif_scanlinesize = TIFFScanlineSize(tif); |
| result = 1; /* NB: pseudo tag */ |
| break; |
| default: |
| result = (*sp->vsetparent)(tif, tag, ap); |
| } |
| return (result); |
| } |
| |
| static int |
| PixarLogVGetField(TIFF* tif, uint32 tag, va_list ap) |
| { |
| PixarLogState *sp = (PixarLogState *)tif->tif_data; |
| |
| switch (tag) { |
| case TIFFTAG_PIXARLOGQUALITY: |
| *va_arg(ap, int*) = sp->quality; |
| break; |
| case TIFFTAG_PIXARLOGDATAFMT: |
| *va_arg(ap, int*) = sp->user_datafmt; |
| break; |
| default: |
| return (*sp->vgetparent)(tif, tag, ap); |
| } |
| return (1); |
| } |
| |
| static const TIFFField pixarlogFields[] = { |
| {TIFFTAG_PIXARLOGDATAFMT, 0, 0, TIFF_ANY, 0, TIFF_SETGET_INT, TIFF_SETGET_UNDEFINED, FIELD_PSEUDO, FALSE, FALSE, "", NULL}, |
| {TIFFTAG_PIXARLOGQUALITY, 0, 0, TIFF_ANY, 0, TIFF_SETGET_INT, TIFF_SETGET_UNDEFINED, FIELD_PSEUDO, FALSE, FALSE, "", NULL} |
| }; |
| |
| int |
| TIFFInitPixarLog(TIFF* tif, int scheme) |
| { |
| static const char module[] = "TIFFInitPixarLog"; |
| |
| PixarLogState* sp; |
| |
| assert(scheme == COMPRESSION_PIXARLOG); |
| |
| /* |
| * Merge codec-specific tag information. |
| */ |
| if (!_TIFFMergeFields(tif, pixarlogFields, |
| TIFFArrayCount(pixarlogFields))) { |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "Merging PixarLog codec-specific tags failed"); |
| return 0; |
| } |
| |
| /* |
| * Allocate state block so tag methods have storage to record values. |
| */ |
| tif->tif_data = (uint8*) _TIFFmalloc(sizeof (PixarLogState)); |
| if (tif->tif_data == NULL) |
| goto bad; |
| sp = (PixarLogState*) tif->tif_data; |
| _TIFFmemset(sp, 0, sizeof (*sp)); |
| sp->stream.data_type = Z_BINARY; |
| sp->user_datafmt = PIXARLOGDATAFMT_UNKNOWN; |
| |
| /* |
| * Install codec methods. |
| */ |
| tif->tif_fixuptags = PixarLogFixupTags; |
| tif->tif_setupdecode = PixarLogSetupDecode; |
| tif->tif_predecode = PixarLogPreDecode; |
| tif->tif_decoderow = PixarLogDecode; |
| tif->tif_decodestrip = PixarLogDecode; |
| tif->tif_decodetile = PixarLogDecode; |
| tif->tif_setupencode = PixarLogSetupEncode; |
| tif->tif_preencode = PixarLogPreEncode; |
| tif->tif_postencode = PixarLogPostEncode; |
| tif->tif_encoderow = PixarLogEncode; |
| tif->tif_encodestrip = PixarLogEncode; |
| tif->tif_encodetile = PixarLogEncode; |
| tif->tif_close = PixarLogClose; |
| tif->tif_cleanup = PixarLogCleanup; |
| |
| /* Override SetField so we can handle our private pseudo-tag */ |
| sp->vgetparent = tif->tif_tagmethods.vgetfield; |
| tif->tif_tagmethods.vgetfield = PixarLogVGetField; /* hook for codec tags */ |
| sp->vsetparent = tif->tif_tagmethods.vsetfield; |
| tif->tif_tagmethods.vsetfield = PixarLogVSetField; /* hook for codec tags */ |
| |
| /* Default values for codec-specific fields */ |
| sp->quality = Z_DEFAULT_COMPRESSION; /* default comp. level */ |
| sp->state = 0; |
| |
| /* we don't wish to use the predictor, |
| * the default is none, which predictor value 1 |
| */ |
| (void) TIFFPredictorInit(tif); |
| |
| /* |
| * build the companding tables |
| */ |
| PixarLogMakeTables(sp); |
| |
| return (1); |
| bad: |
| TIFFErrorExt(tif->tif_clientdata, module, |
| "No space for PixarLog state block"); |
| return (0); |
| } |
| #endif /* PIXARLOG_SUPPORT */ |
| |
| /* vim: set ts=8 sts=8 sw=8 noet: */ |
| /* |
| * Local Variables: |
| * mode: c |
| * c-basic-offset: 8 |
| * fill-column: 78 |
| * End: |
| */ |