引言
年初開始進(jìn)入3D打印行業(yè)，受命以Cura為基礎(chǔ)，研發(fā)一款自主的3D打印切片軟件。
自主研發(fā)要取其長處，補(bǔ)其不足，首先自然是要搞清楚Cura到底做了什么，讀Cura的代碼是必需的。我一向都覺得比起自己寫代碼來，讀別人的代碼是一個漫又而痛苦的過程，讀者的思想與寫者總有偏差，往往又無法驗證自己的猜想是否正確，只嘆人腦不是電腦，無法把眼前的代碼從頭到尾執(zhí)行一遍。不知道各位資深程序會有什么辦法，我的辦法是“翻譯”，看著別人寫的代碼，加上自己的理解之后，按自己的喜好重新寫出來，看一段翻譯一段，等全部翻譯完成，理論上作者的思路也明白了，同時還有了一份功能一模一樣的代碼，自己的理解是否正確，也可以通過執(zhí)行“翻譯”出來的代碼驗證。

計劃總是美好的，中間的工程確總是充滿變數(shù)，之間的曲折折疊不說。經(jīng)過若干次推倒重寫，勉強(qiáng)算是有了一份自己的切片軟件，又經(jīng)過了半年的推敲摸索以及打印經(jīng)驗積累，一個還算另自己滿意的切片軟件最終誕生。起名Pango，先觀大略。

Pango的介紹和說明先按下不表，以會有機(jī)會另外發(fā)文詳述。

在Pango的開發(fā)過程中，我對于Cura的理解也日益深入。時至今日，我也有了信心可以把我的這些經(jīng)驗、理解和心得分享出來，供大家參詳一二，若能對后來的Cura研究者有所助益，那是再好不過。

Cura的架構(gòu)

Cura是一個python語言實現(xiàn)，使用wxpython圖形界面框架的3D打印切片界面軟件，說它是界面軟件是因為Cura本身并不會進(jìn)行實際的切片操作。實際的切片工作是由另外一個C++語言實現(xiàn)的CuraEngine命令行軟件來具體負(fù)責(zé)的，用戶在Cura提供的界面上的絕大多數(shù)操作如加載模型、平穩(wěn)旋轉(zhuǎn)縮放、參數(shù)設(shè)置等最終會轉(zhuǎn)換成并執(zhí)行一條CuraEngine命令；CuraEngine把輸入的STL、DAE或OBJ模型文件切片輸出成gcode字符串返回給Cura；Cura再把gcode在3D界面上可視化成路徑展現(xiàn)給用戶。

我主要參考的代碼是CuraEngine，本文主要篇幅也會放在CuraEngine上。而Pango的界面代碼就主要靠我自己發(fā)揮了。

Cura和CuraEngine都可以Github上找到，地址：
https://github.com/daid/Cura
https://github.com/Ultimaker/CuraEngine

我所參考的版本是15.04，15.06之后Cura和CuraEngine都有較大的改動，但核心思想沒變。所以本文分析的代碼也到15.04為止。

言歸正傳，下面我們將開始一步一步揭開CuraEngine把一個模型文件轉(zhuǎn)換成為gcode的過程。
切片流程概述
從總體上講，CuraEngine的切片分為五個步驟：

步驟一：模型載入
有一點(diǎn)3D編程經(jīng)驗的人都知道，計算機(jī)中的3D模型大多是以三角形面組合成的表面所包裹的空間來表示的。三角形作為3D模型的基本單元，有結(jié)構(gòu)簡單，通用性強(qiáng)，可組合成任意面的特點(diǎn)；空間坐標(biāo)中只要三個點(diǎn)就可以表示一個唯一的三角形，兩點(diǎn)只能表示一條直線，而再多的直線也無法組成一個平面；空間中的任意三個不共線的點(diǎn)都可以組成一個三角形，而四個點(diǎn)所組成的四邊形就必需要求四點(diǎn)共面；任意的表面都可以拆解成三角形，一個四邊形可以拆解成兩個三角形，但一個三角形確沒有辦法用四邊形組合而成。計算機(jī)所擅長的事情就是把簡單的事情不斷重復(fù)，而三角形正是因為這些特性，成為了計算機(jī)3D世界的基石。

CuraEngine內(nèi)部也是用三角形組合來表示模型的，不過同樣一個三角形組合，確有無窮多種數(shù)據(jù)結(jié)構(gòu)來進(jìn)行存儲，CuraEngine切片的第一步，就是從外部讀入模型數(shù)據(jù)，轉(zhuǎn)換成以CuraEngine內(nèi)部的數(shù)據(jù)結(jié)構(gòu)所表示的三角形組合。

有了三角形組合還不夠，CuraEngine在載入模型階段還要對三角形進(jìn)行關(guān)聯(lián)。兩個三角形共有一條邊的，就可以判斷它們?yōu)橄噜徣切?，一個三角形有三條邊，所以最多可以有三個相鄰三角形，一般而言，如果模型是封閉的，那它的每一個三角形都會有三個相鄰三角形。

有了三角形的相鄰關(guān)系，就可以大幅提高下一個步驟分層過程的處理速度。Cura之所以成為當(dāng)前市場切片速度最快的軟件，這是其中最顯著的優(yōu)化之一。

模型載入更詳細(xì)的過程會另文分析，敬請期待。

步驟二：分層
如果把模型放在XY平面上，Z軸對應(yīng)的就是模型高度。我們把XY平面抬高一定高度再與模型的表面相交，就可以得到模型在這個高度上層切片。所謂的分層就是每隔一定高度就用一個XY平面去和模型相交作層切片，全部切完后就可以得到模型在每一個高度上的輪廓線。就像是切土豆片一樣，把一個圓的或不圓異或不管什么奇形怪狀的土豆用菜刀一刀一刀切開，最后就能得到一盤薄如紙片的土豆片，當(dāng)然那還得你的刀功要足夠好才行。

分層本質(zhì)上就是一個把3D模型轉(zhuǎn)化為一系列2D平面的過程，自此之后的所有操作就都是在2D圖形的基礎(chǔ)上進(jìn)行了。

