CN110426047B - 基于成像理论的大气偏振模式的建立及其畸变程度分析方法 - Google Patents
基于成像理论的大气偏振模式的建立及其畸变程度分析方法 Download PDFInfo
- Publication number
- CN110426047B CN110426047B CN201910818304.1A CN201910818304A CN110426047B CN 110426047 B CN110426047 B CN 110426047B CN 201910818304 A CN201910818304 A CN 201910818304A CN 110426047 B CN110426047 B CN 110426047B
- Authority
- CN
- China
- Prior art keywords
- formula
- imaging
- mode
- angle
- obtaining
- 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
Links
Images
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/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Abstract
本发明公开了一种基于成像理论的大气偏振模式建模方法,其步骤包括:1根据瑞利模型,构建大气偏振模式;2分析大气偏振模式在相似成像理论和实际采集的非相似成像理论情况下的差异;3分析获取大气偏振模式的偏振成像系统的光路传输特性;4根据偏振成像系统的光路传输特性,构建一种基于成像理论的大气偏振模式建模方法;5根据该模型计算出大气偏振模式经过偏振成像系统后的畸变程度。本发明能够实现不同成像系统、不同太阳位置以及不同观测位置等条件下,实际天空中的大气偏振模式的合理描述,对偏振数据解算精度有所帮助,也为建立更准确的可解析的大气偏振模式表征模型提供了新思路。
Description
技术领域
本发明涉及一种基于成像理论的大气偏振模式的建模及其畸变程度分析方法,属于光学遥感研究领域。
背景技术
太阳光进入大气层后,在传输的过程中受到大气的散射、吸收等作用发生了偏振,进一步形成了具有稳定分布的天空偏振模式。这种天空偏振模式蕴含的偏振信息可以作为不同动物导航信息来源,例如沙蚁,蝗虫,蜜蜂等生物能够利用其独特的复眼结构感知天空偏振模式,为其提供准确的罗盘信息,从而实现导航能力。仿生偏振光导航是一种以生物高度敏感的偏振视觉感知为基础的自主导航方法,这种方法通过对大气偏振模式的检测和演算来实现导航信息的获取。
为了获取全天域的大气偏振模式,采集系统前加装了鱼眼镜头。然而目前的大气偏振模式表征模型并没有充分地描述经过成像系统畸变后的这一过程。尤其是忽略了鱼眼镜头对偏振特性光路的改变,这在很大程度上造成了大气偏振模型的不精确,从而导致理论模型与实际大气偏振模式的采集结果相差甚大。
发明内容
本发明为克服现有技术存在的不足之处,提供一种基于成像理论的大气偏振模式的建立及其畸变程度分析方法,以期能在大气偏振模式中更准确的提取可用于偏振光导航的信息,从而实现更精确的自主导航。
本发明为解决技术问题采取如下技术方案:
本发明一种基于成像理论的大气偏振模式建立方法的特点是按如下步骤进行:
步骤1:根据瑞利模型,构建大气偏振模式;
步骤1.1:以观测者位置为原点O,以向着天顶点的方向为Z轴,以正东方向为X轴,以正北方向为Y轴,建立球形空间坐标系OXYZ;
令天空球面上任意一点为其中,r是天空球面的半径,θT是点T的天顶角,是点T的方位角;令太阳在天空球面上的空间位置为其中,θs为太阳的天顶角,为太阳的方位角,令太阳的高度角hs=90°-θs;则利用式(1)得到天空球面中任意一点处的偏振角αT:
步骤1.2:利用式(2)得到散射光的偏振度P:
P=Pmaxsin2γT/(1+cos2γT) (2)
式(2)中,Pmax表示最大偏振度;γT为散射光在点T处的散射角度,并有:
步骤2:将天空球面近似为“物在无限远处”,从而利用式(4)得到物体在无限远时成像高度的理想像高度y0:
y0=f tanω (4)
式(4)中,f代表鱼眼镜头焦距,ω代表镜头的半视场角;
步骤3:分析各个大气偏振模式下的光路传输特性;
步骤3.1:利用式(5)得到等距成像模式下成像高度y0与镜头的半视场角ω的关系式:
y0=fω (5)
步骤3.2:利用式(6)得到等立体角成像模式下成像高度y0与镜头的半视场角ω的关系式:
步骤3.3:利用式(7)得到体视成像模式下成像高度y0与镜头的半视场角ω的关系式:
步骤4:各个成像模式下构建大气偏振模式:
步骤4.1:在等距成像模式下构建大气偏振模式:
在所述球形空间坐标系OXYZ中,令T′代表鱼眼镜头等距成像模式下,点T投影到二维平面OXY上的点,令点T的高度角hT=90°-θT,其中θT等于点T所在镜头的半视场角ω;
步骤4.2:由式(1)和式(8)得到式(9):
由式(3)和式(8)得到式(10):
步骤4.4:由式(1)和式(11)得到式(12):
由式(3)和式(11)得到式(13):
步骤4.6:由式(1)和式(14得到式(15):
由式(3)和式(14)得到式(16):
本发明所述的基于成像理论的大气偏振模式建立方法的畸变程度分析方法的特点是按如下步骤进行:
步骤5.1:利用式(17)和式(18)得到等距成像模式下的径向放大率βr和切向放大率βt:
βr=f (17)
步骤5.2:利用式(19)和式(20)得到等立体角成像模式下的径向放大率βr和切向放大率βt:
步骤5.3:利用式(21)和式(22)得到体视成像模式下的径向放大率βr和切向放大率βt:
与已有技术相比,本发明的有益效果体现在:
1、目前学术界缺少有效的考虑由成像系统引起的大气偏振模式畸变的技术方法,本发明从成像角度上分析了理论模型与实际检测的差异,实现了大气偏振模式与成像系统结合,充分地考虑了由成像系统引起的畸变,能够更合理地描述实际天空中的大气偏振模式分布,丰富了大气偏振模式建模的方法研究,对建立更准确的可解析的大气偏振模式表征模型提供了重要的理论意义,对偏振数据解算的准确性和导航精度有应用价值。
2、本发明考虑了大气分子散射作用引起的天空光偏振特性,光路进入成像系统后的畸变效应,理论模型与实际采集结果的差异等因素,能适用于不同成像系统采集到的大气偏振模式模拟,具有良好的应用适用性。
3、本发明首次定量分析了三种成像模式下大气偏振模式经过成像系统后的畸变程度。大气偏振模式在三种成像系统下的二维表征具有不同的畸变量,且受到天顶角的调控,可为建立更准确的可解析的大气偏振模式表征模型提供理论依据。
附图说明
图1为本发明的方法流程图;
图2为本发明中高度角为0°方位角为90°时等距成像模式下的偏振化方向角;
图3为本发明中高度角为0°方位角为90°时等距成像模式下的偏振度;
图4为本发明中高度角为0°方位角为90°时等立体角成像模式下的偏振化方向角;
图5为本发明中高度角为0°方位角为90°时等立体角成像模式下的偏振度;
图6为本发明中高度角为0°方位角为90°时体视成像模式下的偏振化方向角;
图7为本发明中高度角为0°方位角为90°时体视成像模式下的偏振度;
图8为本发明中大气偏振模式上一微小面元与像平面上的对应关系。
具体实施方式
本实施例中,如图1所示,一种基于成像理论的大气偏振模式建模方法能够实现不同成像系统、不同太阳位置以及不同观测位置等条件下,实际天空中的大气偏振模式的合理描述,对偏振数据解算精度有所帮助,也为建立更准确的可解析的大气偏振模式表征模型提供了新思路。具体的说是按如下过程进行:
步骤1:根据瑞利模型,构建大气偏振模式;
步骤1.1:以观测者位置为原点O,以向着天顶点的方向为Z轴,以正东方向为X轴,以正北方向为Y轴,建立球形空间坐标系OXYZ;
令天空球面上任意一点为其中,r是天空球面的半径,θT是点T的天顶角,是点T的方位角;令太阳在天空球面上的空间位置为其中,θs为太阳的天顶角,为太阳的方位角,令太阳的高度角hs=90°-θs;则利用式(1)得到天空球面中任意一点处的偏振角αT:
步骤1.2:利用式(2)得到散射光的偏振度P:
P=Pmaxsin2γT/(1+cos2γT) (2)
式(2)中,Pmax表示最大偏振度(理论值为1,晴朗天空下实际值约为0.7);γT为散射光在点T处的散射角度,并有:
步骤2:将天空球面近似为“物在无限远处”,从而利用式(4)得到物体在无限远时成像高度的理想像高度y0:
y0=f tanω (4)
在普通光学系统,成像遵循相似成像理论,即当物体在近距离处,成像高度与物体实际高度遵循式(5):
y0=βy (5)
式(4)中,f代表鱼眼镜头焦距,ω代表镜头的半视场角,在式(5)中,β代表横向放大率,为一定值,y代表实际物体高度,f代表鱼眼镜头焦距。在实际通过拍摄天空成像来获取大气偏振模式时,天空近似为“物在无限远处”,因此,大气偏振模式成像应该遵循式(4);
分析式(4)可知,当ω→90°时,tanω→∞,y0→∞。因此,大气偏振模式在成像时,如果仍然按照相似成像理论成像,那么当视场达到90°时,会出现成像平面无穷大的局面。所以,需要引用“非相似成像”理论。通过牺牲成像的相似性,进而达到成像平面尺寸的理想化。实际采集大气偏振模式时,通过鱼眼镜头的畸变实现非相似成像,从而达到成像平面尺寸的理想化。因此在构建大气偏振模式表征模型时需要分析鱼眼镜头的畸变原理,进一步在成像理论的基础上建立大气偏振模式二维表征模型。
步骤3:分析各个成像模式下大气偏振模式下的光路传输特性;
在实际偏振成像系统的鱼眼镜头成像中,具体成像模式可分为等距成像模式,等立体角成像模式,体视成像模式等。
步骤3.1:利用式(6)得到等距成像模式下成像高度y0与镜头的半视场角ω的关系式:
y0=fω (6)
步骤3.2:利用式(7)得到等立体角成像模式下成像高度y0与镜头的半视场角ω的关系式:
步骤3.3:利用式(8)得到体视成像模式下成像高度y0与镜头的半视场角ω的关系式:
步骤4:各个成像模式下构建大气偏振模式:
步骤4.1:在等距成像模式下构建大气偏振模式:
在球形空间坐标系OXYZ中,令T′代表鱼眼镜头等距成像模式下,点T投影到二维平面OXY上的点,令点T的高度角hT=90°-θT,其中θT等于点T所在镜头的半视场角ω;
步骤4.2:由式(1)和式(10)得到式(11):
由式(3)和式(10)得到式(12):
图2、图3分别为高度角为0°方位角为90°时等距成像模式下的偏振化方向角与高度角为0°方位角为90°时等距成像模式下的偏振度;
步骤4.4:由式(1)和式(13)得到式(14):
由式(3)和式(13)得到式(15):
图4、图5分别为高度角为0°方位角为90°时等立体角成像模式下的偏振化方向角与高度角为0°方位角为90°时等立体角成像模式下的偏振度;
步骤4.6:由式(1)和式(16)得到式(17):
由式(3)和式(16)得到式(18):
图6、图7分别为高度角为0°方位角为90°时体视成像模式下的偏振化方向角与高度角为0°方位角为90°时体视成像模式下的偏振度;
本实施例中,一种基于成像理论的大气偏振模式建立方法的畸变程度分析方法是按如下步骤进行:
图8为大气偏振模式上一微小面元PQVS与像平面上P'Q'V'S'的对应关系,其中XOY为成像平面,光从Z轴射入,经过成像系统后成像在XOY面上。
取天球上一微小面元PQVS,按照等距成像模式,等立体角成像模式,体视成像模式依次成像到XOY面上,对应P'Q'V'S'。径向放大率βr和切向放大率βt是衡量鱼眼镜头成像下产生畸变量的大小,表述形式如式(19);分别计算PQ=Rdθ、P'Q'=dy0、 利用式(19),解算出三种成像模式下大气偏振模式畸变量的大小。
步骤5.1:利用式(20)和式(21)得到等距成像模式下的径向放大率βr和切向放大率βt:
根据等距成像模式下,y0=fθ,进一步解得:
βr=f (20)
步骤5.2:利用式(22)和式(23)得到等立体角成像模式下的径向放大率βr和切向放大率βt:
步骤5.3:利用式(24)和式(25)得到体视成像模式下的径向放大率βr和切向放大率βt:
根据大气偏振模式的径向放大率βr和切向放大率βt的表达式可以看出,大气偏振模式在三种成像系统下具有不同的畸变量,且受到角度θ的调控。θ越大,产生的畸变量会越大,即大气偏振模式在经过三种成像系统在天顶区域的畸变量最小,地平附近的畸变量最大。
Claims (2)
1.一种基于成像理论的大气偏振模式建立方法,其特征是按如下步骤进行:
步骤1:根据瑞利模型,构建大气偏振模式;
步骤1.1:以观测者位置为原点O,以向着天顶点的方向为Z轴,以正东方向为X轴,以正北方向为Y轴,建立球形空间坐标系OXYZ;
令天空球面上任意一点为其中,r是天空球面的半径,θT是点T的天顶角,是点T的方位角;令太阳在天空球面上的空间位置为其中,θs为太阳的天顶角,为太阳的方位角,令太阳的高度角hs=90°-θs;则利用式(1)得到天空球面中任意一点处的偏振角αT:
步骤1.2:利用式(2)得到散射光的偏振度P:
P=Pmaxsin2γT/(1+cos2γT) (2)
式(2)中,Pmax表示最大偏振度;γT为散射光在点T处的散射角度,并有:
步骤2:将天空球面近似为“物在无限远处”,从而利用式(4)得到物体在无限远时成像高度的理想像高度y0:
y0=f tanω (4)
式(4)中,f代表鱼眼镜头焦距,ω代表镜头的半视场角;
步骤3:分析各个大气偏振模式下的光路传输特性;
步骤3.1:利用式(5)得到等距成像模式下成像高度y0与镜头的半视场角ω的关系式:
y0=fω (5)
步骤3.2:利用式(6)得到等立体角成像模式下成像高度y0与镜头的半视场角ω的关系式:
步骤3.3:利用式(7)得到体视成像模式下成像高度y0与镜头的半视场角ω的关系式:
步骤4:各个成像模式下构建大气偏振模式:
步骤4.1:在等距成像模式下构建大气偏振模式:
在所述球形空间坐标系OXYZ中,令T′代表鱼眼镜头等距成像模式下,点T投影到二维平面OXY上的点,令点T的高度角hT=90°-θT,其中θT等于点T所在镜头的半视场角ω;
步骤4.2:由式(1)和式(8)得到式(9):
由式(3)和式(8)得到式(10):
步骤4.4:由式(1)和式(11)得到式(12):
由式(3)和式(11)得到式(13):
步骤4.6:由式(1)和式(14)得到式(15):
由式(3)和式(14)得到式(16):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910818304.1A CN110426047B (zh) | 2019-08-30 | 2019-08-30 | 基于成像理论的大气偏振模式的建立及其畸变程度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910818304.1A CN110426047B (zh) | 2019-08-30 | 2019-08-30 | 基于成像理论的大气偏振模式的建立及其畸变程度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110426047A CN110426047A (zh) | 2019-11-08 |
CN110426047B true CN110426047B (zh) | 2021-05-04 |
Family
ID=68418342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910818304.1A Active CN110426047B (zh) | 2019-08-30 | 2019-08-30 | 基于成像理论的大气偏振模式的建立及其畸变程度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110426047B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539458A (zh) * | 2009-04-30 | 2009-09-23 | 合肥工业大学 | 全自动大气偏振模式图像获取系统及系统控制方式 |
CN103544679A (zh) * | 2012-07-09 | 2014-01-29 | 中国电力科学研究院 | 一种全天空成像仪地基云图畸变校正方法 |
CN104713555A (zh) * | 2015-03-03 | 2015-06-17 | 南昌大学 | 应用全天域中性点辅助定向的车辆自主导航方法 |
CN105209868A (zh) * | 2013-01-16 | 2015-12-30 | 莱顿大学 | 用于启用宽视场偏振测量的装置 |
CN105608706A (zh) * | 2016-01-18 | 2016-05-25 | 中国人民解放军国防科学技术大学 | 一种偏振视觉传感器结构设计与几何标定方法 |
CN108225563A (zh) * | 2017-12-26 | 2018-06-29 | 北京航空航天大学 | 野外环境天空偏振建模对准测量装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105628622A (zh) * | 2014-11-07 | 2016-06-01 | 田佳聪 | 一种基于三相机偏振成像系统 |
-
2019
- 2019-08-30 CN CN201910818304.1A patent/CN110426047B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539458A (zh) * | 2009-04-30 | 2009-09-23 | 合肥工业大学 | 全自动大气偏振模式图像获取系统及系统控制方式 |
CN103544679A (zh) * | 2012-07-09 | 2014-01-29 | 中国电力科学研究院 | 一种全天空成像仪地基云图畸变校正方法 |
CN105209868A (zh) * | 2013-01-16 | 2015-12-30 | 莱顿大学 | 用于启用宽视场偏振测量的装置 |
CN104713555A (zh) * | 2015-03-03 | 2015-06-17 | 南昌大学 | 应用全天域中性点辅助定向的车辆自主导航方法 |
CN105608706A (zh) * | 2016-01-18 | 2016-05-25 | 中国人民解放军国防科学技术大学 | 一种偏振视觉传感器结构设计与几何标定方法 |
CN108225563A (zh) * | 2017-12-26 | 2018-06-29 | 北京航空航天大学 | 野外环境天空偏振建模对准测量装置 |
Non-Patent Citations (6)
Title |
---|
A generic camera model and calibration method for conventional, wide-angle, and;Kannala, J.; Brandt, S.S.;《IEEE T Pattern Anal》;20060831;第28卷(第8期);第1335–1340页 * |
Digital All-Sky Polarization Imaging of Partly Cloudy Skies;Nathan J. Pust,Joseph A. Shaw;《Applied Optics》;20081231;第1-28页 * |
基于偏振成像技术的空间大气偏振模式分布获取;张楠等;《激光与光电子学进展》;20151231;第1-8页 * |
基于复球面映射的大气偏振模式表征与分析;吴良海等;《仪器仪表学报》;20110430;第32卷(第4期);第870-876页 * |
基于大气偏振模式对称性检测的航向角获取方法;田柳等;《信 号 处 理》;20111130;第27卷(第11期);第1701-1705页 * |
眼相机恒星法检校技术研究;原玉磊;《中国博士学位论文全文数据库 信息科技辑》;20130615;第10-12页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110426047A (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Real-time multi-target localization from unmanned aerial vehicles | |
Pang et al. | LiCHy: The CAF’s LiDAR, CCD and hyperspectral integrated airborne observation system | |
Hakala et al. | Acquisition of bidirectional reflectance factor dataset using a micro unmanned aerial vehicle and a consumer camera | |
Li et al. | Quick image-processing method of UAV without control points data in earthquake disaster area | |
CN108387206B (zh) | 一种基于地平线与偏振光的载体三维姿态获取方法 | |
CN113340277B (zh) | 一种基于无人机倾斜摄影的高精度定位方法 | |
CN109900274B (zh) | 一种图像匹配方法及系统 | |
CN104913780A (zh) | 集成gnss和ccd天顶筒的高精度垂线偏差快速测量方法 | |
Karpenko et al. | Visual navigation of the UAVs on the basis of 3D natural landmarks | |
Li et al. | Retrieving 2-D leaf angle distributions for deciduous trees from terrestrial laser scanner data | |
Konovalenko et al. | New approaches to the integration of navigation systems for autonomous unmanned vehicles (UAV) | |
Doneus et al. | Automated archiving of archaeological aerial images | |
CN112529957A (zh) | 确定摄像设备位姿的方法和装置、存储介质、电子设备 | |
Han et al. | A novel orientation method for polarized light compass under tilted conditions | |
CN117036300A (zh) | 基于点云-rgb异源图像多级配准映射的路面裂缝识别方法 | |
Alves Júnior et al. | High accuracy mapping with cartographic assessment for a fixed-wing remotely piloted aircraft system | |
Zhou et al. | Automatic orthorectification and mosaicking of oblique images from a zoom lens aerial camera | |
Pérez-Álvarez et al. | Mining exploration with UAV, low-cost thermal cameras and GIS tools—application to the specific case of the complex sulfides hosted in Carbonates of Udías (Cantabria, Spain) | |
Liu et al. | Attitude calculation method based on full-sky atmospheric polarization mode | |
CN112098926B (zh) | 一种利用无人机平台的智能测角训练样本生成方法 | |
CN110426047B (zh) | 基于成像理论的大气偏振模式的建立及其畸变程度分析方法 | |
CN112950671B (zh) | 一种无人机对运动目标实时高精度参数测量方法 | |
Fan et al. | Bio-inspired multisensor navigation system based on the skylight compass and visual place recognition for unmanned aerial vehicles | |
CN1556374A (zh) | 机载推帚式多维成像装置 | |
Mohamad Azmi et al. | 3D data fusion using unmanned aerial vehicle (UAV) photogrammetry and terrestrial laser scanner (TLS) |
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 |