尹小奎， 李奇敏， 葉仲泉， 蔣恒恒

（1.重慶大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，重慶 400044；2. 重慶大學(xué)機(jī)械工程學(xué)院，重慶 400044）

在飛機(jī)、汽車(chē)、船舶以及家用電器等的計(jì)算機(jī)輔助設(shè)計(jì)中經(jīng)常遇到許多由二次曲線(xiàn)弧與二次曲面所表示的形狀，為描述這些形狀，Gordon和Riesenfeld于1974年提出了B樣條曲線(xiàn)曲面，較成功地解決了曲線(xiàn)曲面局部控制問(wèn)題，并在參數(shù)連續(xù)性基礎(chǔ)上解決了連接問(wèn)題。由于B樣條方法不能精確表示圓錐及初等解析曲面，Verspril于1975年首次提出有理B樣條方法。后來(lái)由于Pieg和Tille等的工作，終于使非均勻有理B樣條(NURBS)方法成為現(xiàn)代曲面造型中最為廣泛流行的技術(shù)。隨著生產(chǎn)的發(fā)展，人們對(duì)于產(chǎn)品的外形也有了要求，其中之一就是外形的光順性。如果曲線(xiàn)、曲面不光順，就不能滿(mǎn)足產(chǎn)品的設(shè)計(jì)要求，也不便于加工。因此，NURBS曲線(xiàn)、曲面的光順處理就成為 CAD/CAM 中非常重要的一個(gè)問(wèn)題。

由于曲線(xiàn)、曲面的光順性涉及幾何外形的美觀性，難免受到主觀因素的影響，而且在不同的實(shí)際問(wèn)題中，對(duì)其的要求也不同。對(duì)于判斷曲線(xiàn)、曲面的光順性的光順準(zhǔn)則一般可分為兩類(lèi)：一類(lèi)是整體光順準(zhǔn)則，如在能量法中采用的光順準(zhǔn)則為“彈性模板應(yīng)變能較小”，而最小二乘法采用的是“剪力躍度較小”。另一類(lèi)是局部光順準(zhǔn)則，如根據(jù)曲面上的關(guān)鍵曲線(xiàn)（如 u，v方向的參數(shù)線(xiàn)，或者曲面與平行于坐標(biāo)平面的一系列平面的截面線(xiàn)等）是否光順，以及曲面的曲率（主曲率、高斯曲率、平均曲率等）的變化是否均勻等來(lái)判定。浙江大學(xué)的馬利莊和石教英給出了如下的曲面光順準(zhǔn)則[1]：

關(guān)鍵曲線(xiàn)（如飛機(jī)或船舶曲面的骨架線(xiàn)）光順；

網(wǎng)格線(xiàn)無(wú)多余拐點(diǎn)（或平點(diǎn)）及變撓點(diǎn)；

主曲率（低次曲面）在節(jié)點(diǎn)處的躍度和足夠??；

彈性薄板的應(yīng)變能較??；

高斯曲率變化均勻。

此外，在汽車(chē)工業(yè)中，還常用平行光照射到車(chē)身上來(lái)檢查曲面的光順性。對(duì)于帶有特征的曲線(xiàn)、曲面，如在復(fù)雜產(chǎn)品設(shè)計(jì)中需要設(shè)計(jì)出的曲面帶有邊界輪廓，表面棱線(xiàn)等特征，這些帶有關(guān)鍵信息的特征應(yīng)該在光順過(guò)程中被保留。本文在總結(jié)了曲面光順中的幾種常用方法的特點(diǎn)，針對(duì)這些方法在光順中保特征性的不足，結(jié)合各向異性小波在表達(dá)高維信息的優(yōu)勢(shì)，提出將其融入曲面的多分辨率分析，以達(dá)到保特征的目的。

1 基于能量的曲面光順?lè)椒?/h2>

1.1 最小二乘法

