CN106407506B - 一种道路三维线形建模方法及系统 - Google Patents

一种道路三维线形建模方法及系统 Download PDF

Info

Publication number
CN106407506B
CN106407506B CN201610717294.9A CN201610717294A CN106407506B CN 106407506 B CN106407506 B CN 106407506B CN 201610717294 A CN201610717294 A CN 201610717294A CN 106407506 B CN106407506 B CN 106407506B
Authority
CN
China
Prior art keywords
constraint
control point
order parameter
road
parameters
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.)
Expired - Fee Related
Application number
CN201610717294.9A
Other languages
English (en)
Other versions
CN106407506A (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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN201610717294.9A priority Critical patent/CN106407506B/zh
Publication of CN106407506A publication Critical patent/CN106407506A/zh
Application granted granted Critical
Publication of CN106407506B publication Critical patent/CN106407506B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Road Paving Structures (AREA)

Abstract

本发明提供种道路三维线形建模方法和系统,该方法包括:将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;根据所述多项式表达式,获取道路中线上任意点的三维坐标。该方法用统一的样条曲线来替代原有组合设计的三维线形,遵循道路线形的三维本质特性,并且可以大大简化原有设计流程。

Description

一种道路三维线形建模方法及系统
技术领域
本发明涉及公路几何设计领域,更具体地,涉及一种道路三维线形建模方法及系统。
背景技术
现有的公路三维选线方法包括收集和分析线路区域内的地质资料,设定多个可选路线,实地勘测后确定终选路线,这种线路设计的选线方式无法看到实地立体图像,依赖于选线人员的实际经验和技术水平,不适合工程要求较高的线路。
有些技术方案提供了三维线形设计,将道路线形设计分为两个阶段,首先是平曲线设计,然后是竖曲线设计,最后根据设计指标之间的平衡、视觉良好以及行车舒适等原则,将平曲线与竖曲线进行组合,最终得到公路中线的三维线形。
然而,公路中线本质是一条连续的空间曲线,传统组合设计公路线形的方法违背了公路线形的三维本质特性,是一种拼凑的思想。这种组合设计的公路线形在组合线形的节点处会出现线形几何不连续(曲率和挠率突变),从而导致汽车行驶时速度、加速度以及加速度的突变,影响行车安全和舒适性。
发明内容
本发明提供一种克服上述问题或者至少部分地解决上述问题的道路三维线形建模方法及系统。
第一方面,本发明提供一种道路三维线形建模方法,包括:步骤1,将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;步骤2,获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;步骤3,利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;步骤4,根据所述多项式表达式,获取道路中线上任意点的三维坐标。
步骤3还包括对相邻两控制点之间用多项式来连接,通过改变控制点的位置或者在两个控制点之间增加控制点,来调整两个控制点之间的路线。
步骤1的将道路按照平纵组合划分为多个路段中的多个路段包括直线与直坡、缓和曲线与直坡、圆曲线与直坡、直线与竖曲线、缓和曲线与竖曲线以及圆曲线与竖曲线六种组合方式。
步骤3中,控制点c和c+1的第p段样条多项式表达式为:
其中,u∈[c,c+1),p=c+1;p代表控制点c和控制点c+1之间的样条多项式;表示第p段位置u处的横坐标,表示第p段位置u处的纵坐标,表示第p段位置u处的竖坐标;apx,bpx,cpx,dpx,epx,fpx为第p段样条X多项式表达式的系数,apy,bpy,cpy,dpy,epy,fpy为第p段样条Y多项式表达式的系数,apz,bpz,cpz,dpz,epz,fpz为第p段样条Y多项式表达式的系数。
步骤4进一步包括:
根据相邻控制点处的坐标和切线限制,建立相邻控制点处的坐标方程、一阶参数约束方程和二阶参数约束方程;
根据所述坐标方程、一阶参数约束方程和二阶参数约束方程,建立表达式系数和控制点间参数约束值之间的方程组;
根据所述方程组,确定所述样条某点的坐标值。
步骤3进一步包括:
根据切线向量约束,获取一阶参数约束;
根据曲率、单位法向量约束以及上述一阶参数约束,获取二阶参数约束;
基于对曲率变化率和挠率的约束,根据所述一阶参数约束和二阶参数约束,获取三阶参数约束;
根据所述一阶参数约束、二阶参数约束和三阶参数约束,基于挠率变化率约束,获取四阶参数约束。
根据本发明的另一方面,提供一种道路三维线形建模系统,包括:
道路分段模块,用于将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;
几何约束获取模块,用于获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;
参数约束获取模块,用于利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;
坐标获取模块,用于根据所述多项式表达式,获取道路中线上任意点的三维坐标。
坐标获取模块进一步用于:
根据相邻控制点处的坐标和切线限制,建立相邻控制点处的坐标方程、一阶参数约束方程和二阶参数约束方程;
根据所述坐标方程、一阶参数约束方程和二阶参数约束方程,建立表达式系数和控制点间参数约束值之间的方程组;
根据所述方程组,确定所述样条某点的坐标值。
参数约束获取模块进一步用于:
根据切线向量约束,获取一阶参数约束;
根据曲率、单位法向量约束以及上述一阶参数约束,获取二阶参数约束;
基于对曲率变化率和挠率的约束,根据所述一阶参数约束和二阶参数约束,获取三阶参数约束;
根据所述一阶参数约束、二阶参数约束和三阶参数约束,基于挠率变化率约束,获取四阶参数约束
由上述技术方案可知,本发明提出的基于控制点处几何不变量约束的道路三维线形建模型方法,该方法用统一的样条曲线来替代原有组合设计的三维线形,遵循道路线形的三维本质特性,并且可以大大简化原有设计流程。
采用曲率和挠率两个空间曲线不变量作为几何指标来指导设计,不仅减少设计指标,而且可以建立起三维框架下道路几何要素和汽车动力学、运动学之间的关系。
本发明提供了控制点处一阶(切向量)、二阶(曲率)、三阶(曲率变化率或挠率)和四阶(挠率变化率)几何约束转化成参数约束的方法,利用五次多项式来连接满足控制点处曲率约束的道路线形,所构建的道路线形满足道路线形设计的基本要求—G2连续,使得行车安全和舒适。
附图说明
图1为根据本发明实施例的一种道路三维线形建模方法的流程图。
图2是根据本发明实施例的一种道路三维线形建模系统的结构图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1为根据本发明实施例的一种道路三维线形建模方法的流程图。
如图1所示的一种道路三维线形建模方法,包括:
S1,将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;
S2,获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;
S3,利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;
S4,根据所述多项式表达式,获取道路中线上任意点的三维坐标。
作为一种优选实施例,S3还包括对相邻两控制点之间用多项式来连接,通过改变控制点的位置或者在两个控制点之间增加控制点,来调整两个控制点之间的路线。
作为一种优选实施例,S1的将道路按照平纵组合划分为多个路段中的多个路段包括直线与直坡、缓和曲线与直坡、圆曲线与直坡、直线与竖曲线、缓和曲线与竖曲线以及圆曲线与竖曲线六种组合方式。
作为一种优选实施例,S3中,控制点c和c+1的第p段样条多项式表达式为:
其中,u∈[c,c+1),p=c+1;p代表控制点c和控制点c+1之间的样条多项式;表示第p段位置u处的横坐标,表示第p段位置u处的纵坐标,表示第p段位置u处的竖坐标;apx,bpx,cpx,dpx,epx,fpx为第p段样条X多项式表达式的系数,apy,bpy,cpy,dpy,epy,fpy为第p段样条Y多项式表达式的系数,apz,bpz,cpz,dpz,epz,fpz为第p段样条Y多项式表达式的系数。
作为一种优选实施例,S4进一步包括:
根据相邻控制点处的坐标和切线限制,建立相邻控制点处的坐标方程、一阶参数约束方程和二阶参数约束方程;
根据所述坐标方程、一阶参数约束方程和二阶参数约束方程,建立表达式系数和控制点间参数约束值之间的方程组;
根据所述方程组,确定所述样条某点的坐标值。
作为一种优选实施例,S3进一步包括:
根据切线向量约束,获取一阶参数约束;
根据曲率、单位法向量约束以及上述一阶参数约束,获取二阶参数约束;
基于对曲率变化率和挠率的约束,根据所述一阶参数约束和二阶参数约束,获取三阶参数约束;
根据所述一阶参数约束、二阶参数约束和三阶参数约束,基于挠率变化率约束,获取四阶参数约束。
本发明提出的基于控制点处几何不变量约束的道路三维线形建模型方法,该方法用统一的样条曲线来替代原有组合设计的三维线形,遵循道路线形的三维本质特性,并且可以大大简化原有设计流程。
采用曲率和挠率两个空间曲线不变量作为几何指标来指导设计,不仅减少设计指标,而且可以建立起三维框架下道路几何要素和汽车动力学、运动学之间的关系。
本发明提供了控制点处一阶(切向量)、二阶(曲率)、三阶(曲率变化率或挠率)和四阶(挠率变化率)几何约束转化成参数约束的方法,利用五次多项式来连接满足控制点处曲率约束的道路线形,所构建的道路线形满足道路线形设计的基本要求—G2连续,使得行车安全和舒适。
下面对本发明进行详细说明。
总的来说,本发明基于控制点处几何不变量的约束,提出了满足控制点处几何约束的道路中线建模方法。通过该方法进行设计时,不采用平纵分离的方式进行设计,而是用统一的三维解析模型来表示道路线形。所构建的道路中线满足公路设计的基本连续性要求—G2连续,特殊路段设计要求较高时可以通过局部调整达到G3连续。根据控制点处的约束条件求得两控制点间道路线形的多项式表达式,根据该多项式表达式可以求得中线上任意点的三维坐标,供施工放样参照。
在道路线形设计中,对于给定的3D地形,道路设计师会根据所给定地形选取道路必经的点,这些点称之为控制点。控制点处的三维坐标、曲率κ、挠率τ以及曲率变化率κ'、挠率变化率τ'可以根据设计需要来定义,其中曲率、挠率、曲率变化率、挠率变化率称之为几何约束。
为满足控制点处的几何约束,可以使用分段多项式来表示整条线路。对相邻两控制点之间用多项式来连接,如果需要对两个控制点之间的路线进行调整,可以通过改变控制点的位置,或者在两个控制点之间增加控制点。
控制点处给定的几何约束(曲率κ、挠率τ、曲率变化率κ'、挠率变化率τ')直接用于相邻控制点之间样条多项式系数的求解,计算量大而且复杂。以下说明如何利用微分计算将控制点处几何约束转换成便于控制点间多项式系数计算的参数约束
控制点处几何约束转化
1)一阶几何约束转化参数约束:
假定路线上控制点P的坐标(x,y,z),Frenet框架(T,N,B)已给定,以控制点P为例来讨论此点处的一阶约束。
在定义一阶导数参数约束时,用到切线向量公式:
一阶参数约束可以由上述公式变换得到:
式2中T(u)已知,大小可控,所以一阶约束随之可以进行不同的取值,但方向始终与给定切向量方向保持一致。
2)二阶几何约束转化参数约束
二阶参数约束可以由二阶几何约束(曲率和单位法向量)得到
为了能够更好的定义法向向量,对进行解析计算。根据链式法则,切线的一阶导数可以用下式表示:
对式5中右边第二项运用商的求导法则:
的一阶导数计算如下:
将式7的推导结果代入到式5中,得到:
式8可以通过Frenet公式和曲率建立起关系。为了进一步简化这种关系,假定是垂直的,以至于上式分子中值为0,最后可以得到二阶参数约束关于曲率和单位法向量的表达式:
3)三阶几何约束转化参数约束
a.挠率和三阶导数约束
讨论三阶约束时,我们首先考虑挠率,挠率的导数表达形式如下:
上式表明了挠率与一阶、二阶、三阶参数约束之间的关系,已经由单位切向量、单位法向量以及曲率所定义,所以三阶参数限制也可以确定。
将上式进行扩展,得到挠率的参数表达式为:
上式分子合并同类项,三阶参数约束可以分离出来:
由于一、二阶约束在前面已经被定义过,所以上式中除了x”',y”',z”'外其他全可以当作常数。式13可以简化成:
τ(u)=a0x”'+b0y”'+c0z”'
a0,b0,c0分别为关于一阶,二阶导数约束的常数,从而构成了τ(u)关于x”',y”',z”'的线性方程。
b.曲率变化率和三阶参数约束
控制点处的三阶参数约束也可以由曲率的一阶导数转化得来,如下式:
由于之前假定了一阶参数约束和二阶参数约束相互垂直,所示上式分子第二项为零,公式14可以简化成:
可以计算得到,现在用常数来替换,
式16的分子可以建立起三阶参数约束的线形方程的:
上式结合15式,得到最后结果:
现在,对的约束和对挠率τ的约束构成了三阶参数约束的线性方程:
4)四阶几何约束转化参数约束
最后,讨论挠率变化率τ'对构造路线的四阶约束。τ'的表达式如下:
由等式20可以看出,挠率变化率τ'可以是关于一阶,二阶,三阶,四阶导数的方程,由于一阶,二阶,三阶参数约束值由之前步骤可以计算出。因此,四阶参数约束可以分离并且解出来。公式20可以转换成:
一阶,二阶,三阶参数约束可以再一次看作常数,所以上式可以写作四阶参数约束的线性方程,如下:
κ,κ',τ和τ'这些几何约束可以根据设计需要来定义,通过计算得到满足参数限制的四阶几何约束。
满足控制点处约束的道路空间线形
N={0,1,2,…,l-1}代表路线上的必须经过的n个有序控制点。其中c∈C,c=0代表n个有序点的第一个点,c=l-1代表第n个节点。(xc,yc,zc)为控制点集C中的任意一点笛卡尔坐标。
X={x0,x1,…,xn-1}代表控制点集中n个点的横坐标,Y={y0,y1,…,yn-1}代表n个点的纵坐标,Z={z0,z1,…,zn-1}代表n个点的竖坐标。
根据公路线形设计的连续性要求,公路三维设计时为了保证行车舒适性和安全性,至少保证道路中线的曲率连续。控制点处几何约束转换成二阶参数约束,需要用五次样条来连接控制点。假定路线上有n个控制点,所以需要n-1段分段多项式。
连接控制点c和c+1的第p段样条多项式方程如下:
其中u∈[c,c+1),p=c+1。p代表控制点c和控制点c+1之间的样条多项式。例如p=2代表控制点1和2之间的第二段样条多项式。
表示第p段位置u处的横坐标
表示第p段位置u处的纵坐标
表示第p段位置u处的竖坐标。
apx,bpx,cpx,dpx,epx,fpx:第p段样条X多项式表达式的系数
apy,bpy,cpy,dpy,epy,fpy:第p段样条Y多项式表达式的系数
apz,bpz,cpz,dpz,epz,fpz:第p段样条Y多项式表达式的系数
第c个控制点与c+1个控制点我们用pc和pc+1来表示,对于控制点之间的样条多项式的X,Y,Z都有6个未知系数,分别需要6个方程来求解。以X表达式为例来说明系数求解。
根据控制点处的已知约束条件,可以知道求解系数所需要的方程。
点C与点C+1点处的坐标已知:2个方程
点C和点C+1处切线限制转化成一阶参数约束2个方程
点C和点C+1处切线限制转化成二阶参数约束2个方程
利用上述六个已知条件,可以建立起X表达式系数与控制点间参数约束值之间的方程组如下:
从而推得X表达式的系数为:
Y,Z多项式表达式系数也可以用同样的方法推得:
例子:
选取二级公路一段使用五次样条曲线进行设计计算。首先按照传统平纵组合的设计方法将公路划分为直线与直坡、缓和曲线与直坡、圆曲线与直坡、直线与竖曲线、缓和曲线与竖曲线以及圆曲线与竖曲线六种组合方式。
公路的起终点和相邻两组合路段的连接点作为计算控制点,用于新路线表达式计算。结合公路中心线上任意一点坐标计算公式:
计算控制点处的三维坐标、Frenet标架(T,N,B)以及曲率。将这些计算结果(控制点处的三维坐标、Frenet标架(T,N,B)以及曲率)作为控制点处的已知的几何约束,利用上述提出的几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,最终求得满足控制点处几何约束公路中线的分段多项式表达式。
根据平竖曲线设计表按平纵组合不同将实验段划分为12段,则控制点有13个,为方便计算将该实验段起点坐标设为(0,0,0),K5+429.3~K5+479.302直线段其实方位角设为90度。控制点处的几何约束计算结果如下:
表1几何约束
对控制点处的几何约束利用上述的一阶,二阶转化公式,求得控制点处的参数约束如下:
表2参数约束
控制点 一阶参数约束 二阶参数约束
K5+429.300 (0,-1,0.003) (0,0,0)
K5+479.302 (0,-1,0.003) (0,0,0)
K5+486.21 (0.00034,-1,0.003) (-0.000099,-0.000034,0)
K5+569.302 (0.058,-0.998,0.0085) (-0.0013,-0.000075,-0.000067)
K5+628.995 (0.135,-0.991,0.013) (-0.0013,-0.00018,-0.000067)
K5+693.79 (0.188,-0.982,0.0168) (-0.000355,-0.000069,-0.000067)
K5+718.995 (0.192,-0.981,0.0168) (0,0,0)
K6+035.486 (0.192,-0.981,0.0168) (0,0,0)
K6+125.486 (0.277,-0.961,0.0168) (-0.00187,-0.000546,0)
K6+204.909 (0.421,-0.907,0.0168) (-0.00176,-0.000828,0)
K6+284.909 (0.491,-0.871,0.0168) (0,0,0)
K6+361.89 (0.554,-0.832,0.0168) (-0.0016,-0.00107,0)
K6+374.905 (0.577,-0.817,0.0197) (-0.0018,-0.00131,-0.00022)
表2中的参数约束结合公式4-27,4-28,4-29建立起每组合路段X,Y,Z系数方程组,分别解方程组,最终求出每个组合路段X,Y,Z的多项式表达式系数,最终结果如表3所示:
表3二级公路多项式表达式子
本发明提出了基于控制点处几何不变量约束的道路三维线形建模型方法,该方法用统一的样条曲线来替代原有组合设计的三维线形,遵循道路线形的三维本质特性,并且可以大大简化原有设计流程。采用曲率和挠率两个空间曲线不变量作为几何指标来指导设计不仅减少设计指标,而且可以建立起三维框架下道路几何要素和汽车动力学、运动学之间的关系。
本发明介绍了控制点处一阶(切向量),二阶(曲率),三阶(曲率变化率或挠率),四阶(挠率变化率)几何约束转化成参数约束的方法,利用五次多项式来连接满足控制点处曲率约束的道路线形,所构建的道路线形满足道路线形设计的基本要求—G2连续,使得行车安全和舒适。
图2是根据本发明实施例的一种道路三维线形建模系统的结构图。
如图2所示的一种道路三维线形建模系统,包括:
道路分段模块21,用于将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;
几何约束获取模块22,用于获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;
参数约束获取模块23,用于利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;
坐标获取模块24,用于根据所述多项式表达式,获取道路中线上任意点的三维坐标。
作为一种优选实施例,坐标获取模块24进一步用于:
根据相邻控制点处的坐标和切线限制,建立相邻控制点处的坐标方程、一阶参数约束方程和二阶参数约束方程;
根据所述坐标方程、一阶参数约束方程和二阶参数约束方程,建立表达式系数和控制点间参数约束值之间的方程组;
根据所述方程组,确定所述样条某点的坐标值。
作为一种优选实施例,参数约束获取模块23进一步用于:
根据切线向量约束,获取一阶参数约束;
根据曲率、单位法向量约束以及上述一阶参数约束,获取二阶参数约束;
基于对曲率变化率和挠率的约束,根据所述一阶参数约束和二阶参数约束,获取三阶参数约束;
根据所述一阶参数约束、二阶参数约束和三阶参数约束,基于挠率变化率约束,获取四阶参数约束。
由于本发明的一种道路三维线形建模系统和一种道路三维线形建模方法是一一对应的,因此对一种道路三维线形建模系统不再详述。
本领域普通技术人员可以理解:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (8)

