Nothing Special   »   [go: up one dir, main page]

CN111489585A - Vehicle and pedestrian collision avoidance method based on edge calculation - Google Patents

Vehicle and pedestrian collision avoidance method based on edge calculation Download PDF

Info

Publication number
CN111489585A
CN111489585A CN202010143534.5A CN202010143534A CN111489585A CN 111489585 A CN111489585 A CN 111489585A CN 202010143534 A CN202010143534 A CN 202010143534A CN 111489585 A CN111489585 A CN 111489585A
Authority
CN
China
Prior art keywords
vehicle
pedestrian
pedestrians
collision
noise
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
CN202010143534.5A
Other languages
Chinese (zh)
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.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202010143534.5A priority Critical patent/CN111489585A/en
Publication of CN111489585A publication Critical patent/CN111489585A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/16Anti-collision systems
    • G08G1/164Centralised systems, e.g. external to vehicles
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/06Testing, supervising or monitoring using simulated traffic
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/30Services specially adapted for particular environments, situations or purposes
    • H04W4/40Services specially adapted for particular environments, situations or purposes for vehicles, e.g. vehicle-to-pedestrians [V2P]
    • H04W4/44Services specially adapted for particular environments, situations or purposes for vehicles, e.g. vehicle-to-pedestrians [V2P] for communication between vehicles and infrastructures, e.g. vehicle-to-cloud [V2C] or vehicle-to-home [V2H]

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种基于边缘计算的车辆行人碰撞避免方法,属于物联网领域,在边缘计算环境下,通过边缘服务器实现车辆和行人设备的通信,减小通信延迟,采用卡尔曼滤波来预测车辆的位置,利用行人设备上的传感器来预测行人轨迹,克服了视距的限制,提高了碰撞检测精度,适用于大多数非视距(NLOS)场景,并根据车辆与行人之间的相对运动关系,建立了基于碰撞时间(TTC)的碰撞风险模型,仿真结果表明,本发明所提出的碰撞避免算法可以有效降低车辆行人的碰撞风险。

Figure 202010143534

The invention discloses a vehicle-pedestrian collision avoidance method based on edge computing, which belongs to the field of Internet of Things. In the edge computing environment, communication between vehicles and pedestrian equipment is realized through an edge server, communication delay is reduced, and Kalman filtering is used to predict vehicles The position of the pedestrian device is used to predict the pedestrian trajectory, which overcomes the limitation of the line of sight and improves the collision detection accuracy. , a collision risk model based on time to collision (TTC) is established, and the simulation results show that the collision avoidance algorithm proposed in the present invention can effectively reduce the collision risk of vehicles and pedestrians.

Figure 202010143534

Description

一种基于边缘计算的车辆行人碰撞避免方法A vehicle-pedestrian collision avoidance method based on edge computing

技术领域technical field

本发明涉及物联网技术领域,特别涉及一种基于边缘计算的车辆行人碰撞避免方法。The invention relates to the technical field of the Internet of Things, in particular to a vehicle-pedestrian collision avoidance method based on edge computing.

背景技术Background technique

2017年Hassan Artail在Avoiding car-pedestrian collisions using a VANETto cellular communication framework一文中提出了一种利用蜂窝技术实现行人和车辆通信的方案,从而来检测可能发生的车辆行人碰撞。该技术方案中采用了远程云服务器做为中间服务器来间接实现车辆和行人的通信,克服了视距的限制,但是其系统延迟很高并且难以确定,这对于对延迟极其敏感的行人安全算法程序而言无疑是致命的缺陷;并且,该技术方案中的行人和车辆位置预测算法都极为粗糙,大部分情况下会将行人和车辆可能同时出现的位置全部划为非安全区域,误警概率较高。In 2017, Hassan Artail proposed a scheme to use cellular technology to realize pedestrian and vehicle communication in Avoiding car-pedestrian collisions using a VANET to cellular communication framework to detect possible vehicle-pedestrian collisions. In this technical solution, a remote cloud server is used as an intermediate server to indirectly realize the communication between vehicles and pedestrians, which overcomes the limitation of sight distance, but the system delay is very high and difficult to determine, which is extremely sensitive to the delay of pedestrian safety algorithm programs. It is undoubtedly a fatal flaw; in addition, the pedestrian and vehicle position prediction algorithms in this technical solution are extremely rough. In most cases, the positions where pedestrians and vehicles may appear at the same time are all classified as unsafe areas, and the probability of false alarms is relatively high. high.

发明内容SUMMARY OF THE INVENTION

本发明的技术目的是提供一种基于边缘计算的车辆行人碰撞避免方法,该方法不仅降低了车辆与行人之间的通信延迟,还提高了车辆行人碰撞事故的预警准确率。The technical purpose of the present invention is to provide a vehicle-pedestrian collision avoidance method based on edge computing, which not only reduces the communication delay between vehicles and pedestrians, but also improves the early warning accuracy of vehicle-pedestrian collision accidents.

本发明的上述技术目的是通过以下技术方案得以实现的:The above-mentioned technical purpose of the present invention is achieved through the following technical solutions:

一种基于边缘计算的车辆行人碰撞避免方法,包括以下步骤:A vehicle-pedestrian collision avoidance method based on edge computing, comprising the following steps:

Step1:建立行人设备到车辆通信的网络模型:在十字路口或者道路边缘布置路侧单元(RSU)/基站、与所述路侧单元(RSU)/基站连接的移动边缘计算(MEC)服务器,所述行人设备为具备GPS模块的智能手机,所述行人设备和车辆能在所述网络模型覆盖范围内实现通信;Step1: Establish a network model for pedestrian device-to-vehicle communication: Arrange roadside units (RSUs)/base stations at intersections or road edges, and mobile edge computing (MEC) servers connected to the roadside units (RSUs)/base stations, so The pedestrian device is a smart phone with a GPS module, and the pedestrian device and the vehicle can communicate within the coverage of the network model;

在这种模式下MEC服务器可以连接到路边单元(RSU)/基站上,车辆通过LTE-V2X实现和RSU/基站的通信,行人则通过智能手机和RSU/基站进行通信;计算任务可以在车辆或者MEC服务器上进行;这样的系统既利用了传统的车联网中车辆的资源,并且在计算任务繁重时通过在MEC服务器上运行计算任务缓解了车辆上计算资源的压力,同时也利用了智能手机上的位置服务等资源,让智能手机参与到了行人安全保护的应用中;In this mode, the MEC server can be connected to the roadside unit (RSU)/base station, the vehicle communicates with the RSU/base station through LTE-V2X, and the pedestrian communicates with the RSU/base station through a smartphone; computing tasks can be performed in the vehicle Or on the MEC server; such a system not only utilizes the resources of the vehicle in the traditional Internet of Vehicles, but also relieves the pressure of the computing resources on the vehicle by running the computing task on the MEC server when the computing task is heavy, and also utilizes the smart phone. resources such as location services on the Internet, allowing smartphones to participate in the application of pedestrian safety protection;

Step2:根据车辆和行人设备当前信息和历史信息,利用卡尔曼滤波器来对车辆和行人未来位置分别进行预测;Step2: According to the current information and historical information of vehicles and pedestrian equipment, use Kalman filter to predict the future positions of vehicles and pedestrians respectively;

Step3:根据Step2中预测到的车辆和行人路径,根据相互运动关系,运行车辆行人碰撞避免算法,对可能发生的碰撞进行预警。Step3: According to the predicted vehicle and pedestrian paths in Step2, according to the mutual motion relationship, run the vehicle-pedestrian collision avoidance algorithm to give early warning of possible collisions.

进一步地,所述Step2中预测车辆未来位置的具体步骤为:Further, the specific steps of predicting the future position of the vehicle in the Step2 are:

假设汽车具有恒定转率和速度,选取车辆的状态向量为:Assuming that the car has a constant rotation rate and speed, the state vector of the selected vehicle is:

X=[x y v θ ω]T X=[xyv θ ω] T

θ为偏航角,是追踪的目标车辆在当前车辆坐标系下与x轴的夹角,逆时针方向为正,取值范围是[0,2π),ω是偏航角速度;行驶的汽车由于道路和交规等的限制其运动状态一般是直行或者转弯,这俩种运动状态都可以用此种状态向量来表示,当ω≠0时是非直线行驶,此时模型的状态转移函数为:θ is the yaw angle, which is the angle between the tracked target vehicle and the x-axis in the current vehicle coordinate system, the counterclockwise direction is positive, and the value range is [0, 2π), and ω is the yaw angular velocity; The motion states of roads and traffic regulations are generally straight or turning. Both motion states can be represented by this state vector. When ω≠0, it is non-straight travel. At this time, the state transfer function of the model is:

Figure RE-GDA0002536505580000021
Figure RE-GDA0002536505580000021

当ω=0时为直线行驶,此时模型的状态转移函数变为:When ω=0, it is a straight line, and the state transition function of the model becomes:

Figure RE-GDA0002536505580000022
Figure RE-GDA0002536505580000022

对于(1)和(2)的多元函数,使用多元泰勒级数进行展开,忽略高阶级数,只考虑利用雅克比矩阵进行线性化,在此模型中,对各个元素求偏导数可以得到雅可比矩阵JA,当ω≠0时:For the multivariate functions of (1) and (2), the multivariate Taylor series is used for expansion, the higher order number is ignored, and only the Jacobian matrix is used for linearization. In this model, the partial derivative of each element can be obtained. Jacobian Matrix J A , when ω≠0:

Figure RE-GDA0002536505580000023
Figure RE-GDA0002536505580000023

当ω=0时,雅可比矩阵为:When ω=0, the Jacobian matrix is:

Figure RE-GDA0002536505580000024
Figure RE-GDA0002536505580000024

