曲面样条函数在地图数字化中应用浅析.docVIP

曲面样条函数在地图数字化中的应用浅析 张佰慧1 杨明强2 1、平顶山工业职业技术学院 2、平顶山工学院 河南 平顶山 467001 摘要：本文通过对曲面样条函数在进行曲面拟合的分析，总结了实际应用中参数的选取对曲面拟合结果的影响程度，提出了实际应用的步骤。 关键词：样条函数 地图数字化 应用浅析 1.引言 在ARC/INFO中，通常采用TIC点来建立地图坐标与地理坐标之间的联系，特别是发生旋转变形后的地图坐标，如图－1所示。 图－1 用TIC点实现地图坐标与地理坐标的转换示意图 地图坐标与地理坐标之间的转换通常有两种实现方式：一是投影公式正算与反算，二是通过实际控制点集（TIC）。 通过投影公式正算与反算，其应用有很大的局限性，因为它需要知道地理坐标与地图坐标的解析变换关系，而这种解析关系是多样的也是复杂的，在此不做讨论。而通过实际控制点集是一种简易快速的坐标变换方法，它不需要源坐标和目的坐标的解析关系，只要知道有限的地物特征点（实际控制点）即可。由于曲面样条函数在2.5维空间表面拟合有优良的表现，通常把它用来实现地理坐标与地图坐标转换。 地图精度校正，以及不同比例尺、不同投影系的地图综合处理等等，这些涉及到的坐标变换问题也可以通过曲面样条函数来解决。因此，曲面样条函数在“坐标换算”中有着重要的应用。 曲面样条函数具有诸多优点，就用于曲面样条函数拟合的原始数据点而言，它们不必按规则排列，只要提供不少于三个的任意点（控制点或采样点），就可以得到任意阶可微的光滑曲面。显然，这将有利于得到光滑的曲面，不至于因噪声（数字化误差等）的不均匀分布产生的局部畸变，同时可以通过参数控制实现灵活的方式进行曲面拟合，以满足实际需要。2.应用原理 曲面样条函数可看成是无限大平板在多个受力点的作用下发生的弯曲变形。作为平板来讲，有其物理特性，如挠度等。在不同的物理特性情况下，即便是相同的受力情况，也会产生不同的曲面变形。因此，通过不同拟合参数模拟其物理特性可灵活地实现满足不同要求的曲面拟合。例如在数字化地图校正中，一般要求以平差的方式对图形进行校正，此时可以假定曲面为一刚性平板。 曲面样条函数的表达式为： （1） 式中，――待定系数 ――调节曲面曲率大小的经验参数，它视实际情况适当选取。一般地，需要曲面减少畸变、曲面平滑，应使，对于畸变较大的曲面甚至取。 待定系数可以通过下列方程组求得： （2） 式中，，是关于点的弹性系数，＝，＝0。 （2）式的矩阵形式为： （3） 其中 这是一个对称方程组，为保证解的稳定性，可用对称方程组的HOUSEHOLDER变换法求解。 如上所述，曲面样条函数在进行曲面拟合时，有两个可以调控参数和。 为拟合曲面曲率调节的经验参数，在地图校正的实际应用中，当越大时，插值点的取值受最近的控制点（已知点）的影响越大，这样可以表示一些“细节信息”，同时由于不能过多考虑周围的控制点（已知点），直接的后果是容易导致畸变，如不相交的图形发生相交等；当越小时，由于能够较好地考虑周围控制点（已知点）对插值的影响，曲面的变化率较低，比较平缓光滑，但同时消弱了“高频信息”。以下是规则的格网在不同参数下的坐标变换的比较，如图－2。 (a) 变换之前格网 (b) 变换之后格网 图－2 不同参数下的曲面样条函数的坐标变换示意图 是一个关于弹性系数的参数，根据的不同取值，可以将曲面拟合分为三种应用方式：a.柔性拟合（如生成DEM、建立地理坐标与地图坐标联系等）；b.弹性拟合（如光滑表面）；b.刚性拟合（如数字化地图的校正），三种拟合方式的曲面的一些性质如表－1： 表－1 不同拟合方式的比较表 拟合方式 柔性拟合 弹性拟合 刚性拟合 曲面形状 弯曲 弯曲 近似平面 通过控制点 通过 不通过 不通过 在实际应用中，当＝0时称为柔性拟合曲面，其特点是曲面通过控制点（已知点），为弯曲变形性最大的曲面，如图－3（a）；当＝时称为刚性拟合曲面，其特点是曲面为一（近似）平面，曲面不通过控制点（已知点），实际控制点的值到该曲面的距离平方和一般为最小，类似如最小二乘法的结果，如图－3（b）；介于两者之间的称为弹性拟合曲面，其特点是不通过控制点（已知点），曲面类似如在力的作用下发生弹性形变。 (a) 柔性拟和曲面变换 （b）刚性拟和曲面变换 图－3 不同弹性系数的变换示意图 3.应用步骤 以下是以数字化地图校正为例，应用曲面样条函数的几个主要步骤： ①