1.一种道路三维线形建模方法,其特征在于,包括:
步骤1,将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;
步骤2,获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;
步骤3,利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;
步骤4,根据所述多项式表达式,获取道路中线上任意点的三维坐标;
其中,步骤3中,控制点c和c+1的第p段样条多项式表达式为:
其中,u∈[c,c+1),p=c+1;p代表控制点c和控制点c+1之间的样条多项式;表示第p段位置u处的横坐标,表示第p段位置u处的纵坐标,表示第p段位置u处的竖坐标;apx,bpx,cpx,dpx,epx,fpx为第p段样条X多项式表达式的系数,apy,bpy,cpy,dpy,epy,fpy为第p段样条Y多项式表达式的系数,apz,bpz,cpz,dpz,epz,fpz为第p段样条Y多项式表达式的系数。
2.根据权利要求1所述的方法,其特征在于,步骤3还包括对相邻两控制点之间用多项式来连接,通过改变控制点的位置或者在两个控制点之间增加控制点,来调整两个控制点之间的路线。
3.根据权利要求1所述的方法,其特征在于,步骤1的将道路按照平纵组合划分为多个路段中的多个路段包括直线与直坡、缓和曲线与直坡、圆曲线与直坡、直线与竖曲线、缓和曲线与竖曲线以及圆曲线与竖曲线六种组合方式。
4.根据权利要求1所述的方法,其特征在于,步骤4进一步包括:
根据相邻控制点处的坐标和切线限制,建立相邻控制点处的坐标方程、一阶参数约束方程和二阶参数约束方程;
根据所述坐标方程、一阶参数约束方程和二阶参数约束方程,建立表达式系数和控制点间参数约束值之间的方程组;
根据所述方程组,确定所述样条某点的坐标值。
5.根据权利要求1所述的方法,其特征在于,步骤3进一步包括:
根据所述控制点处的切线向量约束,获取一阶参数约束;
根据所述控制点处的曲率、单位法向量约束以及上述一阶参数约束,获取二阶参数约束;
基于对所述控制点处的曲率变化率和挠率的约束,根据所述一阶参数约束和二阶参数约束,获取三阶参数约束;
根据所述一阶参数约束、二阶参数约束和三阶参数约束,基于所述控制点处的挠率变化率约束,获取四阶参数约束。
6.一种道路三维线形建模系统,其特征在于,包括:
道路分段模块,用于将道路按照平纵组合划分为多个路段,将道路的起终点和相邻路段的连接点作为控制点;
几何约束获取模块,用于获取控制点处的三维坐标、Frenet标架和曲率,将所述控制点处的三维坐标、Frenet标架和曲率作为控制点处的几何约束;
参数约束获取模块,用于利用几何约束转化参数约束的方法,将控制点处的几何约束转化成参数约束,获取满足控制点处几何约束公路中线的分段多项式表达式;
坐标获取模块,用于根据所述多项式表达式,获取道路中线上任意点的三维坐标;
其中,所述参数约束获取模块具体获取控制点c和c+1的第p段样条多项式表达式为:
其中,u∈[c,c+1),p=c+1;p代表控制点c和控制点c+1之间的样条多项式;表示第p段位置u处的横坐标,表示第p段位置u处的纵坐标,表示第p段位置u处的竖坐标;apx,bpx,cpx,dpx,epx,fpx为第p段样条X多项式表达式的系数,apy,bpy,cpy,dpy,epy,fpy为第p段样条Y多项式表达式的系数,apz,bpz,cpz,dpz,epz,fpz为第p段样条Y多项式表达式的系数。
7.根据权利要求6所述的系统,其特征在于,坐标获取模块进一步用于:
根据相邻控制点处的坐标和切线限制,建立相邻控制点处的坐标方程、一阶参数约束方程和二阶参数约束方程;
根据所述坐标方程、一阶参数约束方程和二阶参数约束方程,建立表达式系数和控制点间参数约束值之间的方程组;
根据所述方程组,确定所述样条某点的坐标值。
8.根据权利要求6所述的系统,其特征在于,参数约束获取模块进一步用于:
根据所述控制点处的切线向量约束,获取一阶参数约束;
根据所述控制点处的曲率、单位法向量约束以及上述一阶参数约束,获取二阶参数约束;
基于对所述控制点处的曲率变化率和挠率的约束,根据所述一阶参数约束和二阶参数约束,获取三阶参数约束;
根据所述一阶参数约束、二阶参数约束和三阶参数约束,基于所述控制点处的挠率变化率约束,获取四阶参数约束。
CN201610717294.9A 2016-08-24 2016-08-24 一种道路三维线形建模方法及系统 Expired - Fee Related CN106407506B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610717294.9A CN106407506B (zh) 2016-08-24 2016-08-24 一种道路三维线形建模方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610717294.9A CN106407506B (zh) 2016-08-24 2016-08-24 一种道路三维线形建模方法及系统

