DGtal  0.9.3
lower-integer-convex-hull.cpp
Go to the documentation of this file.
1 
42 #include "DGtal/base/Common.h"
45 #include "DGtal/arithmetic/LatticePolytope2D.h"
46 #include "DGtal/io/boards/Board2D.h"
47 #include "DGtal/shapes/Shapes.h"
48 #include "DGtal/helpers/StdDefs.h"
50 
52 
53 using namespace DGtal;
54 
56 
57 void usage( int, char** argv )
58 {
59  std::cerr << "Usage: " << argv[ 0 ] << " <a> <b> <c>" << std::endl;
60  std::cerr << "\t - Cuts a square by ax+by <= c. Displays everything in files "
61  << "lower-integer-convex-hull*.eps" << std::endl;
62 }
63 
67 int main( int argc, char** argv )
68 {
69  if ( argc < 4 )
70  {
71  usage( argc, argv );
72  return 0;
73  }
74 
76  using namespace Z2i;
77  typedef LatticePolytope2D<Space> CIP;
79 
81  CIP cip;
82  cip.push_front( Point( -10, -10 ) );
83  cip.push_front( Point( -10, 10 ) );
84  cip.push_front( Point( 10, 10 ) );
85  cip.push_front( Point( 10, -10 ) );
86  Domain domain = cip.boundingBoxDomain();
87  Board2D board;
88  board << domain
89  << CustomStyle( cip.className(),
91  << cip;
92  board.saveEPS( "lower-integer-convex-hull.eps" );
93  board.clear();
95 
96  int a = atoi( argv[ 1 ] );
97  int b = atoi( argv[ 2 ] );
98  int c = atoi( argv[ 3 ] );
99 
101  typedef LatticePolytope2D<Z2>::HalfSpace HalfSpace;
102  HalfSpace hs( Vector( a, b ), c );
103  cip.cut( hs );
104  DigitalSet aSet( domain );
106  board << domain
107  << CustomStyle( aSet.className(),
109  << SetMode( Point().className(), "Grid" )
110  << aSet
111  << CustomStyle( cip.className(),
113  << cip;
114  board.saveEPS( "lower-integer-convex-hull-cut.eps" );
116 
118  std::cout << "Number of vertices = " << cip.size() << std::endl;
119  std::cout << "Area = " << (((double)cip.twiceArea())/2.0) << std::endl;
120  std::cout << "Number of interior points = " << cip.numberInteriorPoints() << std::endl;
121  std::cout << "Number of boundary points = " << cip.numberBoundaryPoints() << std::endl;
123  return 0;
124 }
125 
DigitalSetSelector< Domain, BIG_DS+HIGH_BEL_DS >::Type DigitalSet
Definition: StdDefs.h:100
void push_front(const Point &K)
const Domain domain(Point(1, 2), Point(6, 5))
Custom style class redefining the pen color and the fill color. You may use Board2D::Color::None for ...
Definition: Board2D.h:278
static const Color None
Definition: Color.h:387
void usage(int, char **argv)
void clear(const DGtal::Color &color=DGtal::Color::None)
Definition: Board.cpp:152
void saveEPS(const char *filename, PageSize size=Board::BoundingBox, double margin=10.0) const
Definition: Board.cpp:805
Space::Point Point
Definition: StdDefs.h:95
DGtal is the top-level namespace which contains all DGtal functions and types.
Aim: A half-space specified by a vector N and a constant c. The half-space is the set ...
Modifier class in a Board2D stream. Useful to choose your own mode for a given class. Realizes the concept CDrawableWithBoard2D.
Definition: Board2D.h:247
static const Color Red
Definition: Color.h:391
std::string className() const
static const Color Green
Definition: Color.h:392
Space::Vector Vector
Definition: StdDefs.h:96
Aim: Represents a 2D polytope, i.e. a convex polygon, in the two-dimensional digital plane...
int main(int argc, char **argv)
static void makeSetFromPointPredicate(DigitalSet &aSet, const PointPredicate &aPP)
Aim: This class specializes a &#39;Board&#39; class so as to display DGtal objects more naturally (with <<)...
Definition: Board2D.h:70