專利名稱：曲面網(wǎng)格化方法
技術領域：
本發(fā)明涉及一種曲面重建方法，尤其是關于一種曲面網(wǎng)格化方法。
背景技術：
在產(chǎn)品的開發(fā)及制造過程中，幾何造型技術已使用得相當廣泛，但是，由于種種原因， 仍有許多產(chǎn)品并非由CAD模型描述，設計和制造者面對的是實物樣件。為了適應先進制造技 術的發(fā)展，需要通過一定途徑，將這些實物轉化為CAD模型，使之能利用CAD、 CAM、 RPM、 PDM及CMIS先進技術進行處理或管理。目前，與這種從實物樣件獲取產(chǎn)品數(shù)學模型技術相關 的技術，已發(fā)展成為CAD/CAM中的一個相對獨立的范疇，稱為"反求工程"(Reverse Engineering)。通過反求工程復現(xiàn)實物的CAD模型，簡稱模型重建。
目前較為成熟的模型重建技術是通過構建曲面來實現(xiàn)的，在構建曲面的過程中，根據(jù)實 物模型的坐標測量數(shù)據(jù)構建曲面、對曲面進行三角網(wǎng)格化以逼近實物面在實踐中被廣泛應用 。常用的曲面三角網(wǎng)格化方法有迭代法和迪式(Delaunay)算法，然而這兩種方法都存在不 足之處用迭代法直接算出曲面的控制點，再根據(jù)控制點坐標擬合三角形的速度非常慢，且 擬合出來的三角形大小不一致，速度和精度都不盡人意；用Dela皿ay算法所獲得的三角網(wǎng)格 比較規(guī)則，但Dela皿ay算法致力于提高剖分質量，而導致剖分生成的三角面片的數(shù)量巨大， 所以速度緩慢。

發(fā)明內容
鑒于以上內容，有必要提供一種曲面網(wǎng)格化方法，能夠根據(jù)用戶實際需求快速生成均勻 的三角形網(wǎng)格曲面，提高曲面重建的速度和精度。
一種曲面網(wǎng)格化方法，該方法包括以下步驟(A)接收來自數(shù)據(jù)庫的曲面對象及用戶 設置的曲面網(wǎng)格化后曲面上的三角形網(wǎng)點數(shù)目；(B)根據(jù)曲面的面積及曲面網(wǎng)格化所需的 三角形網(wǎng)點數(shù)目計算劃分參數(shù)平面時水平方向及垂直方向的步長；(C)獲取曲面上的封閉 邊界曲線及各封閉邊界曲線的長度；(D)根據(jù)各封閉邊界曲線的長度及劃分參數(shù)平面時水 平方向及垂直方向的步長計算出曲面網(wǎng)格化時曲面上各封閉邊界曲線被分得的網(wǎng)點數(shù)；(E )由各封閉邊界曲線的參數(shù)方程得到各被分得的網(wǎng)點的二維坐標值；(F)將上述二維坐標 值所對應的參數(shù)平面內的點用線段依次連接，得到與曲面上各封閉邊界曲線相應的參數(shù)平面 內的封閉邊界曲線，并根據(jù)劃分參數(shù)平面時水平方向及垂直方向的步長將參數(shù)平面內最大封
閉邊界曲線的包圍盒用水平直線及垂直直線劃分成若干個小包圍盒；(G)獲取水平、垂直 直線與參數(shù)平面內各封閉邊界曲線的交點、各封閉邊界曲線上的頂點，以及各小包圍盒頂點 中落入相鄰封閉邊界曲線之間區(qū)域的點，將所述交點及頂點加入相應的小包圍盒內；(H) 將參數(shù)平面內各小包圍盒內的上述交點、頂點及各小包圍盒頂點中落入相鄰封閉邊界曲線之 間區(qū)域的點連接成三角形，并調整所有三角形的頂點順序為逆時針；(I)根據(jù)曲面的參數(shù) 方程將所有三角形頂點的二維坐標映射回三維空間曲面中的三維點坐標，并根據(jù)三角形的三 維點坐標計算每個三角形所在平面的法向量；及(J)根據(jù)三角形信息，包括三角形頂點三 維坐標、三角形頂點的連接方式、及各三角形所在平面的法向量輸出三角網(wǎng)格化后的曲面。
相較于現(xiàn)有技術，本發(fā)明所提供的曲面網(wǎng)格化方法將曲面網(wǎng)格劃分問題轉化為對參數(shù)平 面的平面網(wǎng)格劃分問題，再通過映射將參數(shù)平面網(wǎng)格化得到的三角形轉化為三維空間的三角 形，曲面網(wǎng)格化得到的三角形分布均勻、數(shù)量可根據(jù)用戶實際需求靈活設定。