上世紀(jì) 60年代初，最小二乘法曾是最有影響的一種光順?lè)椒?。最小二乘法將樣條的剪力躍度平方和作為目標(biāo)函數(shù)進(jìn)行優(yōu)化，使其達(dá)到最小值來(lái)實(shí)現(xiàn)光順，因而可以使樣條在光順后的曲率變化比較均勻。朱心雄等人將其推廣到對(duì)空間參數(shù)三次B樣條曲線(xiàn)的光順[2]。此外，用拉格朗日乘數(shù)法對(duì)帶約束的光順目標(biāo)函數(shù)進(jìn)行優(yōu)化設(shè)計(jì)，可以得到保形的光順逼近曲線(xiàn)[3]。但是最小二乘法是整體光順?lè)?，需要將曲線(xiàn)、曲面的所有控制頂點(diǎn)作為未知量，計(jì)算量大、計(jì)算速度慢。

1.2 能量法

能量法的基本思想是將光順問(wèn)題轉(zhuǎn)化為包含型值點(diǎn)偏差和光順兩部分加權(quán)平均的目標(biāo)函數(shù)的極小化問(wèn)題。能量法的關(guān)鍵是能量函數(shù)（光順準(zhǔn)則）E(V)的選取，通常是通過(guò)對(duì)彈性薄板應(yīng)變能進(jìn)行簡(jiǎn)化或近似得到的。在選取能量函數(shù)時(shí)，需要考慮其是否具有明確的物理意義或幾何意義，以及實(shí)際問(wèn)題中對(duì)光順性的具體要求等。如針對(duì)參數(shù)曲面，Lott和Pullin[4]提出的一種薄板應(yīng)變能的簡(jiǎn)化形為曲面的主曲率）對(duì)曲面光順時(shí)會(huì)使曲面趨于平坦。Moreton[5]提出的采用

能量法具有很強(qiáng)的物理背景，光順效果好，但計(jì)算量大，計(jì)算速度慢。對(duì)此，穆國(guó)旺等[6]提出分片能量法，大大地減少了存儲(chǔ)量和計(jì)算量。此外，將能量法與選點(diǎn)法相結(jié)合進(jìn)行曲線(xiàn)曲面光順[7]，一次性修改所有對(duì)壞點(diǎn)有影響的控制頂點(diǎn)，減少了計(jì)算量，也使得曲線(xiàn)曲面的形狀改變盡可能小。將能量法與最小二乘法相結(jié)合[8]，在對(duì) B樣條曲線(xiàn)光順時(shí)可以使曲線(xiàn)的光順性及逼近性提高，使曲線(xiàn)的曲率變化更加均勻。

1.3 選點(diǎn)修改法

選點(diǎn)修改法是一種局部光順?lè)?，適用于少數(shù)型值點(diǎn)處不光順的情況。該方法充分利用NURBS曲面可以通過(guò)調(diào)整控制頂點(diǎn)來(lái)達(dá)到修改曲面形狀的特點(diǎn)，逐次找出不符合光順準(zhǔn)則的型值點(diǎn)（稱(chēng)為“壞點(diǎn)”），并依次對(duì)“壞點(diǎn)”進(jìn)行修改。特別是對(duì)于曲面上“壞點(diǎn)”的選取較為復(fù)雜，需借助曲面的等曲率線(xiàn)、等照度線(xiàn)或反射線(xiàn)等曲面分析工具，結(jié)合交互的方法來(lái)判斷曲面在哪些區(qū)域不光順。目前，“壞點(diǎn)”的修改方法主要有Kejellander方法和Farin的節(jié)點(diǎn)刪除與插入法、圓率法、基樣條方法等。選點(diǎn)修改法的優(yōu)點(diǎn)是計(jì)算速度快，修改能力強(qiáng)，具有局部性。但當(dāng)連續(xù)出現(xiàn)多個(gè)壞點(diǎn)時(shí)，光順效果往往比較差。胡海艷等[9]應(yīng)用廣義逆矩陣使節(jié)點(diǎn)消去，并通過(guò)升階和最小二乘逼近等步驟得到光順的目標(biāo)曲線(xiàn)。該方法減少了Farin光順?lè)椒ǖ挠?jì)算量，但不能直接用于樣條曲面的光順。將選點(diǎn)修改法與小波多分辨率分析方法結(jié)合，Martin Bertram提出的單節(jié)點(diǎn)小波[10]，刪除“壞點(diǎn)”后根據(jù)小波系數(shù)以重新插入控制頂點(diǎn)，可以在保存幾何細(xì)節(jié)的基礎(chǔ)上使曲面更逼近原曲面，也更光順。

