首页 > 代码库 > cura-engine学习(3)
cura-engine学习(3)
目前,我们还在函数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) const 2 { 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 else 44 openPolygonList.add(poly); 45 }
这一段讲述了把刚才求得的线段们连接成多边形的故事。可以看到,首先选取一个线段,然后通过它所在的三角面,找到和其相邻的三角面,在找到相邻三角面上的线段,查看其是否匹配,依次进行下去,最后会回到最初的那个三角面上,连接结束。
然而这种方法并不能完成所有的情况,比如几何体有一些缺陷什么的,后面的代码是针对各种情况所作的处理,这些目前还很难懂,因为在这之前要搞懂Polygon这个基础类到底是如何构建的。咱们下回分解。
cura-engine学习(3)