Function main

examples/spherical_function3d.cpp:52–114

--------------------------- Main procedure --------------------------*/

  --------------------------*/
int main() {

  CImgList<unsigned char> object_colors;
  CImgList<float> object_opacities;

  // Define a 3D centered box.
  CImg<float> object_vertices = CImg<float>(3,8,1,1, // Define the 8 vertices of the cube
                                            -1,-1,-1, // (x0,y0,z0)
                                            1,-1,-1, // (x1,y1,z1)
                                            1,1,-1, // ...
                                            -1,1,-1,
                                            -1,-1,1,
                                            1,-1,1,
                                            1,1,1, // (x6,y6,z6)
                                            -1,1,1).transpose(); // (x7,y7,z7)

  // Define the 12 segments of the cube.
  CImgList<unsigned int> object_primitives =
    CImg<unsigned int>(1,12*2,1,1, 0,1, 1,2, 2,3, 3,0, 4,5, 5,6, 6,7, 7,4, 0,4, 1,5, 2,6, 3,7).get_split('y',12);

  object_colors.insert(object_primitives.size(),CImg<unsigned char>::vector(32,64,255));
  object_opacities.insert(object_primitives.size(),CImg<float>::vector(0.3f));

  // Define the spherical function's vertices.
  CImgList<float> spherical_vertices;
  const float a = 1;
  const unsigned int na = 132, nb = 132;
  for (unsigned int v = 0; v<na; ++v)
    for (unsigned int u = 0; u<nb; ++u) {
      const float
	alpha = (float)(u*2*cimg::PI/nb),
	beta = (float)(-cimg::PI/2 + v*cimg::PI/na),
	x = (float)((a*std::cos(beta))*std::cos(alpha)),
	y = (float)((a*std::cos(beta))*std::sin(alpha)),
	z = (float)(a*std::sin(beta)),
	altitude = cimg::abs(1 + 0.8f*std::cos(3*alpha)*std::sin(4*beta));
      spherical_vertices.insert(CImg<float>::vector(altitude*x,altitude*y,altitude*z));
    }

  // Define the spherical function's mesh.
  CImgList<unsigned int> spherical_primitives;
  for (unsigned int vv = 0; vv<na - 1; ++vv)
    for (unsigned int uu = 0; uu<nb; ++uu) {
      const int nv = (vv + 1)%na, nu = (uu + 1)%nb;
      spherical_primitives.insert(CImg<unsigned int>::vector(nb*vv + nu,nb*nv + uu,nb*vv + uu));
      spherical_primitives.insert(CImg<unsigned int>::vector(nb*vv + nu,nb*nv + nu,nb*nv + uu));
      object_colors.insert(CImg<>::vector(0,255,255));
      object_colors.insert(CImg<>::vector(100,200,255));
      object_opacities.insert(2,CImg<>::vector(1));
    }

  // Merge 3D objects together.
  object_vertices.append_object3d(object_primitives,spherical_vertices>'x',spherical_primitives);
  char title[4096] = { 0 };
  std::sprintf(title,"3D Spherical Function (%u vertices, %u primitives)",
               object_vertices.width(),object_primitives.size());
  CImgDisplay disp(640,480,title,0);
  CImg<unsigned char>(disp.width(),disp.height(),1,3,220).