Loading [MathJax]/extensions/TeX/AMSsymbols.js
DGtal 2.0.0
toricdomainvolumetric.cpp File Reference
DGtal Examples
#include <iostream>
#include <iomanip>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/io/colormaps/GrayscaleColorMap.h"
#include "DGtal/io/colormaps/HueShadeColorMap.h"
#include "DGtal/io/colormaps/TickedColorMap.h"
#include "DGtal/io/colormaps/GradientColorMap.h"
#include "DGtal/io/boards/Board2D.h"
#include "DGtal/images/ImageSelector.h"
#include "DGtal/images/SimpleThresholdForegroundPredicate.h"
#include "DGtal/geometry/volumes/distance/DistanceTransformation.h"
#include <boost/algorithm/minmax_element.hpp>
Include dependency graph for toricdomainvolumetric.cpp:

Go to the source code of this file.

Functions

int main ()

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author
David Coeurjolly (david.nosp@m..coe.nosp@m.urjol.nosp@m.ly@l.nosp@m.iris..nosp@m.cnrs.nosp@m..fr ) Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), CNRS, France
Date
2016/10/17

An example file named toricdomainvolumetric. The aim of this example is to demonstrate the volumetric analysis functions on toric domains.

This file is part of the DGtal library.

Definition in file toricdomainvolumetric.cpp.

Function Documentation

â—† main()

int main ( void )

[DTDef]

[DTDef]

[DTPredicate]

[DTPredicate]

[DTComputeToric]

[DTComputeToric]

[DTComputePartialToric]

[DTComputePartialToric]

[DTColormaps]

[DTColormaps]

Definition at line 61 of file toricdomainvolumetric.cpp.