在前面模型載入階段我說到了CuraEngine埋了一個三角形關(guān)聯(lián)的伏筆，作用是什么，現(xiàn)在就可以揭曉了。我們知道，兩個平面相交，得到的是一條直線，一個平面和一個三角形相交，就得到一條線段。當(dāng)然也有可能什么也得不到，平臺平行啦，三角形的三個點(diǎn)都在平面的同一面之類，這些我們可以不管，我們現(xiàn)在只關(guān)心和平面有交集的那些三角形即可。我們把一個平面和所有的三角形都相交了一遍，得到了許許多多的線段，但我們需要的是2D圖形，三角形是2D圖形，四邊形，任意多邊形都是2D圖形，但線段不是。所以我們就要把這些線段試著連成一個多邊形，那么問題來了，要把這些線段連起來，只能兩個兩個地去試，看看它們是不是共端點(diǎn)了，粗算一下，每一層都是平方級的復(fù)雜度，再算上層數(shù)，那就是三次方級了。但現(xiàn)在，我們知道了三角形的關(guān)聯(lián)關(guān)系，兩個關(guān)聯(lián)的三角形，如果都與一個平面相交，那它們的交線一定也是關(guān)聯(lián)的，這一下，每一條線段只需要判斷三它與它相鄰三角形，看看與這個平面有沒有交線即可，一下子就把問題的復(fù)雜度降了一個次元。速度自然可以有質(zhì)的提升。

分層更詳細(xì)的過程會另文分析，敬請期待。

步驟三：劃分組件

經(jīng)過分層之后，我們得到了一疊2D平面圖形。接下來需要做的事情就是對每一層的平面圖形進(jìn)行跑馬圈地，標(biāo)記出哪里是外墻、內(nèi)墻、填充、上下表面、支撐等等。

3D打印在每一層是以組件單位，所謂組件指的就是每一層2D平面圖形里可以連通的區(qū)域，比如左圖就可以拆分為黃綠藍(lán)三個組件。而打印的順序就每打印完一個組件，接著會挑選一個離上一個組件最近的組件作為下一個進(jìn)行打印，如此循環(huán)直至一層的組件全部打印完成；接著會Z軸上升，重復(fù)上述步驟打印下一層的所有組件。

至于每一個組件怎么打印，就和我們手工畫畫一樣，先打邊線再對邊線內(nèi)部填充。邊線可以打印多層，最外層的邊線稱為外墻，其它的統(tǒng)稱為內(nèi)墻，CuraEngine之所以要對內(nèi)外墻進(jìn)行區(qū)分，是為了可以為它們定制不同的打印參數(shù)：外墻會被人觀察到，所以可以采用低速以提高表面質(zhì)量，內(nèi)墻只是起增加強(qiáng)度的作用，可以稍稍加快打印速度以節(jié)省時間。這些都可以在Cura界面的高級選項里進(jìn)行配置。

有一點(diǎn)值得注意的是，也是我半年打印的經(jīng)驗，由于FDM擠出裝置的特性所至，擠出機(jī)的擠出通過影響的只是加熱腔里的熔絲壓力，間接決定了噴頭的擠出速度，而加熱腔本身對于壓力就有一個緩沖作用，所以擠出速度的突變并不會使得噴頭的擠出速度立即跟著變化，而是有一個延遲，這一點(diǎn)在遠(yuǎn)端送絲的機(jī)器上更為明顯。而恰恰我們公司的主打產(chǎn)品F3CL就是遠(yuǎn)端送絲，在Pango中考慮到這個問題，并加上了特殊處理，事實證明的確對打印質(zhì)量有一定的提升。具體辦法是什么，我先賣個關(guān)子，會Pango的專文里進(jìn)行講解。

內(nèi)外墻標(biāo)記完之后就是填充和上下表面的標(biāo)記了，填充有一個填充率，0%填充率就是無填充，100%就是打成一個密實的平面，所以上下表面就是填充率為100%的填充。中間的填充率自然介于兩者之間，就像一張漁網(wǎng)，填充率越高網(wǎng)眼越細(xì)。

軟件會先把內(nèi)墻里面的部分統(tǒng)統(tǒng)標(biāo)記成填充，之后再進(jìn)一步判斷其中有哪些部分要轉(zhuǎn)換成為上下表面。至于是哪些部分，在設(shè)置里會有一個上下表面層數(shù)的設(shè)置，它代表了模型的上下與空氣接觸的表面有幾層，這里就會用到這個參數(shù)，CuraEngine會把當(dāng)前層上下n層（上下表面層數(shù)）取出來與當(dāng)前層進(jìn)行比較，凡是當(dāng)前層有而上下n層沒有的部分就會被劃歸到表皮。而原來的填充區(qū)域在割除被劃到表皮的部分后剩下的部分就是最終的填充區(qū)域。

CuraEngine在處理過程中大量用到了2D圖形運(yùn)算操作，有關(guān)2D圖形的運(yùn)算，有很多人研究，也被做成許多成熟的庫可以調(diào)用。CuraEngine的作者拿來主義，選取了一個他認(rèn)為比較好用的庫，叫ClipperLib的庫直接內(nèi)嵌到軟件之中，ClipperLib所使用的2D圖形算法也很著名，叫Vatti's clipping algorithm，很復(fù)雜，我也沒有完全搞懂，有興趣的讀者要是搞懂了可以多多交流。ClipperLib的網(wǎng)址是：http://www.angusj.com/delphi/clipper.php

這里我先簡單介紹一下CuraEngine所用到的幾種2D圖形的運(yùn)算，都是由ClipperLib實現(xiàn)的：交、并、差、偏移。與集合操作類似先看圖：

圖形相交
二元圖形操作，最終結(jié)果為兩個圖形共同包含的區(qū)域。記作：A * B
圖形相并
二元圖形操作，最終結(jié)果為兩個圖形其中的一個或兩者所包含的區(qū)域。記作：A + B
圖形相減
二元圖形操作，最終結(jié)果為屬于前者但不屬于后者的區(qū)域。記作：A - B
圖形偏移（外擴(kuò)）
一元圖形操作，最終結(jié)果為圖形區(qū)域的邊界向外擴(kuò)展指定的距離。
圖形偏移（內(nèi)縮）
一元圖形操作，最終結(jié)果為圖形區(qū)域的邊界向內(nèi)收縮指定的距離。內(nèi)縮與外擴(kuò)互為逆運(yùn)算。
這些就是CuraEngine所用到的2D圖形操作，運(yùn)算不多，確可以做許許多多的事情，比如上面所說的上下表面計算，就可以用數(shù)學(xué)公式來表示：


表面(i) = [填充(i) - 層(i + n)] + [填充(i) - 層(i - n)]




