DGtal  1.5.beta
Board3DTo2D.ih
1 /**
2  * This program is free software: you can redistribute it and/or modify
3  * it under the terms of the GNU Lesser General Public License as
4  * published by the Free Software Foundation, either version 3 of the
5  * License, or (at your option) any later version.
6  *
7  * This program is distributed in the hope that it will be useful,
8  * but WITHOUT ANY WARRANTY; without even the implied warranty of
9  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
10  * GNU General Public License for more details.
11  *
12  * You should have received a copy of the GNU General Public License
13  * along with this program. If not, see <http://www.gnu.org/licenses/>.
14  *
15  **/
16 
17 /**
18  * @file Board3DTo2D.ih
19  * @author Martial Tola <http://liris.cnrs.fr/martial.tola/>
20  * @date mercredi 22 juin 2011
21  *
22  * @brief
23  *
24  * Implementation of inline methods defined in Board3DTo2D.h
25  *
26  * This file is part of the DGtal library.
27  */
28 
29 ///////////////////////////////////////////////////////////////////////////////
30 // IMPLEMENTATION of inline methods.
31 ///////////////////////////////////////////////////////////////////////////////
32 
33 #include <limits>
34 
35 // Cairo includes
36 #if defined(__clang__)
37 #pragma clang diagnostic push
38 #pragma clang diagnostic ignored "-Wdocumentation"
39 #endif
40 #include <cairo.h>
41 #include <cairo-pdf.h>
42 #include <cairo-ps.h>
43 #include <cairo-svg.h>
44 #if defined(__clang__)
45 #pragma clang diagnostic pop
46 #endif
47 // Cairo includes
48 //////////////////////////////////////////////////////////////////////////////
49 #include <cstdlib>
50 #include "DGtal/io/CDrawableWithDisplay3D.h"
51 #include "DGtal/io/boards/CDrawableWithBoard3DTo2D.h"
52 #include "DGtal/io/Color.h"
53 #include "DGtal/io/boards/Board3DTo2DFactory.h"
54 
55 //////////////////////////////////////////////////////////////////////////////
56 
57 
58 ///////////////////////////////////////////////////////////////////////////////
59 // //
60 // //
61 ///////////////////////////////////////////////////////////////////////////////
62 
63 
64 ///////////////////////////////////////////////////////////////////////////////
65 // Implementation of inline methods //
66 
67 /**
68  * Set the default color for future drawing.
69  * @param aColor: a DGtal::Color (allow to set a trasnparency value).
70  *
71  **/
72 template < typename Space, typename KSpace>
73 inline
74 DGtal::Board3DTo2D<Space, KSpace> &
75 DGtal::Board3DTo2D<Space, KSpace>::operator<<(const DGtal::Color & aColor)
76 {
77  myDefaultColor=aColor;
78  return *this;
79 }
80 
81 
82 /**
83  * Draws the drawable [object] in this board. It should satisfy
84  * the concept CDrawableWithBoard3DTo2D, which requires for instance a
85  * method setStyle( Board3DTo2D & ).
86  *
87  * @param object any drawable object.
88  * @return a reference on 'this'.
89  */
90 template < typename Space, typename KSpace>
91 template <typename TDrawableWithBoard3DTo2D>
92 inline
93 DGtal::Board3DTo2D<Space, KSpace> &
94 DGtal::Board3DTo2D<Space, KSpace>::operator<<( const TDrawableWithBoard3DTo2D & object )
95 {
96  BOOST_CONCEPT_ASSERT((concepts::CDrawableWithBoard3DTo2D< TDrawableWithBoard3DTo2D, Space, KSpace>));
97 
98  DGtal::Board3DTo2DFactory<Space,KSpace>::draw(*this, object);
99  return *this;
100 }
101 
102 // //
103 ///////////////////////////////////////////////////////////////////////////////
104 
105 
106 ///////////////////////////////////////////////////////////////////////////////
107 // Standard services - public :
108 
109 /*!
110  * \brief Constructor.
111  */
112 template < typename Space, typename KSpace>
113 DGtal::Board3DTo2D<Space, KSpace>::Board3DTo2D()
114 {
115  init();
116 }
117 // //
118 ///////////////////////////////////////////////////////////////////////////////
119 
120 ///////////////////////////////////////////////////////////////////////////////
121 // Interface - public :
122 
123 /**
124  * Writes/Displays the object on an output stream.
125  * @param out the output stream where the object is written.
126  */
127 template < typename Space, typename KSpace>
128 void
129 DGtal::Board3DTo2D<Space, KSpace>::selfDisplay ( std::ostream & out ) const
130 {
131  out << "[Board3DTo2D]";
132 }
133 
134 /**
135  * Checks the validity/consistency of the object.
136  * @return 'true' if the object is valid, 'false' otherwise.
137  */
138 template < typename Space, typename KSpace>
139 bool
140 DGtal::Board3DTo2D<Space, KSpace>::isValid() const
141 {
142  return true;
143 }
144 
145 
146 /**
147  * Transpose a 4x4 matrix.
148  * @param tmat destination matrix.
149  * @param mat source matrix.
150  */
151 template < typename Space, typename KSpace>
152 void
153 DGtal::Board3DTo2D<Space, KSpace>::TransposeMt(double tmat[16], double mat[16])
154 {
155  tmat[0] = mat[0]; tmat[1] = mat[4]; tmat[2] = mat[8]; tmat[3] = mat[12];
156  tmat[4] = mat[1]; tmat[5] = mat[5]; tmat[6] = mat[9]; tmat[7] = mat[13];
157  tmat[8] = mat[2]; tmat[9] = mat[6]; tmat[10] = mat[10]; tmat[11] = mat[14];
158  tmat[12] = mat[3]; tmat[13] = mat[7]; tmat[14] = mat[11]; tmat[15] = mat[15];
159 }
160 
161 /**
162  * Multiply a 3d vector by a 4x4 matrix.
163  * @param v destination vector.
164  * @param mat source matrix.
165  * @param b source vector.
166  */
167 template < typename Space, typename KSpace>
168 void
169 DGtal::Board3DTo2D<Space, KSpace>::MulMt(double v[4], double mat[16], double b[4])
170 {
171  v[0] = mat[0] * b[0] + mat[1] * b[1] + mat[2] * b[2] + mat[3] * b[3];
172  v[1] = mat[4] * b[0] + mat[5] * b[1] + mat[6] * b[2] + mat[7] * b[3];
173  v[2] = mat[8] * b[0] + mat[9] * b[1] + mat[10] * b[2] + mat[11] * b[3];
174  v[3] = mat[12] * b[0] + mat[13] * b[1] + mat[14] * b[2] + mat[15] * b[3];
175 }
176 
177 /**
178  * Compute 4x4 LookAt matrix.
179  * @param mat destination matrix.
180  * @param eyex x position of eye.
181  * @param eyey y position of eye.
182  * @param eyez z position of eye.
183  * @param dirx x direction of eye.
184  * @param diry y direction of eye.
185  * @param dirz z director of eye.
186  * @param upx x coordinate of up-vector.
187  * @param upy y coordinate of up-vector.
188  * @param upz z coordinate of up-vector.
189  */
190 template < typename Space, typename KSpace>
191 void
192 DGtal::Board3DTo2D<Space, KSpace>::LookAtMt(double mat[16],
193  double eyex, double eyey, double eyez,
194  double dirx, double diry, double dirz,
195  double upx, double upy, double upz)
196 {
197  double up[3]; up[0]= upx; up[1]= upy; up[2]= upz;
198 
199  double z[3]; z[0]= -dirx; z[1]= -diry; z[2]= -dirz; Display3D<Space,KSpace>::normalize(z);
200  double x[3]; Display3D<Space,KSpace>::cross (x, up, z); Display3D<Space,KSpace>::normalize(x);
201  double y[3]; Display3D<Space,KSpace>::cross (y, z, x); Display3D<Space,KSpace>::normalize(y);
202 
203  double m[16];
204  m[0] = x[0]; m[1] = x[1]; m[2] = x[2]; m[3] = 0;
205  m[4] = y[0]; m[5] = y[1]; m[6] = y[2]; m[7] = 0;
206  m[8] = z[0]; m[9] = z[1]; m[10] = z[2]; m[11] = 0;
207  m[12] = 0; m[13] = 0; m[14] = 0; m[15] = 1;
208 
209  double e[4]; e[0]= -eyex; e[1]= -eyey; e[2]= -eyez; e[3]= 1;
210  double eyePrime[4]; MulMt(eyePrime, m, e);
211 
212  TransposeMt(mat, m);
213  mat[12] = eyePrime[0]; mat[13] = eyePrime[1]; mat[14] = eyePrime[2];
214 }
215 
216 /**
217  * Precompute 4x4 projection matrix for 3D->2D projection.
218  */
219 template < typename Space, typename KSpace>
220 void DGtal::Board3DTo2D<Space, KSpace>::precompute_projection_matrix()
221 {
222  // Projection: from qglviewer
223  /*const double f = 1.0/tan(fieldOfView()/2.0);
224  projectionMatrix_[0] = f/aspectRatio();
225  projectionMatrix_[5] = f;
226  projectionMatrix_[10] = (ZNear + ZFar) / (ZNear - ZFar);
227  projectionMatrix_[11] = -1.0;
228  projectionMatrix_[14] = 2.0 * ZNear * ZFar / (ZNear - ZFar);
229  projectionMatrix_[15] = 0.0;
230  // same as gluPerspective( 180.0*fieldOfView()/M_PI, aspectRatio(), zNear(), zFar() );*/
231 
232  double fieldOfView = M_PI/4.;
233  double f = 1.0/tan(fieldOfView/2.0);
234  double aspectRatio = (double)Viewport[2]/Viewport[3];
235 
236  double Projection[16] = { f/aspectRatio, 0.00, 0.00, 0.00,
237  0.00, f, 0.00, 0.00,
238  0.00, 0.00, (ZNear + ZFar) / (ZNear - ZFar), -1.00,
239  0.00, 0.00, 2.0 * ZNear * ZFar / (ZNear - ZFar), 0.00 };
240 
241  double Modelview[16];
242  LookAtMt(Modelview,
243  camera_position[0], camera_position[1], camera_position[2],
244  camera_direction[0], camera_direction[1], camera_direction[2],
245  camera_upVector[0], camera_upVector[1], camera_upVector[2]);
246 
247  for (unsigned short m=0; m<4; ++m)
248  {
249  for (unsigned short l=0; l<4; ++l)
250  {
251  double sum = 0.0;
252  for (unsigned short k=0; k<4; ++k)
253  sum += Projection[l+4*k]*Modelview[k+4*m];
254  matrix[l+4*m] = sum;
255  }
256  }
257 }
258 
259 /**
260  * Project a 3d point (3D->2D).
261  * @param x3d x position of the 3d point.
262  * @param y3d y position of the 3d point.
263  * @param z3d z position of the 3d point.
264  * @param x2d x destination projection position of the 2d point.
265  * @param y2d y destination projection position of the 2d point.
266  */
267 template < typename Space, typename KSpace>
268 void DGtal::Board3DTo2D<Space, KSpace>::project(double x3d, double y3d, double z3d, double &x2d, double &y2d)
269 {
270  double v[4], vs[4];
271  v[0]=x3d; v[1]=y3d; v[2]=z3d; v[3]=1.0;
272 
273  vs[0]=matrix[0 ]*v[0] + matrix[4 ]*v[1] + matrix[8 ]*v[2] + matrix[12 ]*v[3];
274  vs[1]=matrix[1 ]*v[0] + matrix[5 ]*v[1] + matrix[9 ]*v[2] + matrix[13 ]*v[3];
275  vs[2]=matrix[2 ]*v[0] + matrix[6 ]*v[1] + matrix[10]*v[2] + matrix[14 ]*v[3];
276  vs[3]=matrix[3 ]*v[0] + matrix[7 ]*v[1] + matrix[11]*v[2] + matrix[15 ]*v[3];
277 
278  vs[0] /= vs[3];
279  vs[1] /= vs[3];
280  vs[2] /= vs[3];
281 
282  vs[0] = vs[0] * 0.5 + 0.5;
283  vs[1] = vs[1] * 0.5 + 0.5;
284  vs[2] = vs[2] * 0.5 + 0.5;
285 
286  vs[0] = vs[0] * Viewport[2] + Viewport[0];
287  vs[1] = vs[1] * Viewport[3] + Viewport[1];
288 
289  x2d = vs[0];
290  y2d = Viewport[3]-vs[1];
291 }
292 
293 /**
294  * Save a Cairo image.
295  * @param filename filename of the image to save.
296  * @param type type of the image to save (CairoPDF, CairoPNG, CairoPS, CairoEPS, CairoSVG).
297  * @param bWidth width of the image to save.
298  * @param bHeight height of the image to save.
299  */
300 template < typename Space, typename KSpace>
301 void
302 DGtal::Board3DTo2D<Space, KSpace>::saveCairo(const char *filename, CairoType type, int bWidth, int bHeight)
303 {
304  for(unsigned int i =0; i< Board3DTo2D<Space, KSpace>::myClippingPlaneList.size(); i++)
305  trace.info() << "-> ClippingPlane not implemented in Board3DTo2D" << std::endl;
306 
307  Viewport[0] = 0; Viewport[1] = 0; Viewport[2] = bWidth; Viewport[3] = bHeight;
308  precompute_projection_matrix();
309 
310  cairo_surface_t *surface;
311  cairo_t *cr;
312 
313  switch (type)
314  {
315  case CairoPDF:
316  surface = cairo_pdf_surface_create (filename, Viewport[2], Viewport[3]); break;
317  case CairoPS:
318  surface = cairo_ps_surface_create (filename, Viewport[2], Viewport[3]); break;
319  case CairoEPS:
320  surface = cairo_ps_surface_create (filename, Viewport[2], Viewport[3]);
321  cairo_ps_surface_set_eps(surface, true); break;
322  case CairoSVG:
323  surface = cairo_svg_surface_create (filename, Viewport[2], Viewport[3]); break;
324  case CairoPNG:
325  default:
326  surface = cairo_image_surface_create (CAIRO_FORMAT_ARGB32, Viewport[2], Viewport[3]);
327  }
328 
329  cr = cairo_create (surface);
330 
331  // Draw the shapes
332  // myBallSetList
333  for(unsigned int i=0; i<Board3DTo2D<Space, KSpace>::myBallSetList.size(); i++)
334  {
335  for (typename std::vector<typename Board3DTo2D<Space, KSpace>::BallD3D>::iterator s_it = Board3DTo2D<Space, KSpace>::myBallSetList.at(i).begin();
336  s_it != Board3DTo2D<Space, KSpace>::myBallSetList.at(i).end();
337  ++s_it)
338  {
339  {
340  cairo_save (cr);
341 
342  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
343  s_it->color.blue()/255.0, s_it->color.alpha()/255.0);
344  cairo_set_line_width (cr, 1.); // arbitraire car non set
345 
346  double x1, y1, x2, y2, x3, y3, x4, y4;
347  //double width=s_it->size/120.; // arbitraire
348  double width=0.05; // arbitraire
349  // TODO:
350  /*double distCam =sqrt((camera_position[0]-centerS.x)*(camera_position[0]-centerS.x)+
351  (camera_position[1]-centerS.y)*(camera_position[1]-centerS.y)+
352  (camera_position[2]-centerS.z)*(camera_position[2]-centerS.z));*/
353 
354  //z+
355  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]+width, x1, y1);
356  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
357  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]+width, x3, y3);
358  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]+width, x4, y4);
359  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
360  //z-
361  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]-width, x1, y1);
362  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]-width, x2, y2);
363  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]-width, x3, y3);
364  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
365  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
366  //x+
367  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]+width, x1, y1);
368  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
369  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]-width, x3, y3);
370  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
371  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
372  //x-
373  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]+width, x1, y1);
374  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
375  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]-width, x3, y3);
376  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
377  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
378  //y+
379  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]+width, x1, y1);
380  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
381  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]-width, x3, y3);
382  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]-width, x4, y4);
383  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
384  //y-
385  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]+width, x1, y1);
386  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]+width, x2, y2);
387  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]-width, x3, y3);
388  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
389  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
390 
391  cairo_restore (cr);
392  }
393  }
394  }
395 
396  // myLineSetList
397  for(unsigned int i=0; i<Board3DTo2D<Space, KSpace>::myLineSetList.size(); i++)
398  {
399  for (typename std::vector<typename Board3DTo2D<Space, KSpace>::LineD3D>::iterator s_it = Board3DTo2D<Space, KSpace>::myLineSetList.at(i).begin();
400  s_it != Board3DTo2D<Space, KSpace>::myLineSetList.at(i).end();
401  ++s_it)
402  {
403  {
404  cairo_save (cr);
405 
406  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
407  s_it->color.blue()/255.0, s_it->color.alpha()/255.0);
408 
409  double x1, y1;
410  double x2, y2;
411  project(s_it->point1[0], s_it->point1[1], s_it->point1[2], x1, y1);
412  project(s_it->point2[0], s_it->point2[1], s_it->point2[2], x2, y2);
413  cairo_move_to (cr, x1, y1);
414  cairo_line_to (cr, x2, y2);
415 
416  //cairo_set_line_width (cr, s_it->width);
417  cairo_set_line_width (cr, 1.); // arbitraire car non set
418 
419  cairo_stroke (cr);
420 
421  cairo_restore (cr);
422  }
423  }
424  }
425 
426  // myCubeSetList
427  for(typename Board3DTo2D<Space, KSpace>::CubesMap::const_iterator it = Board3DTo2D<Space, KSpace>::myCubesMap.begin();
428  it != Board3DTo2D<Space, KSpace>::myCubesMap.end(); it++)
429  {
430  for (typename std::vector< typename Board3DTo2D<Space, KSpace>::CubeD3D>::const_iterator s_it = it->second.begin();
431  s_it != it->second.end(); ++s_it)
432  {
433  {
434  cairo_save (cr);
435 
436  if (Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode")
437  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
438  s_it->color.blue()/255.0, s_it->color.alpha()/(255.0*1.75)); // *1.75 arbitraire
439  else
440  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
441  s_it->color.blue()/255.0, s_it->color.alpha()/(255.0*0.75)); // *0.75 arbitraire
442 
443  cairo_set_line_width (cr, 1.); // arbitraire car non set
444 
445  double x1, y1, x2, y2, x3, y3, x4, y4;
446  double width=s_it->width;
447 
448  //z+
449  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]+width, x1, y1);
450  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
451  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]+width, x3, y3);
452  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]+width, x4, y4);
453  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
454  //z-
455  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]-width, x1, y1);
456  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]-width, x2, y2);
457  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]-width, x3, y3);
458  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
459  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
460  //x+
461  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]+width, x1, y1);
462  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
463  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]-width, x3, y3);
464  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
465  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
466  //x-
467  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]+width, x1, y1);
468  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
469  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]-width, x3, y3);
470  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
471  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
472  //y+
473  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]+width, x1, y1);
474  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]+width, x2, y2);
475  project(s_it->center[0]+width, s_it->center[1]+width, s_it->center[2]-width, x3, y3);
476  project(s_it->center[0]-width, s_it->center[1]+width, s_it->center[2]-width, x4, y4);
477  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
478  //y-
479  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]+width, x1, y1);
480  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]+width, x2, y2);
481  project(s_it->center[0]+width, s_it->center[1]-width, s_it->center[2]-width, x3, y3);
482  project(s_it->center[0]-width, s_it->center[1]-width, s_it->center[2]-width, x4, y4);
483  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
484 
485  cairo_restore (cr);
486  }
487  }
488  }
489 
490  for(typename Display3D<Space, KSpace>::QuadsMap::iterator it = myQuadsMap.begin(); it != myQuadsMap.end(); it++)
491  {
492  for (typename std::vector< typename Board3DTo2D<Space, KSpace>::QuadD3D>::const_iterator s_it = it->second.begin();
493  s_it != it->second.end(); ++s_it)
494  {
495  {
496  cairo_save (cr);
497 
498  if (Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode")
499  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
500  s_it->color.blue()/255.0, s_it->color.alpha()/(255.0*1.75)); // *1.75 arbitraire
501  else
502  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
503  s_it->color.blue()/255.0, s_it->color.alpha()/(255.0*0.75)); // *0.75 arbitraire
504 
505  cairo_set_line_width (cr, 1.); // arbitraire car non set
506 
507  double x1, y1, x2, y2, x3, y3, x4, y4;
508 
509  project(s_it->point1[0], s_it->point1[1], s_it->point1[2], x1, y1);
510  project(s_it->point2[0], s_it->point2[1], s_it->point2[2], x2, y2);
511  project(s_it->point3[0], s_it->point3[1], s_it->point3[2], x3, y3);
512  project(s_it->point4[0], s_it->point4[1], s_it->point4[2], x4, y4);
513  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
514  cairo_restore (cr);
515  }
516  }
517  }
518 
519  // Drawing all Khalimsky Space Cells
520 
521  // Prism (from updateList)
522  for (typename std::vector<typename Board3DTo2D<Space, KSpace>::QuadD3D>::iterator s_it =Board3DTo2D<Space, KSpace>:: myPrismList.begin();
523  s_it != Board3DTo2D<Space, KSpace>::myPrismList.end();
524  ++s_it)
525  {
526  {
527  cairo_save (cr);
528 
529  if (Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode")
530  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
531  s_it->color.blue()/255.0, s_it->color.alpha()/(255.0*3.75)); // *3.75 arbitraire
532  else
533  cairo_set_source_rgba (cr, s_it->color.red()/255.0, s_it->color.green()/255.0,
534  s_it->color.blue()/255.0, s_it->color.alpha()/(255.0*0.75)); // *0.75 arbitraire
535 
536  cairo_set_line_width (cr, 1.); // arbitraire car non set
537 
538  double x1, y1, x2, y2, x3, y3, x4, y4;
539 
540  project(s_it->point1[0], s_it->point1[1], s_it->point1[2], x1, y1);
541  project(s_it->point2[0], s_it->point2[1], s_it->point2[2], x2, y2);
542  project(s_it->point3[0], s_it->point3[1], s_it->point3[2], x3, y3);
543  project(s_it->point4[0], s_it->point4[1], s_it->point4[2], x4, y4);
544  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr); Board3DTo2D<Space, KSpace>::myModes["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
545 
546  cairo_restore (cr);
547  }
548  }
549 
550  /*
551  // KSLinel
552  for(unsigned int i=0; i< myKSLinelList.size();i++)
553  {
554  //if (!myKSLinelList.at(i).isSigned) // for the moment, same for Signed or NonSigned
555  {
556  {
557  cairo_save (cr);
558 
559  cairo_set_source_rgba (cr, myKSLinelList.at(i).R/255.0, myKSLinelList.at(i).G/255.0, myKSLinelList.at(i).B/255.0, myKSLinelList.at(i).T/255.0);
560  cairo_set_line_width (cr, 4.); // arbitraire car non set
561 
562  double x1, y1, x2, y2;
563 
564  project(myKSLinelList.at(i).x1, myKSLinelList.at(i).y1, myKSLinelList.at(i).z1, x1, y1);
565  project(myKSLinelList.at(i).x2, myKSLinelList.at(i).y2, myKSLinelList.at(i).z2, x2, y2);
566  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_close_path (cr); cairo_stroke (cr);
567 
568  cairo_restore (cr);
569  }
570  }
571  }
572 
573 
574  // KSBallel
575  for(std::vector<ballD3D>::iterator it=myKSBallelList.begin(); it != myKSBallelList.end(); it++)
576  {
577  //if (!it->isSigned) // for the moment, same for Signed or NonSigned
578  {
579  {
580  cairo_save (cr);
581 
582  cairo_set_source_rgba (cr, it->R/255.0, it->G/255.0, it->B/255.0, it->T/255.0);
583  cairo_set_line_width (cr, 1.); // arbitraire car non set
584 
585  double x1, y1, x2, y2, x3, y3, x4, y4;
586  //double width=it->size/120.; // arbitraire
587  double width=0.04; // arbitraire
588 
589  //z+
590  project(it->x-width, it->y+width, it->z+width, x1, y1);
591  project(it->x+width, it->y+width, it->z+width, x2, y2);
592  project(it->x+width, it->y-width, it->z+width, x3, y3);
593  project(it->x-width, it->y-width, it->z+width, x4, y4);
594  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
595  //z-
596  project(it->x-width, it->y+width, it->z-width, x1, y1);
597  project(it->x+width, it->y+width, it->z-width, x2, y2);
598  project(it->x+width, it->y-width, it->z-width, x3, y3);
599  project(it->x-width, it->y-width, it->z-width, x4, y4);
600  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
601  //x+
602  project(it->x+width, it->y-width, it->z+width, x1, y1);
603  project(it->x+width, it->y+width, it->z+width, x2, y2);
604  project(it->x+width, it->y+width, it->z-width, x3, y3);
605  project(it->x+width, it->y-width, it->z-width, x4, y4);
606  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
607  //x-
608  project(it->x-width, it->y-width, it->z+width, x1, y1);
609  project(it->x-width, it->y+width, it->z+width, x2, y2);
610  project(it->x-width, it->y+width, it->z-width, x3, y3);
611  project(it->x-width, it->y-width, it->z-width, x4, y4);
612  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
613  //y+
614  project(it->x-width, it->y+width, it->z+width, x1, y1);
615  project(it->x+width, it->y+width, it->z+width, x2, y2);
616  project(it->x+width, it->y+width, it->z-width, x3, y3);
617  project(it->x-width, it->y+width, it->z-width, x4, y4);
618  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
619  //y-
620  project(it->x-width, it->y-width, it->z+width, x1, y1);
621  project(it->x+width, it->y-width, it->z+width, x2, y2);
622  project(it->x+width, it->y-width, it->z-width, x3, y3);
623  project(it->x-width, it->y-width, it->z-width, x4, y4);
624  cairo_move_to (cr, x1, y1); cairo_line_to (cr, x2, y2); cairo_line_to (cr, x3, y3); cairo_line_to (cr, x4, y4); cairo_line_to (cr, x1, y1); cairo_close_path (cr);Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]=="SolidMode"?cairo_fill (cr):cairo_stroke (cr);
625 
626  cairo_restore (cr);
627  }
628  }
629  }
630  */
631 
632  if (type==CairoPNG)
633  cairo_surface_write_to_png (surface, filename);
634 
635  cairo_destroy (cr);
636  cairo_surface_destroy (surface);
637 }
638 
639 /*!
640  * \brief init function (should be in Constructor).
641  */
642 template < typename Space, typename KSpace>
643 void
644 DGtal::Board3DTo2D<Space, KSpace>::init()
645 {
646  Board3DTo2D<Space, KSpace>::createNewCubeList();
647 
648  std::vector<typename Board3DTo2D<Space, KSpace>::LineD3D> listeLine;
649  Board3DTo2D<Space, KSpace>::myLineSetList.push_back(listeLine);
650 
651  std::vector<typename Board3DTo2D<Space, KSpace>::BallD3D> listeBall;
652  Board3DTo2D<Space, KSpace>::myBallSetList.push_back(listeBall);
653 
654  Board3DTo2D<Space, KSpace>::myCurrentFillColor = DGtal::Color (220, 220, 220);
655  Board3DTo2D<Space, KSpace>::myCurrentLineColor = DGtal::Color (22, 22, 222, 50);
656 
657  Board3DTo2D<Space, KSpace>::myDefaultColor= DGtal::Color(255, 255, 255);
658 
659  //
660 
661  camera_position[0] = 5.000000; camera_position[1] = 5.000000; camera_position[2] = 29.893368;
662  camera_direction[0] = 0.000000; camera_direction[1] = 0.000000; camera_direction[2] = -1.000000;
663  camera_upVector[0] = 0.000000; camera_upVector[1] = 1.000000; camera_upVector[2] = 0.000000;
664 
665  ZNear = 0.001;
666  ZFar = 100.0;
667 
668  Board3DTo2D<Space, KSpace>::myModes ["Board3DTo2D"]="SolidMode";
669 }
670 
671 // //
672 ///////////////////////////////////////////////////////////////////////////////