计算处理噪声带来的不确定性Q;在所述网络模型中,噪声的引入主要来源于两处:直线加速度和偏航角加速度噪声,所述直线加速度和偏航角加速度噪声会影响状态量(x,y,v,θ,ω),所述直线加速度和偏航角加速度对状态的影响的表达式如下:Calculate and deal with the uncertainty Q caused by noise; in the network model, the introduction of noise mainly comes from two sources: linear acceleration and yaw angular acceleration noise, the linear acceleration and yaw angular acceleration noise will affect the state quantity (x, y, v, θ, ω), the expression of the influence of the linear acceleration and yaw angular acceleration on the state is as follows:

Figure RE-GDA0002536505580000031
Figure RE-GDA0002536505580000031

其中μa

Figure RE-GDA0002536505580000036
为直线上和转角上的加速度,Q就是处理噪声的协方差矩阵,其表达式为:where μ a ,
Figure RE-GDA0002536505580000036
is the acceleration on the straight line and on the corner, Q is the covariance matrix for dealing with noise, and its expression is:

Q=E[n·nT]=E[GμμTGT]=G·E[μμT]·GT (4)Q=E[n·n T ]=E[GμμT G T ]=G·E[μμ T ]·G T ( 4)

其中in

Figure RE-GDA0002536505580000032
Figure RE-GDA0002536505580000032

所以,在CTRV模型中的处理噪声的协方差矩阵Q的计算公式就是:Therefore, the calculation formula of the covariance matrix Q of processing noise in the CTRV model is:

Figure RE-GDA0002536505580000033
Figure RE-GDA0002536505580000033

建立量测方程时选取激光雷达和雷达这俩个传感器收集的数据来进行传感器融合预测车辆轨迹;When establishing the measurement equation, the data collected by the lidar and radar sensors are selected to perform sensor fusion to predict the vehicle trajectory;

首先是激光雷达,其的数据处于笛卡尔坐标系,可检测到车辆位置,但没有速度信息,其测量值为目标车辆的坐标(x,y);因此,激光雷达的测量模型仍然是线性的,其测量矩阵为:The first is lidar, whose data is in the Cartesian coordinate system, which can detect the vehicle position, but has no speed information, and its measurement value is the coordinates (x, y) of the target vehicle; therefore, the measurement model of lidar is still linear. , and its measurement matrix is:

Figure RE-GDA0002536505580000034
Figure RE-GDA0002536505580000034

相应的量测方程为:The corresponding measurement equation is:

Figure RE-GDA0002536505580000035
Figure RE-GDA0002536505580000035

其次是雷达,其测量的目标车辆数据处于极坐标系下,雷达的预测映射到测量空间是非线性的[16],其表达式为:The second is the radar, whose measured target vehicle data is in the polar coordinate system, and the radar's prediction is mapped to the measurement space is nonlinear [16] , and its expression is:

Figure RE-GDA0002536505580000041
Figure RE-GDA0002536505580000041

此时使用h(x)来表示这样一个非线性映射,那么在求解卡尔曼增益时也要将该非线性过程使用泰勒公式将其线性化,参照预测过程,求解h(x)的雅可比矩阵:At this time, h(x) is used to represent such a nonlinear mapping, then when solving the Kalman gain, the nonlinear process should also be linearized using Taylor's formula, and the Jacobian matrix of h(x) should be solved with reference to the prediction process. :

Figure RE-GDA0002536505580000042
Figure RE-GDA0002536505580000042

预测和修正模型中的未知参数JH,JA,Q,R(处理噪声,一般由厂商提供)都已求解完毕,之后就能根据所述网络模型对数据进行滤波处理,从而实现车辆未来位置的预测。The unknown parameters J H , JA , Q, R (processing noise, generally provided by the manufacturer) in the prediction and correction model have been solved, and then the data can be filtered according to the network model, so as to realize the future position of the vehicle Prediction.

进一步地,所述Step2中确定行人未来位置的具体步骤为:Further, the specific steps of determining the future position of the pedestrian in the Step2 are:

根据所述智能手机的GPS模块收集的信息来确定行人的位置,但是GPS模块的动态定位中还存在很多误差,例如卫星钟差等;为了消除误差,提高定位精度,将卡尔曼滤波方法应用于GPS动态滤波的“当前”统计模型;The location of the pedestrian is determined according to the information collected by the GPS module of the smartphone, but there are still many errors in the dynamic positioning of the GPS module, such as satellite clock errors, etc. In order to eliminate the errors and improve the positioning accuracy, the Kalman filter method is applied to "Current" statistical model of GPS dynamic filtering;

直接从GPS接收机输出的定位结果入手,将GPS定位的各种误差源在各个方向上造成的总误差等效为一个当前均值和一个符合一阶马尔科夫过程的有色噪声的和;根据“当前”统计模型,行人的状态向量选取为:Starting directly from the positioning results output by the GPS receiver, the total error caused by various error sources of GPS positioning in all directions is equivalent to the sum of a current mean and a colored noise that conforms to the first-order Markov process; according to " The current "statistical model", the pedestrian's state vector is selected as:

Figure RE-GDA0002536505580000043
Figure RE-GDA0002536505580000043

系统的状态方程可定为:The state equation of the system can be defined as:

Figure RE-GDA0002536505580000044
Figure RE-GDA0002536505580000044

其中,A为系统状态转移矩阵:Among them, A is the system state transition matrix:

A=diag[Ax,Ay,Az]A=diag[A x , A y , A z ]

Figure RE-GDA0002536505580000051
Figure RE-GDA0002536505580000051

Ay,Az依此类推;A y , A z and so on;

Figure RE-GDA0002536505580000052
Figure RE-GDA0002536505580000052

W(t)为:W(t) is:

=[0 0 Wax Wx 0 0 Way Wy 0 0 Waz Wz]T =[0 0 W ax W x 0 0 W ay W y 0 0 W az W z ] T

为系统噪声矢量;is the system noise vector;

τax,τay,τaz分别为加速度相关时间常数;τx,τy,τz分别为对应马尔科夫过程的相关时间常数;Wax,Way,Waz分别为

Figure RE-GDA0002536505580000053
的高斯白噪声;
Figure RE-GDA0002536505580000054
是三个坐标轴上的“当前”加速度均值;τ ax , τ ay , τ az are the acceleration-related time constants respectively; τ x , τ y , τ z are the relevant time constants corresponding to the Markov process; W ax , W ay , and W az are respectively
Figure RE-GDA0002536505580000053
Gaussian white noise;
Figure RE-GDA0002536505580000054
is the "current" mean acceleration on the three axes;

建立量测方程时取GPS接收机在3个坐标轴方向上的定位结果Lx,Ly,Lz为测量值,则:When establishing the measurement equation, take the positioning results L x , Ly , and L z of the GPS receiver in the directions of the three coordinate axes as the measurement values, then:

Lx=x+∈xlx L x = x + ∈ x + ω lx

Ly=y+∈yly L y =y+∈ yly

Lz=z+∈zlz L z = z+∈ zlz

写成矩阵形式:Written in matrix form:

L=[Lx,Ly,Lz]T L=[L x , Ly , L z ] T

V为观测噪声V=[ωlx,ωly,ωlz]T V is the observation noise V=[ω lx , ω ly , ω lz ] T

H为观测方程H is the observation equation

Figure RE-GDA0002536505580000055
Figure RE-GDA0002536505580000055

则相应的量测方程为Then the corresponding measurement equation is

L=HX+V (7)L=HX+V (7)

到此我们用GPS系统的观测方程建立了系统的测量方程,用“当前”统计模型和GPS系统的误差模型建立了系统的状态转移方程,下面根据上述系统方程和观测方程,设采样周期为T,通过离散化处理[18],建立离散卡尔曼滤波方程如下:So far, we have established the measurement equation of the system with the observation equation of the GPS system, and established the state transition equation of the system with the "current" statistical model and the error model of the GPS system. According to the above system equation and observation equation, let the sampling period be T. , through the discretization process [18] , the discrete Kalman filter equation is established as follows:

Figure RE-GDA0002536505580000061
Figure RE-GDA0002536505580000061

Figure RE-GDA0002536505580000062
Figure RE-GDA0002536505580000062

K(k+1)=P(k+1,k)HT(k+1)[H(k+1)P(k+1)HT(k+1)+R(k+1)]-1 K(k+1)=P(k+1,k)H T (k+1)[H(k+1)P(k+1)H T (k+1)+R(k+1)] -1

(9)(9)

P(k+1,k)=Φ(k+1,k)P(k)ΦT(k+1,k)+Q(k) (10)P(k+1, k)=Φ(k+1, k)P(k)Φ T (k+1, k)+Q(k) (10)

P(k+1)=[I-K(k+1)H(k+1)]P(k+1,k) (11)P(k+1)=[I-K(k+1)H(k+1)]P(k+1,k) (11)

其中Φl(k+1,k)=where Φ l (k+1, k)=

diag[Φlx(k+1,k),Φly(k+1,k),Φlz(k+1,k)]diag[Φ lx (k+1, k), Φ ly (k+1, k), Φ lz (k+1, k)]

Figure RE-GDA0002536505580000063
Figure RE-GDA0002536505580000063

Φly,Φlz依此类推;Φ ly , Φ lz and so on;

式(10)中Φ(k+1,k)为系统转移矩阵A的离散化矩阵:In formula (10), Φ(k+1, k) is the discretization matrix of the system transition matrix A:

Φ(k+1,k)=diag[Φx(k+1,k),Φy(k+1,k),Φz(k+1,k)]Φ(k+1, k)=diag[Φ x (k+1, k), Φ y (k+1, k), Φ z (k+1, k)]

Figure RE-GDA0002536505580000064
Figure RE-GDA0002536505580000064

Φy(k+1,k),Φz(k+1,k)依此类推,只需将(12)中τax,τx分别换成τay,τy和τaz,τz即可;Φ y (k+1, k), Φ z (k+1, k) and so on, just replace τ ax , τ x in (12) with τ ay , τ y and τ az , τ z is Can;

Q(k)为系统噪声协方差阵Q的离散化矩阵Q(k) is the discretization matrix of the system noise covariance matrix Q

