MatLabWriter.hpp

#ifndef MATLABWRITER_HPP
#define MATLABWRITER_HPP

#include <iostream>
#include "tiffio.h"
#include "map_grid.hpp"
#include <string>

template<unsigned int dim, typename T>
class MatLabWriter
{
    template<typename grid, typename Mem> int operator<<(grid data, std::string file)
    {
        // Grid can have several properties we can save only scalar fields
        // for each properties save one scalar fields

        // Tiff files can be saved up to 5D

        if (dim > 5)
        {
            std::cerr << "Error Tiff writer support until 5D images" << "\n";
        }

        // Open the tiff image

        uint32 width;
        uint32 height;
        TIFF *tif = TIFFOpen(file.c_str(),"w");

        // if the file is open

        if(tif)
        {
            // set width and height for 2D

            width = data.getGrid().size(0);
            height = data.getGrid().size(1);

            TIFFSetField(tif,TIFFTAG_IMAGEWIDTH, &width);
            TIFFSetField(tif,TIFFTAG_IMAGELENGTH, &height);

            // Create the tiff line, in case the grid is CPU, we have only
            // one property and is a scalar, we can directly copy the line

            typename boost::fusion::result_of::at<T::type,0>::type first_element_type;

            if (typeid(grid).name() == "grid_cpu" && T::num_prop == 1 && boost::is_array<first_element_type>::type::value == true)
            {
                // Get the grid key iterator

                grid_key_dx_iterator<dim> key = data.getIterator();

                // write all lines

                for(int i = 0; i < height ; i++)
                {
                    // select the correct lines

                    key.set(1,i);

                    // we have only one scalar properties, get the buffer pointer
                    void * buf = &data.template get<0>(key);

                    TIFFWriteScanline(tif,buf,i, 0);
                }
            }
            else
            {
                // we have to create the a scan line for each properties and index array
                // each property and index array became a channel

                // Get the grid key iterator

                grid_key_dx_iterator<dim> key = data.getIterator();

                // count how many properties and how many indexes we have

                const int n_prp = total_prop<T>;

                // write all lines

                for(int i = 0; i < height ; i++)
                {
                    // select the correct lines

                    key.set(1,i);

                    // we have only one scalar properties, get the buffer pointer
                    void * buf = &data.template get<0>(key);

                    TIFFWriteScanline(tif,buf,i, 0);
                }
            }
        }
    }

};

#endif