填充(i) = 填充(i) - 表面(i)


其中，i為當(dāng)前層號，n為上下表面層數(shù)（可以不一樣）。多簡單，數(shù)學(xué)就是這么任性！
同樣的，組件里面內(nèi)外墻，填充怎么劃分，只用一個內(nèi)縮運(yùn)算就可以搞定：


外墻 = 組件.offset(-線寬)




內(nèi)墻1 = 組件.offset(-線寬 * 2)

...

內(nèi)墻n = 組件.offset(-線寬 * (n + 1))




填充 = 組件.offset(-線寬 * (n + 2))



如果模型無需支撐，那組件劃分到這里就可以收工。否則，接下就是計算支撐的時間了。

我用CuraEngine半年下來覺得它最大的不足就是在支撐上，這也是我在Pango投入最大精力要改進(jìn)的地方，這里就先簡單介紹一下CuraEngine所用的支撐算法。

CuraEngine首先把整個打印空間在XY平臺上劃分成為50um*50um的網(wǎng)格，每個網(wǎng)格的中心點(diǎn)再延Z軸向上作一條直線，這條直線可能會與組成3D模型的三角形相交，三角形與直線的交點(diǎn)以及這個三角形的傾斜度會被記錄到網(wǎng)格里面。

現(xiàn)在每個網(wǎng)格里會記錄下一串被稱為支撐點(diǎn)的列表，每個支撐點(diǎn)包含一個高度和一個傾斜度信息。接下來會對每個網(wǎng)格的支撐點(diǎn)列表按照高度從低到高排序。根據(jù)這些信息就可以判斷模型上任意一個點(diǎn)是否需要支撐了，怎么判斷，我們看圖說話：

讓我們從底面開始延著一條網(wǎng)格中心往上走，起始我們是在模型外部，當(dāng)遇到第一個支撐點(diǎn)的時候，就從模型外部進(jìn)行了模型內(nèi)部，我們稱這個支撐點(diǎn)為進(jìn)點(diǎn)。

繼續(xù)向上，遇到了第二個支撐點(diǎn)，從模型內(nèi)部又退到了模型外部，我們稱這個支撐點(diǎn)為出點(diǎn)。

接著向上，我們可以發(fā)現(xiàn)，進(jìn)點(diǎn)與出點(diǎn)總是交替出現(xiàn)的。

利用這個規(guī)律，對于模型上任何一個點(diǎn)，我們只要找到這個點(diǎn)所對應(yīng)的網(wǎng)格，再找到這個網(wǎng)格里在這個點(diǎn)以上最近的一個支撐點(diǎn)，我們就可以得到兩個信息：這個點(diǎn)之上是否有模型懸空；這個點(diǎn)上面的懸空點(diǎn)的面的傾斜度是多少。

Cura界面的專家設(shè)置里面有支撐角度的設(shè)置，如果一個點(diǎn)處于模型懸空部分以下，并且懸空點(diǎn)傾斜度大于支撐角度，那這個點(diǎn)就是需要支撐的。所一個平臺上所有的需要支撐的點(diǎn)連接起來圍成的2D圖形就是支撐區(qū)域。

CuraEngine所使用的支撐算法比較粗糙，但勝在速度很快。先不說網(wǎng)格化后失去了精度，通過傾斜角度來判斷，模型下方一旦傾斜角發(fā)生了突變，像左圖這種從負(fù)45
度一下突變成正45度，傾斜角判斷無能為力，除非把它改大到60度，這樣的話，整個模型都會被過度支撐。這樣矯枉過正，既不科學(xué)，也浪費(fèi)材料和打印時間，還會對模型表面質(zhì)量帶來不好的影響。

科學(xué)的支撐算法應(yīng)該是找到模型局部最低點(diǎn)進(jìn)行支撐，最低點(diǎn)以上不一定需要支撐。因為FDM材料本身的粘性，使得材料的走線可以有一部分懸空而不坍塌，這個效果被稱為Overhang，只要上層材料的懸空距離小于一定的值，它就不需要支撐，這個距離以我的經(jīng)驗應(yīng)該在1/4到1/2線寬之間。我在Pango中就基于這個思路重新實現(xiàn)了支撐的算法，結(jié)果雖然

速度不如Cura的支撐算法那么快，但效果非常好，該撐的地方撐，不該撐的地方也不會多此一舉。

Pango的支撐算法我會在以后專文介紹。順帶一說，CuraEngine在下半年做了很大的改動，其中之一就是拋棄了之前的支撐算法，而新的算法也和我上面所講的思想異曲同工。我要聲明的是Pango的支撐算法和CuraEngine誰也沒有抄誰，我的算法是自己拍腦袋想出來的。算是英雄所見略同吧。

支撐范圍確定之后，也和組件一樣，可以有外墻、內(nèi)墻、填充、表面。依樣畫葫蘆即可。CuraEngine對于支撐，只會生成外墻和填充，Pango更多。

組件和支撐就是CuraEngine在這一步所生成的結(jié)果，這一步可以說是整個切片過程的核心，更詳細(xì)的過程會另文分析，敬請期待。

步驟四：路徑生成

地圈好了，就該在里面種菜了。這一步路徑生成就要開始規(guī)劃噴頭在不同的組件中怎么運(yùn)動。路徑按大類來分，有輪廓和填充兩種。

輪廓很簡單，沿著2D圖形的邊線走一圈即可。前一步所生成的外墻、內(nèi)墻都屬于輪廓，可以直接把它們的圖形以設(shè)置里的線寬轉(zhuǎn)換為輪廓路徑。

填充稍微要復(fù)雜一些，2D圖形指定的只是填充的邊界，而生成的路徑則是在邊界的范圍內(nèi)的條紋或網(wǎng)格結(jié)構(gòu)，就像窗簾或者漁網(wǎng)，如左圖。這兩種就最基本的結(jié)構(gòu)，當(dāng)然也許你還可以想出其它花式的填充，比如蜂窩狀或者S型，這些在新的Cura或者別的什么切片軟件里可能會實現(xiàn)，但我打印下來還是這兩種基本結(jié)構(gòu)更讓人放心。

