LongvolReader.h

 
#pragma once
 
#if defined(LongvolReader_RECURSES)
#error Recursive header files inclusion detected in LongvolReader.h
#else // defined(LongvolReader_RECURSES)
#define LongvolReader_RECURSES
 
#if !defined LongvolReader_h
#define LongvolReader_h
 
// Inclusions
#include <iostream>
#include <sstream>
#include <string>
#include <cstdio>
#include "DGtal/base/Common.h"
#include <boost/static_assert.hpp>
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/base/CUnaryFunctor.h"
 
 
namespace DGtal
{
  
  // template class LongvolReader
  template <typename TImageContainer,
  typename TFunctor= functors::Cast< typename TImageContainer::Value > >
  struct LongvolReader
  {
    // ----------------------- Standard services ------------------------------
    
    typedef TImageContainer ImageContainer;
    typedef typename TImageContainer::Value Value;
    typedef TFunctor Functor;
    
    BOOST_CONCEPT_ASSERT((  concepts::CUnaryFunctor<TFunctor, DGtal::uint64_t, Value > )) ;
    BOOST_STATIC_ASSERT(ImageContainer::Domain::dimension == 3);
    
    
    static ImageContainer importLongvol(const std::string & filename,
                                        const Functor & aFunctor =  Functor());
    
    
    
  private:
    
    template <typename Word>
    static
    std::stringstream& read_word( std::stringstream & fin, Word& aValue )
    {
      aValue = 0;
      char c;
      for (unsigned int size = 0; size < sizeof( Word ); ++size)
      {
        fin.get( c ) ;
        unsigned char cc=static_cast<unsigned char>(c);
        aValue |= (cc << (8 * size));
      }
      return fin;
    }
    
    
    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 LongvolReader
  
  
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/io/readers/LongvolReader.ih"
 
//                                                                           //
 
#endif // !defined LongvolReader_h
 
#undef LongvolReader_RECURSES
#endif // else defined(LongvolReader_RECURSES)