一种基于POI数据的城市商业中心识别方法及系统与流程

本发明涉及到城市地理学和城市规划技术领域，具体涉及一种基于poi数据，结合道路交点和网络核密度估计法，识别城市商业中心的方法及系统。

背景技术：

作为现代城市最重要的功能区之一，城市商业中心的界定及其空间结构历来是城市地理学研究的重要内容。然而在复杂的地理环境中，城市商业中心是一个模糊的空间对象，通过影像或土地利用现状数据判读和边界数字化获取的边界带有模糊不确定性。关于如何定量化地从高阶模糊地理现象中识别出城市商业中心一直广受关注。目前已有的城市商业中心识别方法通常可以分为两大类：(1)基于问卷调查的方法和(2)基于地理数据的方法。利用第一种方法从市民感知中识别的城市商业中心，受调查对象的主观认知影响大，调查周期耗时长。相对而言，采用第二种方法识别商业中心更加客观和可行。

网络地图兴趣点(pointofinterest,poi)作为一种地理大数据，与传统的地理数据相比，其现势性、准确性、共享性和多分类的特点为研究员提供了更多的研究价值。例如，liu和long利用openstreetmap(osm)和兴趣点(poi)数据，提出了一种基于矢量元胞自动机模型的地块识别的方法；jiang等开发了一种利用兴趣点(poi)数据进行土地利用分类的方法；zhang等基于层次语义认知方法，提出了一种利用高分辨率(hsc)卫星图像和兴趣点(poi)数据的城市功能区划分方法。一些学者根据商业类poi点的地理分布特征，对城市商业中心空间进行了一系列的探索。例如，chen等基于核密度估计法提取出了城市商业中心，并分析了不同零售业态对商业聚集区位选择的差异性；yu等针对网络距离计算问题的特殊性，提出了一种基于网络配置的cbd划分方法。

核密度估计法作为空间密度分析方法的一种，常被用于探测空间的“热点区”，且在poi数据分析中得到了广泛的应用。平面核密度估计法将城市空间环境视为均质的，没有考虑到城市的发展限制于道路网络布局的事实，在城市空间中是网络距离而不是欧氏距离影响着社会和经济现象。因此，一些学者提出了利用网络核密度估计法代替平面核密度估计法来计算道路路段周围的密度分布，剖析道路交通事故和城市空间布局。路段长度的确定是网络核密度估计法研究的关键技术，目前道路路段的划分方法主要有以下两种：(1)等长度道路路段和(2)不等长度道路路段。

虽然网络核密度估计法可以弥补平面核密度估计法的缺陷，但是无论是等长度路段还是不等长度路段的划分都存在一定的缺陷。当路段长度太小，凸显出的局部变化过多；而当路段长度太大，很多局部变化会被忽略。虽然经验法是等长度路段划分的常见方法，但是由于城市的特殊性，经验值并不可能适用于所有城市，存在局限性。利用不等长度道路路段得到的热点位置是多变的，在空间分析中这种变化被称为可塑性面积单元问题。且不管是传统的网络核密度估计还是一些改进的网络核密度估计，其算法过程都是相当复杂。

技术实现要素：

发明目的：针对现有技术中存在的问题，本发明目的在于一种基于poi数据的城市商业中心识别方法及系统，利用道路交点代替道路路段构建密度面，在顾及道路网络对城市商业中心发展作用的同时解决网络核密度估计法中路段划分的问题。

技术方案：为实现上述发明目的，本发明所述的一种基于poi数据的城市商业中心识别方法，包括如下步骤：

(1)获取指定区域的网络地图poi数据，根据城市建设用地分类标准和网络地图poi数据分类标准之间的对应关系提取出商业类poi；

(2)将所指定区域的道路网络数据转换为道路网格数据，提取道路交点；

(3)识别不同网格内的所有商业类poi，将与道路交点相连网格内的商业poi作为道路交点邻域内的商业poi点；

(4)根据每个道路交点邻域内的商业poi点分布情况计算各道路交点的商业核密度值；

(5)根据各带宽内道路交点的分布并结合道路交点的商业核密度值计算各道路交点的商业-道路交点核密度值；

(6)基于商业-道路交点核密度值构建商业-道路交点核密度面；

(7)结合核密度阈值及城市商业中心面积指标，识别出城市商业中心。

作为优选，所述步骤(4)中和步骤(5)中采用高斯核函数计算商业核密度和商业-道路交点核密度。

作为优选，所述步骤(6)中采用克里金插值法基于核密度值构建核密度面。

作为优选，所述步骤(7)中基于商业-道路交点核密度面生成生成密度等值线，利用三倍标准差法识别出候选城市商业中心，再结合商业中心面积指标，从候选城市商业中心中提取出城市商业中心。