CuraEngine在專家設(shè)置里可以對填充類型進(jìn)行選擇，里面除了條紋和網(wǎng)格外還有一個自動選項，默認(rèn)就是自動。自動模式會根據(jù)當(dāng)前的填充率進(jìn)行切換，當(dāng)填充率小于20%就用條紋填充，否則使用網(wǎng)格填充。因為網(wǎng)格結(jié)構(gòu)雖然更為合理，但它有一個問題，就是交點(diǎn)的地方會打兩次。填充率越高，交點(diǎn)越密，對打印質(zhì)量的影響會越大。我們知道，表面就是100%的填充，如果表面用網(wǎng)格打，不但無法打密實，表面還會坑坑洼洼，所以100%填充只能用條紋打，這就是CuraEngine推薦自動模式的原因。

至于填充率，就反映在線與線的間距上。100%填充率間距為0；0%填充率間距無限大，一根線條也不會有。

每個組件獨(dú)立的路徑生成好了，還要確定打印的先后順序。順序先好了可以少走彎路，打印速度和質(zhì)量都會有提升。路徑的順序以先近后遠(yuǎn)為基本原則：每打印完一條路徑，當(dāng)前位置是上一條路徑的終點(diǎn)；在當(dāng)前層里剩下還沒打印的路徑中挑選一條起點(diǎn)離當(dāng)前位置最近的一條路徑開打。路徑的起點(diǎn)可以是路徑中的任意一個點(diǎn)，程序會自行判斷。而路徑的終點(diǎn)有兩種可能：對于直線，圖形只有兩個點(diǎn)，終點(diǎn)就是除起點(diǎn)之外的那個點(diǎn)；對于輪廓，終點(diǎn)就是起點(diǎn)，因為輪廓是一個封閉圖形，從它的起點(diǎn)開始沿任意方向走一圈，最后還會回到起點(diǎn)。CuraEngine對路徑選擇做了一個估值，除了考慮到先近后遠(yuǎn)外，還順便參考了下一個點(diǎn)相對于當(dāng)前點(diǎn)的方向，它的物理意義就是減少噴頭轉(zhuǎn)彎。賽車在直道上開得一定比彎道快，不是么。

路徑的順序也確定了，還有一個問題需要考慮：如果前后兩條路徑首尾相連，那直接走就是了，但大多數(shù)情況并非如此，前一條路徑的終點(diǎn)往往和后一條路徑起點(diǎn)之間有一段距離。這時候去往下一點(diǎn)的路上要小心了，肯定不能繼續(xù)擠出材料，否則輕則拉絲，重則模型面目全非。這段路噴頭就需要空走，即噴頭只移動，不吐絲，那只要把擠出機(jī)停下來不轉(zhuǎn)就行了嗎？也不行，因為前面分析過，擠出機(jī)的速度要傳導(dǎo)到噴嘴，有一個延遲，不是你說停它就立即停下來的。這是FDM打印的通病，解決辦法就是回抽。所謂回抽，就是在空走之前先讓擠出機(jī)高速反轉(zhuǎn)一段材料，這樣就可以瞬間把加熱腔里的材料抽光，再移動過去，中間就不會擠出材料，到了下一個點(diǎn)，在打印之前，先把剛才抽回去的絲再按一樣的長度放回來，繼續(xù)打印?；爻榭梢院芎玫亟鉀Q空走拉絲的問題，但是它很慢，以抽一次0.5秒來算的話，如果打印一個表面，0.4線寬，10厘米的距離至少回抽25下，10幾秒鐘的時間一層，幾百上千層打下來，光回抽所用的時間就是幾個小時，是可忍孰不可忍！

CuraEngine給我們提供了解決方案就是Comb，也就是繞路。我們先來看，是不是所有的回抽都是必需的呢？不回抽會拉絲是肯定的，但如果需要空走的路徑本來就要打印的，那拉絲又有何妨。按這個思路，就可以給每個組件設(shè)定一個邊界，只要路徑的起點(diǎn)和終點(diǎn)都在這個邊界之內(nèi)的，空走都不回抽。這樣可以解決80%的問題，但如果是左圖這樣的情況就行不通。

紅色是起點(diǎn)，綠色是終點(diǎn)，直接走過去會走出邊界的范圍。這時我們就要繞一點(diǎn)路，走一條曲線到達(dá)我們的目的地。這就是Comb所做的事情，在Cura專家設(shè)置里面可以對Comb進(jìn)行設(shè)置，選擇開啟、關(guān)閉還有表面不Comb。Comb可以大幅節(jié)省打印時間，但是同一個地方打印多次對模型質(zhì)量還是會有細(xì)微的影響，個中利弊，交給用戶自己判斷。

Comb的調(diào)整是個細(xì)致活，Pango花了相當(dāng)多的時間來微調(diào)Comb功能以求達(dá)到更好的效果，過程繁瑣，不再贅述。

至此路徑生成完成，更詳細(xì)的過程另文分析，敬請期待。

步驟五：gcode生成
路徑都生成好了，還需要翻譯對打印機(jī)可以實別的gcode代碼才行。這一步花樣不多，按部就班即可。

先讓打印機(jī)做一些準(zhǔn)備工作：歸零、加熱噴頭和平臺、抬高噴頭、擠一小段絲、風(fēng)扇設(shè)置。

從下到上一層一層打印，每層打印之前先用G0抬高Z坐標(biāo)到相應(yīng)位置。

按照路徑，每個點(diǎn)生成一條gcode。其中空走G0；邊擠邊走用G1，Cura的設(shè)置里有絲材的直徑、線寬，可以算出走這些距離需要擠出多少材料；G0和G1的速度也都在設(shè)置里可以調(diào)整。

若需回抽，用G1生成一條E軸倒退的代碼。在下一條G1執(zhí)行之前，再用G1生成一條相應(yīng)的E軸前進(jìn)的代碼。

所有層都打完后讓打印機(jī)做一些收尾工作：關(guān)閉加熱、XY歸零、電機(jī)釋放。

生成gcode的過程中，CuraEngine也會模擬一遍打印過程，用來計算出打印所需要的時間和材料長度，這些也會寫在gcode的注釋里供用戶參考。

gcode生成不用另文詳細(xì)分析，但是gcode的說明還是可以專文分析一下，敬請期待。

待續(xù)

寫了這么多，Cura的切片流程也只能講個大概，也算是個提綱，希望對大家有所幫助。我計劃對于上面的第一個步驟再專文分析。除此之外，還有Cura界面部分以及Cura與CuraEngine的通訊也可以講講。之后就是我半年創(chuàng)作，自我感覺良好到覺得可以超越Cura的Pango，也是不說不快的。