Figure RE-GDA0002536505580000065
Figure RE-GDA0002536505580000065

R(k)为系统的观测噪声协方差阵R的离散化矩阵R(k) is the discretization matrix of the observation noise covariance matrix R of the system

Figure RE-GDA0002536505580000066
Figure RE-GDA0002536505580000066

Q(k),R(k)均为正定阵;Q(k) and R(k) are both positive definite matrices;

至此预测和修正模型中的未知参数都已求解完毕,接下来按照常规的卡尔曼滤波程序进行数据的相应处理即可。So far, the unknown parameters in the prediction and correction models have been solved, and then the data can be processed according to the conventional Kalman filter procedure.

进一步地,所述Step3中车辆行人碰撞避免算法的具体步骤为:Further, the specific steps of the vehicle-pedestrian collision avoidance algorithm in the Step 3 are:

Step3.1:车辆和行人设备周期性的通过自身传感器收集各自的状态信息,并通过RSU/基站将其上传至MEC服务器,在MEC服务器上根据行人和车辆的状态信息运行相应的卡尔曼滤波器预测出行人和车辆的轨迹;Step3.1: Vehicles and pedestrian equipment periodically collect their own state information through their own sensors, and upload it to the MEC server through the RSU/base station, and run the corresponding Kalman filter on the MEC server according to the state information of pedestrians and vehicles Predict the trajectories of pedestrians and vehicles;

所述Step3.1中,行人的运动需看成是非线性的,因为行人的运动有很大不确定性,而EKF便是预测非线性运动物体位置的滤波器;相对于行人,在道路上行驶的车辆由于道路的限制,其运动大致可以分为俩类,一类是直行,偏航角速度为零,可以将其看成是线性运动;一类是转弯此时偏航角速度不为零是非线性运动;In the above Step3.1, the movement of pedestrians needs to be regarded as nonlinear, because the movement of pedestrians has great uncertainty, and EKF is a filter for predicting the position of nonlinear moving objects; compared with pedestrians, driving on the road Due to the restrictions of the road, the motion of the vehicle can be roughly divided into two categories, one is straight, and the yaw angular velocity is zero, which can be regarded as a linear motion; the other is when the yaw angular velocity is not zero when turning is nonlinear sports;

Step3.2:运行碰撞避免算法计算出碰撞时间(TTC)并根据TTC的值判断车辆行人碰撞的可能性,来警告车辆驾驶员和行人;Step3.2: Run the collision avoidance algorithm to calculate the time to collision (TTC) and judge the possibility of vehicle-pedestrian collision according to the value of TTC to warn the vehicle driver and pedestrians;

所述TTC(time to collision)为:The TTC (time to collision) is:

Figure RE-GDA0002536505580000071
Figure RE-GDA0002536505580000071

L1,L2分别是车辆和行人到碰撞点的轨迹长度,V1,V2分别是车辆和行人当前的速度,tv,tp分别为车辆和行人到边缘服务器的通信延迟,根据此TTC的值来判断是否发生碰撞;L 1 , L 2 are the trajectory lengths of vehicles and pedestrians to the collision point respectively, V 1 , V 2 are the current speeds of vehicles and pedestrians, respectively, t v , t p are the communication delays from vehicles and pedestrians to the edge server, respectively, according to this The value of TTC to determine whether a collision occurs;

当TTC=0时,即如果当前车辆和行人不采取措施而保持当前状态行驶,则碰撞可能发生,此时,将警告信息发送给车辆,在事故发生前通知车辆采取措施,从而避免碰撞。When TTC=0, that is, if the current vehicle and pedestrian keep the current state without taking measures, a collision may occur. At this time, the warning information is sent to the vehicle, and the vehicle is notified to take measures before the accident occurs, so as to avoid the collision.

进一步地,所述Step3.2中,由于算法中存在着无法避免的误差,因此,设定TTCR(time to collision range)和阈值ξ,当Further, in the Step 3.2, due to the unavoidable error in the algorithm, set TTCR (time to collision range) and threshold ξ, when

TTCR≤ξTTCR≤ξ

就进行碰撞预警,从而消除轨迹预测产生的误差和车长的影响。Collision warning is carried out to eliminate the error caused by trajectory prediction and the influence of vehicle length.

本发明的有益效果是:The beneficial effects of the present invention are:

其一、本发明通过MEC服务器建立了当前网络环境下车辆和行人可行的通信方式,有效降低了通信延迟,车辆到MEC服务器延迟在10ms以下,行人设备到服务器的延迟在50ms以下,提高了车辆行人碰撞事故预警的准确率,有效降低了雨雾道路盲区等非视距情形下车辆和行人碰撞事故发生的概率。First, the present invention establishes a feasible communication method between vehicles and pedestrians in the current network environment through the MEC server, effectively reducing the communication delay. The accuracy of pedestrian collision early warning effectively reduces the probability of collision between vehicles and pedestrians in non-line-of-sight situations such as blind spots on rainy and foggy roads.

其二、本发明通过MEC服务器建立了当前网络环境下车辆和行人可行的通信方式,引入MEC极大的缓解了云服务器的负担,提高了系统的效率,增大了系统的吞吐量。Second, the present invention establishes a feasible communication mode between vehicles and pedestrians in the current network environment through the MEC server, and the introduction of the MEC greatly relieves the burden of the cloud server, improves the efficiency of the system, and increases the throughput of the system.

其三、本发明引入了行人设备——智能手机,克服了非视距限制,有效扩大了应用场景,将行人随身携带的可随时随地提供位置服务的设备加入到了车联网中,显著地提升了非机动道路使用者的安全,并且,引入智能手机势必对当前行人安全研究以车辆为重心的现状产生积极影响。Third, the present invention introduces a pedestrian device-smartphone, which overcomes the limitation of non-line-of-sight distance, effectively expands the application scenarios, and adds the device that pedestrians carry with them that can provide location services anytime and anywhere into the Internet of Vehicles, which significantly improves the performance of the Internet of Vehicles. The safety of non-motorized road users, and the introduction of smartphones is bound to have a positive impact on the current vehicle-focused status of pedestrian safety research.

其四、本发明还采用了卡尔曼滤波对车辆和行人未来的位置进行了预测,进一步的提高了车辆行人碰撞事故预警的准确率。Fourth, the present invention also uses Kalman filtering to predict the future positions of vehicles and pedestrians, which further improves the accuracy of vehicle-pedestrian collision accident warning.

附图说明Description of drawings

图1是本发明的系统模型图;1 is a system model diagram of the present invention;

图2是车速对碰撞预警成功率的影响数据结果图;Figure 2 is a graph showing the effect of vehicle speed on the success rate of collision warning;

图3是行人速度对碰撞预警成功率的影响数据结果图;Figure 3 is a graph of the impact data of pedestrian speed on the success rate of collision warning;

图4是对比试验数据结果图。Figure 4 is a graph of the results of the comparative test data.

具体实施方式Detailed ways

为了更详细的描述本发明和便于本领域人员的理解,下面结合附图以及实施例对本发明做进一步的描述,本部分的实施例用于解释说明本发明,便于理解的目的,不以此来限制本发明。In order to describe the present invention in more detail and facilitate the understanding of those skilled in the art, the present invention will be further described below with reference to the accompanying drawings and embodiments. Limit the invention.

实施例:Example:

一种基于边缘计算的车辆行人碰撞避免方法,如图1所示,包括以下步骤:A vehicle-pedestrian collision avoidance method based on edge computing, as shown in Figure 1, includes the following steps:

Step1:建立行人设备到车辆通信的网络模型:在十字路口或者道路边缘布置路侧单元(RSU)/基站、与所述路侧单元(RSU)/基站连接的移动边缘计算(MEC) 服务器,所述行人设备为具备GPS模块的智能手机,所述行人设备和车辆能在所述网络模型覆盖范围内实现通信;Step1: Establish a network model for pedestrian device-to-vehicle communication: Arrange roadside units (RSUs)/base stations at intersections or road edges, and mobile edge computing (MEC) servers connected to the roadside units (RSUs)/base stations, so The pedestrian device is a smart phone with a GPS module, and the pedestrian device and the vehicle can communicate within the coverage of the network model;

在这种模式下MEC服务器可以连接到路边单元(RSU)/基站上,车辆通过 LTE-V2X实现和RSU/基站的通信,行人则通过智能手机和RSU/基站进行通信;计算任务可以在车辆或者MEC服务器上进行;这样的系统既利用了传统的车联网中车辆的资源,并且在计算任务繁重时通过在MEC服务器上运行计算任务缓解了车辆上计算资源的压力,同时也利用了智能手机上的位置服务等资源,让智能手机参与到了行人安全保护的应用中;In this mode, the MEC server can be connected to the roadside unit (RSU)/base station, the vehicle communicates with the RSU/base station through LTE-V2X, and the pedestrian communicates with the RSU/base station through a smartphone; computing tasks can be performed in the vehicle Or on the MEC server; such a system not only utilizes the resources of the vehicle in the traditional Internet of Vehicles, but also relieves the pressure of the computing resources on the vehicle by running the computing task on the MEC server when the computing task is heavy, and also utilizes the smart phone. resources such as location services on the Internet, allowing smartphones to participate in the application of pedestrian safety protection;

Step2:根据车辆和行人设备当前信息和历史信息,利用卡尔曼滤波器来对车辆和行人未来位置分别进行预测;Step2: According to the current information and historical information of vehicles and pedestrian equipment, use Kalman filter to predict the future positions of vehicles and pedestrians respectively;

Step3:根据Step2中预测到的车辆和行人路径,根据相互运动关系,运行车辆行人碰撞避免算法,对可能发生的碰撞进行预警。Step3: According to the predicted vehicle and pedestrian paths in Step2, according to the mutual motion relationship, run the vehicle-pedestrian collision avoidance algorithm to give early warning of possible collisions.

