jjzjj

【地图匹配(ST-matching)】GPS 轨迹数据预处理——T-Driver数据集【持续更新中】

★ Skyman·Curry_ 2023-12-02 原文

一、数据介绍

1.数据链接

https://download.csdn.net/download/skyman30/87432266?spm=1001.2014.3001.5503https://download.csdn.net/download/skyman30/87432266?spm=1001.2014.3001.5503

2.数据内容

此数据集包含北京市出租车从2008年2月2日到2008年2月8日的GPS轨迹数据,其中共包含10357辆出租车的数据,其中每个文件由出租车ID,时间、经度、纬度构成。该数据集中的轨迹点总数约为1500 万条,轨迹的总距离达到900万公里。其中连续两个轨迹坐标点的平均采样间隔约为 177秒,距离约为623米。这个数据集的每个文件,由出租车ID命名,包含一个出租车的轨迹。以下是一个文件的样例:

id            time            longitude         latitude

1,2008-02-02 15:36:08,116.51172,39.92123

1,2008-02-02 15:46:08,116.51135,39.93883

1,2008-02-02 15:46:08,116.51135,39.93883

1,2008-02-02 15:56:08,116.51627,39.91034 

二、数据预处理

1.合并文件 

首先将所有的轨迹数据合并到一个文件中方便处理

import os

path = '数据所在文件夹路径'
files = [f for f in os.listdir(path) if f.endswith('.txt')]

#保存到merged_file.txt文件中
with open('merged_file.txt', 'w') as outfile:
    for file in files:
        with open(os.path.join(path, file)) as infile:
            outfile.write(infile.read() + '\n')

2.读取文件

利用pandas读取合并后的文件

import pandas as pd
gps_data = pd.read_csv('合并文件路径',names=['id','time','longitude','latitude'])
  • 读取结果

3.排序

按照id号和时间进行排序

gps_data=gps_data.sort_values(by=['id','time']).reset_index(drop=True)
  •  排序后的结果

 4.统计原始数据量

 一共有17662984条gps数据

raw_length = len(gps_data) #值为17662984

5.加入时间戳

 

6.计算采样点之间的时间差

 

7.计算采样点之间的经纬度距离(Haversine距离计算) 

 

8.计算采样点之间的速度

 

9.按照日期划分数据集

三、数据清洗 

1.清除重复数据

2.清除超出地理坐标范围的轨迹

3.删除微小轨迹


【后面内容待更改】 

5. 数据去重

由于数据中重复的数据是无效数据,因此需要去除

gps_data_drop_dup = gps_data.drop_duplicates().reset_index(drop=True)
  •  去重之后的结果

 6.统计去重后的数据量

去重之后共有16325487条数据,去除了1337497条数据(7.5%的数据) 

7.经纬度坐标分布

  • 查看数据中最大与最小的经纬度
max_latitude = max(gps_data_1.latitude)
max_longitude =max(gps_data_1.longitude)
min_latitude =min(gps_data_1.latitude)
min_longitude =min(gps_data_1.longitude)

#运行结果
#(96.06767, 0.0, 255.3, 0.0)

可以由此推断出经纬度坐标数据中存在异常,因此进一步查看情况 

  • 查看经纬度的数据分布 
#查看经度数据分布
gps_data_1.longitude.describe()

#查看纬度分布
gps_data_1.latitude.describe()

 由此推断部分坐标点超出北京坐标范围,因此需要将这部分数据进行删除,为了方便后续研究,这里的范围采取北京是五环的坐标范围。

8.去除超出范围的数据

#如果坐标范围限定在北京市五环内那么范围是116.17 - 116.62,39.83 - 40.05,对gps数据进行筛选
gps_data_2 = gps_data_1[(gps_data_1['latitude']>39.83)&
                        (gps_data_1['latitude']<40.05)&
                        (gps_data_1['longitude']>116.17)&
                        (gps_data_1['longitude']<116.62)]
gps_data_2 = gps_data_2.reset_index(drop=True)
  •  对坐标异常值清洗后的数据

 其中共清除4629807条数据(占26%) 

9.保存数据

#保存数据清理之后的gps数据
gps_data_2.to_csv('文件路径',index=False)

【由于以上处理方式会误删很多点,因此没有往后继续做】 


