ImageRGBLookupReader.hpp

/* ***** BEGIN LICENSE BLOCK *****
 * FW4SPL - Copyright (C) IRCAD, 2009-2018.
 * Distributed under the terms of the GNU Lesser General Public License (LGPL) as
 * published by the Free Software Foundation.
 * ****** END LICENSE BLOCK ****** */

#pragma once

#include "fwDcmtkIO/config.hpp"

#include <fwMedData/DicomSeries.hpp>

#include <dcmtk/config/osconfig.h>
#include <dcmtk/dcmdata/dcdeftag.h>
#include <dcmtk/dcmdata/dcfilefo.h>
#include <dcmtk/dcmimgle/dcmimage.h>
#include <dcmtk/dcmnet/diutil.h>

namespace fwDcmtkIO
{
namespace reader
{
namespace rgblookup
{

class FWDCMTKIO_CLASS_API ImageRGBLookupReader
{
public:
    typedef ::fwMedData::DicomSeries::DicomContainerType DicomContainerType;

    template< typename T, typename U >
    static void fillImageBuffer(unsigned int rows, unsigned int columns,
                                unsigned int depth,
                                DicomContainerType& instances,
                                void* destination, const T* redLookup,
                                const T* greenLookup,
                                const T* blueLookup)
    {
        U* tempoBuffer = ::fwDcmtkIO::reader::main::ImageReader::createTemporaryBuffer<U>(
            rows, columns, depth, instances);

        T* arrayBuffer = static_cast< T* >(destination);
        unsigned short x, y, z;
        unsigned int frameSize = columns * rows;
        x = y = z = 0;
        for (x = 0; x < columns; ++x)
        {
            for (y = 0; y < rows; ++y)
            {
                unsigned int yshift = x * rows;

                for (z = 0; z < depth; ++z)
                {
                    unsigned int position = y + yshift + z * frameSize;
                    U value               = tempoBuffer[position];

                    arrayBuffer[position*3]   = static_cast< T >((redLookup[value]/(double)0xffff)*256);
                    arrayBuffer[position*3+1] = static_cast< T >((greenLookup[value]/(double)0xffff)*256);
                    arrayBuffer[position*3+2] = static_cast< T >((blueLookup[value]/(double)0xffff)*256);
                }
            }
        }
        delete tempoBuffer;
    }
};

} //rgblookup
} //reader
} //fwDcmtkIO