tutorial-examples/2DSliceImageFromVol.cpp

Example of tutorial 5: ImageAdapter: some simple image manipulations

See also

Tutorial "Image -> ImageAdapter: some simple image manipulations"

Result of slice image number 20 on Z axis.

#include <iostream>
#include <DGtal/base/Common.h>
#include "DGtal/io/readers/GenericReader.h"
#include "DGtal/io/writers/GenericWriter.h"

#include <DGtal/helpers/StdDefs.h>
#include <DGtal/images/ImageContainerBySTLVector.h>
#include <DGtal/images/ConstImageAdapter.h>


#include "DGtal/kernel/BasicPointFunctors.h"

using namespace std;
using namespace DGtal;
using namespace Z2i;



int main(int argc, char ** argv)
{
    
  typedef ImageContainerBySTLVector<Z2i::Domain, unsigned char>  Image2D;
  typedef ImageContainerBySTLVector<Z3i::Domain, unsigned char>  Image3D;

  if(argc < 3)
    {
      trace.error() << "You need to indicate the volumetric image name and slice number as parameters." << std::endl;
      trace.error() << std::endl;
      return 1;
    }
  std::string filename(argv[1]);  
  std::string outputFileName = "sliceImage.pgm";

  unsigned int numSlice = atoi(argv[2]);

  trace.beginBlock("Loading file");
  Image3D  image3d = GenericReader< Image3D >::import ( filename );
  
   
  functors::Identity id;  
  typedef ConstImageAdapter<Image3D, Image2D::Domain, functors::Projector<Z3i::Space>,
                            Image3D::Value,  functors::Identity >  SliceImageAdapter;
  functors::Projector<Z2i::Space >  proj(2);
  Z2i::Domain domain2D(proj(image3d.domain().lowerBound()),
                              proj(image3d.domain().upperBound()));

  DGtal::functors::Projector<Z3i::Space> aSliceFunctor(numSlice); 
  aSliceFunctor.initAddOneDim(2);
  SliceImageAdapter sliceImageZ(image3d, domain2D, aSliceFunctor, id);
  
  trace.endBlock();

  trace.beginBlock("Exporting...");
  sliceImageZ >> outputFileName; 
  trace.endBlock();
  return 0;
}