基于非匀质空间下Voronoi图的城市影响范围研究.pdf
《基于非匀质空间下Voronoi图的城市影响范围研究.pdf》由会员分享,可在线阅读,更多相关《基于非匀质空间下Voronoi图的城市影响范围研究.pdf(6页珍藏版)》请在咨信网上搜索。
1、2097-3012(2023)03-0410-06 Journal of Spatio-temporal Information 时空信息学报 收稿日期:2022-04-01;修订日期:2023-06-01 作者简介:张晓贺,研究方向为地理信息系统与遥感图像处理。E-mail:nwu_ 基于非匀质空间下 Voronoi 图的城市影响范围研究 张晓贺 山西省测绘地理信息院,太原 030001 摘 要:为解决城市影响范围研究中忽略非匀质空间的问题,本文以光栅扫描算法为基础,提出了一种非匀质空间下面状加权 Voronoi 图生成算法。首先,对面状目标和栅格空间赋予不同的权重;其次,采用 33 邻域模
2、板进行十字交叉 4 次光栅扫描,在距离变换中按栅格像元传递权重对距离进行分解,完成非匀质空间下面状加权Voronoi 图的生成;最后,以山西省 11 个地级市为研究区,进行了验证。结果表明,改进后的算法实现了非匀质空间面状加权 Voronoi 图的生成;相比点状加权 V 图,城市之间权重差异较大,距离较近时,面状加权 V 图更适合城市影响范围划分。关键词:加权;Voronoi 图;非匀质空间;城市影响范围 引用格式:张晓贺.2023.基于非匀质空间下 Voronoi 图的城市影响范围研究.时空信息学报,30(3):410-415 Zhang X H.2023.Research on urban
3、 influence range based on Voronoi diagram in inhomogeneous space.Journal of Spatio-temporal Information,30(3):410-415,doi:10.20117/j.jsti.202303013 1 引 言 城市影响范围是城市在社会、经济、政治、文化等因素共同作用下对周围区域所及的最大地域范围(顾朝林等,1999)。确定城市影响范围,是城市规划的一项重要工作,有利于明确城市间的空间作用关系,促进城市合理布局。常用的城市影响范围研究方法有赖利模型、引力模型、潜力模型、断裂点模型、缓冲区分析、Vor
4、onoi 图(简称 V 图)、加权 V 图等。其中,赖利模型、引力模型、潜力模型、断裂点模型侧重理论分析,不能实现城市影响范围的具体划分(邓羽等,2013;常小刚和闫浩文,2015);缓冲区分析不能对空间进行完全划分(黄杏元和徐寿成,1998);V 图是对空间的一种完全划分,可以实现城市影响范围的具体划分(赵春燕等,2010;赵佳星等,2020)。V 图生成算法可分为矢量和栅格两类。常用的矢量算法有分治算法、插入算法、扫描线算法等。矢量算法计算复杂,不利于海量数据的处理(李成名和陈军,1998),且只能处理点和线段(Okabe 等,2000)。常见的栅格生成算法有扩张算法、距离变 换算法、层次
5、算法、细分算法、确定归属算法、光栅扫描算法(胡鹏等,2006)等。其中,综合考虑精度和效率因素,光栅扫描算法更具有优势,不需要考虑空间目标数量、分布和形状(Fabbri 等,2008)。传统光栅扫描算法,通过 33(或者 55)邻域模板在正反两次扫描中将一个栅格的信息传递给邻近的栅格,生成普通 V 图(Chen,1999)。目前虽然有一些改进算法(Xu 等,2014;江威等,2015;刘青平等,2019),但还不能满足非匀质空间和面状加权。李佳田等(2016)提出了非理想平面梯度 V 图构建算法,只适合简单的点发生元;田松等(2014)基于结晶法提出了顾及障碍物的加权V 图算法,谢顺平等(20
6、10)提出了一种基于网络分析的加权 V 面域图构建算法,可以将交通路网纳入 V 图的生成中,但都未扩展到整个非匀质空间。基于此,本文改进了光栅扫描算法,研究了非匀质空间下面状加权 V 图的生成。以山西省 11 个地级市为研究区,选取城市常住人口、辖区面积、国内生产总值(GDP)等 18 项指标,采用主成分分析法构建了城市综合实力评价指标体系,用山西省坡度图构建非匀质空间,生成了其地级市影响范围 张晓贺:基于非匀质空间下 Voronoi 图的城市影响范围研究 411 图,并进行了对比分析。以期为促进山西省城镇科学规划提供参考。2 非匀质空间面状加权 V 图构建方法 2.1 方法原理 本文针对面状
7、加权目标,改进了传统光栅扫描算法距离变换中的传递信息;针对非匀质空间,由于每个栅格像元传递权重不同,单个像元与目标点之间有众多不同传递权重的其他像元,其到目标点的最短路径会沿着传递权重最高的像元延伸,这样就导致最短路径不一定是直线,不能直接通过横纵坐标值来计算欧氏距离,所以本文按栅格像元传递权重对距离进行分解,用距离累加值代替传统的欧氏距离;为了提高归属精度,采用了“十字交叉”光栅扫描。在光栅扫描的时候使用 33 邻域模板(图 1),通过“十字交叉”4 次扫描实现距离变换,如图 2所示:“左下右上”纵向扫描时使用了模板(a);“右上左下”纵向扫描时使用了模板(b);“左上右下”横向扫描时使用了
8、模板(c);“右下左上”横向扫描时使用了模板(d)。图 1 八邻域示意图 Fig.1 8-neighborhood sketch 图 2 扫描模板 Fig.2 Scanning template 在距离变换的时候,所传递的信息为距离累加值最小时对应父结点行列号与最近目标点的行列号、目标编号、目标权重及距离累加值:12/()kknDDdww,minknnDD(1)式中,knD 为第 n 次扫描 k 邻域计算的距离累加值;kD 为 k 邻域传递过来的距离累加值;d 为距离递增值,如图 3 所示,四角取值2,上下左右取值1;w1为 k 邻域传递过来的目标权重;w2为该栅格像元在非匀质空间中的传递权重
9、。Dn为该栅格像元第 n次扫描计算的最小距离累加值。图 3 距离递增图 Fig.3 Incremental distance graph 四次交叉扫描结束后,连通最短路径生成非匀质空间下面状加权V图。由于采用了光栅扫描方式,改进后的算法时间复杂度同样只受栅格化后像元数量的影响,不受目标数量、分布和形状的影响。2.2 算法描述 输入:(1)发生元,即面状目标矢量数据,按城市综合实力评价指标赋予权重;(2)非匀质空间,即栅格传递权重图,本文按坡度对栅格赋予传递权重。输出:非匀质空间下面状加权V图。算法:(1)将面状目标和非匀质空间组合栅格化为MN大小的3通道图像,通道1记录最近城市编号,通道2记录
10、最近城市权重,通道3记录栅格所代表空间的传递权重;(2)定 义 用 于 标 注 的 结 构 体 数 组mark TMNT.r,T.c,T.x,T.y,T.n,T.w,T.d,分别记录父结点行列号(r,c)、最近目标点的行列号(x,y)、城市编号 n、城市权重 w,以及距离累加值 d;(3)“左下右上”纵向扫描,按式(1)计算4邻域传递过来的最小距离累加值,并用Tij记录最小值对应父结点行列号与最近目标点的行列号、城市编号、城市权重,以及距离累加值;(4)同步骤(3)继续完成“十字交叉”4次扫描并连通最短路径;(5)用TMN.n 生成面状加权V图。3 城市综合实力评价指标体系 3.1 选择评价指
11、标 城市综合实力评价是一个多维度的概念,因此412 Journal of Spatio-temporal Information 时空信息学报 2023,30(3)要遵循全面性、可比性、可操作性等原则,选取不同角度的指标进行评价。如表1所示,本文从城市规模水平、城市经济水平、社会发展水平、基础设施水平四个方面选取了18项指标来构建山西省11个地级市评价指标体系,评价指标数据来自2020年山西省统计年鉴(山西省统计局,2020)。表 1 城市综合实力评价指标体系 Tab.1 Evaluation index system for urban comprehensive strength 一级指标
12、 二级指标 城市规模水平 X1 城市常住人口/万人 X2 辖区面积/km2 城市经济水平 X3 GDP/万元 X4 人均 GDP/元 X5 一般公共预算收入/万元 X6 一般公共预算支出/万元 X7 第二产业总产值/亿元 X8 第三产业总产值/亿元 X9 社会消费品零售总额/万元 X10 固定资产投资(不含农户)/万元 社会发展水平 X11 住户存款余额/万元 X12 全体居民人均可支配收入/元 X13 教育支出/万元 X14 研究与试验发展经费/亿元 基础设施水平 X15 卫生机构床位数/张 X16 普通高等学校数量/所 X17 互联网宽带接入用户数/万户 X18 公路通车里程/km 3.2
13、 数据处理(1)标准化处理。对原始数据进行标准分数处理,首先计算数据的平均值和标准差,然后计算数据的标准化值。(2)主成分分析。先通过正交变换将一组可能 存在相关性的变量数据转换为一组线性不相关的变量,求得标准化后数据的协方差矩阵,计算特征值和相对应的特征向量;然后以主成分对总方差累计贡献率大于85%的标准提取 n 个主成分。由表2知,前两个因子的累计贡献率已达92.061%。表 2 总方差解释表 Tab.2 Explanation of total variance 初始特征值 提取载荷平方和 成分 总计 方差百分比 累计/%总计 方差百分比 累计/%1 12.314 68.409 68.4
14、09 12.314 68.409 68.409 2 4.257 23.652 92.061 4.257 23.652 92.061 3 0.771 4.286 96.347 (3)求得11个地级市的综合实力分数:211nmijjikkjSAPW (2)式中,Ak为第k个主成分贡献率;Pij为i市第j项指标的z-score标准化值,Wj为指标j在主成分k上的得分系数。为了消除负值,便于后期V图研究,张晓贺:基于非匀质空间下 Voronoi 图的城市影响范围研究 413 利用式(2)对城市得分进行线性映射,映射区间A,B取值为1,10,计算结果作为城市权重,如表3所示。其中:min1maxmin(
15、)iiSSwBAASS (3)式中,A为映射区间最小值;B为映射区间最大值。表 3 城市综合实力最终得分 Tab.3 Final score for urban comprehensive strength 城市 综合实力得分 映射后得分 太原 1.96 10 大同 0.07 4 阳泉 1.08 1 长治 0.02 4.27 晋城 0.49 2.76 朔州 0.83 1.75 晋中 0.03 4.28 运城 0.35 5.24 忻州 0.25 3.45 临汾 0.21 4.81 吕梁 0.15 4.64 4 山西省城市影响范围 以2021年山西省11个地级市建城区范围作为V图的发生元,并将表3
16、中映射后得分作为影响权重赋予相应的发生元;因为城市影响主要通过人类活动实现,而人类活动会受地形影响,坡度越大,人类活动越受限制,即栅格像元的传递能力和坡度成反比,所以取山西省坡度图作为构建非匀质空间要素,按照式(4)给每个栅格像元赋予传递权重,坡度线性映射区间A,B取值为1,50:min2maxmin1()iiGGwBAAGG (4)式中,w2i为栅格像元i的传递权重;Gi为栅格像元对应的坡度值();A为映射区间最小值;B为映射区间最大值。其中,山西省坡度区间为0 45.264,如图4所示。以改进后的面状加权V图算法为基础,利用ArcEngine开发相应模块,生成山西省11个地级市影响范围图。
17、作为对比,匀质空间下取欧氏距离,分别构建点状和面状加权V图。如图5所示,通过生成的V图分析,总体上,太原作为省会城市,位于山西省中部,行政和区位 图 4 山西省坡度图 Fig.4 Slope map of Shanxi Province 优势得天独厚,其综合水平名列榜首,影响范围几乎囊括了整个山西省中部;运城、临汾、长治得分区别不大,占据了晋南大部分地区,晋城得分相对较低,划分影响范围相对较小;大同得分虽然不高,但由于独处晋北,影响范围囊括了整个山西北部,添加地形因素后,划分的影响范围增大;忻州、朔州、阳泉、晋中、吕梁影响范围受限于太原,尤其是晋中得分虽然比大同高,但因为离太原较近,其影响范围
18、要远小于大同;添加地形因素后,由于山脉阻隔,朔州和忻州划分的影响范围增大。叠加对比图5(a)、(b),点状加权V图(呈弧形)不仅不能准确反映城市范围特征,且会将权重小的城市(阳泉、晋中)凸出区域划分到权重大的城市(太原)影响范围内(图6中点状加权V错分区域),可见当城市之间权重差异较大、距离较近时,面状加权V图(呈不规则线段)比点状加权V图更能精确划分城市影响范围;结合图4,对比图5(b)、(c),山西省山区较多,均质空间加权V图不能反映地形特点,非匀质空间加权V图沿地形分布,城市影响范围划分受地形因素影响较大,尤其 414 Journal of Spatio-temporal Informa
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 非匀质 空间 Voronoi 城市 影响 范围 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。