二、数据预处理【第二版】

三、地图匹配

1.获取路网数据 

路网数据从openstreatmap获取,利用osmnx包下载,获取之后保存为shapefile文件以便地图匹配时需要 

import osmnx as ox

#获取路网数据
beijing_road=ox.graph_from_bbox(40.05,39.83,116.62,116.17,network_type='drive')
#保存为shapefile文件
ox.save_graph_shapefile(beijing_road,'文件路径')

2.可视化路网

#可视化路网
ox.plot_graph(beijing_road,figsize=(15,15),show=False,close=False,node_size=4)

3.将gps数据投影到路网上

  •  先投影一部分数据看看
import matplotlib.pyplot as plt
#1.先将gps点的经纬度提取出来
latitude = gps_data_2.latitude.to_list()
longitude = gps_data_2.longitude.to_list()

#投影
fig,ax = ox.plot_graph(beijing_road,figsize=(15,15),show=False,close=False,node_size=4)
ax.scatter(longitude[:100000],latitude[:100000],s=0.5,alpha=1,c='red')#投影10万个坐标点
plt.show()

  • 投影全部坐标点 

  

很明显,坐标点大部分都偏离了路网 ,由第一个图看出,路网坐标系和gps轨迹点的坐标系应该是一致的,不然会发生大的偏移,因次这里只需要进行地图匹配即可


