/*=========================================================================

  Program:   Advanced Normalization Tools
  Module:    $RCSfile: itkVectorImageFileReader.hxx,v $
  Language:  C++
  Date:      $Date: 2008/11/15 23:46:06 $
  Version:   $Revision: 1.16 $

  Copyright (c) ConsortiumOfANTS. All rights reserved.
  See accompanying COPYING.txt or
 http://sourceforge.net/projects/advants/files/ANTS/ANTSCopyright.txt for details.

     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef _itkVectorImageFileReader_hxx
#define _itkVectorImageFileReader_hxx
#include "itkVectorImageFileReader.h"

#include "itkObjectFactory.h"
#include "itkImageIOFactory.h"
#include "itkConvertPixelBuffer.h"
#include "itkImageRegion.h"
#include "itkPixelTraits.h"
#include "itkVectorImage.h"
#include "itkImageRegionIterator.h"

#include <itksys/SystemTools.hxx>
#include <fstream>

namespace itk
{
template <class TImage, class TVectorImage, class ConvertPixelTraits>
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::VectorImageFileReader()
{
  m_ImageIO = 0;
  m_FileName = "";
  m_UserSpecifiedImageIO = false;
  m_UseAvantsNamingConvention = true;
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::~VectorImageFileReader()
{
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::PrintSelf(std::ostream& os, Indent indent) const
{
  Superclass::PrintSelf(os, indent);

  if( m_ImageIO )
    {
    os << indent << "ImageIO: \n";
    m_ImageIO->Print(os, indent.GetNextIndent() );
    }
  else
    {
    os << indent << "ImageIO: (null)" << "\n";
    }

  os << indent << "UserSpecifiedImageIO flag: " << m_UserSpecifiedImageIO << "\n";
  os << indent << "m_FileName: " << m_FileName << "\n";
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::SetImageIO( ImageIOBase * imageIO)
{
  itkDebugMacro("setting ImageIO to " << imageIO );
  if( this->m_ImageIO != imageIO )
    {
    this->m_ImageIO = imageIO;
    this->Modified();
    }
  m_UserSpecifiedImageIO = true;
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::GenerateOutputInformation(void)
{
  typename TVectorImage::Pointer output = this->GetOutput();

  itkDebugMacro(<< "Reading file for GenerateOutputInformation()" << m_FileName);

  // Check to see if we can read the file given the name or prefix
  //
  if( m_FileName == "" )
    {
    throw VectorImageFileReaderException(__FILE__, __LINE__, "FileName must be specified", ITK_LOCATION);
    }

  // Test if the files exist and if it can be open.
  // and exception will be thrown otherwise.
  //
  std::string tmpFileName = this->m_FileName;

  // Test if the file exists and if it can be opened.
  // An exception will be thrown otherwise.
  // We catch the exception because some ImageIO's may not actually
  // open a file. Still reports file error if no ImageIO is loaded.

  try
    {
    m_ExceptionMessage = "";
    this->TestFileExistanceAndReadability();
    }
  catch( itk::ExceptionObject & err )
    {
    m_ExceptionMessage = err.GetDescription();
    }

  std::string::size_type pos = this->m_FileName.rfind( "." );
  std::string            extension( this->m_FileName, pos, this->m_FileName.length() - 1 );
  std::string            filename = std::string( this->m_FileName, 0, pos );

  std::string gzExtension( "" );
  if( extension == std::string( ".gz" ) )
    {
    gzExtension = extension;
    std::string::size_type pos2 = filename.rfind( "." );
    extension = std::string( filename, pos2, filename.length() - 1 );
    filename = std::string( this->m_FileName, 0, pos2 );
    }
//   unsigned int dimension = itk::GetVectorDimension
//      <VectorImagePixelType>::VectorDimension;
// Assume that the first image read contains all the information to generate
//  the output image.
  for( unsigned int i = 0; i <= 1; i++ )
    {
    this->m_FileName = filename;

    if( this->m_UseAvantsNamingConvention )
      {
      switch( i )
        {
        case 0:
          {
          this->m_FileName += std::string( "xvec" );
          }
          break;
        case 1:
          {
          this->m_FileName += std::string( "yvec" );
          }
          break;
        case 2:
          {
          this->m_FileName += std::string( "zvec" );
          }
          break;
        default:
          {
          this->m_FileName += std::string( "you_are_screwed_vec" );
          }
          break;
        }
      }
    else
      {
      std::ostringstream buf;
      buf << i;
      this->m_FileName += ( std::string( "." )  + std::string( buf.str().c_str() ) );
      }
    this->m_FileName += extension;
    if( !gzExtension.empty() )
      {
      this->m_FileName += std::string( ".gz" );
      }

    if( i == 0 )
      {
      itkDebugMacro( << "Generating output information from the file " << this->m_FileName );

      if( m_UserSpecifiedImageIO == false )  // try creating via factory
        {
        m_ImageIO = ImageIOFactory::CreateImageIO( m_FileName.c_str(), ImageIOFactory::ReadMode );
        }

      if( m_ImageIO.IsNull() )
        {
        std::ostringstream msg;
        msg << " Could not create IO object for file "
            << m_FileName.c_str() << std::endl;
        if( m_ExceptionMessage.size() )
          {
          msg << m_ExceptionMessage;
          }
        else
          {
          msg << "  Tried to create one of the following:" << std::endl;
          std::list<LightObject::Pointer> allobjects =
            ObjectFactoryBase::CreateAllInstance("itkImageIOBase");
          for( std::list<LightObject::Pointer>::iterator i = allobjects.begin();
               i != allobjects.end(); ++i )
            {
            ImageIOBase* io = dynamic_cast<ImageIOBase *>(i->GetPointer() );
            msg << "    " << io->GetNameOfClass() << std::endl;
            }
          msg << "  You probably failed to set a file suffix, or" << std::endl;
          msg << "    set the suffix to an unsupported type." << std::endl;
          }
        VectorImageFileReaderException e(__FILE__, __LINE__, msg.str().c_str(), ITK_LOCATION);
        throw e;
        return;
        }

      // Got to allocate space for the image. Determine the characteristics of
      // the image.
      //
      m_ImageIO->SetFileName(m_FileName.c_str() );
      m_ImageIO->ReadImageInformation();

      typename TVectorImage::SizeType dimSize;
      double spacing[TVectorImage::ImageDimension];
      double origin[TVectorImage::ImageDimension];
      typename TVectorImage::DirectionType direction;
      std::vector<double> axis;
      for( unsigned int k = 0; k < TImage::ImageDimension; k++ )
        {
        if( k < m_ImageIO->GetNumberOfDimensions() )
          {
          dimSize[k] = m_ImageIO->GetDimensions(k);
          spacing[k] = m_ImageIO->GetSpacing(k);
          origin[k]  = m_ImageIO->GetOrigin(k);
          // Please note: direction cosines are stored as columns of the
          // direction matrix
          axis = m_ImageIO->GetDirection(k);
          for( unsigned j = 0; j < TImage::ImageDimension; j++ )
            {
            if( j < m_ImageIO->GetNumberOfDimensions() )
              {
              direction[j][k] = axis[j];
              }
            else
              {
              direction[j][k] = 0.0;
              }
            }
          }
        else
          {
          // Number of dimensions in the output is more than number of dimensions
          // in the ImageIO object (the file).  Use default values for the size,
          // spacing, origin and direction for the final (degenerate) dimensions.
          dimSize[i] = 1;
          spacing[i] = 1.0;
          origin[i] = 0.0;
          for( unsigned j = 0; j < TImage::ImageDimension; j++ )
            {
            if( i == j )
              {
              direction[j][k] = 1.0;
              }
            else
              {
              direction[j][k] = 0.0;
              }
            }
          }
        }

      output->SetSpacing( spacing );     // Set the image spacing
      output->SetOrigin( origin );       // Set the image origin
      output->SetDirection( direction ); // Set the image direction cosines

      // Copy MetaDataDictionary from instantiated reader to output image.
      output->SetMetaDataDictionary(m_ImageIO->GetMetaDataDictionary() );
      this->SetMetaDataDictionary(m_ImageIO->GetMetaDataDictionary() );

      this->m_Image->SetSpacing( spacing );
      this->m_Image->SetOrigin( origin );
      this->m_Image->SetDirection( direction );
      this->m_Image->SetMetaDataDictionary(m_ImageIO->GetMetaDataDictionary() );

      typedef typename TVectorImage::IndexType IndexType;

      IndexType start;
      start.Fill(0);

      VectorImageRegionType region;
      region.SetSize(dimSize);
      region.SetIndex(start);

      ImageRegionType imageregion;
      imageregion.SetSize(dimSize);
      imageregion.SetIndex(start);

      // If a VectorImage, this requires us to set the
      // VectorLength before allocate
      // if( strcmp( output->GetNameOfClass(), "VectorImage" ) == 0 )
      //  {
      //  typedef typename TImage::AccessorFunctorType AccessorFunctorType;
      //  AccessorFunctorType::SetVectorLength( output, m_ImageIO->GetNumberOfComponents() );
      //  }

      output->SetLargestPossibleRegion( region );
      this->m_Image->SetLargestPossibleRegion( imageregion );
      }
    }
  this->m_FileName = tmpFileName;
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::TestFileExistanceAndReadability()
{
  std::string tmpFileName = this->m_FileName;

  std::string::size_type pos = this->m_FileName.rfind( "." );
  std::string            extension( this->m_FileName, pos, this->m_FileName.length() - 1 );
  std::string            filename = std::string( this->m_FileName, 0, pos );

  std::string gzExtension( "" );

  if( extension == std::string( ".gz" ) )
    {
    gzExtension = extension;
    std::string::size_type pos2 = filename.rfind( "." );
    extension = std::string( filename, pos2, filename.length() - 1 );
    filename = std::string( this->m_FileName, 0, pos2 );
    }

  unsigned int dimension = itk::GetVectorDimension
    <VectorImagePixelType>::VectorDimension;
  for( unsigned int i = 0; i < dimension; i++ )
    {
    this->m_FileName = filename;

    if( this->m_UseAvantsNamingConvention )
      {
      switch( i )
        {
        case 0:
          {
          this->m_FileName += std::string( "xvec" );
          }
          break;
        case 1:
          {
          this->m_FileName += std::string( "yvec" );
          }
          break;
        case 2:
          {
          this->m_FileName += std::string( "zvec" );
          }
          break;
        default:
          {
          this->m_FileName += std::string( "you_are_screwed_vec" );
          }
          break;
        }
      }
    else
      {
      std::ostringstream buf;
      buf << i;
      this->m_FileName += ( std::string( "." )  + std::string( buf.str().c_str() ) );
      }
    this->m_FileName += extension;
    if( !gzExtension.empty() )
      {
      this->m_FileName += std::string( ".gz" );
      }

    itkDebugMacro( << "Checking for the file " << this->m_FileName );

    // Test if the file exists.
    if( !itksys::SystemTools::FileExists( m_FileName.c_str() ) )
      {
      VectorImageFileReaderException e(__FILE__, __LINE__);
      std::ostringstream             msg;
      msg << "The file doesn't exists. "
          << std::endl << "Filename = " << m_FileName
          << std::endl;
      e.SetDescription(msg.str().c_str() );
      throw e;
      return;
      }

    // Test if the file can be open for reading access.
    std::ifstream readTester;
    readTester.open( m_FileName.c_str() );
    if( readTester.fail() )
      {
      readTester.close();
      std::ostringstream msg;
      msg << "The file couldn't be opened for reading. "
          << std::endl << "Filename: " << m_FileName
          << std::endl;
      VectorImageFileReaderException e(__FILE__, __LINE__, msg.str().c_str(), ITK_LOCATION);
      throw e;
      return;
      }
    readTester.close();
    }
  this->m_FileName = tmpFileName;
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::EnlargeOutputRequestedRegion(DataObject *output)
{
  typename TVectorImage::Pointer out = dynamic_cast<TVectorImage *>(output);

  // the ImageIO object cannot stream, then set the RequestedRegion to the
  // LargestPossibleRegion
  if( !m_ImageIO->CanStreamRead() )
    {
    if( out )
      {
      out->SetRequestedRegion( out->GetLargestPossibleRegion() );
      }
    else
      {
      throw VectorImageFileReaderException(__FILE__, __LINE__,
                                           "Invalid output object type");
      }
    }
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::GenerateData()
{
  typename TVectorImage::Pointer output = this->GetOutput();

  // allocate the output buffer
  output->SetBufferedRegion( output->GetRequestedRegion() );
  output->Allocate();

  this->m_Image->SetBufferedRegion( output->GetRequestedRegion() );
  this->m_Image->Allocate();

  // Test if the file exist and if it can be open.
  // and exception will be thrown otherwise.
  try
    {
    m_ExceptionMessage = "";
    this->TestFileExistanceAndReadability();
    }
  catch( itk::ExceptionObject & err )
    {
    m_ExceptionMessage = err.GetDescription();
    }

  std::string tmpFileName = this->m_FileName;

  std::string::size_type pos = this->m_FileName.rfind( "." );
  std::string            extension( this->m_FileName, pos, this->m_FileName.length() - 1 );
  std::string            filename = std::string( this->m_FileName, 0, pos );

  std::string gzExtension( "" );
  if( extension == std::string( ".gz" ) )
    {
    gzExtension = extension;
    std::string::size_type pos2 = filename.rfind( "." );
    extension = std::string( filename, pos2, filename.length() - 1 );
    filename = std::string( this->m_FileName, 0, pos2 );
    }

  unsigned int dimension = itk::GetVectorDimension
    <VectorImagePixelType>::VectorDimension;
  for( unsigned int i = 0; i < dimension; i++ )
    {
    this->m_FileName = filename;

    if( this->m_UseAvantsNamingConvention )
      {
      switch( i )
        {
        case 0:
          {
          this->m_FileName += std::string( "xvec" );
          }
          break;
        case 1:
          {
          this->m_FileName += std::string( "yvec" );
          }
          break;
        case 2:
          {
          this->m_FileName += std::string( "zvec" );
          }
          break;
        default:
          {
          this->m_FileName += std::string( "you_are_screwed_vec" );
          }
          break;
        }
      }
    else
      {
      std::ostringstream buf;
      buf << i;
      this->m_FileName += ( std::string( "." )  + std::string( buf.str().c_str() ) );
      }
    this->m_FileName += extension;
    if( !gzExtension.empty() )
      {
      this->m_FileName += std::string( ".gz" );
      }

    itkDebugMacro( << "Reading image buffer from the file " << this->m_FileName );

    // Tell the ImageIO to read the file
    m_ImageIO->SetFileName(m_FileName.c_str() );

    this->m_FileName = tmpFileName;
    ImageIORegion ioRegion(TImage::ImageDimension);

    ImageIORegion::SizeType  ioSize = ioRegion.GetSize();
    ImageIORegion::IndexType ioStart = ioRegion.GetIndex();

    typename TImage::SizeType dimSize;
    for( unsigned int j = 0; j < TImage::ImageDimension; j++ )
      {
      if( j < m_ImageIO->GetNumberOfDimensions() )
        {
        dimSize[j] = m_ImageIO->GetDimensions(j);
        }
      else
        {
        // Number of dimensions in the output is more than number of dimensions
        // in the ImageIO object (the file).  Use default values for the size,
        // spacing, and origin for the final (degenerate) dimensions.
        dimSize[j] = 1;
        }
      }
    for( unsigned int j = 0; j < dimSize.GetSizeDimension(); ++j )
      {
      ioSize[j] = dimSize[j];
      }

    typedef typename TImage::IndexType IndexType;
    IndexType start;
    start.Fill(0);
    for( unsigned int j = 0; j < start.GetIndexDimension(); ++j )
      {
      ioStart[j] = start[j];
      }

    ioRegion.SetSize(ioSize);
    ioRegion.SetIndex(ioStart);

    itkDebugMacro(<< "ioRegion: " << ioRegion);

    m_ImageIO->SetIORegion(ioRegion);

    char *loadBuffer = 0;
    // the size of the buffer is computed based on the actual number of
    // pixels to be read and the actual size of the pixels to be read
    // (as opposed to the sizes of the output)
    try
      {
      if(  /** FIXME */
