MagickReader.ih

#include <cstdlib>

// Interface - public :

template <typename TImageContainer>
inline
TImageContainer 
DGtal::MagickReader<TImageContainer>::importImage(const std::string & filename, bool topbotomOrder ) throw(DGtal::IOException)
{
  int w,h;
  const Magick::PixelPacket *cacheRead;
  DGtal::IOException dgtalio;

  Magick::InitializeMagick(NULL);
  
  Magick::Image img;

  try 
    {
      img.read ( filename );
    }
  catch( ... )
    {
      std::cout << "MagickReader : can't open " << filename << endl;
      throw dgtalio;
    }
  
  w = img.size().width();
  h = img.size().height();

  //@todo Check Magick::Types according to ValueType
  img.type ( Magick::TrueColorType );
  img.modifyImage();
  cacheRead = img.getConstPixels ( 0, 0, w, h );
  
  //@todo check ImageMagick errors.
  //trace.error() << "MagickReader : can't open "<< filename<<endl;
  
  typename TImageContainer::Point firstPoint;
  typename TImageContainer::Point lastPoint;
      
  firstPoint = TImageContainer::Point::zero;
  lastPoint[0] = w-1;
  lastPoint[1] = h-1;

  typename TImageContainer::Domain domain(firstPoint,lastPoint);
  TImageContainer image(domain);
  typename TImageContainer::Value val;

  //We scan the file
  typename TImageContainer::Domain::ConstIterator it = domain.begin(),
    itend=domain.end();
  
  for(; it != itend; ++it)
    {
      const Magick::PixelPacket *pixel = cacheRead + w * ((topbotomOrder)? ( h - 1 - (*it)[1] ): (*it)[1] ) + (*it)[0];
      val = (pixel->red + pixel->green + pixel->blue) % 256;
      // std::cout<<" valread= "<< (int)val << std::endl;
      image.setValue( (*it), val  );
    }

  return image;
}