Publications (2)

Publication Number Publication Date
CN106407506A CN106407506A (zh) 2017-02-15
CN106407506B true CN106407506B (zh) 2019-08-27

Family

ID=58004656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610717294.9A Expired - Fee Related CN106407506B (zh) 2016-08-24 2016-08-24 一种道路三维线形建模方法及系统

Country Status (1)

Country Link
CN (1) CN106407506B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107909611A (zh) * 2017-10-01 2018-04-13 华南理工大学 一种利用微分几何理论提取空间曲线曲率特征的方法
CN108763775B (zh) * 2018-05-30 2022-03-18 东南大学 一种基于限制因素对高速公路进行立体改扩建的定线方法及其系统
CN112445207B (zh) * 2019-08-15 2023-07-11 纳恩博(北京)科技有限公司 位移的确定方法、装置、存储介质及电子装置
CN112146667B (zh) * 2020-09-29 2022-10-14 广州小鹏自动驾驶科技有限公司 一种车辆过渡轨迹的生成方法和装置
CN112257020B (zh) * 2020-11-04 2023-09-15 中国矿业大学 一种确定折线卷筒折线区设计圆心角的方法
CN113379873B (zh) * 2021-08-11 2021-11-09 北京赛目科技有限公司 一种道路曲线确定方法、装置、电子设备
CN114863385B (zh) * 2022-03-23 2023-04-07 禾多科技(北京)有限公司 道路曲面信息生成方法、装置、设备和计算机可读介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1391005A (zh) * 2001-12-30 2003-01-15 吴国雄 山区公路平面线形自动设计仿真系统
CN101814102A (zh) * 2010-03-16 2010-08-25 中交第二公路勘察设计研究院有限公司 一种基于高分辨卫星影像的公路测设方法
US8204680B1 (en) * 2007-07-05 2012-06-19 Navteq B.V. Method of operating a navigation system to provide road curvature
CN103310699A (zh) * 2012-03-16 2013-09-18 北京四维图新科技股份有限公司 一种提取道路的线形参数的方法
CN104732066A (zh) * 2015-02-16 2015-06-24 同济大学 路径约束条件下车辆行为时空演化的建模方法及其应用
CN104794286A (zh) * 2015-04-22 2015-07-22 杭州映墨科技有限公司 一种用于虚拟现实过山车场景的轨道构建和运动计算方法
CN105760392A (zh) * 2014-12-17 2016-07-13 高德软件有限公司 一种基于电子地图的数据处理方法和装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1391005A (zh) * 2001-12-30 2003-01-15 吴国雄 山区公路平面线形自动设计仿真系统
US8204680B1 (en) * 2007-07-05 2012-06-19 Navteq B.V. Method of operating a navigation system to provide road curvature
CN101814102A (zh) * 2010-03-16 2010-08-25 中交第二公路勘察设计研究院有限公司 一种基于高分辨卫星影像的公路测设方法
CN103310699A (zh) * 2012-03-16 2013-09-18 北京四维图新科技股份有限公司 一种提取道路的线形参数的方法
CN105760392A (zh) * 2014-12-17 2016-07-13 高德软件有限公司 一种基于电子地图的数据处理方法和装置
CN104732066A (zh) * 2015-02-16 2015-06-24 同济大学 路径约束条件下车辆行为时空演化的建模方法及其应用
CN104794286A (zh) * 2015-04-22 2015-07-22 杭州映墨科技有限公司 一种用于虚拟现实过山车场景的轨道构建和运动计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
空间曲线不变量约束下的公路几何线形三维设计理论与方法;葛婷;《中国博士学位论文全文数据库 工程科技II辑》;20160415(第04期);正文第87第1段-第111页第3段 *

