CN1299220C - 数字道路地图自动拼接方法 - Google Patents
数字道路地图自动拼接方法 Download PDFInfo
- Publication number
- CN1299220C CN1299220C CNB2004100183066A CN200410018306A CN1299220C CN 1299220 C CN1299220 C CN 1299220C CN B2004100183066 A CNB2004100183066 A CN B2004100183066A CN 200410018306 A CN200410018306 A CN 200410018306A CN 1299220 C CN1299220 C CN 1299220C
- Authority
- CN
- China
- Prior art keywords
- road
- vector
- intersection
- vectors
- map
- 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
Links
Images
Landscapes
- Processing Or Creating Images (AREA)
- Instructional Devices (AREA)
- Image Processing (AREA)
- Traffic Control Systems (AREA)
Abstract
一种数字道路地图自动拼接方法,通过路口中心点的坐标计算得到两个待拼接数字道路地图的坐标变换关系,用得到的坐标变化关系将待拼接的两幅地图配准到同一坐标系中并计算其重合部分,找出地图中重合的路口矢量对、道路矢量对并确立路口矢量和道路矢量的拓扑关系,将重合的道路拼接为新的道路矢量,将重合的道路路口合并为新的道路路口矢量,并将新生成的道路矢量和道路路口矢量添加到新建的数字道路地图中,将非重合的道路矢量和非重合的道路路口矢量直接添加到新建的数字道路地图中,完成数字道路地图的自动拼接。本发明的操作步骤简单,提高了数字道路地图拼接的精度和速度,降低对计算机硬件系统性能的依赖。
Description
技术领域
本发明涉及一种数字道路地图自动拼接方法,用于交通设备的自动导航。属于先进制造与自动化领域。
背景技术
数字道路地图在智能交通、车载导航、车辆定位和城市规划等应用中有极其重要的位置。目前数字道路地图的制作有很大一部分是通过对纸面地图扫描,进而用计算机数字化处理以得到数字道路地图。然而,受扫描仪尺寸的限制,采用扫描纸面地图得到数字道路地图就不可避免地涉及到地图的拼接问题。数字道路地图的拼接方法有两种:一是手工拼接,该方法通过操作人员的判断来确定待拼接图的对应位置关系,其拼接质量完全由操作人员的经验控制掌握,人为因素很大,这是主要的误差来源。二是自动拼接,其优点是拼接精度和速度的大量提高。
数字道路地图自动拼接的最关键部分是两幅待拼接图上数字道路之间的拓扑连接关系。因此,两幅待拼接地图的配准和确定它们上面数字道路拓扑连接关系是方案设计的核心。传统的拼接方法是先扫描得到光栅图对光栅图进行拼接,然后再对拼接后的光栅图进行数字道路矢量化处理。这类拼接方法的不利之处在于:由于拼接后的光栅图很大,因此在数字道路矢量化的处理时间、处理速度以及所使用的计算机内存容量上受计算机硬件性能的约束很大,同时也会增加数字道路矢量化算法的空间复杂度和时间复杂度。(R.Kasturi,et.al,“A System for Interpretation of Line Drawings”,IEEE Trans.PAMI,Vol12,No.10,Oct.1990,pp978-992.)
发明内容
本发明的目的在于针对现有技术的不足,提出一种新的数字道路地图自动拼接方法,解决拼接精度、速度的问题和对计算机硬件性能的依赖。
为实现这样的目的,本发明的技术方案中,利用二维平面的三点定位技术,通过对数字道路地图中路口中心点的坐标计算得到两个待拼接数字道路地图的坐标变换关系,使用得到的坐标变化关系将待拼接的两幅地图配准到同一坐标系中并计算出它们的重合部分,找出地图中重合的路口矢量对、道路矢量对并确立路口矢量和道路矢量的拓扑关系,将重合的道路拼接为一个新的道路矢量,将重合的道路路口合并为一个新的道路路口矢量,并将新生成的道路矢量和新生成的道路路口矢量添加到新建的数字道路地图中,将非重合的道路矢量和非重合的道路路口矢量直接添加到新建的数字道路地图中,完成数字道路地图的自动拼接。
本发明主要包括以下几个具体步骤:
1.地图的配准。
本发明利用二维平面的三点定位技术,首先在两幅待拼接数字道路地图的重合区域中,选取三对不在同一条直线上的道路的路口,运用最小二乘法计算标准化的道路路口中心坐标,对三对道路路口中心坐标建立方程组,从而将地图的配准问题抽象成在已知平面上求坐标变换的几何模型。解该方程组得到两幅需要拼接的数字道路地图的坐标对应转换关系。根据该坐标转换关系和待拼接图的尺寸按高度将拼接图像分为三部分,自左到右逐行处理,分别读取待拼接地图内容并转换到新的坐标系中同时确定待拼接地图的重合区域。
2.确定需要连接的道路对及它们的拓扑关系。
根据两幅待拼接数字道路地图在同一坐标系统中的位置和它们的重合区域范围,采用广度优先算法沿数字道路网搜索查找出所有位于重合部位的道路路口矢量和道路矢量。根据道路路口矢量的坐标将重合的道路路口矢量配对成需要合并的道路路口矢量对。根据道路矢量的坐标以及该道路矢量的方向角将重合的道路矢量配对成需要拼接的道路对。根据路口矢量所连接的道路矢量确立道路路口矢量和道路矢量的拓扑关系。
3.道路对的拼接。
对于需要拼接的道路对,它的拼接存在以下三种情况:(1)当两条道路完全重合时,将其中的一条道路矢量作为新的道路矢量并将它连接到相应的两个道路路口矢量中,删除另外一条道路矢量;(2)当两条道路部分重合时,删除一条道路上重合的道路段,并将剩余道路段连接到另外一条道路矢量上作为一条新的道路矢量,最后将新生成的道路矢量连接到相应的两个道路路口矢量中;(3)当两条道路完全不重合时,将这两条道路矢量直接连接为一条新的道路矢量,再将连接后的新的道路矢量连接到相应的两个道路路口矢量中。依次对所有道路对做以上拼接直到所有道路对都已经拼接为止。
4.道路路口对的合并。
新建一个道路路口矢量作为合并后的道路路口。对道路路口对所连接的所有道路进行依次判断,如果属于已经合并的道路对则将该道路对对应的新生成的道路矢量添加到该新建的道路路口矢量中;如果不属于已经合并的道路对则将该道路直接添加到该新建的道路路口矢量中。为道路路口矢量中的参数依次赋值。对所有道路路口对做以上合并直到所有道路路口对都已经合并为止。最后将新生成的道路路口矢量和新生成的道路矢量添加到新建的数字道路地图中,完成数字道路地图的自动拼接。
利用本发明的拼接方法实现数字道路地图的拼接,可简化操作步骤,提高拼接精度和速度,降低对计算机硬件系统性能的依赖。
附图说明
图1为本发明实施例中待拼接的数字道路地图A。
图2为本发明实施例中待拼接的数字道路地图B。
图3为采用本发明方法拼接后的数字道路地图C。
图1、图2、图3中:a、b、c所标识的圆圈中的道路路口是用于拼接的标志点。
图4为本发明方法中采用的待拼接图像的8种图像相对关系。
其中:Dx表示两副图像x方向上的相对关系,Dy表示两副图像y方向上的相对关系。
图5为本发明实施例中待拼接的两幅图像相对关系。
其中:Wa表示图像A的宽度,Ha表示图像A的高度,Wb表示图像B的宽度,Hb表示图像B的高度,Dx表示两副图像x方向上的相对关系,Dy表示两副图像y方向上的相对关系。
图6为计算道路位置和方向角的示意图。
图6中,a是道路路口中心点,a1、a2、a3分别标识道路路口所连接的三条道路。
图7为道路合并示意图。
图7中,A和B分别表示两条需要合并的道路。a、b、b’、c、c’、d分别标识道路上的点,其中b、b’和c、c’是道路上同一个点分别在两条道路段上的位置。
具体实施方式
为了更好地理解本发明的技术方案,以下结合实施例和附图作进一步的详细描述。
1.地图的配准。
图1图2分别为本发明实施例中待拼接的两幅数字道路地图A和B。
本发明利用二维平面的三点定位技术,首先在两幅待拼接数字道路地图的重合区域中,选取三对不在同一条直线上的道路的路口(附图1、2中的a、b、c点),运用最小二乘法计算标准化的道路路口矢量的路口中心坐标,对三对道路路口中心坐标建立方程组,从而将地图的配准问题抽象成在已知平面上求坐标变换的几何模型。解该方程组得到两幅需要拼接的数字道路地图的坐标对应转换关系:
分析两幅图像的坐标关系,由两幅图像的大小、位置,其对应关系有8种情况,如附图4所示。根据Dx和Dy的正负号及待拼接图A和图B的尺寸关系,按高度将拼接图像分为仅有图A、既有图A又有图B和仅有图B的三部分,对这三个部分分别自左到右逐行处理(如图5)。设A图左上角为原点(0,0),分以下三步:如果仅有图A部分,即y∈「0,Dy」读取图A的[0,0]→「Wa,Dy」部分将其转换到新的坐标系中;如果既有图A又有图B,即y∈(Dy,Ha),则以Dx为界先读取图A的0,Dy→Dx,Ha部分将其转换到新的坐标系中,然后读取图B的[0,0]→「Wb,Ha-Dy」部分将其转换到新的坐标系中,同时标记该区域为重合区域;如果仅有图B部分,即y∈Ha,Dy+Hb,读取图B的0,Ha-Dy→Wb,Hb部分并将其转换到新的坐标系中。
循环执行以上三步,直到将两幅待拼接地图全部转换到同一个坐标系中并标记出重合区域。其余七种图像对应关系,也是根据此法进行,只是其所分区间和内容有所变化而已。
2.确定需要连接的道路对及它们的拓扑关系。
根据两幅待拼接数字道路地图在同一坐标系统中的位置和它们的重合区域范围,调用广度优先算法(BFS)沿数字道路网搜索查找出所有位于重合部位的道路路口矢量和道路矢量。根据道路路口矢量的坐标找出附图1中道路路口(a,b,c)在附图2中对应的道路路口(a,b,c),将重合的道路路口矢量配对成需要合并的道路路口矢量对。根据道路矢量的坐标以及该道路矢量的方向角(公式2)找出附图1中道路路口a所连接的每条道路(如附图6中的a1、a2、a3)在附图2中对应的道路,将重合的道路矢量配对成需要合并的道路对。
根据路口和道路的连接关系确立道路路口矢量和道路矢量的拓扑关系。
3.道路对的拼接。
对于需要拼接的道路对,它的拼接存在以下三种情况(如附图7):(1)当两条道路完全重合时,即点a、b、b’重合、点c、c’、d重合,将其中的一条道路矢量(abcd)作为新的道路矢量并将它连接到相应的两个道路路口矢量中,删除另外一条道路矢量(b’c’d);(2)当两条道路部分重合时,即点b、b’重合、点c、c’重合,删除一条道路上重合的道路段(b’c’),并将剩余道路段(c’d)连接到另外一条道路矢量(abc)上作为一条新的道路矢量(abcc’d),最后将新生成的道路矢量连接到相应的两个道路路口矢量中;(3)当两条道路完全不重合时,即点a、b、b’、d、c、c’都不重合,将这两条道路矢量直接连接为一条新的道路矢量(a、b、b’、c、c’、d),再将连接后的新的道路矢量连接到相应的两个道路路口矢量中。其基本算法如下:
(1)对道路A,从端点a开始扫描并判断该点是否也在道路B上。
(2)标记第一个既属于道路A也属于道路B的点为b和b’。
(3)对道路B从b’的开始删除直到最后一个既属于道路A也属于道路B的点c’。
(4)连接道路点c和道路点c’。
(5)将拼接后的道路(abcd)添加的对应的道路路口中。
依次对所有道路对做以上拼接直到所有道路对都已经拼接为止。
4.道路路口对的合并。
在附图3中新建道路路口矢量(a,b,c)作为合并后的道路路口。对道路路口对所连接的所有道路(如附图6中的a1、a2、a3)进行依次判断,如果属于已经合并的道路对则将该道路对对应的新生成的道路矢量添加到该新建的道路路口矢量中;如果不属于已经合并的道路对则将该道路直接添加到该新建的道路路口矢量中。对所有道路路口对做以上合并直到所有道路路口对都已经合并为止。最后将新生成的道路路口矢量和新生成的道路矢量添加到新建的数字道路地图中,完成数字道路地图的自动拼接。
Claims (1)
1、一种数字道路地图自动拼接方法,其特征在于包括如下具体步骤:
1)地图的配准:首先在两幅待拼接地图的重合区域中选取三对不在同一条直线上的道路路口,运用最小二乘法计算路口中心坐标,对三对道路路口中心坐标建立方程组确定两数字地图的坐标转换关系,根据坐标转换关系和待拼接地图尺寸将待拼接地图转换到新坐标系中,同时确定待拼接地图的重合区域;
2)确定需要连接的道路对及路口矢量和道路矢量的拓扑关系:根据两幅待拼接数字道路地图在新坐标系中的位置和重合区域,采用广度优先算法找出所有位于重合区域的路口矢量和道路矢量,根据路口矢量坐标将重合的路口矢量配对成需合并路口矢量对,根据道路矢量坐标以及道路矢量方向角将重合的道路矢量配对成需拼接道路对,根据路口所连接的道路确立路口矢量和道路矢量的拓扑关系;
3)道路对的拼接:对需要拼接的道路对分三种情况处理:(1)当两条道路完全重合时,将其中的一条道路矢量作为新的道路矢量并连接到相应的两个路口矢量中,删除另外一条道路矢量;(2)当两条道路部分重合时,删除一条道路上重合的道路段,将剩余道路段连到另外一条道路矢量上作为一条新的道路矢量,将新生成的道路矢量连接到相应的两个路口矢量中;(3)当两条道路完全不重合时,将这两条道路矢量直接连接为一条新的道路矢量,将新道路矢量连接到相应的两个路口矢量中;依次对所有道路对做以上拼接直到所有道路对都已经拼接为止;
4)道路路口对的合并:新建一个路口矢量作为合并后的道路路口,对道路路口对所连的所有道路进行依次判断,如属于已经合并的道路对则将该道路对对应的新道路矢量添加到新建的上述路口矢量中,如不属于已经合并的道路对则将该道路直接添加到新建的上述路口矢量中;对所有道路路口对做以上合并直到所有道路路口对都已经合并为止,最后将新生成的道路路口矢量和新生成的道路矢量添加到新建的数字道路地图中,完成数字道路地图的自动拼接。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100183066A CN1299220C (zh) | 2004-05-13 | 2004-05-13 | 数字道路地图自动拼接方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100183066A CN1299220C (zh) | 2004-05-13 | 2004-05-13 | 数字道路地图自动拼接方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1570917A CN1570917A (zh) | 2005-01-26 |
| CN1299220C true CN1299220C (zh) | 2007-02-07 |
Family
ID=34479454
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2004100183066A Expired - Fee Related CN1299220C (zh) | 2004-05-13 | 2004-05-13 | 数字道路地图自动拼接方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1299220C (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630463B (zh) * | 2009-07-30 | 2010-12-01 | 上海交通大学 | 路网栅格数字地图自动矢量化方法 |
| DE102016205436A1 (de) * | 2015-11-25 | 2017-06-01 | Volkswagen Aktiengesellschaft | Verfahren und System zum Erstellen einer digitalen Karte |
| CN107838926B (zh) * | 2017-10-18 | 2020-11-17 | 歌尔科技有限公司 | 一种自动接送机器人 |
| CN109993813B (zh) * | 2017-12-29 | 2023-06-27 | 长城汽车股份有限公司 | 创建地图的方法、装置、车辆及可读存储介质 |
| CN108711356B (zh) * | 2018-02-08 | 2020-07-07 | 中国人民解放军61540部队 | 矢量地理pdf地图制作中地理目标与符号图形配准方法 |
| CN109192054B (zh) * | 2018-07-27 | 2020-04-28 | 阿里巴巴集团控股有限公司 | 一种地图区域合并的数据处理方法及装置 |
| CN114526753B (zh) * | 2022-04-24 | 2022-07-15 | 深圳依时货拉拉科技有限公司 | 跨路口交规的关联方法及装置、计算机设备及可读存储介质 |
| CN115797365B (zh) * | 2022-11-16 | 2025-09-26 | 武汉中海庭数据技术有限公司 | 高精度地图区域拼接单元的接边处平齐处理方法及系统 |
| CN116067389A (zh) * | 2022-12-31 | 2023-05-05 | 武汉中海庭数据技术有限公司 | 一种高精度地图规划道路修正方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1190480A (zh) * | 1994-11-04 | 1998-08-12 | 洛克系德·马丁有限公司 | 构造数据嵌拼地图的装置和方法 |
| CN1356668A (zh) * | 2001-12-30 | 2002-07-03 | 上海交通大学 | 数字道路地图自动生成滚球法 |
| JP2003135073A (ja) * | 2001-11-02 | 2003-05-13 | Inst Of Physical & Chemical Res | MAPキナーゼp38αの新規なスプライシング変異体 |
-
2004
- 2004-05-13 CN CNB2004100183066A patent/CN1299220C/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1190480A (zh) * | 1994-11-04 | 1998-08-12 | 洛克系德·马丁有限公司 | 构造数据嵌拼地图的装置和方法 |
| JP2003135073A (ja) * | 2001-11-02 | 2003-05-13 | Inst Of Physical & Chemical Res | MAPキナーゼp38αの新規なスプライシング変異体 |
| CN1356668A (zh) * | 2001-12-30 | 2002-07-03 | 上海交通大学 | 数字道路地图自动生成滚球法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1570917A (zh) | 2005-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114018239A (zh) | 一种三维车道地图构建方法、装置、设备及存储介质 | |
| CN102968634B (zh) | 一种主方向约束下的停车场结构提取方法 | |
| CN1299220C (zh) | 数字道路地图自动拼接方法 | |
| CN107885790B (zh) | 一种道路空间网络多因子自动更新方法 | |
| CN1854685A (zh) | 导航装置 | |
| CN110136159A (zh) | 面向高分辨率遥感影像的线段提取方法 | |
| CN111368016B (zh) | 多阈值约束的轨道交通控制保护区空间信息提取分析方法 | |
| CN101630463A (zh) | 路网栅格数字地图自动矢量化方法 | |
| CN109448046B (zh) | 一种基于多描述子的半自动道路中心线快速提取方法 | |
| CN110688958A (zh) | 基于GoogLeNet神经网络的交叉路口识别方法 | |
| CN116860906B (zh) | 轨迹生成方法、装置、设备及存储介质 | |
| CN111127530A (zh) | 一种多期道路三维激光点云精确配准方法 | |
| CN109101743A (zh) | 一种高精度路网模型的构建方法 | |
| CN115752432A (zh) | 无人机采集道路交通图中虚线车道线自动提取方法及系统 | |
| CN107871317A (zh) | 一种基于图像处理技术的砂浆饱满度检测方法 | |
| CN111681250B (zh) | 一种基于激光栅格地图的分割方法 | |
| JPH0561407A (ja) | 図形処理装置 | |
| CN1399223A (zh) | 房屋数字地图自动生成方法 | |
| CN111626971A (zh) | 具有图像语义感知的智慧城市cim实时成像方法 | |
| CN111238500B (zh) | 道路地图区域的道路线段的地图生成方法、装置及系统 | |
| CN1141669C (zh) | 数字道路地图自动生成滚球法 | |
| CN115439792A (zh) | 一种基于人工智能的监控方法和系统 | |
| CN116955517B (zh) | 路口区域检测方法、装置、设备及计算机可读存储介质 | |
| CN111462195B (zh) | 基于主线约束的非规则角度方向代价聚合路径确定方法 | |
| CN111737508B (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 | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070207 |