// m_ImageIO->GetComponentTypeInfo()
//           != typeid(ITK_TYPENAME ConvertPixelTraits::ComponentType)
//           ||
        (m_ImageIO->GetNumberOfComponents()
         != ConvertPixelTraits::GetNumberOfComponents() ) )
        {
        // the pixel types don't match so a type conversion needs to be
        // performed
        itkDebugMacro(<< "Buffer conversion required from: "
                      <<  " FIXME m_ImageIO->GetComponentTypeInfo().name() "
                      << " to: "
                      << typeid(ITK_TYPENAME ConvertPixelTraits::ComponentType).name() );

        loadBuffer = new char[m_ImageIO->GetImageSizeInBytes()];
        m_ImageIO->Read( static_cast<void *>(loadBuffer) );

        // See note below as to why the buffered region is needed and
        // not actualIOregion
        this->DoConvertBuffer(static_cast<void *>(loadBuffer),
                              output->GetBufferedRegion().GetNumberOfPixels() );
        }
      else // a type conversion is not necessary
        {
        itkDebugMacro(<< "No buffer conversion required.");

        ImagePixelType *buffer =
          this->m_Image->GetPixelContainer()->GetBufferPointer();
        m_ImageIO->Read( buffer );
        }
      }
    catch( ... )
      {
      // if an exception is thrown catch it

      if( loadBuffer )
        {
        // clean up
        delete [] loadBuffer;
        loadBuffer = 0;
        }

      // then rethrow
      throw;
      }

    // clean up
    if( loadBuffer )
      {
      delete [] loadBuffer;
      loadBuffer = 0;
      }

    ImageRegionIterator<TVectorImage> Id( this->GetOutput(),
                                          this->GetOutput()->GetLargestPossibleRegion() );
    ImageRegionIterator<TImage> It( this->m_Image,
                                    this->m_Image->GetLargestPossibleRegion() );

    Id.GoToBegin();
    It.GoToBegin();
    while( !Id.IsAtEnd() || !It.IsAtEnd() )
      {
      VectorImagePixelType V = Id.Get();
      V[i] = static_cast<typename VectorImagePixelType::ValueType>( It.Get() );
      Id.Set( V );
      ++Id;
      ++It;
      }
    }
}

