斜轴抛物线插值的改进算法与近似算法-中国科学院.PDF

第 34 卷 第 12 期 武 汉 大 学 学 报 ·信 息 科 学 版 Vol . 34 No . 12 2009 年 12 月 Geomatics and Information Science of Wuhan U niver sit y Dec . 2009 文章编号 :167 18860 (2009) 文献标志码 :A 斜轴抛物线插值的改进算法与近似算法 李云锦1 ,2 　钟耳顺1 ,2 　黄跃峰1 ,2 ( 1 　中国科学院地理科学与资源研究所 ,北京市大屯路甲 11 号 ,100 10 1) (2 　中国科学院研究生院 ,北京市玉泉路甲 19 号 , 100049) 摘 　要 :为提高加权平均法求插值点的效率 ,用相邻斜轴抛物线顶点的切矢量构造 Hermit e 插值曲线 ,避免了 开方运算 。为提高切矢量计算效率 ,用长度比值点确定的切线方向作为切矢量近似方向 ,用弦线在切线上的 投影长度作为切矢量模 。实验表明 ,所提出的算法与斜轴抛物线加权平均算法相比,插值曲线形状相近而效 率更高 。 关键词 : 插值 ;斜轴抛物线 ; Hermit e 中图法分类号 :P208 　　常用的曲线插值方法有线性插值 、二次插值 的参数包括抛物线系数 、坐标系变换系数 。如图 ( ) 抛物线插值 、样条插值等 。文献[ 1]提出了斜轴 1 所示 , P1 、P2 、P3 为 X O Y 平面直角坐标系中非 抛物线插值 ,插值曲线光滑且与常用算法相比能 共线的 3 点, 将坐标原点 O 平移至 P2 , 令平移后 更好地表现原始曲线的局部特征 ,但是 ,该算法的 的坐标系为 X ′P2 Y′; 将 X ′P2 Y′绕 P2 逆时针旋转 计算量较大 。求斜轴抛物线需解三次方程 ,一般 α 角度, 令旋转后的坐标系为 U P V ; 设 P 、P 在 2 1 3 采用牛顿迭代法 。为减少迭代次数 ,文献[2 ,3 ]分 ( ) ( ) U P2 V 中的坐标分别为 u1 , v1 、u3 , v3 , 并用 v 别提出了改进算法 。然而 ,上述算法均采用加权 2 = k u 表示斜轴抛物线, 则有 : 平均法求插值点 ,运算量较大 。使用 Her mit e 插 2 v3 u3