进一步地,所述Step2中预测车辆未来位置的具体步骤为:Further, the specific steps of predicting the future position of the vehicle in the Step2 are:

假设汽车具有恒定转率和速度,选取车辆的状态向量为:Assuming that the car has a constant rotation rate and speed, the state vector of the selected vehicle is:

X=[x y v θ ω]T X=[xyv θ ω] T

θ为偏航角,是追踪的目标车辆在当前车辆坐标系下与x轴的夹角,逆时针方向为正,取值范围是[0,2π),ω是偏航角速度;行驶的汽车由于道路和交规等的限制其运动状态一般是直行或者转弯,这俩种运动状态都可以用此种状态向量来表示,当ω≠0时是非直线行驶,此时模型的状态转移函数为:θ is the yaw angle, which is the angle between the tracked target vehicle and the x-axis in the current vehicle coordinate system, the counterclockwise direction is positive, and the value range is [0, 2π), and ω is the yaw angular velocity; The motion states of roads and traffic regulations are generally straight or turning. Both motion states can be represented by this state vector. When ω≠0, it is non-straight travel. At this time, the state transfer function of the model is:

Figure RE-GDA0002536505580000091
Figure RE-GDA0002536505580000091

当ω=0时为直线行驶,此时模型的状态转移函数变为:When ω=0, it is a straight line, and the state transition function of the model becomes:

Figure RE-GDA0002536505580000092
Figure RE-GDA0002536505580000092

对于(1)和(2)的多元函数,使用多元泰勒级数进行展开,忽略高阶级数,只考虑利用雅克比矩阵进行线性化,在此模型中,对各个元素求偏导数可以得到雅可比矩阵JA,当ω≠0时:For the multivariate functions of (1) and (2), the multivariate Taylor series is used for expansion, the higher order number is ignored, and only the Jacobian matrix is used for linearization. In this model, the partial derivative of each element can be obtained. Jacobian Matrix J A , when ω≠0:

Figure RE-GDA0002536505580000093
Figure RE-GDA0002536505580000093

当ω=0时,雅可比矩阵为:When ω=0, the Jacobian matrix is:

Figure RE-GDA0002536505580000101
Figure RE-GDA0002536505580000101

计算处理噪声带来的不确定性Q;在所述网络模型中,噪声的引入主要来源于两处:直线加速度和偏航角加速度噪声,所述直线加速度和偏航角加速度噪声会影响状态量(x,y,v,θ,ω),所述直线加速度和偏航角加速度对状态的影响的表达式如下:Calculate and deal with the uncertainty Q caused by noise; in the network model, the introduction of noise mainly comes from two sources: linear acceleration and yaw angular acceleration noise, the linear acceleration and yaw angular acceleration noise will affect the state quantity (x, y, v, θ, ω), the expression of the influence of the linear acceleration and yaw angular acceleration on the state is as follows:

Figure RE-GDA0002536505580000102
Figure RE-GDA0002536505580000102

其中μa

Figure RE-GDA0002536505580000105
为直线上和转角上的加速度,Q就是处理噪声的协方差矩阵,其表达式为:where μ a ,
Figure RE-GDA0002536505580000105
is the acceleration on the straight line and on the corner, Q is the covariance matrix for dealing with noise, and its expression is:

Q=E[n·nT]=E[GμμTGT]=G·E[μμT]·GT (4)Q=E[n·n T ]=E[GμμT G T ]=G·E[μμ T ]·G T ( 4)

其中in

Figure RE-GDA0002536505580000103
Figure RE-GDA0002536505580000103

所以,在CTRV模型中的处理噪声的协方差矩阵Q的计算公式就是:Therefore, the calculation formula of the covariance matrix Q of processing noise in the CTRV model is:

Figure RE-GDA0002536505580000104
Figure RE-GDA0002536505580000104

建立量测方程时选取激光雷达和雷达这俩个传感器收集的数据来进行传感器融合预测车辆轨迹;When establishing the measurement equation, the data collected by the lidar and radar sensors are selected to perform sensor fusion to predict the vehicle trajectory;

首先是激光雷达,其的数据处于笛卡尔坐标系,可检测到车辆位置,但没有速度信息,其测量值为目标车辆的坐标(x,y);因此,激光雷达的测量模型仍然是线性的,其测量矩阵为:The first is lidar, whose data is in the Cartesian coordinate system, which can detect the vehicle position, but has no speed information, and its measurement value is the coordinates (x, y) of the target vehicle; therefore, the measurement model of lidar is still linear. , and its measurement matrix is:

Figure RE-GDA0002536505580000111
Figure RE-GDA0002536505580000111

相应的量测方程为:The corresponding measurement equation is:

Figure RE-GDA0002536505580000112
Figure RE-GDA0002536505580000112

其次是雷达,其测量的目标车辆数据处于极坐标系下,雷达的预测映射到测量空间是非线性的[16],其表达式为:The second is the radar, whose measured target vehicle data is in the polar coordinate system, and the radar's prediction is mapped to the measurement space is nonlinear [16] , and its expression is:

Figure RE-GDA0002536505580000113
Figure RE-GDA0002536505580000113

此时使用h(x)来表示这样一个非线性映射,那么在求解卡尔曼增益时也要将该非线性过程使用泰勒公式将其线性化,参照预测过程,求解h(x)的雅可比矩阵:At this time, h(x) is used to represent such a nonlinear mapping, then when solving the Kalman gain, the nonlinear process should also be linearized using Taylor's formula, and the Jacobian matrix of h(x) should be solved with reference to the prediction process. :

Figure RE-GDA0002536505580000114
Figure RE-GDA0002536505580000114

预测和修正模型中的未知参数JH,JA,Q,R(处理噪声,一般由厂商提供)都已求解完毕,之后就能根据所述网络模型对数据进行滤波处理,从而实现车辆未来位置的预测。The unknown parameters J H , JA , Q, R (processing noise, generally provided by the manufacturer) in the prediction and correction model have been solved, and then the data can be filtered according to the network model, so as to realize the future position of the vehicle Prediction.

进一步地,所述Step2中确定行人未来位置的具体步骤为:Further, the specific steps of determining the future position of the pedestrian in the Step2 are:

根据所述智能手机的GPS模块收集的信息来确定行人的位置,但是GPS模块的动态定位中还存在很多误差,例如卫星钟差等;为了消除误差,提高定位精度,将卡尔曼滤波方法应用于GPS动态滤波的“当前”统计模型;The location of the pedestrian is determined according to the information collected by the GPS module of the smartphone, but there are still many errors in the dynamic positioning of the GPS module, such as satellite clock errors, etc. In order to eliminate the errors and improve the positioning accuracy, the Kalman filter method is applied to "Current" statistical model of GPS dynamic filtering;

直接从GPS接收机输出的定位结果入手,将GPS定位的各种误差源在各个方向上造成的总误差等效为一个当前均值和一个符合一阶马尔科夫过程的有色噪声的和;根据“当前”统计模型,行人的状态向量选取为:Starting directly from the positioning results output by the GPS receiver, the total error caused by various error sources of GPS positioning in all directions is equivalent to the sum of a current mean and a colored noise that conforms to the first-order Markov process; according to " The current "statistical model", the pedestrian's state vector is selected as:

Figure RE-GDA0002536505580000121
Figure RE-GDA0002536505580000121

系统的状态方程可定为:The state equation of the system can be defined as:

Figure RE-GDA0002536505580000122
Figure RE-GDA0002536505580000122

其中,A为系统状态转移矩阵:Among them, A is the system state transition matrix:

A=diag[Ax,Ay,Az]A=diag[A x , A y , A z ]

Figure RE-GDA0002536505580000123
Figure RE-GDA0002536505580000123

Ay,Az依此类推;A y , A z and so on;

Figure RE-GDA0002536505580000124
Figure RE-GDA0002536505580000124

W(t)为:W(t) is:

=[0 0 Wax Wx 0 0 Way Wy 0 0 Waz Wz]T =[0 0 W ax W x 0 0 W ay W y 0 0 W az W z ] T

为系统噪声矢量;is the system noise vector;

τax,τay,τaz分别为加速度相关时间常数;τx,τy,τz分别为对应马尔科夫过程的相关时间常数;Wax,Way,Waz分别为

Figure RE-GDA0002536505580000125
的高斯白噪声;
Figure RE-GDA0002536505580000126
是三个坐标轴上的“当前”加速度均值;τ ax , τ ay , τ az are the acceleration-related time constants respectively; τ x , τ y , τ z are the relevant time constants corresponding to the Markov process; W ax , W ay , and W az are respectively
Figure RE-GDA0002536505580000125
Gaussian white noise;
Figure RE-GDA0002536505580000126
is the "current" mean acceleration on the three axes;

建立量测方程时取GPS接收机在3个坐标轴方向上的定位结果Lx,Ly,Lz为测量值,则:When establishing the measurement equation, take the positioning results L x , Ly , and L z of the GPS receiver in the directions of the three coordinate axes as the measurement values, then:

Lx=x+∈xlx L x = x + ∈ x + ω lx

Ly=y+∈yly L y =y+∈ yly

Lz=z+∈zlz L z = z+∈ zlz

写成矩阵形式:Written in matrix form:

L=[Lx,Ly,Lz]T L=[L x , Ly , L z ] T

V为观测噪声V=[ωlx,ωly,ωlz]T V is the observation noise V=[ω lx , ω ly , ω lz ] T

H为观测方程H is the observation equation

Figure RE-GDA0002536505580000131
Figure RE-GDA0002536505580000131

则相应的量测方程为 L=HX+V (7)Then the corresponding measurement equation is L=HX+V (7)

到此我们用GPS系统的观测方程建立了系统的测量方程,用“当前”统计模型和GPS系统的误差模型建立了系统的状态转移方程,下面根据上述系统方程和观测方程,设采样周期为T,通过离散化处理[18],建立离散卡尔曼滤波方程如下:So far, we have established the measurement equation of the system with the observation equation of the GPS system, and established the state transition equation of the system with the "current" statistical model and the error model of the GPS system. According to the above system equation and observation equation, let the sampling period be T. , through the discretization process [18] , the discrete Kalman filter equation is established as follows:

