CN108917749A - 一种基于偏振成像提取太阳子午线的方法 - Google Patents
一种基于偏振成像提取太阳子午线的方法 Download PDFInfo
- Publication number
- CN108917749A CN108917749A CN201810758227.0A CN201810758227A CN108917749A CN 108917749 A CN108917749 A CN 108917749A CN 201810758227 A CN201810758227 A CN 201810758227A CN 108917749 A CN108917749 A CN 108917749A
- Authority
- CN
- China
- Prior art keywords
- polarization
- straight line
- sun
- coordinate system
- meridian
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
- Eyeglasses (AREA)
Abstract
本发明提供了一种基于偏振成像提取太阳子午线的方法,属于天文导航技术领域,适用于从天空偏振模式图中提取有效的导航信息。现有太阳子午线提取方法主要针对晴朗天气做的研究,不适用于非晴朗天气,适用范围小。本发明根据太阳子午线上的偏振方位角为90度的特点,利用边缘检测和霍夫变换算法对90度偏振方位角所在直线进行检测识别,从而提取太阳子午线,该方法不易受其他因素干扰,适用于不同天气状况。
Description
技术领域
本发明属于天文导航领域,是一种基于偏振成像提取太阳子午线的方法,适用于从天空偏振模式图中提取有效的导航信息。
背景技术
仿生偏振光导航已成为目前导航技术领域研究的新热点。生物研究发现,昆虫导航是从整个(或部分)具有稳定分布的偏振模式中提取太阳子午线,再判断太阳子午线与体轴的夹角,用以确定航向。由此可知,在仿生天空偏振光导航中,太阳子午线位置信息的获取是偏振光导航应用的关键。
目前,从天空偏振模式图中提取太阳子午线的方法十分有限。田柳等利用大气偏振模式存在对称性分布的规律,通过对称性检测获得大气偏振模式对称轴位置,从而实现导航参考基准太阳子午线方向的获取。任建斌等利用模糊聚类分析的方法对偏振数据进行分析,分别就偏振度和E-矢量方向两组数据进行数据分析,确定类中心,将其作为所求取的太阳位置和太阳子午线位置。
上述两种方法主要针对在晴朗天气下获得的天空偏振模式图做的研究,在非晴朗天气下,前者偏振模式图的对称性很难检测,而后者聚类的对象受到干扰不易获取,两者都有很大的使用局限性,适用范围小。因此,提出一种能适用于从不同天气条件的天空偏振模式图中提取太阳子午线的方法。
发明内容
本发明要解决的技术问题是现有从天空偏振模式图中提取太阳子午线的方法不适用于非晴朗天气,适用范围小。本发明根据天空偏振模式图中太阳子午线上的偏振方位角为90度的特点,利用边缘检测和霍夫变换算法对90度偏振方位角所在直线进行检测识别,从而提取太阳子午线,该方法不易受其他因素干扰,适用于不同天气状况。
本发明的技术方案:
一种基于偏振成像提取太阳子午线的方法,对偏振方位角模式图进行处理,根据太阳子午线方向的偏振方位角为90度的特点,利用边缘检测和霍夫变换原理对偏振模式图中的90度偏振方位角所在直线进行检测,即可得到太阳子午线。具体步骤如下:
步骤一、利用全天空偏振模式测试系统采集数据,并进行处理得到偏振方位角模式图;
步骤二、依据90度偏振方位角所在直线为太阳子午线方向,以及偏振方位角的分布关于子午线逆对称的规律,对偏振方位角模式图进行图像二值化,设置阈值(90-a,90+a),a在1~6之间取值,偏振方位角在(90-a,90+a)区间的设置为黑色,其余值设置为白色,由此得到二值化图像;
步骤三、利用Canny算子采用双阈值法对二值化图像进行边缘检测,再利用霍夫变换识别二值化图像中的90度偏振方位角所在的直线:
1)将二值化图像的x-y直角坐标系中的直线上的每个点(x,y)转化为ρ-θ极坐标系内相应的曲线,转化关系如下:
ρ=xcosθ+ysinθ,(ρ≥0,0≤θ≤2π) (1)
其中,ρ为直角坐标系中的直线与坐标原点的距离,θ为直线的倾斜角度;
同一条直线上的点对应的曲线相交于一点,极坐标系内每一个交点都对应着原直角坐标系的一条直线;如附图中图1所示,三条曲线相交的交点所表示的ρ和θ即为原直角坐标系中三个点所在的直线位置;
2)检测极坐标系中的峰值点:峰值点为经过最多次数的曲线相交得到的交点;
3)直线连接:将中间有断开的直线线段连接成完整的长直线段,即可得到太阳子午线;
步骤四、确定太阳子午线位置:在极坐标系中可以确定峰值点的坐标(ρ1,θ1)、(ρ2,θ2)、…、(ρk,θk),峰值点被检测的次数累加为n1,n2,…,nk,则太阳子午线的位置为
θ=(n1*θ1+n2*θ2+…+nk*θk)/(n1+n2+…+nk) (2)。
本发明的有益效果:利用Canny算子进行边界检测和霍夫变换算法对90度偏振方位角所在直线进行识别,从而提取出太阳子午线,得到子午线位置参数θ,其误差角度在2度以下,图像式偏振成像与Canny算子边界检测和霍夫变换算法的结合适用于不同天气下的太阳子午线的提取,满足利用偏振模式仿生导航的需要。
附图说明
图1为霍夫变换检测原理图。
图2为测试偏振方位角模式图。
图3为太阳子午线提取流程图。
图4(a)为二值化图像;图4(b)为极坐标系图。
图5为直线检测结果图。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式:
一种基于偏振成像提取太阳子午线的方法,具体步骤如下:
步骤一、利用全天空偏振模式测试系统采集数据,并进行处理得到偏振方位角模式图;
步骤二、如图3所示,设置偏振方位角选取范围(90-a,90+a),其中a取2,则偏振方位角的选取范围是(88,92),将如图2所示的测试偏振方位角模式图导入程序,得到如图4(a)所示的二值化图像,其中黑色区域即为检测到的(88,92)偏振方位角区域,此区域关于90度偏振方位角对称。
步骤三、利用Canny算子对图4(a)进行边缘检测做预处理,再利用霍夫变换原理对此二值化图像的90度偏振方位角所在直线进行检测识别:首先将二值化图像从x-y直角坐标系转化为ρ-θ极坐标系,再对图像进行峰值点检测,检测出密集的交点,得到如图4(b)所示的极坐标系图像,图中的亮点即为峰值点的位置;最后,由极坐标系的亮点返回原来直角坐标系的边界直线,如图5所示的直线即为边界直线。
步骤四、确定太阳子午线位置:由极坐标图可知亮点的坐标为θ1=67°,θ2=64°,被检测的次数分别为1335,1313,代入公式(2)中,即可得到子午线的位置信息θ=(1335*67°+1313*64°)/(1335+1313)=65.512°。
Claims (1)
1.一种基于偏振成像提取太阳子午线的方法,其特征在于,具体步骤如下:
步骤一、利用全天空偏振模式测试系统采集数据,并进行处理得到偏振方位角模式图;
步骤二、依据90度偏振方位角所在直线为太阳子午线方向,以及偏振方位角的分布关于子午线逆对称的规律,对偏振方位角模式图进行图像二值化,设置阈值(90-a,90+a),a在1~6之间取值,偏振方位角在(90-a,90+a)区间的设置为黑色,其余值设置为白色,由此得到二值化图像;
步骤三、利用Canny算子采用双阈值法对二值化图像进行边缘检测,再利用霍夫变换识别二值化图像中的90度偏振方位角所在的直线:
1)将二值化图像的x-y直角坐标系中的直线上的每个点(x,y)转化为ρ-θ极坐标系内相应的曲线,转化关系如下:
ρ=xcosθ+ysinθ,(ρ≥0,0≤θ≤2π) (1)
其中,ρ为直角坐标系中的直线与坐标原点的距离,θ为直线的倾斜角度;
同一条直线上的点对应的曲线相交于一点,极坐标系内每一个交点都对应着原直角坐标系的一条直线;
2)检测极坐标系中的峰值点:峰值点为经过最多次数的曲线相交得到的交点;
3)直线连接:将中间有断开的直线线段连接成完整的长直线段,即可得到太阳子午线;
步骤四、确定太阳子午线位置:在极坐标系中可以确定峰值点的坐标(ρ1,θ1)、(ρ2,θ2)、…、(ρk,θk),峰值点被检测的次数累加为n1,n2,…,nk,则太阳子午线的位置为
θ=(n1*θ1+n2*θ2+…+nk*θk)/(n1+n2+…+nk) (2)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810758227.0A CN108917749B (zh) | 2018-07-11 | 2018-07-11 | 一种基于偏振成像提取太阳子午线的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810758227.0A CN108917749B (zh) | 2018-07-11 | 2018-07-11 | 一种基于偏振成像提取太阳子午线的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108917749A true CN108917749A (zh) | 2018-11-30 |
CN108917749B CN108917749B (zh) | 2021-05-11 |
Family
ID=64412273
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810758227.0A Active CN108917749B (zh) | 2018-07-11 | 2018-07-11 | 一种基于偏振成像提取太阳子午线的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108917749B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103759727A (zh) * | 2014-01-10 | 2014-04-30 | 大连理工大学 | 一种基于天空偏振光分布模式的导航定位方法 |
CN103913167A (zh) * | 2014-04-11 | 2014-07-09 | 中北大学 | 利用自然光偏振模式确定大气层内飞行器空间姿态的方法 |
CN106643704A (zh) * | 2017-01-16 | 2017-05-10 | 中国人民解放军国防科学技术大学 | 基于大气偏振模式的太阳方位获取方法 |
-
2018
- 2018-07-11 CN CN201810758227.0A patent/CN108917749B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103759727A (zh) * | 2014-01-10 | 2014-04-30 | 大连理工大学 | 一种基于天空偏振光分布模式的导航定位方法 |
CN103913167A (zh) * | 2014-04-11 | 2014-07-09 | 中北大学 | 利用自然光偏振模式确定大气层内飞行器空间姿态的方法 |
CN106643704A (zh) * | 2017-01-16 | 2017-05-10 | 中国人民解放军国防科学技术大学 | 基于大气偏振模式的太阳方位获取方法 |
Non-Patent Citations (5)
Title |
---|
CHEN FAN 等: "Integrated Polarized Skylight Sensor and MIMU With a Metric Map for Urban Ground Navigation", 《IEEE SENSORS JOURNAL》 * |
YUJIE WANG 等: "Bionic Orientation and Visual Enhancement With a Novel Polarization Camera", 《IEEE SENSORS JOURNAL》 * |
卢皓 等: "基于偏振成像的方位角度解算算法的设计与验证", 《清华大学学报 (自然科学版)》 * |
吴川: "基于大气偏振模式的导航定位方法研究", 《中国优秀硕士学位论文全文数据库》 * |
崔岩 等: "东海海域的天空光偏振及子午线研究", 《激光与光电子学进展》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108917749B (zh) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108305264B (zh) | 一种基于图像处理的无人机精确着陆方法 | |
CN104637073B (zh) | 一种基于太阳照射阴影补偿的带状地下结构探测方法 | |
CN103727930B (zh) | 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法 | |
CN107677274B (zh) | 基于双目视觉的无人机自主着陆导航信息实时解算方法 | |
CN104200461B (zh) | 基于互信息图像选块和sift特征的遥感图像配准方法 | |
CN110147095A (zh) | 基于路标信息与多传感器数据融合的机器人重定位方法 | |
CN106056619A (zh) | 基于梯度约束Radon变换的无人机视觉电线巡检方法 | |
JP2009020014A (ja) | 自己位置推定装置 | |
CN110197157B (zh) | 一种基于历史裂纹数据的路面裂纹生长检测方法 | |
CN101593351A (zh) | 基于距离变换和刚性变换参数估计的眼底图像配准方法 | |
CN110009680B (zh) | 基于圆特征及异面特征点的单目图像位置、姿态测量方法 | |
CN105307115A (zh) | 一种基于行动机器人的分布式视觉定位系统及方法 | |
CN103632384B (zh) | 组合式标记点及标记点中心的快速提取方法 | |
CN104207801A (zh) | 一种超声检测图像三维标定方法 | |
CN102095370B (zh) | 三x组合标记的检测识别方法 | |
CN105427304A (zh) | 基于多种特征联合的目标sar图像和光学图像配准方法 | |
CN105938554A (zh) | 基于图像自动判读的尖轨伸缩位移监测方法及系统 | |
CN106295491A (zh) | 车道直线检测方法及装置 | |
CN102944226A (zh) | 基于明暗区域配对的陨石坑检测方法 | |
CN104315998A (zh) | 一种基于深度图像和方位角的门开度判断方法 | |
Huang et al. | An object-based change detection approach by integrating intensity and texture differences | |
CN104457691B (zh) | 一种建筑物主体高程信息获取方法 | |
CN108875454A (zh) | 交通标志识别方法、装置和车辆 | |
CN108917749A (zh) | 一种基于偏振成像提取太阳子午线的方法 | |
CN105844614B (zh) | 一种基于校对机器人角度的视觉指北方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |