一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法技术

本发明专利技术涉及一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法。通过提取高空间分辨率遥感影像的加权均值视觉特征和出租车轨迹上下车点数据的时序统计特征进行城市功能区的识别。在提取遥感影像的视觉特征时，通过计算各基本分析单元的均值权重，能有效缓解样本的不均衡问题，提高精度；综合使用高分遥感影像的视觉特征和出租车轨迹数据的统计特征，充分发挥二者优势，从自然属性和社会属性两方面表征城市用地的属性信息，能有效对城市功能区进行分类，与仅使用视觉特征或者时序统计特征的方法相比，本发明专利技术提出的方法得到精确度、召回率和F1

【技术实现步骤摘要】
一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法


[0001]本专利技术属于城市地理学
，特别是涉及一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法。

技术介绍

[0002]城市用地是城市规划区范围内赋以一定用途和功能的土地的统称，包含自然属性、社会属性、经济属性和法律属性等四大用地属性。城市空间结构主要包括城市自然要素和社会要素的空间分布模式以及要素之间的相互关系。城市空间结构形态的研究在西方社会出现较早，作为系统的理论研究则是在城市规划实践的推动下逐步产生和发展起来的。1923年，伯吉斯以芝加哥为例，基于社会生态学的入侵和继承的概念，创立了城市发展与土地使用空间组织方式的同心圆模式。 1932年，巴布科引入交通轴线的辐射作用，将同心圆模式修正为星状环形模式，使之更接近单中心的中小规模城市的真实状况。而我国在20世纪80年代前，城市空间结构与形态只有一些零星的研究成果，真正意义上的研究于80年代之后才开始。改革开放以来，我国经济发展迅猛，城市规模不断扩张，城市人口暴涨，出现了一系列城市病，如何合理地规划和管理城市空间结构成为一个亟待解决的问题。
[0003]当前，各级土地管理部门编制土地利用现状图主要依据遥感影像解译和建设用地审批信息的内业分析、外业调查核实结果等进行制作，人力物力消耗较高，更新周期较长。其中，使用遥感影像进行土地利用类型分类的研究，起初关注于提取基于像元的光谱特征和纹理特征，之后学者们加入了目标地物的结构信息，基于面向对象的思想进行遥感影像分类，近来又引入了图像的上下文语义信息，从场景理解的角度进行图像分类。然而，上述研究仅反映了城市的自然地理环境，即城市构成要素的空间组成和空间配置，大部分的人类活动都是在建筑物内进行，仅依靠遥感影像无法判断室内活动区域的社会经济功能属性。近些年来，出租车轨迹数据和社交媒体数据等城市居民时空行为数据的涌现，提供了大样本量的城市居民连续型经济、社会活动信息，不仅可以用来进行个体活动模式的研究，也为宏观形态的城市土地利用类型的空间分布模式提供了社会经济环境信息。出租车轨迹数据凭借其数据量大、覆盖率高、成本低、更新速度快等优点，目前已广泛应用于交通状态的评估、移动模式探测、热点区域提取、道路网络更新、土地利用分类等研究领域。

技术实现思路

