CN102930244B - 一种基于像素点可区分度的自然路标及着陆区域确定方法 - Google Patents

一种基于像素点可区分度的自然路标及着陆区域确定方法 Download PDF

Info

Publication number
CN102930244B
CN102930244B CN201210352611.3A CN201210352611A CN102930244B CN 102930244 B CN102930244 B CN 102930244B CN 201210352611 A CN201210352611 A CN 201210352611A CN 102930244 B CN102930244 B CN 102930244B
Authority
CN
China
Prior art keywords
pixel
differentiate
touchdown area
image
prime
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
CN201210352611.3A
Other languages
English (en)
Other versions
CN102930244A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201210352611.3A priority Critical patent/CN102930244B/zh
Publication of CN102930244A publication Critical patent/CN102930244A/zh
Application granted granted Critical
Publication of CN102930244B publication Critical patent/CN102930244B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

本发明涉及一种基于像素点可区分度的自然路标及着陆区域确定方法,属于深空探测邻域。该方法基于图像中平坦区域中的像素点与其邻域中的像素点的可区分度较低,而陨石坑,岩石等区域中的像素点与其邻域中的像素点的可区分度较高的特点,通过计算图像中每个像素点的可区分度,将可区分度低于给定阈值的像素点作为着陆区域的候选点,将可区分度高于给定阈值的像素点作为自然路标候选点,将给定半径内的自然路标候选点的集合作为自然路标;通过对自然路标候选点以及图像边界上的像素进行形态学膨胀,以最后被膨胀操作覆盖的着陆区域候选点作为圆心,以该圆心到距离其最近的自然路标候选点的距离为半径,得到最大可着陆圆。本发明可以快速可靠的确定可着陆区域。

Description