62{
63
65 using namespace std;
66 using namespace DGtal;
67 using namespace Z2i;
68
69 Point a ( 0, 0 );
70 Point b ( 32, 16);
71
72 //Input image with unsigned char values
73 typedef ImageSelector<Domain, unsigned int>::Type Image;
74 Image image ( Domain(a, b ));
75
76 //We fill the image with the 128 value
77 for ( Image::Iterator it = image.begin(), itend = image.end();it != itend; ++it)
78 (*it)=128;
79
80 //We add 3 seeds with 0 values.
81 image.setValue(Point(16,2), 0);
82 image.setValue(Point(2,11), 0);
83 image.setValue(Point(30,15), 0);
85
86 trace.beginBlock ( "Example toricdomainvolumetric" );
87 //Input shape output
88 typedef GrayscaleColorMap<Image::Value> Gray;
89 Board2D board;
90 board.setUnit ( LibBoard::Board::UCentimeter );
91 Display2DFactory::drawImage<Gray>(board, image, (unsigned int)0, (unsigned int)129);
92 board.saveSVG("toric-inputShape.svg");
93
95 //Point Predicate from random seed image
96 typedef functors::SimpleThresholdForegroundPredicate<Image> PointPredicate;
97 PointPredicate predicate(image,0);
99
101 typedef DistanceTransformation<Space, PointPredicate, L2Metric> DTL2;
102 typedef DistanceTransformation<Space, PointPredicate, L2Metric> DTL2Toric;
103
104 //Regular 2D domain
105 DTL2 dtL2(image.domain(), predicate, l2Metric);
106 //Full toric 2D domain
107 DTL2Toric dtL2Toric(image.domain(), predicate, l2Metric, {{true, true}} );
109
111 typedef DistanceTransformation<Space, PointPredicate, L2Metric> DTL2ToricX;
112 typedef DistanceTransformation<Space, PointPredicate, L2Metric> DTL2ToricY;
113
114 // 2D domain that is periodic along the first dimension.
115 DTL2ToricX dtL2ToricX( image.domain(), predicate, l2Metric, {{true, false}} );
116 // 2D domain that is periodic along the second dimension.
117 DTL2ToricY dtL2ToricY( image.domain(), predicate, l2Metric, {{false, true}} );
119
120 //We compute the maximum DT value on the L2 map
121 const DTL2::Value maxv2 = * (boost::first_max_element(dtL2.constRange().begin(), dtL2.constRange().end()));
122 const DTL2Toric::Value maxvtoric = * (boost::first_max_element(dtL2Toric.constRange().begin(), dtL2Toric.constRange().end()));
123 const DTL2ToricX::Value maxvtoricX = * (boost::first_max_element(dtL2ToricX.constRange().begin(), dtL2ToricX.constRange().end()));
124 const DTL2ToricY::Value maxvtoricY = * (boost::first_max_element(dtL2ToricY.constRange().begin(), dtL2ToricY.constRange().end()));
125
126 // Color map based on the maximal value for all maps (in order to compare results with similar colors).
127 const auto maxvall = std::max( { maxv2, maxvtoric, maxvtoricX, maxvtoricY } );
128
130 //Colormap used for the SVG output
131 typedef HueShadeColorMap<DTL2::Value, 1> HueTwice;
133
134 trace.warning() << "DT maxValue= " << maxv2 << endl;
135 board.clear();
136 Display2DFactory::drawImage<HueTwice>(board, dtL2, 0.0, maxvall + 1);
137 board.saveSVG ( "toric-example-DT-L2.svg" );
138
139 trace.warning() << "Full toric maxValue= " << maxvtoric << endl;
140 board.clear();
141 Display2DFactory::drawImage<HueTwice>(board, dtL2Toric, 0.0, maxvall + 1);
142 board.saveSVG ( "toric-example-DT-L2-toric.svg" );
143
144 trace.warning() << "1th dimension periodic maxValue= " << maxvtoricX << endl;
145 board.clear();
146 Display2DFactory::drawImage<HueTwice>(board, dtL2ToricX, 0.0, maxvall + 1);
147 board.saveSVG ( "toric-example-DT-L2-toricX.svg" );
148
149 trace.warning() << "2nd dimension periodic maxValue= " << maxvtoricY << endl;
150 board.clear();
151 Display2DFactory::drawImage<HueTwice>(board, dtL2ToricY, 0.0, maxvall + 1);
152 board.saveSVG ( "toric-example-DT-L2-toricY.svg" );
153
154 //Explicit export with ticked colormap
155 //We compute the maximum DT value on the L2 map
156 TickedColorMap<double, GradientColorMap<double> > ticked(0.0,maxv2, Color::White);
157 ticked.addRegularTicks(3, 0.5);
158 ticked.finalize();
159 ticked.colormap()->addColor( Color::Red );
160 ticked.colormap()->addColor( Color::Black );
161 board.clear();
162 for ( auto it = dtL2.domain().begin(), itend = dtL2.domain().end();it != itend; ++it)
163 {
164 board<< CustomStyle((*it).className(),new CustomColors(ticked(dtL2(*it)),ticked(dtL2(*it))));
165 board << *it;
166 }
167 board.saveSVG("toric-example-DT-L2-ticked.svg");
168
169 board.clear();
170 for ( auto it = dtL2Toric.domain().begin(), itend = dtL2Toric.domain().end();it != itend; ++it)
171 {
172 board<< CustomStyle((*it).className(),new CustomColors(ticked(dtL2Toric(*it)),ticked(dtL2Toric(*it))));
173 board << *it;
174 }
175 board.saveSVG("toric-example-DT-L2-ticked-toric.svg");
176
177 board.clear();
178 for ( auto it = dtL2ToricX.domain().begin(), itend = dtL2ToricX.domain().end();it != itend; ++it)
179 {
180 board<< CustomStyle((*it).className(),new CustomColors(ticked(dtL2ToricX(*it)),ticked(dtL2ToricX(*it))));
181 board << *it;
182 }
183 board.saveSVG("toric-example-DT-L2-ticked-toricX.svg");
184
185 board.clear();
186 for ( auto it = dtL2ToricY.domain().begin(), itend = dtL2ToricY.domain().end();it != itend; ++it)
187 {
188 board<< CustomStyle((*it).className(),new CustomColors(ticked(dtL2ToricY(*it)),ticked(dtL2ToricY(*it))));
189 board << *it;
190 }
191 board.saveSVG("toric-example-DT-L2-ticked-toricY.svg");
192
193 //Voronoi vector output
194 board.clear();
195 board << dtL2.domain();
196 for ( auto it = dtL2.domain().begin(), itend = dtL2.domain().end();it != itend; ++it)
197 if ( dtL2.getVoronoiSite(*it) != *it )
198 Display2DFactory::draw(board,dtL2.getVoronoiSite(*it) - (*it), (*it));
199 board.saveSVG("toric-example-Voro-L2.svg");
200
201 board.clear();
202 board << dtL2Toric.domain();
203 for ( auto it = dtL2Toric.domain().begin(), itend = dtL2Toric.domain().end();it != itend; ++it)
204 if ( dtL2Toric.getVoronoiSite(*it) != *it )
205 Display2DFactory::draw(board, dtL2Toric.getVoronoiSite(*it) - (*it), (*it));
206 board.saveSVG("toric-example-Voro-L2-toric.svg");
207
208 board.clear();
209 board << dtL2Toric.domain();
210 for ( auto it = dtL2Toric.domain().begin(), itend = dtL2Toric.domain().end();it != itend; ++it)
211 if ( dtL2Toric.getVoronoiSite(*it) != *it )
212 Display2DFactory::draw(board, dtL2Toric.projectPoint(dtL2Toric.getVoronoiSite(*it)) - (*it), (*it));
213 board.saveSVG("toric-example-Voro-L2-toric-projected.svg");
214
215 board.clear();
216 board << dtL2ToricX.domain();
217 for ( auto it = dtL2ToricX.domain().begin(), itend = dtL2ToricX.domain().end();it != itend; ++it)
218 if ( dtL2ToricX.getVoronoiSite(*it) != *it )
219 Display2DFactory::draw(board, dtL2ToricX.getVoronoiSite(*it) - (*it), (*it));
220 board.saveSVG("toric-example-Voro-L2-toricX.svg");
221
222 board.clear();
223 board << dtL2ToricY.domain();
224 for ( auto it = dtL2ToricY.domain().begin(), itend = dtL2ToricY.domain().end();it != itend; ++it)
225 if ( dtL2ToricY.getVoronoiSite(*it) != *it )
226 Display2DFactory::draw(board, dtL2ToricY.getVoronoiSite(*it) - (*it), (*it));
227 board.saveSVG("toric-example-Voro-L2-toricY.svg");
228
229 trace.endBlock();
230 return 0;
231}
DGtal::Board2D
Aim: This class specializes a 'Board' class so as to display DGtal objects more naturally (with <<)....
Definition Board2D.h:71
DGtal::Color::Red
static const Color Red
Definition Color.h:425
DGtal::Color::Black
static const Color Black
Definition Color.h:422
DGtal::Color::White
static const Color White
Definition Color.h:424
DGtal::DistanceTransformation
Aim: Implementation of the linear in time distance transformation for separable metrics.
Definition DistanceTransformation.h:99
DGtal::GrayscaleColorMap
Aim: This class template may be used to (linearly) convert scalar values in a given range into gray l...
Definition GrayscaleColorMap.h:94
DGtal::HueShadeColorMap
Aim: This class template may be used to (linearly) convert scalar values in a given range into a colo...
Definition HueShadeColorMap.h:91
DGtal::TickedColorMap
Aim: This class adapts any colormap to add "ticks" in the colormap colors.
Definition TickedColorMap.h:79
DGtal::functors::SimpleThresholdForegroundPredicate
Aim: Define a simple Foreground predicate thresholding image values given a single thresold....
Definition SimpleThresholdForegroundPredicate.h:66
LibBoard::Board::clear
void clear(const DGtal::Color &color=DGtal::Color::None)
Definition Board.cpp:151
LibBoard::Board::saveSVG
void saveSVG(const char *filename, PageSize size=Board::BoundingBox, double margin=10.0) const
Definition Board.cpp:1011
LibBoard::Board::setUnit
void setUnit(Unit unit)
Definition Board.cpp:239
DGtal::Z2i
Z2i this namespace gathers the standard of types for 2D imagery.
DGtal
DGtal is the top-level namespace which contains all DGtal functions and types.
Definition ClosedIntegerHalfPlane.h:49
DGtal::trace
Trace trace
std
STL namespace.
DGtal::CustomColors
Custom style class redefining the pen color and the fill color. You may use Board2D::Color::None for ...
Definition Board2D.h:279
DGtal::Display2DFactory::drawImage
static void drawImage(DGtal::Board2D &board, const Image &i, const typename Image::Value &minV, const typename Image::Value &maxV)
DGtal::Display2DFactory::draw
static void draw(DGtal::Board2D &board, const DGtal::DiscreteExteriorCalculus< dimEmbedded, dimAmbient, TLinearAlgebraBackend, TInteger > &calculus)
DGtal::ImageSelector::Type
ImageContainerBySTLVector< Domain, Value > Type
Definition ImageSelector.h:84
Point
MyPointD Point
Definition testClone2.cpp:381
image
Image image(domain)
Image
ImageContainerBySTLVector< Domain, Value > Image
Definition testSimpleRandomAccessRangeFromPoint.cpp:45
Domain
HyperRectDomain< Space > Domain
Definition testSimpleRandomAccessRangeFromPoint.cpp:44

References DGtal::Display2DFactory::drawImage(), image(), LibBoard::Board::saveSVG(), LibBoard::Board::setUnit(), DGtal::trace, and LibBoard::Board::UCentimeter.