本发明所述的一种基于poi数据的城市商业中心识别系统，包括：

商业类poi提取模块，用于获取指定区域的网络地图poi数据，根据城市建设用地分类标准和网络地图poi数据分类标准之间的对应关系提取出商业类poi；

道路交点提取模块，用于将所指定区域的道路网络数据转换为道路网格数据，提取道路交点；

道路交点领域确定模块，用于识别不同网格内的所有商业类poi，将与道路交点相连网格内的商业poi作为道路交点邻域内的商业poi点；

商业核密度计算模块，用于根据每个道路交点邻域内的商业poi点分布情况计算各道路交点的商业核密度值；

商业-道路交点核密度计算模块，用于根据各带宽内道路交点的分布并结合道路交点的商业核密度值计算各道路交点的商业-道路交点核密度值；

核密度面构建模块，用于基于商业-道路交点核密度值构建商业-道路交点核密度面；

以及，城市商业中心界定模块，用于结合核密度阈值及城市商业中心面积指标，识别出城市商业中心。

本发明所述的一种基于poi数据的城市商业中心识别系统，至少包括一台计算机设备，所述计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的城市商业中心识别方法。

有益效果：本发明可以解决城市商业中心范围的模糊不确定性，能够更客观、科学的反映城市商业中心范围与商业网络点和交通通达性之间的关联，有助于进一步分析城市商业中心的空间结构，为城市规划实施评估提供科学的基础。本发明的市级城市商业中心的识别数量和吻合精度优于平面核密度估计法，且与网络核密度估计法相比，本发明降低的算法的复杂度。并且利用本发明改进后的核密度估计法还可以提前探测出道路发达但是商业并不繁华的潜在商业中心。

附图说明

图1是本发明实施例的方法流程示意图。

图2是本发明实施例中道路交通点邻域的构建示意图。

图3是本发明实施例中构建的南京市商业-道路交点核密度面。

图4是本发明实施例中商业-道路交点核密度面提取结果图。

具体实施方式

下面结合附图和具体实施例对本发明的技术方案做进一步说明。为了便于说明，本发明实施例中的网络地图poi数据来源于高德poi数据，可以理解的是，本发明并不限于该数据来源。

本发明实施例以高德网络地图poi数据为研究对象，提供了一种基于poi数据的城市商业中心识别方法，流程如图1所示，主要包括如下步骤：

s1：获取指定区域的网络地图poi数据，根据城市建设用地分类标准和网络地图poi数据分类标准之间的对应关系提取出商业类poi。

网络地图poi数据是一种以点形式表达地理实体的数据，为迎合公众的兴趣，更关注生活服务类的分类，如餐饮、酒店、休闲娱乐、医院等。而在城市规划中，规划人员根据用地的属性，将地理实体划分为商业服务业设施用地(b)、居住用地(r)、工业用地(m)等8大类。因此，为了从网络地图poi数据中提取出城市商业中心，必须对poi数据进行重分类，识别出商业类poi点。对比城市建设用地和网络地图poi数据的分类标准，商业类poi与商业服务业设施用地(b)之间主要存在以下两种关系：

(1)“包含关系”：某一类poi数据完全属于商业服务业设施用地(b)，这类poi可以直接被划分为商业类poi点。

(2)“介于关系”：某一类poi数据的一部分属于商业服务业设施用地(b)，另一部分属于其它类型的用地。所以需要从此类poi数据中抽取出商业类poi点。

本实施例中选择南京市作为研究区域，结合以上的分析，得到南京市高德商业poi重分类表：

表1.基于建设用地的南京市高德商业poi重分类表

s2：将所指定区域的道路网络数据转换为道路网格数据，提取道路交点，即每个道路网格的顶点。

s3：识别不同网格内的所有商业类poi，将与道路交点相连网格内的商业poi作为道路交点邻域内的商业poi点。

传统平面核密度估计法的邻域一般由带宽决定，且常用的经验法选择带宽长度。在本实施例中，将与道路交点邻近的道路网格视为邻域。各道路交点的商业核密度受其邻域内商业pois分布的影响，靠近道路交点的商业pois越多，其商业核密度越高。图2中的黑色加粗边框即为道路交点s的邻域边界。

具体地，本步骤中首先识别每个网格的顶点，并建立二维矩阵，通过网格顶点的id号将道路交点和道路网格关联起来，每个道路交点的邻域即为与之相连的网格。根据每个商业poi点与不同格网的空间关系(包含关系)，再识别不同网格内的所有商业类poi点，则与道路交点相连网格内的商业poi即为道路交点领域内的商业poi点。