2 升階光順?lè)?/h2>

升階[11]是 B樣條方法配套技術(shù)中一項(xiàng)重要技術(shù)，通過(guò)增加B樣條基函數(shù)的次數(shù)來(lái)提高曲面的光順性，使曲線(xiàn)、曲面更為光順。B樣條曲線(xiàn)升階的主要方法有：Prautzsch的升階方法、Cohen等的升階方法、Piegl等的升階方法、秦開(kāi)懷的升階方法[12]等。采用升階的方法，即是將B樣條函數(shù)的次數(shù)升高，特別是對(duì)于數(shù)據(jù)點(diǎn)個(gè)數(shù)變化較大的部分，得到的擬合曲線(xiàn)也更光順[13]，但加大了數(shù)據(jù)量。

3 保特征光順?lè)椒?/h2>

為滿(mǎn)足產(chǎn)品的設(shè)計(jì)意圖，要求所設(shè)計(jì)曲面的光順性不僅僅要美觀，同時(shí)也要滿(mǎn)足功能上的需求，這就需要在光順過(guò)程中，對(duì)曲面特征進(jìn)行保留，即采用保特征光順?lè)椒ā５壳皩?duì)這類(lèi)NURBS曲面的保特征光順問(wèn)題研究還很少。

另一種對(duì)曲面進(jìn)行光順的方式是將其轉(zhuǎn)化為網(wǎng)格線(xiàn)的光順性問(wèn)題來(lái)進(jìn)行處理。網(wǎng)格光順的目的是要從原始網(wǎng)格模型中去除噪聲，并在局部和整體的體積變形最小化的基礎(chǔ)上保持模型的幾何形狀特征。近年來(lái)，一些學(xué)者致力于這方面的研究，提出了不少好的保特征的網(wǎng)格光順?biāo)惴?。Klaus Hildebrandt等[14]融合平均曲率流與各向異性方法定義了各向異性平均曲率向量，并提出了去除噪聲保特征的光順?biāo)惴āＴ撍惴軌虮Ａ舨⒓訌?qiáng)曲面上的非線(xiàn)性幾何特征，如曲面上的彎曲部分、特征線(xiàn)等。Xu Guoliang[15]通過(guò)用有限元離散化的方法求解偏微分方程——曲面熱流方程(Heat equation），提出的保特征光順?biāo)惴芙鉀Q曲面光順不夠或過(guò)于光順而丟失細(xì)節(jié)的問(wèn)題。目前，這些保特征光順?lè)椒ㄖ饕杏诨谥髑?、平均曲率等幾何流的各向異性融合方法。這些方法雖然能夠增強(qiáng)曲面尖銳度，但是對(duì)多尺度形狀特征并存的曲面光順還需進(jìn)一步探討。

4 基于各向同性小波的曲面光順?biāo)惴?/h2>

近年來(lái)，小波分析這種新的時(shí)頻分析數(shù)學(xué)方法在NURBS曲線(xiàn)曲面光順中得到廣泛的運(yùn)用。利用小波的分解與重構(gòu)算法，可以使B樣條曲線(xiàn)曲面具有很好的多分辨特性。因此，可用于B樣條曲線(xiàn)曲面的幾何造型、光順處理以及數(shù)據(jù)壓縮等。采用不同的基函數(shù)對(duì)曲線(xiàn)曲面進(jìn)行多分辨率分析也會(huì)得到不同的結(jié)果。Finkelstein[16]提出了在閉區(qū)間上基于半正交 B樣條小波構(gòu)造的多分辨率曲線(xiàn)編輯系統(tǒng)。孫延奎、朱心雄等[17]給出了準(zhǔn)均勻雙二次 B樣條曲面小波分解與重構(gòu)的基本概念，描述了曲面小波分解與重構(gòu)的快速算法并論述了準(zhǔn)均勻 B樣條曲線(xiàn)曲面的多分辨分析理論。趙罡等[18]基于多分辨率技術(shù)，對(duì)曲線(xiàn)、曲面的光順和編輯進(jìn)行了研究，實(shí)現(xiàn)了在誤差控制下準(zhǔn)均勻B樣條曲線(xiàn)的光順。