Figure RE-GDA0002536505580000132
Figure RE-GDA0002536505580000132

Figure RE-GDA0002536505580000133
Figure RE-GDA0002536505580000133

K(k+1)=P(k+1,k)HT(k+1)[H(k+1)P(k+1)HT(k+1)+R(k+1)]-1 K(k+1)=P(k+1,k)H T (k+1)[H(k+1)P(k+1)H T (k+1)+R(k+1)] -1

(9) (9)

P(k+1,k)=Φ(k+1,k)P(k)ΦT(k+1,k)+Q(k) (10)P(k+1, k)=Φ(k+1, k)P(k)Φ T (k+1, k)+Q(k) (10)

P(k+1)=[I-K(k+1)H(k+1)]P(k+1,k) (11)P(k+1)=[I-K(k+1)H(k+1)]P(k+1,k) (11)

其中Φl(k+1,k)=where Φ l (k+1, k)=

diag[Φlx(k+1,k),Φly(k+1,k),Φlz(k+1,k)]diag[Φ lx (k+1, k), Φ ly (k+1, k), Φ lz (k+1, k)]

Figure RE-GDA0002536505580000134
Figure RE-GDA0002536505580000134

Φly,Φlz依此类推;Φ ly , Φ lz and so on;

式(10)中Φ(k+1,k)为系统转移矩阵A的离散化矩阵:In formula (10), Φ(k+1, k) is the discretization matrix of the system transition matrix A:

Φ(k+1,k)=diag[Φx(k+1,k),Φy(k+1,k),Φz(k+1,k)]Φ(k+1, k)=diag[Φ x (k+1, k), Φ y (k+1, k), Φ z (k+1, k)]

Figure RE-GDA0002536505580000135
Figure RE-GDA0002536505580000135

Φy(k+1,k),Φz(k+1,k)依此类推,只需将(12)中τax,τx分别换成τay,τy和τaz,τz即可;Φ y (k+1, k), Φ z (k+1, k) and so on, just replace τ ax , τ x in (12) with τ ay , τ y and τ az , τ z is Can;

Q(k)为系统噪声协方差阵Q的离散化矩阵Q(k) is the discretization matrix of the system noise covariance matrix Q

Figure RE-GDA0002536505580000141
Figure RE-GDA0002536505580000141

R(k)为系统的观测噪声协方差阵R的离散化矩阵R(k) is the discretization matrix of the observation noise covariance matrix R of the system

Figure RE-GDA0002536505580000142
Figure RE-GDA0002536505580000142

Q(k),R(k)均为正定阵;Q(k) and R(k) are both positive definite matrices;

至此预测和修正模型中的未知参数都已求解完毕,接下来按照常规的卡尔曼滤波程序进行数据的相应处理即可。So far, the unknown parameters in the prediction and correction models have been solved, and then the data can be processed according to the conventional Kalman filter procedure.

进一步地,所述Step3中车辆行人碰撞避免算法的具体步骤为:Further, the specific steps of the vehicle-pedestrian collision avoidance algorithm in the Step 3 are:

Step3.1:车辆和行人设备周期性的通过自身传感器收集各自的状态信息,并通过RSU/基站将其上传至MEC服务器,在MEC服务器上根据行人和车辆的状态信息运行相应的卡尔曼滤波器预测出行人和车辆的轨迹;Step3.1: Vehicles and pedestrian equipment periodically collect their own state information through their own sensors, and upload it to the MEC server through the RSU/base station, and run the corresponding Kalman filter on the MEC server according to the state information of pedestrians and vehicles Predict the trajectories of pedestrians and vehicles;

所述Step3.1中,行人的运动需看成是非线性的,因为行人的运动有很大不确定性,而EKF便是预测非线性运动物体位置的滤波器;相对于行人,在道路上行驶的车辆由于道路的限制,其运动大致可以分为俩类,一类是直行,偏航角速度为零,可以将其看成是线性运动;一类是转弯此时偏航角速度不为零是非线性运动;In the above Step3.1, the movement of pedestrians needs to be regarded as nonlinear, because the movement of pedestrians has great uncertainty, and EKF is a filter for predicting the position of nonlinear moving objects; compared with pedestrians, driving on the road Due to the restrictions of the road, the motion of the vehicle can be roughly divided into two categories, one is straight, and the yaw angular velocity is zero, which can be regarded as a linear motion; the other is when the yaw angular velocity is not zero when turning is nonlinear sports;

Step3.2:运行碰撞避免算法计算出碰撞时间(TTC)并根据TTC的值判断车辆行人碰撞的可能性,来警告车辆驾驶员和行人;Step3.2: Run the collision avoidance algorithm to calculate the time to collision (TTC) and judge the possibility of vehicle-pedestrian collision according to the value of TTC to warn the vehicle driver and pedestrians;

所述TTC(time to collision)为:The TTC (time to collision) is:

Figure RE-GDA0002536505580000143
Figure RE-GDA0002536505580000143

L1,L2分别是车辆和行人到碰撞点的轨迹长度,V1,V2分别是车辆和行人当前的速度,tv,tp分别为车辆和行人到边缘服务器的通信延迟,根据此TTC的值来判断是否发生碰撞;L 1 , L 2 are the trajectory lengths of vehicles and pedestrians to the collision point respectively, V 1 , V 2 are the current speeds of vehicles and pedestrians, respectively, t v , t p are the communication delays from vehicles and pedestrians to the edge server, respectively, according to this The value of TTC to determine whether a collision occurs;

当TTC=0时,即如果当前车辆和行人不采取措施而保持当前状态行驶,则碰撞可能发生,此时,将警告信息发送给车辆,在事故发生前通知车辆采取措施,从而避免碰撞。When TTC=0, that is, if the current vehicle and pedestrian keep the current state without taking measures, a collision may occur. At this time, the warning information is sent to the vehicle, and the vehicle is notified to take measures before the accident occurs, so as to avoid the collision.

进一步地,所述Step3.2中,由于算法中存在着无法避免的误差,因此,设定TTCR(time to collision range)和阈值ξ,当Further, in the Step 3.2, due to the unavoidable error in the algorithm, set TTCR (time to collision range) and threshold ξ, when

TTCR≤ξTTCR≤ξ

就进行碰撞预警,从而消除轨迹预测产生的误差和车长的影响。Collision warning is carried out to eliminate the error caused by trajectory prediction and the influence of vehicle length.

从图2中可以看出,在TTC=0时车速越高预警准确率越低,因为车速越高,误差和延迟对其影响越大。当阈值ξ增加时,首先预警准确率会随之增加,并且达到最高值,这是因为阈值增大后会逐渐降低系统误差的影响。到达峰值之后预警准确率开始下降,因为时间差即阈值越大轨迹预测越不准确。并且阈值增大到一定程度之后,存在车辆和行人中一方或者双方已经通过碰撞点的情况,这种情况出现的概率会随着阈值的增大而增大。As can be seen from Figure 2, when TTC=0, the higher the vehicle speed, the lower the warning accuracy, because the higher the vehicle speed, the greater the impact of errors and delays on it. When the threshold ξ increases, the early warning accuracy rate will first increase and reach the highest value, because the increase of the threshold will gradually reduce the influence of the system error. After reaching the peak, the early warning accuracy begins to decline, because the time difference, that is, the larger the threshold, the less accurate the trajectory prediction is. And after the threshold is increased to a certain extent, there is a situation that one or both of the vehicle and the pedestrian have passed the collision point, and the probability of such a situation will increase with the increase of the threshold.

而且可以看到当车速越高时预警准确率曲线会越“陡峭”,车速越低曲线越平滑。因为车速越高越对延迟越敏感,车速越低轨迹预测越准确,相应的准确率较高。Moreover, it can be seen that the curve of the warning accuracy rate will be “steeper” when the vehicle speed is higher, and the curve will be smoother when the vehicle speed is lower. Because the higher the vehicle speed, the more sensitive it is to the delay, the lower the vehicle speed, the more accurate the trajectory prediction is, and the corresponding accuracy rate is higher.

由图3可以看出,行人速度和车辆速度对系统的影响大致相同,曲线走势也符合速度越快越“陡峭”,速度越慢越“平缓”的趋势。当TTC=0时速度和预警成功率成反相关。不同之处在于阈值ξ相同时,本部分的预警准确率一般是高于 5.3中的准确率,这是因为行人速度相较于车辆来说比较缓慢,对系统的影响不会像车辆那么明显。It can be seen from Figure 3 that the pedestrian speed and vehicle speed have roughly the same impact on the system, and the trend of the curve also conforms to the trend that the faster the speed is, the more "steep" it is, and the slower the speed is, the more "smooth". When TTC=0, speed and early warning success rate are inversely correlated. The difference is that when the threshold ξ is the same, the early warning accuracy rate in this part is generally higher than the accuracy rate in 5.3, because the pedestrian speed is slower than the vehicle, and the impact on the system will not be as obvious as the vehicle.

图4是对比试验结果图,本发明中碰撞点是根据车辆行人轨迹进行计算的,并且每个行人和车辆都有其各自的轨迹,所以车辆行人总数会影响系统性能。本部分就研究行人车辆总数对系统产生的影响。并与Avoiding car-pedestrian collisions using a VANETto cellular communication framework中的方法进行对比。从图4可以看到的是随着行人车辆总数增加,预警准确率会有所降低。因为此时会出现行人和车辆不会发生碰撞但只是因为距离较为接近就会进行预警的情况,这是由于有阈值ξ的影响系统有可能将此种情况判定为碰撞发生。Figure 4 is a comparison test result diagram. In the present invention, the collision point is calculated according to the trajectory of the vehicle and pedestrian, and each pedestrian and vehicle have their own trajectory, so the total number of vehicle pedestrians will affect the system performance. This part studies the impact of the total number of pedestrians and vehicles on the system. And contrast with the approach in Avoiding car-pedestrian collisions using a VANETto cellular communication framework. It can be seen from Figure 4 that as the total number of pedestrians and vehicles increases, the warning accuracy rate will decrease. Because at this time, there will be a situation where pedestrians and vehicles will not collide, but an early warning will be given only because the distance is relatively close. This is due to the influence of the threshold ξ. The system may determine this situation as a collision.