s4：根据每个道路交点邻域内的商业poi点分布情况计算各道路交点的商业核密度值。

商业核密度是评价空间商业密度的一个值，它受邻域划分的影响。每个道路交点的商业核密度km可以按照公式1和2计算得到，公式如下：

其中，r是道路交点到周围邻域内商业poi点的最远距离(带宽)，di是邻域内第i个商业点到道路交点的距离，k(di，r)是高斯核函数。

其中，km是第m个道路交点的商业核密度(概率密度估计值)，n是道路交点邻域内商业点的数量。当km值越大，表示在道路交点处聚集了越多的商业pois。

s5：根据各带宽内道路交点的分布并结合道路交点的商业核密度值计算各道路交点的商业-道路交点核密度值。

事实上，km仍然是描述商业密度的一个值，但它是基于道路交点计算的。在本实施例中，计划不仅从商业密度的角度来辨识商业中心区，并且加入道路网络作为一个约束条件，基于网络来识别商业中心区。然后，将km值作为重要的参数来计算道路交点的核密度是一个综合性指标，同时表示了道路网络和商业pois的聚集度。利用公式3和4计算道路交点s的商业-道路交点核密度，公式如下：

其中，r是实验选择的带宽，di是带宽内的第i个道路交点和道路交点m之间的距离，km是道路交点m的商业核密度，n是带宽内道路交点的数量，是道路交点的核密度。符合距离衰减规律，当商业poi点越靠近道路交点，道路交点的核密度较大。此外，与km成正比例关系。

s6：基于商业-道路交点核密度值构建商业-道路交点核密度面。

是一个受商业核密度和道路交点核密度共同约束的指标，本实施例采用arcgis中的克里金插值法，基于核密度值构建核密度面，并将此面定义为商业-道路交点核密度面。图3为利用以上6个步骤，带宽设为300米构建的南京市商业-道路交点核密度面。图中，颜色越深的栅格面表示其商业-道路交点核密度值越大，成为商业中心面的可能性越高。

s7：结合核密度阈值及城市商业中心面积指标，识别出城市商业中心。

为了能进一步界定南京市城市商业中心范围，基于图3生成密度等值线，然后利用三倍标准差法识别出候选城市商业中心，再结合每个城市的商业中心面积指标，从候选城市商业中心中提取出城市市级商业中心。

在《南京市商业网点规划建设管理办法》中规定，市级商业中心区的面积必须大于30万平方米，因此必须剔除面积小于30万平方米的候选商业中心区，保留下的热点区即为城市商业中心的范围。图4为南京市市级商业中心的识别结果。图中，深黑色区域是利用本实验方法提取的市级商业中心，浅色区域是规划中的先验区，深灰色区域为实验提取的非市级商业中心。

为验证本发明方法的有效性及优势，将识别结果与每个城市的市级商业系统规划图进行对比，并计算每个商业面的f1-score指标，评价商业中心区的吻合度，提取的商业中心区与先验范围越一致，则f1-score值较高。公式如下：

其中acomputed是利用研究方法识别出的城市商业中心区的范围，acomparative是先验范围，aoverlap是两个范围间的重叠面积。

表2.提取结果的精度对比

如表2所示，将本发明方法与平面核密度估计法的识别结果对比可以发现，本发明的方法精度更高，可以识别出潜在城市商业中心，仙林商业中心。

本发明另一实施例公开的一种基于poi数据的城市商业中心识别系统，包括：商业类poi提取模块，用于获取指定区域的网络地图poi数据，根据城市建设用地分类标准和网络地图poi数据分类标准之间的对应关系提取出商业类poi；道路交点提取模块，用于将所指定区域的道路网络数据转换为道路网格数据，提取道路交点；道路交点领域确定模块，用于识别不同网格内的所有商业类poi，将与道路交点相连网格内的商业poi作为道路交点邻域内的商业poi点；商业核密度计算模块，用于根据每个道路交点邻域内的商业poi点分布情况计算各道路交点的商业核密度值；商业-道路交点核密度计算模块，用于根据各带宽内道路交点的分布并结合道路交点的商业核密度值计算各道路交点的商业-道路交点核密度值；核密度面构建模块，用于基于商业-道路交点核密度值构建商业-道路交点核密度面；以及，城市商业中心界定模块，用于结合核密度阈值及城市商业中心面积指标，识别出城市商业中心。该系统实施例与上述方法实施例属于相同的发明构思，具体实现细节可参考上述方法实施例，此处不再赘述。

基于相同的发明构思，本发明实施例还公开了一种基于poi数据的城市商业中心识别系统，至少包括一台计算机设备，所述计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的城市商业中心识别方法。本发明实施例未详细说明的部分均为现有技术。