基于doogles-pecker算法的化简质量研究.docxVIP

基于doogles-pecker算法的化简质量研究 在当今快速通信的信息社会中，为了快速处理和更新信息数据，有必要简化信息数据。Douglas-Peucker算法作为一种代表性的矢量线要素化简算法,在地理信息处理中发挥着重要作用。该算法可称作“数学最优”,对原线要素视觉表示最好,在分级结构、多尺度处理、图像处理、计算几何等多领域广泛应用,使其成为最经典的线要素化简方法。 运行Douglas-Peucker算法时只需要指定一个参数值(又称为阈值)。阈值的选择直接影响到化简结果和线要素的最终形态。在实际应用时,线要素化简的最佳阈值通常通过试验来确定。该阈值确定的关键在于对化简结果的评价,人们从视觉感受或经验角度定性分析及评价化简结果的质量,平衡线要素表示质量与数据量间的关系,确定最好的线要素化简结果,进而确定阈值。该评价过程具有很大的不确定性,不同的人可能对同样的化简过程得出差异较大的化简结果。化简过程可以从几何、艺术、传输、技术、表示法等多个方面加以评价,是一个灵活的、多层次的智能加工过程。化简结果是否达到要求,也需要从几何约束、拓扑约束、结构约束、Gestalt约束等多个方面加以衡量。从数学角度就已经提出了许多针对线要素复杂性的测量方法,如应用于单个线要素的“单属性测量”,及应用于线要素化简前后几何对比的“位移或比较测量”,按几何属性还可以细分为长度测量、密度测量、角度测量及弧度测量等。可以说,线要素化简的评价结果“没有最好,只有更好”,利用尽可能少的数据量表示质量尽可能高的线要素信息是判断一个化简算法好坏的最高准则。 依据具体的Douglas-Peucker算法实验结果数据,从分析化简算法中阈值与评价化简结果的几个重要几何指标的关系着手,利用曲线拟合法,建立阈值与几何指标的函数关系,通过分析函数曲线特征,确定最优的化简算法阈值。所用试验数据为经过处理的Shape格式全国主要道路线要素数据,采用WGS84坐标系,数据单位为度。Shape数据仅含有简单线要素对象(Polyline),即道路与单路径(Single Path)线要素对象间属一一对应关系。全部数据共含22 629条线性道路,道路总长计约为502 741 km,总共由366 813个基本点要素构成,每条线要素平均长度约22 km,平均由16.2个点要素构成。 1 道路几何属性信息 利用迭代方法,规定Douglas-Peucker算法有意义的阈值取值范围和增长步长,依次生成均匀布满整个阈值取值域的线要素化简结果,记录各次化简时所需统计的度量数据值。 以下是算法执行的伪函数。 /* Douglas-Peucker算法阈值迭代执行函数 ptLines-线要素集; dmin,dmax-化简阈值的最小、最大值 */ AnalyseDP(ptLines,dmin,dmax,dlnterval) { /*阈值的迭代*/ for(double dVal=dmin,dValdmax,daVl+=dlnterval) { /*遍历线要素集中所有线要素*/ for each ptLine in ptLines { 统计原线要素ptLine的几何属性信息; /*线要素化简算法/* Douglas-Peucker(ptLine,dVal); 统计化简后线要素ptLine相应的几何属性信息; } 统计线要素集ptLines化简前后的几何属性信息; 输出线要素集ptLines化简前后的几何属性信息; } } 此处所统计的与阈值取值相关的化简前后道路的几何属性信息包括:道路总长度、构成全部道路总点数、道路的平均长度、构成每条道路的平均点数等信息。 因试验用Shape数据由以度为单位的地理坐标值构成,程序运行所取的阈值dThresh范围相对较小,取值区间为[0.000 1, 0.070 0],所取步长为0.000 3°,由此共取得235组等阈值间隔统计数据。 2 do膜拟合优度函数的确定 基于最小二乘法曲线拟合是指:对给定的一组数据(xi,f(xi))(i=0,1,…,m),要求在函数类φ=span{φ1,…,φn}中指定一个函数y=s*(x),使误差平方和最小,即 ∥δ∥22=m∑i=0δ2i=m∑i=0(s*(xi)-f(xi))2=mins∈φm∑i=0(s(xi)-f(xi))2∥δ∥22=∑i=0mδ2i=∑i=0m(s?(xi)?f(xi))2=mins∈φ∑i=0m(s(xi)?f(xi))2 其中,s(x)=a0φ0(x)+a1φ1(x)+…+anφn(x), (nm)。 对于函数φ类,根据实际情况可以选择不同的函数类型。最常见的函数有:指数函数y=cebx,线性函数y=ax+b,对数函数y=clnx+b,多项式函数y=b+c0x+c1x2+c2x3+…+cixi+1,以及幂函数y=c