Also Published As

Publication number Publication date
CN106407506A (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN106407506B (zh) 一种道路三维线形建模方法及系统
CN108088456B (zh) 一种具有时间一致性的无人驾驶车辆局部路径规划方法
CN113448335A (zh) 路径规划方法和装置、车辆和可读存储介质
CN109322221A (zh) 一种采用四次Bezier曲线分段构造公路三维空间线形的方法
EP2172748A2 (en) Creating geometry for advanced drive assistance systems
CN109241674A (zh) 一种智能网联车辆编队的多时延稳定性分析方法
CN102096919B (zh) 一种基于双向加权聚合的实时立体匹配方法
CN110017832A (zh) 一种基于Gauss解群优选的短弧初轨确定方法
CN103559797A (zh) 一种相位通用型干道双向绿波协调控制方法
CN105761270A (zh) 一种基于外极线距离变换的树型滤波立体匹配方法
CN114906173B (zh) 一种基于两点预瞄驾驶员模型的自动驾驶决策方法
CN113886764A (zh) 一种基于Frenet坐标系的智能车辆多场景轨迹规划方法
Patton Development of vehicle dynamics tools for motorsports
CN109166173B (zh) 基于bim的多精度三维测绘数据融合方法
CN112991421A (zh) 一种机器人视觉立体匹配方法
CN106570916A (zh) 用于公路三维线形设计的ph曲线生成方法和装置
CN115268504B (zh) 一种大型无人机仿地飞行路径规划方法
CN105096328B (zh) 一种基于ios‑pso的摄像机自标定方法
Bianco et al. Velocity planning for autonomous vehicles
CN113741179B (zh) 一种面向异构车辆的统一运动规划方法和系统
CN116295398A (zh) 一种基于协同聚集和分支偏差的路径优化方法
CN106251307A (zh) 基于 Extended Smoothlets (ExSmoothlets)变换的图像近似方法
CN102930176B (zh) 大地水准面模型九网格五样条插值方法
CN116331256A (zh) 分布式驱动电动汽车换道轨迹规划方法、设备及存储介质
CN113467480B (zh) 一种用于无人驾驶方程式的全局路径规划算法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190827

Termination date: 20200824