但是準(zhǔn)均勻小波光順?biāo)惴ú荒芴幚磉吔缂s束的問(wèn)題，不適合進(jìn)行曲線(xiàn)曲面局部光順。對(duì)此，吳維勇[19]提出將能量法與小波光順?lè)ńY(jié)合，構(gòu)造分層能量算法，以達(dá)到靈活高效地進(jìn)行曲線(xiàn)局部光順。然而此算法無(wú)法直接應(yīng)用于任意B樣條曲線(xiàn)，需要先對(duì)原始曲線(xiàn)進(jìn)行準(zhǔn)均勻B樣條曲線(xiàn)逼近，這增加了計(jì)算量和計(jì)算誤差。此外，由于準(zhǔn)均勻小波光順?biāo)惴ㄒ蟊还忭?B樣條曲線(xiàn)的控制頂點(diǎn)數(shù)必須為2j+k（ j為整數(shù)，k為B樣條曲線(xiàn)的次數(shù)）。對(duì)此，可以直接插入節(jié)點(diǎn)來(lái)滿(mǎn)足小波光順對(duì)控制點(diǎn)個(gè)數(shù)的要求[20]。但是這個(gè)方法并沒(méi)有從本質(zhì)上做出改變。紀(jì)小剛等[21]根據(jù)多分辨分析的思想，將 m個(gè)控制點(diǎn)的曲線(xiàn)用 n(n

由于準(zhǔn)均勻 B樣條曲線(xiàn)曲面的邊界處采用了重節(jié)點(diǎn)技術(shù)，與定義域內(nèi)部的小波表達(dá)形式不一致，導(dǎo)致在小波重構(gòu)時(shí)計(jì)算量較大。對(duì)此，趙罡等[22]從幾何概念出發(fā)論述了基于小波的均勻三次 B樣條曲線(xiàn)曲面多分辨表示的原理及其實(shí)現(xiàn)。與準(zhǔn)均勻B樣條曲線(xiàn)曲面采用的半正交小波不同的是，均勻B樣條曲線(xiàn)曲面采用的是雙正交小波，這使其計(jì)算速度更快。但是其小波分解后的邊界控制和誤差控制還有待于進(jìn)一步研究。

為了進(jìn)一步減少小波方法的運(yùn)算量，出現(xiàn)了基于離散范數(shù)l2的雙正交非均勻B樣條小波基等的分解和重構(gòu)算法。這種雙正交小波只需要尺度函數(shù)和小波線(xiàn)性獨(dú)立，且重構(gòu)矩陣滿(mǎn)足=0，因此小波分解和重構(gòu)不需要積分運(yùn)算，計(jì)算較簡(jiǎn)便。根據(jù)這一特點(diǎn)，吳福鳴等[23]提出利用雙正交非均勻 B樣條小波對(duì)非均勻 B樣條曲面進(jìn)行光順的方法。但是這種光順?lè)椒ㄎ纯紤]曲面邊界的約束誤差精度的控制，較適用于對(duì)誤差精度要求不高的情況。

為了用小波實(shí)現(xiàn)對(duì)曲線(xiàn)、曲面的某一個(gè)局部進(jìn)行光順操作，Lyche[24]給出了非均勻小波的構(gòu)造。李奇敏等[25]將非均勻 B樣條小波分解用于NURBS曲面光順，通過(guò)整體光順度量確定B樣條曲面的光順?lè)较?，?shí)現(xiàn)曲面的整體光順。同時(shí)根據(jù)曲面在節(jié)點(diǎn)處的光順性選出曲面上的壞點(diǎn)，通過(guò)局部小波分解實(shí)現(xiàn)曲面的局部光順，并通過(guò)添加邊界約束的方式實(shí)現(xiàn)了曲線(xiàn)的保邊界光順。但對(duì)于曲面其算法只能保證單參數(shù)方向分解時(shí)邊界不發(fā)生變化，因此不具有真正的保邊界性。如圖1(a)所示為原始曲面，圖1(b)所示為經(jīng)過(guò)非均勻小波方法光順后所得到的曲面，圖1(c)所示是經(jīng)過(guò)升階光順后的曲面，可以看出光順的效果要差些。

