VolReader.h

#pragma once

#if defined(VolReader_RECURSES)
#error Recursive header files inclusion detected in VolReader.h
#else // defined(VolReader_RECURSES)

#define VolReader_RECURSES

#if !defined VolReader_h

#define VolReader_h

// Inclusions
#include <iostream>
#include <sstream>
#include <string>
#include <cstdio>
#include "DGtal/base/Common.h"
#include "DGtal/base/CUnaryFunctor.h"

namespace DGtal
{

  // template class VolReader
  template <typename TImageContainer,  
            typename TFunctor = functors::Cast< typename TImageContainer::Value > >
  struct VolReader
  {
    // ----------------------- Standard services ------------------------------

    typedef TImageContainer ImageContainer;
    typedef typename TImageContainer::Value Value;
    typedef TFunctor Functor;
    
    BOOST_CONCEPT_ASSERT((  concepts::CUnaryFunctor<TFunctor, unsigned char, Value > )) ;    

    BOOST_STATIC_ASSERT(ImageContainer::Domain::dimension == 3);


    static ImageContainer importVol(const std::string & filename, 
                                    const Functor & aFunctor =  Functor()) throw(DGtal::IOException);
    
  private:

    typedef unsigned char voxel;
    /* In recent C++, we should use a std::map, but we prefer (badly) code it
       by hand for compatibility with old compilers.
       At this time, there is a limit of 30 fields in header :-} */
    struct HeaderField {
      HeaderField( const char *t, const char *v ) :
        type( strdup(t) ), value( strdup(v) ) {}
      ~HeaderField() {
        free( type );
        free( value );
      }
      HeaderField( const HeaderField &h ) :
        type( strdup(h.type) ), value( strdup(h.value) ) {};
      HeaderField() : type(NULL), value(NULL) {};
      const HeaderField &operator = (const HeaderField &h) {
        free( type );
        free( value );
        if (h.type != NULL) {
          type = strdup( h.type );
          value = strdup( h.value );
        }
        return *this;
      }
      char *type;
      char *value;
    };


    static const char *getHeaderValue( const char *type, const HeaderField * header );

    static     int getHeaderValueAsInt( const char *type, int *dest , const HeaderField * header); 
    
    static int getHeaderField( const char *type, const HeaderField * header ) ;
    
    static const char *requiredHeaders[];
   
  }; // end of class VolReader


} // namespace DGtal


// Includes inline functions.
#include "DGtal/io/readers/VolReader.ih"

//                                                                           //

#endif // !defined VolReader_h

#undef VolReader_RECURSES
#endif // else defined(VolReader_RECURSES)