CN111220150B - Sun vector calculation method based on underwater polarization distribution mode - Google Patents
Sun vector calculation method based on underwater polarization distribution mode Download PDFInfo
- Publication number
- CN111220150B CN111220150B CN201911252040.4A CN201911252040A CN111220150B CN 111220150 B CN111220150 B CN 111220150B CN 201911252040 A CN201911252040 A CN 201911252040A CN 111220150 B CN111220150 B CN 111220150B
- Authority
- CN
- China
- Prior art keywords
- vector
- polarization
- underwater
- angle
- refraction
- 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
- 239000013598 vector Substances 0.000 title claims abstract description 121
- 230000010287 polarization Effects 0.000 title claims abstract description 94
- 238000004364 calculation method Methods 0.000 title claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 26
- 238000000034 method Methods 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 2
- 230000009471 action Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 2
- 239000011664 nicotinic acid Substances 0.000 description 2
- 241000256837 Apidae Species 0.000 description 1
- 241000288673 Chiroptera Species 0.000 description 1
- 241001327682 Oncorhynchus mykiss irideus Species 0.000 description 1
- 241000736128 Solenopsis invicta Species 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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/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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention relates to a sun vector calculation method based on an underwater polarization distribution mode. Firstly, acquiring underwater polarization azimuth angle information in at least two arbitrary observation directions under an underwater celestial coordinate system; then, inverting the atmospheric incident light propagation direction in the underwater observation direction based on the light refraction law; then, determining an atmospheric polarization E-vector corresponding to the polarization E-vector observed underwater by utilizing a Fresnel refraction formula; and finally, obtaining a solar vector by utilizing the vertical relation between the atmospheric polarization E-vector and the solar vector according to the Rayleigh scattering model. The invention establishes a solar vector calculation basic model based on the polarization distribution mode in the underwater Snell window under the horizontal attitude, and can be popularized and applied to underwater application scenes of point source type and image type polarization sensors.
Description
Technical Field
The invention relates to a method for calculating a solar vector based on an underwater polarization distribution mode, which considers the change of refraction on an underwater E-vector, realizes the calculation of the solar vector by utilizing the polarization distribution mode in an underwater Snell window and belongs to the field of underwater polarized light navigation.
Background
At present, the information available for underwater autonomous navigation is very deficient. Inertial navigation is a navigation mode widely applied at present, has the advantages of full autonomy, strong anti-interference performance and the like, but errors of the inertial navigation can be accumulated over time, so that other navigation modes are needed for assisting correction. Because electromagnetic waves are seriously attenuated in water, the traditional radio and satellite navigation cannot be applied underwater; on the one hand, in unknown unfamiliar sea areas, prior geophysical field information is lacked, so that navigation based on the geophysical field, such as geomagnetism and terrain matching, is difficult to realize; on the other hand, no preset underwater transponder is arranged in a strange sea area, so that the underwater sound positioning and other means cannot be applied. The outstanding requirement of current underwater navigation is to realize auxiliary correction of inertial navigation by means of error-free accumulation and independent of prior information.
Biological studies have shown that organisms such as solenopsis invicta, honeybees, bats, rainbow trout, etc. can navigate through sensitive atmospheric polarization distribution patterns. The sunlight generates Rayleigh scattering with particles such as atmospheric molecules in the process of passing through the atmosphere, and a stable and regular polarized light distribution mode is presented on the celestial sphere. Spatial information is contained in the optical field, and the position information of the sun can be calculated through sensing of atmospheric polarization distribution, so that functions of navigation positioning and the like are realized. The method for acquiring the navigation information by sensing the natural environment has the advantages of autonomy, stability, no error accumulation and the like.
The current bionic polarization navigation research mainly focuses on the application scene of the atmosphere, the polarization degree and the polarization E-vector are two parameters of the atmospheric polarization distribution mode, and the key point of the polarization navigation research is how to extract spatial information from the perceived polarization distribution mode, namely the polarization degree or the polarization E-vector distribution mode. For example, chinese patent CN201410652332.8 and chinese patent CN201510303533.1 respectively propose a method for resolving a sun vector based on atmospheric polarization degree distribution characteristics, and chinese patent CN201710027484.2 propose a method for obtaining a sun azimuth by using an atmospheric polarization E-vector. The intrinsic models of the methods are all based on a single Rayleigh scattering model, and the solar vector is solved by utilizing the geometric relationship of the polarization E-vector, the solar perpendicularity and the like. In the above patent, the calculation of the solar vector is a model established only based on the atmospheric scattering effect, and the change of the water body refraction effect on the polarization E-vector direction is not considered, so the underwater polarization distribution does not completely conform to the single Rayleigh scattering model, and therefore, the above patent method is not suitable for the calculation of the solar vector based on the underwater polarization distribution.
Disclosure of Invention
In order to solve the technical problem, the invention provides a solar vector calculation method based on a polarization distribution mode in an underwater Snell window on the basis of the research on the underwater distribution mode of a single Rayleigh scattering model under the refraction action. The underwater polarized light mainly has two sources, namely, a polarized distribution mode formed in a Snell window after the atmospheric polarized light is refracted; the other source is the polarization distribution formed by the sunlight under the scattering action of the water medium after penetrating into the water, and exists inside and outside the Snell window, but the polarized light intensity of the source is weaker, and can be ignored compared with the polarized light after refraction in the Snell window, and only the polarized light of the source exists outside the Snell window. Because the Snell window is transmitted light, the light intensity is stronger, so the polarization distribution mode of the part is more added to the acquisition of navigation information from the perspective of the sensor. The influence of the refraction effect of water on the underwater distribution of atmospheric polarized light mainly has two aspects, namely the refraction changes the propagation direction of non-vertical incident light, and the refraction effect can influence the polarization E-vector angle. Therefore, in the method, the acquired polarization distribution mode in the underwater Snell window can be used for resolving the sun vector so as to obtain navigation information.
The technical scheme adopted by the invention for solving the technical problems is as follows: a sun vector resolving method based on an underwater polarization distribution mode comprises the following implementation steps:
firstly, acquiring any n observation directions under an underwater celestial coordinate system through a horizontally placed polarization sensorUnderwater polarization E-vector information onWherein k is 1, 2.., n; n is more than or equal to 2;
step (2) combining the underwater observation direction of the step (1)Reflecting the incident angle i of the light in the underwater observation direction based on the law of refractionk;
And (3) combining a Fresnel refraction formula, and obtaining the polarization E-vector of the underwater observation obtained in the step (1)And the light incidence angle i calculated in the step (2)kInverting its corresponding atmospheric polarization E-vector
Step (4) according to the atmospheric polarization E-vectorOrthogonal relation to sun vector, using multiple atmospheric polarization E-vectorsAnd resolving to obtain a sun vector s.
Further, the underwater observation direction of the step (1)The specific requirements are as follows:
for any observation directionThe angle from zenith to observation, i.e. zenith angle, is ζk(ii) a Determining zenith angle according to Snell window visual angle thetaUnder the coordinate system of the northeast China,expressed as:
wherein,to observe the directionAnd defining the north direction as the azimuth angle zero position and the north east as the positive direction of the corresponding azimuth angle.
Further, the step (2) is based onThe law of refraction of light is inverted to the incident angle i of light in the underwater observation directionkThe method is concretely realized as follows:
observing direction under waterAngle of refraction r of lightkThe incident angle i of the incident light on water is obtained by the law of refraction of lightkComprises the following steps:
wherein n isair,nwaterRefractive indices of air and water, respectively; angle of refraction r of refracted rays in underwater observation directionkZenith angle zeta with observation directionkEqual, while the surface refraction does not change the azimuth of the ray propagation.
Further, the combination of Fresnel refraction formula in the step (3) and the polarization E-vector of the underwater observation obtained in the step (1)And the light incidence angle i calculated in the step (2)kInverting its corresponding atmospheric polarization E-vector
Setting the included angle between a certain underwater polarization E-vector and the local meridian, namely the underwater polarization azimuth angle is chiwater(ii) a The underwater E-vector is formed by refracting a polarized light in the atmosphere, and the polarization azimuth angle of the polarized E-vector in the atmosphere is Xair(ii) a Then there is the following relationship:
wherein σ⊥,σ||Refractive indices representing a pair of orthogonal components of the E-vector normal and parallel to the plane of incidence, respectively, from FresnelThe equation for the molar refraction yields:
therefore, the method comprises the following steps:
the E-vector is then expressed in the northeast coordinate system as:
eair=vcosχair+usinχair
wherein v, u respectively represent the tangential unit vectors of the longitude direction and the latitude direction of the observation point, namely:
so that any underwater observation directionThe waterborne polarization E-vector corresponding to the upper light ray is as follows:
further, the step (4) utilizes the atmospheric polarization E-vectorPerpendicular to the sun vector, using n atmospheric polarization E-vectorsThe sun vector s is obtained by resolving, and the specific method is as follows:
the sun vector can be obtained by randomly selecting two of n atmospheric polarization E-vectors:
namely, it isWherein,representing the E-vector obtained at n inversionsAny two different amounts of; wherein k is 1,2, 3.., n; to fully utilize the acquired n polarization E-vectors, let:
where m ∈ {1,2, 3.,. n-1} then:
E1×E2=(λ1m+1+λ1m+2+…+λmn)s'
therefore, the method comprises the following steps:
in the above formula, s' is the ambiguous sun vector, and since the z direction of the sun vector is positive in the northeast sky coordinate system, the zenith vector z is introduced as [ 001 ═ c]TAnd eliminating ambiguity through the following formula to obtain the sun vector under the northeast coordinate system as:
s=sign[s'·z]s'
wherein sign [ a · b ] represents the sign of the inner product of vectors a and b.
Has the advantages that:
compared with the prior art, the invention has the following advantages: the existing solar vector calculating method based on the polarization distribution mode is used in an atmospheric scene, and the solar vector calculating method based on the underwater polarization distribution is not available. The underwater polarization distribution is different from an atmospheric polarization distribution mode under the refraction action of water, the method can realize the solar vector calculation based on the polarization distribution mode in the underwater Snell window, and the application scene of the bionic polarized light navigation is expanded from the atmosphere to the underwater.
Drawings
FIG. 1 is a flow chart of a solar vector calculation method based on an underwater polarization distribution mode according to the present invention;
FIG. 2 is an underwater polarization distribution three-dimensional spatial coordinate system according to the present invention;
FIG. 3 is a diagram of the principle of refraction of polarized light on water.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, rather than all embodiments, and all other embodiments obtained by a person skilled in the art based on the embodiments of the present invention belong to the protection scope of the present invention without creative efforts.
According to an embodiment of the present invention, as shown in fig. 1, a solar vector solution method based on an underwater polarization distribution mode of the present invention includes the following steps:
step 1, obtaining any n observation directions under an underwater celestial coordinate system by utilizing a polarization sensorUnderwater polarization E-vector information onWherein k is 1, 2.., n; n.gtoreq.2, e.g.FIG. 2 shows a unit hemisphere as an underwater celestial coordinate system, which represents an upward observation field from an observation point O under water; the circle shown by the dot-dash line in FIG. 2 is the Snell window, and if the view angle of the Snell window is θ, the zenith angle is determinedTherefore, under the coordinate system of the northeast China,expressed as:
wherein,to observe the directionAnd defining the north direction as the azimuth angle zero position and the north east as the positive direction of the corresponding azimuth angle.
Step 2, inverting the incident angle i of the light in the underwater observation direction based on the light refraction lawkUnderwater observation directionAngle of refraction r of lightkAngle of incidence i of incident light on water from the law of refraction of lightkComprises the following steps:
wherein n isair,nwaterRefractive indices of air and water, respectively; angle of refraction r of refracted rays in underwater observation directionkZenith angle zeta with observation directionkEqual, while the water surface refraction does not change the azimuth of the light propagation.
Step 3, combining with a Fresnel refraction formula,polarization E-vector observed from underwaterAnd the light incidence angle i calculated in the step (2)kInverting its corresponding atmospheric polarization E-vector
As shown in FIG. 2, the P point is any observed point in the celestial coordinate system, and a right-hand rectangular coordinate system is established by passing the P point, wherein x isPAxis down the meridian tangent, yPAxial direction is positive east along tangential direction of latitude, zPThe axis is upward along the normal direction of the spherical surface of the point P; the thick line segment of the double arrows passing through the P point represents an E-vector in the underwater observation direction, the vector is a spherical tangent, and the included angle between the vector and the local meridian is the polarization azimuth angle chi of the E-vectorwater(ii) a As shown in FIG. 3, the atmospheric polarization E-vector can be decomposed into two orthogonal components, E, perpendicular and parallel to the plane of incidence⊥And E||(ii) a Atmospheric polarization E-vector EairPolarization E-vector E formed underwater after refractionwater(ii) a Since the refractive action of water is not uniform for the refractive indices of the two orthogonal components, the following relationship holds:
T⊥,T||respectively representing underwater polarization E-vector EwaterOrthogonal components at right angles, parallel to the plane of incidence; sigma⊥,σ||Respectively represents the refractive indexes of orthogonal components of the E-vector which are vertical to and parallel to the incident surface, and the refractive indexes are respectively shown as follows according to a Fresnel refractive formula:
let the polarization azimuth angle of the polarization E-vector in the atmosphere be chiairThen, the following relationship exists:
and because the zenith angle is geometrically equal to the refraction angle, namely ζ ═ r knows:
therefore, the method comprises the following steps:
the E-vector is then expressed in the northeast coordinate system as:
eair=vcosχair+usinχair
wherein v, u respectively represent the tangential unit vector passing through the observation point in the warp direction and the weft direction, namely:
so that the direction is observed underwaterThe waterborne polarization E-vector corresponding to the upper light ray is as follows:
step 4, the specific method for obtaining the sun vector by utilizing the vertical relation between the atmospheric polarization E-vector and the sun vector is as follows:
the sun vector can be obtained by randomly selecting two of n atmospheric polarization E-vectors:
namely, it isWhereinRepresenting the E-vector obtained at n inversionsTwo different E-vectors are arbitrarily selected;
to fully utilize the n polarization E-vectors obtained by the polarization sensor, let:
where m ∈ {1,2,3,.., n-1}, then:
E1×E2=(λ1m+1+λ1m+2+…+λmn)s'
therefore, the method comprises the following steps:
in the above formula, s' is the ambiguous sun vector, and since the z direction of the sun vector is positive in the northeast sky coordinate system, the zenith vector z is introduced as [ 001 ═ c]TAnd disambiguation by the following equationAnd obtaining the sun vector under the northeast coordinate system as follows:
s=sign[s'·z]s'
wherein sign [ a · b ] represents the sign of the inner product of vectors a and b.
Although illustrative embodiments of the present invention have been described above to facilitate the understanding of the present invention by those skilled in the art, it should be understood that the present invention is not limited to the scope of the embodiments, but various changes may be apparent to those skilled in the art, and it is intended that all inventive concepts utilizing the inventive concepts set forth herein be protected without departing from the spirit and scope of the present invention as defined and limited by the appended claims.
Claims (4)
1. A sun vector resolving method based on an underwater polarization distribution mode is characterized by comprising the following steps:
firstly, acquiring any n observation directions under an underwater celestial coordinate system through a horizontally placed polarization sensorUnderwater polarization E-vector ofWherein k is 1, 2.., n; n is more than or equal to 2;
step (2) combining the underwater observation direction of the step (1)Reflecting the incident angle i of the light in the underwater observation direction based on the law of refractionk;
And (3) combining a Fresnel refraction formula, and obtaining the polarization E-vector of the underwater observation obtained in the step (1)And the light incidence angle i calculated in the step (2)kInverting its corresponding atmospheric polarization E-vector
The polarization E-vector of the underwater observation obtained in the step (1) is combined with the Fresnel refraction formula in the step (3)And the light incidence angle i calculated in the step (2)kInverting its corresponding atmospheric polarization E-vector
Setting the included angle between a certain underwater polarization E-vector and the local meridian, namely the underwater polarization azimuth angle is chiwater(ii) a The underwater polarization E-vector is formed by refracting a certain polarized light in the atmosphere, and the polarization azimuth angle of the atmospheric polarization E-vector is Xair(ii) a Then there is the following relationship:
wherein σ⊥,σ||The refractive indices, which represent a pair of orthogonal components of the E-vector, perpendicular and parallel to the plane of incidence, respectively, are given by the Fresnel refractive equation:
therefore, the method comprises the following steps:
the E-vector is then expressed in the northeast coordinate system as:
eair=v cosχair+u sinχair
wherein v, u respectively represent the tangential unit vector passing through the observation point in the warp direction and the weft direction, namely:
so that any underwater observation directionThe corresponding polarization E-vector of the upper ray is:
2. The underwater polarization distribution mode-based solar vector calculation method according to claim 1, wherein:
for any observation directionFrom zenith to observerThe angle of the direction, i.e. zenith angle, being ζk(ii) a Determining zenith angle according to Snell window visual angle thetaUnder the coordinate system of the northeast China,expressed as:
3. The underwater polarization distribution mode-based solar vector calculation method according to claim 1, wherein:
the step (2) is used for inverting the incident angle i of the light in the underwater observation direction based on the light refraction lawkThe method is concretely realized as follows:
observing direction under waterAngle of refraction r of refracted raykThe incident angle i of the incident light on water is obtained by the law of refraction of lightkComprises the following steps:
wherein n isair,nwaterRefractive indices of air and water, respectively;angle of refraction r of refracted rays in underwater observation directionkZenith angle zeta with observation directionkEqual, while the surface refraction does not change the azimuth of the ray propagation.
4. The underwater polarization distribution mode-based solar vector calculation method according to claim 1, wherein:
the step (4) of utilizing the atmospheric polarization E-vectorPerpendicular to the sun vector, using n atmospheric polarization E-vectorsThe sun vector s is obtained by resolving, and the specific method is as follows:
the sun vector can be obtained by randomly selecting two of n atmospheric polarization E-vectors:
namely, it isWherein,representing the E-vector obtained at n inversionsAny two different amounts of; wherein k is 1,2, 3.., n; to fully utilize the acquired n polarization E-vectors, let:
where m ∈ {1,2, 3.,. n-1} then:
E1×E2=(λ1m+1+λ1m+2+…+λmn)s'
therefore, the method comprises the following steps:
in the above formula, s' is the ambiguous sun vector, and since the z direction of the sun vector is positive in the northeast sky coordinate system, the zenith vector z is introduced as [ 001 ═ c]TAnd eliminating ambiguity through the following formula to obtain the sun vector under the northeast coordinate system as:
s=sign[s'·z]s'
wherein sign [ a · b ] represents the sign of the inner product of vectors a and b.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911252040.4A CN111220150B (en) | 2019-12-09 | 2019-12-09 | Sun vector calculation method based on underwater polarization distribution mode |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911252040.4A CN111220150B (en) | 2019-12-09 | 2019-12-09 | Sun vector calculation method based on underwater polarization distribution mode |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111220150A CN111220150A (en) | 2020-06-02 |
CN111220150B true CN111220150B (en) | 2021-09-14 |
Family
ID=70808128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911252040.4A Active CN111220150B (en) | 2019-12-09 | 2019-12-09 | Sun vector calculation method based on underwater polarization distribution mode |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111220150B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112419410B (en) * | 2020-11-20 | 2021-10-19 | 北京航空航天大学 | Horizontal attitude determination method based on underwater Snell window edge identification |
CN112461191B (en) * | 2020-11-20 | 2022-04-05 | 北京航空航天大学 | Sun height calculating method based on underwater refraction and scattering coupling polarization degree |
CN113819904B (en) * | 2021-11-22 | 2022-02-22 | 北京航空航天大学 | polarization/VIO three-dimensional attitude determination method based on zenith vector |
CN114937075B (en) * | 2022-06-06 | 2024-05-14 | 北京航空航天大学 | Underwater polarized light field autonomous orientation method based on three-dimensional solar meridian plane fitting |
CN115014354B (en) * | 2022-06-13 | 2024-04-09 | 北京航空航天大学 | Compass calculation method for self-adaptive fusion of underwater solar light spots and polarized fields |
CN114877898B (en) * | 2022-07-11 | 2022-10-18 | 北京航空航天大学 | Sun dynamic tracking method based on underwater polarization attitude and refraction coupling inversion |
CN114910080B (en) * | 2022-07-15 | 2022-10-21 | 北京航空航天大学 | Three-dimensional navigation attitude determination method based on underwater downlink radiation light intensity and polarized light field |
CN116222550B (en) * | 2023-05-08 | 2023-07-07 | 北京航空航天大学 | Underwater polarized sun calculation method with depth adaptability |
CN116222580B (en) * | 2023-05-09 | 2023-07-25 | 北京航空航天大学 | Underwater polarization orientation method based on cross-medium refraction interference compensation correction |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181145A (en) * | 2015-11-13 | 2015-12-23 | 北京航空航天大学 | Method for simulating underwater polarization field of wave water surface transmitted light |
CN106643704A (en) * | 2017-01-16 | 2017-05-10 | 中国人民解放军国防科学技术大学 | Solar azimuth acquisition method based on atmospheric polarization modes |
CN106767752A (en) * | 2016-11-25 | 2017-05-31 | 北京航空航天大学 | A kind of Combinated navigation method based on polarization information |
CN108225336A (en) * | 2018-01-24 | 2018-06-29 | 北京航空航天大学 | A kind of polarization independent combined navigation method based on confidence level |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7630077B2 (en) * | 2006-10-13 | 2009-12-08 | University Of Massachusetts | System and method for imaging through an irregular water surface |
-
2019
- 2019-12-09 CN CN201911252040.4A patent/CN111220150B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181145A (en) * | 2015-11-13 | 2015-12-23 | 北京航空航天大学 | Method for simulating underwater polarization field of wave water surface transmitted light |
CN106767752A (en) * | 2016-11-25 | 2017-05-31 | 北京航空航天大学 | A kind of Combinated navigation method based on polarization information |
CN106643704A (en) * | 2017-01-16 | 2017-05-10 | 中国人民解放军国防科学技术大学 | Solar azimuth acquisition method based on atmospheric polarization modes |
CN108225336A (en) * | 2018-01-24 | 2018-06-29 | 北京航空航天大学 | A kind of polarization independent combined navigation method based on confidence level |
Non-Patent Citations (2)
Title |
---|
Polarization Patterns of Transmitted Celestial Light under Wavy Water Surfaces;Guanhua Zhou et al.;《REMOTE SENSING》;20170329;第9卷(第4期);正文第1-21页 * |
水下偏振光导航技术;王光辉等;《舰船科学技术》;20110731;第33卷(第7期);第79-82页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111220150A (en) | 2020-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111220150B (en) | Sun vector calculation method based on underwater polarization distribution mode | |
CN106679645B (en) | Real time navigation apparatus based on multi-direction polarised light | |
CN109459027B (en) | Navigation method based on polarization-geomagnetic vector tight combination | |
CN104880191B (en) | A kind of polarization auxiliary navigation method based on solar vector | |
CN107063170B (en) | Course angle estimation method based on atmospheric polarization angle mould formula under complex environment | |
CN101344391A (en) | Lunar vehicle pose self-confirming method based on full-function sun-compass | |
CN102589544A (en) | Three-dimensional attitude acquisition method based on space characteristics of atmospheric polarization mode | |
US11396354B2 (en) | Covert underwater navigation via polarimetry | |
CN112461191B (en) | Sun height calculating method based on underwater refraction and scattering coupling polarization degree | |
CN114877898B (en) | Sun dynamic tracking method based on underwater polarization attitude and refraction coupling inversion | |
CN114910080B (en) | Three-dimensional navigation attitude determination method based on underwater downlink radiation light intensity and polarized light field | |
CN103837150A (en) | Method for performing rapid celestial fix through CCD (charge coupled device) zenith telescope on ground | |
CN115014354B (en) | Compass calculation method for self-adaptive fusion of underwater solar light spots and polarized fields | |
Zhang et al. | A self-contained interactive iteration positioning and orientation coupled navigation method based on skylight polarization | |
CN110887476B (en) | Autonomous course and attitude determination method based on polarization-astronomical included angle information observation | |
Song et al. | HAUD: A high-accuracy underwater dataset for visual-inertial odometry | |
CN114937075B (en) | Underwater polarized light field autonomous orientation method based on three-dimensional solar meridian plane fitting | |
CN106018402B (en) | A kind of the visibility detection system and method for UAV system refractive and reflective panorama stereoscopic camera | |
CN115164872B (en) | Autonomous positioning method based on time sequence polarized light field | |
CN112419410B (en) | Horizontal attitude determination method based on underwater Snell window edge identification | |
Hu et al. | Underwater downwelling radiance fields enable three-dimensional attitude and heading determination | |
CN109974692A (en) | A kind of hidden environment astronomical positioning system and method based on neutrino signal | |
Hu et al. | Underwater autonomous orientation using submarine light intensity gradient | |
CN106767766B (en) | A kind of sky polarization localization method based on single neutral point model | |
Dupeyroux et al. | Bio-inspired celestial compass yields new opportunities for urban localization |
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 |