testObjectBorder.cpp File Reference

#include <iostream>
#include <iterator>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/domains/DomainPredicate.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/kernel/sets/DigitalSetSelector.h"
#include "DGtal/kernel/sets/DigitalSetConverter.h"
#include "DGtal/topology/MetricAdjacency.h"
#include "DGtal/topology/DomainMetricAdjacency.h"
#include "DGtal/topology/DomainAdjacency.h"
#include "DGtal/topology/DigitalTopology.h"
#include "DGtal/topology/Object.h"
#include "DGtal/graph/Expander.h"
#include "DGtal/io/boards/Board2D.h"

Go to the source code of this file.

Functions
bool	testObjectBorder ()
bool	testBoard2D ()
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

2010/09/17

Functions for testing class ObjectBorder.

This file is part of the DGtal library.

Definition in file testObjectBorder.cpp.

Function Documentation

main()

int main

(

void

)

Definition at line 309 of file testObjectBorder.cpp.

{
    bool res = testObjectBorder()
               && testBoard2D();
    return res ? 0 : 1;
}

References testBoard2D(), and testObjectBorder().

testBoard2D()

bool testBoard2D

(

)

Simple test of Board2D. Illustrates the border extraction of a simple 2D object considering different topologies.

Definition at line 232 of file testObjectBorder.cpp.

{
    trace.beginBlock ( "Testing Board2D with Object Borders in 2D ..." );
 
    //typedef int Integer;                // choose your digital line here.
    typedef SpaceND<2> Z2;          // Z^2
    typedef Z2::Point Point;
    typedef MetricAdjacency<Z2, 1> Adj4; // 4-adjacency type
    typedef MetricAdjacency<Z2, 2> Adj8; // 8-adjacency type
    typedef DigitalTopology< Adj8, Adj4 > DT8_4; //8,4 topology type
    typedef HyperRectDomain< Z2 > Domain;
    typedef Domain::ConstIterator DomainConstIterator;
    typedef DigitalSetSelector < Domain, BIG_DS + HIGH_BEL_DS >::Type DigitalSet;
    typedef Object<DT8_4, DigitalSet> ObjectType;
 
 
    Point p1 ( -20, -10 );
    Point p2 ( 20, 10 );
    Domain domain ( p1, p2 );
 
    Adj4 adj4;                          // instance of 4-adjacency
    Adj8 adj8;                          // instance of 8-adjacency
    DT8_4 dt8_4 ( adj8, adj4, JORDAN_DT );
 
    Point c ( 0, 0 );
 
    //We construct a simple 3-bubbles set
    DigitalSet bubble_set ( domain );
    for ( DomainConstIterator it = domain.begin(); it != domain.end(); ++it )
    {
        int x = ( *it ) [0];
        int y = ( *it ) [1];
        if ( ( x*x + y*y < 82 ) ||
                ( ( x - 14 ) * ( x - 14 ) + ( y + 1 ) * ( y + 1 ) < 17 ) ||
                ( ( x + 14 ) * ( x + 14 ) + ( y - 1 ) * ( y - 1 ) < 17 ) )
            bubble_set.insertNew ( *it );
    }
 
    ObjectType bubble ( dt8_4, bubble_set );
 
    //Connectedness Check
    if (bubble.computeConnectedness() == CONNECTED)
        trace.info() << "The object is (8,4)connected." << endl;
    else
        trace.info() << "The object is not (8,4)connected." << endl;
 
    //Border Computation
    ObjectType bubbleBorder = bubble.border();
    if (bubbleBorder.computeConnectedness() == CONNECTED)
        trace.info() << "The object (8,4) border is connected." << endl;
    else
        trace.info() << "The object (8,4) border is not connected." << endl;
 
    //Board Export
    Board2D board;
    board.setUnit ( Board::UCentimeter );
 
    board << SetMode( domain.className(), "Grid" )
    << CustomStyle ( domain.className(), new MyDrawStyleCustomGreen )
    << domain
    << CustomStyle ( bubble_set.className(), new MyDrawStyleCustomRed )
    << bubble_set;
    board.saveSVG ( "bubble-set-dgtalboard.svg" );
 
    board << SetMode( bubbleBorder.className(), "DrawAdjacencies" )
    << CustomStyle ( bubbleBorder.className(), new MyDrawStyleCustomBlue )
    << bubbleBorder;
    board.saveSVG ( "bubble-object-border-dgtalboard.svg" );
    board.clear();
 
    trace.endBlock();
    return true;
}

References DGtal::HyperRectDomain< TSpace >::begin(), DGtal::Trace::beginBlock(), DGtal::HyperRectDomain< TSpace >::className(), DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::className(), LibBoard::Board::clear(), DGtal::CONNECTED, domain, DGtal::HyperRectDomain< TSpace >::end(), DGtal::Trace::endBlock(), DGtal::Trace::info(), DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::insertNew(), DGtal::JORDAN_DT, LibBoard::Board::saveSVG(), LibBoard::Board::setUnit(), DGtal::trace, and LibBoard::Board::UCentimeter.

Referenced by main().

testObjectBorder()

bool testObjectBorder

(

)

Simple test to illustrate the border extraction of a simple 2D object considering different topologies.

Definition at line 113 of file testObjectBorder.cpp.