一种基于像素点可区分度的自然路标及着陆区域确定方法
技术领域
本发明属于深空探测邻域,涉及到图像处理的相关技术。
背景技术
自然路标检测和着陆区域确定是深空探测领域中的重要问题。在深空探测过程中,由于着陆器距离地球非常遥远,信号传递具有较大的延时,因此需要着陆器能够自主的确定自身的姿态和位置以及自主的找到可着陆区域。
自然路标是指环境中已有的,并非人工设置的,能够用以标识不同环境的参考对象或模式。自然路标检测对于着陆器确定自身的姿态和位置,进行自主导航有着重要的作用。现有的自然路标检测方法多是针对行星表面的各种地形(如陨石坑,岩石等,如图1所示)的图像表观特征进行自然路标检测。但是这类方法受到光照以及具体地形的影响,无法取得很好的结果。
而自动的着陆区域确定对于探测器来说也是非常重要。在着陆器着陆阶段,由于距离地球非常遥远,依靠地面通讯不能满足实时安全着陆的需求,因此需要探测器具有较强的自主确定着陆区域的能力。现有的自主着陆区域确定方法主要是通过检测行星表面的陨石坑,岩石等障碍物,通过避开障碍物来确定着陆区域。也有根据行星表面的纹理特征进行着陆区域的确定方法。通过检测障碍物进行着陆区域确定的方法的问题是障碍物的检测很难做到精确和鲁棒,同时受光照条件影响较大。现有的纹理分析方法是通过统计区域内的对比度等特征来进行,也易受光照等条件的影响。
我们通过观察发现,行星表面的障碍物,不管是陨石坑(图1(a)),岩石(图1(b)),或者裂纹(图1(c)),都具有以下特点:在这些障碍物边界上的点(如图1中的1,2,3,4),与其邻域内的其他点的可区分程度很高,既很容易将这些点与周围其他的点区分开来。而平坦区域内的点(如图1中的5,6,7),与其周围的点比较相似,基本无法与周围点区分开来。我们基于这一观察,将自然路标检测与着陆区域确定统一在像素点的可区分度框架下进行。
发明内容
本发明的目的是针对上述自然路标检测与着陆区域确定方法存在的不足而提出一种基于像素点可区分度的自然路标检测及着陆区域确定方法。
本发明的主要内容为:采用摄像机在着陆器下降过程中获取着陆器下方的图像信息,在每一帧获取到的图像中,计算每个像素点的可区分度,既该像素点与其邻域内的像素点的可区分程度。设定阈值,可区分度高于给定阈值的像素点作为自然路标候选点,将一定半径内的自然路标候选点的集合作为自然路标;可区分度低于给定阈值的像素点作为可安全着陆区域的候选点,通过对自然路标候选点及图像边界上的像素点进行形态学膨胀操作,以最后被膨胀操作覆盖的着陆区域候选点作为圆心,以该圆心到距离其最近的自然路标候选点的距离为半径,得到最大可着陆圆。
本发明的目的是通过下述技术方案实现的。
一种基于像素点可区分度的自然路标检测及着陆区域确定方法,具体实现步骤如下:
步骤一:图像获取
摄像机安装在着陆器的下方,在着陆器下降过程中可以连续的获得地面的图像。
步骤二:计算像素点的可区分度
设s(p,p′)为像素p与像素p′的相似程度,设p(s(p,p′))表示像素p与像素p′相似的概率,定义像素的可区分度为该像素与其邻域内像素相似概率分布的熵,既
distinguishbility p = - Σ p ′ ∈ N p p ( s ( p , p ′ ) ) log p ( s ( p , p ′ ) ) - - - ( 1 )
其中Np为像素p的邻域。熵越大,则可区分度越低。
步骤三:基于像素点的可区分度确定可着陆区域
将可区分度高和可区分度中等的像素所覆盖的区域都看作是障碍物区域,可区分度低的像素所覆盖的区域看作是可着陆区域。将图像边界上的点也看作是可区分度高的像素点,既也看作障碍区域。每次对障碍物区域的像素进行形态学膨胀,将障碍物周围的像素也设置为障碍物区域,随着膨胀的进行,可着陆区域越来越小,最后被膨胀操作覆盖的可着陆区域就是最大可着陆区域的圆心。最大可着陆圆的半径可通过计算与该点距离最近的膨胀之前的高可区分度或中等区分度的像素点的距离得到
步骤四:基于高区分度像素的自然路标检测
具有高可区分度的像素能够与其周围的像素区分开来,可以作为自然路标。我们将一定半径内的具有高可区分度的像素的集合作为自然路标。
本发明的优点
本发明与其它自然路标及着陆区域确定方法相比,具有以下几个方面的优势:
(1)将陨石坑,岩石等各种自然路标的检测以及可着陆区域的检测统一在可区分度的框架下进行处理,可同时得到自然路标和着陆区域。
(2)通过对非自然路标候选点进行形态学的膨胀操作来得到最大可着陆圆的圆心及半径,可以快速可靠的确定可着陆区域。
附图说明
图1行星表面各种地形;
图2基于像素点可区分度的自然路标及着陆区域确定方法流程图;
图3像素点与其邻域内像素的相似概率分布图;
图4最大可着陆圆确定方法示意图。
具体实施方式
下面结合附图和实施例对本发明做详细说明。
本发明提出的基于像素点可区分度的自然路标检测与着陆区域确定方法的流程图如图2所示,具体实施步骤如下:
步骤一:图像获取
摄像机安装在着陆器的下方,在着陆器下降过程中可以连续的获得地面的图像。
步骤二:计算像素点的可区分度
设p表示图像中的一个像素,其坐标为(px,py),其像素值为I(px,py),则像素点p的可区分度为p与其邻域中的其他像素的可区分程度。设像素p的邻域为Np,Np为以像素p为中心的一个矩形,其长宽可根据实际情况进行设定,如果其长宽设定为图像的宽度和高度,则是在整副图像中进行该像素的可区分度,设置邻域Np的长度和宽度分别为图像宽度和高度的三分之一。定义s(p,p′)为像素p与像素p′的相似程度,此处,我们采用SSD来计算两个像素的相似程度,SSD值越小,表明两个像素越相似,设计算SSD时采用的邻域大小为(2n+1)×(2m+1),则
s ( p , p ′ ) = Σ i = - n n Σ j = - m m [ I ( px + i , py + j ) - I ( p ′ x + i , p ′ y + j ) ] 2 - - - ( 2 )
其中(px,py)为像素p在图像中的坐标,(p′x,p′y)为像素p′在图像中的坐标。设p(s(p,p′))表示像素p与像素p′相似的概率,则
p ( s ( p , p ′ ) ) = exp { - s ( p , p ′ ) } Σ p ′ ∈ N p exp { - s ( p , p ′ ) } - - - ( 3 )
定义像素的可区分度为该像素与其邻域内像素相似概率分布的熵,既
distinguishbility p = - Σ p ′ ∈ N p p ( s ( p , p ′ ) ) log p ( s ( p , p ′ ) ) - - - ( 4 )
对于图1中标记的像素点计算与其邻域内像素的相似概率分布,如图3所示,可以看出,分布有三种情况:
1.点状分布,如图3(1-2)所示,这些像素一般是图像中的角点,其相似概率分布的熵较小,可区分度高。
2.线状分布,如图3(3-4)所示,这些像素一般是图像中的边缘,其相似概率分布的熵居中,可区分度中等。
3.面状分布,如图3(5-7)所示,这些像素一般属于平坦区域,其相似概率分布的熵较大,可区分度低。
步骤三:基于像素点的可区分度确定可着陆区域
我们将可区分度高和可区分度中等的像素所覆盖的区域都看作是障碍物区域,可区分度低的像素所覆盖的区域看作是可着陆区域。可着陆区域应满足以下两个要求,一是面积足够大,以便着陆器着陆;二是距离障碍物足够远,以确保着陆过程的安全。由于图像边界外面的区域是未知区域,因此理想的着陆区域应该距离图像边界也足够远。因此我们将图像边界上的点的可区分度也设置为高,既边界上的点也看作是可区分度高的像素点。如图4所示,图中红色的点障碍物区域,边界上的点也认为是障碍物区域。绿色部分为可区分度低的像素,可作为着陆区域,我们的目标是找到一个最大的圆,使得该圆仅包含绿色区域作为可着陆区域。具体做法是每次对障碍物区域的像素进行形态学膨胀,将障碍物周围的像素也设置为障碍物区域,如图4(b),(c),(d)所示,可以看到,随着膨胀的进行,绿色区域越来越小,最后变红的绿色区域就是最大可着陆区域的圆心。具体流程为:
1)依次对障碍物区域中的像素(高可区分度和中等区分度的像素)点进行形态学膨胀,每次膨胀一个像素;既将障碍物区域周围一个像素的区域也设置为障碍物区域。
2)判断可区分度低的像素的个数,如果数目大于1,则转第一步;如果数目等于1,则转第三步,如果数目为0,则转第四步。
3)如果可区分度低的像素点的数目等于1,则该像素点是最大可着陆圆的圆心,最大可着陆圆的半径可通过计算与该点距离最近的膨胀之前的高可区分度或中等区分度的像素点的距离得到,如图4(d)所示。
4)如果可区分度低的像素点的数目等于0,则选取此次膨胀之前,未被覆盖的可区分度低的像素点作为可着陆区域的圆心,如有多个这样的像素点,则取最接近图像中心的一个作为可着陆区域的圆心,半径可通过计算与该点距离最近的膨胀之前的高可区分度或中等区分度的像素点的距离得到。
步骤四:基于高区分度像素的自然路标检测
自然路标是指环境中已有的,并非人工设置的,能够用以标识不同环境的参考对象或模式。由自然路标的定义可以发现,具有高可区分度的像素符合自然路标的定义,既具有高可区分度的像素能够与其周围的像素区分开来,可以作为自然路标。我们将一定半径内的具有高可区分度的像素的集合作为自然路标。由于这些点的可区分度很高,因此在后续的图像中可以比较有保证的跟踪到这些像素点。

