Logo Search packages:      
Sourcecode: openexr version File versions  Download package

testCompression.cpp

///////////////////////////////////////////////////////////////////////////
//
// 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.
//
///////////////////////////////////////////////////////////////////////////


#include <tmpDir.h>

#include <ImfOutputFile.h>
#include <ImfInputFile.h>
#include <ImfChannelList.h>
#include <ImfArray.h>
#include <ImathRandom.h>
#include <half.h>
#include <compareFloat.h>

#include <stdio.h>
#include <assert.h>

using namespace std;
using namespace Imath;
using namespace Imf;

namespace {

void
fillPixels1 (Array2D<unsigned int> &pi,
           Array2D<half> &ph,
           Array2D<float> &pf,
           int width,
           int height)
{
    cout << "only zeroes" << endl;

    for (int y = 0; y < height; ++y)
      for (int x = 0; x < width; ++x)
      {
          pi[y][x] = 0;
          ph[y][x] = 0;
          pf[y][x] = 0;
      }
}


void
fillPixels2 (Array2D<unsigned int> &pi,
           Array2D<half> &ph,
           Array2D<float> &pf,
           int width,
           int height)
{
    cout << "pattern 1" << endl;

    for (int y = 0; y < height; ++y)
      for (int x = 0; x < width; ++x)
      {
          pi[y][x] = (x + y) & 1;
          ph[y][x] = pi[y][x];
          pf[y][x] = pi[y][x];
      }
}


void
fillPixels3 (Array2D<unsigned int> &pi,
           Array2D<half> &ph,
           Array2D<float> &pf,
           int width,
           int height)
{
    cout << "pattern 2" << endl;

    for (int y = 0; y < height; ++y)
      for (int x = 0; x < width; ++x)
      {
          pi[y][x] = x % 100 + 100 * (y % 100);
          ph[y][x] = sin (double (x)) + sin (y * 0.5);
          pf[y][x] = sin (double (y)) + sin (x * 0.5);
      }
}


void
fillPixels4 (Array2D<unsigned int> &pi,
           Array2D<half> &ph,
           Array2D<float> &pf,
           int width,
           int height)
{
    cout << "random bits" << endl;

    //
    // Use of a union to extract the bit pattern from a float, as is
    // done below, works only if int and float have the same size.
    //

    assert (sizeof (int) == sizeof (float));

    Rand48 rand;

    for (int y = 0; y < height; ++y)
      for (int x = 0; x < width; ++x)
      {
          pi[y][x] = rand.nexti();

          ph[y][x].setBits (rand.nexti());

          union {int i; float f;} u;
          u.i = rand.nexti();

          pf[y][x] = u.f;
      }
}


void
writeRead (const Array2D<unsigned int> &pi1,
         const Array2D<half> &ph1,
         const Array2D<float> &pf1,
         const char fileName[],
         int width,
         int height,
         int xOffset,
         int yOffset,
         Compression comp,
         int xs,
         int ys)
{
    //
    // Write the pixel data in pi1, ph1 and ph2 to an
    // image file using the specified compression type
    // and subsampling rates.  Read the pixel data back
    // from the file and verify that the data did not
    // change.
    //

    cout << "compression " << comp <<
          ", x sampling " << xs <<
          ", y sampling " << ys <<
          ":" << flush;


    Header hdr ((Box2i (V2i (0, 0),             // display window
                    V2i (width - 1, height -1))),
            (Box2i (V2i (xOffset, yOffset),           // data window
                    V2i (xOffset + width - 1, yOffset + height - 1))));

    hdr.compression() = comp;

    hdr.channels().insert ("I",                 // name
                     Channel (UINT, // type
                            xs,           // xSampling
                            ys)           // ySampling
                    );

    hdr.channels().insert ("H",                 // name
                     Channel (HALF, // type
                            xs,           // xSampling
                            ys)           // ySampling
                    );

    hdr.channels().insert ("F",                 // name
                     Channel (FLOAT,      // type
                            xs,           // xSampling
                            ys)           // ySampling
                    );

    {
      FrameBuffer fb; 

      fb.insert ("I",                                 // name
               Slice (UINT,                           // type
                    (char *) &pi1[-yOffset / ys][-xOffset / xs], // base
                    sizeof (pi1[0][0]),               // xStride
                    sizeof (pi1[0][0]) * (width / xs),      // yStride
                    xs,                         // xSampling
                    ys)                         // ySampling
              );
      
      fb.insert ("H",                                 // name
               Slice (HALF,                           // type
                    (char *) &ph1[-yOffset / ys][-xOffset / xs], // base
                    sizeof (ph1[0][0]),               // xStride
                    sizeof (ph1[0][0]) * (width / xs),      // yStride
                    xs,                         // xSampling
                    ys)                         // ySampling
              );
      
      fb.insert ("F",                                 // name
               Slice (FLOAT,                    // type
                    (char *) &pf1[-yOffset / ys][-xOffset / xs], // base
                    sizeof (pf1[0][0]),               // xStride
                    sizeof (pf1[0][0]) * (width / xs),      // yStride
                    xs,                         // xSampling
                    ys)                         // ySampling
              );
      
      cout << " writing" << flush;

      remove (fileName);
      OutputFile out (fileName, hdr);
      out.setFrameBuffer (fb);
      out.writePixels (height);
    }

    {
      cout << " reading" << flush;

      InputFile in (fileName);
      
      const Box2i &dw = in.header().dataWindow();
      int w = dw.max.x - dw.min.x + 1;
      int h = dw.max.y - dw.min.y + 1;
      int dx = dw.min.x;
      int dy = dw.min.y;

      Array2D<unsigned int> pi2 (h / ys, w / xs);
      Array2D<half>         ph2 (h / ys, w / xs);
      Array2D<float>        pf2 (h / ys, w / xs);

      FrameBuffer fb;

      {
          int xs = in.header().channels()["I"].xSampling;
          int ys = in.header().channels()["I"].ySampling;

          fb.insert ("I",                             // name
                   Slice (UINT,                       // type
                        (char *) &pi2[-dy / ys][-dx / xs], // base
                        sizeof (pi2[0][0]),           // xStride
                        sizeof (pi2[0][0]) * (w / xs),      // yStride
                        xs,                     // xSampling
                        ys)                     // ySampling
                  );
      }
          
      {
          int xs = in.header().channels()["H"].xSampling;
          int ys = in.header().channels()["H"].ySampling;

          fb.insert ("H",                             // name
                   Slice (HALF,                       // type
                        (char *) &ph2[-dy / ys][-dx / xs], // base
                        sizeof (ph2[0][0]),           // xStride
                        sizeof (ph2[0][0]) * (w / xs),      // yStride
                        xs,                     // xSampling
                        ys)                     // ySampling
                  );
      }
          
      {
          int xs = in.header().channels()["F"].xSampling;
          int ys = in.header().channels()["F"].ySampling;

          fb.insert ("F",                             // name
                   Slice (FLOAT,                      // type
                        (char *) &pf2[-dy / ys][-dx / xs], // base
                        sizeof (pf2[0][0]),           // xStride
                        sizeof (pf2[0][0]) * (w / xs),      // yStride
                        xs,                     // xSampling
                        ys)                     // ySampling
                  );
      }
      
      in.setFrameBuffer (fb);
      in.readPixels (dw.min.y, dw.max.y);

      cout << " comparing" << flush;

      assert (in.header().displayWindow() == hdr.displayWindow());
      assert (in.header().dataWindow() == hdr.dataWindow());
      assert (in.header().pixelAspectRatio() == hdr.pixelAspectRatio());
      assert (in.header().screenWindowCenter() == hdr.screenWindowCenter());
      assert (in.header().screenWindowWidth() == hdr.screenWindowWidth());
      assert (in.header().lineOrder() == hdr.lineOrder());
      assert (in.header().compression() == hdr.compression());

      ChannelList::ConstIterator hi = hdr.channels().begin();
      ChannelList::ConstIterator ii = in.header().channels().begin();

      while (hi != hdr.channels().end())
      {
          assert (!strcmp (hi.name(), ii.name()));
          assert (hi.channel().type == ii.channel().type);
          assert (hi.channel().xSampling == ii.channel().xSampling);
          assert (hi.channel().ySampling == ii.channel().ySampling);

          ++hi;
          ++ii;
      }

      assert (ii == in.header().channels().end());

      for (int y = 0; y < h / ys; ++y)
      {
          for (int x = 0; x < w / xs; ++x)
          {
            assert (pi1[y][x] == pi2[y][x]);
            assert (ph1[y][x].bits() == ph2[y][x].bits());
            assert (equivalent (pf1[y][x], pf2[y][x], comp));
          }
      }
    }

    remove (fileName);
    cout << endl;
}


void
writeRead (const Array2D<unsigned int> &pi,
         const Array2D<half> &ph,
         const Array2D<float> &pf,
         int w,
         int h,
         int dx,
         int dy)
{
    const char *filename = IMF_TMP_DIR "imf_test_comp.exr";

    for (int xs = 1; xs <= 2; ++xs)
    {
      for (int ys = 1; ys <= 2; ++ys)
      {
          for (int comp = 0; comp < NUM_COMPRESSION_METHODS; ++comp)
          {
            writeRead (pi, ph, pf,
                     filename,
                     w  * xs, h  * ys,
                     dx * xs, dy * ys,
                     Compression (comp),
                     xs, ys);
          }
      }
    }
}

} // namespace