圖1 不同方法對(duì)曲面進(jìn)行光順

此外，利用細(xì)分和提升模式構(gòu)造的第2代小波也具有計(jì)算簡(jiǎn)單，運(yùn)算速度快等特點(diǎn)。Martin[26]構(gòu)造了雙正交提升B樣條小波用于科學(xué)可視化。吳維勇在基于二元提升 B樣條小波給出了復(fù)雜拓?fù)淝嫘〔ü忭標(biāo)惴?，突破了能量法只能?duì)四邊域參數(shù)曲面進(jìn)行光順的限制，并指出適當(dāng)修改提升小波規(guī)則，可以保留尖角棱邊等不連續(xù)的特征[27]。

上述這些小波的光順?biāo)惴ㄖ饕m用于具有各向同性奇異性(Isotropic Singularities)的對(duì)象，而對(duì)于具有各向異性奇異性特征(Anisotropic Singularities)的曲面，上述小波光順?biāo)惴ň皇呛芎玫墓ぞ?。因?yàn)楦飨蛲孕〔ㄊ强煞蛛x的和各向同性的，難以有效地表達(dá)圖像、圖形中帶有關(guān)鍵信息的邊緣、曲線(xiàn)等更高維的特征，也增加了圖像、圖形表達(dá)的冗余度和處理時(shí)間。如圖2(a)所示為依次帶有s形、矩形、三角形類(lèi)特征的不光順曲面，圖2(b)所示為經(jīng)過(guò)小波光順后的曲面，其中黑色部分為光順中的改變，可以看出為各特征已經(jīng)沒(méi)有了原有的形狀。

圖2 帶特征曲面的光順

5 基于各向異性小波的 NURBS曲面光順

5.1 各向異性小波

各向異性小波理論是近幾年出現(xiàn)的一種新的數(shù)學(xué)分析方法，它為圖像提供了新的稀疏表示方法，能夠有效地捕捉信號(hào)中的高維奇異特征，從而保留幾何結(jié)構(gòu)等重要信息。各向異性小波包括 Ridgelet、Curvelet、Contourlet等。

1998至1999年，D.L Donoho和E.J.Candes提出了脊波變換(Ridgelet Transform)[28]。其本質(zhì)上是通過(guò)對(duì)小波基函數(shù)添加一個(gè)表征方向的參數(shù)得到的，因此具有很強(qiáng)的方向選擇和辨識(shí)的能力，可以非常有效地表示信號(hào)中具有方向性的奇異特征，如圖像的線(xiàn)性輪廓等。為處理以奇異邊界為曲線(xiàn)的對(duì)象，Candes在脊波理論的基礎(chǔ)上，利用剖分的方法構(gòu)造了單尺度脊波（Monoscale Ridgelet）。但是，其獲得好的逼近性能也是以計(jì)算復(fù)雜度增加為代價(jià)的。

接著，Candes和Donoho在1999年提出第1代曲波變換[29]（Curvelet transform），可以在所有可能的尺度S≥0上進(jìn)行分解。但是，第1代曲波的數(shù)字實(shí)現(xiàn)比較復(fù)雜，需要子帶分解、平滑分塊、正規(guī)化和脊波分析等一系列步驟，其金字塔的分解也帶來(lái)了巨大的數(shù)據(jù)冗余量。因此Candes等人又提出了實(shí)現(xiàn)更簡(jiǎn)單、更便于理解的快速曲波變換算法，即第2代曲波。

在繼承了Curvelet變換的各向異性尺度關(guān)系的基礎(chǔ)上，Donoho和Vetterli于 2002年提出了Contourlet變換[30]，也稱(chēng)為金字塔方向?yàn)V波器組(Pyramid Direction Filter Bank, PDFB)。Contourlet變換的優(yōu)點(diǎn)在于能夠在使用少量系數(shù)的基礎(chǔ)上有效地表示平滑輪廓，而平滑輪廓正是自然圖像中的重要特征。目前，Contourlet變換已經(jīng)在圖像重建、壓縮、融合等領(lǐng)域有了初步的應(yīng)用。

