CN106441310B - 一种基于cmos的太阳方位角计算方法 - Google Patents

一种基于cmos的太阳方位角计算方法 Download PDF

Info

Publication number
CN106441310B
CN106441310B CN201611085818.3A CN201611085818A CN106441310B CN 106441310 B CN106441310 B CN 106441310B CN 201611085818 A CN201611085818 A CN 201611085818A CN 106441310 B CN106441310 B CN 106441310B
Authority
CN
China
Prior art keywords
polarization
cmos
azimuth
coordinate system
channel
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.)
Active
Application number
CN201611085818.3A
Other languages
English (en)
Other versions
CN106441310A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201611085818.3A priority Critical patent/CN106441310B/zh
Publication of CN106441310A publication Critical patent/CN106441310A/zh
Application granted granted Critical
Publication of CN106441310B publication Critical patent/CN106441310B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00Measuring angles

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Solid State Image Pick-Up Elements (AREA)

Abstract

本发明涉及一种基于CMOS的太阳方位角计算方法,首先,利用三通道CMOS偏振传感器获取各通道偏振图像测量信息,并对偏振图像做图像预处理;其次,利用在实验室内标定好的参数对图像预处理过的偏振图像测量信息补偿,获取偏振方位角;最后,充分利用基于CMOS的偏振传感器采样点多的优点,利用基于原始偏振图像信息得到偏振方位角的分布区间,选择偏振方位角计算式,获取太阳方位角。本发明具有精度高、计算量小的优点,可用于基于CMOS的偏振导航传感器的太阳方位角计算方法。

Description

