[發(fā)明專利]一種基于三角形域上L曲面和W曲面的遞歸曲面構(gòu)造方法

無效

• 申請?zhí)枺?/td>: 200710031024.3
• 申請日: 2007-10-24
• 公開（公告）號: CN101149839A
• 公開（公告）日: 2008-03-26
• 發(fā)明（設(shè)計(jì)）人: 俞晏林;曹穎;羅笑南
• 申請（專利權(quán)）人: 中山大學(xué)
• 主分類號: G06T11/20
• 分類號: G06T11/20;G06T17/30;G06F17/50
• 代理公司: 廣州市深研專利事務(wù)所
• 代理人: 陳雅平

1.一種基于三角形域上L曲面和W曲面的遞歸曲面構(gòu)造方法，其特征在于，它的主要步驟包括：

1)首先定義三角參數(shù)域上的一般遞歸形式曲面表示方法：

定義符號：J₀＝(1，0，0)，J₁＝(0，1，0)，J₂＝(0，0，1)，$a=(a_0,a_1,a_2)\∈Z_{+}^3,$設(shè){P_a||a|＝n}為三角域上特征網(wǎng)的頂點(diǎn)，$a\∈Z_{+}^3,$L_a^T((x)，M_a^T(x)，N_a^T(x)為三角域上相應(yīng)節(jié)點(diǎn)處的調(diào)配函數(shù)，具體要求如下：

L_a^T(x)應(yīng)該滿足1≤a₀≤n，0≤a₁，a₂≤n；

M_a^T(x)應(yīng)該滿足1≤a₀≤n，0≤a₁，a₂≤n；

N_a^T(x)應(yīng)該滿足1≤a₀≤n，0≤a₁，a₂≤n；

且$L_{a+J_0}^T\left(x\right)+M_{a+J_1}^T\left(x\right)+N_{a+J_2}^T\left(x\right)\≡1,r=1,\mathrm{\Λ},n,\left|a\right|=n-r;$

定義三角參數(shù)域上一般遞歸形式曲面的表示方法為：

$\left\{\begin{array}{cc}{d}_a^0\left(x\right)=P_a& \left|a\right|=n\\ d_a^m\left(x\right)=L_a^{n-m+1}\left(x\right)d_{a-J_0}^{m-1}+M_a^{n-m+1}\left(x\right)d_{a-J_1}^{m-1}+N_a^{n-m+1}\left(x\right)d_{a-J_2}^{m-1},m=1,\mathrm{\Λ},n& \left|a\right|=m\end{array}\right.$

其中若a_i-1＜0或者＞0，則令$d_{a-J_i}^{m-1}=0,i=\mathrm{0,1,2};$

2)定義三角域上的n次L曲面

設(shè)：$L_a^r\left(x\right)=a_a^{r}x+b_a^r,$$M_a^r\left(x\right)=c_a^{r}x+d_a^r,$$N_a^r\left(x\right)=e_a^{r}x+f_a^r$

(且滿足$a_a^r\≠0,$$c_a^r\≠0,$$e_a^r\≠0;$r＝1，2，Λn，|a|＝n-r)

根據(jù)1)定義三角域上的n次L曲面為：n次遞歸曲線必須滿足

$\left\{\begin{array}{c}\frac{b_a^r}{a_a^{r+1}}=\frac{b_a^{r+1}}{a_a^{r+1}}\\ \frac{1-b_a^{r+1}}{a_a^{r+1}}=\frac{1-b_{a+J_0}^r}{b_{a+J_0}^r}\end{array}\right.,r=\mathrm{1,2},\mathrm{\Λn};\left|a\right|=n-r$

$\left\{\begin{array}{c}\frac{d_a^r}{c_a^r}=\frac{d_a^{r+1}}{c_a^{r+1}}\\ \frac{1-d_{a+J_1}^r}{c_{a+J_1}^r}=\frac{1-d_{a_0}^{r+1}}{c_{a_0}^{r+1}}\end{array}\right.,r=\mathrm{1,2},\mathrm{\Λn};\left|a\right|=n-r$

$\left\{\begin{array}{c}\frac{f_a^r}{e_a^{r+1}}=\frac{f_a^{r+1}}{e_a^{r+1}}\\ \frac{1-f_{a+J_1}^r}{e_a^{r+1}}=\frac{1-f_a^{r+1}}{e_a^{r+1}}\end{array}\right.,r=\mathrm{1,2},\mathrm{\Λn};\left|a\right|=n-r$

則為三角域上的n次L曲面；

3)定義W曲面

若n次L曲面繼續(xù)滿足

$0\≤L_a^r\left(x\right),M_a^r\left(x\right),N_a^r\left(x\right)\≤1,r=\mathrm{1,2},\mathrm{\Λn};\left|a\right|=n-r$

則為n次W曲面；

4)轉(zhuǎn)換成Bernstein-Bezier曲面

令：

$\left\{\begin{array}{c}{L}_a^r\left(x\right)=u,\\ M_a^r\left(x\right)=v,\\ N_a^r\left(x\right)=w\end{array}\right.r=\mathrm{1,2},\mathrm{\Λ},n;\left|a\right|=n-r$

則，得到Bernstein-Bezier曲面

$P(u,v,w)=\underset{i+j+k=n}{\mathrm{\Σ}}{P}_{\mathrm{ijk}}{B}_{\mathrm{ijk}}^n(u,v,w),(u,v,w\≥0,u+v+w=1)$

其中$B_{\mathrm{ijk}}^n(u,v,w)=\frac{n!}{i!j!k!}{u}^i{v}^j{w}^k,(i+j+k=n,u+v+w=1),$P_ijk為三角形域網(wǎng)的節(jié)點(diǎn)；

5)三角域上Bezier曲面連續(xù)拼接

如果給定兩片三角域上n次Bezier曲面：

$\left\{\begin{array}{c}P(u,v,w)=\underset{i+j+k}{\mathrm{\Σ}}{P}_{\mathrm{ijk}}(u,v,w)\\ Q(u,v,w)=\underset{i+j+k=n}{\mathrm{\Σ}}{Q}_{\mathrm{ijk}}(u,v,w)\end{array}\right.$

則Q(u，v，w)與P(u，v，w)在拼接線上達(dá)到幾何連續(xù)的條件為：

$\left\{\begin{array}{cc}{Q}_{\mathrm{ij}0}=P_{\mathrm{ij}0},& i+j=0\\ Q_{\mathrm{ij}1}=a_1{P}_{(i+1)j0}+a_2{P}_{i(j+1)0}+a_3{P}_{\mathrm{ij}1},& i+j=n-1;a_1+a_2+a_3=1\end{array}\right.$

然后依次連接所有劃分的三角域曲面即得到所要的曲面構(gòu)造。