目前，我们还在函数prepareModel中徘徊，因为这函数实在是太长了，近乎包含了整个数据处理过程。通过前面两篇，几何图形已经被导入到内存中，并且由一个simplemodel变成了一个optimizedmodel。现在程序拿起了屠刀，开始了切片过程。

　　通常来讲，在3D打印中所说的slicer过程包含了切割几何体，连接多边形，生成Gcode的所有过程，但在此程序中，slicer步仅为生成切片的切割操作。每一个切包含若干个闭合平面多边形。这些多边形的产生分两步。第一，所有的三角形面片被切割成为一条线段；第二，将每层的线段连接起来形成闭合多边形。前篇所述的点-面对应关系可以帮助这一过程快速进行，因为两个相邻的面上的切割线段很有可能就是相邻线段。这段程序还可以帮助修补几何图形中的空洞，因此STL模型不需要那么完美。它也可以处理错误的法向量，所以它可以翻转线段来满足端点到端点的连接。slicer步骤之后我们得到了一系列闭合多边形，然后会用到clipper来处理这些多边形（clipper只可以处理2D闭合多边形）。这一段来自cura-engine github主页的介绍。

　　切割过程代码在slicer.h和slicer.cpp这两个文件中。首先我们看其中的class slicer

1  SlicerSegment project2D(Point3& p0, Point3& p1, Point3& p2, int32_t z) const2     {3         SlicerSegment seg;4         seg.start.X = p0.x + int64_t(p1.x - p0.x) * int64_t(z - p0.z) / int64_t(p1.z - p0.z);5         seg.start.Y = p0.y + int64_t(p1.y - p0.y) * int64_t(z - p0.z) / int64_t(p1.z - p0.z);6         seg.end.X = p0.x + int64_t(p2.x - p0.x) * int64_t(z - p0.z) / int64_t(p2.z - p0.z);7         seg.end.Y = p0.y + int64_t(p2.y - p0.y) * int64_t(z - p0.z) / int64_t(p2.z - p0.z);8         return seg;9     }

　　class slicer中的这段代码是计算一个三角形p0,p1,p2和平面z的相交线的过程。就算公式也比较简单粗暴，如下

　　画得虽丑，但科学是美丽的,咳咳.。。。函数project2D所表达的就是上面的公式了，至于六种情况的分类，我们可以在构造函数slicer（）里面找到。没错，就是那六排if-else

Slicer(OptimizedVolume* ov, int32_t initial, int32_t thickness, bool keepNoneClosed, bool extensiveStitching)   把优化过的几何体ov按照厚度（thickness）切割成片片即layer。函数的前大半部分计算了三角面上被切割得到的线段，这些线段只有起点和终点，之间没有什么联系。后小部分使用makepolygon函数将线段们连接成多边形。

1  for(unsigned int layerNr=0; layerNr<layers.size(); layerNr++)2     {3         layers[layerNr].makePolygons(ov, keepNoneClosed, extensiveStitching);4     }

　　那么makepolygon函数又如何呢，这是一个200多行的大函数。本篇只能学习一小部分，欲窥其全貌，要等到之后我们学习了class Polygons和class Polygon这些类之后才能领悟。本篇先看这一段：

 1  Polygons openPolygonList; 2      3     for(unsigned int startSegment=0; startSegment < segmentList.size(); startSegment++) 4     { 5         if (segmentList[startSegment].addedToPolygon) 6             continue; 7          8         Polygon poly; 9         poly.add(segmentList[startSegment].start);10         11         unsigned int segmentIndex = startSegment;12         bool canClose;13         while(true)14         {15             canClose = false;16             segmentList[segmentIndex].addedToPolygon = true;17             Point p0 = segmentList[segmentIndex].end;18             poly.add(p0);19             int nextIndex = -1;20             OptimizedFace* face = &ov->faces[segmentList[segmentIndex].faceIndex];21             for(unsigned int i=0;i<3;i++)22             {23                 if (face->touching[i] > -1 && faceToSegmentIndex.find(face->touching[i]) != faceToSegmentIndex.end())24                 {25                     Point p1 = segmentList[faceToSegmentIndex[face->touching[i]]].start;26                     Point diff = p0 - p1;27                     if (shorterThen(diff, MM2INT(0.01)))28                     {29                         if (faceToSegmentIndex[face->touching[i]] == static_cast<int>(startSegment))30                             canClose = true;31                         if (segmentList[faceToSegmentIndex[face->touching[i]]].addedToPolygon)32                             continue;33                         nextIndex = faceToSegmentIndex[face->touching[i]];34                     }35                 }36             }37             if (nextIndex == -1)38                 break;39             segmentIndex = nextIndex;40         }41         if (canClose)42             polygonList.add(poly);43         else44             openPolygonList.add(poly);45     }

　　这一段讲述了把刚才求得的线段们连接成多边形的故事。可以看到，首先选取一个线段，然后通过它所在的三角面，找到和其相邻的三角面，在找到相邻三角面上的线段，查看其是否匹配，依次进行下去，最后会回到最初的那个三角面上，连接结束。

　　然而这种方法并不能完成所有的情况，比如几何体有一些缺陷什么的，后面的代码是针对各种情况所作的处理，这些目前还很难懂，因为在这之前要搞懂Polygon这个基础类到底是如何构建的。咱们下回分解。

延伸阅读

• cura-engine学习（1） 2017-06-05 16:17:39.0
• cura-engine学习（1） 2017-03-01 11:39:25.053

学习是年轻人改变自己的最好方式