引言
年初開始進(jìn)入3D打印行業(yè)，受命以Cura為基礎(chǔ)，研發(fā)一款自主的3D打印切片軟件。
自主研發(fā)要取其長處，補(bǔ)其不足，首先自然是要搞清楚Cura到底做了什么，讀Cura的代碼是必需的。我一向都覺得比起自己寫代碼來，讀別人的代碼是一個漫又而痛苦的過程，讀者的思想與寫者總有偏差，往往又無法驗證自己的猜想是否正確，只嘆人腦不是電腦，無法把眼前的代碼從頭到尾執(zhí)行一遍。不知道各位資深程序會有什么辦法，我的辦法是“翻譯”，看著別人寫的代碼，加上自己的理解之后，按自己的喜好重新寫出來，看一段翻譯一段，等全部翻譯完成，理論上作者的思路也明白了，同時還有了一份功能一模一樣的代碼，自己的理解是否正確，也可以通過執(zhí)行“翻譯”出來的代碼驗證。

計劃總是美好的，中間的工程確總是充滿變數(shù)，之間的曲折折疊不說。經(jīng)過若干次推倒重寫，勉強(qiáng)算是有了一份自己的切片軟件，又經(jīng)過了半年的推敲摸索以及打印經(jīng)驗積累，一個還算另自己滿意的切片軟件最終誕生。起名Pango，先觀大略。

Pango的介紹和說明先按下不表，以會有機(jī)會另外發(fā)文詳述。

在Pango的開發(fā)過程中，我對于Cura的理解也日益深入。時至今日，我也有了信心可以把我的這些經(jīng)驗、理解和心得分享出來，供大家參詳一二，若能對后來的Cura研究者有所助益，那是再好不過。

Cura的架構(gòu)

Cura是一個python語言實現(xiàn)，使用wxpython圖形界面框架的3D打印切片界面軟件，說它是界面軟件是因為Cura本身并不會進(jìn)行實際的切片操作。實際的切片工作是由另外一個C++語言實現(xiàn)的CuraEngine命令行軟件來具體負(fù)責(zé)的，用戶在Cura提供的界面上的絕大多數(shù)操作如加載模型、平穩(wěn)旋轉(zhuǎn)縮放、參數(shù)設(shè)置等最終會轉(zhuǎn)換成并執(zhí)行一條CuraEngine命令；CuraEngine把輸入的STL、DAE或OBJ模型文件切片輸出成gcode字符串返回給Cura；Cura再把gcode在3D界面上可視化成路徑展現(xiàn)給用戶。

我主要參考的代碼是CuraEngine，本文主要篇幅也會放在CuraEngine上。而Pango的界面代碼就主要靠我自己發(fā)揮了。

Cura和CuraEngine都可以Github上找到，地址：
https://github.com/daid/Cura
https://github.com/Ultimaker/CuraEngine

我所參考的版本是15.04，15.06之后Cura和CuraEngine都有較大的改動，但核心思想沒變。所以本文分析的代碼也到15.04為止。

言歸正傳，下面我們將開始一步一步揭開CuraEngine把一個模型文件轉(zhuǎn)換成為gcode的過程。
切片流程概述
從總體上講，CuraEngine的切片分為五個步驟：

步驟一：模型載入
有一點(diǎn)3D編程經(jīng)驗的人都知道，計算機(jī)中的3D模型大多是以三角形面組合成的表面所包裹的空間來表示的。三角形作為3D模型的基本單元，有結(jié)構(gòu)簡單，通用性強(qiáng)，可組合成任意面的特點(diǎn)；空間坐標(biāo)中只要三個點(diǎn)就可以表示一個唯一的三角形，兩點(diǎn)只能表示一條直線，而再多的直線也無法組成一個平面；空間中的任意三個不共線的點(diǎn)都可以組成一個三角形，而四個點(diǎn)所組成的四邊形就必需要求四點(diǎn)共面；任意的表面都可以拆解成三角形，一個四邊形可以拆解成兩個三角形，但一個三角形確沒有辦法用四邊形組合而成。計算機(jī)所擅長的事情就是把簡單的事情不斷重復(fù)，而三角形正是因為這些特性，成為了計算機(jī)3D世界的基石。

CuraEngine內(nèi)部也是用三角形組合來表示模型的，不過同樣一個三角形組合，確有無窮多種數(shù)據(jù)結(jié)構(gòu)來進(jìn)行存儲，CuraEngine切片的第一步，就是從外部讀入模型數(shù)據(jù)，轉(zhuǎn)換成以CuraEngine內(nèi)部的數(shù)據(jù)結(jié)構(gòu)所表示的三角形組合。

有了三角形組合還不夠，CuraEngine在載入模型階段還要對三角形進(jìn)行關(guān)聯(lián)。兩個三角形共有一條邊的，就可以判斷它們?yōu)橄噜徣切?，一個三角形有三條邊，所以最多可以有三個相鄰三角形，一般而言，如果模型是封閉的，那它的每一個三角形都會有三個相鄰三角形。

有了三角形的相鄰關(guān)系，就可以大幅提高下一個步驟分層過程的處理速度。Cura之所以成為當(dāng)前市場切片速度最快的軟件，這是其中最顯著的優(yōu)化之一。

模型載入更詳細(xì)的過程會另文分析，敬請期待。

步驟二：分層
如果把模型放在XY平面上，Z軸對應(yīng)的就是模型高度。我們把XY平面抬高一定高度再與模型的表面相交，就可以得到模型在這個高度上層切片。所謂的分層就是每隔一定高度就用一個XY平面去和模型相交作層切片，全部切完后就可以得到模型在每一個高度上的輪廓線。就像是切土豆片一樣，把一個圓的或不圓異或不管什么奇形怪狀的土豆用菜刀一刀一刀切開，最后就能得到一盤薄如紙片的土豆片，當(dāng)然那還得你的刀功要足夠好才行。

分層本質(zhì)上就是一個把3D模型轉(zhuǎn)化為一系列2D平面的過程，自此之后的所有操作就都是在2D圖形的基礎(chǔ)上進(jìn)行了。