圖l是本發(fā)明曲面網(wǎng)格化方法較佳實施例的應用環(huán)境圖。
圖2是本發(fā)明曲面網(wǎng)格化方法較佳實施例的主流程圖。
圖3是圖2中步驟S102的細化流程圖。
圖4是圖3中步驟S203的細化流程圖。
圖5是圖2中步驟S104的細化流程圖。
圖6是圖2中步驟S 106中獲取內點的細化流程圖。
圖7是圖2中步驟S107的細化流程圖。
圖8是一個曲面對象的三維圖示。
圖9是圖8中曲面對象在參數(shù)平面中的相應圖示，及執(zhí)行圖2中步驟S105劃分參數(shù)平面的 示意圖。
圖1 O是執(zhí)行圖2中步驟S 107得到的點的示意圖。
圖11是執(zhí)行圖5中步驟S303分割參數(shù)平面時各封閉邊界曲線分得的網(wǎng)點數(shù)的示意圖。 圖12和圖13舉例說明圖7中連接中間三角形和邊界三角形的步驟。 圖14是在參數(shù)平面內生成邊界三角形的圖示。
圖15是在參數(shù)平面內生成所有三角形包括邊界三角形和中間三角形的圖示，其中相較于 圖14多出的三角形為中間三角形。
圖16是圖15的參數(shù)平面內生成所有三角形映射回三維空間后得到的網(wǎng)格化后均勻分布的 三角形的示意圖。
具體實施例方式
如圖1所示，是本發(fā)明曲面網(wǎng)格化方法較佳實施例的應用環(huán)境圖。該應用環(huán)境圖主要包 括數(shù)據(jù)庫IO、與數(shù)據(jù)庫10連接的計算機20以及與計算機20連接的顯示設備30。
數(shù)據(jù)庫10用于儲存曲面對象。每一個曲面對象包括該曲面的參數(shù)方程，從曲面參數(shù)方程 可得到曲面上三維點與參數(shù)平面中二維點的對應關系，曲面的面積，曲面上封閉邊界曲線的 參數(shù)方程；由封閉邊界曲線的參數(shù)方程可以得到封閉邊界曲線的長度。于本實施例中，參數(shù) 平面的水平方向用U表示，垂直方向用V表示，U、 V的取值范圍為0-1。
計算機20用于接收來自數(shù)據(jù)庫10的曲面對象及用戶設置的曲面網(wǎng)格化后的三角形網(wǎng)點數(shù) 目并執(zhí)行曲面三角網(wǎng)格化程序200，將曲面網(wǎng)格劃分問題轉化為對參數(shù)平面的平面網(wǎng)格劃分 問題，再通過映射將參數(shù)平面網(wǎng)格化得到的三角形轉化為三維空間的三角形，得到均勻分布 的三角形網(wǎng)格曲面。
顯示設備30提供一顯示界面用于顯示三角網(wǎng)格化前的曲面及三角網(wǎng)格化后的曲面。 曲面三角網(wǎng)格化程序200包括接收模塊210、計算模塊220、獲取模塊230、網(wǎng)格化模塊 240及輸出模塊250。
接收模塊210用于接收來自數(shù)據(jù)庫10的曲面對象及用戶設置的曲面網(wǎng)格化后曲面上的三 角形網(wǎng)點數(shù)目。
計算模塊220用于根據(jù)曲面的參數(shù)方程計算曲面面積，并根據(jù)曲面的面積及曲面網(wǎng)格化 所需的三角形網(wǎng)點數(shù)目計算劃分參數(shù)平面時U方向及V方向的步長。
獲取模塊230用于獲取曲面上的封閉邊界曲線及各封閉邊界曲線的長度。
計算模塊220還用于根據(jù)封閉邊界曲線的參數(shù)方程計算封閉邊界曲線的長度，并根據(jù)各 封閉邊界曲線的長度及劃分參數(shù)平面時U方向及V方向的步長計算出曲面網(wǎng)格化時曲面上各封 閉邊界曲線被分得的網(wǎng)點數(shù)。
接下來，獲取模塊230由各封閉邊界曲線的參數(shù)方程得到各被分得的網(wǎng)點的二維UV坐標值。
網(wǎng)格化模塊240用于將上述二維UV坐標值所對應的參數(shù)平面內的點用線段依次連接，得 到與曲面上各封閉邊界曲線相應的參數(shù)平面內的封閉邊界曲線，并根據(jù)劃分參數(shù)平面時U方 向及V方向的步長將參數(shù)平面內最大封閉邊界曲線的包圍盒用若干平行于U方向的直線(以下 簡稱U直線)及若干平行于V方向的直線(以下簡稱V直線)劃分成若干個小包圍盒。
接下來，獲取模塊230獲取U、 V直線與參數(shù)平面內各封閉邊界曲線的交點、各邊界曲線 上的頂點，以及各小包圍盒頂點中落入相鄰封閉邊界曲線間區(qū)域內的點(以下簡稱內點)。
網(wǎng)格化模塊240將參數(shù)平面內的上述交點、頂點及內點連接成三角形，并調整所有三角 形的頂點連接順序為逆時針。
計算模塊220根據(jù)曲面的參數(shù)方程將所有三角形頂點的二維UV坐標映射回三維空間曲面 中的三維點坐標，并根據(jù)三角形的三維點坐標計算每個三角形所在平面的法向量。
輸出模塊250根據(jù)三角形信息，包括三角形頂點三維坐標、三角形頂點的連接順序、及 各三角形所在平面的法向量輸出三角網(wǎng)格化后的曲面。
如圖2所示，是本發(fā)明曲面網(wǎng)格化方法較佳實施例的主流程圖。首先，接收模塊210接收 用戶自數(shù)據(jù)庫10導入的一個曲面對象，該曲面對象包括該曲面的參數(shù)方程(步驟SIOO)。接 下來，接收模塊210接收用戶設置的曲面網(wǎng)格化所需的三角形網(wǎng)點數(shù)目(步驟SIOI)。計算 模塊220根據(jù)該曲面的面積和所設置的三角形網(wǎng)點數(shù)目計算分割參數(shù)平面時U方向及V方向的 步長(具體說明請參見圖3)(步驟S102)。獲取模塊230從曲面對象中獲取封閉邊界曲線( loops)的參數(shù)方程，計算模塊220根據(jù)該封閉邊界曲線的參數(shù)方程得到該封閉邊界曲線的長 度(步驟S103)。計算模塊220根據(jù)封閉邊界曲線的長度及上述U方向及V方向的步長中的較 小值計算出在曲面三角網(wǎng)格化過程中各封閉邊界曲線被分得的網(wǎng)點數(shù)，并根據(jù)該封閉邊界曲 線的參數(shù)方程得到這些網(wǎng)點的二維UV坐標(步驟S104)(具體介紹請參見圖5)。網(wǎng)格化模 塊240將在這些二維點用線段依次連接，得到參數(shù)平面中相應的封閉邊界曲線(如圖9中的四 條封閉邊界曲線分別與圖8中三維空間中曲面的四條封閉邊界曲線對應)，并根據(jù)計算所得 步長，用U、 V直線將參數(shù)平面中最大封閉邊界曲線的包圍盒等分成若干個小包圍盒(步驟 S105)(如圖9所示)。獲取模塊230獲取U、 V直線與各封閉邊界曲線的交點(以下簡稱交點 )、小包圍盒頂點中落入相鄰封閉邊界曲線間區(qū)域(如圖ll中的loopO， loopl， loop2， 1o叩3之間的區(qū)域)內的點(以下簡稱內點)(具體介紹請參見圖6)，以及封閉邊界曲線的 頂點(以下簡稱頂點)，并將所述交點和頂點加入對應的小包圍盒內，上述交點、內點和頂 點稱作有效點(步驟S106)(如圖10所示)。網(wǎng)格化模塊240將參數(shù)平面中各小包圍盒內的 上述有效點連接成三角形，并調整所有三角形頂點連接順序為順時針(具體介紹請參見圖7 )(步驟S107)。計算模塊220根據(jù)曲面的參數(shù)方程將所有三角形頂點的二維UV坐標映射回 三維空間曲面中的三維點坐標，并根據(jù)三角形的三維點坐標計算每個三角形所在平面的法向 量(步驟S108)。輸出模塊250根據(jù)三角形信息，包括三角形頂點三維坐標、三角形頂點的 連接順序、及各三角形所在平面的法向量輸出三角網(wǎng)格化后的曲面(步驟S109)。
如圖3所示，是圖2中步驟S102的細化流程圖。計算模塊220根據(jù)曲面的參數(shù)方程計算得 到曲面的面積Area (步驟S200)。計算模塊220根據(jù)曲面的面積及三角形網(wǎng)點數(shù)目N計算劃分
參數(shù)平面時U方向的理論步長Len. UO及V方向的理論步長Len. V0，
Len. U0=Len. V0=sqrt (Area/N)步驟S201 )。接收模塊210接收用戶所設置的劃分參數(shù)平面時 U、 V方向的測試步長Len. test (步驟S202)。計算模塊220以所述測試步長Len. test為間 距分別沿U、 V方向取點，得到U方向的點與V方向點各自所對應的三維點，并計算U方向兩點 間的距離與相應三維點兩點間的距離的比例關系Ru及V方向兩點間的距離與相應三維兩點間 點距離的比例關系Rv (步驟S203)(具體介紹請參見圖4)。計算模塊220根據(jù)所述U方向及V 方向的理論步長及相應方向的比例關系計算劃分參數(shù)平面時U方向的實際步長Len. U及V方向 的實際步長Len. V， Len. U= Len.UO*Ru， Len. V=Len. VO*Rv (步驟S204)。
圖4以具體例子說明圖3中的步驟S203。首先，計算模塊220在參數(shù)平面中任意取一點A ( 步驟S2031)。接下來，計算模塊220以點A為基點沿U方向、V方向分別取測試步長Len. test的距離，得到第二點B及第三點C (步驟S2032)。計算模塊220根據(jù)點A、 B、 C在UV參數(shù) 平面的UV坐標及曲面的參數(shù)方程計算得到點A、 B、 C所對應的三維空間點A1、 Bl、 Cl的三維 坐標(步驟S2033)。計算模塊220計算點A1、 B1的距離D1，點A1、 C1的距離D2 (步驟S2034 )。計算模塊220計算得到UV參數(shù)平面中U方向點A、 B之間的距離Len. test與A、 B點相應的 三維空間點A1、 Bl之間的距離Dl的比例關系R『Len. test/Dl，同理計算得到Rv二 Len. test/D2 (步驟S2035)。本實施例以沿U方向、V方向各取一個點舉例說明UV參數(shù)平面中U方 向、V方向點與點間的距離與相應三維空間點與三維空間點的距離的比例關系，沿U方向、V 方向取點越多，通過平均值法計算得到的比例關系更準確，在此不再贅述。
如圖5所示，是圖2中步驟S104的細化流程圖。獲取模塊230根據(jù)曲面的參數(shù)方程得到曲 面上的封閉邊界曲線(loops)(步驟S301)。獲取模塊230獲取分割參數(shù)平面時U方向及V方 向的實際步長Len. U， Len. V中的較小值(如若Len. U〈Len. V則取Len. U)，并根據(jù)該封 閉邊界曲線的參數(shù)方程得到該封閉邊界曲線的長度L (步驟S302)。計算模塊220根據(jù)各封閉 邊界曲線的長度L及所述較小值(如Len. U)計算得到曲面網(wǎng)格化時各封閉邊界曲線分得的 網(wǎng)點數(shù)n， n=L/Len. U (步驟S303)(如圖ll所示)。獲取模塊230根據(jù)各封閉邊界曲線的參 數(shù)方程得到該封閉邊界曲線分得的各網(wǎng)點的二維UV坐標(步驟S304)。
如圖6所示，是圖2中步驟S106中獲取內點的細化流程圖。獲取模塊230獲取參數(shù)平面中 的封閉邊界曲線(步驟S401)。獲取模塊230獲取一條U直線或V直線(步驟S402)。獲取模 塊230獲取該U直線或V直線與參數(shù)平面中封閉邊界曲線的第一個交點(步驟S403)。獲取模 塊230繼續(xù)沿U方向或V方向査找，得到該U直線或V直線與參數(shù)平面中封閉邊界曲線的第二個 交點(步驟S404)。獲取模塊230獲取第一個交點與第二個交點之間的U方向線段或V方向線
段上的小包圍盒的頂點，該頂點即為內點(步驟S405)。獲取模塊230繼續(xù)沿U方向或V方向 査找，判斷該U直線或V直線與參數(shù)平面中封閉邊界曲線是否還有新的交點(步驟S406)。若 有新的交點，則獲取模塊230以該新的交點為第一個交點(步驟S407)，重復步驟S404至步 驟S406。若無新的交點，則獲取模塊230判斷是否還有U、 V直線未査找(步驟S408)。若還 有U、 V直線未査找，則返回步驟S402;若所有U、 V直線都査找完畢，則流程結束。
如圖7所示，是圖2中步驟S107的細化流程圖。獲取模塊230獲取參數(shù)平面內的一個小包 圍盒(步驟S501)。獲取模塊230判斷該小包圍盒內有效點中是否有交點和頂點(步驟S502 )。若該小包圍盒的有效點中無交點和頂點，則該小包圍盒內的有效點均為內點，即該小包 圍盒的四個頂點均落入相鄰封閉邊界曲線之間的區(qū)域，則網(wǎng)格化模塊240將該小包圍盒的四 個頂點按照逆時針順序連接成兩個三角形，如圖12中所示的小包圍盒box4的四個頂點P、 Q、 I、 O均落入曲面邊界內，則網(wǎng)格化模塊240將頂點P、 Q、 I、 O按照逆時針順序連接組成兩個 三角形OQP、 OIQ，如圖13所示(步驟S503)，之后流程進入步驟S512。
在步驟S502，若獲取模塊230判斷該小包圍盒內有效點中有交點或/和頂點，則獲取模塊 230依次査找該小包圍盒四條邊得到各邊上的有效點，如圖12所示boxl的邊IJ、邊IB上無 交點，邊BA上有交點E、 F，邊AJ上無交點，boxl有邊界曲線落入該boxl內的頂點M，則獲取 模塊230獲取點M、 E、 F; box2的邊AB上有交點F、 E，邊BC上無交點，邊CD上有交點H、內點 D，邊DA上有內點D、交點G，則獲取模塊230獲取點F、 E、 H、 D、 G; box3的邊DC上有內點D、 交點H，邊CP上有交點N、內點P、邊PO上有內點P、 0、邊OD上有內點D、 0，則獲取模塊230獲 取點D、 H、 N、 P、 0 (步驟S504)。
計算模塊220將各邊上的有效點按照査找時各有效點到該邊上出現(xiàn)的第一個小包圍盒頂 點的距離從小到大進行排序，如査找boxl各邊時按照JI、 IB、 BA、 AJ依次査找，有效點E、 F到邊BA出現(xiàn)的第一個頂點B的距離lBEl〈lBFl，則boxl內有效點M、 E、 F的順序為M、 E、 F; 按照與査找boxl各邊時同樣的査找順序，box2的邊AB上有效點F、 E到邊AB出現(xiàn)的的第一個頂 點A的距離IAFI 〈 IAE I ，邊CD上有效點H、 D到邊CD上出現(xiàn)的第一個頂點C的距離| CH | 〈 | CD | ，邊 DG上有效點D、 G到邊DG出現(xiàn)的第一個頂點D的距離0〈lDGl ，故box2內有效點的排序順序為F、 E、 H、 D、 G;同理得到box3內有效點的排序順序為D、 H、 N、 P、 0 (步驟S505)。
網(wǎng)格化模塊240將上述排序后的有效點按照排序順序連接成一個多邊形，如圖13所示， 連接boxl內的有效點得到多邊形MEF，連接box2內的有效點得到多邊形FEHDG，連接box3內的 有效點得到多邊形DHNPO (步驟S506)。
獲取模塊230按照有效點的排序順序査找該多邊形的邊得到該多邊形的第一個最小內角
，如依有效點的排序順序査找多邊形MEF各邊得到第一個最小內角ZFME，循環(huán)多邊形FEHDG 各邊得到第一個最小內角ZEHD，依有效點的排序順序査找多邊形D麗PO得到第一個最小內角 ZNP0 (步驟S507)。網(wǎng)格化模塊240連接該多邊形相應兩個頂點得到以該最小內角為內角的 一個邊界三角形，如圖13所示，連接多邊形MEF的頂點E、 F得到邊界三角形AMEF;連接多邊 形FEHDG的頂點E、 D得到三角形AEHD;連接多邊形DHNPO的頂點N、 O得到三角形ANPO (步驟 S508)。
網(wǎng)格化模塊240判斷除去該第一最小內角的頂點后該多邊形剩余有效點是否組成多邊形 (步驟S509)。若剩余有效點不足以組成多邊形，如圖13所示，多邊形MEF除去頂點M僅剩下 有效點E、 F，兩個點E、 F不足以組成多邊形，則流程進入步驟S512。
若剩余有效點可以組成多邊形，如圖13所示，多邊形FEHDG除去ZEHD的頂點H剩下有效 點E、 D、 G、 F，多邊形DHNPO除去ANPO的頂點P剩下有效點H、 N、 0、 D，則流程進入步驟 S510，網(wǎng)格化模塊240判斷剩余有效點組成的多邊形是否為三角形，若不為三角形，如有效 點E、 D、 G、 F組成多邊形EDGF，有效點H、 N、 0、 D組成多邊形HNOD，則對該由剩余有效點組 成的多邊形自步驟S507開始重復直到剩下有效點組成的多邊形為三角形。
在步驟S510，若網(wǎng)格化模塊240判斷剩余有效點組成的多邊形為三角形，則網(wǎng)格化模塊 240調整該小包圍盒內有效點組成的所有邊界三角形的頂點連接順序為逆時針(步驟S511) 。接下來，網(wǎng)格化模塊240判斷是否參數(shù)平面內的所有小包圍盒都査找完畢(步驟S512)。 若還有小包圍盒未査找完畢則對待査找的小包圍盒自步驟S501開始重復流程直到參數(shù)平面內 的所有有效點都被連接成邊界三角形和中間三角形；若所有小包圍盒都査找完畢，則流程結 束。
如圖14是在參數(shù)平面內生成邊界三角形的圖示，圖15是在參數(shù)平面內生成所有三角形的 圖示，其中相較于圖14多出的三角形為依據(jù)圖7中步驟S503直接連接內點組成的三角形；圖 16是圖15的參數(shù)平面內生成所有三角形映射回三維空間后得到的網(wǎng)格化后均勻分布的三角形 的示意圖。
權利要求
權利要求1一種曲面網(wǎng)格化方法，其特征在于，該方法包括以下步驟(a)接收來自數(shù)據(jù)庫的曲面對象及用戶設置的曲面網(wǎng)格化所需的三角形網(wǎng)點數(shù)目，所述曲面對象包括曲面的參數(shù)方程及曲面上各封閉邊界曲線的參數(shù)方程；(b)根據(jù)曲面的參數(shù)方程計算曲面面積，并根據(jù)曲面面積及曲面網(wǎng)格化所需的三角形網(wǎng)點數(shù)目計算劃分參數(shù)平面時水平方向及垂直方向的實際步長；(c)根據(jù)各封閉邊界曲線的參數(shù)方程計算得到曲面上各封閉邊界曲線的長度；(d)根據(jù)曲面上各封閉邊界曲線的長度及劃分參數(shù)平面時水平方向及垂直方向的步長計算出曲面網(wǎng)格化時曲面上各封閉邊界曲線被分得的網(wǎng)點數(shù)；(e)由各封閉邊界曲線的參數(shù)方程得到曲面上各封閉邊界曲線被分得的網(wǎng)點的二維坐標值；(f)將上述二維坐標值所對應的參數(shù)平面內的點用線段依次連接，得到與曲面上各封閉邊界曲線相應的參數(shù)平面內的封閉邊界曲線，并根據(jù)劃分參數(shù)平面時水平方向及垂直方向的步長將參數(shù)平面內最大封閉邊界曲線的包圍盒用若干水平直線及若干垂直直線劃分成若干個小包圍盒；(g)獲取所述若干水平、垂直直線與參數(shù)平面內各封閉邊界曲線的交點、各封閉邊界曲線上的頂點，以及各小包圍盒頂點中落入?yún)?shù)平面內相鄰封閉邊界曲線之間區(qū)域的點，并將所述交點及頂點加入相應的小包圍盒內；(h)將參數(shù)平面內各小包圍盒內的上述交點、頂點及各小包圍盒頂點中落入?yún)?shù)平面內相鄰封閉邊界曲線之間區(qū)域的點連接成三角形，并調整所有三角形的頂點連接順序為逆時針；(i)根據(jù)曲面的參數(shù)方程將所有三角形頂點的二維坐標映射回曲面中的三維點坐標，并根據(jù)三角形的三維點坐標計算每個三角形所在平面的法向量；及(j)根據(jù)各三角形頂點的三維坐標、各三角形頂點的連接順序及各三角形所在平面的法向量，輸出三角網(wǎng)格化后的曲面。
2.如權利要求l所述的曲面網(wǎng)格化方法，其特征在于，其中步驟(b )包括以下步驟根據(jù)曲面面積及三角形網(wǎng)點數(shù)目計算劃分參數(shù)平面時水平方向及垂直方向的理論步長 接收用戶設置的劃分參數(shù)平面的測試步長；根據(jù)所述理論步長及測試步長計算參數(shù)平面內水平方向兩點之間距離與相應兩三維點 之間距離的比例關系、垂直方向兩點之間距離與相應兩三維點之間距離的比例關系；及根據(jù)所述理論步長及所述比例關系分別計算劃分參數(shù)平面時水平方向及垂直方向的實 際步長。
3.如權利要求2所述的曲面網(wǎng)格化方法，其特征在于，其中步驟(d )包括以下步驟獲取劃分參數(shù)平面時水平方向及垂直方向的實際步長中的較小值；及 根據(jù)曲面上各封閉邊界曲線的長度及所述較小值計算曲面網(wǎng)格化時曲面上各封閉邊界曲線被分得的網(wǎng)點數(shù)。
4.如權利要求l所述的曲面網(wǎng)格化方法，其特征在于，其中步驟(g)中獲取各小包圍盒頂點中落入?yún)?shù)平面內相鄰封閉邊界曲線之間區(qū)域的點包括以下步驟 獲取一條水平直線或垂直直線；獲取該水平直線或垂直直線與參數(shù)平面內封閉邊界曲線的第一個交點； 繼續(xù)沿水平方向或垂直方向査找得到該水平直線或垂直直線與所述參數(shù)平面內封閉邊界曲線的第二個交點；獲取第一個交點與第二個交點之間的線段上的小包圍盒的頂點； 判斷該水平直線或垂直直線與所述參數(shù)平面內封閉邊界曲線是否還有新的交點； 若有新的交點，則以該新的交點為第一個交點，返回査找該水平直線或垂直直線與所述參數(shù)平面內封閉邊界曲線的第二個交點的步驟；若無新的交點，則判斷參數(shù)平面內是否還有其它水平直線或垂直直線未査找；及 若還有未査找的水平直線或垂直直線，則返回獲取水平直線或垂直直線的步驟。
5.如權利要求l所述的曲面網(wǎng)格化方法，其特征在于，其中步驟(h)包括以下步驟(hl)獲取參數(shù)平面內的一個小包圍盒；(h2)判斷該小包圍盒中是否有水平、垂直直線與參數(shù)平面內封閉邊界曲線的交點或 參數(shù)平面內封閉邊界曲線上的頂點；(h3)若該小包圍盒內無水平、垂直直線與參數(shù)平面內封閉邊界曲線的交點及參數(shù)平 面內封閉邊界曲線上的頂點，則將該小包圍盒落入?yún)?shù)平面內相鄰封閉邊界曲線之間區(qū)域的 四個頂點按照逆時針順序連接成兩個三角形；(h4)判斷參數(shù)平面內的小包圍盒是否都査找完畢；(h5)若還有小包圍盒未査找則返回獲取小包圍盒的步驟。
6.如權利要求5所述的曲面網(wǎng)格化方法，其特征在于，在步驟(h2 )若判斷該小包圍盒內有效點中有水平、垂直直線與封閉邊界曲線的交點或封閉邊界曲線上 的頂點，則執(zhí)行以下步驟(kl)査找該小包圍盒的四條邊得到各邊上的水平、垂直直線與封閉邊界曲線的交點 、該小包圍盒落入封閉邊界曲線邊界內區(qū)域的頂點及封閉邊界曲線上的頂點；(k2)將各邊上的上述水平、垂直直線與封閉邊界曲線的交點、該小包圍盒落入封閉 邊界曲線邊界內區(qū)域的頂點及封閉邊界曲線上的頂點按照査找時各點到該邊上出現(xiàn)的第一個 小包圍盒頂點的距離從小到大的順序進行排序；(k3)將該小包圍盒內的上述水平、垂直直線與封閉邊界曲線的交點、該小包圍盒落 入封閉邊界曲線邊界內區(qū)域的頂點及封閉邊界曲線上的頂點按照排序順序連接成一個多邊形(k4)按照排序順序査找該多邊形的邊，得到該多邊形的第一個最小內角； (k5 )連接多邊形相應兩個頂點得到以該第一個最小內角為內角的一個三角形； (k6)判斷該多邊形除去該第一個最小內角的頂點后，剩余的頂點是否組成一個多邊形；(k7)若該多邊形剩余的頂點不足以組成一個多邊形，則判斷是否所有小包圍盒都査 找完畢；(k8)若還有小包圍盒未査找完畢，則返回獲取小包盒的步驟。
7.如權利要求6所述的曲面網(wǎng)格化方法，其特征在于，在步驟(k6 )若判斷得到該多邊形除去該第一個最小內角的頂點后，剩余的頂點組成一個多邊形，則執(zhí) 行以下步驟(ml )判斷該多邊形剩余頂點組成的多邊形是否為一個三角形；(m2)若該多邊形剩余頂點組成的多邊形為一個三角形，則調整該包圍盒內水平、垂 直直線與封閉邊界曲線的交點、該小包圍盒落入封閉邊界曲線邊界內區(qū)域的頂點及封閉邊界 曲線上的頂點組成的所有三角形的頂點連接順序為逆時針；及 (m3)進入判斷是否所有小包圍盒都査找完畢的步驟。
8.如權利要求7所述的曲面網(wǎng)格化方法，其特征在于，在步驟(ml )若判斷得到該多邊形剩余頂點組成的多邊形不為一個三角形，則返回步驟(k4)。
全文摘要
本發(fā)明提供一種曲面網(wǎng)格化方法。該方法將曲面網(wǎng)格劃分問題轉化為對參數(shù)平面的平面網(wǎng)格劃分問題，再通過映射將參數(shù)平面網(wǎng)格化得到的三角形轉化為三維空間的三角形，曲面網(wǎng)格化得到的三角形分布均勻、數(shù)量可根據(jù)用戶實際需求靈活設定。
文檔編號G06T11/00GK101383047SQ200710201590
公開日2009年3月11日 申請日期2007年9月3日 優(yōu)先權日2007年9月3日
發(fā)明者吳新元, 張旨光, 敏 王 申請人:鴻富錦精密工業(yè)(深圳)有限公司;鴻海精密工業(yè)股份有限公司