Claims (1)

1.一种基于像素点可区分度的自然路标及着陆区域确定方法,其特征在于:一种基于像素点可区分度的自然路标检测及着陆区域确定方法,具体实现步骤如下:
步骤一:图像获取
摄像机安装在着陆器的下方,在着陆器下降过程中连续的获得地面的图像;
步骤二:计算步骤一获取的图像中的像素点的可区分度
设p表示图像中的一个像素点,其坐标为(px,py),其像素值为I(px,py),则像素点p的可区分度为像素点p与其邻域中的其他像素点的可区分程度;设像素点p的邻域为Np,Np为以像素点p为中心的一个矩形,其长宽可根据实际情况进行设定,如果其长宽设定为图像的宽度和高度,则在整副图像中对该像素点进行可区分度区分,设置邻域Np的长度和宽度分别为图像宽度和高度的三分之一,定义s(p,p')为像素点p与像素点p'的相似程度,采用SSD来计算两个像素的相似程度,SSD值越小,表明两个像素越相似,设计算SSD时采用的邻域大小为(2n+1)×(2m+1),则
s ( p , p ′ ) = Σ i = - n n Σ j = - m m [ I ( p x + i , p y + j ) - I ( p ′ x + i , p ′ y + j ) ] 2 - - - ( 2 )
其中(px,py)为像素点p在图像中的坐标,(p'x,p'y)为像素点p'在图像中的坐标,设p(s(p,p'))表示像素点p与像素点p'相似的概率,则
p ( s ( p , p ′ ) ) = exp { - s ( p , p ′ ) } Σ p ′ ∈ N p exp { - s ( p , p ′ ) } - - - ( 3 )
定义像素点的可区分度为该像素点与其邻域内像素点相似概率分布的熵,既
distinguishbility p = - Σ p ′ ∈ N p p ( s ( p , p ′ ) ) log p ( s ( p , p ′ ) ) - - - ( 4 )
其中,熵越大,像素点的可区分度越低,因之得出以下关于像素点可区分度的三种分布状态:
(1)点状分布,这些像素点是图像中的角点,其相似概率分布的熵较小,为可区分度高像素点;
(2)线状分布,这些像素点是图像中的边缘,其相似概率分布的熵居中,为可区分度中等像素点;
(3)面状分布,这些像素点是图像中平坦区域,其相似概率分布的熵较大,为可区分度低像素点;
步骤三:基于步骤二计算的像素点的可区分度确定可着陆区域
将可区分度高和可区分度中等的像素点所覆盖的区域都看作是障碍物区域,可区分度低的像素点所覆盖的区域看作是可着陆区域;将图像边界上的点也看作是可区分度高的像素点,既也看作障碍物区域;每次对障碍物区域的像素点进行形态学膨胀,将障碍物周围的像素点也设置为障碍物区域,随着膨胀的进行,可着陆区域越来越小,最后被膨胀操作覆盖的可着陆区域就是最大可着陆区域的圆心;最大可着陆圆的半径通过计算与圆心距离最近的膨胀之前的可区分度高或可区分度中等的像素点的距离得到,其最大可着陆圆的圆心及半径的确定步骤为:
1)依次对障碍物区域中的像素点即可区分度高和可区分度中等的像素点进行形态学膨胀,每次膨胀一个像素点;既将障碍物区域周围一个像素点的区域也设置为障碍物区域;
2)判断可区分度低的像素点的个数,如果数目大于1,则转第一步;如果数目等于1,则转第三步,如果数目为0,则转第四步;
3)如果可区分度低的像素点的数目等于1,则该像素点是最大可着陆圆的圆心,最大可着陆圆的半径通过计算与该像素点距离最近的膨胀之前的可区分度高或可区分度中等的像素点的距离得到;
4)如果可区分度低的像素点的数目等于0,则选取此次膨胀之前,未被覆盖的可区分度低的像素点作为可着陆区域的圆心,如有多个这样的像素点,则取最接近图像中心的一个作为可着陆区域的圆心,半径可通过计算与该像素点距离最近的膨胀之前的高可区分度或可区分度中等的像素点的距离得到;
步骤四:基于步骤二得到的高区分度像素点进行自然路标检测
具有高可区分度的像素点能够与其周围的像素点区分开来,符合作为路标的要求,作为自然路标;将一定半径内的具有高可区分度的像素点集合作为自然路标。
CN201210352611.3A 2012-09-20 2012-09-20 一种基于像素点可区分度的自然路标及着陆区域确定方法 Expired - Fee Related CN102930244B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210352611.3A CN102930244B (zh) 2012-09-20 2012-09-20 一种基于像素点可区分度的自然路标及着陆区域确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210352611.3A CN102930244B (zh) 2012-09-20 2012-09-20 一种基于像素点可区分度的自然路标及着陆区域确定方法