在前面模型載入階段我說到了CuraEngine埋了一個三角形關(guān)聯(lián)的伏筆，作用是什么，現(xiàn)在就可以揭曉了。我們知道，兩個平面相交，得到的是一條直線，一個平面和一個三角形相交，就得到一條線段。當(dāng)然也有可能什么也得不到，平臺平行啦，三角形的三個點(diǎn)都在平面的同一面之類，這些我們可以不管，我們現(xiàn)在只關(guān)心和平面有交集的那些三角形即可。我們把一個平面和所有的三角形都相交了一遍，得到了許許多多的線段，但我們需要的是2D圖形，三角形是2D圖形，四邊形，任意多邊形都是2D圖形，但線段不是。所以我們就要把這些線段試著連成一個多邊形，那么問題來了，要把這些線段連起來，只能兩個兩個地去試，看看它們是不是共端點(diǎn)了，粗算一下，每一層都是平方級的復(fù)雜度，再算上層數(shù)，那就是三次方級了。但現(xiàn)在，我們知道了三角形的關(guān)聯(lián)關(guān)系，兩個關(guān)聯(lián)的三角形，如果都與一個平面相交，那它們的交線一定也是關(guān)聯(lián)的，這一下，每一條線段只需要判斷三它與它相鄰三角形，看看與這個平面有沒有交線即可，一下子就把問題的復(fù)雜度降了一個次元。速度自然可以有質(zhì)的提升。

分層更詳細(xì)的過程會另文分析，敬請期待。

步驟三：劃分組件

經(jīng)過分層之后，我們得到了一疊2D平面圖形。接下來需要做的事情就是對每一層的平面圖形進(jìn)行跑馬圈地，標(biāo)記出哪里是外墻、內(nèi)墻、填充、上下表面、支撐等等。

3D打印在每一層是以組件單位，所謂組件指的就是每一層2D平面圖形里可以連通的區(qū)域，比如左圖就可以拆分為黃綠藍(lán)三個組件。而打印的順序就每打印完一個組件，接著會挑選一個離上一個組件最近的組件作為下一個進(jìn)行打印，如此循環(huán)直至一層的組件全部打印完成；接著會Z軸上升，重復(fù)上述步驟打印下一層的所有組件。

至于每一個組件怎么打印，就和我們手工畫畫一樣，先打邊線再對邊線內(nèi)部填充。邊線可以打印多層，最外層的邊線稱為外墻，其它的統(tǒng)稱為內(nèi)墻，CuraEngine之所以要對內(nèi)外墻進(jìn)行區(qū)分，是為了可以為它們定制不同的打印參數(shù)：外墻會被人觀察到，所以可以采用低速以提高表面質(zhì)量，內(nèi)墻只是起增加強(qiáng)度的作用，可以稍稍加快打印速度以節(jié)省時間。這些都可以在Cura界面的高級選項里進(jìn)行配置。

有一點(diǎn)值得注意的是，也是我半年打印的經(jīng)驗，由于FDM擠出裝置的特性所至，擠出機(jī)的擠出通過影響的只是加熱腔里的熔絲壓力，間接決定了噴頭的擠出速度，而加熱腔本身對于壓力就有一個緩沖作用，所以擠出速度的突變并不會使得噴頭的擠出速度立即跟著變化，而是有一個延遲，這一點(diǎn)在遠(yuǎn)端送絲的機(jī)器上更為明顯。而恰恰我們公司的主打產(chǎn)品F3CL就是遠(yuǎn)端送絲，在Pango中考慮到這個問題，并加上了特殊處理，事實證明的確對打印質(zhì)量有一定的提升。具體辦法是什么，我先賣個關(guān)子，會Pango的專文里進(jìn)行講解。

內(nèi)外墻標(biāo)記完之后就是填充和上下表面的標(biāo)記了，填充有一個填充率，0%填充率就是無填充，100%就是打成一個密實的平面，所以上下表面就是填充率為100%的填充。中間的填充率自然介于兩者之間，就像一張漁網(wǎng)，填充率越高網(wǎng)眼越細(xì)。

軟件會先把內(nèi)墻里面的部分統(tǒng)統(tǒng)標(biāo)記成填充，之后再進(jìn)一步判斷其中有哪些部分要轉(zhuǎn)換成為上下表面。至于是哪些部分，在設(shè)置里會有一個上下表面層數(shù)的設(shè)置，它代表了模型的上下與空氣接觸的表面有幾層，這里就會用到這個參數(shù)，CuraEngine會把當(dāng)前層上下n層（上下表面層數(shù)）取出來與當(dāng)前層進(jìn)行比較，凡是當(dāng)前層有而上下n層沒有的部分就會被劃歸到表皮。而原來的填充區(qū)域在割除被劃到表皮的部分后剩下的部分就是最終的填充區(qū)域。

CuraEngine在處理過程中大量用到了2D圖形運(yùn)算操作，有關(guān)2D圖形的運(yùn)算，有很多人研究，也被做成許多成熟的庫可以調(diào)用。CuraEngine的作者拿來主義，選取了一個他認(rèn)為比較好用的庫，叫ClipperLib的庫直接內(nèi)嵌到軟件之中，ClipperLib所使用的2D圖形算法也很著名，叫Vatti's clipping algorithm，很復(fù)雜，我也沒有完全搞懂，有興趣的讀者要是搞懂了可以多多交流。ClipperLib的網(wǎng)址是：http://www.angusj.com/delphi/clipper.php

這里我先簡單介紹一下CuraEngine所用到的幾種2D圖形的運(yùn)算，都是由ClipperLib實現(xiàn)的：交、并、差、偏移。與集合操作類似先看圖：

圖形相交
二元圖形操作，最終結(jié)果為兩個圖形共同包含的區(qū)域。記作：A * B
圖形相并
二元圖形操作，最終結(jié)果為兩個圖形其中的一個或兩者所包含的區(qū)域。記作：A + B
圖形相減
二元圖形操作，最終結(jié)果為屬于前者但不屬于后者的區(qū)域。記作：A - B
圖形偏移（外擴(kuò)）
一元圖形操作，最終結(jié)果為圖形區(qū)域的邊界向外擴(kuò)展指定的距離。
圖形偏移（內(nèi)縮）
一元圖形操作，最終結(jié)果為圖形區(qū)域的邊界向內(nèi)收縮指定的距離。內(nèi)縮與外擴(kuò)互為逆運(yùn)算。
這些就是CuraEngine所用到的2D圖形操作，運(yùn)算不多，確可以做許許多多的事情，比如上面所說的上下表面計算，就可以用數(shù)學(xué)公式來表示：


表面(i) = [填充(i) - 層(i + n)] + [填充(i) - 層(i - n)]




