一种基于出租车gps的道路网络提取方法
【技术领域】
[0001] 本发明涉及道路网络提取方法，特别是一种基于出租车GPS的道路网络提取方 法。
【背景技术】
[0002] 中国城市交通发展迅猛，道路建设日新月异。国内尚缺乏可供交通规划师或交通 研究者用于交通分析的权威道路网络数据发布，因此，国内交通规划或研究的现状路网往 往采用规划部门提供的早先勘测的路网信息，或者规划团队花费大量的人力物力进行实地 调研以更新与矫正路网信息。前者获得的路网信息严重滞后；后者则也要花费大量的时间， 获得的数据缺乏了时效性。
【发明内容】
[0003] 发明目的：本发明的目的是提供一种低成本、高效率的基于出租车GPS的道路网 络提取方法。
[0004] 本发明所述的基于出租车GPS的道路网络提取方法，包括如下的步骤：
[0005] S1 :提取一个工作日的出租车GPS数据库：数据库包括8列，即车辆编号、触发事 件、运营状态、GPS时间、GPS经度、GPS炜度、GPS速度和GPS方向；
[0006] S2 :限定识别城市边界：提取城市行政区域边界最西端、最东端、最北端和最南端 的节点坐标，节点坐标包含经度与炜度；设置最西端的经度坐标为城市边界的经度最小值 L0N_、设置最东端的经度坐标为城市边界的经度最大值L0N_、设置最北端的炜度坐标为 城市边界的炜度最大值LAT_、设置最南端的炜度坐标为城市边界的炜度最小值LAT_，四 个边界值围合形成的区域即为观测范围；剔除步骤S1中获得的出租车GPS数据库中不满足 式⑴的数据行；
[0007] L0UGPS 经度 α〇ΝΜΧ，并且 LAUGPS 炜度〈LAI^x (1)
[0008] S3 :数据清洗：将步骤S2处理后的出租车GPS数据库依次按照GPS时间、车辆编 号两个变量升序排列，使得新形成的数据库中同一车辆的数据行连在一起并按时间先后排 列；删除非道路数据行与位置异常数据行，得到有效出租车GPS数据库；设定全局迭代变量 iteration为1 ;其中，删除非道路数据行是删除符合式（2)的数据行：
[0009] 判断条件1 :运营状态=2或者3,并且GPS速度=0 (2)
[0010] 删除位置异常数据行是删除DXiX+1满足式（3)时的第x+1行数据：
[0011] 判断条件 2 :?χ,χ+1>νΜΧ · (Tx+1-Tx) (3)
[0012] Dx,x+1 如式（4)所示：
[0014] 其中，DXiX+1为删除了非道路数据行之后得到的数据库中的相邻行两个数据点之间 的距离，V_为城市道路的最大限速值，r为地球的半径，Τχ、Φχ、λχ*删除了非道路数据行 之后得到的数据库中的第X行GPS时间、GPS炜度和GPS经度，Τχ+1、Φχ+1、λ x+1为删除了非 道路数据行之后得到的数据库中的第χ+l行GPS时间、GPS炜度和GPS经度；
[0015] S4 :观测范围的划分：将步骤S2确定的观测范围切割为m行η列的矩形地块，其 中m与η由式（5)与式（6)得到，切割边界的炜度latai和经度Ion ^分别由式（7)与式 (8)得到：
[0020] S5 :交叉口的判别：对步骤S4中得到的矩形地块进行交叉口识别，也即将落在一 个地块中的有效出租车GPS数据库的所有经炜度坐标组成子数据库，对所有子数据库，依 据经度与炜度的跨度，均匀划分为10X10个检测单元，然后根据交叉口判断条件判断每个 检测单元中是否存在交叉口以及交叉口的位置，将所有识别的交叉口的位置记录在交叉口 检测列表中，最后对全局迭代变量iteration进行赋值与判别；其中，交叉口判断条件为： 检测单元中包含的数据行超过40行，并且这些数据行所包含的不重复的"GPS方位"值的个 数大于40 ;
[0021] S6 :基于交叉口的路段识别：将步骤S5得到的交叉口检测列表转化为待检测交叉 口点对表，每个待检测交叉口点对由起始交叉口与终止交叉口组成；对每一组待检测交叉 口点对，判断起始交叉口与终止交叉口之间是否存在道路，若存在，则将起始交叉口与终止 交叉口的坐标记录在路段检测列表中，并删除重复误判路段；
[0022] S7 :输出道路属性表：输出步骤S5得到的交叉口检测列表以及步骤S6得到的路 段检测列表。
[0023] 进一步，所述步骤S5中的交叉口的位置包括交叉口经度和交叉口炜度，交叉口经 度为检测单元中所有数据行的经度的平均值，交叉口炜度为检测单元中所有数据行的炜度 的平均值。
[0024] 进一步，所述步骤S5中对全局迭代变量iteration进行赋值与判别是对全局迭代 变量iteration加1并重新赋值，判断如果全局迭代变量iteration不等于3,则回到步骤 S4,如果全局迭代变量iteration等于3,则进入步骤S6。
[0025] 进一步，所述步骤S6包括以下的子步骤：
[0026] S6. 1 :交叉口检测列表转化为待检测交叉口点对表：根据式（4)计算交叉口检测 列表中的任一交叉口与其余交叉口之间的距离，筛选出该距离小于2km的交叉口点对组成 待检测交叉口点对表，该点对表包含5个变量，即起始交叉口经度lonstart、起始交叉口炜度 latjstiii't、终止父叉口经度l〇nend、终止父叉口韩度latend、两父叉口距1? Ds e;其中，根据式（9) 来计算，" S "代表起始交叉口，" e "代表终止交叉口：
[0028] S6. 2 :对每一组待检测交叉口点对，判断起始交叉口与终止交叉口之间是否存在 道路：首先，限定道路检测区域，该区域由起始交叉口与终止交叉口坐标连线向两侧各偏移 20m的两条直线和分别经过起始交叉口、终止交叉口的这两条直线的垂线围合而成；然后 将检测区域沿交叉口连线方向平均分成个矩形子块，若少于20个连续矩形 子块内不含有任何有效出租车GPS数据库经炜度坐标点，则表明两交叉口之间存在道路， 将该交叉口点对的5个变量记录到道路检测列表中；
[0029] S6. 3 :删除重复误判路段：根据式（10)计算步骤S6. 2获得的道路检测列表中终 止交叉口与起始交叉口连线所形成的角度Θ，其中atan2是已将象限纳入考量的反正切函 数；将道路检测列表先后按照角度Θ与起始交叉口经度lonstart排序，使得相同起始交叉口 的数据行连在一起并按照角度值大小排序；对比相同起始交叉口的相邻数据行的角度Θ， 如果二者差别小于40°，则删除较大距离Dse的数据行；然后将所有数据行中的角度Θ介 于270°与360°之间的值替换为Θ-360°，重新将道路检测列表先后按照角度Θ与起始 交叉口经度l〇nstart排序，对比相同起始交叉口的相邻数据行的角度Θ，如果二者差别小于 40°，则删除较大距离Ds e的数据行；
[0030] Θ = atan2 (latend_latstart，lonend_lonstart) (10) 〇
[0031] 有益效果：与现有技术相比，本发明具有如下的有益效果：
[0032] 1)提取的路网信息更新、更精确：传统通过规划部门提供的现状城市CAD路网提 取道路信息的方法，很大程度上依赖规划部门路网的更新时间；中国绝大多数城市规划部 门路网更新往往与城市规划同步，几年进行一次，势必导致提取的路网信息滞后；而本发明 方法可以根据近期的出租车GPS数据提取，只要有出租车经过的道路均可以被识别，大大 提尚了路网?目息的时效性和精确性；
[0033] 2)运算速度快：传统的道路网络提取方法，在滞后的现状城市CAD基础上，结合大 量的人力实地踏勘才能获得最新的路网信息，整个过程耗时耗力；调查过程可以持续几周 甚至更长时间；而本发明方法利用近期的出租车GPS数据，可在几个小时内提取详细的最 新的路网信息，大大节省了工时，提高了交通相关规划的效率；
[0034] 3)计算机自动运算，减少了人力投入，降低了误差：该发明方法可以完全由计算 机运算得到，既节省了人力调查的投入，又防止了人力参与而产生的人为误差，也大大节省 了数据更新的成本。
【附图说明】
[0035] 图1为本发明的方法流程图；
[0036] 图2为本发明的观测范围划分的示意图；
[0037] 图3为本发明的道路检测区域的示意图。
【具体实施方式】
[0038] 下面结合【具体实施方式】对本发明的技术方案作进一步的介绍。
[0039] 本发明提供的一种基于出租车GPS的道路网络提取方法，如图1所示，包括如下的 步骤：
[0040] S1 :提取一个工作日的出租车GPS数据库：数据库包括8列，即车辆编号、触发事 件、运营状态、GPS时间、GPS经度、GPS炜度、GPS速度和GPS方向，变量的提取格式要求与 含义如式（1)所示：
[0042] S2 :限定识别城市边界：提取城市行政区域边界最西端、最东端、最北端和最南端 的节点坐标，节点坐标包含经度与炜度；设置最西端的经度坐标为城市边界的经度最小值 L0N_、设置最东端的经度坐标为城市边界的经度最大值L0N_、设置最北端的炜度坐标为 城市边界的炜度最大值LAT_、设置最南端的炜度坐标为城市边界的炜度最小值LAT_，四 个边界值围合形成的区域即为观测范围；剔除步骤S1中获得的出租车GPS数据库中不满足 式⑵的数据行；
[0043] L0UGPS 经度 α〇ΝΜΧ，并且 LAUGPS 炜度〈LAI^x (2)
[0044] S3 :数据清洗：将步骤S2处理后的出租车GPS数据库依次按照GPS时间、车辆编 号两个变量升序排列，使得新形成的数据库中同一车辆的数据行连在一起并按时间先后排 列；删除非道路数据行与位置异常数据行，得到有效出租车GPS数据库；设定全局迭代变量 iteration为1 ;其中，删除非道路数据行是删除符合式（3)的数据行：
[0045] 判断条件1 :运营状态=2或者3,并且GPS速度=0 (3)删除位置异常数据行是 删除Dx，x+1满足式（4)时的第x+1行数据：
[0046] 判断条件 2 :?χ,χ+1>νΜΧ · (Tx+1-Tx) (4)DX,X+1 如式（5)所示：
[0048] 其中，DXiX+1为删除了非道路数据行之后得到的数据库中的相邻行两个数据点之间 的距离，V_为城市道路的最大限速值，r为地球的半径，Τχ、Φχ、λχ*删除了非道路数据行 之后得到的数据库中的第X行GPS时间、GPS炜度和GPS经度，Τχ+1、Φχ+1、λ χ+1为删除了非 道路数据行之后得到的数据库中的