template <class TImage, class TVectorImage, class ConvertPixelTraits>
void
VectorImageFileReader<TImage, TVectorImage, ConvertPixelTraits>
::DoConvertBuffer(void* inputData,
                  unsigned long numberOfPixels)
{
  // get the pointer to the destination buffer
  ImagePixelType *imageData =
    this->m_Image->GetPixelContainer()->GetBufferPointer();

  // TODO:
  // Pass down the PixelType (RGB, VECTOR, etc.) so that any vector to
  // scalar conversion be type specific. i.e. RGB to scalar would use
  // a formula to convert to luminance, VECTOR to scalar would use
  // vector magnitude.

// Create a macro as this code is a bit lengthy and repetitive
// if the ImageIO pixel type is typeid(type) then use the ConvertPixelBuffer
// class to convert the data block to TImage's pixel type
// see DefaultConvertPixelTraits and ConvertPixelBuffer

// The first else if block applies only to images of type itk::VectorImage
// VectorImage needs to copy out the buffer differently.. The buffer is of
// type InternalPixelType, but each pixel is really 'k' consecutive pixels.

#define ITK_CONVERT_BUFFER_IF_BLOCK(type)               \
  else if( true /** FIXME  m_ImageIO->GetComponentTypeInfo() == typeid(type) )  */ )   \
    {                                                   \
    if( strcmp( this->GetOutput()->GetNameOfClass(), "VectorImage" ) == 0 ) \
      { \
      ConvertPixelBuffer<                                 \
        type,                                             \
        ImagePixelType,                             \
        ConvertPixelTraits                                \
        >                                                 \
      ::ConvertVectorImage(                             \
        static_cast<type *>(inputData),                  \
        m_ImageIO->GetNumberOfComponents(),             \
        imageData,                                     \
        numberOfPixels);                                \
      } \
    else \
      { \
      ConvertPixelBuffer<                                 \
        type,                                             \
        ImagePixelType,                             \
        ConvertPixelTraits                                \
        >                                                 \
      ::Convert(                                        \
        static_cast<type *>(inputData),                  \
        m_ImageIO->GetNumberOfComponents(),             \
        imageData,                                     \
        numberOfPixels);                                \
      } \
    }

  if( 0 )
    {
    }
  ITK_CONVERT_BUFFER_IF_BLOCK(unsigned char)
  ITK_CONVERT_BUFFER_IF_BLOCK(char)
  ITK_CONVERT_BUFFER_IF_BLOCK(unsigned short)
  ITK_CONVERT_BUFFER_IF_BLOCK( short)
  ITK_CONVERT_BUFFER_IF_BLOCK(unsigned int)
  ITK_CONVERT_BUFFER_IF_BLOCK( int)
  ITK_CONVERT_BUFFER_IF_BLOCK(unsigned long)
  ITK_CONVERT_BUFFER_IF_BLOCK( long)
  ITK_CONVERT_BUFFER_IF_BLOCK(float)
  ITK_CONVERT_BUFFER_IF_BLOCK( double)
  else
    {
    VectorImageFileReaderException e(__FILE__, __LINE__);
    std::ostringstream             msg;
    msg << "Couldn't convert component type: "
        << std::endl << "    "
        << m_ImageIO->GetComponentTypeAsString(m_ImageIO->GetComponentType() )
        << std::endl << "to one of: "
        << std::endl << "    " << typeid(unsigned char).name()
        << std::endl << "    " << typeid(char).name()
        << std::endl << "    " << typeid(unsigned short).name()
        << std::endl << "    " << typeid(short).name()
        << std::endl << "    " << typeid(unsigned int).name()
        << std::endl << "    " << typeid(int).name()
        << std::endl << "    " << typeid(unsigned long).name()
        << std::endl << "    " << typeid(long).name()
        << std::endl << "    " << typeid(float).name()
        << std::endl << "    " << typeid(double).name()
        << std::endl;
    e.SetDescription(msg.str().c_str() );
    e.SetLocation(ITK_LOCATION);
    throw e;
    return;
    }
#undef ITK_CONVERT_BUFFER_IF_BLOCK
}
} // namespace ITK

#endif