基于分形理论的北京形态结构分析.docVIP

基于分形理论的北京形态结构分析 作者：肖汉，李志鹏，徐金泽 指导老师：明冬萍 （中国地质大学（北京）信息工程学院，北京 海淀 100083） 摘要 北京城市形态结构分析对北京市的长远规划具有重要意义。本文基于分形理论，运用GIS空间分析等方法，对北京城市形态结构及变化进行了定量的分析与评价，预测了北京市未来城市发展范围。首先，利用TM遥感影像监督分类，提取1992、1999和2006年北京市建筑用地面积，然后以城市中心为圆心，取几十个半径递增的同心圆进行剪裁，得出各个同心圆内的城市建筑用地面积。接着通过反复试验确定双标度区并得到城市建筑用地面积折线图，并对面积进行基于乘幂关系的函数拟合，拟合优度R2≥0.995。并得到分维数D以及北京市范围及其变化。在得到这些数据后，本文进行了深入的分析以解释这些现象，同时也对北京市与其他大城市的城市形态进行对比。最后，本文对北京市的城市建筑用地的未来发展以及可能出现的情况作出预测。联系人：肖汉，联系电话E-Mail：520043801@。Area-Perimeter Relation)法、盒子计数法(Box Counting Method)以及面积-半径关系(Area-Radius Relation)法三种。目前看来，面积-半径关系有助于从动态的角度刻画城市生长与形态[1]。本次研究就是采用这一方法，做出同心圆编号-面积折线图并拟合出乘幂函数求出分维数。这种方法在作图时也可以采用面积-半径双对数折线拟合线性函数，斜率即为分维数。这种方法给出的维数一般是从统计自相似的角度定义的，所得到的维数俗称半径维数[2]。 本次研究应用数理模型——分形理论对不同年份北京市的城市形态的发展进行评估，利用提取不同半径范围内的城市建筑面积拟合函数，首先通过函数方程判定城市分形特征，即确定城市是否存在分形。其次通过函数图像进行评价和预测，对不同年份计算出的分维指数进行分析，同时也对城市未来发展的变化及趋势做出判断。 2 基于分形理论的城市形态研究方法 2.1 城市形态分形理论 分形理论是当今世界十分风靡和活跃的新理论、新学科。分形的概念最早是由美籍数学家曼德布罗特(B.B.Mandelbort)在1967年首先提出的。自相似原则和迭代生成原则是分形理论的重要原则，它表征分形在通常的几何变换下具有不变性,即标度无关性，分形理论以分维作为的定量表征和基本参数[3]。Batty在1985年开创了分形城市形态研究[4]，他和Longley研究了城市边界和城市土地利用形态分形[5][6]。 地学中的分形方法就是用某种尺度n对空间现象或特征进行空间度量，根据尺度变化得到相应的测度A(n)计算结果（在本次研究中为城市建筑用地面积）。改变尺度n（在本次研究中为同心圆编号），测度A(n)也会随之改变。如果测度、尺度之间服从如下标度不变规律：A(kn)∝k±αA(n)，则可认为该空间现象具有分形性质，即其空间结构特征不随尺度改变而变化。从数学上，可以得出乘幂函数A(n)∝n±α满足上述泛函方程。其中，k为尺度比，α为标度指数。α通常是分维的函数，或者就是分维本身，即有α=D，这里D为分维[7]。 2.2 Smeed模型推导 一个几何体的长度L、面积S和体积V之间存在如下测度关系：L1/1∝S1/2∝V1/3∝M1/D，这里M为广义体积(可代表L、S、V=1，23)，如果某种测度对应于分形，则D可被推广到分维。对于城市来说，半径r相当于线度L，即有L∝r，但面积未必是2维。假定城市形态为D维，则由上式可知r1/1∝S(r)1/D，从而S(r)=krDk为系数，D为城市用地形态的维数。又设从城市中心向外半径为r的圆域面积为A(r)，则显然有A(r)=πr2，对上式求导，则有dA(r)/dr=2πr再对式S(r)=krD求导，可得dS(r)/dr=DkrD-1。用式dA(r)/dr=2πr两边除式dS(r)/dr=DkrD-1两边，可得ρ(r)=dS(r)/dA(r)∝rD-2。 当D2时，ρ(r)∝1/r2-D，r变大则ρ(r)变小，城市利用密度从中心向外递减； 当D=2时，ρ(r)=Dk/2π，城市利用密度从中心到外围没有差异； 当D2时，ρ(r)∝rD-2，r变大则ρ(r)变大，城市利用密度从中心向外围递增[8]。 2.3 研究方法 图1 本次研究流程图 3 数据源及预处理 3.1 数据源 本次研究采用的数据为1992年9月17日、1999年7月1日和2006年9月6日的北京市TM卫片（图2），由于成图时段均为夏季，减少了植被对于城区建筑用地的混分情况。经多次试验，本次研究RGB分别采用5、4、3.img格式文件，这样能够最大限度地减少波段之间的相关性，提取到最多的遥感信息。 a b c 图2 北京市TM影