一种基于CMOS的太阳方位角计算方法
技术领域
本发明涉及一种基于CMOS的太阳方位角计算方法,可用于提升无人机或地面车辆用偏振光导航系统的太阳方位角解算精度。
背景技术
大气偏振模式蕴含着太阳方位信息,可用于导航目的。本发明提出一种基于CMOS的太阳方位角计算方法,通过测量天空偏振光的方向和强度,得到太阳在模块坐标系(m系)下的方位角信息,结合天文年历信息,可以进一步获得载体的航向信息。
传统偏振导航方法受所用传感器器性能局限,未能充分利用大气偏振信息。一种基于六通道光电传感器的偏振方位角确定方法,专利号:201410046227,受限于光电式传感器空间分辨率低,所得偏振方位角精度有待提高。利用天空偏振模式分布规律计算导航方向角的方法,专利号:201010541166.6,提出一种基于偏振分布模式的导航定向方法,利用多个采样点的偏振信息得到太阳子午线所在位置,从而得到导航所需方向角信息,采样点多易受外界环境影响,且在线计算困难。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种基于CMOS的太阳方位角计算方法,利用载体上携带的偏振传感器得到天空偏振图像并进行图像预处理,然后利用在实验室内标定好的参数对图像预处理过的偏振图像测量信息补偿,获取偏振方位角。最后,充分利用基于CMOS的偏振导航传感器采样点多的优点,得到太阳方位角。
本发明的技术解决方案为:一种基于CMOS的太阳方位角计算方法,其特征实现步骤如下:
(1)利用三通道CMOS偏振传感器同时获取偏振图像数据Img1、Img2和Img3,采用中值滤波对图像Img1、Img2和Img3预处理,然后将三幅处理后的图像以图像中心为原点裁剪为10×10像素的方块,最后分别将三个10×10的像素矩阵求均值得到像素灰度均值S1、S2和S3
(2)利用在室内标定好的CMOS偏振传感器参数对像素灰度均值S1、S2和S3进行补偿,求取模块坐标系(m系)下的偏振方位角φ';
(3)由步骤(1)所得的原始偏振图像信息Img1,对图像Img1依据模块坐标系(m系)分割图像,根据图像光强信息判断太阳在模块坐标系(m系)中的区间位置,由步骤(2)所得偏振方位角φ'得到模块坐标系(m系)下的太阳方位角φ;
所述步骤(2)具体实现如下:
首先,三通道CMOS偏振传感器每个通道像素灰度均值可表示为:
Sn=K1nI(1+K2ndcos2(φ'-φn-Δφn))
其中,n∈{1,2,3}表示第n个通道,Sn表示第n个通道的像素灰度均值,Sn∈{S1,S2,S3},I表示相对光强功率,d表示偏振度,φ'表示偏振方位角,φn为第n个通道的主检偏方向与模块坐标系(m系)x轴正方向夹角,且K1n是第n个通道的CMOS器件电气性能有关的常数,K2n是第n个通道的与偏振片消光比有关的常数,Δφn是第n个通道的由于偏振片实际检偏方向与预设方向的小角度误差;
然后,求取偏振方位角φ':
将标定好的偏振传感器参数K1n、K2n和Δφn带入像素灰度均值S1、S2和S3中,可得:
对P1和P2求解可得偏振方位角φ'计算公式:
其中,P1、P2、l1、l2、n1、n2均为中间变量,P1是像素灰度均值S2与S1的比值,P2是像素灰度均值S3与S1的比值,l1、l2、n1、n2满足下式:
所述步骤(3)具体实现如下:
如图2所示,偏振传感器水平安装在载体上(图中点o处),镜头垂直对准天顶点t,以偏振传感器o点为原点建立模块坐标系M,模块坐标系x轴与y轴所在平面为传感器安装平面,且x轴与图3中通道1中实线所表示的主检偏方向相同,z轴方向服从右手定则指向天顶点t。将模块坐标系原点平移至原始偏振图像Img1中心,设原始偏振图像Img1在x轴上方的像素灰度值和为I1,x轴下方的像素灰度值和为I2。如果I1>I2,说明太阳在图像中位于象限1和2中,如果I1<I2,说明太阳在图像中位于象限3或4中。
通过判断x轴上方像素灰度值和和x轴下方像素灰度值和的大小关系,求取模块坐标系(m系)下太阳方位角φ:
其中,φ表示模块坐标系(m系)下的太阳方位角,φ∈(-π,π]。
本发明与现有技术相比的优点在于:
(1)本发明提出了一种基于CMOS的太阳方位角方法,通过偏振传感器得到模块系下偏振图像,充分利用CMOS图像传感器采样点多的优点,得到模块坐标系(m系)下的太阳方位角,为载体提供快速、精确的导航信息。
(2)本发明采用的基于CMOS的偏振传感器为自主导航传感器,不受外界干扰,可以在复杂环境下实现无源、无辐射导航工作,且误差不随时间的积累,在应用过程中具有隐蔽性和自主性。
附图说明
图1为本发明的设计流程图;
图2为本发明涉及偏振方位角示意图;
图3为本发明涉及三通道CMOS偏振导航传感器结构示意图。
具体实施方式
如图1所示,本发明的具体实现步骤如下:
1、利用三通道CMOS偏振传感器获取三幅原始偏振图像Img1、Img2和Img3,然后对原始偏振图像进行图像预处理,分别得到三幅图像的中心裁剪所得区域的像素灰度均值S1、S2和S3
如图2所示,偏振传感器水平安装在载体上,镜头沿着z轴垂直对准天顶点t,以偏振传感器中心o点为原点建立模块坐标系M,模块坐标系x轴与y轴所在水平面为传感器安装平面,且x轴与图3中通道1中实线所表示的主检偏方向相同,z轴方向服从右手定则指向天顶点t。点s表示太阳在水平面xoy中的投影点,有EE'过原点o,且EE'⊥so,则∠xoE即是偏振方位角φ',取值范围∠sox即是太阳方位角φ,取值范围(-π,π]。
对CMOS偏振传感器同时获取偏振图像数据Img1、Img2和Img3做图像预处理,首先将三幅图分别做中值滤波,然后将三幅处理后的图像以图像中心为原点裁剪为10×10像素的方块,最后分别将三个10×10的像素矩阵求均值得到像素灰度均值S1、S2和S3,本实施例中,选用像素10位的CMOS偏振传感器,取S1=338.87,S2=298.75和S3=322.17;
2、利用在室内标定好的CMOS偏振传感器参数对各通道像素灰度均值S1、S2和S3进行补偿,求取模块坐标系(m系)下的偏振方位角φ',具体实现如下:
(1)三通道CMOS偏振传感器各通道像素灰度均值可表示为:
Sn=K1nI(1+K2ndcos2(φ'-φn-Δφn))
其中,在本实施例中,n∈{1,2,3}表示第n个通道,Sn表示第n个通道的像素灰度均值,I表示相对光强功率且I=1,d表示偏振度,取值范围d∈[0,1],φ'表示偏振方位角,φn为第n个通道的主检偏方向与模块坐标系(m系)x轴正方向夹角,取φ1=0,K1n是CMOS器件电气性能有关的常数,一般取K11=1,K12∈[0.7,1.3],K13∈[0.7,1.3],K2n是与偏振片消光比有关的常数,一般取K21=1,K22∈[0.85,1.1],K23∈[0.85,1.1],Δφn是由于偏振片实际检偏方向与预设方向的小角度误差,与安装精度有关,取Δφn∈[-6°,6°],K1n、K2n和Δφn由于是器件误差,具体数值以实际标定结果为准。在本实例中,取K1n为K11=1,K12=0.9100,K13=1.0600,取K2n为K21=1,K22=1.0085,K23=1.0057,Δφ1=0°,Δφ2=0.3760°,Δφ3=1.7674°。
(2)求取偏振方位角φ':
将标定好的偏振传感器参数K1n、K2n和Δφn带入各通道像素灰度均值S1、S2和S3中,可得:
其中,P1、P2、l1、l2、n1、n2均为中间变量,P1是像素灰度均值S2与S1的比值,P2是像素灰度均值S3与S1的比值,l1、l2、n1、n2满足下式:
带入后,可得偏振方位角φ':
3、根据通道1图像光强信息,求取模块坐标系(m系)下太阳方位角,具体实现如下:
如图3所示,为三通道CMOS偏振导航传感器结构示意图,这里将步骤1中所述模块坐标系原点o平移至图中通道一中心处,图中实线表示为各个通道预设的主检偏方向,各主检偏方向与x轴正方向夹角即为φn,图中虚线表示实际偏振片安装方向,各通道实线与虚线的夹角即安装方向与各主检偏方向夹角为Δφn。设原始偏振图像Img1在x轴上方的像素灰度值和为I1,x轴下方的像素灰度值和为I2。如果I1>I2,说明太阳在图像中位于象限1和2中,如果I1<I2,说明太阳在图像中位于象限3或4中。本实例中,由于I1=335.7,I2=476.5,所以可以根据下式:
得到模块坐标系下太阳方位角φ=-68.298°。

Claims (2)

1.一种基于CMOS的太阳方位角计算方法,其特征实现步骤如下:
(1)利用三通道CMOS偏振传感器同时获取偏振图像数据Img1、Img2和Img3,采用中值滤波对图像Img1、Img2和Img3预处理,然后将三幅处理后的图像以图像中心为原点裁剪为10×10像素的方块,最后分别将三个10×10的像素矩阵求均值得到像素灰度均值S1、S2和S3
(2)利用在室内标定好的CMOS偏振传感器参数对像素灰度均值S1、S2和S3进行补偿,求取模块坐标系下的偏振方位角φ';所述偏振方位角的求取方法如下:
1)三通道CMOS偏振传感器每个通道预处理结果可表示为:
Sn=K1nI(1+K2ndcos2(φ'-φn-Δφn))
其中,n∈{1,2,3}表示第n个通道,Sn表示第n个通道的像素灰度均值,Sn∈{S1,S2,S3},I表示相对光强功率,d表示偏振度,φ'表示偏振方位角,φn为第n个通道的主检偏方向与模块坐标系x轴正方向夹角,且K1n是第n个通道CMOS器件电气性能有关的常数,K2n是第n个通道的与偏振片消光比有关的常数,Δφn是第n个通道由于偏振片实际检偏方向与预设方向的小角度误差;
2)求取模块坐标系下偏振方位角φ':
将室内标定好的CMOS偏振传感器参数K1n、K2n和Δφn带入像素灰度均值S1、S2和S3中,可得:
对P1和P2求解可得偏振方位角φ'计算公式:
其中,P1、P2、l1、l2、n1、n2均为中间变量,P1是像素灰度均值S2与S1的比值,P2是像素灰度均值S3与S1的比值,l1、l2、n1、n2满足下式:
(3)由步骤(1)所得的原始偏振图像信息Img1,对图像Img1依据模块坐标系分割图像,根据图像光强信息判断太阳在模块坐标系中的区间位置,由步骤(2)所得偏振方位角φ'得到模块坐标系下的太阳方位角φ。
2.根据权利要求1所述的基于CMOS的太阳方位角计算方法,其特征在于:所述步骤(3)具体实现如下:
选择CMOS偏振传感器的第一个通道的预设检偏方向为参考方向,CMOS偏振传感器水平安装在载体上,镜头垂直对准天顶点t,以CMOS偏振传感器o点为原点建立模块坐标系M,模块坐标系x轴与y轴所在平面为CMOS偏振传感器安装平面,x轴与参考方向相同,z轴方向服从右手定则指向天顶点t,将模块坐标系原点平移至原始偏振图像Img1中心,设原始偏振图像Img1在x轴上方的像素灰度值和为I1,x轴下方的像素灰度值和为I2;如果I1>I2,说明太阳在图像中位于象限1和2中,如果I1<I2,说明太阳在图像中位于象限3或4中;
通过判断x轴上方像素灰度值和和x轴下方像素灰度值和的大小关系,求取模块坐标系下太阳方位角φ:
其中,φ表示模块坐标系下太阳方位角,φ∈(-π,π]。
CN201611085818.3A 2016-11-30 2016-11-30 一种基于cmos的太阳方位角计算方法 Active CN106441310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611085818.3A CN106441310B (zh) 2016-11-30 2016-11-30 一种基于cmos的太阳方位角计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611085818.3A CN106441310B (zh) 2016-11-30 2016-11-30 一种基于cmos的太阳方位角计算方法

Publications (2)

Publication Number Publication Date
CN106441310A CN106441310A (zh) 2017-02-22
CN106441310B true CN106441310B (zh) 2019-06-04

Family

ID=58222806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611085818.3A Active CN106441310B (zh) 2016-11-30 2016-11-30 一种基于cmos的太阳方位角计算方法

Country Status (1)

Country Link
CN (1) CN106441310B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020241130A1 (ja) * 2019-05-29 2020-12-03 古野電気株式会社 情報処理システム、方法、及びプログラム
CN111121781A (zh) * 2019-12-18 2020-05-08 天津大学 基于地基云图的太阳定位和太阳光干扰判别方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102538783A (zh) * 2012-02-14 2012-07-04 北京大学 基于遥感天空偏振模式图的仿生导航方法及导航定位系统
CN103759725A (zh) * 2014-02-10 2014-04-30 北京航空航天大学 一种基于六通道光电传感器的偏振方位角确定方法
CN104359454A (zh) * 2014-11-17 2015-02-18 中北大学 基于大气偏振光的太阳空间位置获取方法
CN105157667A (zh) * 2015-06-04 2015-12-16 北京航空航天大学 一种基于大气偏振信息的太阳高度角计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8390696B2 (en) * 2009-01-06 2013-03-05 Panasonic Corporation Apparatus for detecting direction of image pickup device and moving body comprising same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102538783A (zh) * 2012-02-14 2012-07-04 北京大学 基于遥感天空偏振模式图的仿生导航方法及导航定位系统
CN103759725A (zh) * 2014-02-10 2014-04-30 北京航空航天大学 一种基于六通道光电传感器的偏振方位角确定方法
CN104359454A (zh) * 2014-11-17 2015-02-18 中北大学 基于大气偏振光的太阳空间位置获取方法
CN105157667A (zh) * 2015-06-04 2015-12-16 北京航空航天大学 一种基于大气偏振信息的太阳高度角计算方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于Rayleigh大气偏振模式的太阳空间位置优化计算;王飞等;《光子学报》;20141231;第43卷(第12期);正文第1-3节
基于偏振成像的方位角度解算算法的设计与验证;卢皓等;《清华大学学报(自然科学版)》;20141231;第54卷(第11期);全文
基于多偏振敏感器的太阳矢量测量方法;杨中光等;《光子学报》;20130131;第42卷(第1期);正文第1-2节
基于大气偏振光学的太阳跟踪方法;王晨光等;《科学技术与工程》;20161031;第16卷(第29期);正文第1-2节

Also Published As

Publication number Publication date
CN106441310A (zh) 2017-02-22

Similar Documents

Publication Publication Date Title
CN108613628B (zh) 一种基于双目视觉的架空输电线路弧垂测量方法
JP5357902B2 (ja) 撮像装置向き検出装置および当該装置を備える移動体
CN110842940A (zh) 一种建筑测量机器人多传感器融合三维建模方法及系统
CN105716542B (zh) 一种基于柔性特征点的三维数据拼接方法
CN111412916B (zh) 一种基于大气偏振光场的天文航海船位计算方法
CN108868268A (zh) 基于点到面距离和互相关熵配准的无人车位姿估计方法
CN106127787B (zh) 一种基于逆投影变换的相机标定方法
Stürzl et al. A fisheye camera system for polarisation detection on UAVs
CN104200086A (zh) 宽基线可见光相机位姿估计方法
CN104794743A (zh) 一种车载激光移动测量系统彩色点云生产方法
CN106679676B (zh) 一种单视场多功能光学敏感器及实现方法
CN106873004B (zh) 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法
CN110246194A (zh) 一种像机与惯性测量单元旋转关系快速标定方法
CN107453811B (zh) 一种基于可见光视觉通信的无人机协同slam的方法
CN113627473A (zh) 基于多模态传感器的水面无人艇环境信息融合感知方法
CN110537197A (zh) 图像处理装置、生长调查图像创建系统以及程序
CN106441310B (zh) 一种基于cmos的太阳方位角计算方法
CN108225335A (zh) 一种用于多目偏振视觉的航向角求解方法
CN113819904B (zh) 一种基于天顶点矢量的偏振/vio三维姿态确定方法
CN113379848A (zh) 基于双目ptz相机的目标定位方法
CN111307140A (zh) 一种用于多云天气条件下的大气偏振光定向方法
Ye et al. Extrinsic calibration of a monocular camera and a single line scanning Lidar
CN117310627A (zh) 一种应用于车路协同路侧感知系统的联合标定方法
CN107705272A (zh) 一种空间影像的高精度几何校正方法
CN115100290B (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Guo Lei

Inventor after: Yang Jian

Inventor after: Wang Gang

Inventor after: Niu Ben

Inventor before: Guo Lei

Inventor before: Wang Gang

Inventor before: Yang Jian

Inventor before: Niu Ben

GR01 Patent grant
GR01 Patent grant