并且轨迹增多单个车辆或者行人可能会收到多个预警,此时预警总数会增加,当某一车辆收到准确预警信息进行减速成功避免碰撞后,其收到的其他预警信息都会作废,即单个车辆或行人的预警成功次数小于等于1,也就是预警成功次数的上限较小,而只会随着总数的增大而越来越大并且上限很高,所以预警成功率的变化趋势是随着行人车辆总数的增大而逐渐变小。和Avoiding car-pedestrian collisions using a VANET tocellular communication framework中的算法相比本文中的算法预警成功率更高。And the trajectory increases. A single vehicle or pedestrian may receive multiple warnings. At this time, the total number of warnings will increase. When a vehicle receives accurate warning information to decelerate and successfully avoid a collision, other warning information it receives will be invalid. The number of successful early warnings of vehicles or pedestrians is less than or equal to 1, that is, the upper limit of the number of successful early warnings is small, but will only increase with the increase of the total number and the upper limit is very high, so the change trend of the early warning success rate is as follows: The total number of pedestrians and vehicles increases and gradually decreases. Compared with the algorithm in Avoiding car-pedestrian collisions using a VANET tocellular communication framework, the algorithm in this paper has a higher early warning success rate.

以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。The specific embodiments of the present invention have been described in detail above in conjunction with the accompanying drawings, but the present invention is not limited to the above-mentioned embodiments, and can also be made within the scope of knowledge possessed by those of ordinary skill in the art without departing from the spirit of the present invention. Various changes.

Claims (5)