Publications (2)

Publication Number Publication Date
CN102930244A CN102930244A (zh) 2013-02-13
CN102930244B true CN102930244B (zh) 2016-04-06

Family

ID=47645041

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210352611.3A Expired - Fee Related CN102930244B (zh) 2012-09-20 2012-09-20 一种基于像素点可区分度的自然路标及着陆区域确定方法

Country Status (1)

Country Link
CN (1) CN102930244B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101539422A (zh) * 2009-04-22 2009-09-23 北京航空航天大学 一种单目视觉实时测距方法
CN102173313A (zh) * 2010-12-24 2011-09-07 北京控制工程研究所 一种软着陆接力避障方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101539422A (zh) * 2009-04-22 2009-09-23 北京航空航天大学 一种单目视觉实时测距方法
CN102173313A (zh) * 2010-12-24 2011-09-07 北京控制工程研究所 一种软着陆接力避障方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种基于形态学的探测器安全着陆点选择方法;邵巍等;《宇航学报》;20080331;第29卷(第2期);659-663 *

Also Published As

Publication number Publication date
CN102930244A (zh) 2013-02-13

Similar Documents

Publication Publication Date Title
US11288521B2 (en) Automated road edge boundary detection
CN104637073B (zh) 一种基于太阳照射阴影补偿的带状地下结构探测方法
CN111551958A (zh) 一种面向矿区无人驾驶的高精地图制作方法
CN104700399A (zh) 一种基于高分遥感影像的大变形滑坡位移场标定方法
CN102542551B (zh) 极地冰盖边缘浮冰自动变化检测技术
CN104850834A (zh) 基于三维激光雷达的道路边界检测方法
US20220363263A1 (en) Automated bump and/or depression detection in a roadway
CN104536009A (zh) 一种激光红外复合的地面建筑物识别及导航方法
US9239960B2 (en) Three-dimensional object detection device
CN102081801A (zh) 多特征自适应融合船舶跟踪和航迹检测方法
CN115774444B (zh) 一种基于稀疏导航地图的路径规划优化方法
CN104134209A (zh) 一种视觉导航中的特征提取与匹配方法及***
CN104166163A (zh) 基于三维大数据量地震数据体的断层曲面自动提取方法
CN104574365A (zh) 障碍物检测装置及方法
CN109241855B (zh) 基于立体视觉的智能车辆可行驶区域探测方法
CN115728781A (zh) 一种基于激光雷达点云的小障碍物检测方法及装置
Lee et al. Creation of river terrain data using region growing method based on point cloud data from UAV photography
CN111721302A (zh) 不规则小行星表面复杂地形特征识别与感知方法
CN103456027B (zh) 一种机场空间关系约束下的时敏目标检测定位方法
CN114332658B (zh) 基于无人机巡检的铁路工务设备与周边环境隐患排查方法
CN102930244B (zh) 一种基于像素点可区分度的自然路标及着陆区域确定方法
Fang et al. Extraction 3D road boundaries from mobile laser scanning point clouds
CN104700428B (zh) 一种位于无人机上的沙漠区域道路检测方法
CN115797900B (zh) 基于单目视觉的车路姿态感知方法
Gong et al. Roof-cut guided localization for building change detection from imagery and footprint map

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160406

Termination date: 20160920

CF01 Termination of patent right due to non-payment of annual fee