itkMRCImageIO.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkMRCImageIO.h,v $
  Language:  C++
  Date:      $Date: 2010-06-23 16:11:17 $
  Version:   $Revision: 1.5 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm 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 __itkMRCImageIO_h
#define __itkMRCImageIO_h

#include "itkStreamingImageIOBase.h"
#include "itkMRCHeaderObject.h"
#include <numeric>

namespace itk
{

class ITK_EXPORT MRCImageIO 
  : public StreamingImageIOBase
{
public:
  typedef MRCImageIO           Self;
  typedef StreamingImageIOBase Superclass;
  typedef SmartPointer<Self>   Pointer;

  itkNewMacro(Self);

  itkTypeMacro(MRCImageIO, StreamingImageIOBase);


  // we don't use this method
  virtual void WriteImageInformation( void ) {};

  //-------- This part of the interface deals with reading data. ------ 

  // See super class for documentation
  virtual bool CanReadFile(const char*);

  // See super class for documentation
  virtual void ReadImageInformation();

  // See super class for documentation
  virtual void Read(void* buffer);

  // -------- This part of the interfaces deals with writing data. -----

  virtual bool CanWriteFile(const char*);


  // see super class for documentation
  virtual void Write(const void* buffer);

  static const char * m_MetaDataHeaderName;

protected:
  MRCImageIO();
  // ~MRCImageIO(); // default works
  void PrintSelf(std::ostream& os, Indent indent) const;

  virtual SizeType GetHeaderSize( void ) const;

private:

  MRCImageIO(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented
  
  // internal methods to update the min and max in the header based on
  // the data, in the image buffer to be written
  template<typename TPixelType>
  void UpdateHeaderWithMinMaxMean( const TPixelType * bufferBegin)
    {
    typedef const TPixelType*  ConstPixelPointer;

    ConstPixelPointer bufferEnd = bufferBegin + m_IORegion.GetNumberOfPixels();

    // this could be replaced with std::min_element and
    // std::max_element, but that is slighlty less efficient
    std::pair<ConstPixelPointer, ConstPixelPointer> mm =
      min_max_element(bufferBegin, bufferEnd);

    double mean = std::accumulate( bufferBegin, bufferEnd, double(0.0) ) /
      std::distance( bufferBegin, bufferEnd );

    m_MRCHeader->m_Header.amin = float(*mm.first);
    m_MRCHeader->m_Header.amax = float(*mm.second);
    m_MRCHeader->m_Header.amean = float(mean);
    }

  void UpdateHeaderWithMinMaxMean( const void * bufferBegin );

  // internal methods to update the header object from the ImageIO's
  // set member variables
  void UpdateHeaderFromImageIO( void );
  
  // reimplemented
  void InternalReadImageInformation(std::ifstream& is);
 
  virtual void WriteImageInformation( const void * bufferBegin );


  MRCHeaderObject::Pointer m_MRCHeader;
};


} // namespace itk


#endif