填充(i) = 填充(i) - 表面(i)


其中，i為當(dāng)前層號，n為上下表面層數(shù)（可以不一樣）。多簡單，數(shù)學(xué)就是這么任性！
同樣的，組件里面內(nèi)外墻，填充怎么劃分，只用一個內(nèi)縮運(yùn)算就可以搞定：


外墻 = 組件.offset(-線寬)




內(nèi)墻1 = 組件.offset(-線寬 * 2)

...

內(nèi)墻n = 組件.offset(-線寬 * (n + 1))




填充 = 組件.offset(-線寬 * (n + 2))



如果模型無需支撐，那組件劃分到這里就可以收工。否則，接下就是計算支撐的時間了。

我用CuraEngine半年下來覺得它最大的不足就是在支撐上，這也是我在Pango投入最大精力要改進(jìn)的地方，這里就先簡單介紹一下CuraEngine所用的支撐算法。

CuraEngine首先把整個打印空間在XY平臺上劃分成為50um*50um的網(wǎng)格，每個網(wǎng)格的中心點(diǎn)再延Z軸向上作一條直線，這條直線可能會與組成3D模型的三角形相交，三角形與直線的交點(diǎn)以及這個三角形的傾斜度會被記錄到網(wǎng)格里面。

現(xiàn)在每個網(wǎng)格里會記錄下一串被稱為支撐點(diǎn)的列表，每個支撐點(diǎn)包含一個高度和一個傾斜度信息。接下來會對每個網(wǎng)格的支撐點(diǎn)列表按照高度從低到高排序。根據(jù)這些信息就可以判斷模型上任意一個點(diǎn)是否需要支撐了，怎么判斷，我們看圖說話：

讓我們從底面開始延著一條網(wǎng)格中心往上走，起始我們是在模型外部，當(dāng)遇到第一個支撐點(diǎn)的時候，就從模型外部進(jìn)行了模型內(nèi)部，我們稱這個支撐點(diǎn)為進(jìn)點(diǎn)。

繼續(xù)向上，遇到了第二個支撐點(diǎn)，從模型內(nèi)部又退到了模型外部，我們稱這個支撐點(diǎn)為出點(diǎn)。

接著向上，我們可以發(fā)現(xiàn)，進(jìn)點(diǎn)與出點(diǎn)總是交替出現(xiàn)的。

利用這個規(guī)律，對于模型上任何一個點(diǎn)，我們只要找到這個點(diǎn)所對應(yīng)的網(wǎng)格，再找到這個網(wǎng)格里在這個點(diǎn)以上最近的一個支撐點(diǎn)，我們就可以得到兩個信息：這個點(diǎn)之上是否有模型懸空；這個點(diǎn)上面的懸空點(diǎn)的面的傾斜度是多少。

Cura界面的專家設(shè)置里面有支撐角度的設(shè)置，如果一個點(diǎn)處于模型懸空部分以下，并且懸空點(diǎn)傾斜度大于支撐角度，那這個點(diǎn)就是需要支撐的。所一個平臺上所有的需要支撐的點(diǎn)連接起來圍成的2D圖形就是支撐區(qū)域。

CuraEngine所使用的支撐算法比較粗糙，但勝在速度很快。先不說網(wǎng)格化后失去了精度，通過傾斜角度來判斷，模型下方一旦傾斜角發(fā)生了突變，像左圖這種從負(fù)45
度一下突變成正45度，傾斜角判斷無能為力，除非把它改大到60度，這樣的話，整個模型都會被過度支撐。這樣矯枉過正，既不科學(xué)，也浪費(fèi)材料和打印時間，還會對模型表面質(zhì)量帶來不好的影響。

科學(xué)的支撐算法應(yīng)該是找到模型局部最低點(diǎn)進(jìn)行支撐，最低點(diǎn)以上不一定需要支撐。因為FDM材料本身的粘性，使得材料的走線可以有一部分懸空而不坍塌，這個效果被稱為Overhang，只要上層材料的懸空距離小于一定的值，它就不需要支撐，這個距離以我的經(jīng)驗應(yīng)該在1/4到1/2線寬之間。我在Pango中就基于這個思路重新實現(xiàn)了支撐的算法，結(jié)果雖然

速度不如Cura的支撐算法那么快，但效果非常好，該撐的地方撐，不該撐的地方也不會多此一舉。

Pango的支撐算法我會在以后專文介紹。順帶一說，CuraEngine在下半年做了很大的改動，其中之一就是拋棄了之前的支撐算法，而新的算法也和我上面所講的思想異曲同工。我要聲明的是Pango的支撐算法和CuraEngine誰也沒有抄誰，我的算法是自己拍腦袋想出來的。算是英雄所見略同吧。

支撐范圍確定之后，也和組件一樣，可以有外墻、內(nèi)墻、填充、表面。依樣畫葫蘆即可。CuraEngine對于支撐，只會生成外墻和填充，Pango更多。

組件和支撐就是CuraEngine在這一步所生成的結(jié)果，這一步可以說是整個切片過程的核心，更詳細(xì)的過程會另文分析，敬請期待。

步驟四：路徑生成

地圈好了，就該在里面種菜了。這一步路徑生成就要開始規(guī)劃噴頭在不同的組件中怎么運(yùn)動。路徑按大類來分，有輪廓和填充兩種。

輪廓很簡單，沿著2D圖形的邊線走一圈即可。前一步所生成的外墻、內(nèi)墻都屬于輪廓，可以直接把它們的圖形以設(shè)置里的線寬轉(zhuǎn)換為輪廓路徑。

填充稍微要復(fù)雜一些，2D圖形指定的只是填充的邊界，而生成的路徑則是在邊界的范圍內(nèi)的條紋或網(wǎng)格結(jié)構(gòu)，就像窗簾或者漁網(wǎng)，如左圖。這兩種就最基本的結(jié)構(gòu)，當(dāng)然也許你還可以想出其它花式的填充，比如蜂窩狀或者S型，這些在新的Cura或者別的什么切片軟件里可能會實現(xiàn)，但我打印下來還是這兩種基本結(jié)構(gòu)更讓人放心。

