有理参数曲线的C-1分段根式弧角重新参数化算法研究_顾特.pdf
《有理参数曲线的C-1分段根式弧角重新参数化算法研究_顾特.pdf》由会员分享,可在线阅读,更多相关《有理参数曲线的C-1分段根式弧角重新参数化算法研究_顾特.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 22卷 第 2期2023年 2月Vol.22 No.2Feb.2023软 件 导 刊Software Guide有理参数曲线的C1分段根式弧角重新参数化算法研究顾特(广西民族大学 人工智能学院,广西 南宁530006)摘要:曲线和曲面的参数化与重新参数化是计算机辅助几何设计研究中的基本问题。针对一类角速度函数在单位区间上存在零点的有理曲线参数表示,研究其近似弧角重新参数化,通过分段根式变换复合分段Mbius变换并结合C1连续的弧角新参数化算法,提出一种基于C1分段根式变换的重新参数化算法。实验表明,该算法能极大提升有理参数表示的角速度均匀度。关键词:有理曲线参数表示;近似弧角重新参数化;C
2、1分段根式变换DOI:10.11907/rjdk.221239开 放 科 学(资 源 服 务)标 识 码(OSID):中图分类号:TP301.6 文献标识码:A文章编号:1672-7800(2023)002-0100-08Research on C1 Piecewise Radical Arc Angle Re-parameterization Algorithm for Rational Parametric CurvesGU Te(School of Artificial Intelligence,Guangxi University for Nationalities,Nanning 5
3、30006,China)Abstract:Parameterization and re-parameterization of curves and surfaces is a basic problem in computer-aided geometric design.The approximate arc angle re-parameterization is studied for a class of rational curve parameter representation of angular velocity function with zero on unit in
4、terval.Based on the piecewise radical transformation and the piecewise Mbius transformation combined with the C1 arc angle new parameterization algorithm,a re-parameterization algorithm based on the C1 piecewise radical transformation was proposed.Experiments show that this algorithm can greatly imp
5、rove the angular velocity uniformity represented by the rational parameters.Key Words:rational parameter representation;approximate arc angle re-parameterization;C1 piecewise radical transformation0 引 言曲线曲面的参数表示常被应用于计算机图形学与计算机辅助几何设计领域,其优点为表达方便、易于显示。对于参数曲线而言,有些参数具有几何意义,主要包括弧长参数、弦长参数和弧角参数等。其中,曲线的弧长参数能
6、够均匀划分曲线的长度,在实际生活中有广泛应用1-4;弧角参数则能更好地反映曲线在各个位置上曲率的变化。尽管所体现的性质有所不同,但是在研究弧长和弧角参数的过程中所采用的思想有很多相似之处。在几何设计和造型中,有理参数表示具有简单、直观、易于实现的优点,因此研究有理弧长参数和有理弧角参数具有重要意义。关于有理弧长参数,Gil等5证明了任意平面曲线都具有弧长参数,但是这些参数表达式不一定是有理的;Farouki等6证明了在平面曲线中,只有直线具有有理弧长参数化;Salkkalis等7证明了在Rn空间中,仅有直线具有有理弧长参数。因此,求得一般曲线的有理弧长参数表示几乎是不可能的。于是研究者们转而寻
7、找一般曲线的有理近似弧长参数。例如,Farouki8提出一种当次数不变的情况下,基于Mbius参数变换计算B线条的最优参数化方法;Costantini等9提出用一种分段的Mbius变换近似弧长参数化变换,并且证明了分段区间越小分段越多,得到的参数变换越接近于弧长参数变换;Liang等10提出一种C1连续的分段有理在参数化算法,该算法是基于实验结果实现的。收稿日期:2022-03-06基金项目:广西科技计划项目(GUIKE-AD18126010)作者简介:顾特(1993-),男,广西民族大学人工智能学院硕士研究生,研究方向为高性能计算。第 2 期顾特:有理参数曲线的C1分段根式弧角重新参数化算法
8、研究关于有理弧角参数,最早由Patterson等11提出一种曲率参数,该参数是由参数曲线的弧长参数和曲率公式推导得到的。利用曲率参数作图能将更多的点集中在曲率更大的地方,从而反映曲线在各个位置的弯曲程度。与弧长参数表示相似,一般来说曲率参数没有解析表达式,因此对一般曲线无法精确求出其曲率参数表示。为此,Patterson等给出了一种数值近似曲率参数化的算法,但该算法高度依赖于作图点的数目,如果作图点数目发生改变,则必须重新计算其曲率参数化,因而效率较低。此后,Yang等12-16提出一套曲线均匀弧角参数化理论,该理论首先定义了一个角速度均匀度函数,用衡量曲线的某种参数表示接近弧角参数表示的程度
9、,同时指出如果一条参数曲线某一参数表达式的角速度是均匀的(也就是角速度函数等于不为0的常数),那么随着参数变化曲线的弯曲程度也是均匀的,这也就是标准弧角参数,其与曲率参数在作图时有着相同的特性。同样,弧角参数也难以求得有理的解析表达式。因此,该理论提出了计算曲线有理近似弧角参数表示的多种算法,其主要思想是通过分段的Mbius变换近似非有理的弧角参数变换。此后,刘振华等17根据曲线均匀弧角参数化理论设计出软件包ImUp+,这个软件包可以对原来曲线参数进行优化而得到新的参数,利用新参数绘制的曲线比原来的曲线更加光滑。然而,如果给出的参数角速度函数含零点,这些重新参数化的算法无法很好地均匀化零点附近
10、的角速度,这也是为什么ImUp+17要求处理的参数角速度函数不含有零点。为了完善曲线均匀弧角参数化理论,优化角速度函数含有零点参数表达式的角速度均匀度,本文针对角速度包含零点的参数曲线,考虑其 C1连续的根式近似弧角参数化。通过引入C1连续性的约束计算得到角速度不包含零点的C1连续参数表示,然后使用文献 16 中C1连续的优化弧角重新参数化算法,对已经得到的参数表示进行有理的重新参数化,以得到曲线的近似弧角参数表示。计算实例表明本文方法可以在C1连续性的约束条件下有效解决角速度包含零点的弧角重新参数化问题。1 研究背景本章阐述关于角速度均匀的变量用于衡量近似弧角参数接近弧角参数的程度,给出分段
11、Mbius变换的定义,其将在 2.1 节求取曲线角速度函数不含零点的参数以及2.2节优化角速度均匀度时被用到。对平面曲线的某一参数表示为:p=p(x(t),y(t)(1)其角速度函数为:p(t)=x(t)y(t)-x(t)y(t)2 x(t)2+y(t)2(2)定义其角速度均匀度为15:up=11+2p2p=2p012p(t)dt(3)其中:p=01p(t)dt,p2=01(p(t)-p)2dt(4)设t=r(s)为 0,1 上的单调变换,将p与r复合可以得到原曲线的一个新的参数表示p r,该过程称为重新参数化,可以证明重新参数化后的角速度函数为15:p r(s)=p(r(s)r(s)(5)其
12、角速度均匀度为:p r(s)=p(r(s)r(s)(6)定义 113 若T=t0,.,ti,.,tN,S=s0,.,si,.,sN,=1,.,i,.,N,其 中ti-1 ti,si-1 si,t0=s0=0,tN=sN=1,0 i 0。式(10)是一个微分方程,通过求解该方程,即可得到想要的变换(z)。另一方面,参数表示的角速度函数可能有多个零点,而式(10)所得到的结果只针对其中某个零点,为了能够处理每一个零点,考虑(z)为分段变换。通过调整式(10)解的形式,可以使得(z)在 ti-1012023 年软 件 导 刊,ti+上保持连续,从而p (z)在 ti-,ti+上连续,但是(z)在区间
13、端点ti-和ti+由于左右极限不连续,p (z)将出现跳跃间断点。为了使得变换后角速度函数在0,1上连续,再对z做分段Mbius变换mc,使得变换后p mc(s)在0,1连续。定理1 设p为曲线的一个参数表示,p的角速度函数在(0,1)仅有N个零点t1,.,ti,.,tN,i为零点ti的重数,那么参数表达式p mc的角速度函数在(0,1)上连续且不含零点,其中为分段根式变换,表示为:(z)=i(z)z ti-,ti+z,其他 (11)其中:i(z)=-(-iz+iti)1i+1+ti,z t1-,ti(iz-iti)1i+1+ti,z ti,ti+(12)此处,0 0,而p (z)=p(z)(
14、z),当t ti,p(t)0,因此对除zi=-1(ti)以外的点,p (z)0。现 在 考 虑zi=-1(ti)的 情 况,当t ti-,ti+时,p(t)可以表示为:p(t)=|t-ti|i(t),(t)0(17)则:p (zi)=p(zi)(zi)=|-diz+diti|ii+1di|-diz+diti|-ii+1i+1(ti)=di(ti)i+1 0(18)综上可知,p (z)在(0,1)上不含零点。以下证明p mc(s)在(0,1)上连续且不含零点,首先证 明p mc(s)不 含 零 点。由 于mi(s)=(1-i)i(1-i)s?+(1-s?)i2 0,故p mc=p (mc)mc
15、0。下 证p mc(s)连 续:由 于p mc=p(mc)(mc)mc,而和mc分段点相同(即ti-和ti+),且在分段内是连续的,因此p mc在除端点ti-和ti+外连续。下面考虑p mc在端点ti-的连续性:lims (ti-)+(z(s)mc(s)=i(-i(ti-)+iti)-ii+1i+1 1-2i2i=1i+11-2i2i(19)lims (ti-)-(z(s)mc(s)=1 2i-11-2i-1=2i-11-2i-1因为2i=12i-11-2i-1(i+1)+1,所以:lims (ti-d)-(z(s)mc(s)=lims (ti-d)+(z(s)mc(s)=(z(si)mc(s
16、i)(20)故p mc在端点ti-处连续,同理可证p mc在端点ti+处也连续,因此p mc在(0,1)上没有零点且连续。利用定理 1中的变换对于角速度函数含有零点的参 102第 2 期顾特:有理参数曲线的C1分段根式弧角重新参数化算法研究数表示重新参数化后,可以得到角速度函数不含零点且连续的参数表示,这一过程称为去零点。2.2重新参数化后的角速度均匀度优化2.1节中去零点的过程并未对参数表示的角速度均匀度进行优化,以下将改进文献 16 中的变换,对定理 1中得到的参数表示p*=p mc进行重新参数化,继续对角速度均匀度进行优化,从而得到原参数曲线的近似弧角参数化。定理216 若p为曲线的参数
17、表达式,p在 0,1 上C1连续,mo是由序列S,R,确定的分段Mbius变换,其中S、R表示为:S=sj:p(sj)=0 sj=0 sj=1,sj-1 up,且p mo在 0,1 上连续。定理2要求p在 0,1 上C1连续,而p*仅在分段内满足C1连续的条件,在分段点ti-和ti+处则为C0连续。若要使用定理2对p*进行优化的重新参数化,则需对连续条件进行弱化。定理2中选取分段点的目的是将p拆分为单调的分段,虽然p*在分段点处不连续,但仍可以根据其左右导数的符号确定是否将该分段点加入S中。即:若lims (ti-)+p*(s)lims (ti-)-p*(s)0,则 将ti-并 入S;若lim
18、s (ti+)+p*(s)lims (ti+)-p*(s)0成立,则将ti+并入S。在定理2中式(21)用于确认分段点,保证分段点内部角速度函数单调性一致,式(22)是求得基于式(21)的分段Mbius变换最优的。在保证分段内角速度函数单调性一致的条件下,分段数越多,定理2给出的分段Mbius变换能更好地近似弧角参数变换。文献 16 还给出了一种加细划分的算法,对式(21)给出的分段做进一步划分,这种划分方式是基于式(23)对分段 si-1,si进行求解:p*(si-si-1)is-si-1si-si-1+(1-i)(1-s-si-1si-si-1)-2p*(1-2i)=0(23)将得到的新解
19、并入S得到新的分段S1,令S=S1,重复定理2的过程就能得到分段数更多的分段Mbius变换,反复执行这一过程就能得到近似弧角参数变换m*。由于p*在分段点处不连续,若要使用式(23)进行加细划分,也需要对连续条件进行弱化,令:F(s)=p*(si-si-1)is-si-1si-si-1+(1-i)(1-s-si-1si-si-1)-2p*(1-2i)(24)若lims (ti-)+F(s)lims (ti-)-F(s)0则ti+并 入S;若lims (ti+)+F(s)lims (ti+)-F(s)0成立,则将ti+并入S。2.3基于分段根式变换的近似弧角参数变换m*可 与mc复 合 为 一
20、分 段 Mbius 变 换mf,即mf=mc m*。最后得到*=mf,即为C1连续的近似弧角参数变换。定 理 3 设m1为 分 段 Mbius 变 换,其 参 数 序 列 为T1,S1,1;m2为相容的分段 Mbius 变换,其参数序列为S2,R2,2,变换m3=m1 m2,那么m3也是分段 Mbius 变换且由T3,R3,3确定,其中:T1=t10,t11,.,t1n,S1=s10,s11,.,s1n,1=11,12,.,1nS2=s20,s21,.,s2m,R2=r20,r21,.,r1m,2=21,22,.,2ms10=s20=r10=r20=0,s1n=s2n=r1m=r2m=1T3=
21、s3 k:t3 k=t1 i t3 k=m-11 i(s2 j),s1 i-1 s2 j s1 i+1,t3 k t3 k+1,k=0.m+n+1 R3=r3 k:r3 k=r2 j r3 k=m2 j(s1 i),s2 j-1 s1 i s2 j+1,r3 k r3 k+1,k=0.m+n+13=31,32,.,3 m+n+1(25)当 t3 k-1,t3 k t1 i-1,t1 i,r3 k-1,r3 k r2 j-1,r2 j时:3 k=(2s2 j-1+2s2 j-2s1 i-1-2s1 i)t3 k+(-t1 i-1-t1 i)s2 j+2s1 i-1t1 i-s2 j-1t1 i-
22、t1 i-1(s2 j-1-2s1 i)2 j+(s1 i-2s2 j+s1 i-1)t3 k+(t1 i-1+t1 i)s2 j-t1 i-1s1 i-s1 i-1t1 i)1 i-(s2 j+s2 j-1-2s1 i)2 j-s2 j+s1 i-1)(t3 k-t1 i)/(2s2 j-1+2s2 j-2s1 i-1-2s1 i)t3 k+(-2t1 i-1+2t3 k-1-2t1 i)s2 j+(-2t3 k-1+4t1 i)s1 i-1-2s2 j-1t1 i+(2s2 j-1-2s1 i)t3 k-1-2t1 i-1(s2 j-1-2s1 i-1)2 j+(s1 i-2s2 j+s1
23、 i-1)t3 k+(2t1 i-1-2t3 k-1+2t1 i)s2 j+(t3 k-1-2t1 i)s1 i-1-(2t1 i-1-t3 k-1)s1 i)1 i-(t3 k+t3 k-1-2t1 i)(s2 j-1+s2 j-2s1 i-1)2 j-s2 j+s1 i-1)(26)1032023 年软 件 导 刊当 t3 k-1,t3 k t1 i-1,t1 i,r3 k-1,r3 k r2 j-1,r2 j时:m1=m1 i(s)=t1 i-1+(t1 i-t1 i-1)(1-1 i)(s-s1 i-1s1 i-s1 i-1)(1-1 i)(s-s1 i-1s1 i-s1 i-1)+(
24、1-s-s1 i-1s1 i-s1 i-1)1 i(27)m2=m2 j(r)=s2 j-1+(s2 j-s2 j-1)(1-2 j)(r-r2 j-1r2 j-r2 j-1)(1-2 j)(r-r2 j-1r2 j-r2 j-1)+(1-r-r2 j-1r2 j-r2 j-1)2 j(28)从而有m3(r)=m1 m2=Kr+MNr+L,其中:K=-(1 i-1)(s2 j-1+s2 j-2s1 i-1)2 j-s2 j+2s1 i-1)t1 i-(s2 j-1+s2 j-2s1 i)2 j-s2 j+s1 i)1 it1 i (29)N=(2s1 i-2s2 j-1-2s2 j+2s1 i
25、-1)2 j+2s2 j-s1 i-s1 i-1)1 i+(s2 j-1+s2 j-2s1 i-1)2 j-s2 j+s1 i-1(30)M,L为常数,那么:Kr+MNr+L=KN+M-KLNNr+L(31)又由于:m3(r)=t3 k-1+(t3 k-t3 k-1)(1-3 k)(r-r3 k-1r3 k-r3 k-1)(1-3 k)(r-r3 k-1r3 k-r3 k-1)+(1-r-r3 k-1r3 k-r3 k-1)3 k=3 kt3 k-1+3 kt3 k-t3 k23 k-1+A(23 k-1)r-3 kt3 k-1-3 kt3 k+t3 k-1(32)这里A,-3 kt3 k-1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有理 参数 曲线 分段 根式 重新 算法 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。