CN109269514A - 确定运动轨迹的方法和装置 - Google Patents

确定运动轨迹的方法和装置 Download PDF

Info

Publication number
CN109269514A
CN109269514A CN201710586941.1A CN201710586941A CN109269514A CN 109269514 A CN109269514 A CN 109269514A CN 201710586941 A CN201710586941 A CN 201710586941A CN 109269514 A CN109269514 A CN 109269514A
Authority
CN
China
Prior art keywords
point
filtered
initial trace
road
trace point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710586941.1A
Other languages
English (en)
Inventor
杨力
孙未未
吴昊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fudan University
Huawei Technologies Co Ltd
Original Assignee
Fudan University
Huawei Technologies Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fudan University, Huawei Technologies Co Ltd filed Critical Fudan University
Priority to CN201710586941.1A priority Critical patent/CN109269514A/zh
Priority to PCT/CN2018/094740 priority patent/WO2019015485A1/zh
Priority to EP18834716.5A priority patent/EP3640596A4/en
Publication of CN109269514A publication Critical patent/CN109269514A/zh
Priority to US16/746,666 priority patent/US20200151885A1/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/393—Trajectory determination or predictive tracking, e.g. Kalman filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/215—Motion-based segmentation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/14—Receivers specially adapted for specific applications
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40—Correcting position, velocity or attitude
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/207—Analysis of motion for motion estimation over a hierarchy of resolutions
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00—Network arrangements or protocols for supporting network services or applications
    • H04L67/50—Network services
    • H04L67/535—Tracking the activity of the user
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/28—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
    • G01C21/30—Map- or contour-matching
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/50—Determining position whereby the position solution is constrained to lie upon a particular curve or surface, e.g. for locomotives on railway tracks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10016—Video; Image sequence

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Automation & Control Theory (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Navigation (AREA)
  • Image Processing (AREA)

Abstract

本申请提供了一种确定运动轨迹的方法和装置。该方法包括：获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该移动终端经过第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻；将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理用于删除和/或更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置；根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹。本申请实施例的确定运动轨迹的方法和装置，针对基站定位这样的高噪声的数据获取方式的地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。

Description

确定运动轨迹的方法和装置
技术领域
本申请涉及信息技术领域，并且更具体地，涉及确定运动轨迹的方法和装置。
背景技术
地图匹配是一种移动轨迹数据的误差修正技术，该技术以模式识别理论为依据，可以基于车辆始终行驶在道路上的假设，其基本思想是结合车辆定位的轨迹与数字地图中的道路网络，将定位方法所测得的车辆位置信息与导航系统的电子地图数据比较和匹配，找到车辆所在的路段。定位技术结合地图匹配技术，可以有效减少由于定位采样带来的误差，极大地提高了车辆定位精度，现已广泛应用于全球定位系统(global positioningsystem，GPS)导航、交通流分析等领域。
现有的地图匹配技术中，基于隐马尔可夫模型(Hidden Markov Model,HMM)的地图匹配具有最好的准确率和鲁棒性。该方法主要还是针对小误差的轨迹设计的，例如GPS轨迹，其误差在10m的级别，匹配准确率能够达到比较理想的结果。
目前，由于基站定位数据获取比较简单，并且用户手机十分普及，且无需GPS定位那样需要主动开启GPS开关，获得的数据量也较大，蕴含着丰富的信息，因此提出一种基于基站定位数据的地图匹配技术十分必要。若采用现有HMM地图匹配算法，对于基站轨迹数据这种大误差的轨迹，其误差在100m的级别，匹配的准确率则得不到理想的结果。
发明内容
本申请提供了一种确定运动轨迹的方法和装置，能够过滤与修正地图匹配过程中可能起了负面作用的轨迹点。
第一方面，提供了一种确定运动轨迹的方法，该方法包括：获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该移动终端经过该n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻，n为正整数；将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理用于删除该n个原始轨迹点P1至Pn中的至少一个原始轨迹点，和/或更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置；根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹，该道路地图包括多条道路。
因此，本申请实施例的确定运动轨迹的方法，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
应理解，服务器可以获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该n个原始轨迹点P1至Pn可以为基站采集的该若干轨迹点中的全部轨迹点或者部分轨迹点，例如，服务器可以获取该基站采集的全部轨迹点，该n个原始轨迹点P1至Pn即为该全部轨迹点；或者，服务器可以对该全部轨迹点进行处理，该n个原始轨迹点P1至Pn可以为服务器处理原全部轨迹点得到的其中部分轨迹点。
应理解，基站采集移动终端在运动过程中轨迹点，可以采用相同的采样时间，例如，采用相同的时间间隔，获取每个时刻的轨迹点的位置，或者也可以采用不同的采样时间，即在不同的时间间隔，采集移动终端移动的轨迹点的位置。
可选地，作为一个实施例，该方法中的过滤处理可以包括差值过滤、速度过滤、角度过滤、扭曲度过滤、切尾均值过滤和过近点过滤中的至少一个，其中，速度过滤、角度过滤、扭曲度过滤和过近点过滤可以用于删除该n个原始轨迹点P1至Pn中的至少一个原始轨迹点，差值过滤和切尾均值过滤可以用于更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置。
结合第一方面，在第一方面的一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行差值过滤处理，该差值过滤处理用于将位置重合的点进行重新排列，即更新重合点的位置。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：在该n个原始轨迹点P1至Pn中获取a个待过滤点Px至Px+a-1，该a个待过滤点Px至Px+a-1位于第一位置，x和a为正整数；在该n个原始轨迹点P1至Pn中获取b个待过滤点Px+a至Px+a+b-1，该b个待过滤点Px+a至Px+a+b-1位于第二位置，b为大于1的正整数；在该n个原始轨迹点P1至Pn中获取c个待过滤点Px+a+b至Px+a+b+c-1，该c个待过滤点Px+a+b至Px+a+b+c-1位于第三位置，c为正整数，x+a+b+c-1小于或者等于n；在该第一位置到该第二位置的直线上均匀排列待过滤点至在该第二位置到该第三位置的直线上均匀排列待过滤点至 表示向下取整。
应理解，对于a大于1时，该a个待过滤点Px至Px+a-1在该第一位置上完全重合；该b个待过滤点Px+a至Px+a+b-1在第二位置上完全重合；对于c大于1时，该c个待过滤点Px+a+b至Px+a+b+c-1在第三位置上完全重合。
应理解，在该第一位置到该第二位置的直线上均匀排列待过滤点至包括：在第一位置上为待过滤点在第二位置上为待过滤点将第一位置与第二位置之间的直线按照距离等分，依次排列待过滤点至之间的每个待过滤点。
应理解，在该第二位置到该第三位置的直线上均匀排列待过滤点至包括：在第二位置上为待过滤点在第三位置上为待过滤点将第二位置与第三位置之间的直线按照距离等分，依次排列待过滤点至之间的每个待过滤点。
因此，通过差值过滤，可以将重合的轨迹点重新排列，更新部分轨迹点的位置，使得轨迹点分布更加均匀。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行速度过滤处理，该速度过滤处理用于过滤速度较大的点，例如，可以过滤掉速度大于或者等于预设速度的点。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，a为小于或者等于n的正整数；若该第一待过滤点Pa的速度大于或者等于预设速度，删除该第一待过滤点Pa。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，包括：在该n个原始轨迹点P1至Pn中获取该第一待过滤点Pa和第二待过滤点Pb，b＝a+1或者b＝a-1，a小于n；确定从该第一待过滤点Pa到该第二待过滤点Pb的平均速度为该第一待过滤点Pa的速度。
应理解，预设速度可以根据实际情况进行设置，例如，可以根据移动终端在整个运动过程中的平均速度设置，或者是根据移动终端在某一段时间或距离内的速度进行设置，本申请实施例并不限于此。
因此，通过速度过滤，可以将部分误差速度过大的轨迹点过滤掉，这些轨迹的速度不符合常理，与其它轨迹点的速度相差过大，速度过滤可以使得剩余数据更加合理和准备。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行角度过滤处理，该角度过滤用于删除部分角度大于或者等于预设角度的角的顶点。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：在该n个原始轨迹点P1至Pn中获取4个待过滤点Pc、Pc+1、Pc+2和Pc+3，c为正整数，c+3小于或者等于n；确定第一角度∠PcPc+1Pc+2和第二角度∠Pc+1Pc+ 2Pc+3；若该第一角度小于或者等于第一预设角度且该第二角度小于或者等于第二预设角度，删除该4个待过滤点Pc、Pc+1、Pc+2和Pc+3中的待过滤点Pc+1或者待过滤点Pc+2。
应理解，预设角度可以根据实际情况进行设置，其中，第一预设角度和第二预设角度可以设置为相同或者不同，本申请实施例并不限于此。
考虑到移动终端在道路上运动时，其轨迹一般应该为较平滑的线段，如果出现了较小的锐角的情况，则很有可能是采样误差造成的，因此，可以通过角度过滤把误差较大的轨迹点过滤掉。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行扭曲度过滤处理，该扭曲度过滤处理可以删除部分扭曲度过大的点。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：在该n个原始轨迹点P1至Pn中，获取2α+1个待过滤点Pf-α至Pf+α，f和α为正整数，f大于α；确定该2α+1个待过滤点Pf-α至Pf+α中待过滤点Pg与相邻的待过滤点Pg+1之间的距离，g取f-α至f+α-1，获得2α个距离；确定该2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离为第一距离；确定该2α个距离之和与该第一距离的商为扭曲度；当该扭曲度大于或等于预设扭曲度时，删除该2α+1个待过滤点Pf-α至Pf+α中第α+1个待过滤点Pf。
具体地，服务器可以在该n个原始轨迹点P1至Pn中遍历每一个点，这里以任意一个轨迹点Pf为例。以该待过滤点Pf为中心，设置半窗大小为α，即以待过滤点Pf为中心，在n个原始轨迹点P1至Pn中取Pf前后各α个待过滤点，共获取2α+1个待过滤点Pf-α至Pf+α。确定该2α+1个待过滤点Pf-α至Pf+α中两两待过滤点之间的距离，并将所有距离求和。再确定该2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离，将该距离称为第一距离，则该待过滤点Pf的扭曲度等于上述2α个距离之和与该第一距离的商。
应理解，根据预设扭曲度ζ，若tort(Pf)≥ζ,则认为Pf的扭曲度过大，即该Pf附近的若干个点组成的折线段比较扭曲，置信度较低，需要将Pf过滤；否则可以保留Pf。
可选地，预设扭曲度ζ可以根据实际情况进行设置，一般可以将扭曲度ζ设置为满足ζ≥1。
由于移动终端在路网上做运动时，其轨迹应该为较平滑的线段，也就是对应的扭曲度较小，即趋近于1，如果出现了扭曲严重的情况，则很有可能是采样误差造成的，因此需要执行过滤，通过该扭曲度过滤能过滤掉一部分大误差点。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行切尾均值过滤处理，该切尾均值过滤处理可以更新部分或全部轨迹点的坐标。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：在该n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点中每个待过滤点的经度的平均值，以及h个待过滤点中每个待过滤点的纬度的平均值，k和h为小于2β+1的正整数；将该经度的平均值和该纬度的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第β+1个待过滤点Pl的经度和纬度。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：建立直角坐标系，该直角坐标系包括水平方向的X轴和垂直方向的Y轴；确定该n个原始轨迹点P1至Pn中每个原始轨迹点的坐标；在该n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点中每个待过滤点的X轴坐标的平均值，以及h个待过滤点中每个待过滤点的Y轴坐标的平均值，k和h为小于2β+1的正整数；将该X轴坐标的平均值和该Y轴坐标的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第β+1个待过滤点Pl的X轴坐标和Y轴坐标。
应理解，在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点以及h个待过滤点，该k个待过滤点以及h个待过滤点可以为连续的，也可以为不连续的，k与h的取值可以相等，也可以不相等，该该k个待过滤点以及h个待过滤点可以为相同的点，也可以不相同的点。
可选地，对于在该2β+1个待过滤点Pl-β至Pl+β中确定k个待过滤点，可以将该2β+1个待过滤点中每个待过滤点按照第一坐标值的大小进行排列，该第一坐标可以指上述经度值或者X轴坐标，去掉部分最大和最小值，仅取中间大小的k个待过滤点；类似的，再将每个待过滤点按照第二坐标值的大小进行排列，该第二坐标可以指上述纬度值或者Y轴坐标，去掉部分最大和最小值，仅取中间大小的h个待过滤点。
应理解，可以对该n个原始轨迹点P1至Pn中每个轨迹点都进行该切尾均值过滤处理，或者，也可以对部分轨迹点进行切尾均值过滤处理，例如，仅对部分坐标值变化较大的点进行切尾均值过滤处理。可选地，可以设置差值阈值，在坐标值变化量小于该差值阈值时，确定坐标值变化较小，而不进行切尾均值过滤，但本申请实施例并不限于此。
因此，切尾均值过滤可以通过联系前后若干个点的位置来将当前点的坐标进行一定的调整，使得轨迹的形状更加趋于平滑，对后续地图匹配起到了帮助的作用。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行过近点过滤处理，该过近点过滤处理用于删除距离较近的轨迹点。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：在该n个原始轨迹点P1至Pn中获取第y个待过滤点Py和第y+1个待过滤点Py+1；若该第y个待过滤点Py和该第y+1个待过滤点Py+1之间的距离小于或者等于预设距离，删除第y个待过滤点Py和/或该第y+1个待过滤点Py+1。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，该将该n个原始轨迹点P1至Pn进行过滤处理，包括：在该n个原始轨迹点P1至Pn中获取第y-1个待过滤点Py-1、第y个待过滤点Py和第y+1个待过滤点Py+1；若该第y-1个待过滤点Py-1和该第y个待过滤点Py之间的距离小于或者等于预设距离，且该第y个待过滤点Py和该第y+1个待过滤点Py+1之间的距离也小于或者等于预设距离，则删除第y个待过滤点P y；否则保留该待过滤点Py。
因此，过近点过滤处理可以将由于移动终端的移动速度过慢或静止而造成的定位误差以及对地图匹配的影响给排除掉，解决了移动终端停留导致定位点漂移的问题。例如，该移动终端可以为车辆，车辆行驶过程中可能发生停车或减速，过近点过滤处理可以排除车辆停车或减速过程造成的定位误差。
应理解，由于上述过滤处理可以包括上述差值过滤、速度过滤、角度过滤、扭曲度过滤、切尾均值过滤和过近点过滤中的至少一个，因此，在该过滤处理包括至少两个过滤处理过程时，对于任意一次的非首次过滤处理中的n个原始轨迹点P1至Pn，可以为上一次过滤处理后剩余的轨迹点进行重新编号获得。
应理解，考虑到计算复杂度和精确度，可以任意选择差值过滤、速度过滤、角度过滤、扭曲度过滤、切尾均值过滤和过近点过滤中至少一种过滤方式。例如，对于角度过滤和扭曲度过滤，由于过程类似，可以仅选择其中一种执行。再例如，由于差值过滤较复杂且原始轨迹点完全重叠可能性较小，可以选择不执行差值过滤。再例如，考虑到计算效果，如果确定执行差值过滤，可以将该差值过滤设置为第一次过滤，之后再进行其他过滤过程。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，获取该道路地图，该道路地图中的该多条道路中的每条道路仅包括两个端点，不包括任何岔路。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，将该n个原始轨迹点P1至Pn进行过滤处理，包括：将该n个原始轨迹点P1至Pn进行过滤处理获得m个目标轨迹点，m为小于或者等于n的正整数。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹，包括：将该m个目标轨迹点与该道路地图中的位于道路上的m个点相匹配，该m个点按照时间顺序排列；根据该m个点占用的道路，确定该移动终端在该至少一条道路上的运动轨迹。
结合第一方面及其上述实现方式，在第一方面的另一种实现方式中，根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹，包括：获取该m个点中在时间上连续的第一点和第二点，该第一点的时间早于该第二点的时间；若该第一点和该第二点位于同一条道路，确定该移动终端的轨迹沿该同一条道路由该第一点运动至该第二点；若该第一点位于第一道路，该第二点位于第二道路，该第一道路与该第二道路不同且相连，确定该移动终端的轨迹沿该第一道路上的该第一点运动至该第二道路上的该第二点；若该第一点位于第一道路，该第二点位于第二道路，该第一道路与该第二道路不同且不相连，确定该第一道路的终点到该第二道路的起点之间沿着该道路地图中的道路的最短路径，并确定该移动终端的轨迹沿该第一道路上的该第一点运动至该第一道路的终点，经过该最短路径达到该第二道路的起点，经过该第二道路的起点运动至该第二道路上的该第二点。
因此，本申请实施例的确定运动轨迹的方法，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
第二方面，提供了一种确定运动轨迹的装置，用于执行上述第一方面或第一方面的任意可能的实现方式中的方法。具体地，该装置包括用于执行上述第一方面或第一方面的任意可能的实现方式中的方法的单元。
第三方面，提供了一种确定运动轨迹的装置，包括：存储单元和处理器，该存储单元用于存储指令，该处理器用于执行该存储器存储的指令，并且当该处理器执行该存储器存储的指令时，该执行使得该处理器执行第一方面或第一方面的任意可能的实现方式中的方法。
第四方面，提供了一种计算机可读介质，用于存储计算机程序，该计算机程序包括用于执行第一方面或第一方面的任意可能的实现方式中的方法的指令。
第五方面，提供了一种包括指令的计算机程序产品，当计算机运行该计算机程序产品的该指令时，该计算机执行上述第一方面或第一方面的任意可能的实现方式中的确定运动轨迹的方法。具体地，该计算机程序产品可以运行于上述第三方面的确定运动轨迹的装置上。
附图说明
图1是根据本申请实施例的确定运动轨迹的方法的示意性流程图。
图2是根据本申请实施例的差值过滤的示意图。
图3是根据本申请实施例的速度过滤的示意图。
图4是根据本申请实施例的角度过滤的示意图。
图5是根据本申请实施例的扭曲度过滤的示意图。
图6是根据本申请实施例的切尾均值过滤的示意图。
图7是根据本申请实施例的过近点过滤的示意图。
图8是根据本申请实施例的道路地图匹配方法的示意图。
图9是根据本申请实施例的道路地图匹配方法的另一示意图。
图10是根据本申请实施例的确定运动轨迹的装置的示意性框图。
图11是根据本申请实施例的确定运动轨迹的装置的另一示意性框图。
具体实施方式
下面将结合附图，对本申请中的技术方案进行描述。
应理解，本申请实施例的技术方案可以应用于各种通信系统，例如：全球移动通讯(Global System of Mobile communication，GSM)系统、码分多址(Code DivisionMultiple Access，CDMA)系统、宽带码分多址(Wideband Code Division MultipleAccess，WCDMA)系统、通用分组无线业务(General Packet Radio Service，GPRS)、长期演进(Long Term Evolution，LTE)系统、LTE频分双工(Frequency Division Duplex，FDD)系统、LTE时分双工(Time Division Duplex，TDD)、通用移动通信系统(Universal MobileTelecommunication System，UMTS)、或全球互联微波接入(Worldwide Interoperabilityfor Microwave Access，WiMAX)通信系统等。
在本申请实施例中，移动终端可以包括但不限于移动台(Mobile Station，MS)、移动终端(Mobile Terminal)、移动电话(Mobile Telephone)、用户设备(User Equipment，UE)、手机(handset)及便携设备(portable equipment)、车辆(vehicle)等，该移动终端可以经无线接入网(Radio Access Network，RAN)与一个或多个核心网进行通信，例如，移动终端可以是车辆，还可以是便携式、袖珍式、手持式、计算机内置的或者车载的移动装置等。
本申请实施例所还涉及到基站，该基站可以采集到移动终端在移动过程中的运动轨迹，可选地，该基站可以复用现有的基站，即可以为一种部署在无线接入网中用以为终端设备提供无线通信功能的装置。该基站可以包括各种形式的宏基站，微基站，中继站，接入点等。在采用不同的无线接入技术的系统中，具有基站功能的设备的名称可能会有所不同。例如在LTE网络中，称为演进的节点B(Evolved NodeB，eNB或eNodeB)，在第三代(3rdGeneration，3G)网络中，称为节点B(Node B)等等
图1示出了根据本申请实施例的确定运动轨迹的方法100的示意性流程图，该方法100可以由服务器执行，例如定位服务器，该定位服务器可以获取基站采集到的移动终端的移动轨迹的数据，并将该数据进行相关处理，可选的，该定位服务器还可以为该基站，本申请实施例并不限于此。为了便于描述，下面以该方法100由服务器执行为例进行说明。
如图1所示，该方法100包括：S110，获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该移动终端经过该n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻，n为正整数。
在本申请实施例中，移动终端在移动过程中，可以由基站采集移动终端的运动轨迹，获取若干轨迹点的位置，在S110中，服务器可以获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该n个原始轨迹点P1至Pn可以为基站采集的该若干轨迹点中的全部轨迹点或者部分轨迹点，例如，服务器可以获取该基站采集的全部轨迹点，该n个原始轨迹点P1至Pn即为该全部轨迹点；或者，服务器可以对该全部轨迹点进行处理，该n个原始轨迹点P1至Pn可以为服务器处理原全部轨迹点得到的其中部分轨迹点，本申请实施例并不限于此。
应理解，移动终端的运动过程，随着时间的迁移，移动终端可能位于不同的位置，基站采集移动终端运动过程中的轨迹点，对应每个轨迹点包括移动终端经过该轨迹点的时刻，服务器获取移动终端的n个原始轨迹点P1至Pn可以按照移动终端经过该原始轨迹点的时刻进行排列，例如，该n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻。
应理解，基站采集移动终端在运动过程中轨迹点，可以采用相同的采样时间，例如，采用相同的时间间隔，获取每个时刻的轨迹点的位置，或者也可以采用不同的采样时间，即在不同的时间间隔，采集移动终端移动的轨迹点的位置。
如图1所示，该方法100还包括：S120，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以用于删除该n个原始轨迹点P1至Pn中的至少一个原始轨迹点，和/或更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置。
可选地，作为一个实施例，该过滤处理可以包括差值过滤、速度过滤、角度过滤、扭曲度过滤、切尾均值过滤和过近点过滤中的至少一个，其中，速度过滤、角度过滤、扭曲度过滤和过近点过滤可以用于删除该n个原始轨迹点P1至Pn中的至少一个原始轨迹点，差值过滤和切尾均值过滤可以用于更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置。
可选地，作为一个实施例，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以为差值过滤。具体地，服务器在该n个原始轨迹点P1至Pn中获取a个待过滤点Px至Px+a-1，该a个待过滤点Px至Px+a-1位于第一位置，x和a为正整数，即对于a大于1时，该a个待过滤点Px至Px+a-1在该第一位置上完全重合。再在该n个原始轨迹点P1至Pn中获取b个待过滤点Px+a至Px+a+b-1，该b个待过滤点Px+a至Px+a+b-1位于第二位置，b为大于1的正整数，即该b个待过滤点Px+a至Px+a+b-1在第二位置上完全重合。再在该n个原始轨迹点P1至Pn中获取c个待过滤点Px+a+b至Px+a+b+c-1，该c个待过滤点Px+a+b至Px+a+b+c-1位于第三位置，c为正整数，x+a+b+c-1小于或者等于n，即对于c大于1时，该c个待过滤点Px+a+b至Px+a+b+c-1在第三位置上完全重合。连接第一位置到第二位置的直线，在该第一位置到该第二位置的直线上均匀排列待过滤点至其中，在第一位置上为待过滤点在第二位置上为待过滤点将第一位置与第二位置之间的直线按照距离等分，依次排列待过滤点至之间的每个待过滤点。同样的，连接第二位置到第三位置的直线，在该第二位置到该第三位置的直线上均匀排列待过滤点至其中，在第二位置上为待过滤点在第三位置上为待过滤点将第二位置与第三位置之间的直线按照距离等分，依次排列待过滤点至之间的每个待过滤点，表示向下取整。
例如，图2示出了根据本申请实施例的差值过滤的示意图。如图2所示，取n个原始轨迹点中位于第一位置的a个待过滤点Px至Px+a-1，令a＝3，x＝11，第一位置为A点，即该a个待过滤点为位于A点的P11、P12和P13；取n个原始轨迹点中位于第二位置的b个待过滤点Px+a至Px+a+b-1，令b＝5，第二位置为B点，即该b个待过滤点为位于B点的P14、P15、P16、P17和P18；取n个原始轨迹点中位于第三位置的c个待过滤点Px+a+b至Px+a+b+c-1，令c＝2，第三位置为C点，即该c个待过滤点为位于C点的P19和P20。连接A点与B点得到线段AB，根据a＝3，b＝5可知，在线段AB之间均匀排列P12至P16共5个待过滤点，则把线段AB分为距离相等的4段，得到如图2所示的A点、D点、E点、F点和B点5个点，满足AD＝DE＝EF＝FB，排列结果为：在A点处为P12，在D点处为P'13，E点处为P'14，F点处为P'15，B点处为P16，其中，P'13为P13移动后的位置，P'14为P14移动后的位置，P'15为P15移动后的位置。同样的，连接B点和C点得到线段AB，根据b＝5，c＝2可知，在线段BC之间均匀排列P16至P20共5个待过滤点，则把线段BC分为距离相等的4段，得到如图2所示的B点、G点、H点、K点和C点5个点，满足BG＝GH＝HK＝KC，排列结果为：在B点处为P16，在G点处为P'17，H点处为P'18，K点处为P'19，C点处为P20，其中，P'17为P17移动后的位置，P'18为P18移动后的位置，P'19为P19移动后的位置。
因此，通过差值过滤，可以将重合的点重新排列，更新部分轨迹点的位置，使得轨迹点分布更加均匀。
可选地，作为一个实施例，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以为速度过滤。具体地，服务器获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，该第一待过滤点Pa可以依次取n个原始轨迹点P1至Pn中的每一个轨迹点；若该第一待过滤点Pa的速度大于或者等于预设速度，删除该第一待过滤点Pa。其中，确定第一待过滤点Pa的速度包括：服务器可以在该n个原始轨迹点P1至Pn中获取该第一待过滤点Pa和第二待过滤点Pb，b＝a+1或者b＝a-1；确定从该第一待过滤点Pa到该第二待过滤点Pb的平均速度为该第一待过滤点Pa的速度，但本申请实施例并不限于此。
例如，图3示出了根据本申请实施例的速度过滤的示意图。如图3所示，第一待过滤点Pa可以依次取n个原始轨迹点P1至Pn中的每一个轨迹点，这里以第一待过滤点为P9为例进行说明。确定P9的速度，可以通过确定从P8到P9的平均速度为P9的速度，或者，确定从P9到P10的平均速度为P9的速度。若该P9的速度大于或者等于预设速度，则如图3所示，可以将P9删除；若该P9的速度小于预设速度，则保留P9。
应理解，预设速度可以根据实际情况进行设置，例如，可以根据移动终端在整个运动过程中的平均速度设置，或者是根据移动终端在某一段时间或距离内的速度进行设置，本申请实施例并不限于此。
因此，通过速度过滤，可以将部分误差速度过大的轨迹点过滤掉，这些轨迹的速度不符合常理，与其它轨迹点的速度相差过大，速度过滤可以使得剩余数据更加合理和准备。
可选地，作为一个实施例，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以为角度过滤。具体地，该角度过滤可以用于过滤到大于或者等于预设角度的角的顶点，例如，服务器可以在该n个原始轨迹点P1至Pn中获取4个待过滤点Pc、Pc+1、Pc+2和Pc+3，其中，c为正整数，c+3小于等于n；确定第一角度∠PcPc+1Pc+2和第二角度∠Pc+1Pc+2Pc+3；若该第一角度小于或者等于第一预设角度且该第二角度小于或者等于第二预设角度，删除该4个待过滤点Pc、Pc+1、Pc+2和Pc+3中的待过滤点Pc+1或者待过滤点Pc+2，其中，该预设角度包括该第一预设角度和该第二预设角度。
例如，图4示出了根据本申请实施例的角度过滤的示意图。如图4所示，取n个原始轨迹点P1至Pn中获取4个待过滤点，以P13、P14、P15和P16为例。根据图4可知，测量∠P13P14P15和∠P14P15P16的大小分别为∠P13P14P15＝60°和∠P14P15P16＝55°，假设第一预设角度和第二预设角度均设置为80°，则∠P13P14P15和∠P14P15P16都小于预设角度，因此可以将P14或P15删除，例如图4可以将P15删除。
应理解，预设角度可以根据实际情况进行设置，其中，第一预设角度和第二预设角度可以设置为相同或者不同，本申请实施例并不限于此。
考虑到移动终端在道路上运动时，其轨迹一般应该为较平滑的线段，如果出现了较小的锐角的情况，则很有可能是采样误差造成的，因此，可以通过角度过滤把误差较大的轨迹点过滤掉。
可选地，作为一个实施例，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以为扭曲度过滤。具体地，服务器在该n个原始轨迹点P1至Pn中，获取2α+1个待过滤点Pf-α至Pf+α，f和α为正整数，f大于α。确定该2α+1个待过滤点Pf-α至Pf+α中任意一个待过滤点Pg与相邻的待过滤点Pg+1之间的距离，g分别取f-α至f+α-1中的每一个，一共可以获得2α个距离；确定该2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离为第一距离；确定该2α个距离之和与该第一距离的商为扭曲度；当该扭曲度大于或等于预设扭曲度时，删除该2α+1个待过滤点Pf-α至Pf+α中第α+1个待过滤点Pf。
具体地，服务器可以在该n个原始轨迹点P1至Pn中遍历每一个点，这里以任意一个轨迹点Pf为例。以该待过滤点Pf为中心，设置半窗大小为α，即以待过滤点Pf为中心，在n个原始轨迹点P1至Pn中取Pf前后各α个待过滤点，共获取2α+1个待过滤点Pf-α至Pf+α。确定该2α+1个待过滤点Pf-α至Pf+α中两两待过滤点之间的距离，并将所有距离求和。再确定该2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离，将该距离称为第一距离，则该待过滤点Pf的扭曲度等于上述2α个距离之和与该第一距离的商，即可以通过下面的公式计算该待过滤点Pf的扭曲度tort(Pf)：
其中，dist(A,B)表示计算AB两点之间的距离。
应理解，根据预设扭曲度ζ，若tort(Pf)≥ζ,则认为Pf的扭曲度过大，即该Pf附近的若干个点组成的折线段比较扭曲，置信度较低，需要将Pf过滤；否则可以保留Pf。
可选地，预设扭曲度ζ可以根据实际情况进行设置，一般可以将扭曲度ζ设置为满足ζ≥1，但本申请实施例并不限于此。
例如，图5示出了根据本申请实施例的扭曲度过滤的示意图。如图5所示，服务器在该n个原始轨迹点P1至Pn中遍历每个轨迹点，这里以待过滤点P24为例进行说明。计算P24的扭曲度，这里将半窗大小设置为α＝2，即共取5个待过滤点，分别为点P22至点P26。根据图5可知，待过滤点P22至P26中每相邻的两个点的距离分别为：P22P23＝2，P23P24＝4.5，P24P25＝4，P25P26＝3；该5个待过滤点中第一个待过滤点P22和最后一个待过滤点P26之间的距离为：P22P26＝10。因此，根据上述的公式可知，待过滤点P24的扭曲度为：
假设将预设扭曲度ζ设置为1.2，则1.35>1.2，因此，如图5所示，需要将待过滤点P24删除。
由于移动终端在路网上做运动时，其轨迹应该为较平滑的线段，也就是对应的扭曲度较小，即趋近于1，如果出现了扭曲严重的情况，则很有可能是采样误差造成的，因此需要执行过滤，通过该扭曲度过滤能过滤掉一部分大误差点。
可选地，作为一个实施例，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以为切尾均值过滤。具体地，服务器在该n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；在该2β+1个待过滤点Pl-β至Pl+β中，每个待过滤点的位置可以通过至少两个坐标值来表示，这里以每个待过滤点的位置通过两个坐标值来表示为例进行说明，该两个坐标分别为第一坐标和第二坐标。确定k个待过滤点中每个待过滤点的第一坐标的平均值，以及h个待过滤点中每个待过滤点的第二坐标的平均值，k和h为小于2β+1的正整数；类似的，对于多于两个坐标值，还可以继续求其它坐标值的平均值。将该第一坐标的平均值和该第二坐标的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第l个待过滤点Pl的第一坐标值和第二坐标值。
应理解，在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点以及h个待过滤点，该k个待过滤点以及h个待过滤点可以为连续的，也可以为不连续的，k与h的取值可以相等，也可以不相等，该该k个待过滤点以及h个待过滤点可以为相同的点，也可以不相同的点。
可选地，对于在该2β+1个待过滤点Pl-β至Pl+β中确定k个待过滤点，可以将该2β+1个待过滤点中每个待过滤点按照第一坐标值的大小进行排列，去掉部分最大和最小值，仅取中间大小的k个待过滤点；类似的，再将每个待过滤点按照第二坐标值的大小进行排列，去掉部分最大和最小值，仅取中间大小的h个待过滤点。
可选地，该2β+1个待过滤点Pl-β至Pl+β中每个待过滤点的至少两个坐标值可以为经度坐标和纬度坐标。具体地，在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点Pm至Pm+k-1中每个待过滤点的经度的平均值，以及h个待过滤点Pm至Pm+h-1中每个待过滤点的纬度的平均值，将该经度的平均值和该纬度的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第l个待过滤点Pl的经度和纬度。
可选地，该2β+1个待过滤点Pl-β至Pl+β中每个待过滤点的至少两个坐标值可以指直角坐标系中的X轴坐标和Y轴坐标，其中，该直角坐标系的原点(0,0)可以设置在任意一点上，即以任意一点作为参考地建立该直角坐标系。具体地，假设将n个原始轨迹点P1至Pn中第一轨迹点P1的位置作为原点建立直角坐标系，以水平方向为X轴，垂直方向为Y轴，对于获取的2β+1个待过滤点Pl-β至Pl+β中每个待过滤点都可以确定对应的X轴和Y轴。确定k个待过滤点Pm至Pm+k-1中每个待过滤点的X轴坐标的平均值，以及h个待过滤点Pm至Pm+h-1中每个待过滤点的Y轴坐标的平均值，将该X轴的平均值和该Y轴坐标的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第l个待过滤点Pl的X轴坐标和Y轴坐标。
可选地，该2β+1个待过滤点Pl-β至Pl+β中每个待过滤点的至少两个坐标值还可以为其它坐标系中的坐标值，在此不再一一赘述。
下面以直角坐标系的坐标值举例说明。例如，图6示出了根据本申请实施例的切尾均值过滤的示意图。如图6所示，服务器获取2β+1个待过滤点Pl-β至Pl+β分别为P15至P19，即β＝2，l＝17，每个待过滤点的坐标值如图6所示，其中，假设该直角坐标系取P14的位置为原点，水平方向为X轴，垂直方向为Y轴。对于X轴坐标，将该5个待过滤点的X轴坐标按照大小顺序排列为(0,2,3,4,6)，这里将切尾值设置为1，即取k＝3，去掉1个最大值和1个最小值，剩下的3个坐标为(2,3,4)；同样的，对于Y轴坐标，将该5个待过滤点的Y轴坐标按照大小顺序排列为(-0.5,0,0.2,0.5,2)，这里同样可以将切尾值设置为1，即取h＝3，去掉1个最大值和1个最小值，剩下的3个坐标为(0,0.2,0.5)。因此，计算获得X轴平均值Xnew和Y轴平均值Ynew为：
因此，将Xnew和Ynew作为点P17的新的X轴和Y轴的坐标值(3,0.23)，即如图6所示的P'17为更新后的位置。
应理解，可以对该n个原始轨迹点P1至Pn中每个轨迹点都进行该切尾均值过滤处理，或者，也可以对部分轨迹点进行切尾均值过滤处理，例如，仅对部分坐标值变化较大的点进行切尾均值过滤处理。可选地，可以设置差值阈值，在坐标值变化量小于该差值阈值时，确定坐标值变化较小，而不进行切尾均值过滤，但本申请实施例并不限于此。
以上面的实施例为例，由于原来的轨迹点P17的坐标为(3,2)，经过切尾均值处理后变为P'17＝(3,0.23)，假设设置差值阈值为0.5。由于X轴坐标值的变化量为0，这里以Y轴坐标值变化为例。根据对该P17采取切尾均值过滤的计算方法获得的Y轴坐标值0.23可知，变化量3-0.23＝2.77，且2.77>0.5，因此，可以确定执行切尾均值过滤，更新P17的坐标值为P'17＝(3,0.23)。假设这里计算获得的Y轴坐标不是0.23，而是1.9，由于1.9与原坐标值2相差为0.1，0.1小于0.5，因此，可以不对该轨迹点P17采取切尾均值过滤，仍然使用P17的原坐标值。
应理解，对于任意待过滤点Pl的至少两个坐标值，可以确定将部分或全部坐标值进行更新，本申请实施例并不限于此。
因此，切尾均值过滤可以通过联系前后若干个点的位置来将当前点的坐标进行一定的调整，使得轨迹的形状更加趋于平滑，对后续地图匹配起到了帮助的作用。
可选地，作为一个实施例，将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理可以为过近点过滤。具体地，服务器在该n个原始轨迹点P1至Pn中获取第y个待过滤点Py和第y+1个待过滤点Py+1；若该第y个待过滤点Py和该第y+1个待过滤点Py+1之间的距离小于或者等于预设距离，删除第y个待过滤点Py和/或该第y+1个待过滤点Py+1；否则保留该两个待过滤点Py和Py+1。
可选地，服务器还可以在该n个原始轨迹点P1至Pn中获取第y-1个待过滤点Py-1、第y个待过滤点Py和第y+1个待过滤点Py+1；若该第y-1个待过滤点Py-1和该第y个待过滤点Py之间的距离小于或者等于预设距离，且该第y个待过滤点Py和该第y+1个待过滤点Py+1之间的距离也小于或者等于预设距离，则删除第y个待过滤点Py；否则保留该待过滤点Py。
应理解，该预设距离可以根据实际情况设置，例如，可以根据移动终端的移动速度和采样时间设置该预设距离，或者，还可以根据n个原始轨迹点P1至Pn中每相邻两点之间的距离的平均值进行设置，本申请实施例并不限于此。
例如，图7示出了根据本申请实施例的过近点过滤的示意图。如图7所示，第y个待过滤点Py可以依次取n个原始轨迹点P1至Pn中的每一个轨迹点，这里以待过滤点P15至P23为例进行说明。分别计算P15至P23中每相邻两个待过滤点之间的距离，假设其中仅P16P17、P17P18、P18P19、P19P20和P20P21的距离小于或者等于预设距离。以P16P17为例，可以将P16和P17两点均删除，则类似的需要将P16至P21全部删除；或者，也可以将P16或P17删除，则删除P16对应的删除P16至P20，删除P17对应的删除P17至P21，即如图7所示，剩下的点为P15、P16、P22和P23；或者，将点P16至P21中两端距离均小于或者等于预设距离的待过滤点删除，即删除P17至P20。
因此，过近点过滤处理可以将由于移动终端的移动速度过慢或静止而造成的定位误差以及对地图匹配的影响给排除掉，解决了移动终端停留导致定位点漂移的问题。例如，该移动终端可以为车辆，车辆行驶过程中可能发生停车或减速，过近点过滤处理可以排除车辆停车或减速过程造成的定位误差。
应理解，由于在S120中的过滤处理可以包括上述差值过滤、速度过滤、角度过滤、扭曲度过滤、切尾均值过滤和过近点过滤中的至少一个，因此，在该过滤处理包括至少两个过滤处理过程时，对于任意一次的非首次过滤处理中的n个原始轨迹点P1至Pn，可以为上一次过滤处理后剩余的轨迹点进行重新编号获得。
例如，以7个原始轨迹点P1至P7为例，假设经过一次处理后删除了点P3和P5，并将P6更新为P'6，则剩余P1、P2、P4、P'6和P7，重新编号获得P1、P2、P3(原P4)、P4(原P'6)和P5(原P7)，即作为下一次处理时原始轨迹点，即下一次的n个原始轨迹点变为5个原始轨迹点P1至P5。
应理解，考虑到计算复杂度和精确度，可以任意选择差值过滤、速度过滤、角度过滤、扭曲度过滤、切尾均值过滤和过近点过滤中至少一种过滤方式。例如，对于角度过滤和扭曲度过滤，由于过程类似，可以仅选择其中一种执行。再例如，由于差值过滤较复杂且原始轨迹点完全重叠可能性较小，可以选择不执行差值过滤。再例如，考虑到计算效果，如果确定执行差值过滤，可以将该差值过滤设置为第一次过滤，之后再进行其他过滤过程。
可选地，考虑到最终的计算效果，可以依次进行速度过滤、角度过滤、切尾均值过滤和过近点过滤。或者，也可以依次进行差值过滤、速度过滤、、扭曲度过滤、切尾均值过滤和过近点过滤。
如图1所示，该方法100还包括：S130，根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹，该道路地图包括多条道路。
在本申请实施例中，服务器获取道路地图，该道路地图中包括多条道路。具体地，本申请实施例中的道路地图可以使用现有的各种地图原数据，例如，可以使用全球最大的开源地图项目OpenStreetMap，简记为OSM。OSM官网下载的地图原文件为XML格式；或者，也可以使用其他地图格式，本申请实施例并不限于此。
在本申请实施例中，道路地图中可以包括多条道路，该道路可以分为许多种类，例如可以包括：高速公路(motorway)、主干线(trunk)、住宅区内路线(residential)以及一些无法归类的(unclassified)等等。其中，如果移动终端为车辆，那么可以按照车辆是否可以通过，将该多条道路分为可以通过和不可以通过的，例如，residential等级的道路一般车辆无法通过。
应理解，为了便于将移动终端的运动轨迹与道路地图相匹配，可以将原始的地图中的道路进行切割，获取包括多条道路的道路地图，例如根据道路与道路直接的交叉点进行切割，使得获得的道路地图中的每条道路为最短道路，即使得每条道路上只有首尾两个节点，道路中间不包括任何岔路。可选地，每条道路可以是双向的，也可以是单向的，本申请实施例并不限于此。
可选地，下面以OSM官网下载的XML格式的地图原文件为例，进行道路切割获取道路地图。在原地图中，可以用node标识地图中的点，该点可以为地图中的任意一个点，每个点可以有对应的经纬度信息；还可以用way标识原地图中的折线，可以为每一段way设置k值，当k值为“highway”时，表示该way对应标识了一条道路。
应理解，按照该地图中的道路级别，可以将道路等级从高到底分为motorway、trunk、primary、secondary、tertiary、unclassified、residential与service，由于最后两个等级(residential与service)是车辆无法通过的，这里假设以移动终端为车辆为例，因此在构建并切割原地图获取道路地图时，可以忽略这两个等级的道路。因此，这里仅将表示车辆可以通过的道路保留，即way标识地图中车辆可以通过的道路。
应理解，在原地图中可以用node标识任意一个点，为了获取道路地图，这里仅保留道路上的节点，即保留原地图中用node标识的道路上的节点，即当且仅当两个以上的way同时经过同一个node的情况下，保留该node，该node为道路网络中的节点，而其余的node都只是用来描述路段的形状，这里可以忽略不计。
在本申请实施例中，对每条way需要进行切割，当其中的node被多条way经过时，当前way需要被该node分割成两条路段，这样不停地分割直至一条路段保证只有首尾两个点是被多条way所经过的。具体地，可以通过下面的算法1获得切割后的道路地图。算法1描述了道路地图构建算法的流程，其中，N与W分别代表从OSM所读到的node与way；V为切割后获得的道路地图的顶点集，E为切割后获得的道路地图的边集；E’为一个临时边集，用来记录当前w被切割成的路段集合。split(E,n)函数的行为指：对边集合E中的每条边，检查是否有道路e经过n，如果存在e经过n，则将e根据n切割成两条边e1和e2。
应理解，服务器获取包括多条道路的道路地图，将经过上述过滤处理得到的目标轨迹点中每个目标轨迹点与道路地图中的道路相匹配，获得对应的至少一条道路。具体地，对n个原始轨迹点P1至Pn进行该过滤处理得到m个目标轨迹点，m小于等于n，将该m个目标轨迹点与道路地图进行道路匹配，获得m条道路，其中，该m条道路中可以存在相同的道路。
可选地，可以通过多种算法将m个目标轨迹点与道路例如进行匹配，例如，考虑准确率和鲁棒性，可以采用基于隐马尔可夫模型(Hidden Markov Model，HMM)进行道路匹配，HMM的地图匹配具有最好的准确率和鲁棒性。具体地，HMM是马尔可夫链的一种，它的状态不能直接观察到，但能通过观测序列观察到，假设节点O1,O2,…,On为观测到的序列，HMM假设每个观测到的序列由对应的隐状态产生，对应有H1,H2,…,Hn为隐状态序列，每个隐状态都与前一个隐状态相关，而不与更前面的隐状态相关。HMM地图匹配算法将每个观测点对应的真实路段作为隐状态，移动终端看作在隐状态中进行转移(即路段的转移)，每一个隐状态生成一个观测到的状态(即实际获取的目标轨迹点位置)。通过求出给定观测状态下的最大概率的隐状态序列来得到轨迹匹配后的路径(即路段序列)。
具体地，对于道路地图匹配过程，要对两个概率进行建模。第一个为放射概率，也叫做测量概率，表示的是从空间测量来看，一条路段匹配该状态的概率。图8示出了根据本申请实施例的道路地图匹配方法的示意图，如图8所示，对于获得的过滤后的m个目标轨迹点Z1至Zm，取任意一个目标轨迹点Zi，对于候选的道路ra，其放射概率记为P(Zi|ra)，该概率描述了目标轨迹点Zi在道路ra上的似然。直观上地有，候选路段离采样点距离越远，则实际在该路段运动的可能性也越小，它的放射概率也就越小。假定定位数据误差服从高斯分布，因此有：
其中，如图8所示，Xi,a为Zi在ra上的投影点。
第二个需要建模的为转移概率。对于任意一个目标轨迹点Zi有一系列候选的匹配道路路段，例如道路ra,rb…。转移概率表示的是从当前状态转移到下一个状态的概率。直观上来看，实际上的行驶路径一般都会尽可能地接近最短路径，因此道路上的距离与最短距离相差越大，则其转移概率也越低。
对于任意一个目标轨迹点Zi，其候选的匹配道路为路段ra；目标轨迹点Zi的下一个目标轨迹点Zi+1，其候选的匹配道路为路段rb，记Xi+1,b为目标轨迹点Zi+1在路段rb上的投影点，则从Zi匹配的ra到Zi+1匹配的rb上的转移概率为：
其中δi＝|dist(Zi,Zi+1)-distG(Xi,a,Xi+1,b)|；dist(Zi,Zi+1)为目标轨迹点Zi和Zi+1之间的实际距离；distG(Xi,a,Xi+1,b)为X i,a和Xi+1,b在道路上的距离，即图8中从X i,a到Xi+1,b的灰色折线的长度。
对于一条待匹配的轨迹，即m个目标轨迹点Z＝Z1→Z2→…→Zm，找出其对应的匹配后的m个道路路段R*＝r1,r2,…,rm应该满足：
R*＝arg maxR P(R|Z)
对P(R|Z)展开，利用马尔可夫链的一阶马尔可夫性质，有：
P(R|Z)∝P(Z|R)P(R)
＝P(Z1,Z2,…,Zm|r1,r2,…,rm)P(r1,r2,…,rm)
＝∏P(Zi|ri)∏P(rj|ri)
其中P(zi|ri)为放射概率，P(rj|ri)为转移概率。使用Viterbi算法即可求出一个状态序列R，使得整个序列的后验概率最大，该序列R即对应了m条与m个目标轨迹点对应的道路。
在本申请实施例中，根据m个目标轨迹点获取m条道路；根据该m条道路确定移动终端在道路地图上的移动轨迹。具体地，对于获得的m条道路，可以通过多种算法确定移动终端在道路地图上的移动轨迹。例如，可以通过下述方法确定，这里以相邻的目标轨迹点Zi和Zi+1为例，确定移动终端对应目标轨迹点Zi和Zi+1的实际行驶路径，其中，Zi匹配后的道路记为ri，Zi+1匹配后的道路记为ri+1；服务器确定的移动终端从Zi移动到Zi+1的实际道路路径标记为Ri。
具体地，若ri＝ri+1，即道路ri和道路ri+1为同一条道路时，Ri＝{ri}，即移动终端一直在道路ri上移动。若ri≠ri+1，且ri和ri+1相邻，即道路ri和道路ri+1不是同一条道路，但道路ri和道路ri+1的首尾相连时，Ri＝{ri,ri+1}，即移动终端的运动轨迹为从道路ri至道路ri+1。若ri≠ri+1，且ri和ri+1也不相邻，即道路ri和道路ri+1不是同一条道路，且道路ri和道路ri+1之间不连接时，确定道路ri的末端，标记为ri,s，确定道路ri+1的首端，标记为ri+1,e，根据ri,s和ri+1,e，确定从ri,s到ri+1,e的最短道路路径，标记为r’，则Ri＝{ri}+r’+{ri+1}，例如，如图9所示，根据ri,s和ri+1,e的位置可知，从ri,s到ri+1,e的最短道路路径即为图中r’，因此，移动终端的运动轨迹则为从道路ri至最短道路路径r’，再到道路ri+1，例如图9中灰色线段所示的路径。
以此类推，将m个目标轨迹点对应的m条道路根据上述方式进行匹配，获得连续的至少一条道路即为移动终端在道路地图上的运动轨迹。
因此，本申请实施例的确定运动轨迹的方法，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
应理解，在本申请的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。
另外，本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
上文中结合图1至9，详细描述了根据本申请实施例的确定运动轨迹的方法，下面将结合图10至图11，描述根据本申请实施例的确定运动轨迹的装置。
如图10所示，根据本申请实施例的确定运动轨迹的装置200包括：获取单元210、处理单元220和确定单元230。
该获取单元210用于：获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该移动终端经过该n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻，n为正整数。
该处理单元220用于：将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理用于删除该n个原始轨迹点P1至Pn中的至少一个原始轨迹点，和/或更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置。
该确定单元230用于：根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹，该道路地图包括多条道路。
因此，本申请实施例的确定运动轨迹的装置，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
可选地，该处理单元220具体用于：在该n个原始轨迹点P1至Pn中获取a个待过滤点Px至Px+a-1，该a个待过滤点Px至Px+a-1位于第一位置，x和a为正整数；在该n个原始轨迹点P1至Pn中获取b个待过滤点Px+a至Px+a+b-1，该b个待过滤点Px+a至Px+a+b-1位于第二位置，b为大于1的正整数；在该n个原始轨迹点P1至Pn中获取c个待过滤点Px+a+b至Px+a+b+c-1，该c个待过滤点Px+a+b至Px+a+b+c-1位于第三位置，c为正整数，x+a+b+c-1小于或者等于n；在该第一位置到该第二位置的直线上均匀排列待过滤点至在该第二位置到该第三位置的直线上均匀排列待过滤点至 表示向下取整。
可选地，该处理单元220具体用于：获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，a为小于或者等于n的正整数；若该第一待过滤点Pa的速度大于或者等于预设速度，删除该第一待过滤点Pa。
可选地，该处理单元220具体用于：在该n个原始轨迹点P1至Pn中获取该第一待过滤点Pa和第二待过滤点Pb，b＝a+1或者b＝a-1，a小于n；确定从该第一待过滤点Pa到该第二待过滤点Pb的平均速度为该第一待过滤点Pa的速度。
可选地，该处理单元220具体用于：在该n个原始轨迹点P1至Pn中获取4个待过滤点Pc、Pc+1、Pc+2和Pc+3，c为正整数，c+3小于或者等于n；确定第一角度∠PcPc+1Pc+2和第二角度∠Pc+1Pc+2Pc+3；若该第一角度小于或者等于第一预设角度且该第二角度小于或者等于第二预设角度，删除该4个待过滤点Pc、Pc+1、Pc+2和Pc+3中的待过滤点Pc+1或者待过滤点Pc+2。
可选地，该处理单元220具体用于：在该n个原始轨迹点P1至Pn中，获取2α+1个待过滤点Pf-α至Pf+α，f和α为正整数，f大于α；确定该2α+1个待过滤点Pf-α至Pf+α中待过滤点Pg与相邻的待过滤点Pg+1之间的距离，g取f-α至f+α-1，获得2α个距离；确定该2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离为第一距离；确定该2α个距离之和与该第一距离的商为扭曲度；当该扭曲度大于或等于预设扭曲度时，删除该2α+1个待过滤点Pf-α至Pf+α中第α+1个待过滤点Pf。
可选地，该处理单元220具体用于：在该n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点中每个待过滤点的经度的平均值，以及h个待过滤点中每个待过滤点的纬度的平均值，k和h为小于2β+1的正整数；将该经度的平均值和该纬度的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第β+1个待过滤点Pl的经度和纬度。
可选地，该处理单元220具体用于：在该n个原始轨迹点P1至Pn中获取第y个待过滤点Py和第y+1个待过滤点Py+1；若该第y个待过滤点Py和该第y+1个待过滤点Py+1之间的距离小于或者等于预设距离，删除第y个待过滤点Py和/或该第y+1个待过滤点Py+1。
应理解，根据本申请实施例的确定运动轨迹的装置200可对应于执行本申请实施例中的方法100，并且该装置200中的各个单元的上述和其它操作和/或功能分别为了实现图1至图9中的各个方法中服务器的相应流程，为了简洁，在此不再赘述。
因此，本申请实施例的确定运动轨迹的装置，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
图11示出了根据本申请实施例的确定运动轨迹的装置300的示意性框图，如图11所示，该装置300包括：处理器310，可选地，该装置300还包括存储器320，存储器320与处理器310相连。其中，处理器310和存储器320之间可以通过内部连接通路互相通信，传递和/或控制数据信号，该存储器320可以用于存储指令，该处理器310用于执行该存储器320存储的指令，该处理器310用于：获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，该移动终端经过该n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻，n为正整数；将该n个原始轨迹点P1至Pn进行过滤处理，该过滤处理用于删除该n个原始轨迹点P1至Pn中的至少一个原始轨迹点，和/或更新该n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置；根据对该n个原始轨迹点P1至Pn进行该过滤处理得到的目标轨迹点，确定该移动终端在道路地图中的至少一条道路上的运动轨迹，该道路地图包括多条道路。
因此，本申请实施例的确定运动轨迹的装置，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
可选地，作为一个实施例，该处理器310用于：在该n个原始轨迹点P1至Pn中获取a个待过滤点Px至Px+a-1，该a个待过滤点Px至Px+a-1位于第一位置，x和a为正整数；在该n个原始轨迹点P1至Pn中获取b个待过滤点Px+a至Px+a+b-1，该b个待过滤点Px+a至Px+a+b-1位于第二位置，b为大于1的正整数；在该n个原始轨迹点P1至Pn中获取c个待过滤点Px+a+b至Px+a+b+c-1，该c个待过滤点Px+a+b至Px+a+b+c-1位于第三位置，c为正整数，x+a+b+c-1小于或者等于n；在该第一位置到该第二位置的直线上均匀排列待过滤点至在该第二位置到该第三位置的直线上均匀排列待过滤点至 表示向下取整。
可选地，作为一个实施例，该处理器310用于：获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，a为小于或者等于n的正整数；若该第一待过滤点Pa的速度大于或者等于预设速度，删除该第一待过滤点Pa。
可选地，作为一个实施例，该处理器310用于：在该n个原始轨迹点P1至Pn中获取该第一待过滤点Pa和第二待过滤点Pb，b＝a+1或者b＝a-1，a小于n；确定从该第一待过滤点Pa到该第二待过滤点Pb的平均速度为该第一待过滤点Pa的速度。
可选地，作为一个实施例，该处理器310用于：在该n个原始轨迹点P1至Pn中获取4个待过滤点Pc、Pc+1、Pc+2和Pc+3，c为正整数，c+3小于或者等于n；确定第一角度∠PcPc+1Pc+2和第二角度∠Pc+1Pc+2Pc+3；若该第一角度小于或者等于第一预设角度且该第二角度小于或者等于第二预设角度，删除该4个待过滤点Pc、Pc+1、Pc+2和Pc+3中的待过滤点Pc+1或者待过滤点Pc+2。
可选地，作为一个实施例，该处理器310用于：在该n个原始轨迹点P1至Pn中，获取2α+1个待过滤点Pf-α至Pf+α，f和α为正整数，f大于α；确定该2α+1个待过滤点Pf-α至Pf+α中待过滤点Pg与相邻的待过滤点Pg+1之间的距离，g取f-α至f+α-1，获得2α个距离；确定该2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离为第一距离；确定该2α个距离之和与该第一距离的商为扭曲度；当该扭曲度大于或等于预设扭曲度时，删除该2α+1个待过滤点Pf-α至Pf+α中第α+1个待过滤点Pf。
可选地，作为一个实施例，该处理器310用于：在该n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；在该2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点中每个待过滤点的经度的平均值，以及h个待过滤点中每个待过滤点的纬度的平均值，k和h为小于2β+1的正整数；将该经度的平均值和该纬度的平均值分别确定为该2β+1个待过滤点Pl-β至Pl+β中第β+1个待过滤点Pl的经度和纬度。
可选地，作为一个实施例，该处理器310用于：在该n个原始轨迹点P1至Pn中获取第y个待过滤点Py和第y+1个待过滤点Py+1；若该第y个待过滤点Py和该第y+1个待过滤点Py+1之间的距离小于或者等于预设距离，删除第y个待过滤点Py和/或该第y+1个待过滤点Py+1。
应理解，根据本申请实施例的确定运动轨迹的装置300可对应于本申请实施例中的装置200，并可以对应于执行根据本申请实施例的方法100中的相应主体，并且装置300中的各个单元的上述和其它操作和/或功能分别为了实现图1至图9中的各个方法中服务器的相应流程，为了简洁，在此不再赘述。
因此，本申请实施例的确定运动轨迹的装置，通过对获取移动终端在运动过程中的n个原始轨迹点进行至少一种过滤处理获得m个目标轨迹点，再对该m个目标轨迹点进行道路地图匹配，进而获得移动终端的运动轨迹，能够过滤或改变原始的轨迹点中误差较大以及起到了负面作用的轨迹点，使得经过处理的目标轨迹点在进行地图匹配时更加准确，尤其针对基站定位这样的高噪声的数据获取方式，以及快速定位这种特殊的情况，使得后续地图匹配过程更加准确，获得的运动轨迹更加符合实际情况。
应注意，本申请上述方法实施例可以应用于处理器中，或者由处理器实现。处理器可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法实施例的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(FieldProgrammable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
可以理解，本申请实施例中的存储器可以是易失性存储器或非易失性存储器，或可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(Read-Only Memory，ROM)、可编程只读存储器(Programmable ROM，PROM)、可擦除可编程只读存储器(Erasable PROM，EPROM)、电可擦除可编程只读存储器(Electrically EPROM，EEPROM)或闪存。易失性存储器可以是随机存取存储器(Random Access Memory，RAM)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(Static RAM，SRAM)、动态随机存取存储器(Dynamic RAM，DRAM)、同步动态随机存取存储器(Synchronous DRAM，SDRAM)、双倍数据速率同步动态随机存取存储器(Double Data RateSDRAM，DDR SDRAM)、增强型同步动态随机存取存储器(Enhanced SDRAM，ESDRAM)、同步连接动态随机存取存储器(Synchlink DRAM，SLDRAM)和直接内存总线随机存取存储器(DirectRambus RAM，DR RAM)。应注意，本文描述的系统和方法的存储器旨在包括但不限于这些和任意其它适合类型的存储器。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。

Claims (17)

1.一种确定运动轨迹的方法，其特征在于，包括：
获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，所述移动终端经过所述n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻，n为正整数；
将所述n个原始轨迹点P1至Pn进行过滤处理，所述过滤处理用于删除所述n个原始轨迹点P1至Pn中的至少一个原始轨迹点，和/或更新所述n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置；
根据对所述n个原始轨迹点P1至Pn进行所述过滤处理得到的目标轨迹点，确定所述移动终端在道路地图中的至少一条道路上的运动轨迹，所述道路地图包括多条道路。
2.根据权利要求1所述的方法，其特征在于，所述将所述n个原始轨迹点P1至Pn进行过滤处理，包括：
在所述n个原始轨迹点P1至Pn中获取a个待过滤点Px至Px+a-1，所述a个待过滤点Px至Px+a-1位于第一位置，x和a为正整数；
在所述n个原始轨迹点P1至Pn中获取b个待过滤点Px+a至Px+a+b-1，所述b个待过滤点Px+a至Px+a+b-1位于第二位置，b为大于1的正整数；
在所述n个原始轨迹点P1至Pn中获取c个待过滤点Px+a+b至Px+a+b+c-1，所述c个待过滤点Px+a+b至Px+a+b+c-1位于第三位置，c为正整数，x+a+b+c-1小于或者等于n；
在所述第一位置到所述第二位置的直线上均匀排列待过滤点至
在所述第二位置到所述第三位置的直线上均匀排列待过滤点至 表示向下取整。
3.根据权利要求1所述的方法，其特征在于，所述将所述n个原始轨迹点P1至Pn进行过滤处理，包括：
获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，a为小于或者等于n的正整数；
若所述第一待过滤点Pa的速度大于或者等于预设速度，删除所述第一待过滤点Pa。
4.根据权利要求3所述的方法，其特征在于，所述获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，包括：
在所述n个原始轨迹点P1至Pn中获取所述第一待过滤点Pa和第二待过滤点Pb，b＝a+1或者b＝a-1，a小于n；
确定从所述第一待过滤点Pa到所述第二待过滤点Pb的平均速度为所述第一待过滤点Pa的速度。
5.根据权利要求1所述的方法，其特征在于，所述将所述n个原始轨迹点P1至Pn进行过滤处理，包括：
在所述n个原始轨迹点P1至Pn中获取4个待过滤点Pc、Pc+1、Pc+2和Pc+3，c为正整数，c+3小于或者等于n；
确定第一角度∠PcPc+1Pc+2和第二角度∠Pc+1Pc+2Pc+3；
若所述第一角度小于或者等于第一预设角度且所述第二角度小于或者等于第二预设角度，删除所述4个待过滤点Pc、Pc+1、Pc+2和Pc+3中的待过滤点Pc+1或者待过滤点Pc+2。
6.根据权利要求1所述的方法，其特征在于，所述将所述n个原始轨迹点P1至Pn进行过滤处理，包括：
在所述n个原始轨迹点P1至Pn中，获取2α+1个待过滤点Pf-α至Pf+α，f和α为正整数，f大于α；
确定所述2α+1个待过滤点Pf-α至Pf+α中待过滤点Pg与相邻的待过滤点Pg+1之间的距离，g取f-α至f+α-1，获得2α个距离；
确定所述2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离为第一距离；
确定所述2α个距离之和与所述第一距离的商为扭曲度；
当所述扭曲度大于或等于预设扭曲度时，删除所述2α+1个待过滤点Pf-α至Pf+α中第α+1个待过滤点Pf。
7.根据权利要求1所述的方法，其特征在于，所述将所述n个原始轨迹点P1至Pn进行过滤处理，包括：
在所述n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；
在所述2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点中每个待过滤点的经度的平均值，以及h个待过滤点中每个待过滤点的纬度的平均值，k和h为小于2β+1的正整数；
将所述经度的平均值和所述纬度的平均值分别确定为所述2β+1个待过滤点Pl-β至Pl+β中第β+1个待过滤点Pl的经度和纬度。
8.根据权利要求1所述的方法，其特征在于，所述将所述n个原始轨迹点P1至Pn进行过滤处理，包括：
在所述n个原始轨迹点P1至Pn中获取第y个待过滤点Py和第y+1个待过滤点Py+1；
若所述第y个待过滤点Py和所述第y+1个待过滤点Py+1之间的距离小于或者等于预设距离，删除第y个待过滤点Py和/或所述第y+1个待过滤点Py+1。
9.一种确定运动轨迹的装置，其特征在于，包括：
获取单元，用于获取移动终端在运动过程中的n个原始轨迹点P1至Pn的位置，其中，所述移动终端经过所述n个原始轨迹点P1至Pn中第i个原始轨迹点Pi的时刻早于经过第i+1个原始轨迹点Pi+1的时刻，n为正整数；
处理单元，用于将所述n个原始轨迹点P1至Pn进行过滤处理，所述过滤处理用于删除所述n个原始轨迹点P1至Pn中的至少一个原始轨迹点，和/或更新所述n个原始轨迹点P1至Pn中的至少一个原始轨迹点的位置；
确定单元，用于根据对所述n个原始轨迹点P1至Pn进行所述过滤处理得到的目标轨迹点，确定所述移动终端在道路地图中的至少一条道路上的运动轨迹，所述道路地图包括多条道路。
10.根据权利要求9所述的装置，其特征在于，所述处理单元具体用于：
在所述n个原始轨迹点P1至Pn中获取a个待过滤点Px至Px+a-1，所述a个待过滤点Px至Px+a-1位于第一位置，x和a为正整数；
在所述n个原始轨迹点P1至Pn中获取b个待过滤点Px+a至Px+a+b-1，所述b个待过滤点Px+a至Px+a+b-1位于第二位置，b为大于1的正整数；
在所述n个原始轨迹点P1至Pn中获取c个待过滤点Px+a+b至Px+a+b+c-1，所述c个待过滤点Px+a+b至Px+a+b+c-1位于第三位置，c为正整数，x+a+b+c-1小于或者等于n；
在所述第一位置到所述第二位置的直线上均匀排列待过滤点至
在所述第二位置到所述第三位置的直线上均匀排列待过滤点至 表示向下取整。
11.根据权利要求9所述的装置，其特征在于，所述处理单元具体用于：
获取n个原始轨迹点P1至Pn中的第一待过滤点Pa的速度，a为小于或者等于n的正整数；
若所述第一待过滤点Pa的速度大于或者等于预设速度，删除所述第一待过滤点Pa。
12.根据权利要求11所述的装置，其特征在于，所述处理单元具体用于：
在所述n个原始轨迹点P1至Pn中获取所述第一待过滤点Pa和第二待过滤点Pb，b＝a+1或者b＝a-1，a小于n；
确定从所述第一待过滤点Pa到所述第二待过滤点Pb的平均速度为所述第一待过滤点Pa的速度。
13.根据权利要求9所述的装置，其特征在于，所述处理单元具体用于：
在所述n个原始轨迹点P1至Pn中获取4个待过滤点Pc、Pc+1、Pc+2和Pc+3，c为正整数，c+3小于或者等于n；
确定第一角度∠PcPc+1Pc+2和第二角度∠Pc+1Pc+2Pc+3；
若所述第一角度小于或者等于第一预设角度且所述第二角度小于或者等于第二预设角度，删除所述4个待过滤点Pc、Pc+1、Pc+2和Pc+3中的待过滤点Pc+1或者待过滤点Pc+2。
14.根据权利要求9所述的装置，其特征在于，所述处理单元具体用于：
在所述n个原始轨迹点P1至Pn中，获取2α+1个待过滤点Pf-α至Pf+α，f和α为正整数，f大于α；
确定所述2α+1个待过滤点Pf-α至Pf+α中待过滤点Pg与相邻的待过滤点Pg+1之间的距离，g取f-α至f+α-1，获得2α个距离；
确定所述2α+1个待过滤点Pf-α至Pf+α中第一个待过滤点Pf-α与最后一个待过滤点Pf+α之间的距离为第一距离；
确定所述2α个距离之和与所述第一距离的商为扭曲度；
当所述扭曲度大于或等于预设扭曲度时，删除所述2α+1个待过滤点Pf-α至Pf+α中第α+1个待过滤点Pf。
15.根据权利要求9所述的装置，其特征在于，所述处理单元具体用于：
在所述n个原始轨迹点P1至Pn中，获取2β+1个待过滤点Pl-β至Pl+β，l和β为正整数，l大于β；
在所述2β+1个待过滤点Pl-β至Pl+β中，确定k个待过滤点中每个待过滤点的经度的平均值，以及h个待过滤点中每个待过滤点的纬度的平均值，k和h为小于2β+1的正整数；
将所述经度的平均值和所述纬度的平均值分别确定为所述2β+1个待过滤点Pl-β至Pl+β中第β+1个待过滤点Pl的经度和纬度。
16.根据权利要求9所述的装置，其特征在于，所述处理单元具体用于：
在所述n个原始轨迹点P1至Pn中获取第y个待过滤点Py和第y+1个待过滤点Py+1；
若所述第y个待过滤点Py和所述第y+1个待过滤点Py+1之间的距离小于或者等于预设距离，删除第y个待过滤点Py和/或所述第y+1个待过滤点Py+1。
17.一种计算机可读介质，用于存储计算机程序，该计算机程序包括用于执行权利要求1至8中任一项所述的方法的指令。
CN201710586941.1A 2017-07-18 2017-07-18 确定运动轨迹的方法和装置 Pending CN109269514A (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201710586941.1A CN109269514A (zh) 2017-07-18 2017-07-18 确定运动轨迹的方法和装置
PCT/CN2018/094740 WO2019015485A1 (zh) 2017-07-18 2018-07-06 确定运动轨迹的方法和装置
EP18834716.5A EP3640596A4 (en) 2017-07-18 2018-07-06 METHOD AND DEVICE FOR DETERMINING A TRAIL OF TRAVEL
US16/746,666 US20200151885A1 (en) 2017-07-18 2020-01-17 Method and apparatus for determining moving track

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710586941.1A CN109269514A (zh) 2017-07-18 2017-07-18 确定运动轨迹的方法和装置

Publications (1)

Publication Number Publication Date
CN109269514A true CN109269514A (zh) 2019-01-25

Family

ID=65015376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710586941.1A Pending CN109269514A (zh) 2017-07-18 2017-07-18 确定运动轨迹的方法和装置

Country Status (4)

Country Link
US (1) US20200151885A1 (zh)
EP (1) EP3640596A4 (zh)
CN (1) CN109269514A (zh)
WO (1) WO2019015485A1 (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110650437A (zh) * 2019-11-27 2020-01-03 苏宁云计算有限公司 电子围栏上的围栏点删除方法、装置及计算机设备
CN111158368A (zh) * 2019-12-31 2020-05-15 深圳市优必选科技股份有限公司 一种双足机器人及其轨迹跟随方法和装置
CN111368014A (zh) * 2019-12-23 2020-07-03 广东小天才科技有限公司 一种运动轨迹的生成方法、终端设备及存储介质
CN111679297A (zh) * 2020-05-08 2020-09-18 四川超影科技有限公司 一种gps定位轨迹的噪声点漂移去除方法
CN111723304A (zh) * 2020-01-03 2020-09-29 腾讯科技（深圳）有限公司 一种轨迹点识别方法和相关装置
CN111735461A (zh) * 2020-06-10 2020-10-02 腾讯科技（深圳）有限公司 行驶轨迹的处理方法、装置及电子设备
CN111857113A (zh) * 2019-04-12 2020-10-30 北京地平线机器人技术研发有限公司 可移动设备的定位方法及定位装置
CN112985430A (zh) * 2019-12-13 2021-06-18 百度在线网络技术（北京）有限公司 地图匹配的容灾方法、装置、设备及存储介质
CN113408973A (zh) * 2020-03-17 2021-09-17 北京京东振世信息技术有限公司 生成轨迹数据的方法和装置
CN113495939A (zh) * 2020-04-08 2021-10-12 中移智行网络科技有限公司 基于道格拉斯-普克算法的车辆轨迹纠偏方法和装置
CN113539050A (zh) * 2020-04-20 2021-10-22 华为技术有限公司 数据处理方法、装置及设备
CN115454144A (zh) * 2022-10-28 2022-12-09 中国电子科技集团公司第二十八研究所 一种动目标飞行轨迹平滑方法及系统
CN116430423A (zh) * 2023-06-13 2023-07-14 广州悦跑信息科技有限公司 一种运动数据中卫星导航定位轨迹点坐标方法
CN118033703A (zh) * 2024-04-12 2024-05-14 华能信息技术有限公司 一种基于地图融合的无人驾驶矿车定位方法及系统

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110858990B (zh) * 2018-08-22 2021-04-20 华为技术有限公司 一种网络重选的方法及装置
US11218887B2 (en) * 2019-10-15 2022-01-04 Cisco Technology, Inc. Switching of networks for a mobile device using location based predictive algorithm
CN110991525B (zh) * 2019-11-29 2022-08-16 西安交通大学 基于运营商轨迹数据的伴随模式匹配方法
CN111615061B (zh) * 2020-05-09 2022-02-15 国家计算机网络与信息安全管理中心山东分中心 移动终端轨迹数据的去噪方法及装置
CN112487114B (zh) * 2020-11-10 2021-08-31 电子科技大学 一种城市交叉口地图匹配方法
CN112365088B (zh) * 2020-11-28 2024-04-26 北京梧桐车联科技有限责任公司 行程关键点的确定方法、装置、设备及可读存储介质
CN112527932B (zh) * 2020-12-04 2023-09-26 北京百度网讯科技有限公司 道路数据处理的方法、装置、设备及存储介质
CN112541852B (zh) * 2020-12-24 2024-04-12 南方科技大学 城市人流监控方法、装置、电子设备及存储介质
CN112907686B (zh) * 2021-02-09 2021-12-14 青海师范大学 用于矢量轨迹压缩的余弦垂距判别方法、装置和设备
US20220341742A1 (en) * 2021-04-27 2022-10-27 Via Transportation, Inc. Systems and methods for route reconstruction
CN113157848A (zh) * 2021-05-06 2021-07-23 清华大学 航路确定方法及装置、电子设备和存储介质
CN113806426A (zh) * 2021-11-16 2021-12-17 亿海蓝（北京）数据技术股份公司 轨迹数据的处理方法、装置、可读存储介质和电子设备
CN114838735B (zh) * 2022-03-21 2024-05-10 江苏盛海智能科技有限公司 一种基于movebase的路径循迹方法及终端
KR20240015458A (ko) * 2022-07-27 2024-02-05 홍익대학교 산학협력단 지구위치데이터 맵 매칭방법 및 이를 위한 컴퓨팅 장치

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101458434A (zh) * 2009-01-08 2009-06-17 浙江大学 精确测量和预测乒乓球轨迹系统与系统运行方法
US20090198391A1 (en) * 2008-02-05 2009-08-06 Ajith Kuttannair Kumar System, method and computer software code for obtaining information for routing a powered system and adjusting a route in accordance with relevant information
CN102298152A (zh) * 2010-06-23 2011-12-28 爱信艾达株式会社 轨迹信息生成装置和方法以及计算机可读存储介质
CN105143909A (zh) * 2012-06-26 2015-12-09 多伦多大学理事会 用于无线电地图的动态生成的系统、方法和计算机程序
CN106650771A (zh) * 2016-09-29 2017-05-10 百度在线网络技术（北京）有限公司 基于聚类分析的轨迹去噪方法以及装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4348441B2 (ja) * 2007-01-22 2009-10-21 国立大学法人 大阪教育大学 位置検出装置、位置検出方法、データ判定装置、データ判定方法、コンピュータプログラム及び記憶媒体
CN103037507A (zh) * 2012-12-17 2013-04-10 浙江鸿程计算机系统有限公司 一种基于Cell-ID定位技术的地图匹配方法
CN104115180B (zh) * 2013-02-21 2017-06-09 索尼公司 信息处理设备、信息处理方法和程序
CN104596507B (zh) * 2015-02-09 2017-10-03 成都小步创想畅联科技有限公司 一种移动终端出行轨迹的确定方法
CN105628033B (zh) * 2016-02-26 2019-04-02 广西鑫朗通信技术有限公司 一种基于道路连通关系的地图匹配方法
CN105890599A (zh) * 2016-04-12 2016-08-24 上海牵挂网络科技有限公司 一种基于地理位置足迹噪点的过滤算法
CN106231671B (zh) * 2016-08-09 2019-05-17 南京掌控网络科技有限公司 一种移动设备的移动轨迹优化方法
CN106483537A (zh) * 2016-09-22 2017-03-08 深圳市时代经纬科技有限公司 一种卫星定位轨迹的优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090198391A1 (en) * 2008-02-05 2009-08-06 Ajith Kuttannair Kumar System, method and computer software code for obtaining information for routing a powered system and adjusting a route in accordance with relevant information
CN101458434A (zh) * 2009-01-08 2009-06-17 浙江大学 精确测量和预测乒乓球轨迹系统与系统运行方法
CN102298152A (zh) * 2010-06-23 2011-12-28 爱信艾达株式会社 轨迹信息生成装置和方法以及计算机可读存储介质
CN105143909A (zh) * 2012-06-26 2015-12-09 多伦多大学理事会 用于无线电地图的动态生成的系统、方法和计算机程序
CN106650771A (zh) * 2016-09-29 2017-05-10 百度在线网络技术（北京）有限公司 基于聚类分析的轨迹去噪方法以及装置

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111857113B (zh) * 2019-04-12 2023-11-03 北京地平线机器人技术研发有限公司 可移动设备的定位方法及定位装置
CN111857113A (zh) * 2019-04-12 2020-10-30 北京地平线机器人技术研发有限公司 可移动设备的定位方法及定位装置
CN110650437B (zh) * 2019-11-27 2020-02-18 苏宁云计算有限公司 电子围栏上的围栏点删除方法、装置及计算机设备
CN110650437A (zh) * 2019-11-27 2020-01-03 苏宁云计算有限公司 电子围栏上的围栏点删除方法、装置及计算机设备
CN112985430B (zh) * 2019-12-13 2024-06-25 百度在线网络技术（北京）有限公司 地图匹配的容灾方法、装置、设备及存储介质
CN112985430A (zh) * 2019-12-13 2021-06-18 百度在线网络技术（北京）有限公司 地图匹配的容灾方法、装置、设备及存储介质
CN111368014A (zh) * 2019-12-23 2020-07-03 广东小天才科技有限公司 一种运动轨迹的生成方法、终端设备及存储介质
CN111368014B (zh) * 2019-12-23 2024-04-19 广东小天才科技有限公司 一种运动轨迹的生成方法、终端设备及存储介质
CN111158368A (zh) * 2019-12-31 2020-05-15 深圳市优必选科技股份有限公司 一种双足机器人及其轨迹跟随方法和装置
CN111158368B (zh) * 2019-12-31 2024-02-02 深圳市优必选科技股份有限公司 一种双足机器人及其轨迹跟随方法和装置
CN111723304A (zh) * 2020-01-03 2020-09-29 腾讯科技（深圳）有限公司 一种轨迹点识别方法和相关装置
CN113408973A (zh) * 2020-03-17 2021-09-17 北京京东振世信息技术有限公司 生成轨迹数据的方法和装置
CN113495939A (zh) * 2020-04-08 2021-10-12 中移智行网络科技有限公司 基于道格拉斯-普克算法的车辆轨迹纠偏方法和装置
CN113539050B (zh) * 2020-04-20 2022-09-23 华为技术有限公司 数据处理方法、装置及设备
CN113539050A (zh) * 2020-04-20 2021-10-22 华为技术有限公司 数据处理方法、装置及设备
CN111679297A (zh) * 2020-05-08 2020-09-18 四川超影科技有限公司 一种gps定位轨迹的噪声点漂移去除方法
CN111735461B (zh) * 2020-06-10 2023-11-17 腾讯科技（深圳）有限公司 行驶轨迹的处理方法、装置及电子设备
CN111735461A (zh) * 2020-06-10 2020-10-02 腾讯科技（深圳）有限公司 行驶轨迹的处理方法、装置及电子设备
CN115454144A (zh) * 2022-10-28 2022-12-09 中国电子科技集团公司第二十八研究所 一种动目标飞行轨迹平滑方法及系统
CN115454144B (zh) * 2022-10-28 2023-02-24 中国电子科技集团公司第二十八研究所 一种动目标飞行轨迹平滑方法及系统
CN116430423A (zh) * 2023-06-13 2023-07-14 广州悦跑信息科技有限公司 一种运动数据中卫星导航定位轨迹点坐标方法
CN116430423B (zh) * 2023-06-13 2023-08-29 广州悦跑信息科技有限公司 一种运动数据中卫星导航定位轨迹点坐标修正方法
CN118033703A (zh) * 2024-04-12 2024-05-14 华能信息技术有限公司 一种基于地图融合的无人驾驶矿车定位方法及系统

Also Published As

Publication number Publication date
EP3640596A4 (en) 2020-09-16
WO2019015485A1 (zh) 2019-01-24
US20200151885A1 (en) 2020-05-14
EP3640596A1 (en) 2020-04-22

Similar Documents

Publication Publication Date Title
CN109269514A (zh) 确定运动轨迹的方法和装置
US10323948B2 (en) GPS data repair
CN106610294B (zh) 定位方法及装置
CN105403228B (zh) 一种运动距离的确定方法和装置
CN107643085A (zh) 一种路径推荐方法与装置
CN105223548B (zh) 用于基于预期路线来向移动设备提供无线发射机历书信息的方法和装置
CN106610981B (zh) 电子地图中道路信息的验证更新方法及系统
CN106339716A (zh) 一种基于加权欧氏距离的移动轨迹相似度匹配方法
CN101836081A (zh) 处理定位数据的方法
JP2005283575A (ja) 移動先予測装置および移動先予測方法
CN108362293A (zh) 一种基于关键点技术的车辆轨迹匹配方法
CN105738932A (zh) 轨迹的生成方法和装置
CN105352520A (zh) 一种导航路线修正方法和装置
CN107330085B (zh) 大数据环境下固定传感器错误位置的判断识别和矫正方法
CN105180951B (zh) 对车辆的路线规划
CN107871400A (zh) 路网信息更新方法及装置
CN111664848A (zh) 一种多模态室内定位导航方法及系统
CN103968850B (zh) 更新兴趣点引导信息的方法及装置
CN102997927A (zh) 信息采集和处理方法及装置
CN105091888A (zh) 导航装置及更新其图资的方法
CN104101357B (zh) 导航系统以及在导航系统上显示影像地图的方法
CN108225349A (zh) 导航中车辆位置信息修正方法和装置
CN102087113B (zh) 一种车载机自主导航方法
CN106887184B (zh) 路线更新方法及装置
CN106525060A (zh) 一种信息更正方法、智能终端、服务器和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190125

RJ01 Rejection of invention patent application after publication