///////////////////////////////////////////////////////////////////////////
//
// Copyright (c) 2004, Industrial Light & Magic, a division of Lucas
// Digital Ltd. LLC
//
// All rights reserved.
//
// 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
// from this software without specific prior written permission.
//
// 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.
//
///////////////////////////////////////////////////////////////////////////
//-----------------------------------------------------------------------------
//
// Miscellaneous stuff related to tiled files
//
//-----------------------------------------------------------------------------
#include <ImfTiledMisc.h>
#include "Iex.h"
#include <ImfMisc.h>
#include <ImfChannelList.h>
#include <ImfTileDescription.h>
#include <algorithm>
#include "ImfNamespace.h"
OPENEXR_IMF_INTERNAL_NAMESPACE_SOURCE_ENTER
using IMATH_NAMESPACE::Box2i;
using IMATH_NAMESPACE::V2i;
int
levelSize (int min, int max, int l, LevelRoundingMode rmode)
{
if (l < 0)
throw IEX_NAMESPACE::ArgExc ("Argument not in valid range.");
int a = max - min + 1;
int b = (1 << l);
int size = a / b;
if (rmode == ROUND_UP && size * b < a)
size += 1;
return std::max (size, 1);
}
Box2i
dataWindowForLevel (const TileDescription &tileDesc,
int minX, int maxX,
int minY, int maxY,
int lx, int ly)
{
V2i levelMin = V2i (minX, minY);
V2i levelMax = levelMin +
V2i (levelSize (minX, maxX, lx, tileDesc.roundingMode) - 1,
levelSize (minY, maxY, ly, tileDesc.roundingMode) - 1);
return Box2i(levelMin, levelMax);
}
Box2i
dataWindowForTile (const TileDescription &tileDesc,
int minX, int maxX,
int minY, int maxY,
int dx, int dy,
int lx, int ly)
{
V2i tileMin = V2i (minX + dx * tileDesc.xSize,
minY + dy * tileDesc.ySize);
V2i tileMax = tileMin + V2i (tileDesc.xSize - 1, tileDesc.ySize - 1);
V2i levelMax = dataWindowForLevel
(tileDesc, minX, maxX, minY, maxY, lx, ly).max;
tileMax = V2i (std::min (tileMax[0], levelMax[0]),
std::min (tileMax[1], levelMax[1]));
return Box2i (tileMin, tileMax);
}
size_t
calculateBytesPerPixel (const Header &header)
{
const ChannelList &channels = header.channels();
size_t bytesPerPixel = 0;
for (ChannelList::ConstIterator c = channels.begin();
c != channels.end();
++c)
{
bytesPerPixel += pixelTypeSize (c.channel().type);
}
return bytesPerPixel;
}
void
calculateBytesPerLine (const Header &header,
char* sampleCountBase,
int sampleCountXStride,
int sampleCountYStride,
int minX, int maxX,
int minY, int maxY,
std::vector<int>& xOffsets,
std::vector<int>& yOffsets,
std::vector<Int64>& bytesPerLine)
{
const ChannelList &channels = header.channels();
int pos = 0;
for (ChannelList::ConstIterator c = channels.begin();
c != channels.end();
++c, ++pos)
{
int xOffset = xOffsets[pos];
int yOffset = yOffsets[pos];
int i = 0;
for (int y = minY - yOffset; y <= maxY - yOffset; y++, i++)
for (int x = minX - xOffset; x <= maxX - xOffset; x++)
{
bytesPerLine[i] += sampleCount(sampleCountBase,
sampleCountXStride,
sampleCountYStride,
x, y)
* pixelTypeSize (c.channel().type);
}
}
}
namespace {
int
floorLog2 (int x)
{
//
// For x > 0, floorLog2(y) returns floor(log(x)/log(2)).
//
int y = 0;
while (x > 1)
{
y += 1;
x >>= 1;
}
return y;
}
int
ceilLog2 (int x)
{
//
// For x > 0, ceilLog2(y) returns ceil(log(x)/log(2)).
//
int y = 0;
int r = 0;
while (x > 1)
{
if (x & 1)
r = 1;
y += 1;
x >>= 1;
}
return y + r;
}
int
roundLog2 (int x, LevelRoundingMode rmode)
{
return (rmode == ROUND_DOWN)? floorLog2 (x): ceilLog2 (x);
}
int
calculateNumXLevels (const TileDescription& tileDesc,
int minX, int maxX,
int minY, int maxY)
{
int num = 0;
switch (tileDesc.mode)
{
case ONE_LEVEL:
num = 1;
break;
case MIPMAP_LEVELS:
{
int w = maxX - minX + 1;
int h = maxY - minY + 1;
num = roundLog2 (std::max (w, h), tileDesc.roundingMode) + 1;
}
break;
case RIPMAP_LEVELS:
{
int w = maxX - minX + 1;
num = roundLog2 (w, tileDesc.roundingMode) + 1;
}
break;
default:
throw IEX_NAMESPACE::ArgExc ("Unknown LevelMode format.");
}
return num;
}
int
calculateNumYLevels (const TileDescription& tileDesc,
int minX, int maxX,
int minY, int maxY)
{
int num = 0;
switch (tileDesc.mode)
{
case ONE_LEVEL:
num = 1;
break;
case MIPMAP_LEVELS:
{
int w = maxX - minX + 1;
int h = maxY - minY + 1;
num = roundLog2 (std::max (w, h), tileDesc.roundingMode) + 1;
}
break;
case RIPMAP_LEVELS:
{
int h = maxY - minY + 1;
num = roundLog2 (h, tileDesc.roundingMode) + 1;
}
break;
default:
throw IEX_NAMESPACE::ArgExc ("Unknown LevelMode format.");
}
return num;
}
void
calculateNumTiles (int *numTiles,
int numLevels,
int min, int max,
int size,
LevelRoundingMode rmode)
{
for (int i = 0; i < numLevels; i++)
{
numTiles[i] = (levelSize (min, max, i, rmode) + size - 1) / size;
}
}
} // namespace
void
precalculateTileInfo (const TileDescription& tileDesc,
int minX, int maxX,
int minY, int maxY,
int *&numXTiles, int *&numYTiles,
int &numXLevels, int &numYLevels)
{
numXLevels = calculateNumXLevels(tileDesc, minX, maxX, minY, maxY);
numYLevels = calculateNumYLevels(tileDesc, minX, maxX, minY, maxY);
numXTiles = new int[numXLevels];
numYTiles = new int[numYLevels];
calculateNumTiles (numXTiles,
numXLevels,
minX, maxX,
tileDesc.xSize,
tileDesc.roundingMode);
calculateNumTiles (numYTiles,
numYLevels,
minY, maxY,
tileDesc.ySize,
tileDesc.roundingMode);
}
int
getTiledChunkOffsetTableSize(const Header& header)
{
//
// Save the dataWindow information
//
const Box2i &dataWindow = header.dataWindow();
//
// Precompute level and tile information.
//
int* numXTiles;
int* numYTiles;
int numXLevels;
int numYLevels;
precalculateTileInfo (header.tileDescription(),
dataWindow.min.x, dataWindow.max.x,
dataWindow.min.y, dataWindow.max.y,
numXTiles, numYTiles,
numXLevels, numYLevels);
//
// Calculate lineOffsetSize.
//
int lineOffsetSize = 0;
const TileDescription &desc = header.tileDescription();
switch (desc.mode)
{
case ONE_LEVEL:
case MIPMAP_LEVELS:
for (int i = 0; i < numXLevels; i++)
lineOffsetSize += numXTiles[i] * numYTiles[i];
break;
case RIPMAP_LEVELS:
for (int i = 0; i < numXLevels; i++)
for (int j = 0; j < numYLevels; j++)
lineOffsetSize += numXTiles[i] * numYTiles[j];
break;
case NUM_LEVELMODES :
throw IEX_NAMESPACE::LogicExc("Bad level mode getting chunk offset table size");
}
delete[] numXTiles;
delete[] numYTiles;
return lineOffsetSize;
}
OPENEXR_IMF_INTERNAL_NAMESPACE_SOURCE_EXIT