exampleMaximalSegmentSliceEstimation.cpp File Reference

#include <iostream>
#include "ConfigExamples.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/base/Common.h"
#include "DGtal/geometry/surfaces/estimation/MaximalSegmentSliceEstimation.h"
#include "DGtal/helpers/Shortcuts.h"
#include "DGtal/io/viewers/Viewer3D.h"

Include dependency graph for exampleMaximalSegmentSliceEstimation.cpp:

Go to the source code of this file.

Typedefs
using	KSpace = Z3i::KSpace
using	SH3 = Shortcuts< KSpace >
using	Surface = SH3::DigitalSurface
using	Surfel = Surface::Surfel
using	RealPoint = SH3::RealPoint
using	Point = SH3::Point

Functions
RealPoint	centerSurfel (KSpace const &K, SH3::SCell const &s)
int	main (int argc, char **argv)

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

Jocelyn Meyron (jocel.nosp@m.yn.m.nosp@m.eyron.nosp@m.@lir.nosp@m.is.cn.nosp@m.rs.f.nosp@m.r ) Laboratoire d'InfoRmatique en Image et Systemes d'information - LIRIS (CNRS, UMR 5205), CNRS, France

Date

2020/12/15

An example file that shows how to construct and use a DGtal::MaximalSegmentSliceEstimation on a digital surface.

This file is part of the DGtal library.

Definition in file exampleMaximalSegmentSliceEstimation.cpp.

Typedef Documentation

KSpace

using KSpace = Z3i::KSpace

Definition at line 43 of file exampleMaximalSegmentSliceEstimation.cpp.

Point

using Point = SH3::Point

Definition at line 48 of file exampleMaximalSegmentSliceEstimation.cpp.

RealPoint

using RealPoint = SH3::RealPoint

Definition at line 47 of file exampleMaximalSegmentSliceEstimation.cpp.

SH3

using SH3 = Shortcuts<KSpace>

Definition at line 44 of file exampleMaximalSegmentSliceEstimation.cpp.

Surface

using Surface = SH3::DigitalSurface

Examples

geometry/surfaces/dvcm-2d-curvature.cpp.

Definition at line 45 of file exampleMaximalSegmentSliceEstimation.cpp.

Surfel

using Surfel = Surface::Surfel

Definition at line 46 of file exampleMaximalSegmentSliceEstimation.cpp.

Function Documentation

centerSurfel()

RealPoint centerSurfel	(	KSpace const &	K,
		SH3::SCell const &	s
	)

Definition at line 50 of file exampleMaximalSegmentSliceEstimation.cpp.

{
    auto pointels = SH3::getPrimalVertices(K, s, true);
 
    Point p = K.uCoords(pointels[0]),
          u = K.uCoords(pointels[1]) - p,
          v = K.uCoords(pointels[3]) - p;
 
    static const RealPoint shift(-0.5, -0.5, -0.5);
 
    return p + 0.5 * u + 0.5 * v + shift;
}

References K, and DGtal::KhalimskySpaceND< dim, TInteger >::uCoords().

Referenced by main().

main()

int main	(	int	argc,
		char **	argv
	)

[MaximalSegmentSliceEstimationConstructionUsage]

[MaximalSegmentSliceEstimationConstructionUsage]

Definition at line 65 of file exampleMaximalSegmentSliceEstimation.cpp.

{
    QApplication application(argc, argv);
 
    auto params = SH3::defaultParameters();
    params("polynomial", "ellipsoid")("gridstep", "1.0");
    auto implicit_shape  = SH3::makeImplicitShape3D  ( params );
    auto digitized_shape = SH3::makeDigitizedImplicitShape3D( implicit_shape, params );
    auto K               = SH3::getKSpace( params );
    auto binary_image    = SH3::makeBinaryImage( digitized_shape, params );
    auto surface         = SH3::makeDigitalSurface( binary_image, K, params );
    auto surfels         = SH3::getSurfelRange(surface);
 
    double gridstep = params["gridstep"].as<double>();
 
    Viewer3D<> viewer(K);
    viewer << SetMode3D(Surfel().className(), "Basic");
    viewer.show();
 
    // Instantiation
    using Estimator = MaximalSegmentSliceEstimation<Surface>;
    Estimator estimator;
    estimator.init(gridstep, surfels.begin(), surfels.end());
    estimator.attach(surface);
 
    // Usage
    std::vector<Estimator::Quantity> quantities;
    estimator.eval(surfels.begin(), surfels.end(), std::back_inserter(quantities));
 
    Color fillColor = viewer.getFillColor();
 
    for (std::size_t i = 0; i < surfels.size(); ++i)
    {
        const Surfel& s = surfels[i];
        const Estimator::Quantity& normal = quantities[i];
 
        const RealPoint& n = normal.getNormalized();
        RealPoint origin = centerSurfel(K, s);
 
        viewer.setFillColor(fillColor);
        viewer << s;
 
        viewer.setLineColor(Color::Green);
        viewer.addLine(origin, origin + 1.5 * n, 0.3);
    }
 
    viewer << Viewer3D<>::updateDisplay;
    application.exec();
 
    return 0;
}

References centerSurfel(), DGtal::PointVector< dim, TEuclideanRing, TContainer >::getNormalized(), and K.