DGtal  0.9.3
viewer3D-8-2DSliceImages.cpp
Go to the documentation of this file.
1 
36 #include <iostream>
38 #include "DGtal/base/Common.h"
39 #include "DGtal/io/readers/VolReader.h"
40 #include "DGtal/images/ImageHelper.h"
41 #include "ConfigExamples.h"
42 #include "DGtal/io/viewers/Viewer3D.h"
43 
45 #include "DGtal/io/DrawWithDisplay3DModifier.h"
46 #include "DGtal/io/viewers/DrawWithViewer3DModifier.h"
47 #include "DGtal/io/colormaps/HueShadeColorMap.h"
48 #include "DGtal/io/Color.h"
50 
52 
53 using namespace std;
54 using namespace DGtal;
55 
57 struct hueFct{
58  inline
59  unsigned int operator() (unsigned char aVal) const
60  {
61  HueShadeColorMap<unsigned char> hueShade(0,255);
62  Color col = hueShade((unsigned char)aVal);
63  return (((unsigned int) col.red()) << 16)| (((unsigned int) col.green()) << 8)|((unsigned int) col.blue());
64  }
65 };
68 
69 int main( int argc, char** argv )
70 {
71 
73  QApplication application(argc,argv);
74  typedef Viewer3D<> MyViewer;
75  MyViewer viewer;
76  viewer.show();
77  std::string inputFilename = examplesPath + "samples/lobster.vol";
78  Image3D imageVol = VolReader<Image3D>::importVol(inputFilename);
79 
80 
82  // Extracting the 2D images from the 3D one and from a given dimension.
83  // First image the teenth Z slice (dim=2)
85  Image3D::Value, DGtal::functors::Identity > MySliceImageAdapter;
86 
87  // Define the functor to recover a 2D domain from the 3D one in the Z direction (2):
88  DGtal::functors::Projector<DGtal::Z2i::Space> transTo2DdomainFunctorZ; transTo2DdomainFunctorZ.initRemoveOneDim(2);
89  DGtal::Z2i::Domain domain2DZ(transTo2DdomainFunctorZ(imageVol.domain().lowerBound()),
90  transTo2DdomainFunctorZ(imageVol.domain().upperBound()));
91 
92  // Define the functor to associate 2D coordinates to the 3D one by giving the direction Z (2) and the slide numnber (10):
93  DGtal::functors::Projector<DGtal::Z3i::Space> aSliceFunctorZ(10); aSliceFunctorZ.initAddOneDim(2);
94 
95  // We can now obtain the slice image (a ConstImageAdapter):
96  const auto identityFunctor = DGtal::functors::Identity();
97  MySliceImageAdapter aSliceImageZ(imageVol, domain2DZ, aSliceFunctorZ, identityFunctor );
98 
99  // Second image the fiftieth Y slice (dim=1)
100  // Define the functor to recover a 2D domain from the 3D one in the Y direction (1):
101  DGtal::functors::Projector<DGtal::Z2i::Space> transTo2DdomainFunctorY; transTo2DdomainFunctorY.initRemoveOneDim(1);
102  DGtal::Z2i::Domain domain2DY(transTo2DdomainFunctorY(imageVol.domain().lowerBound()),
103  transTo2DdomainFunctorY(imageVol.domain().upperBound()));
104 
105  // Define the functor to associate 2D coordinates to the 3D one by giving the direction Y (1) and the slide numnber (50):
106  DGtal::functors::Projector<DGtal::Z3i::Space> aSliceFunctorY(50); aSliceFunctorY.initAddOneDim(1);
107 
108  // We can now obtain the slice image (a ConstImageAdapter):
109  MySliceImageAdapter aSliceImageY(imageVol, domain2DY, aSliceFunctorY, identityFunctor );
111 
113  viewer << SetMode3D(aSliceImageZ.className(), "BoundingBox");
114  viewer << MyViewer::updateDisplay;
116 
118  viewer << aSliceImageZ;
119  viewer << aSliceImageY;
121 
122  viewer << SetMode3D(aSliceImageZ.className(), "");
124  viewer << AddTextureImage2DWithFunctor<MySliceImageAdapter, hueFct, Z3i::Space, Z3i::KSpace> (aSliceImageZ, hueFct(), Viewer3D<Z3i::Space, Z3i::KSpace>::RGBMode);
125  viewer << AddTextureImage2DWithFunctor<MySliceImageAdapter, hueFct, Z3i::Space, Z3i::KSpace> (aSliceImageY, hueFct(), Viewer3D<Z3i::Space, Z3i::KSpace>::RGBMode);
127 
128 
130  viewer << DGtal::UpdateImagePosition<Z3i::Space, Z3i::KSpace>(1, MyViewer::yDirection, 0.0, 50.0, 0.0);
131  viewer << DGtal::UpdateImageData<MySliceImageAdapter>(0, aSliceImageZ, 0, 0, 10);
132  viewer << MyViewer::updateDisplay;
134 
135 
137  viewer << DGtal::UpdateImagePosition<Z3i::Space, Z3i::KSpace>(3, MyViewer::yDirection, 500.0, 50.0, 0.0);
138  viewer << DGtal::UpdateImageData<MySliceImageAdapter, hueFct>(2, aSliceImageZ, 500, 0, 10, 0.0, MyViewer::zDirection, hueFct());
139  viewer << MyViewer::updateDisplay;
141 
142 return application.exec();
143 
144  trace.endBlock();
145  return 0;
146 }
147 // //
Trace trace
Definition: Common.h:137
virtual void show()
Overload QWidget method in order to add a call to updateList() method (to ensure that the lists are w...
STL namespace.
double endBlock()
Aim: This class template may be used to (linearly) convert scalar values in a given range into a colo...
void initRemoveOneDim(const Dimension &dimRemoved)
Aim: Functor that maps a point P of dimension i to a point Q of dimension j. The member myDims is an ...
Aim: Define a simple default functor that just returns its argument.
Aim: implements methods to read a "Vol" file format.
Definition: VolReader.h:89
DGtal is the top-level namespace which contains all DGtal functions and types.
void initAddOneDim(const Dimension &newDim)
Structure representing an RGB triple with alpha component.
Definition: Color.h:66
void green(const unsigned char aGreenValue)
Aim: implements a const image adapter with a given domain (i.e. a subdomain) and 2 functors : g for d...
void red(const unsigned char aRedValue)
Modifier class in a Display3D stream. Useful to choose your own mode for a given class. Realizes the concept CDrawableWithDisplay3D.
int main(int argc, char **argv)
[ExampleViewer3D2DImagesExtractImagesColorFct]
void blue(const unsigned char aBlueValue)