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 )
    {
      char* raw = reinterpret_cast<char*>(&aValue);
      for (std::size_t i = 0; i < sizeof(Word); ++i)
          raw[i] = fin.get();
      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)