圖3 基于Contourlet變換和NURBS的重構(gòu)曲面

參考Contourlet在圖像處理中的應(yīng)用，可以把Contourlet用于曲面處理中。我們把準(zhǔn)均勻B樣條曲面控制頂點(diǎn)坐標(biāo)做為變換對(duì)象，對(duì)其進(jìn)行Contourlet變換。圖3(a)是所用原始曲面，圖3(b)是對(duì)控制頂點(diǎn)坐標(biāo)利用小波變換后重構(gòu)得到的曲面，后者邊緣和隆起部分的誤差都較明顯?？梢钥闯?，對(duì)于這種特征比較突出的曲面，小波重構(gòu)方法存在不足。圖3(c)是基于Contourlet變換的NURBS重構(gòu)曲面。與圖3(b)對(duì)比，前者誤差為 0.0046，后者為 0.0034。實(shí)驗(yàn)結(jié)果表明，Contourlet變換的擬合精度要優(yōu)于小波變換，具有較高的實(shí)際應(yīng)用價(jià)值。

綜上所述，這些各向異性小波的出現(xiàn)，為提高曲面多分辨率表達(dá)的效率，減少表達(dá)的冗余，加快曲面顯示、處理的速度等提供了新的思路和方法。

5.2 基于各向異性小波的NURBS曲面光順

在復(fù)雜產(chǎn)品設(shè)計(jì)中，多尺度特征并存是常見(jiàn)的現(xiàn)象。多尺度特征并存曲面的保特征光順是多尺度特征并存曲面建模中的關(guān)鍵。而其目的就是要得到滿(mǎn)足設(shè)計(jì)要求的曲面模型，這就需要解決尺度分析、微分幾何屬性分布規(guī)律提取、特征提取、連續(xù)性、約束優(yōu)化等問(wèn)題。其中微分幾何屬性分布規(guī)律的數(shù)學(xué)描述、全局約束優(yōu)化方程的建立及優(yōu)化算法是實(shí)現(xiàn)曲面保特征光順的難點(diǎn)。目前關(guān)于多尺度特征并存曲面的主要處理方法是把基曲面和各個(gè)局部特征曲面之間通過(guò)剪裁、過(guò)渡和搭接等方法實(shí)現(xiàn)曲面從整體到局部的造型，這樣的造型過(guò)程對(duì)設(shè)計(jì)者要求非常高，而且在很多情況下得不到滿(mǎn)意的結(jié)果。

隨著多分辨率分析方法在圖形學(xué)中的廣泛應(yīng)用，基于小波的多分辨率表達(dá)已經(jīng)被應(yīng)用到參數(shù)曲線(xiàn)，并被推廣到張量基曲面、任意拓?fù)漕?lèi)型曲面、球面及體數(shù)據(jù)。但是已有的多分辨率曲面表達(dá)方法對(duì)多尺度特征并存曲面及其上有曲線(xiàn)類(lèi)特征的曲面表達(dá)能力也不足，相應(yīng)的會(huì)帶來(lái)曲面顯示、壓縮及傳輸方面的數(shù)據(jù)冗余。

基于這種小波多分辨率分析方法，李奇敏提出了平面曲線(xiàn)多尺度分析方法，定義了曲線(xiàn)的局部形狀尺度因子，結(jié)合去除細(xì)節(jié)小波重建方法，分析了曲線(xiàn)上各個(gè)尺度特征的曲率主要分布規(guī)律，實(shí)現(xiàn)了大尺度特征在較低分辨率下分析，小尺度特征在較高分辨率下分析，并利用分析結(jié)果指導(dǎo)曲線(xiàn)的光順和特征提取[31]。但是曲面的情況要比平面曲線(xiàn)復(fù)雜得多，在分析時(shí)既要考慮特征的大小還要考慮特征的方向問(wèn)題。而基于各向同性小波的多分辨率分析技術(shù)對(duì)特征方向的表達(dá)能力有限，為了提取多尺度特征并存曲面各尺度特征微分幾何屬性分布規(guī)律，有必要把基于各向異性小波的多分辨率分析方法和多尺度分析方法結(jié)合起來(lái)，對(duì)曲面各局部特征主要微分幾何信息分布規(guī)律進(jìn)行分析。