1.一种基于边缘计算的车辆行人碰撞避免方法,其特征在于:包括以下步骤:1. a vehicle-pedestrian collision avoidance method based on edge computing, is characterized in that: comprise the following steps: Step1:建立行人设备到车辆通信的网络模型:在十字路口或者道路边缘布置路侧单元(RSU)/基站、与所述路侧单元(RSU)/基站连接的移动边缘计算(MEC)服务器,所述行人设备为具备GPS模块的智能手机,所述行人设备和车辆能在所述网络模型覆盖范围内实现通信;Step1: Establish a network model for pedestrian device-to-vehicle communication: Arrange roadside units (RSUs)/base stations at intersections or road edges, and mobile edge computing (MEC) servers connected to the roadside units (RSUs)/base stations, so The pedestrian device is a smart phone with a GPS module, and the pedestrian device and the vehicle can communicate within the coverage of the network model; Step2:根据车辆和行人设备当前信息和历史信息,利用卡尔曼滤波器来对车辆和行人未来位置分别进行预测;Step2: According to the current information and historical information of vehicles and pedestrian equipment, use Kalman filter to predict the future positions of vehicles and pedestrians respectively; Step3:根据Step2中预测到的车辆和行人路径,根据相互运动关系,运行车辆行人碰撞避免算法,对可能发生的碰撞进行预警。Step3: According to the predicted vehicle and pedestrian paths in Step2, according to the mutual motion relationship, run the vehicle-pedestrian collision avoidance algorithm to give early warning of possible collisions. 2.根据权利要求1所述的一种基于边缘计算的车辆行人碰撞避免方法,其特征在于:所述Step2中预测车辆未来位置的具体步骤为:2. a kind of vehicle-pedestrian collision avoidance method based on edge computing according to claim 1, is characterized in that: the concrete step of predicting the future position of vehicle in described Step2 is: 假设汽车具有恒定转率和速度,选取车辆的状态向量为:Assuming that the car has a constant rotation rate and speed, the state vector of the selected vehicle is: X=[x y v θ ω]T X=[xyv θ ω] T θ为偏航角,是追踪的目标车辆在当前车辆坐标系下与x轴的夹角,逆时针方向为正,取值范围是[0,2π),ω是偏航角速度;行驶的汽车由于道路和交规等的限制其运动状态一般是直行或者转弯,这俩种运动状态都可以用此种状态向量来表示,当ω≠0时是非直线行驶,此时模型的状态转移函数为:θ is the yaw angle, which is the angle between the tracked target vehicle and the x-axis in the current vehicle coordinate system, the counterclockwise direction is positive, and the value range is [0, 2π), and ω is the yaw angular velocity; The motion states of roads and traffic regulations are generally straight or turning. Both motion states can be represented by this state vector. When ω≠0, it is non-straight travel. At this time, the state transfer function of the model is:
Figure RE-FDA0002536505570000011
Figure RE-FDA0002536505570000011
当ω=0时为直线行驶,此时模型的状态转移函数变为:When ω=0, it is a straight line, and the state transition function of the model becomes:
Figure RE-FDA0002536505570000012
Figure RE-FDA0002536505570000012
对于(1)和(2)的多元函数,使用多元泰勒级数进行展开,忽略高阶级数,只考虑利用雅克比矩阵进行线性化,在此模型中,对各个元素求偏导数可以得到雅可比矩阵JA,当ω≠0时:For the multivariate functions of (1) and (2), the multivariate Taylor series is used for expansion, the higher order number is ignored, and only the Jacobian matrix is used for linearization. In this model, the partial derivative of each element can be obtained. Jacobian Matrix J A , when ω≠0:
Figure RE-FDA0002536505570000021
Figure RE-FDA0002536505570000021
当ω=0时,雅可比矩阵为:When ω=0, the Jacobian matrix is:
Figure RE-FDA0002536505570000022
Figure RE-FDA0002536505570000022
计算处理噪声带来的不确定性Q;在所述网络模型中,噪声的引入主要来源于两处:直线加速度和偏航角加速度噪声,所述直线加速度和偏航角加速度噪声会影响状态量(x,y,v,θ,ω),所述直线加速度和偏航角加速度对状态的影响的表达式如下:Calculate and deal with the uncertainty Q caused by noise; in the network model, the introduction of noise mainly comes from two sources: linear acceleration and yaw angular acceleration noise, the linear acceleration and yaw angular acceleration noise will affect the state quantity (x, y, v, θ, ω), the expression of the influence of the linear acceleration and yaw angular acceleration on the state is as follows:
Figure RE-FDA0002536505570000023
Figure RE-FDA0002536505570000023
其中μa
Figure RE-FDA0002536505570000024
为直线上和转角上的加速度,Q就是处理噪声的协方差矩阵,其表达式为:
where μ a ,
Figure RE-FDA0002536505570000024
is the acceleration on the straight line and on the corner, Q is the covariance matrix for dealing with noise, and its expression is:
Q=E[n·nT]=E[GμμTGT]=G·E[μμT]·GT (4)Q=E[n·n T ]=E[GμμT G T ]=G·E[μμ T ]·G T ( 4) 其中in
Figure RE-FDA0002536505570000025
Figure RE-FDA0002536505570000025
所以,在CTRV模型中的处理噪声的协方差矩阵Q的计算公式就是:Therefore, the calculation formula of the covariance matrix Q of processing noise in the CTRV model is:
Figure RE-FDA0002536505570000026
Figure RE-FDA0002536505570000026
Figure RE-FDA0002536505570000031
Figure RE-FDA0002536505570000031
建立量测方程时选取激光雷达和雷达这俩个传感器收集的数据来进行传感器融合预测车辆轨迹;When establishing the measurement equation, the data collected by the lidar and radar sensors are selected to perform sensor fusion to predict the vehicle trajectory; 首先是激光雷达,其的数据处于笛卡尔坐标系,可检测到车辆位置,但没有速度信息,其测量值为目标车辆的坐标(x,y);因此,激光雷达的测量模型仍然是线性的,其测量矩阵为:The first is lidar, whose data is in the Cartesian coordinate system, which can detect the vehicle position, but has no speed information, and its measurement value is the coordinates (x, y) of the target vehicle; therefore, the measurement model of lidar is still linear. , and its measurement matrix is:
Figure RE-FDA0002536505570000032
Figure RE-FDA0002536505570000032
相应的量测方程为:The corresponding measurement equation is:
Figure RE-FDA0002536505570000033
Figure RE-FDA0002536505570000033
其次是雷达,其测量的目标车辆数据处于极坐标系下,雷达的预测映射到测量空间是非线性的[16],其表达式为:The second is the radar, whose measured target vehicle data is in the polar coordinate system, and the radar's prediction is mapped to the measurement space is nonlinear [16] , and its expression is:
Figure RE-FDA0002536505570000034
Figure RE-FDA0002536505570000034
此时使用h(x)来表示这样一个非线性映射,那么在求解卡尔曼增益时也要将该非线性过程使用泰勒公式将其线性化,参照预测过程,求解h(x)的雅可比矩阵:At this time, h(x) is used to represent such a nonlinear mapping, then when solving the Kalman gain, the nonlinear process should also be linearized using Taylor's formula, and the Jacobian matrix of h(x) should be solved with reference to the prediction process. :
Figure RE-FDA0002536505570000035
Figure RE-FDA0002536505570000035
预测和修正模型中的未知参数JH,JA,Q,R(处理噪声,一般由厂商提供)都已求解完毕,之后就能根据所述网络模型对数据进行滤波处理,从而实现车辆未来位置的预测。The unknown parameters J H , JA , Q, R (processing noise, generally provided by the manufacturer) in the prediction and correction model have been solved, and then the data can be filtered according to the network model, so as to realize the future position of the vehicle Prediction.
3.根据权利要求1所述的一种基于边缘计算的车辆行人碰撞避免方法,其特征在于:所述Step2中确定行人未来位置的具体步骤为:3. a kind of vehicle-pedestrian collision avoidance method based on edge computing according to claim 1, is characterized in that: the concrete step of determining the future position of pedestrian in described Step2 is: 根据所述智能手机的GPS模块收集的信息来确定行人的位置,但是GPS模块的动态定位中还存在很多误差,例如卫星钟差等;为了消除误差,提高定位精度,将卡尔曼滤波方法应用于GPS动态滤波的“当前”统计模型;The location of the pedestrian is determined according to the information collected by the GPS module of the smartphone, but there are still many errors in the dynamic positioning of the GPS module, such as satellite clock errors, etc. In order to eliminate the errors and improve the positioning accuracy, the Kalman filter method is applied to "Current" statistical model of GPS dynamic filtering; 直接从GPS接收机输出的定位结果入手,将GPS定位的各种误差源在各个方向上造成的总误差等效为一个当前均值和一个符合一阶马尔科夫过程的有色噪声的和;根据“当前”统计模型,行人的状态向量选取为:Starting directly from the positioning results output by the GPS receiver, the total error caused by various error sources of GPS positioning in all directions is equivalent to the sum of a current mean and a colored noise that conforms to the first-order Markov process; according to " The current "statistical model", the pedestrian's state vector is selected as:
Figure FDA0002399925360000041
Figure FDA0002399925360000041
系统的状态方程可定为:The state equation of the system can be defined as:
Figure FDA0002399925360000042
Figure FDA0002399925360000042
其中,A为系统状态转移矩阵:Among them, A is the system state transition matrix: A=diag[Ax,Ay,Az]A=diag[A x , A y , A z ]
Figure FDA0002399925360000043
Figure FDA0002399925360000043
Ay,Az依此类推;A y , A z and so on;
Figure FDA0002399925360000044
Figure FDA0002399925360000044
W(t)为:W(t) is: =[0 0 Wax Wx 0 0 Way Wy 0 0 Waz Wz]T =[0 0 W ax W x 0 0 W ay W y 0 0 W az W z ] T 为系统噪声矢量;is the system noise vector; τax,τay,τaz分别为加速度相关时间常数;τx,τy,τz分别为对应马尔科夫过程的相关时间常数;Wax,Way,Waz分别为
Figure FDA0002399925360000045
的高斯白噪声;
Figure FDA0002399925360000046
是三个坐标轴上的“当前”加速度均值;
τ ax , τ ay , τ az are the acceleration-related time constants respectively; τ x , τ y , τ z are the relevant time constants corresponding to the Markov process, respectively; W ax , W ay , and W az are respectively
Figure FDA0002399925360000045
Gaussian white noise;
Figure FDA0002399925360000046
is the "current" mean acceleration on the three axes;
建立量测方程时取GPS接收机在3个坐标轴方向上的定位结果Lx,Ly,Lz为测量值,则:When establishing the measurement equation, take the positioning results L x , Ly , and L z of the GPS receiver in the directions of the three coordinate axes as the measurement values, then: Lx=x+∈xlx L x = x + ∈ x + ω lx Ly=y+∈yly L y =y+∈ yly Lz=z+∈zlz L z = z+∈ zlz 写成矩阵形式:Written in matrix form: L=[Lx,Ly,Lz]T L=[L x , Ly , L z ] T V为观测噪声V=[ωlx,ωly,ωlz]T V is the observation noise V=[ω lx , ω ly , ω lz ] T H为观测方程H is the observation equation
Figure FDA0002399925360000051
Figure FDA0002399925360000051
则相应的量测方程为Then the corresponding measurement equation is L=HX+V (7)L=HX+V (7) 到此我们用GPS系统的观测方程建立了系统的测量方程,用“当前”统计模型和GPS系统的误差模型建立了系统的状态转移方程,下面根据上述系统方程和观测方程,设采样周期为T,通过离散化处理[18],建立离散卡尔曼滤波方程如下:So far, we have established the measurement equation of the system with the observation equation of the GPS system, and established the state transition equation of the system with the "current" statistical model and the error model of the GPS system. According to the above system equation and observation equation, let the sampling period be T. , through the discretization process [18] , the discrete Kalman filter equation is established as follows:
Figure FDA0002399925360000052
Figure FDA0002399925360000052
Figure FDA0002399925360000053
Figure FDA0002399925360000053
K(k+1)=P(k+1,k)HT(k+1)[H(k+1)P(k+1)HT(k+1)+R(k+1)]-1 K(k+1)=P(k+1,k)H T (k+1)[H(k+1)P(k+1)H T (k+1)+R(k+1)] -1 (9) (9) P(k+1,k)=Φ(k+1,k)P(k)ΦT(k+1,k)+Q(k) (10)P(k+1, k)=Φ(k+1, k)P(k)Φ T (k+1, k)+Q(k) (10) P(k+1)=[I-K(k+1)H(k+1)]P(k+1,k) (11)P(k+1)=[I-K(k+1)H(k+1)]P(k+1,k) (11) 其中Φl(k+1,k)=where Φ l (k+1, k)= diag[Φlx(k+1,k),Φly(k+1,k),Φlz(k+1,k)]diag[Φ lx (k+1, k), Φ ly (k+1, k), Φ lz (k+1, k)]
Figure FDA0002399925360000061
Figure FDA0002399925360000061
Φly,Φlz依此类推;Φ ly , Φ lz and so on; 式(10)中Φ(k+1,k)为系统转移矩阵A的离散化矩阵:In formula (10), Φ(k+1, k) is the discretization matrix of the system transition matrix A: Φ(k+1,k)=diag[Φx(k+1,k),Φy(k+1,k),Φz(k+1,k)]Φ(k+1, k)=diag[Φ x (k+1, k), Φ y (k+1, k), Φ z (k+1, k)]
Figure FDA0002399925360000062
Figure FDA0002399925360000062
Φy(k+1,k),Φz(k+1,k)依此类推,只需将(12)中τax,τx分别换成τay,τy和τaz,τz即可;Φ y (k+1, k), Φ z (k+1, k) and so on, just replace τ ax , τ x in (12) with τ ay , τ y and τ az , τ z is Can; Q(k)为系统噪声协方差阵Q的离散化矩阵Q(k) is the discretization matrix of the system noise covariance matrix Q
Figure FDA0002399925360000063
Figure FDA0002399925360000063
R(k)为系统的观测噪声协方差阵R的离散化矩阵R(k) is the discretization matrix of the observation noise covariance matrix R of the system
Figure FDA0002399925360000064
Figure FDA0002399925360000064
Q(k),R(k)均为正定阵;Q(k) and R(k) are both positive definite matrices; 至此预测和修正模型中的未知参数都已求解完毕,接下来按照常规的卡尔曼滤波程序进行数据的相应处理即可。So far, the unknown parameters in the prediction and correction models have been solved, and then the data can be processed according to the conventional Kalman filter procedure.
4.根据权利要求1所述的一种基于边缘计算的车辆行人碰撞避免方法,其特征在于:所述Step3中车辆行人碰撞避免算法的具体步骤为:4. a kind of vehicle-pedestrian collision avoidance method based on edge computing according to claim 1, is characterized in that: the concrete steps of vehicle-pedestrian collision avoidance algorithm in described Step3 are: Step3.1:车辆和行人设备周期性的通过自身传感器收集各自的状态信息,并通过RSU/基站将其上传至MEC服务器,在MEC服务器上根据行人和车辆的状态信息运行相应的卡尔曼滤波器预测出行人和车辆的轨迹;Step3.1: Vehicles and pedestrian equipment periodically collect their own state information through their own sensors, and upload it to the MEC server through the RSU/base station, and run the corresponding Kalman filter on the MEC server according to the state information of pedestrians and vehicles Predict the trajectories of pedestrians and vehicles; 所述Step3.1中,行人的运动需看成是非线性的,因为行人的运动有很大不确定性,而EKF便是预测非线性运动物体位置的滤波器;相对于行人,在道路上行驶的车辆由于道路的限制,其运动大致可以分为俩类,一类是直行,偏航角速度为零,可以将其看成是线性运动;一类是转弯此时偏航角速度不为零是非线性运动;In the above Step3.1, the movement of pedestrians needs to be regarded as nonlinear, because the movement of pedestrians has great uncertainty, and EKF is a filter for predicting the position of nonlinear moving objects; compared with pedestrians, driving on the road Due to the restrictions of the road, the motion of the vehicle can be roughly divided into two categories, one is straight, and the yaw angular velocity is zero, which can be regarded as a linear motion; the other is when the yaw angular velocity is not zero when turning is nonlinear sports; Step3.2:运行碰撞避免算法计算出碰撞时间(TTC)并根据TTC的值判断车辆行人碰撞的可能性,来警告车辆驾驶员和行人;Step3.2: Run the collision avoidance algorithm to calculate the time to collision (TTC) and judge the possibility of vehicle-pedestrian collision according to the value of TTC to warn the vehicle driver and pedestrians; 所述TTC(time to collision)为:The TTC (time to collision) is:
Figure FDA0002399925360000071
Figure FDA0002399925360000071
L1,L2分别是车辆和行人到碰撞点的轨迹长度,V1,V2分别是车辆和行人当前的速度,tv,tp分别为车辆和行人到边缘服务器的通信延迟,根据此TTC的值来判断是否发生碰撞;L 1 , L 2 are the trajectory lengths of vehicles and pedestrians to the collision point, respectively, V 1 , V 2 are the current speeds of vehicles and pedestrians, respectively, t v , t p are the communication delays from vehicles and pedestrians to the edge server, respectively, according to this The value of TTC to determine whether a collision occurs; 当TTC=0时,即如果当前车辆和行人不采取措施而保持当前状态行驶,则碰撞可能发生,此时,将警告信息发送给车辆,在事故发生前通知车辆采取措施,从而避免碰撞。When TTC=0, that is, if the current vehicle and pedestrian keep the current state without taking measures, a collision may occur. At this time, the warning information is sent to the vehicle, and the vehicle is notified to take measures before the accident occurs, so as to avoid the collision.
5.根据权利要求4所述的一种基于边缘计算的车辆行人碰撞避免方法,其特征在于:所述Step3.2中,由于算法中存在着无法避免的误差,因此,设定TTCR(time to collisionrange)和阈值ξ,当5. a kind of vehicle-pedestrian collision avoidance method based on edge computing according to claim 4, is characterized in that: in described Step3.2, because there is an unavoidable error in the algorithm, therefore, set TTCR (time to collisionrange) and threshold ξ, when TTCR≤ξTTCR≤ξ 就进行碰撞预警,从而消除轨迹预测产生的误差和车长的影响。Collision warning is carried out to eliminate the error caused by trajectory prediction and the influence of vehicle length.
CN202010143534.5A 2020-03-04 2020-03-04 Vehicle and pedestrian collision avoidance method based on edge calculation Pending CN111489585A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010143534.5A CN111489585A (en) 2020-03-04 2020-03-04 Vehicle and pedestrian collision avoidance method based on edge calculation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010143534.5A CN111489585A (en) 2020-03-04 2020-03-04 Vehicle and pedestrian collision avoidance method based on edge calculation

