DGtal  0.9.4beta
images/extract2DSlicesImagesFrom3D.cpp

Example of 2D slice images extraction from 3D images.

lobsterSliceZ_20.png
Example of resulting 2D slices images of the sample/lobster.vol file.
#include <iostream>
#include <sstream>
#include "DGtal/kernel/BasicPointFunctors.h"
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/io/readers/VolReader.h"
#include "DGtal/io/writers/PGMWriter.h"
#include "DGtal/images/ImageSelector.h"
#include "DGtal/images/ConstImageAdapter.h"
#include "ConfigExamples.h"
using namespace std;
using namespace DGtal;
int main( )
{
Image3D::Value, DGtal::functors::Identity > SliceImageAdapter;
// Importing a 3D image
std::string filename = examplesPath + "samples/lobster.vol";
Image3D image = VolReader<Image3D>::importVol( filename );
DGtal::Z2i::Domain domain(proj(image.domain().lowerBound()),
proj(image.domain().upperBound()));
trace.beginBlock ( "Example extract2DImagesFrom3D" );
// Extracting 2D slices ... and export them in the pgm format.
for (unsigned int i=0; i<30; i+=10){
std::stringstream name;
name << "lobsterSliceZ_" << i << ".pgm";
DGtal::functors::Projector<Z3i::Space> aSliceFunctor(i); aSliceFunctor.initAddOneDim(2);
SliceImageAdapter sliceImageZ(image, domain, aSliceFunctor, idV);
PGMWriter<SliceImageAdapter>::exportPGM(name.str(), sliceImageZ);
}
// trace.endBlock();
return 0;
}
// //