数字绘图中曲线拟合算法研究.docVIP

数字绘图中曲线拟合算法研究 　　摘要：讨论曲线拟合的基本原理，结合一些地形图测绘成果，分析各种曲线拟合算法，并实现加权分段圆弧法拟合曲线，分析拟合效果。 　　关键词：数字测图，曲线拟合，加权分段圆弧法 　　中图分类号：P231.5 文献标识码：A 文章编号： 　　Analysis of Curve Fitting Algorithms for Digital Mapping 　　Abstract:Combined with some topographic map surveying and mapping results to discuss the basic principle of curve fitting and analysis of a variety of curve fitting algorithm ,achieve weighted segmented arc method fitting curve and analysis fitting effect. 　　Key words: digital mapping, curve fitting, weighted segmented arc method 　　 　　1.引言 　　数字测图是以计算机及其软件为核心在外接输入、输出设备的支持下，对地形空间数据进行采集、输入、成图、绘图、输出、管理。作为数字测图过程中主要处理对象的空间数据代表着现实世界地理实体或现象在信息世界中的映射。它具有空间性，属性和时间性。空间数据根据表示对象的不同，具体可分为七种类型，即类型数据、面域数据、网络数据、样本数据、曲面数据、文本数据、符号数据。所有这些不同类型的数据都可以由点、线、面三种不同的图形表示。而曲线的绘制是在实际工作中经常碰到，也是较难处理的问题，如果曲线绘制的方法不够严格，通过数字绘图所得的图形数据往往和实际情况有较大的偏差，甚至与实际情况相悖。 　　2.曲线拟合的基本原理 　　在数字绘图的过程中，往往需要将已知离散采样点连接成光滑曲线，如通过测高点勾画等高线；通过手工数字化采集等高线数据后，重新内插生成等高线；以及通过离散的实测或统计采样点，生成等值区域等。这一过程的实质就是采样点的曲线光滑处理，也称为采样点曲线拟合。 　　曲线拟合是用连续曲线近似地比拟平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。即是用解析表达式逼近离散数据的一种方法。曲线拟合所要完成的任务也就是在一定意义下“最佳”地逼近或拟合已知数据。由于通过测量和观测所得到的数据本身就有误差，所以在绘制曲线的过程中一定要通过型值点（x，y）似乎没有必要，因此如何使得空间内选择函数Y=Ψ（x）尽可能地逼近真实函数Y=f（x），就是曲线拟合所要解决的主要问题。如图一就是使用南方CASS成图软件所绘制的，进行曲线拟合处理前后的图形。 　　 　　图1曲线拟合前后对比 　　3. 实际工作中的问题和简单处理方法 　　在实际工作中，所要绘制的线要素，都可以看成是曲线和直线的组合。如果在计算机绘图时，只考虑直线，不使用曲线的拟合功能，虽然能够保证通过测量所得到的点都能够准确的绘制在线要素上，但是绘制出来线要素并不能准确地描述出实际情况。如果在绘图工作中过于依赖曲线拟合功能，也不能准确地描述出实际情况。图2为某镇地形图的一部分，其中路中间的为隔离带。现将单个隔离带图形提取出来进行研究，图3为采用不同方法绘制出的隔离带图形。 　　 　　 　　图2 　　 　　图3 　　综上所述，在绘图工作中，应当将直线绘制和曲线拟合综合考虑。图2就是综合两种情况考虑的结果，显然绘制的效果要比简单只考虑直线绘制和过于依赖曲线拟合所绘制的效果好。 　　4．拟合算法研究简介 　　常用的曲线拟合算法有张力样条函数，分段三次多项式，分段圆弧和切线方向内插法等。 　　4．1张力样条函数法 　　 设已知特征点序列(xi,yi ),i=1,2,3,…，n，且满足条件x1x2…xn，另外还给出一个张力系数，且σ≠0。现在要求一个具有二阶导数连续的单值函数y=f(x)，使它满足： 　　yi=f(xi) i=1,2,3,…，n， 　　同时还要求f(x)-σ2f(x)必须是连续的，并且在每个区间[xi,xi+1](i=1,2,3,…，n-1)成线性变化，即： 　　 　　式中Hi=Xi+1-Xi。 　　 上式是一个二阶非齐次的常系数性微分方程，它的通解为: 　　 　　式中的Y为它对应的齐次方程f(X)-σ2f(X)=0的通解，即Y=c1e2x+c2e-2x，为它的一个特解，=aX+b。根据通过每一离散点的初始条件，就可解算出c1,c2和a,b的值。最后经过整理可得到单