城市交通网络分形维数的不确定性估计_控制与分析_李畅

  • 格式:pdf
  • 大小:1.63 MB
  • 文档页数:10

城市交通网络分形维数的不确定性估计、控制与分析李畅,李桂娥,朱昱佳华中师范大学城市与环境科学学院 地理过程分析与模拟湖北省重点实验室,武汉 430079摘 要:长度-半径维数模型作为描述城市交通网络复杂不确定性现象的一种分形分维方法,其自身存在的不确定性往往被忽视,且相关研究更是鲜见报道。

故针对该模型在分形维数测算全过程中存在的不确定性问题,本文率先开展了系统剖析、定量估计和质量控制研究。

首先对数据源、矢量化处理、测算中心、尺度选择、以及分维数模型估计等一系列环节进行了不确定性估计与分析,其中首次给出了分形维数在一定置信水平下的不确定性度量区间,并依据误差传播理论对误差的传递和累积进行了描述;然后着重提出了基于LMedS(Least Median of Squares)的质量控制方法。

最后通过对拉萨市的算例实验表明:道路的矢量化过程、测算中心和测算尺度的选择都会导致分维的不确定性;并在对数据质量进行控制的基础上,通过置信区间对长度-半径维数模型的不确定性进行了在一定概率水平下的首次度量;同时结合区域现状对研究结果给出了合乎实际的解释。

本文在描述表征不确定性问题的分形几何和分形维数的基础上,系统地揭示了其自身不确定性的本质,不仅进一步丰富了分形分维理论,为控制其质量奠定理论基础,而且可为城市交通网络分形维数的地学应用提供可靠的科学依据。

关键词:交通网络,长度-半径维数模型,不确定性,LMedS ,质量控制,误差传播中图分类号:N93 文献标志码:A引用格式:李畅, 李桂娥, 朱昱佳. 2017. 城市交通网络分形维数的不确定性估计、控制与分析. 遥感学报, 21(1): 74–83Li C, Li G E and Zhu Y J. 2017. Uncertainty estimation, control, and analysis of fractal dimension for urban traffic network. Journal of Remote Sensing, 21(1): 74–83 [DOI:10.11834/jrs.20176023]1 引 言¢¢¢¢¢¢正是因为客观世界难以找到纯数学的理想分形,但统计意义下的随机分形却是客观存在的,因而分形理论成了描述复杂、非线性空间分布的有力工具(Mandelbrot 和Aizenman ,1979)。

分形理论的精妙在于把一个不确定性的分形几何问题转换成为了一个确定性的维数问题,如Cantor 集、Koch 曲线等,一旦尺度确定,就可以得到确定性的分形维数,例如Cantor 集尺度与对应份数的关系(尺度,份数)为(1/3, 2)、(1/9, 4)、(1/27, 8)、(1/3n, 2n),其对应的Hausdorff 维数为log(2)/log(3)=log(4)/log(9) =log(8)/log(27)==log(2n)/log(3n)=0.6309,取对数后的点对(log3i,log2i)严格¢¢¢¢¢¢成比例,即共线;然而,自然当中存在的统计意义下的分形几何问题,在不同尺度下得到的对应的分维数是不确定的,需做回归分析才可得到一个相对确定的维数,例如经典的海岸线长度问题,以中国海岸线为例,尺度与长度的点对关系(尺度,长度)为(马建华 等,2015):(400,4156)、(200,4776)、(100,5486)、(50,6304)、,然而对其取对数后的点对(log400,log4156)(log50,log6304)并不共线,故需进行回归分析。

剖析上述问题的本质原因在于研究对象自相似性的程度,自然界中事物的非规则本质使得无规分形无法获得绝对意义上的事物的确定性分形特征,并在测算过程中不可避免产生不确定性,且必须进行回归分析,这正是无规分形具有不确定1007-4619(2017)01-0074-10Journal of Remote Sensing 遥感学报收稿日期:2016-01-28;修订日期:2016-06-17;优先数字出版日期:2016-06-24基金项目:国家自然科学基金(编号:41101407,61301278);湖北省自然科学基金(编号:2014CFB377,2010CDZ005);华中师范大学中央高校基本科研项目(编号:CCNU15A02001);武汉市晨光计划(编号:2016070204010137);华中师范大学中央高校基本科研业务费项目(编号:CCNU16JCZX09)第一作者简介:李畅(1982— ),男,副教授,研究方向为地理信息科学、数字摄影测量与遥感的理论、技术和地学应用研究。

E-mail :lcshaka@ ;lichang@性的根源所在,同时这也是激发本研究——分形维数不确定性的理论关键动机。

综上所述,分形理论自身同样存在不确定性,而且已经引起少数学者的思考,当前学者对于分形不确定性的研究已涉及对盒子维数的测量算法(梁东方等,2002)、河流(李静静等,2010)和海岸线(马建华等,2015)长度的测算等方面,且有学者对分形不确定性进行了计算机模拟(李德毅等,2004),相关研究多着眼于测量尺度选择对于分形分维的影响,对本研究也无不启迪,在此基础上我们对其开展了系统性的不确定性研究。

鉴于分形分维理论在地学应用研究中对交通网络分形特征的分析较为经典(B e n g u i g u i和D a o u d,1991;B e n g u i g u i,1992;柏春广和蔡先华,2008;Li 等,2012),因而本文拟以此为研究对象分析常用于描述其分形特征的长度-半径维数的不确定性。