如上所述，要實(shí)現(xiàn)多尺度特征并存曲面保特征光順，首先需要建立基于各向異性小波的曲面多分辨率表達(dá)和分析框架，用各向異性小波基代替各向同性小波基，把分析相關(guān)的Fourier變換、Radon變換等運(yùn)算推廣到張量積曲面和離散曲面，實(shí)現(xiàn)曲面的各項(xiàng)異性小波多分辨率分析。

其次，在曲面光順之前對(duì)曲面形狀即設(shè)計(jì)意圖進(jìn)行分析，提取各局部特征尺度，然后基于各向異性小波的多分辨率分析方法對(duì)曲面進(jìn)行多尺度分析，確定各局部特征存在的尺度及其微分幾何屬性分布規(guī)律，然后用數(shù)學(xué)的語(yǔ)言對(duì)該分布規(guī)律進(jìn)行描述，在所得微分幾何屬性分析基礎(chǔ)上提取特征區(qū)域的邊界、曲線(xiàn)類(lèi)特征的脊線(xiàn)等特征曲線(xiàn)，并分析各特征間其連續(xù)性，以所提取特征曲線(xiàn)和微分幾何屬性分布規(guī)律為約束，利用帶約束的全局優(yōu)化算法實(shí)現(xiàn)曲面保形光順。

6 結(jié)束語(yǔ)

本文在總結(jié)分析了NURBS曲線(xiàn)曲面各種光順?biāo)惴ǖ幕A(chǔ)上，詳細(xì)介紹了小波方法并與其他光順?lè)椒ㄟM(jìn)行了對(duì)比。我們指出了小波光順?biāo)惴ㄔ诙喑叨忍卣鞑⒋媲姹Ｌ卣鞴忭樕系牟蛔悖岢隽藢⒏飨虍愋孕〔ㄈ谌肭娴亩喾直媛史治鲋械乃枷牒蛯?shí)現(xiàn)這種方法的思路。我們將繼續(xù)對(duì)該問(wèn)題進(jìn)行更詳細(xì)的研究。

[1]馬利莊, 石教英. 曲線(xiàn)曲面的幾何光順?biāo)惴╗J]. 計(jì)算機(jī)學(xué)報(bào), 1996, 19: 210-216.

[2]朱心雄. 自由曲線(xiàn)曲面造型技術(shù)[M]. 北京:科學(xué)出版社, 2000: 358-360.

[3]蔣大為, 李安平. B樣條曲線(xiàn)的最小二乘保形光順逼近[J]. 工程數(shù)學(xué)學(xué)報(bào), 2000, 17(1): 125-128.

[4]Lott N J, Pullin D I. Method for fairing B-spline surfaces [J]. Computer-Aided Design, 1988, 20(10):597- 604.

[5]Moreton H P, Sequin C H. Functional optimization for fair surfaces design [J]. Computer Graphics, 1992,26(2): 167-175.

[6]穆國(guó)旺, 朱心雄, 涂侯杰, 等. 曲面光順的分片能量法[J]. 工程圖學(xué)學(xué)報(bào), 1999, (1): 29-34.

[7]Zhang Caiming, Zhang Pifu, Cheng Fuhua. Fairing spline curves and surfaces by minimizing energy [J].Computer-Aided Design, 2001, 33(13): 913-923.

[8]王瑩瑩, 劉德平, 陳建軍, 等. 基于曲率均化的B樣條曲線(xiàn)能量光順?lè)椒╗J]. 機(jī)床與液壓, 2010, 38(3):94- 96.

[9]胡海艷, 吳振遠(yuǎn). 基于廣義逆矩陣的 8樣條曲線(xiàn)的節(jié)點(diǎn)消去與光順[J]. 工程地球物理學(xué)報(bào), 2010, 7(2):106- 110.