Publications (1)

Publication Number Publication Date
CN111489585A true CN111489585A (en) 2020-08-04

Family

ID=71794318

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010143534.5A Pending CN111489585A (en) 2020-03-04 2020-03-04 Vehicle and pedestrian collision avoidance method based on edge calculation

Country Status (1)

Country Link
CN (1) CN111489585A (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112562405A (en) * 2020-11-27 2021-03-26 山东高速建设管理集团有限公司 Radar video intelligent fusion and early warning method and system
CN112634354A (en) * 2020-12-21 2021-04-09 紫清智行科技(北京)有限公司 Road side sensor-based networking automatic driving risk assessment method and device
CN113012469A (en) * 2021-03-16 2021-06-22 浙江亚太机电股份有限公司 Intelligent traffic early warning system based on target recognition
CN113112805A (en) * 2021-04-16 2021-07-13 吉林大学 Intersection monitoring and early warning method based on base station communication and intersection camera positioning
CN113140132A (en) * 2021-04-20 2021-07-20 重庆邮电大学 Pedestrian anti-collision early warning system and method based on 5G V2X mobile intelligent terminal
CN113223325A (en) * 2021-03-26 2021-08-06 南京市德赛西威汽车电子有限公司 Method for safely passing signal-lamp-free intersection
CN113345248A (en) * 2021-04-29 2021-09-03 华设设计集团股份有限公司 Non-signal control intersection safety early warning and control method and system
CN113610263A (en) * 2021-06-18 2021-11-05 广东能源集团科学技术研究院有限公司 Power plant operating vehicle track estimation method and system
CN113727308A (en) * 2021-10-20 2021-11-30 湖北大学 Edge calculation unloading optimization method based on vehicle position prediction
CN113851017A (en) * 2021-08-19 2021-12-28 复旦大学 Pedestrian and vehicle identification and early warning multifunctional system based on road side RSU
CN114120715A (en) * 2020-08-31 2022-03-01 中移(成都)信息通信科技有限公司 Method, device and equipment for collision prevention of low-altitude aircraft and computer storage medium
CN115131990A (en) * 2022-05-20 2022-09-30 兆边(上海)科技有限公司 Early warning method and system for mixed flow collision risk of non-signal control intersection
CN115267868A (en) * 2022-09-27 2022-11-01 腾讯科技(深圳)有限公司 Positioning point processing method and device and computer readable storage medium
CN115457806A (en) * 2022-09-05 2022-12-09 广东轻工职业技术学院 Rear vehicle collision early warning system and method based on Internet of vehicles
CN115830860A (en) * 2022-11-17 2023-03-21 西部科学城智能网联汽车创新中心(重庆)有限公司 Traffic accident prediction method and device
CN116246489A (en) * 2023-02-06 2023-06-09 云控智行科技有限公司 Pedestrian protection method, device and equipment

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114120715A (en) * 2020-08-31 2022-03-01 中移(成都)信息通信科技有限公司 Method, device and equipment for collision prevention of low-altitude aircraft and computer storage medium
CN114120715B (en) * 2020-08-31 2023-02-28 中移(成都)信息通信科技有限公司 Method, device and equipment for collision prevention of low-altitude aircraft and computer storage medium
CN112562405A (en) * 2020-11-27 2021-03-26 山东高速建设管理集团有限公司 Radar video intelligent fusion and early warning method and system
CN112634354A (en) * 2020-12-21 2021-04-09 紫清智行科技(北京)有限公司 Road side sensor-based networking automatic driving risk assessment method and device
CN112634354B (en) * 2020-12-21 2021-08-13 紫清智行科技(北京)有限公司 Road side sensor-based networking automatic driving risk assessment method and device
CN113012469A (en) * 2021-03-16 2021-06-22 浙江亚太机电股份有限公司 Intelligent traffic early warning system based on target recognition
CN113223325A (en) * 2021-03-26 2021-08-06 南京市德赛西威汽车电子有限公司 Method for safely passing signal-lamp-free intersection
CN113112805B (en) * 2021-04-16 2022-08-02 吉林大学 An intersection monitoring and early warning method based on base station communication and intersection camera positioning
CN113112805A (en) * 2021-04-16 2021-07-13 吉林大学 Intersection monitoring and early warning method based on base station communication and intersection camera positioning
CN113140132A (en) * 2021-04-20 2021-07-20 重庆邮电大学 Pedestrian anti-collision early warning system and method based on 5G V2X mobile intelligent terminal
CN113140132B (en) * 2021-04-20 2023-11-03 西安华企众信科技发展有限公司 Pedestrian anti-collision early warning system and method based on 5G V2X mobile intelligent terminal
CN113345248A (en) * 2021-04-29 2021-09-03 华设设计集团股份有限公司 Non-signal control intersection safety early warning and control method and system
CN113610263B (en) * 2021-06-18 2023-06-09 广东能源集团科学技术研究院有限公司 Power plant operation vehicle track estimation method and system
CN113610263A (en) * 2021-06-18 2021-11-05 广东能源集团科学技术研究院有限公司 Power plant operating vehicle track estimation method and system
CN113851017A (en) * 2021-08-19 2021-12-28 复旦大学 Pedestrian and vehicle identification and early warning multifunctional system based on road side RSU
CN113727308B (en) * 2021-10-20 2023-06-30 湖北大学 Edge calculation unloading optimization method based on vehicle position prediction
CN113727308A (en) * 2021-10-20 2021-11-30 湖北大学 Edge calculation unloading optimization method based on vehicle position prediction
CN115131990A (en) * 2022-05-20 2022-09-30 兆边(上海)科技有限公司 Early warning method and system for mixed flow collision risk of non-signal control intersection
CN115457806A (en) * 2022-09-05 2022-12-09 广东轻工职业技术学院 Rear vehicle collision early warning system and method based on Internet of vehicles
CN115267868A (en) * 2022-09-27 2022-11-01 腾讯科技(深圳)有限公司 Positioning point processing method and device and computer readable storage medium
CN115267868B (en) * 2022-09-27 2023-09-19 腾讯科技(深圳)有限公司 Positioning point processing method and device and computer readable storage medium
CN115830860A (en) * 2022-11-17 2023-03-21 西部科学城智能网联汽车创新中心(重庆)有限公司 Traffic accident prediction method and device
CN115830860B (en) * 2022-11-17 2023-12-15 西部科学城智能网联汽车创新中心(重庆)有限公司 Traffic accident prediction method and device
CN116246489A (en) * 2023-02-06 2023-06-09 云控智行科技有限公司 Pedestrian protection method, device and equipment

Similar Documents

Publication Publication Date Title
CN111489585A (en) Vehicle and pedestrian collision avoidance method based on edge calculation
US20230400567A1 (en) Extended Object Tracking Using RADAR
CN110596694B (en) Complex environment radar multi-target tracking and road driving environment prediction method
CN107826104B (en) Method for providing information about a predicted driving intent of a vehicle
CN105206108B (en) A kind of vehicle collision prewarning method based on electronic map
WO2017217265A1 (en) Surrounding environment recognition device
CN110487288A (en) A kind of estimation method and carriage way estimating system of carriage way
CN112106065A (en) Predicting the state and position of an observed vehicle using optical tracking of wheel rotation
US11685406B2 (en) Vehicle control device, vehicle control method, and storage medium
CN109767619B (en) An intelligent network-connected pure electric vehicle driving condition prediction method
CN105682222A (en) Vehicle location positioning information fusion method based on vehicular ad hoc network
Song et al. A fusion strategy for reliable vehicle positioning utilizing RFID and in-vehicle sensors
Ramos et al. Cooperative target tracking in vehicular sensor networks
JP2019128614A (en) Prediction device, prediction method, and program
CN107085938A (en) A fault-tolerant planning method for local trajectory of intelligent driving based on lane lines and GPS following
Mo et al. Simulation and analysis on overtaking safety assistance system based on vehicle-to-vehicle communication
US20220144265A1 (en) Moving Track Prediction Method and Apparatus
JP2023548879A (en) Methods, devices, electronic devices and storage media for determining traffic flow information
CN104395944A (en) Carriageway recognition
CN116022163A (en) Automatic driving vehicle scanning matching and radar attitude estimator based on super local subgraph
KR20200133122A (en) Apparatus and method for preventing vehicle collision
CN117912295A (en) Vehicle data processing method and device, electronic equipment and storage medium
EP4160269A1 (en) Systems and methods for onboard analysis of sensor data for sensor fusion
WO2022078342A1 (en) Dynamic occupancy grid estimation method and apparatus
Song et al. Research on Target Tracking Algorithm Using Millimeter‐Wave Radar on Curved Road

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200804

WD01 Invention patent application deemed withdrawn after publication