回顾以往学者的研究,对于城市交通网络长度-半径维数的不确定性分析仅有少量文献(刘妙龙和黄佩蓓,2003;刘承良等,2013)在研究过程中指出测算中心和尺度选取不同,会对解算交通网络的真实分维值产生影响,但并未对其不确定性进行系统研究和剖析。

基于此,本文尝试对城市交通网络分形分维开展不确定性来源分析、度量及其质量控制的研究,以期丰富分形分维的相关理论与质量控制方法,并为全面、客观分析研究对象的分形分维特征提供可靠依据。

2 分形维数的不确定性分析原理与质量控制本文拟以经典的长度-半径维数模型为例对其计算过程中的不确定性进行分析与度量。

因此,首先回顾一下该模型(Frankhauser,1990; 陈彦光,2008)。

一般地,对于一个几何体,设长度为L,面积为S,体积为V,则有借助广义体积M,则(1)的广义化形式可以表示成:式中,D表示欧氏维数,倘若几何体的某个测度具有分形特征,则其相应的D值为非整数。

如果一个面积为S的区域,其中的交通网络具有分形特征,则其总长L(S)与面积之间应该具有如下关系当所选区域为圆形时,S∝r2,则式(3)可表示为以下形式式中,r为以某枢纽为中心的回旋(转)半径,L(r)为半径为r的区域范围内的道路网络总长度,L1为常系数,幂指数D即分维,Frankhouser将其命名为半径维数(radial dimension),也可称之为长度-半径维数。

从中不难发现,在对城市道路网络进行长度-半径维数测算时,道路网络的数据选择与预处理、测算中心和测算尺度的选择、分形维数模型本身都将成为该模型的不确定性来源对真实分维值的测算产生影响。

所以,本文将从(1)数据源、(2)矢量化处理、(3)测算中心及尺度选择、(4)分维数模型估计等环节对其进行不确定性分析,并基于误差传播的累积效应,在一定概率水平下通过置信区间对长度-半径维数模型的不确定性进行首次度量,同时,基于LMedS(Least Median of Squares)提出相应的质量控制方法。

2.1 分形维数的不确定性分析原理(1)数据源的不确定性。

对于研究学者而言,城市交通道路网的基础数据源主要来自于对纸质或电子地图的矢量化并同时辅助其他工具对数据进行完善。

从数据获取源来看,影响数据质量的因素主要包括地图测量误差、绘制误差,图纸变形误差(柳宗伟等,2002),地图要素本身的密度和复杂程度以及地图综合(王桥等,1998; 艾廷华等,2007; Li,2007)过程中对于地物的取舍程度等。

(2)数据转换的不确定性。

从城市交通道路矢量数据处理方式来看,数字化以及数字化操作人员的熟练程度及操作方式、数据格式转换等同样存在不确定性。

对于同一区域,不同比例尺的地图对分维数会产生影响(陈杰和马素媛,2012),即使是同一比例尺的地图,在不同比例的矢量化环境下,其结果也不相同,从而导致不确定性。

故,矢量化的结果,例如:点、线、面皆存在不确定性,其中点的不确定性可以由误差椭圆来描述;线的不确定性可以由G-带误差模型进行概李畅等:城市交通网络分形维数的不确定性估计、控制与分析75括;面的不确定性则可通过讨论多边形组成顶点的不确定性和基于多边形边界的误差带两种方法对其不确定性进行度量(Shi,2009)。

(3)测算中心和测算尺度的不确定性。

对交通道路网的长度-半径维数进行测算时需要选择测算中心和测算尺度。

从当前学者的研究工作来看,测算中心或为城市交通网络的枢纽或为城市的几何中心亦或是重心,不同的选择所获得结果也不尽相同,故在中心选择上存在认知的不确定性。

即使同是选择城市道路网的交通枢纽中心,也会因研究者不同,而有差别,即相应测算中心都存在一个点位误差(测算点与真实点),可用点位误差椭圆(Shi,2009)来衡量。

对于测算尺度的选取,就目前而言,没有固定的标准,多数学者都是结合研究区域的特征选取合适的尺度进行测算,这也为测算区域交通网络的分维值增加了多样性。

(4)分形维数估计的不确定性度量。

对于分形维数模型,有些不确定性可以通过对数据重新定位、浏览、查错等进行降低,而有些则是无法避免且随处理过程进行误差传播累积。

对于此类不确定性可通过计算分维数的置信区间进行定量估计。

对式(4)作对数变换,则可得到:为了表示方便,在以下公式中,凡涉及点对(ln r,ln L(r))均用(x i,y i)代换,ln L1用a代换,分维数(即回归系数)用b表示,则文中长度-半径对应的点列的回归方程可用下式表示:(a)的最大似然b(回归系数)的最大似然估计值,(x i,y i)表示上述代换后对应的样本点对,εi表示随机误差。

在以下公式中凡涉及到样本的个数均用n表示,为了计算上的方便,特引入下述记号这样,a,b的最大似然估计值就可以表述为为了更有效地描述随机变量之间的相关关系,特引入残差(残差为某样本点观测值与拟合值之差)平方和,记为Qe程值,则¢¢¢将样本中的y i改为Y i(i=1,2,,n)(Y i是指相互独立的随机变量),则若残差平方和Q e服从分布,且令σ2表示均方误差:则这样就可以得到σ2的无偏估计量将道路长度和半径长度标绘在双对数坐标图上,判别其回归效果,若回归效果显著,则可以对分维数作区间估计。