有关【地图匹配(ST-matching)】GPS 轨迹数据预处理——T-Driver数据集【持续更新中】的更多相关文章

  1. ruby - 解析 RDFa、微数据等的最佳方式是什么,使用统一的模式/词汇(例如 schema.org)存储和显示信息 - 2

    我主要使用Ruby来执行此操作,但到目前为止我的攻击计划如下:使用gemsrdf、rdf-rdfa和rdf-microdata或mida来解析给定任何URI的数据。我认为最好映射到像schema.org这样的统一模式,例如使用这个yaml文件,它试图描述数据词汇表和opengraph到schema.org之间的转换:#SchemaXtoschema.orgconversion#data-vocabularyDV:name:namestreet-address:streetAddressregion:addressRegionlocality:addressLocalityphoto:i

  2. ruby - Ruby 有 `Pair` 数据类型吗? - 2

    有时我需要处理键/值数据。我不喜欢使用数组,因为它们在大小上没有限制(很容易不小心添加超过2个项目,而且您最终需要稍后验证大小)。此外,0和1的索引变成了魔数(MagicNumber),并且在传达含义方面做得很差(“当我说0时,我的意思是head...”)。散列也不合适,因为可能会不小心添加额外的条目。我写了下面的类来解决这个问题:classPairattr_accessor:head,:taildefinitialize(h,t)@head,@tail=h,tendend它工作得很好并且解决了问题,但我很想知道:Ruby标准库是否已经带有这样一个类? 最佳

  3. ruby-on-rails - Capybara-webkit 引发 Capybara::Driver::Webkit::WebkitInvalidResponseError - 2

    我在rspec中收到来自webkit驱动程序的以下消息:Capybara::Driver::Webkit::WebkitInvalidResponseError:UnabletoloadURL:http://127.0.0.1:44923/posts几天前它成功了。问题出在save_page方法上。有什么问题吗? 最佳答案 当我的页面出现错误时,我收到过类似的错误消息。您应该通过在测试模式下启动服务器(railss-etest)并自行访问页面来手动检查情况是否如此。 关于ruby-on-

  4. ruby - 我如何添加二进制数据来遏制 POST - 2

    我正在尝试使用Curbgem执行以下POST以解析云curl-XPOST\-H"X-Parse-Application-Id:PARSE_APP_ID"\-H"X-Parse-REST-API-Key:PARSE_API_KEY"\-H"Content-Type:image/jpeg"\--data-binary'@myPicture.jpg'\https://api.parse.com/1/files/pic.jpg用这个:curl=Curl::Easy.new("https://api.parse.com/1/files/lion.jpg")curl.multipart_form_

  5. 世界前沿3D开发引擎HOOPS全面讲解——集3D数据读取、3D图形渲染、3D数据发布于一体的全新3D应用开发工具 - 2

    无论您是想搭建桌面端、WEB端或者移动端APP应用,HOOPSPlatform组件都可以为您提供弹性的3D集成架构,同时,由工业领域3D技术专家组成的HOOPS技术团队也能为您提供技术支持服务。如果您的客户期望有一种在多个平台(桌面/WEB/APP,而且某些客户端是“瘦”客户端)快速、方便地将数据接入到3D应用系统的解决方案,并且当访问数据时,在各个平台上的性能和用户体验保持一致,HOOPSPlatform将帮助您完成。利用HOOPSPlatform,您可以开发在任何环境下的3D基础应用架构。HOOPSPlatform可以帮您打造3D创新型产品,HOOPSSDK包含的技术有:快速且准确的CAD

  6. FOHEART H1数据手套驱动Optitrack光学动捕双手运动(Unity3D) - 2

    本教程将在Unity3D中混合Optitrack与数据手套的数据流,在人体运动的基础上,添加双手手指部分的运动。双手手背的角度仍由Optitrack提供,数据手套提供双手手指的角度。 01  客户端软件分别安装MotiveBody与MotionVenus并校准人体与数据手套。MotiveBodyMotionVenus数据手套使用、校准流程参照:https://gitee.com/foheart_1/foheart-h1-data-summary.git02  数据转发打开MotiveBody软件的Streaming,开始向Unity3D广播数据;MotionVenus中设置->选项选择Unit

  7. 使用canal同步MySQL数据到ES - 2

    文章目录一、概述简介原理模块二、配置Mysql使用版本环境要求1.操作系统2.mysql要求三、配置canal-server离线下载在线下载上传解压修改配置单机配置集群配置分库分表配置1.修改全局配置2.实例配置垂直分库水平分库3.修改group-instance.xml4.启动监听四、配置canal-adapter1修改启动配置2配置映射文件3启动ES数据同步查询所有订阅同步数据同步开关启动4.验证五、配置canal-admin一、概述简介canal是Alibaba旗下的一款开源项目,Java开发。基于数据库增量日志解析,提供增量数据订阅&消费。Git地址:https://github.co

  8. ruby-on-rails - 创建 ruby​​ 数据库时惰性符号绑定(bind)失败 - 2

    我正在尝试在Rails上安装ruby​​,到目前为止一切都已安装,但是当我尝试使用rakedb:create创建数据库时,我收到一个奇怪的错误:dyld:lazysymbolbindingfailed:Symbolnotfound:_mysql_get_client_infoReferencedfrom:/Library/Ruby/Gems/1.8/gems/mysql2-0.3.11/lib/mysql2/mysql2.bundleExpectedin:flatnamespacedyld:Symbolnotfound:_mysql_get_client_infoReferencedf

  9. STM32读取串口传感器数据(颗粒物传感器,主动上传) - 2

    文章目录1.开发板选择*用到的资源2.串口通信(个人理解)3.代码分析(注释比较详细)1.主函数2.串口1配置3.串口2配置以及中断函数4.注意问题5.源码链接1.开发板选择我用的是STM32F103RCT6的板子,不过代码大概在F103系列的板子上都可以运行,我试过在野火103的霸道板上也可以,主要看一下串口对应的引脚一不一样就行了,不一样的就更改一下。*用到的资源keil5软件这里用到了两个串口资源,采集数据一个,串口通信一个,板子对应引脚如下:串口1,TX:PA9,RX:PA10串口2,TX:PA2,RX:PA32.串口通信(个人理解)我就从串口采集传感器数据这个过程说一下我自己的理解,

  10. SPI接收数据异常问题总结 - 2

    SPI接收数据左移一位问题目录SPI接收数据左移一位问题一、问题描述二、问题分析三、探究原理四、经验总结最近在工作在学习调试SPI的过程中遇到一个问题——接收数据整体向左移了一位(1bit)。SPI数据收发是数据交换,因此接收数据时从第二个字节开始才是有效数据,也就是数据整体向右移一个字节(1byte)。请教前辈之后也没有得到解决,通过在网上查阅前人经验终于解决问题,所以写一个避坑经验总结。实际背景:MCU与一款芯片使用spi通信,MCU作为主机,芯片作为从机。这款芯片采用的是它规定的六线SPI,多了两根线:RDY和INT,这样从机就可以主动请求主机给主机发送数据了。一、问题描述根据从机芯片手

随机推荐