[0004]本专利技术针对现有技术的不足，提供一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，通过提取高空间分辨率遥感影像的加权均值视觉特征和出租车轨迹上下车点数据的时序统计特征进行城市功能区的识别。
[0005]为了达到上述目的，本专利技术提供的技术方案是一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，包括以下步骤：
[0006]步骤1，基于OpenStreetMap(OSM)道路数据，建立城市用地的基本分析单元，并注
明基本分析单元的建设用地类型属性，生成土地利用图，包括以下几个子步骤：
[0007]步骤1.1，下载OSM城市底图及相关属性数据的shapefile文件，包括城市道路、建筑物、自然用地、重要地点、铁路、水系等信息，根据区域行政边界矢量图裁剪出研究区域；
[0008]步骤1.2，提取OSM中区域城市道路数据的道路中心线，并参考出租车轨迹点数据聚集产生的轨迹路径，进行拓扑检查与修正，生成区域道路网矢量图层；
[0009]步骤1.3，利用步骤1.2得到的区域道路网矢量图层切割区域行政边界矢量图，生成基本分析单元面图层文件；
[0010]步骤1.4，参考百度在线地图和高德在线地图的城市建设用地实际类型，为相应的基本分析单元添加建设用地类型属性，生成土地利用图，土地利用混合程度过高的区域不在本研究考虑之列；
[0011]步骤2，基于高空间分辨率遥感影像提取加权均值的视觉特征，包括以下几个子步骤：
[0012]步骤2.1，对卫星原始遥感影像数据进行预处理，包括几何校正、图像融合、图像镶嵌、图像裁剪；
[0013]步骤2.2，依次使用棋盘分割、多尺度分割和光谱差分割三种多层次图像分割方法，基于高空间分辨率遥感影像和道路矢量图，得到50184个图像对象；
[0014]步骤2.3，提取各个图像对象包括光谱特征、纹理特征和几何特征在内的19 种视觉特征，得到每个图像对象的19维视觉特征向量；
[0015]步骤2.4，以各基本分析单元为单位，对其包含的所有图像对象的视觉特征值求加权平均值得到各基本分析单元的19维加权平均视觉特征向量；
[0016]步骤3，基于出租车轨迹数据，提取上、下车点数据的时序统计特征，包括以下几个子步骤：
[0017]步骤3.1，对一个月的出租车轨迹数据进行预处理，包括坐标系纠偏、研究区域裁剪、异常点清理；
[0018]步骤3.2，从步骤3.1经过预处理的出租车轨迹数据中提取上、下车点；
[0019]步骤3.3，将步骤3.2提取的上、下车点图层匹配到步骤1.3中的基本单元面图层；
[0020]步骤3.4，提取各基本分析单元内一个月的出租车上、下车点数据的48维时序统计特征；
[0021]步骤4，综合使用高空间分辨率遥感影像的视觉特征和出租车轨迹上下车点数据的时序统计特征进行城市功能区的识别。
[0022]而且，所述步骤1.2中拓扑检查可利用ArcMap软件中的“拓扑”工具进行，拓扑修正可根据出租车的轨迹路径对OSM中提取的区域城市道路数据的道路中心线进行手动修正，出租车轨迹数据由装置在出租车上的GPS设备接收导航卫星发射的信息采集而成，采样时间间隔大于40s，每条记录包含车牌号、时间、经纬度、方向、速度、ACC状态和载客状态等属性信息。
[0023]而且，所述步骤2.4中各基本分析单元的19维加权平均视觉特征向量，计算方法如下：
[0024][0025]式中，为第i个基本分析单元的19维加权平均视觉特征向量，L
ij
为第i 个基本分析单元内第j个图像对象的19维视觉特征向量，m为第i个基本分析单元内图像对象的数量，p
i
为第i个基本分析单元内各图像对象的权重值，计算方法如下：
[0026][0027]式中，n
i
表示第i个基本分析单元中图像对象数目，N表示研究区中总图像对象的数目。
[0028]而且，所述步骤3.3需先将步骤1.3中基本分析单元面图层转化为线图层，并为线图层添加线段起点、终点、中点、方向4种几何属性，然后将点图层根据最邻近距离匹配到线图层，最后根据点图层的方向属性和基本分析单元线图层的方向属性，将轨迹点匹配到线两侧的基本分析单元内，实现点图层与基本分析单元面图层相匹配。
[0029]而且，所述步骤3.4需先统计各基本分析单元一个月内出租车每天24小时每小时内上、下车点的总量，构成每天上车点24维时序统计特征向量和下车点 24维时序统计特征向量，然后再将这些24维时序统计特征分为工作日和周末两类，工作日包括周一到周五5天，周末包括周六和周本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，其特征在于，包括如下步骤：步骤1，基于OpenStreetMap道路数据，建立城市用地的基本分析单元，并注明基本分析单元的建设用地类型属性，生成土地利用图，包括以下几个子步骤：步骤1.1，下载OSM城市底图及相关属性数据的shapefile文件，包括城市道路、建筑物、自然用地、重要地点、铁路、水系信息，根据区域行政边界矢量图裁剪出研究区域；步骤1.2，提取OSM中区域城市道路数据的道路中心线，并参考出租车轨迹点数据聚集产生的轨迹路径，进行拓扑检查与修正，生成区域道路网矢量图层；步骤1.3，利用步骤1.2得到的区域道路网矢量图层切割区域行政边界矢量图，生成基本分析单元面图层文件；步骤1.4，参考城市建设用地实际类型，为相应的基本分析单元添加建设用地类型属性，生成土地利用图；步骤2，基于高空间分辨率遥感影像提取加权均值的视觉特征，包括以下几个子步骤：步骤2.1，对卫星原始遥感影像数据进行预处理；步骤2.2，使用多层次图像分割方法，基于高空间分辨率遥感影像和道路矢量图，得到M个图像对象；步骤2.3，提取各个图像对象包括光谱特征、纹理特征和几何特征在内的λ种视觉特征，得到每个图像对象的λ维视觉特征向量；步骤2.4，以各基本分析单元为单位，对其包含的所有图像对象的视觉特征值求加权平均值，得到各基本分析单元的λ维加权平均视觉特征向量；步骤3，基于出租车轨迹数据，提取上、下车点数据的时序统计特征，包括以下几个子步骤：步骤3.1，对一定时间段内的出租车轨迹数据进行预处理；步骤3.2，从步骤3.1经过预处理的出租车轨迹数据中提取上、下车点；步骤3.3，将步骤3.2提取的上、下车点图层匹配到步骤1.3中的基本单元面图层；步骤3.4，提取各基本分析单元内一定时间段内的出租车上、下车点数据的48维时序统计特征；步骤4，综合使用高空间分辨率遥感影像的视觉特征和出租车轨迹上、下车点数据的时序统计特征进行城市功能区的识别。2.如权利要求1所述的一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，其特征在于：所述步骤1.2中拓扑检查可利用ArcMap软件中的“拓扑”工具进行，拓扑修正可根据出租车的轨迹路径对OSM中提取的区域城市道路数据的道路中心线进行手动修正，出租车轨迹数据由装置在出租车上的GPS设备接收导航卫星发射的信息采集而成，采样时间间隔大于40s，每条记录包含车牌号、时间、经纬度、方向、速度、ACC状态和载客状态属性信息。3.如权利要求1所述的一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，其特征在于：所述步骤2.1中预处理包括几何校正、图像融合、图像镶嵌、图像裁剪。4.如权利要求1所述的一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，其特征在于：所述步骤2.2是依次使用棋盘分割、多尺度分割和光谱差分割三种多层次
图像分割方法，基于高空间分辨率遥感影像和道路矢量图，得到M个图像对象。5.如权利要求1所述的一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，其特征在于：所述步骤2.3中视觉特征包括归一化差分植被指数、归一化差分水体指数、土壤调整植被指数、灰度共生矩阵均值、熵、对比度、相关性、非相似度、均值、标准差、偏度、面积、长宽比、边界指数、紧致度、椭圆拟合、矩形拟合、圆度、形状指数。6.如权利要求1所述的一种综合使用遥感数据和出租车轨迹数据的城市功能区识别方法，其特征在于：所述步骤2.4中各基本分析单元的λ维加权平均视觉特征向量，计...

【专利技术属性】
技术研发人员：秦昆，张晔，
申请(专利权)人：武汉大学，
类型：发明
国别省市：