{
    trace.beginBlock ( "Testing Object Borders in 2D ..." );
 
    typedef SpaceND<2> Z2;          // Z^2
    typedef Z2::Point Point;
    typedef MetricAdjacency<Z2, 1> Adj4; // 4-adjacency type
    typedef MetricAdjacency<Z2, 2> Adj8; // 8-adjacency type
    typedef DigitalTopology< Adj8, Adj4 > DT8_4; //8,4 topology type
    typedef HyperRectDomain< Z2 > Domain;
    typedef Domain::ConstIterator DomainConstIterator;
    typedef DigitalSetSelector < Domain, BIG_DS + HIGH_BEL_DS >::Type DigitalSet;
    typedef Object<DT8_4, DigitalSet> ObjectType;
 
 
    Point p1 ( -20, -10 );
    Point p2 ( 20, 10 );
    Domain domain ( p1, p2 );
 
    Adj4 adj4;                          // instance of 4-adjacency
    Adj8 adj8;                          // instance of 8-adjacency
    DT8_4 dt8_4 ( adj8, adj4, JORDAN_DT );
 
    Point c ( 0, 0 );
 
    //We construct a simple 3-bubbles set
    DigitalSet bubble_set ( domain );
    for ( DomainConstIterator it = domain.begin(); it != domain.end(); ++it )
    {
        int x = ( *it ) [0];
        int y = ( *it ) [1];
        if ( ( x*x + y*y < 82 ) ||
                ( ( x - 14 ) * ( x - 14 ) + ( y + 1 ) * ( y + 1 ) < 17 ) ||
                ( ( x + 14 ) * ( x + 14 ) + ( y - 1 ) * ( y - 1 ) < 17 ) )
            bubble_set.insertNew ( *it );
    }
 
    ObjectType bubble ( dt8_4, bubble_set );
 
    //Connectedness Check
    if (bubble.computeConnectedness() == CONNECTED)
        trace.info() << "The object is (8,4)connected." << endl;
    else
        trace.info() << "The object is not (8,4)connected." << endl;
 
    //Border Computation
    ObjectType bubbleBorder = bubble.border();
    if (bubbleBorder.computeConnectedness() == CONNECTED)
        trace.info() << "The object (8,4) border is connected." << endl;
    else
        trace.info() << "The object (8,4) border is not connected." << endl;
 
    //Board Export
    Board2D board;
    board.setUnit ( Board::UCentimeter );
 
    board << SetMode( domain.className(), "Grid" ) <<  domain << bubble_set;
    board.saveSVG ( "bubble-set.svg" );
 
    board << SetMode( bubbleBorder.className(), "DrawAdjacencies" )
          << CustomStyle ( bubbleBorder.className(), new MyObjectStyleCustom )
          << bubbleBorder;
    board.saveSVG ( "bubble-object-border.svg" );
 
    board.clear();
 
    //the same with the reverse topology
    typedef Object<DT8_4::ReverseTopology, DigitalSet> ObjectType48;
    DT8_4::ReverseTopology dt4_8 = dt8_4.reverseTopology();
 
    ObjectType48 bubble2 ( dt4_8, bubble_set );
 
    //Border Computation
    ObjectType48 bubbleBorder2 = bubble2.border();
    if (bubbleBorder2.computeConnectedness() == CONNECTED)
        trace.info() << "The object (4,8) border is connected." << endl;
    else
        trace.info() << "The object (4,8) border is not connected." << endl;
 
    board << SetMode( domain.className(), "Grid" ) << domain;
    board << bubble_set
          << SetMode( bubbleBorder2.className(), "DrawAdjacencies" )
          << CustomStyle ( bubbleBorder2.className(), new MyObjectStyleCustom )
          << bubbleBorder2;
 
    board.saveSVG ( "bubble-object-border-48.svg" );
 
    //We split the border according to its components
    vector<ObjectType48> borders ( 30 );
    unsigned int nbComponents;
 
    vector<ObjectType48>::iterator it = borders.begin();
    nbComponents = bubbleBorder2.writeComponents ( it );
 
    trace.info() << "The Bubble object has " << nbComponents << " (4,8)-connected components" << endl;
 
    bool flag = true;
    for ( unsigned int k = 0;k < nbComponents ; k++ )
    {
        if ( flag )
            board << SetMode( borders[k].className(), "DrawAdjacencies" ) << CustomStyle ( borders[k].className(), new MyObjectStyleCustom ) << borders[k];
        else
            board << SetMode( borders[k].className(), "DrawAdjacencies" ) << CustomStyle ( borders[k].className(), new MyObjectStyleCustom ) << borders[k];
        flag = !flag;
    }
 
    board.saveSVG ( "bubble-object-color-borders-48.svg" );
    trace.endBlock();
 
    return true;
}

References DGtal::HyperRectDomain< TSpace >::begin(), DGtal::Trace::beginBlock(), DGtal::HyperRectDomain< TSpace >::className(), LibBoard::Board::clear(), DGtal::CONNECTED, domain, DGtal::HyperRectDomain< TSpace >::end(), DGtal::Trace::endBlock(), DGtal::Trace::info(), DGtal::DigitalSetByAssociativeContainer< TDomain, TContainer >::insertNew(), DGtal::JORDAN_DT, LibBoard::Board::saveSVG(), LibBoard::Board::setUnit(), DGtal::trace, and LibBoard::Board::UCentimeter.

Referenced by main().