CuraEngine在專家設(shè)置里可以對填充類型進(jìn)行選擇，里面除了條紋和網(wǎng)格外還有一個自動選項，默認(rèn)就是自動。自動模式會根據(jù)當(dāng)前的填充率進(jìn)行切換，當(dāng)填充率小于20%就用條紋填充，否則使用網(wǎng)格填充。因為網(wǎng)格結(jié)構(gòu)雖然更為合理，但它有一個問題，就是交點(diǎn)的地方會打兩次。填充率越高，交點(diǎn)越密，對打印質(zhì)量的影響會越大。我們知道，表面就是100%的填充，如果表面用網(wǎng)格打，不但無法打密實，表面還會坑坑洼洼，所以100%填充只能用條紋打，這就是CuraEngine推薦自動模式的原因。

至于填充率，就反映在線與線的間距上。100%填充率間距為0；0%填充率間距無限大，一根線條也不會有。

每個組件獨(dú)立的路徑生成好了，還要確定打印的先后順序。順序先好了可以少走彎路，打印速度和質(zhì)量都會有提升。路徑的順序以先近后遠(yuǎn)為基本原則：每打印完一條路徑，當(dāng)前位置是上一條路徑的終點(diǎn)；在當(dāng)前層里剩下還沒打印的路徑中挑選一條起點(diǎn)離當(dāng)前位置最近的一條路徑開打。路徑的起點(diǎn)可以是路徑中的任意一個點(diǎn)，程序會自行判斷。而路徑的終點(diǎn)有兩種可能：對于直線，圖形只有兩個點(diǎn)，終點(diǎn)就是除起點(diǎn)之外的那個點(diǎn)；對于輪廓，終點(diǎn)就是起點(diǎn)，因為輪廓是一個封閉圖形，從它的起點(diǎn)開始沿任意方向走一圈，最后還會回到起點(diǎn)。CuraEngine對路徑選擇做了一個估值，除了考慮到先近后遠(yuǎn)外，還順便參考了下一個點(diǎn)相對于當(dāng)前點(diǎn)的方向，它的物理意義就是減少噴頭轉(zhuǎn)彎。賽車在直道上開得一定比彎道快，不是么。

路徑的順序也確定了，還有一個問題需要考慮：如果前后兩條路徑首尾相連，那直接走就是了，但大多數(shù)情況并非如此，前一條路徑的終點(diǎn)往往和后一條路徑起點(diǎn)之間有一段距離。這時候去往下一點(diǎn)的路上要小心了，肯定不能繼續(xù)擠出材料，否則輕則拉絲，重則模型面目全非。這段路噴頭就需要空走，即噴頭只移動，不吐絲，那只要把擠出機(jī)停下來不轉(zhuǎn)就行了嗎？也不行，因為前面分析過，擠出機(jī)的速度要傳導(dǎo)到噴嘴，有一個延遲，不是你說停它就立即停下來的。這是FDM打印的通病，解決辦法就是回抽。所謂回抽，就是在空走之前先讓擠出機(jī)高速反轉(zhuǎn)一段材料，這樣就可以瞬間把加熱腔里的材料抽光，再移動過去，中間就不會擠出材料，到了下一個點(diǎn)，在打印之前，先把剛才抽回去的絲再按一樣的長度放回來，繼續(xù)打印?；爻榭梢院芎玫亟鉀Q空走拉絲的問題，但是它很慢，以抽一次0.5秒來算的話，如果打印一個表面，0.4線寬，10厘米的距離至少回抽25下，10幾秒鐘的時間一層，幾百上千層打下來，光回抽所用的時間就是幾個小時，是可忍孰不可忍！

CuraEngine給我們提供了解決方案就是Comb，也就是繞路。我們先來看，是不是所有的回抽都是必需的呢？不回抽會拉絲是肯定的，但如果需要空走的路徑本來就要打印的，那拉絲又有何妨。按這個思路，就可以給每個組件設(shè)定一個邊界，只要路徑的起點(diǎn)和終點(diǎn)都在這個邊界之內(nèi)的，空走都不回抽。這樣可以解決80%的問題，但如果是左圖這樣的情況就行不通。

紅色是起點(diǎn)，綠色是終點(diǎn)，直接走過去會走出邊界的范圍。這時我們就要繞一點(diǎn)路，走一條曲線到達(dá)我們的目的地。這就是Comb所做的事情，在Cura專家設(shè)置里面可以對Comb進(jìn)行設(shè)置，選擇開啟、關(guān)閉還有表面不Comb。Comb可以大幅節(jié)省打印時間，但是同一個地方打印多次對模型質(zhì)量還是會有細(xì)微的影響，個中利弊，交給用戶自己判斷。

Comb的調(diào)整是個細(xì)致活，Pango花了相當(dāng)多的時間來微調(diào)Comb功能以求達(dá)到更好的效果，過程繁瑣，不再贅述。

至此路徑生成完成，更詳細(xì)的過程另文分析，敬請期待。

步驟五：gcode生成
路徑都生成好了，還需要翻譯對打印機(jī)可以實別的gcode代碼才行。這一步花樣不多，按部就班即可。

先讓打印機(jī)做一些準(zhǔn)備工作：歸零、加熱噴頭和平臺、抬高噴頭、擠一小段絲、風(fēng)扇設(shè)置。

從下到上一層一層打印，每層打印之前先用G0抬高Z坐標(biāo)到相應(yīng)位置。

按照路徑，每個點(diǎn)生成一條gcode。其中空走G0；邊擠邊走用G1，Cura的設(shè)置里有絲材的直徑、線寬，可以算出走這些距離需要擠出多少材料；G0和G1的速度也都在設(shè)置里可以調(diào)整。

若需回抽，用G1生成一條E軸倒退的代碼。在下一條G1執(zhí)行之前，再用G1生成一條相應(yīng)的E軸前進(jìn)的代碼。

所有層都打完后讓打印機(jī)做一些收尾工作：關(guān)閉加熱、XY歸零、電機(jī)釋放。

生成gcode的過程中，CuraEngine也會模擬一遍打印過程，用來計算出打印所需要的時間和材料長度，這些也會寫在gcode的注釋里供用戶參考。

gcode生成不用另文詳細(xì)分析，但是gcode的說明還是可以專文分析一下，敬請期待。

待續(xù)

寫了這么多，Cura的切片流程也只能講個大概，也算是個提綱，希望對大家有所幫助。我計劃對于上面的第一個步驟再專文分析。除此之外，還有Cura界面部分以及Cura與CuraEngine的通訊也可以講講。之后就是我半年創(chuàng)作，自我感覺良好到覺得可以超越Cura的Pango，也是不說不快的。

酷?。。?！