DGtal  0.9.3beta
tutorial-examples/2DSliceImageFromVol.cpp

Example of tutorial 5: ImageAdapter: some simple image manipulations

See also
Tutorial "Image -> ImageAdapter: some simple image manipulations"
slice20Lobster.png
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)
{
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 );
Image3D::Value, functors::Identity > SliceImageAdapter;
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.beginBlock("Exporting...");
sliceImageZ >> outputFileName;
return 0;
}