CN102944226A - 基于明暗区域配对的陨石坑检测方法 - Google Patents
基于明暗区域配对的陨石坑检测方法 Download PDFInfo
- Publication number
- CN102944226A CN102944226A CN201210508696XA CN201210508696A CN102944226A CN 102944226 A CN102944226 A CN 102944226A CN 201210508696X A CN201210508696X A CN 201210508696XA CN 201210508696 A CN201210508696 A CN 201210508696A CN 102944226 A CN102944226 A CN 102944226A
- Authority
- CN
- China
- Prior art keywords
- bright areas
- shadow region
- centerdot
- rightarrow
- vector
- 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.)
- Granted
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
基于明暗区域配对的陨石坑检测方法。它涉及图像处理技术领域,本发明要解决现有行星着陆段基于路标导航任务中的陨石坑提取方法存在的误提取率高、不易检出形状不规则的陨石坑等问题。步骤如下:基于最大稳定极值区域方法对行星着陆过程中光学相机获得的图像进行初步检测,提取图像中的阴影区域和明亮区域;删除过大和过小的区域;以检测出的陨石坑的阴影区域的矩心为中心,在半径为R的圆内搜索明亮区域,且该明亮区域与阴影区域灰度平均值的差异需大于给定阈值;计算由阴影区域矩心指向明亮区域矩心的矢量;计算该矢量与太阳光线矢量在相机像平面的投影矢量的夹角,若此夹角小于给定阈值,则该阴影区域与明亮区域构成一个陨石坑。
Description
技术领域
本发明涉及图像处理技术领域,具体涉及基于明暗区域配对的陨石坑检测方法。
背景技术
陨石坑是小天体撞击行星表面而形成的环形凹坑,具有独特的地质形态。在行星着陆任务中,可将陨石坑作为导航路标,以提高导航***的精度。在导航路标应用之前,必须将陨石坑从行星着陆器光学相机所拍摄的图像中检测并提取出来。陨石坑的检测是后续的匹配、导航等操作的基础,检测结果直接影响了导航***的性能。
在行星任务中,陨石坑的提取可用于行星着陆段的自主导航。若以行星着陆段自主导航任务为目的,不要求从图像中提取出全部的陨石坑,但要求陨石坑提取算法实时性较高,定位准确。近年来,学者们针对行星着陆这一任务目的对陨石坑检测算法进行了研究,基于边缘信息的陨石坑提取、基于模板匹配的陨石坑提取和基于区域信息的陨石坑检测是具有代表性的陨石坑提取算法。
基于边缘梯度信息和模板匹配的陨石坑检测算法均是以陨石坑具有明显椭圆外廓为假设,但行星表面上也存在着许多并不具有椭圆外廓的陨石坑,以火星表面为例,陨石坑在大风、沙尘暴等长期影响下其唇部被逐渐风蚀,导致边缘不断退化进而失去明显的几何外廓;另外受地表结构变动影响,陨石坑的形态也会出现一定的改变,导致边缘缺失等现象,而这两类算法难以检测此类退化陨石坑。基于区域信息的检测算法利用了陨石坑在光照下形成的特定区域,绕开了椭圆外廓的假设,能够检测出形状不规则的陨石坑。
最大稳定极值区域方法(MSER)是近年新兴的区域特征提取算法,MSER提取的是图像中的区域。在陨石坑提取中,我们期望提取出陨石坑内因光照所形成的明亮与阴影区域信息,并且要求提取算法具有一定的尺度及仿射不变性来应对行星着陆任务中的着陆器位姿变化,而MSER算法恰恰符合我们对陨石坑区域检测的期望,同时MSER算法简单,实时性高,也令其在行星着陆任务中具有潜在应用价值。
发明内容
本发明的目的是为了解决现有行星着陆段基于路标导航任务中的陨石坑提取方法存在的误提取率高、不易检出形状不规则的陨石坑等问题,而提出了基于明暗区域配对的陨石坑检测方法。
本发明的基于明暗区域配对的陨石坑检测方法,步骤如下:
步骤一:采用最大稳定极值区域方法对行星着陆器着陆过程中光学相机拍摄的图像进行初步检测,提取图像中的阴影区域和明亮区域;
步骤二:删除步骤一提取的阴影区域和明亮区域中K和M的区域,得到受光照所形成的n个阴影区域D1、D2…Dn和m个明亮区域L1、L2…Lm,其中,n、m为自然数;
步骤四:在步骤三搜索到的p个明亮区域中,若有j个明亮区域L1…Lj满足average(L1)-average(Dk)>σ1,…,average(Lj)-average(Dk)>σ1,则将明亮区域L1…Lj作为阴影区域Dk的待配对区域;其中,L1…Lj表示j个明亮区域,0≤j≤p,average(Lj)和average(Dk)分别表示明亮区域Lj和阴影区域Dk的图像亮度的平均值,σ1为阈值;
其中,步骤二中所述的K区域为步骤一图像中提取的小于20个像素的阴影区域和明亮区域,M区域为步骤一图像中提取的面积大于整个图像面积75%的阴影区域和明亮区域。
本发明包含以下有益效果:
本发明可用于提高陨石坑提取的正确率,特别适用于形状不规则的陨石坑。本发明与现有技术相比的优点在于:基于最大稳定极值区域方法提取陨石坑受光照所形成的阴暗区域与明亮区域,将明暗区域进行配对,生成阴影区域矩心指向明亮区域矩心的矢量,并通过该矢量与太阳光线的夹角及该矢量的大小判断明暗配对区域是否属于同一陨石坑。该方法简单易实现,实时性高,能够检测出形状不规则的陨石坑,同时具有一定的尺度及仿射不变性,尤其适用于行星着陆任务中。
本发明的原理是:基于最大稳定极值区域方法提取陨石坑受光照所形成的阴暗区域与明亮区域,将明暗区域进行配对并生成指向矢量,并通过该矢量与太阳光线的夹角及该矢量的大小判断明暗配对区域是否属于同一陨石坑。
附图说明
图1为本发明的流程图;
图2为基于明暗区域配对的陨石坑检测方法的检测结果图。
具体实施方式
结合图1说明本实施方式。
具体实施方式一:本实施方式的基于明暗区域配对的陨石坑检测方法,步骤如下:
步骤一:采用最大稳定极值区域方法对行星着陆器着陆过程中光学相机拍摄的图像进行初步检测,提取图像中的阴影区域和明亮区域;
步骤二:删除步骤一提取的阴影区域和明亮区域中K和M的区域,得到受光照所形成的n个阴影区域D1、D2…Dn和m个明亮区域L1、L2…Lm,其中,n、m为自然数;
步骤四:在步骤三搜索到的p个明亮区域中,若有j个明亮区域L1…Lj满足average(L1)-average(Dk)>σ1,…,average(Lj)-average(Dk)>σ1,则将明亮区域L1…Lj作为阴影区域Dk的待配对区域;其中,L1…Lj表示j个明亮区域,0≤j≤p,average(Lj)和average(Dk)分别表示明亮区域Lj和阴影区域Dk的图像亮度的平均值,σ1为阈值;
步骤六:根据星历确定太阳光线矢量,根据探测器自身姿态计算太阳光线矢量在相机像平面的投影矢量
其中,步骤二中所述的K区域为步骤一图像中提取的小于20个像素的阴影区域和明亮区域,M区域为步骤一图像中提取的面积大于整个图像面积75%的阴影区域和明亮区域。
本实施方式的阈值σ1和σ2根据实际情况由以往经验设定。
本实施方式步骤一中所述的最大稳定极值区域方法(MSER),提取图像中的阴影区域和明亮区域具体操作步骤如下:
2)建立连续的像素路径X:(x1,x2,…,xn),对于灰度图像来说,n=0~255;路径内像素间邻域关系按4-邻域建立,S(x)对应的为当前阈值下分割出的区域;在S(x)中寻找闭合的像素路径R(xi,….xk)1≤i<k≤n,当S内其他包含R的闭合路径R~与R均相同,即区域停止扩大时,将R定义为S的极值域;将S(x)遍历图像I(x)便得到了单幅图像的极值域,记为R(I);
3)建立极值域中的稳定区域:
在极值域R(I)中寻找灰度最大值作为阈值集:
I(R)=max(I(x))x∈R,
定义R+Δ与R-Δ分别为极值域R的扩大极小域与缩小极大域,(Δ>0):
其中argmin,argmax分别为求取泛函极小值与极大值函数;为了求取最大稳定极值区域,还需求出面积变化率:
当面积变化率δ(R,Δ)为最小时,对应的区域R就是最大稳定极值区域(MSER);
4)在单幅图像中通常存在多处MSER区域,在提取出MSER区域后,对互相有重叠的MSER区域进行像素点的逐点合并;
5)对合并后的MSER区域,若区域内像素灰度平均值低于整幅图像灰度平均值与标准差之差,则将此MSER区域标记为阴影区域,即完成对阴影区域的提取;若区域内像素灰度平均值高于整幅图像灰度平均值,则将此MSER区域标记为明亮区域。
本实施方式可用于提高陨石坑提取的正确率,特别适用于形状不规则的陨石坑。本实施方式与现有技术相比的优点在于:基于最大稳定极值区域方法提取陨石坑受光照所形成的阴暗区域与明亮区域,将明暗区域进行配对,生成阴影区域矩心指向明亮区域矩心的矢量,并通过该矢量与太阳光线的夹角及该矢量的大小判断明暗配对区域是否属于同一陨石坑。该方法简单易实现,实时性高,能够检测出形状不规则的陨石坑,同时具有一定的尺度及仿射不变性,尤其适用于行星着陆任务中。
本实施方式的原理是:基于最大稳定极值区域方法提取陨石坑受光照所形成的阴暗区域与明亮区域,将明暗区域进行配对并生成指向矢量,并通过该矢量与太阳光线的夹角及该矢量的大小判断明暗配对区域是否属于同一陨石坑。
具体实施方式二:本实施方式与具体实施方式一不同的是:步骤三中所述的
其中,|D|表示阴影区域D内所含像素点的总和。其它步骤和参数与具体是实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤七中所述的
其它步骤和参数与具体实施方式一或二相同。
通过以下试验验证本发明的有益效果:
基于明暗区域配对的陨石坑检测方法,步骤如下:
步骤一:采用最大稳定极值区域方法(MSER)对行星着陆器着陆过程中光学相机拍摄的图像进行初步检测,提取图像中的阴影区域和明亮区域;
步骤二:删除步骤一提取的阴影区域和明亮区域中K和M的区域,得到受光照所形成的n个阴影区域D1、D2…Dn和m个明亮区域L1、L2…Lm,其中,n、m为自然数;
步骤四:在步骤三搜索到的p个明亮区域中,若有j个明亮区域L1…Lj满足average(L1)-average(Dk)>σ1,…,average(Lj)-average(Dk)>σ1,则将明亮区域L1…Lj作为阴影区域Dk的待配对区域;其中,L1…Lj表示j个明亮区域,0≤j≤p,average(Lj)和average(Dk)分别表示明亮区域Lj和阴影区域Dk的图像亮度的平均值,σ1为阈值;
其中,步骤二中所述的K区域为步骤一图像中提取的小于20个像素的阴影区域和明亮区域,M区域为步骤一图像中提取的面积大于整个图像面积75%的阴影区域和明亮区域;
步骤七中所述的
本试验的阈值σ1和σ2根据实际情况由以往经验设定。
本试验的检测结果见图2所示,图2的图像为美国MOC火星观测卫星于2005年10月所拍摄,使用基于明暗区域配对的陨石坑检测方法对该幅图像中的陨石坑进行检测,图中带箭头的线段表示由阴影区域矩心指向明亮区域矩心的矢量。实线方框中阴影区域矩心指向明亮区域矩心的矢量与太阳光线方向的夹角小于给定阈值σ2(此处σ2取为30度),表示该阴影区域与明亮区域构成一个陨石坑,图中共检测出4个陨石坑。虚线方框中阴影区域矩心指向明亮区域矩心的矢量与太阳光线方向的夹角不满足小于给定阈值σ2,表示该阴影区域与明亮区域不构成陨石坑。由图2可以看出,基于明暗区域配对的陨石坑检测方法正确地检测出了4个陨石坑,没有误检测出现。
Claims (3)
1.基于明暗区域配对的陨石坑检测方法,其特征在于它的步骤如下:
步骤一:采用最大稳定极值区域方法对行星着陆器着陆过程中光学相机拍摄的图像进行初步检测,提取图像中的阴影区域和明亮区域;
步骤二:删除步骤一提取的阴影区域和明亮区域中K和M的区域,得到受光照所形成的n个阴影区域D1、D2…Dn和m个明亮区域L1、L2…Lm,其中,n、m为自然数;
步骤四:在步骤三搜索到的p个明亮区域中,若有j个明亮区域L1…Lj满足average(L1)-average(Dk)>σ1,…,average(Lj)-average(Dk)>σ1,则将明亮区域L1…Lj作为阴影区域Dk的待配对区域;其中,L1…Lj表示j个明亮区域,0≤j≤p,average(Lj)和average(Dk)分别表示明亮区域Lj和阴影区域Dk的图像亮度的平均值,σ1为阈值;
步骤六:根据星历确定太阳光线矢量,根据探测器自身姿态计算太阳光线矢量在相机像平面的投影矢量
其中,步骤二中所述的K区域为步骤一图像中提取的小于20个像素的阴影区域和明亮区域,M区域为步骤一图像中提取的面积大于整个图像面积75%的阴影区域和明亮区域。
2.根据权利要求1所述的基于明暗区域配对的陨石坑检测方法,其特征在于步骤三中所述的
其中,|D|表示阴影区域D内所含像素点的总和。
3.根据权利要求1所述的基于明暗区域配对的陨石坑检测方法,其特征在于步骤七中所述的
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210508696.XA CN102944226B (zh) | 2012-12-03 | 2012-12-03 | 基于明暗区域配对的陨石坑检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210508696.XA CN102944226B (zh) | 2012-12-03 | 2012-12-03 | 基于明暗区域配对的陨石坑检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102944226A true CN102944226A (zh) | 2013-02-27 |
CN102944226B CN102944226B (zh) | 2015-02-18 |
Family
ID=47727191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210508696.XA Expired - Fee Related CN102944226B (zh) | 2012-12-03 | 2012-12-03 | 基于明暗区域配对的陨石坑检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102944226B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700115A (zh) * | 2014-07-30 | 2015-06-10 | 北京航空航天大学 | 基于稀疏提升集成分类器的火星探测器软着陆过程中陨石坑的检测方法 |
CN107993234A (zh) * | 2016-10-26 | 2018-05-04 | 中国科学院深圳先进技术研究院 | 一种坑区域的提取方法及装置 |
CN109598243A (zh) * | 2018-12-06 | 2019-04-09 | 山东大学 | 一种月球表面安全着陆区选择方法及*** |
CN111652896A (zh) * | 2020-05-29 | 2020-09-11 | 北京理工大学 | 惯导辅助陨石坑由粗到精检测方法 |
CN111721302A (zh) * | 2020-06-24 | 2020-09-29 | 北京理工大学 | 不规则小行星表面复杂地形特征识别与感知方法 |
CN113077459A (zh) * | 2021-04-28 | 2021-07-06 | 北京的卢深视科技有限公司 | 图像清晰度检测方法、装置、电子设备及存储介质 |
CN114187529A (zh) * | 2021-12-06 | 2022-03-15 | 北京理工大学 | 小天体表面复杂地形特征检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101140624A (zh) * | 2007-10-18 | 2008-03-12 | 清华大学 | 图像匹配方法 |
CN101859382A (zh) * | 2010-06-03 | 2010-10-13 | 复旦大学 | 一种基于最大稳定极值区域的车牌检测与识别的方法 |
-
2012
- 2012-12-03 CN CN201210508696.XA patent/CN102944226B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101140624A (zh) * | 2007-10-18 | 2008-03-12 | 清华大学 | 图像匹配方法 |
CN101859382A (zh) * | 2010-06-03 | 2010-10-13 | 复旦大学 | 一种基于最大稳定极值区域的车牌检测与识别的方法 |
Non-Patent Citations (2)
Title |
---|
ERIK R.URBACH.ETC: "Automatic detection of sub-km craters in high resolution planetary images", 《PLANETARY AND SPACE SCIENCE》, 5 April 2009 (2009-04-05), pages 880 - 887 * |
冯军华等: "行星表面陨石坑检测与匹配方法", 《航空学报》, vol. 31, no. 9, 30 September 2010 (2010-09-30), pages 1858 - 1863 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700115A (zh) * | 2014-07-30 | 2015-06-10 | 北京航空航天大学 | 基于稀疏提升集成分类器的火星探测器软着陆过程中陨石坑的检测方法 |
CN104700115B (zh) * | 2014-07-30 | 2017-12-05 | 北京航空航天大学 | 基于稀疏提升集成分类器的火星探测器软着陆过程中陨石坑的检测方法 |
CN107993234A (zh) * | 2016-10-26 | 2018-05-04 | 中国科学院深圳先进技术研究院 | 一种坑区域的提取方法及装置 |
CN109598243A (zh) * | 2018-12-06 | 2019-04-09 | 山东大学 | 一种月球表面安全着陆区选择方法及*** |
CN111652896A (zh) * | 2020-05-29 | 2020-09-11 | 北京理工大学 | 惯导辅助陨石坑由粗到精检测方法 |
CN111652896B (zh) * | 2020-05-29 | 2023-06-23 | 北京理工大学 | 惯导辅助陨石坑由粗到精检测方法 |
CN111721302A (zh) * | 2020-06-24 | 2020-09-29 | 北京理工大学 | 不规则小行星表面复杂地形特征识别与感知方法 |
CN111721302B (zh) * | 2020-06-24 | 2021-11-09 | 北京理工大学 | 不规则小行星表面复杂地形特征识别与感知方法 |
CN113077459A (zh) * | 2021-04-28 | 2021-07-06 | 北京的卢深视科技有限公司 | 图像清晰度检测方法、装置、电子设备及存储介质 |
CN114187529A (zh) * | 2021-12-06 | 2022-03-15 | 北京理工大学 | 小天体表面复杂地形特征检测方法 |
CN114187529B (zh) * | 2021-12-06 | 2024-05-31 | 北京理工大学 | 小天体表面复杂地形特征检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102944226B (zh) | 2015-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102944226B (zh) | 基于明暗区域配对的陨石坑检测方法 | |
CN111209915B (zh) | 一种基于深度学习的三维图像同步识别和分割方法 | |
CN102768022B (zh) | 采用数码照相技术的隧道围岩变形检测方法 | |
CN103411621B (zh) | 一种面向室内移动机器人的光流场视觉/ins组合导航方法 | |
CN103093459B (zh) | 利用机载LiDAR点云数据辅助影像匹配的方法 | |
CN101509782B (zh) | 一种小型地标捕获定位方法 | |
CN102542551B (zh) | 极地冰盖边缘浮冰自动变化检测技术 | |
CN101750067B (zh) | 一种成像式地球敏感器地球扁率修正方法 | |
CN103593832A (zh) | 一种基于高斯二阶差分特征检测算子的图像拼接方法 | |
CN101714211B (zh) | 高分辨率遥感影像道路中心线的检测的方法 | |
CN104197933B (zh) | 一种望远镜视场中高星等滑动恒星的增强和提取方法 | |
CN102999915B (zh) | 一种基于面积比的陨石坑匹配方法 | |
CN102944238B (zh) | 一种行星探测器接近目标过程中相对位置确定方法 | |
CN105160125A (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN109376208A (zh) | 一种基于智能终端的定位方法、***、存储介质及设备 | |
CN104655135A (zh) | 一种基于地标识别的飞行器视觉导航方法 | |
CN105547287A (zh) | 一种不规则天体导航视线信息提取方法 | |
CN102663778B (zh) | 一种基于多视点视频的目标跟踪方法和*** | |
CN104766341A (zh) | 一种基于dsm和核密度估计的遥感影像建筑物变化检测方法 | |
CN103854290A (zh) | 一种结合骨架特征点和分布场描述子的扩展目标跟踪方法 | |
CN103093214B (zh) | 一种基于车载红外相机的行人检测方法 | |
CN102968795B (zh) | 基于阴影面积与外廓面积比值的陨石坑误匹配判断方法 | |
CN104809720A (zh) | 基于小交叉视场的两相机目标关联方法 | |
TWM600873U (zh) | 用於檢測道路破損之檢測系統 | |
CN114283126B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150218 Termination date: 20151203 |
|
EXPY | Termination of patent right or utility model |