opencv/3rdparty/openexr/IlmImf/ImfB44Compressor.cpp

1070 lines
26 KiB
C++
Raw Normal View History

2012-08-25 04:31:49 +08:00
///////////////////////////////////////////////////////////////////////////
//
// Copyright (c) 2006, Industrial Light & Magic, a division of Lucas
// Digital Ltd. LLC
2012-10-17 07:18:30 +08:00
//
2012-08-25 04:31:49 +08:00
// All rights reserved.
2012-10-17 07:18:30 +08:00
//
2012-08-25 04:31:49 +08:00
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are
// met:
// * Redistributions of source code must retain the above copyright
// notice, this list of conditions and the following disclaimer.
// * Redistributions in binary form must reproduce the above
// copyright notice, this list of conditions and the following disclaimer
// in the documentation and/or other materials provided with the
// distribution.
// * Neither the name of Industrial Light & Magic nor the names of
// its contributors may be used to endorse or promote products derived
2012-10-17 07:18:30 +08:00
// from this software without specific prior written permission.
//
2012-08-25 04:31:49 +08:00
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
// OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
// LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
// DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
// OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
//
///////////////////////////////////////////////////////////////////////////
//-----------------------------------------------------------------------------
//
// class B44Compressor
//
// This compressor is lossy for HALF channels; the compression rate
// is fixed at 32/14 (approximately 2.28). FLOAT and UINT channels
// are not compressed; their data are preserved exactly.
//
// Each HALF channel is split into blocks of 4 by 4 pixels. An
// uncompressed block occupies 32 bytes, which are re-interpreted
// as sixteen 16-bit unsigned integers, t[0] ... t[15]. Compression
// shrinks the block to 14 bytes. The compressed 14-byte block
// contains
//
// - t[0]
//
// - a 6-bit shift value
//
// - 15 densely packed 6-bit values, r[0] ... r[14], which are
// computed by subtracting adjacent pixel values and right-
// shifting the differences according to the stored shift value.
//
// Differences between adjacent pixels are computed according
// to the following diagram:
//
// 0 --------> 1 --------> 2 --------> 3
// | 3 7 11
// |
// | 0
// |
2012-10-17 07:18:30 +08:00
// v
2012-08-25 04:31:49 +08:00
// 4 --------> 5 --------> 6 --------> 7
// | 4 8 12
// |
// | 1
// |
// v
// 8 --------> 9 --------> 10 --------> 11
// | 5 9 13
// |
// | 2
// |
// v
// 12 --------> 13 --------> 14 --------> 15
// 6 10 14
//
// Here
//
// 5 ---------> 6
// 8
//
// means that r[8] is the difference between t[5] and t[6].
//
// - optionally, a 4-by-4 pixel block where all pixels have the
// same value can be treated as a special case, where the
// compressed block contains only 3 instead of 14 bytes:
// t[0], followed by an "impossible" 6-bit shift value and
// two padding bits.
//
// This compressor can handle positive and negative pixel values.
// NaNs and infinities are replaced with zeroes before compression.
//
//-----------------------------------------------------------------------------
#include <ImfB44Compressor.h>
#include <ImfHeader.h>
#include <ImfChannelList.h>
#include <ImfMisc.h>
#include <ImfCheckedArithmetic.h>
#include <ImathFun.h>
#include <ImathBox.h>
#include <Iex.h>
#include <ImfIO.h>
#include <ImfXdr.h>
#include <string.h>
#include <assert.h>
#include <algorithm>
namespace Imf {
using Imath::divp;
using Imath::modp;
using Imath::Box2i;
using Imath::V2i;
using std::min;
namespace {
//
// Lookup tables for
// y = exp (x / 8)
2012-10-17 07:18:30 +08:00
// and
2012-08-25 04:31:49 +08:00
// x = 8 * log (y)
//
#include "b44ExpLogTable.h"
inline void
convertFromLinear (unsigned short s[16])
{
for (int i = 0; i < 16; ++i)
2012-10-17 07:18:30 +08:00
s[i] = expTable[s[i]];
2012-08-25 04:31:49 +08:00
}
inline void
convertToLinear (unsigned short s[16])
{
for (int i = 0; i < 16; ++i)
2012-10-17 07:18:30 +08:00
s[i] = logTable[s[i]];
2012-08-25 04:31:49 +08:00
}
inline int
shiftAndRound (int x, int shift)
{
//
// Compute
//
// y = x * pow (2, -shift),
//
// then round y to the nearest integer.
// In case of a tie, where y is exactly
// halfway between two integers, round
// to the even one.
//
x <<= 1;
int a = (1 << shift) - 1;
shift += 1;
int b = (x >> shift) & 1;
return (x + a + b) >> shift;
}
int
pack (const unsigned short s[16],
unsigned char b[14],
bool optFlatFields,
bool exactMax)
{
//
// Pack a block of 4 by 4 16-bit pixels (32 bytes) into
// either 14 or 3 bytes.
//
//
// Integers s[0] ... s[15] represent floating-point numbers
// in what is essentially a sign-magnitude format. Convert
// s[0] .. s[15] into a new set of integers, t[0] ... t[15],
// such that if t[i] is greater than t[j], the floating-point
// number that corresponds to s[i] is always greater than
// the floating-point number that corresponds to s[j].
//
// Also, replace any bit patterns that represent NaNs or
// infinities with bit patterns that represent floating-point
// zeroes.
//
// bit pattern floating-point bit pattern
// in s[i] value in t[i]
//
// 0x7fff NAN 0x8000
// 0x7ffe NAN 0x8000
// ... ...
// 0x7c01 NAN 0x8000
// 0x7c00 +infinity 0x8000
// 0x7bff +HALF_MAX 0xfbff
// 0x7bfe 0xfbfe
// 0x7bfd 0xfbfd
// ... ...
// 0x0002 +2 * HALF_MIN 0x8002
// 0x0001 +HALF_MIN 0x8001
// 0x0000 +0.0 0x8000
// 0x8000 -0.0 0x7fff
// 0x8001 -HALF_MIN 0x7ffe
// 0x8002 -2 * HALF_MIN 0x7ffd
// ... ...
// 0xfbfd 0x0f02
// 0xfbfe 0x0401
// 0xfbff -HALF_MAX 0x0400
// 0xfc00 -infinity 0x8000
// 0xfc01 NAN 0x8000
// ... ...
// 0xfffe NAN 0x8000
// 0xffff NAN 0x8000
//
unsigned short t[16];
for (int i = 0; i < 16; ++i)
{
2012-10-17 07:18:30 +08:00
if ((s[i] & 0x7c00) == 0x7c00)
t[i] = 0x8000;
else if (s[i] & 0x8000)
t[i] = ~s[i];
else
t[i] = s[i] | 0x8000;
2012-08-25 04:31:49 +08:00
}
2012-10-17 07:18:30 +08:00
2012-08-25 04:31:49 +08:00
//
// Find the maximum, tMax, of t[0] ... t[15].
//
unsigned short tMax = 0;
for (int i = 0; i < 16; ++i)
2012-10-17 07:18:30 +08:00
if (tMax < t[i])
tMax = t[i];
2012-08-25 04:31:49 +08:00
//
// Compute a set of running differences, r[0] ... r[14]:
// Find a shift value such that after rounding off the
// rightmost bits and shifting all differenes are between
// -32 and +31. Then bias the differences so that they
// end up between 0 and 63.
//
int shift = -1;
int d[16];
int r[15];
int rMin;
int rMax;
const int bias = 0x20;
do
{
2012-10-17 07:18:30 +08:00
shift += 1;
//
// Compute absolute differences, d[0] ... d[15],
// between tMax and t[0] ... t[15].
//
// Shift and round the absolute differences.
//
for (int i = 0; i < 16; ++i)
d[i] = shiftAndRound (tMax - t[i], shift);
//
// Convert d[0] .. d[15] into running differences
//
r[ 0] = d[ 0] - d[ 4] + bias;
r[ 1] = d[ 4] - d[ 8] + bias;
r[ 2] = d[ 8] - d[12] + bias;
r[ 3] = d[ 0] - d[ 1] + bias;
r[ 4] = d[ 4] - d[ 5] + bias;
r[ 5] = d[ 8] - d[ 9] + bias;
r[ 6] = d[12] - d[13] + bias;
r[ 7] = d[ 1] - d[ 2] + bias;
r[ 8] = d[ 5] - d[ 6] + bias;
r[ 9] = d[ 9] - d[10] + bias;
r[10] = d[13] - d[14] + bias;
r[11] = d[ 2] - d[ 3] + bias;
r[12] = d[ 6] - d[ 7] + bias;
r[13] = d[10] - d[11] + bias;
r[14] = d[14] - d[15] + bias;
rMin = r[0];
rMax = r[0];
for (int i = 1; i < 15; ++i)
{
if (rMin > r[i])
rMin = r[i];
if (rMax < r[i])
rMax = r[i];
}
2012-08-25 04:31:49 +08:00
}
while (rMin < 0 || rMax > 0x3f);
if (rMin == bias && rMax == bias && optFlatFields)
{
2012-10-17 07:18:30 +08:00
//
// Special case - all pixels have the same value.
// We encode this in 3 instead of 14 bytes by
// storing the value 0xfc in the third output byte,
// which cannot occur in the 14-byte encoding.
//
b[0] = (t[0] >> 8);
b[1] = t[0];
b[2] = 0xfc;
return 3;
2012-08-25 04:31:49 +08:00
}
if (exactMax)
{
2012-10-17 07:18:30 +08:00
//
// Adjust t[0] so that the pixel whose value is equal
// to tMax gets represented as accurately as possible.
//
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
t[0] = tMax - (d[0] << shift);
2012-08-25 04:31:49 +08:00
}
//
// Pack t[0], shift and r[0] ... r[14] into 14 bytes:
//
b[ 0] = (t[0] >> 8);
b[ 1] = t[0];
b[ 2] = (unsigned char) ((shift << 2) | (r[ 0] >> 4));
b[ 3] = (unsigned char) ((r[ 0] << 4) | (r[ 1] >> 2));
b[ 4] = (unsigned char) ((r[ 1] << 6) | r[ 2] );
b[ 5] = (unsigned char) ((r[ 3] << 2) | (r[ 4] >> 4));
b[ 6] = (unsigned char) ((r[ 4] << 4) | (r[ 5] >> 2));
b[ 7] = (unsigned char) ((r[ 5] << 6) | r[ 6] );
b[ 8] = (unsigned char) ((r[ 7] << 2) | (r[ 8] >> 4));
b[ 9] = (unsigned char) ((r[ 8] << 4) | (r[ 9] >> 2));
b[10] = (unsigned char) ((r[ 9] << 6) | r[10] );
b[11] = (unsigned char) ((r[11] << 2) | (r[12] >> 4));
b[12] = (unsigned char) ((r[12] << 4) | (r[13] >> 2));
b[13] = (unsigned char) ((r[13] << 6) | r[14] );
return 14;
}
inline
void
unpack14 (const unsigned char b[14], unsigned short s[16])
{
//
// Unpack a 14-byte block into 4 by 4 16-bit pixels.
//
#if defined (DEBUG)
2012-10-17 07:18:30 +08:00
assert (b[2] != 0xfc);
2012-08-25 04:31:49 +08:00
#endif
s[ 0] = (b[0] << 8) | b[1];
unsigned short shift = (b[ 2] >> 2);
unsigned short bias = (0x20 << shift);
s[ 4] = s[ 0] + ((((b[ 2] << 4) | (b[ 3] >> 4)) & 0x3f) << shift) - bias;
s[ 8] = s[ 4] + ((((b[ 3] << 2) | (b[ 4] >> 6)) & 0x3f) << shift) - bias;
s[12] = s[ 8] + ((b[ 4] & 0x3f) << shift) - bias;
2012-10-17 07:18:30 +08:00
2012-08-25 04:31:49 +08:00
s[ 1] = s[ 0] + ((b[ 5] >> 2) << shift) - bias;
s[ 5] = s[ 4] + ((((b[ 5] << 4) | (b[ 6] >> 4)) & 0x3f) << shift) - bias;
s[ 9] = s[ 8] + ((((b[ 6] << 2) | (b[ 7] >> 6)) & 0x3f) << shift) - bias;
s[13] = s[12] + ((b[ 7] & 0x3f) << shift) - bias;
2012-10-17 07:18:30 +08:00
2012-08-25 04:31:49 +08:00
s[ 2] = s[ 1] + ((b[ 8] >> 2) << shift) - bias;
s[ 6] = s[ 5] + ((((b[ 8] << 4) | (b[ 9] >> 4)) & 0x3f) << shift) - bias;
s[10] = s[ 9] + ((((b[ 9] << 2) | (b[10] >> 6)) & 0x3f) << shift) - bias;
s[14] = s[13] + ((b[10] & 0x3f) << shift) - bias;
2012-10-17 07:18:30 +08:00
2012-08-25 04:31:49 +08:00
s[ 3] = s[ 2] + ((b[11] >> 2) << shift) - bias;
s[ 7] = s[ 6] + ((((b[11] << 4) | (b[12] >> 4)) & 0x3f) << shift) - bias;
s[11] = s[10] + ((((b[12] << 2) | (b[13] >> 6)) & 0x3f) << shift) - bias;
s[15] = s[14] + ((b[13] & 0x3f) << shift) - bias;
for (int i = 0; i < 16; ++i)
{
2012-10-17 07:18:30 +08:00
if (s[i] & 0x8000)
s[i] &= 0x7fff;
else
s[i] = ~s[i];
2012-08-25 04:31:49 +08:00
}
}
inline
void
unpack3 (const unsigned char b[3], unsigned short s[16])
{
//
// Unpack a 3-byte block into 4 by 4 identical 16-bit pixels.
//
#if defined (DEBUG)
2012-10-17 07:18:30 +08:00
assert (b[2] == 0xfc);
2012-08-25 04:31:49 +08:00
#endif
s[0] = (b[0] << 8) | b[1];
if (s[0] & 0x8000)
2012-10-17 07:18:30 +08:00
s[0] &= 0x7fff;
2012-08-25 04:31:49 +08:00
else
2012-10-17 07:18:30 +08:00
s[0] = ~s[0];
2012-08-25 04:31:49 +08:00
for (int i = 1; i < 16; ++i)
2012-10-17 07:18:30 +08:00
s[i] = s[0];
2012-08-25 04:31:49 +08:00
}
void
notEnoughData ()
{
throw Iex::InputExc ("Error decompressing data "
2012-10-17 07:18:30 +08:00
"(input data are shorter than expected).");
2012-08-25 04:31:49 +08:00
}
void
tooMuchData ()
{
throw Iex::InputExc ("Error decompressing data "
2012-10-17 07:18:30 +08:00
"(input data are longer than expected).");
2012-08-25 04:31:49 +08:00
}
} // namespace
struct B44Compressor::ChannelData
{
unsigned short * start;
unsigned short * end;
int nx;
int ny;
int ys;
PixelType type;
bool pLinear;
int size;
};
B44Compressor::B44Compressor
(const Header &hdr,
size_t maxScanLineSize,
size_t numScanLines,
bool optFlatFields)
:
Compressor (hdr),
_maxScanLineSize (maxScanLineSize),
_optFlatFields (optFlatFields),
_format (XDR),
_numScanLines (numScanLines),
_tmpBuffer (0),
_outBuffer (0),
_numChans (0),
_channels (hdr.channels()),
_channelData (0)
{
//
// Allocate buffers for compressed an uncompressed pixel data,
// allocate a set of ChannelData structs to help speed up the
// compress() and uncompress() functions, below, and determine
// if uncompressed pixel data should be in native or Xdr format.
//
_tmpBuffer = new unsigned short
[checkArraySize (uiMult (maxScanLineSize, numScanLines),
sizeof (unsigned short))];
const ChannelList &channels = header().channels();
int numHalfChans = 0;
for (ChannelList::ConstIterator c = channels.begin();
2012-10-17 07:18:30 +08:00
c != channels.end();
++c)
2012-08-25 04:31:49 +08:00
{
2012-10-17 07:18:30 +08:00
assert (pixelTypeSize (c.channel().type) % pixelTypeSize (HALF) == 0);
++_numChans;
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
if (c.channel().type == HALF)
++numHalfChans;
2012-08-25 04:31:49 +08:00
}
//
// Compressed data may be larger than the input data
//
size_t padding = 12 * numHalfChans * (numScanLines + 3) / 4;
_outBuffer = new char
[uiAdd (uiMult (maxScanLineSize, numScanLines), padding)];
_channelData = new ChannelData[_numChans];
int i = 0;
for (ChannelList::ConstIterator c = channels.begin();
2012-10-17 07:18:30 +08:00
c != channels.end();
++c, ++i)
2012-08-25 04:31:49 +08:00
{
2012-10-17 07:18:30 +08:00
_channelData[i].ys = c.channel().ySampling;
_channelData[i].type = c.channel().type;
_channelData[i].pLinear = c.channel().pLinear;
_channelData[i].size =
pixelTypeSize (c.channel().type) / pixelTypeSize (HALF);
2012-08-25 04:31:49 +08:00
}
const Box2i &dataWindow = hdr.dataWindow();
_minX = dataWindow.min.x;
_maxX = dataWindow.max.x;
_maxY = dataWindow.max.y;
//
// We can support uncompressed data in the machine's native
// format only if all image channels are of type HALF.
//
assert (sizeof (unsigned short) == pixelTypeSize (HALF));
if (_numChans == numHalfChans)
2012-10-17 07:18:30 +08:00
_format = NATIVE;
2012-08-25 04:31:49 +08:00
}
B44Compressor::~B44Compressor ()
{
delete [] _tmpBuffer;
delete [] _outBuffer;
delete [] _channelData;
}
int
B44Compressor::numScanLines () const
{
return _numScanLines;
}
Compressor::Format
B44Compressor::format () const
{
return _format;
}
int
B44Compressor::compress (const char *inPtr,
2012-10-17 07:18:30 +08:00
int inSize,
int minY,
const char *&outPtr)
2012-08-25 04:31:49 +08:00
{
return compress (inPtr,
2012-10-17 07:18:30 +08:00
inSize,
Box2i (V2i (_minX, minY),
V2i (_maxX, minY + numScanLines() - 1)),
outPtr);
2012-08-25 04:31:49 +08:00
}
int
B44Compressor::compressTile (const char *inPtr,
2012-10-17 07:18:30 +08:00
int inSize,
Imath::Box2i range,
const char *&outPtr)
2012-08-25 04:31:49 +08:00
{
return compress (inPtr, inSize, range, outPtr);
}
int
B44Compressor::uncompress (const char *inPtr,
2012-10-17 07:18:30 +08:00
int inSize,
int minY,
const char *&outPtr)
2012-08-25 04:31:49 +08:00
{
return uncompress (inPtr,
2012-10-17 07:18:30 +08:00
inSize,
Box2i (V2i (_minX, minY),
V2i (_maxX, minY + numScanLines() - 1)),
outPtr);
2012-08-25 04:31:49 +08:00
}
int
B44Compressor::uncompressTile (const char *inPtr,
2012-10-17 07:18:30 +08:00
int inSize,
Imath::Box2i range,
const char *&outPtr)
2012-08-25 04:31:49 +08:00
{
return uncompress (inPtr, inSize, range, outPtr);
}
int
B44Compressor::compress (const char *inPtr,
2012-10-17 07:18:30 +08:00
int inSize,
Imath::Box2i range,
const char *&outPtr)
2012-08-25 04:31:49 +08:00
{
//
// Compress a block of pixel data: First copy the input pixels
// from the input buffer into _tmpBuffer, rearranging them such
// that blocks of 4x4 pixels of a single channel can be accessed
// conveniently. Then compress each 4x4 block of HALF pixel data
// and append the result to the output buffer. Copy UINT and
// FLOAT data to the output buffer without compressing them.
//
outPtr = _outBuffer;
if (inSize == 0)
{
2012-10-17 07:18:30 +08:00
//
// Special case - empty input buffer.
//
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
return 0;
2012-08-25 04:31:49 +08:00
}
//
// For each channel, detemine how many pixels are stored
// in the input buffer, and where those pixels will be
// placed in _tmpBuffer.
//
int minX = range.min.x;
int maxX = min (range.max.x, _maxX);
int minY = range.min.y;
int maxY = min (range.max.y, _maxY);
2012-10-17 07:18:30 +08:00
2012-08-25 04:31:49 +08:00
unsigned short *tmpBufferEnd = _tmpBuffer;
int i = 0;
for (ChannelList::ConstIterator c = _channels.begin();
2012-10-17 07:18:30 +08:00
c != _channels.end();
++c, ++i)
2012-08-25 04:31:49 +08:00
{
2012-10-17 07:18:30 +08:00
ChannelData &cd = _channelData[i];
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
cd.start = tmpBufferEnd;
cd.end = cd.start;
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
cd.nx = numSamples (c.channel().xSampling, minX, maxX);
cd.ny = numSamples (c.channel().ySampling, minY, maxY);
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
tmpBufferEnd += cd.nx * cd.ny * cd.size;
2012-08-25 04:31:49 +08:00
}
if (_format == XDR)
{
2012-10-17 07:18:30 +08:00
//
// The data in the input buffer are in the machine-independent
// Xdr format. Copy the HALF channels into _tmpBuffer and
// convert them back into native format for compression.
// Copy UINT and FLOAT channels verbatim into _tmpBuffer.
//
for (int y = minY; y <= maxY; ++y)
{
for (int i = 0; i < _numChans; ++i)
{
ChannelData &cd = _channelData[i];
if (modp (y, cd.ys) != 0)
continue;
if (cd.type == HALF)
{
for (int x = cd.nx; x > 0; --x)
{
Xdr::read <CharPtrIO> (inPtr, *cd.end);
++cd.end;
}
}
else
{
int n = cd.nx * cd.size;
memcpy (cd.end, inPtr, n * sizeof (unsigned short));
inPtr += n * sizeof (unsigned short);
cd.end += n;
}
}
}
2012-08-25 04:31:49 +08:00
}
else
{
2012-10-17 07:18:30 +08:00
//
// The input buffer contains only HALF channels, and they
// are in native, machine-dependent format. Copy the pixels
// into _tmpBuffer.
//
for (int y = minY; y <= maxY; ++y)
{
for (int i = 0; i < _numChans; ++i)
{
ChannelData &cd = _channelData[i];
#if defined (DEBUG)
assert (cd.type == HALF);
#endif
if (modp (y, cd.ys) != 0)
continue;
int n = cd.nx * cd.size;
memcpy (cd.end, inPtr, n * sizeof (unsigned short));
inPtr += n * sizeof (unsigned short);
cd.end += n;
}
}
2012-08-25 04:31:49 +08:00
}
//
// The pixels for each channel have been packed into a contiguous
// block in _tmpBuffer. HALF channels are in native format; UINT
// and FLOAT channels are in Xdr format.
//
#if defined (DEBUG)
2012-10-17 07:18:30 +08:00
for (int i = 1; i < _numChans; ++i)
assert (_channelData[i-1].end == _channelData[i].start);
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
assert (_channelData[_numChans-1].end == tmpBufferEnd);
2012-08-25 04:31:49 +08:00
#endif
//
// For each HALF channel, split the data in _tmpBuffer into 4x4
// pixel blocks. Compress each block and append the compressed
// data to the output buffer.
//
// UINT and FLOAT channels are copied from _tmpBuffer into the
// output buffer without further processing.
//
char *outEnd = _outBuffer;
for (int i = 0; i < _numChans; ++i)
{
2012-10-17 07:18:30 +08:00
ChannelData &cd = _channelData[i];
if (cd.type != HALF)
{
//
// UINT or FLOAT channel.
//
int n = cd.nx * cd.ny * cd.size * sizeof (unsigned short);
memcpy (outEnd, cd.start, n);
outEnd += n;
continue;
}
//
// HALF channel
//
for (int y = 0; y < cd.ny; y += 4)
{
//
// Copy the next 4x4 pixel block into array s.
// If the width, cd.nx, or the height, cd.ny, of
// the pixel data in _tmpBuffer is not divisible
// by 4, then pad the data by repeating the
// rightmost column and the bottom row.
//
unsigned short *row0 = cd.start + y * cd.nx;
unsigned short *row1 = row0 + cd.nx;
unsigned short *row2 = row1 + cd.nx;
unsigned short *row3 = row2 + cd.nx;
if (y + 3 >= cd.ny)
{
if (y + 1 >= cd.ny)
row1 = row0;
if (y + 2 >= cd.ny)
row2 = row1;
row3 = row2;
}
for (int x = 0; x < cd.nx; x += 4)
{
unsigned short s[16];
if (x + 3 >= cd.nx)
{
int n = cd.nx - x;
for (int i = 0; i < 4; ++i)
{
int j = min (i, n - 1);
s[i + 0] = row0[j];
s[i + 4] = row1[j];
s[i + 8] = row2[j];
s[i + 12] = row3[j];
}
}
else
{
memcpy (&s[ 0], row0, 4 * sizeof (unsigned short));
memcpy (&s[ 4], row1, 4 * sizeof (unsigned short));
memcpy (&s[ 8], row2, 4 * sizeof (unsigned short));
memcpy (&s[12], row3, 4 * sizeof (unsigned short));
}
row0 += 4;
row1 += 4;
row2 += 4;
row3 += 4;
//
// Compress the contents of array s and append the
// results to the output buffer.
//
if (cd.pLinear)
convertFromLinear (s);
outEnd += pack (s, (unsigned char *) outEnd,
_optFlatFields, !cd.pLinear);
}
}
2012-08-25 04:31:49 +08:00
}
return outEnd - _outBuffer;
}
int
B44Compressor::uncompress (const char *inPtr,
2012-10-17 07:18:30 +08:00
int inSize,
Imath::Box2i range,
const char *&outPtr)
2012-08-25 04:31:49 +08:00
{
//
// This function is the reverse of the compress() function,
// above. First all pixels are moved from the input buffer
// into _tmpBuffer. UINT and FLOAT channels are copied
// verbatim; HALF channels are uncompressed in blocks of
// 4x4 pixels. Then the pixels in _tmpBuffer are copied
// into the output buffer and rearranged such that the data
// for for each scan line form a contiguous block.
//
outPtr = _outBuffer;
if (inSize == 0)
{
2012-10-17 07:18:30 +08:00
return 0;
2012-08-25 04:31:49 +08:00
}
int minX = range.min.x;
int maxX = min (range.max.x, _maxX);
int minY = range.min.y;
int maxY = min (range.max.y, _maxY);
2012-10-17 07:18:30 +08:00
2012-08-25 04:31:49 +08:00
unsigned short *tmpBufferEnd = _tmpBuffer;
int i = 0;
for (ChannelList::ConstIterator c = _channels.begin();
2012-10-17 07:18:30 +08:00
c != _channels.end();
++c, ++i)
2012-08-25 04:31:49 +08:00
{
2012-10-17 07:18:30 +08:00
ChannelData &cd = _channelData[i];
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
cd.start = tmpBufferEnd;
cd.end = cd.start;
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
cd.nx = numSamples (c.channel().xSampling, minX, maxX);
cd.ny = numSamples (c.channel().ySampling, minY, maxY);
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
tmpBufferEnd += cd.nx * cd.ny * cd.size;
2012-08-25 04:31:49 +08:00
}
for (int i = 0; i < _numChans; ++i)
{
2012-10-17 07:18:30 +08:00
ChannelData &cd = _channelData[i];
if (cd.type != HALF)
{
//
// UINT or FLOAT channel.
//
int n = cd.nx * cd.ny * cd.size * sizeof (unsigned short);
if (inSize < n)
notEnoughData();
memcpy (cd.start, inPtr, n);
inPtr += n;
inSize -= n;
continue;
}
//
// HALF channel
//
for (int y = 0; y < cd.ny; y += 4)
{
unsigned short *row0 = cd.start + y * cd.nx;
unsigned short *row1 = row0 + cd.nx;
unsigned short *row2 = row1 + cd.nx;
unsigned short *row3 = row2 + cd.nx;
for (int x = 0; x < cd.nx; x += 4)
{
unsigned short s[16];
if (inSize < 3)
notEnoughData();
if (((const unsigned char *)inPtr)[2] == 0xfc)
{
unpack3 ((const unsigned char *)inPtr, s);
inPtr += 3;
inSize -= 3;
}
else
{
if (inSize < 14)
notEnoughData();
unpack14 ((const unsigned char *)inPtr, s);
inPtr += 14;
inSize -= 14;
}
if (cd.pLinear)
convertToLinear (s);
int n = (x + 3 < cd.nx)?
4 * sizeof (unsigned short) :
(cd.nx - x) * sizeof (unsigned short);
if (y + 3 < cd.ny)
{
memcpy (row0, &s[ 0], n);
memcpy (row1, &s[ 4], n);
memcpy (row2, &s[ 8], n);
memcpy (row3, &s[12], n);
}
else
{
memcpy (row0, &s[ 0], n);
if (y + 1 < cd.ny)
memcpy (row1, &s[ 4], n);
if (y + 2 < cd.ny)
memcpy (row2, &s[ 8], n);
}
row0 += 4;
row1 += 4;
row2 += 4;
row3 += 4;
}
}
2012-08-25 04:31:49 +08:00
}
char *outEnd = _outBuffer;
if (_format == XDR)
{
2012-10-17 07:18:30 +08:00
for (int y = minY; y <= maxY; ++y)
{
for (int i = 0; i < _numChans; ++i)
{
ChannelData &cd = _channelData[i];
if (modp (y, cd.ys) != 0)
continue;
if (cd.type == HALF)
{
for (int x = cd.nx; x > 0; --x)
{
Xdr::write <CharPtrIO> (outEnd, *cd.end);
++cd.end;
}
}
else
{
int n = cd.nx * cd.size;
memcpy (outEnd, cd.end, n * sizeof (unsigned short));
outEnd += n * sizeof (unsigned short);
cd.end += n;
}
}
}
2012-08-25 04:31:49 +08:00
}
else
{
2012-10-17 07:18:30 +08:00
for (int y = minY; y <= maxY; ++y)
{
for (int i = 0; i < _numChans; ++i)
{
ChannelData &cd = _channelData[i];
#if defined (DEBUG)
assert (cd.type == HALF);
#endif
if (modp (y, cd.ys) != 0)
continue;
int n = cd.nx * cd.size;
memcpy (outEnd, cd.end, n * sizeof (unsigned short));
outEnd += n * sizeof (unsigned short);
cd.end += n;
}
}
2012-08-25 04:31:49 +08:00
}
#if defined (DEBUG)
2012-10-17 07:18:30 +08:00
for (int i = 1; i < _numChans; ++i)
assert (_channelData[i-1].end == _channelData[i].start);
2012-08-25 04:31:49 +08:00
2012-10-17 07:18:30 +08:00
assert (_channelData[_numChans-1].end == tmpBufferEnd);
2012-08-25 04:31:49 +08:00
#endif
if (inSize > 0)
2012-10-17 07:18:30 +08:00
tooMuchData();
2012-08-25 04:31:49 +08:00
outPtr = _outBuffer;
return outEnd - _outBuffer;
}
} // namespace Imf