[10]Bertram M. Single-knot wavelets for non-uniform B-splines [J]. Computer Aided Geometric Design,2005, 22(9): 849-864.

[11]施法中. 計(jì)算機(jī)輔助幾何設(shè)計(jì)與非均勻有理 B樣條[M]. 北京: 高等教育出版社, 2001: 285-293.

[12]秦開(kāi)懷. 非均勻B樣條曲線(xiàn)升階的新算法[J]. 計(jì)算機(jī)學(xué)報(bào), 1996, 19(7): 537-542.

[13]曹利新, 游 洪. 基于升階的 B 樣條曲面光順擬合[J]. 機(jī)械制造, 2003, 41(7): 23-25.

[14]Hildebrandt K, Polthier K. Anisotropic filtering of non-linear surface features [C]//Eurographics 2004 Comput Graph Forum, 2004: 391-400.

[15]Xu Guoliang. Surface fairing and featuring by mean curvature motions [J]. Journal of Computational and Applied Mathematics, 2004, 163(1): 295-309.

[16]Finkelstein D H, Salesin. Multiresolution curves [C]//Proceedings of SIGGRAPH. ACM, 1994: 261-268.

[17]孫延奎, 朱心雄, 唐 龍, 等. 準(zhǔn)均勻 B 樣條曲面小波分解的快速算法[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2001(z1): 209-213.

[18]趙 罡, 穆國(guó)旺, 朱心雄. 基于小波的準(zhǔn)均勻 B樣條曲線(xiàn)曲面變分造型[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào), 2002, (1): 61-65.

[19]吳維勇, 王小椿. 自由曲線(xiàn)局部光順的分層能量算法[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào), 2002, 14(10):959-966.

[20]潘洋宇, 姜福祥. 任意控制點(diǎn)曲線(xiàn)小波光順?lè)椒ㄑ芯縖J]. 機(jī)械設(shè)計(jì), 2009, 26(11): 12-13.

[21]紀(jì)小剛, 龔光容. 準(zhǔn)均勻三次B樣條曲線(xiàn)的任意分辨率光順[J]. 機(jī)械科學(xué)與技術(shù), 2009, 28(11): 1477-1481.

[22]趙 罡, 穆國(guó)旺, 閆光榮, 等. 均勻 B 樣條曲線(xiàn)曲面的小波表示[J]. 工程圖學(xué)學(xué)報(bào), 2001, (1): 80-87.

[23]吳福鳴, 潘日晶. 基于雙正交非均勻B樣條小波的曲面光順[J]. 福建師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2010, (1): 33-37.

[24]Lyche T, Mrken K. Spline-wavelets of minimal support [J]. Numerical Methods of Approximation Theory, 1992, (9): 177-194.

[25]李奇敏, 柯映林, 何玉林. 基于非均勻B樣條小波的NURBS曲面光順[J]. 中國(guó)機(jī)械工程, 2007, 18(5):577-582.

[26]Bertram M, Mark A D, Hamann B, el at. Bicubic subdivision-surface wavelets for large-scale isosurface representa-tionand visualization [C]//Proceedings of IEEE Visualization, 2000:389-396&579.

[27]吳維勇. 自由曲線(xiàn)曲面多分辨率造型技術(shù)研究與應(yīng)用[D]. 西安: 西安交通大學(xué), 2003.

[28]Emmanuel J, Candès, David L, et al. Ridgelets:a key to higher-dimensional intermittency [C]//Phil. Trans.R. Soc. Lond. A. 1999: 2495-2509.

[29]焦李成, 譚 山, 劉 芳. 脊波理論: 從脊波變換到 Curvelet變換[J]. 工程數(shù)學(xué)學(xué)報(bào), 2005, 22(5):761-773.

[30]Lounsbery M, De R T, Warren J. Multiresolution analysis for surfaces of arbitrary topological type [J].ACM Transaction on Graphics, 1997, 16(1): 34-73.

[31]李奇敏. 小波在反求工程中的若干應(yīng)用[D]. 杭州:浙江大學(xué), 2006.