1、收稿日期:20220215第 40 卷第 2 期计算机仿真2023 年 2 月文章编号:10069348(2023)02027504城市绿地率空间分异格局及演变规律仿真谢兴宇，高向东(东华理工大学，江西南昌 330000)摘要:为了更好地维持城市绿地系统的稳定，获取城市绿地率演变规律特征，提出城市绿地率空间分异格局及其演变规律模拟方法。预处理遥感数据，基于数据处理结果提取城市绿地率空间分异格局特征;通过城市绿地斑块之间的连通关系构建绿地连通图，采用邻近事件矩阵统计绿地连通空间的排列特征，获取城市绿地率演变规律。将得到的全部特征输入到Matlab 仿真软件中展开仿真分析，经实验测试结果证明，所提
2、方法可以获取满意的模拟结果，演变规律模拟效率与精度较高。关键词:城市绿地率;空间分异格局;演变规律;连通图;邻近事件矩阵中图分类号:TP393文献标识码:BSimulation of Spatial Differentiation Patternand Evolution Law of Urban Green Space ateXIE Xingyu，GAO Xiangdong(East China University of Technology，Jiangxi Nanchang 330000，China)ABSTACT:In order to maintain the stability o
3、f urban greening system and obtain the characteristics of evolutionlaw，this article puts forward a method of simulating the spatial heterogeneity pattern of urban greening rate and its e-volution law At first，the remote sensing data was preprocessed On this basis，the characteristics of spatial heter
4、oge-neity pattern were extracted And then，a connectivity map was constructed through the connectivity relationship be-tween urban green space patches Moreover，the adjacent event matrix was used to obtain the arrangement character-istics of connected space Finally，the evolution law of urban greening
5、rate was obtained All the features were inputinto the MATLAB software for simulation analysis Experimental test results show that the proposed method can obtainsatisfactory simulation results;And the simulation efficiency and accuracy of the evolution law are highKEYWODS:Urban green space rate;Spati
6、al differentiation pattern;Evolution law;Connectivity graph;Adjacent e-vents matrix1引言城市绿地是城市绿化系统和城市景观的重要组成部分，具有十分重要的生态以及经济功能。明确城市绿化用地的概念以及发展过程，了解城市绿化用地对生态环境的作用，获取城市绿地率空间分异格局及其演变规律，可以更好理解城市绿地之间的相互作用和相互联系，为城市绿地管理规划提供更好的决策支持，对城市用地的合理规划以及建设具有十分重要的意义1，2。国内相关专家针对绿地率空间分异格局及其演变规律模拟方面的内容展开了大量研究，例如唐燕等人3 通过超效
7、率 DEA 模型分析产业在设定时间段内的生态水平，同时以地理区域划分依据，获取生态效率演变特征和空间分异格局特征。陈学兄等人4 通过 Fragstats 4.2 计算城市斑块类型水平景观指数，对不同地形起伏程度的分异特征展开分析，最终获取格局分异特征。陈浮等人5 主要通过空间自相关分析方法和多元回归分析方法完成空间格局分异特征以及驱动特征提取。在以上几种方法的基础上，提出一种城市绿地率空间分异格局及其演变规律模拟方法。经大量实验测试证明，所提方法可以以更快的速度和更高的精度完成城市绿地率空间分异格局及其演变规律模拟。5722城市绿地率空间分异格局及其演变规律模拟2.1城市绿地率空间分异格局由于
8、遥感系统空间以及时间等多方面因素的限制，遥感图像中的信息无法得到有效保存，导致获取的遥感数据存在比较大的误差。这不仅会降低数据的质量，同时也会影响城市绿地率空间分异格局分析精度。因此，在获取城市绿地率空间分异格局特征之前，为了有效解决该问题，需要预处理原始城市绿地遥感图像。现阶段，遥感信息的分类6，7 以及分析通常要结合目视解释，在目视解释中，波段选择需要注重以下几方面的因素:1)选择的波段信息量要比较大;2)各个波段之间的相关性要比较小;3)各个地物之间的光谱差异要比较大。通过最佳指数 OIF 相关概念可以获取各个波段之间的相关系数，进而得到最佳指数因子OIF=ni=1Q(i)ni=1|si
9、j|(1)式中，Q(i)代表第 i 个波段的标准差;sij代表波段 i 和波段 j 两者之间的相关系数。选择标准差作为测试指标可以有效描述城市绿地遥感图像灰度均值的离散程度，其中，标准差的取值越大，则说明城市绿地率遥感图像中灰度级分布越分散，这样可以获取更加详细的图像细节信息。选择标准差作为评价指标可以从信息量以及影像空间分辨率两个角度完成遥感图像评价，标准差 std 的具体计算式如下所示std=mi=1nj=1 G(i，j)sijm n(2)式中，G(i，j)代表灰度值;mn 代表城市绿地遥感图像的大小。通过遥感图像的相关系数可以准确反映随机两幅城市绿地率遥感图像之间的相关程度，经过计算得到
10、图像融合前后的相关系数，进而获取图像的变化程度。相关系数的取值越大，则说明图像的保真效果越好。其中，相关系数对应的计算式为D(f，g)=ij(hi，j uf)(gi，j ug)ij(hi，j uf)2(gi，j ug)(3)式中，D(f，g)代表相关系数;hi，j和 gi，j代表两幅图像的灰度值;uf和 ug代表两幅图像像素平均值。为了提升城市绿地率空间分异格局特征提取精度，进一步对遥感图像进行几何校正8，9。遥感图像校正可以划分为多种不同的形式，由于城市绿地内部区域大部分为平摊区域，不需要考虑其它地形因素，所以使用多项式模型校正，因为多项式模型稳定性比较好，操作方式也十分简单，在实际应用过程
11、中利用率也比较高，以下给出详细的操作步骤:1)优先确定各个地面的控制点，获取图像和地图对应控制点的像元坐标，并根据对比分析选取最佳坐标点。经过变换得到图像坐标(x，y)和参考坐标(X，Y)之间的关系式，如式(4)所示x=mi=0nj=0aij(X，Y)y=mi=0nj=0bij(X，Y)(4)上式中，aij和 bij代表多项式系数。对于发生几何形变的遥感图像而言，为了确保后续图像处理结果的准确性，需要借助一次线性方程完成图像校正处理。而针对变形比较严重的图像而言，则需要借助二次或者三次多项式完成校正处理，具体如式(5)x=a0+a1X+a2Y+a3X2+a4XY+a5Y2y=b0+b1X+b2
12、Y+b3X2+b4XY+b5Y2(5)当多项式的项数确定之后，控制点坐标的选取是通过最小二乘回归方式对多项式求解，以下采用式(6)计算均方根误差MSerror=(x x2)+(y y2)(6)式中，MSerror代表均方根误差。在上述步骤的基础上，通过面向技术的分类对象提取方法提取城市绿地率空间分异格局特征10，11，详细的操作步骤如图 1 所示。图 1城市绿地率空间分异格局特征提取流程图2.2城市绿地率空间演变规律模拟城市绿地连通性的相关研究主要依赖于连通图，所以将672连通图作为获取城市绿地率演变规律的基础，同时采用邻近事件矩阵统计绿地连通空间的排列特征，最终获取演变规律。当随机两个绿地之
13、间的最短边缘距离大于选定的距离阈值，则说明两个绿地斑块之间属于连通关系。当绿地之间的最短距离小于或者等于选定的距离，则说明两者是连通的。连通图使用无向完全拓扑图展开描述，通过空间位置描述城市绿地斑块的几何中心点。其中，连通图的建立过程如下所示:1)通过全部绿地的几何中心点构建节点集合，形成连通图，同时使用节点代表绿地斑块以及空间位置。2)根据选定的距离阈值对绿地实行缓冲区分析处理，同时将其表示为多边形。将全部元素和组建线条之间的元素进行对比分析，将符合需求的元素写入边集合中。重复上述操作，直至绿地斑块对应的边全部记录在边集合中。3)依次对各个绿地斑块遍历处理，同时重复步骤(2)，直至全部绿地完
14、成遍历为止。4)通过全部的线以及节点共同构建连通图。其中，连通图中距离阈值的大小需要满足实际研究目的，同时还需要通过城市绿地的特定决定。通过连通图可以准确反映城市绿地连通的分布情况，而绿地连通性的空间属性可以采用空间排列准确反映出来。为了方便对城市绿地之间连通关系的量化处理，在连通图研究的基础上，统计不同类型绿地空间连通排列频率对应的邻近事件矩阵。在图论中，通过连通性可以准确反映连通图结构的连通程度，连通程度越好，则说明图结构就越稳定。以下主要通过贝塔指数以及介数中心性指数完成绿地连通性分析:1)贝塔指数:通过贝塔指数完成连通图非空间关系度量处理。同时，连通性水平也可以利用连通图中边的数量 e
15、 和节点数量 c 两者之间的比值描述，其中贝塔指数对应的计算式为=ec(7)式中，代表贝塔指数。通过贝塔指数可以准确反映城市绿地之间的空间关系和复杂程度。在距离阈值取值相同的情况下，连通图的贝塔指数越高，则说明城市绿地之间的关联性和连通性越强。2)介数中心性指数:通过介数中心指数对连通图中比较重要的绿色斑块实行特征提取，通过介数中心性指数可以准确描述随机一个节点的重要性，取值越大，则说明节点越重要。其中，节点 k 对应的介数中心性指数可以表示为式(8)的形式HB(k)=s，t=1st(k)st(8)式中，HB(k)代表介数中心性指数;st(k)代表最短路径中节点通过的数量。在连通图的基础上获取
16、不同类别城市绿地的邻近事件矩阵，城市绿地率演变规律主要是指城市绿地在时间上的空间动态变化过程。通过图论理论获取城市绿地的动态演变模型，同时经过一系列的分析操作，组建城市绿地率演变图。城市绿地率演变规律演变模型需要满足以下的约束条件:1)连续:当随机两个绿地在两个相邻年份或者相似空间范围内，则说明两块绿地面积的重叠率大于或者等于 95%。2)扩张:随机选择两块绿地作为研究对象，当其中一块城市绿地的空间范围明显小于另外一块空地，则说明两个绿地之间的斑块不存在任何关联或者重叠的面积比较小。根据上述约束条件构建城市绿地演变图，利用图 2 给出城市绿地率演变规律示意图。图 2城市绿地率演变规律示意图城市绿地率演变规律的获取流程如下所示12:1)判断设定年份城市绿地连通图中各个边的变化情况，从设定年份的编号开始，通过城市绿地演变规律判断绿地所属于的演化类型，同时获取在设定年份和绿地对应的绿色斑块13，14。假设斑块处于不连续状态，则说明边相比前一年发生了比较明显的变化，最终会消失。假设绿地斑块是连续的，则说明设定年份的连续边没有发生任何改变，同时也不存在边消失的情况15。2)重复步骤(1)，判断连