栅格成本距离分析法以格网为单元，基于路网栅格数据测度成本栅格数据，再计算各栅格单元到源(中心城市)累计耗费距离的区域可达性，再根据可达性，确定栅格归属的中心城市，最终得到城市影响腹地范围。相对于矢量网络分析法，该方法能比较好的模拟现实。
02 数据准备珠三角地级市矢量图；
珠三角公路网数据(数据来源OSM数据)。
图1｜珠三角地级市矢量图
图2｜珠三角公路网数据
数据获取可参考历史文章：
【粤港澳大湾区专题】01基础地理数据·整理分享
03 操作步骤◐ 第一步
选择【系统工具箱→Conversion Tools→转为栅格→要素转栅格】工具，将广东省路网要素转化为1km×1km的栅格，字段选择公路等级字段。处理后得到的栅格数据仅仅是表示公路等级的栅格，不具备成本意义，需要对其进行重新赋值。
选择【系统工具箱→Spatial Analyst Tools→重分类→重分类】工具，进行参数设置，得到成本栅格数据。
实际运行速度为：高速——100km/h；国道(一级道路)——75km/h；省道(二级道路)——60km/h；县道(三级道路)——40km/h；其他道路(四级道路)——20km/h；
根据速度计算得到不同级别道路的每公里时间成本：高速——36sec；国道(一级道路)——48sec；省道(二级道路)——60sec；县道(三级道路)——90sec；其他道路(四级道路)——180sec，空白处按照600sec。
选择【系统工具箱→Spatial Analyst Tools→距离分析→成本距离】工具，计算得到距离栅格数据。
因为此时得到的时间是以秒为单位，可以通过选择【系统工具箱→Spatial Analyst Tools→地图代数→栅格计算器】工具，计算得到以小时为单位的结果。
主要是根据区域范围内各点在累加成本最小的基础上选择最近的源，从而可以分配各源的点集范围。选择【系统工具箱→Spatial Analyst Tools→距离分析→成本分配】工具，计算得到分配栅格数据。
1.反映了城市之间的联系程度，城市间的交通便捷程度是城市间交流活跃度的重要指标。
2.区域同城化，以广佛同城化为例，从城市空间结构或行政区边界来看，广州市西部区域距离中心较近、地域结构不完整，佛山市中心位置不明确，广佛两地连在一起才能构成一个完整的核心-边缘结构，成为一个完整的地理单元。
3.城市腹地范围与行政区域相差较大。长期以来，中国城市以行政区域内部相对封闭的管理模式为主，难以实现跨市或跨区域的共建发展。随着经济快速发展和城镇化进程加快，城市之间存在空间隔离和行政制度障碍，社会经济等功能联系却不断增强；特别是一些地理位置上相近的城市，在空间融合、产业协作、设施衔接等方面均有较大的发展。导致一些地区在经济联系上与其它地区联系密切。
4.在一些地区出现多个地区的腹地范围，这些区域往往位于城市交界处，而远离所在城市中心区域。在粤港澳大湾区空间一体化的背景下。内部城市由于其距离近、发展关系密切,城市交界地带的建设对城市发展至为敏感,成为多方博弈的热点地区。
编辑批注：阿松宝贝，你可真能扯，像一个科研人了。
不足之处：
1.未引用城市相关经济社会人文数据进行优化，下一步引进社会驱动力等因子，在多因子尺度下进行数学建模是研究重点。
2.城市点选择还需要综合更多影响因素，多重验证城市点选择的科学性，可引用经济数据构建各城市经济重心，以经济重心进行研究，优化模型。
3.源的数据具有单一性，城市点源属性不具备反映城市综合影响力的范围，在进行成本分配时，均匀分配给城市点，忽略了城市吸引力。下一步需对城市赋予权重因子，构建城市吸引力模型进行验证。
作者｜松玄月下
排版｜数读菌
校阅｜数读菌
那今天就到这里结束啦，欢迎留言讨论。文中的图片文字未经许可不要随便“引用”。
如果可以的话，希望能够转发分享，点个在看并且点个赞，给个赞赏~~也欢迎规范转载~
也希望大家和我多留言互动啊！(据说这样可以增加我的推送在你的订阅号里出现的概率)
历史文章推荐：
【数据分享】中国生态服务功能重要性与中国生态系统敏感性
【数据整理】2019年行政区划调整的数据更新与分享
【数据分享】如何获取北京旅游景区AOI数据？
【数据分享】如何获取全国矢量地铁站点与线路数据？
【数据分享】扩张的世界，收缩的城市(分享世界人口密度数据)
【数据分享】地形是如何影响人口与城市分布的？(图多！流量慎点)
【粤港澳大湾区专题】01基础地理数据·整理分享
【长三角一体化专题】01基础地理数据·整理分享
【数据分享·教程】20分钟做一个中国城市经济联系图
【GIS基础】基本空间分析工具：缓冲区、网络分析、相交分析、密度分析
【GIS基础】基于路网数据的最短路径分析(附练习数据)
【GIS基础】填挖方分析
【GIS基础】谁都能学会的地形分析全流程(1)
【GIS基础】谁都能学会的地形分析全流程(2)
【GIS基础】谁都能学会的地形分析全流程(3)
【ENVI实验课】01基于像元二分模型的植被覆盖度反演
【ENVI实验课】02利用ENVI进行农业耕作用地变化监测
【GIS实验课】01山区土壤侵蚀评价及其地形的关系——以皖西大别山为例
【GIS实验课】02 超市选址适宜性评价
【GIS实验课】03 城镇住房选址适宜区评价分析
【GIS实验课】04 GIS数据的批处理方法
【GIS实验课】05 坡度提取的不确定性分析(分享批量重采样、批量掩膜、坡度批量提取代码)
▼需要你的“分享”和“在看”
END>
如需全文转载文章、投稿或者合作
可添加微信
(回复超慢！！！)
(不要添加我问各种问题，我大概率不会的==)
(入群请一定要备注入群)
(添加后会在晚上非工作时间通过，请稍安勿躁)
公众号
微博