void
testCompression ()
{
    try
    {
      cout << "Testing pixel data types, "
            "subsampling and "
            "compression schemes" << endl;

      const int W = 1371;
      const int H = 159;
      const int DX = 17;
      const int DY = 29;

      Array2D<unsigned int> pi (H, W);
      Array2D<half> ph (H, W);
      Array2D<float> pf (H, W);

      //
      // If the following assertion fails, new pixel types have
      // been added to the Imf library; testing code for the new
      // pixel types should be added to this file.
      //

      assert (NUM_PIXELTYPES == 3);

      fillPixels1 (pi, ph, pf, W, H);
      writeRead (pi, ph, pf, W, H, DX, DY);

      fillPixels2 (pi, ph, pf, W, H);
      writeRead (pi, ph, pf, W, H, DX, DY);

      fillPixels3 (pi, ph, pf, W, H);
      writeRead (pi, ph, pf, W, H, DX, DY);

      fillPixels4 (pi, ph, pf, W, H);
      writeRead (pi, ph, pf, W, H, DX, DY);

      cout << "ok\n" << endl;
    }
    catch (const std::exception &e)
    {
      cerr << "ERROR -- caught exception: " << e.what() << endl;
      assert (false);
    }
}

Generated by  Doxygen 1.6.0   Back to index