CN115200543A - Earth sky light polarization basic field observation method based on atmospheric scattering - Google Patents
Earth sky light polarization basic field observation method based on atmospheric scattering Download PDFInfo
- Publication number
- CN115200543A CN115200543A CN202210581758.3A CN202210581758A CN115200543A CN 115200543 A CN115200543 A CN 115200543A CN 202210581758 A CN202210581758 A CN 202210581758A CN 115200543 A CN115200543 A CN 115200543A
- Authority
- CN
- China
- Prior art keywords
- angle
- polarization
- longitude
- time
- light polarization
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 2
- 230000000737 periodic effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 6
- 239000013598 vector Substances 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 241000257303 Hymenoptera Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 241000282414 Homo sapiens Species 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011664 nicotinic acid Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation 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
- G01C1/00—Measuring angles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J4/00—Measuring polarisation of light
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention relates to an earth sky light polarization basic field observation method based on atmospheric scattering, which comprises the following steps: acquiring a time angle and a declination angle of a measuring station at a certain moment based on longitude and latitude coordinates of the measuring station; calculating a solar elevation angle and a solar azimuth angle at each longitude and latitude coordinate at a certain moment of the measuring station based on the hour angle and the declination angle; calculating the real-time polarization degree of the vertical sky at each longitude and latitude coordinate of the world at a certain moment by utilizing a Rayleigh scattering model; and constructing an earth sky light polarization field based on the sun altitude, the sun azimuth and the corresponding real-time polarization degree of each global longitude and latitude coordinate at a certain moment, and finding the global space distribution rule of the earth sky light polarization field. The invention draws the sky polarization field information of continuous time sequence at each longitude and latitude coordinate position of the world, and can find the periodic variation of the world sky light polarization basic field.
Description
Technical Field
The invention relates to a method, a system, equipment and a medium for observing a polarization basic field of earth sky light based on atmospheric scattering, and relates to the field of polarization remote sensing and polarization light navigation.
Background
Sky polarized light information may be used for bio-polarized navigation. Some ants, bees and the like use the polarized light distribution characteristic of the sky to navigate, can quickly position the direction of the nest and select the optimal path to return to the nest. The phenomenon can be found by human beings to invent a bionic polarization navigation instrument for polarization autonomous navigation of robots and the like.
At present, the common use is the local all-sky polarized information of a certain place, and the distribution characteristic of the all-sky polarized light of one place can be obtained. This is a great gap for a large range of even global polarization field information.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide an earth sky light polarization basic field observation method, system, device and medium based on atmospheric scattering, which can construct an earth sky light polarization basic field and obtain a global sky polarization basic field.
In order to achieve the purpose of the invention, the technical scheme provided by the invention is as follows:
in a first aspect, the invention provides an earth sky light polarization basic field observation method based on atmospheric scattering, which includes:
obtaining a polarization degree value of each global longitude and latitude coordinate point in the vertical direction;
and drawing a global sky light polarization field obtained on a global map based on the sun altitude, the sun azimuth angle and the corresponding polarization degree value at each longitude and latitude coordinate at a certain moment, and constructing the global sky light polarization field at the certain moment so as to obtain a spatial distribution rule and a temporal change rule of the global sky light polarization field.
Further, constructing a polarization field of the earth sky light at a certain moment, comprising:
acquiring a time angle and a declination angle of a measuring station at a certain moment based on longitude and latitude coordinates of the measuring station;
calculating a solar altitude angle and a solar azimuth angle at each longitude and latitude coordinate of the measuring station at a certain moment based on the hour angle and the declination angle;
calculating real-time polarization values of vertical sky at each longitude and latitude coordinate of the world at a certain moment by using a Rayleigh scattering model;
and constructing a terrestrial sky light polarization field based on the solar altitude and the solar azimuth at each longitude and latitude coordinate of the whole world at a certain moment and the corresponding real-time polarization value.
Further, obtaining the hour angle and the declination angle of the measurement station at a certain time comprises:
correcting errors of the real solar time T;
based on the corrected true solar time T, the time angle ω and the declination angle δ are calculated.
Further, the formula for correcting the error of the true solar time T is:
T＝t±4min*(120-L LOC )+E
wherein t is a certain time, L LOC The longitude of the measured location is and the corrected time difference is E.
Further, the calculation formula of the hour angle ω and the declination angle δ is:
ω＝(T-12)*15
wherein n is the age day.
Further, the solar altitude angle and the solar azimuth angle at each longitude and latitude coordinate at a certain moment of the measuring station are calculated by the following formulas:
in the formula, h s Is the altitude angle, a s Is the azimuth angle.
In a second aspect, the present invention further provides an earth sky light polarization elementary field observation system based on atmospheric scattering, which includes:
the sun position information calculation unit is configured to obtain a time angle and a declination angle of the measuring station at a certain moment based on the longitude and latitude coordinates of the measuring station, and calculate a sun altitude angle and a sun azimuth angle of each longitude and latitude coordinate of the measuring station at the certain moment based on the time angle and the declination angle;
the polarization degree calculation unit is configured to calculate the real-time polarization degree of the vertical sky at each longitude and latitude coordinate of the world at a certain moment by utilizing a Rayleigh scattering model;
the polarized light field graph drawing unit is configured to construct an earth sky light polarization field based on the sun altitude angle, the sun azimuth angle and the corresponding real-time polarization degree at each longitude and latitude coordinate of the globe at a certain moment, calculate the earth sky light polarization field at each moment and obtain a spatial distribution rule and a time change rule of the earth sky light polarization field.
In a third aspect, the present invention further provides an electronic device, which includes computer program instructions, where the program instructions, when executed by a processor, are configured to implement the method for observing an earth sky light polarization basic field based on atmospheric scattering.
In a fourth aspect, the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium has stored thereon computer program instructions, wherein the program instructions, when executed by a processor, are configured to implement the atmospheric scattering-based earth sky light polarization basic field observation method.
Due to the adoption of the technical scheme, the invention has the following characteristics:
1. the invention utilizes the Rayleigh scattering model to simulate the information of the all-sky polarized field, and has obvious polarized light distribution effect on sunny days.
2. The method calculates the real-time solar altitude and azimuth angle of each longitude and latitude coordinate of the world, further calculates the real-time polarization degree information of the vertical sky of each longitude and latitude coordinate of the world by utilizing a Rayleigh scattering model, draws the sky polarization field information of a certain moment of each longitude and latitude coordinate of the world to obtain the earth sky light polarization field, can obtain the spatial distribution rule of the earth sky light polarization field, further can calculate the earth sky light polarization field of each moment, and can find the time period change of the global sky light polarization basic field, wherein the sky global light polarization basic field has the daily change with a short period of one day and the annual change with a long period of one year.
In conclusion, the method is suitable for calculating the basic field of the sky polarized light field.
Drawings
Various additional advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Like reference numerals refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a flow chart of a method according to an embodiment of the present invention.
FIG. 2 illustrates sky P-point e vector polarization information in accordance with an embodiment of the present invention.
Fig. 3 is a diagram of the effect of global solar altitude distribution according to an embodiment of the present invention (spring equinox, summer solstice, autumn equinox, winter solstice, greenwich mean time 12.
Fig. 4 is a global polarization vertical component effect graph of an embodiment of the present invention (spring equinox, summer equinox, autumn equinox, winter equinox time 06, 00, 09, 12, 00, 15.
FIG. 5 is a diagram illustrating the effect of the degree of polarization of the global polarization field according to the present invention.
FIG. 6 is a graph showing the effect of the polarization degree of the global polarization field according to an embodiment of the present invention.
Detailed Description
It is to be understood that the terminology used herein is for the purpose of describing particular example embodiments only, and is not intended to be limiting. As used herein, the singular forms "a", "an" and "the" may be intended to include the plural forms as well, unless the context clearly indicates otherwise. The terms "comprises," "comprising," "including," and "having" are inclusive and therefore specify the presence of stated features, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, steps, operations, elements, components, and/or groups thereof. The method steps, processes, and operations described herein are not to be construed as necessarily requiring their performance in the particular order discussed or illustrated, unless specifically identified as an order of performance. It should also be understood that additional or alternative steps may be used.
For convenience of description, spatially relative terms, such as "inner", "outer", "lower", "upper", and the like, may be used herein to describe one element or feature's relationship to another element or feature as illustrated in the figures. This spatially relative term is intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures.
The invention provides a method, a system, equipment and a medium for observing an earth sky light polarization basic field based on atmospheric scattering, wherein the method comprises the following steps: and drawing a global map to obtain an earth sky light polarization field based on the sun altitude, the sun azimuth and the corresponding polarization degree value of each latitude and longitude coordinate at a certain moment, and further calculating the earth sky light polarization field at a certain moment to obtain a spatial distribution rule and a time change rule of the earth sky light polarization field. Therefore, the invention can draw the sky polarization field information of each time at each longitude and latitude coordinate of the world, and can find the periodic variation of the world sky light polarization basic field.
Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention can be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
The stokes vector can describe polarized light and consists of four mathematical vectors: i, Q, U and V, wherein the I component represents the overall light intensity, the Q component represents the horizontal and vertical line components, the U component represents the ± 45 ° line components, and the V component describes the left-and right-hand circular polarization.
In the formula I 0° Is the total intensity of light, I 0° 、I 90° 、I +45°° 、I -45°° 、I r And I l Respectively representing the light intensity of linearly polarized light in the directions of 0 degrees, 90 degrees, +45 degrees and 45 degrees, and the light intensity of levorotatory (l) and dextrorotatory (r) circularly polarized light.
The relationship between the incident light intensity and the emergent light intensity can be expressed by the stokes component and the deflection angle:
setting θ to 0 °,60 ° and 120 °, respectively, three calculation formulas can be obtained, and then three stokes components I, Q, U can be obtained.
Further, on the basis of the stokes component, the degree of polarization P of the incident and exit angle can be calculated:
and the polarization azimuth angle of incident light:
for sky scattering dominated by rayleigh scattering, the degree of polarization under the condition that the incident light is natural light is:
where θ is the angle between the scattering direction and the incident direction, i.e. the scattering angle, P max Is the maximum value of the degree of polarization, which is 1.
The first embodiment is as follows: based on the above principle, as shown in fig. 1, the method for observing the earth sky light polarization basic field based on atmospheric scattering provided by this embodiment includes:
s1, determining longitude and latitude of a measuring station, obtaining a time angle and a declination angle by utilizing geographic coordinates (longitude and latitude coordinates) and time of the measuring station based on atmospheric scattering and optical polarization, and further calculating solar altitude angle and solar azimuth angle information of each longitude and latitude coordinate of the measuring station at a certain moment, wherein the method comprises the following steps:
s11, correcting the error of the true sun T.
The calculation of the true solar time T is affected by two errors, namely an error caused by longitude, namely a time difference of 120 degrees from east longitude. Secondly, because of the error caused by the real solar time difference, two errors need to be corrected, and the correction process comprises the following steps:
first, the longitude error is corrected:
T 0 ＝t±4min*(L SI -L LOC ) (9)
where t is a certain time, L LOC For the longitude of the geodetic location, "+" applies to the western hemisphere and "one" to the eastern hemisphere. L is a radical of an alcohol SI Longitude of a standard time zone, L for China SI ＝120°
Next, the error caused by the time difference E is corrected:
E＝0.0172+0.4281 cos Q-7.3515 sin Q-3.3495 cos 2Q-9.3619 sin 2Q (10)
wherein,
and n is the year.
Finally, a real solar time calculation formula can be obtained by combining the two correction results:
T＝t±4min*(120-L LOC )+E (11)
s12, calculating the hour angle omega and the declination angle delta
In the equatorial coordinate system, the position of the sun is determined by both the hour angle ω and the declination angle δ.
The time angle omega is the angular distance from the celestial sphere meridian of the observation station to the time circle of the sun along the equator, and is calculated as shown in formula (7), wherein T is the true solar time;
ω＝(T-12)*15 (7)
the declination angle is the angle between the equatorial plane of the δ earth and the line connecting the sun and the earth's center, and is calculated as shown in equation (8), where n is the age day, i.e., the number of days the test date is 1 month and 1 day apart.
S12, calculating the position of the sun on the celestial sphere based on the time angle omega and the declination angle delta
In the horizon coordinate system, the elevation angle h is usually used s And azimuth angle a s To determine the position of the sun on the celestial sphere.
Height angle h s : the calculation formula of the included angle between the direct solar light and the ground plane where the observation station is located is as follows:
wherein,
the latitude of the observer location.
Azimuth angle a s : the included angle between the projection line of the direct solar ray on the ground plane and the line in the south direction on the ground plane is calculated by the following formula:
according to the calculation formula, the solar altitude and the solar azimuth at any time of a specific place can be obtained.
S2, calculating the polarization degree of the light beam in the observation direction in the sky, and further drawing a global polarization degree graph, wherein the method comprises the following steps:
s21, calculating the polarization degree P and the polarization azimuth phi of the light beam in the observation direction in the sky
As shown in fig. 2, in the horizon coordinate system, S denotes the sun, OP denotes the direction observed in the sky, and Z denotes the zenith. According to Rayleigh scattering theory, the e-vector vibration direction of the light beam to be measured is perpendicular to a plane formed by the sun S, the ground observation point O and the observation direction P point in the sky.
The polarization azimuth angle psi is defined as the included angle between the e-vector vibration direction of a point P in the observation direction in the sky and the meridian (arc ZP) passing the point P, and the e-vector vibration direction of the point P is perpendicular to PS, psi =90 DEG-angle ZPS.
In the spherical triangle ZPS, the cosine theorem of the spherical triangle indicates that:
obtaining:
cosθ＝sinh p ·sinh s +cosh p ·cosh s ·cos(A s -A p ) (15)
where θ is the angle between the incident light direction OS and the observation direction OP, h s Is the altitude of the sun, A s Is the azimuth angle of the sun, h p Is the altitude, A, of an observation point P in the spherical sky p Is the azimuth of the observation point P. The transformation of the coordinates into a true south direction represents an azimuth angle of 0 deg., and a true west direction represents an azimuth angle of 90 deg.. The polarization degree P of the observation direction light beam in the sky can be calculated by the formula.
And S22, forming a polarization field of the earth sky light based on the solar altitude and the solar azimuth of each global position at a certain moment and the polarization degree P of the observation direction light beam in the corresponding sky.
Specifically, a polarization value is arranged in the vertical direction of each coordinate point of the world, and a global polarization degree map is obtained by drawing the global polarization degree map and is used for displaying the earth sky light polarization field. The earth sky light polarization field reveals the distribution rule of the vertical sky polarization information of each global point at a certain moment, and can provide a basis for global polarization navigation and also provide a basis for global polarization remote sensing information correction.
The following describes in detail the application of the earth sky light polarization basic field observation method based on atmospheric scattering according to the present invention by using specific embodiments.
At a certain time of day, a certain position, there is a relatively stable polarization distribution in the sky.
First, the solar altitude at each global position at a certain time is calculated and simulated, as shown in fig. 3.
Secondly, because the polarization degree in the zenith direction has stability and gradual change, global atmospheric polarization degree distribution is simulated for the conditions observed vertically upward from the ground, and the result is shown in fig. 4. The polarization pattern diagram in the sky is closely related to the solar altitude (i.e. the moment), the polarization degrees of different solar altitudes are different, the intensities and directions of polarized light are different, when the solar altitude is low, the polarization intensity in the sky is relatively high, that is, the polarization intensity in the sky is greater than that in the middle of the day in the evening or in the morning, and in different seasons: the distribution of spring minutes, summer solstice, autumn minutes and winter solstice is changed regularly.
FIG. 4 shows four distinct demarcation points in the annual variation of the global polarization field polarization degree. The distribution of the ingredients has regular change in spring equinox, summer solstice, autumn equinox and winter solstice. The central point of the global polarization field reaches the equator in spring time division in March, reaches the return-to-north return line in summer solstice in June, then moves to the south, reaches the equator in autumn equinox, and reaches the return-to-south return line in winter solstice in December. In one year, the global polarization field center moves between the return lines of south and north, crosses the equator twice in spring and autumn, and reaches the farthest positions of south and north in winter and summer, respectively.
Fig. 5 shows the daily change of the polarization degree of the global polarization field, which is taken a month and a day as an example. n represents the day of the year, e.g., 1 month at 1 day, n =1,2 months at 1 day, n =32.UT (Universal Time) refers to coordinated Universal Time, also called Universal Time, universal standard Time and international coordinated Time, which is the standard Time of Greenwich mean. And selecting six time points of 1 month and 1 day on the day, and calculating the theoretical distribution of the polarization degree of the sky in the vertical direction of the top of the sky by using a Rayleigh scattering model. The point at which the degree of polarization is smallest (polarization neutral point) is theoretically the point at which the sun is directly irradiated, and the polarization neutral point changes with the movement of the position of the sun in the diurnal hemisphere as the day and night alternate. From UT =0 to UT =20, the sun moves from east to west, so that the local polarization is the smallest at direct sun rays, while at the morning and evening line the sun's altitude is the smallest and the sky polarization is the largest. Such a stable polarization degree distribution generates a daily change in a period of one day as the earth rotates.
FIG. 6 shows the annual variation of the polarization degree of the global polarization field. Taking the first day of each month as an example, the world standard time UT =12 hours of the theoretical distribution of the degree of polarization of the sky globally. It can be seen that the point with the minimum vertical sky polarization degree is theoretically the point of direct solar radiation, moves between the return lines of south and north, and moves by one year as a period. At month 1, the global polarization field center is in the southern hemisphere. As the month progresses, the central point of minimum polarization gradually moves to the north. In spring time between 3 months and 1 day and 4 months and 1 day, the global polarization field center crosses the equator. During summer, the central point of the global polarization field reaches the return-to-north line in the north-north direction, and then moves to the south. In the autumn time between 9 and 1/10, the global polarization field center again crosses the equator. In winter solstice, the central point of the global polarization field reaches the south-return line from the south to the north, and then moves to the north.
Example two: correspondingly, the embodiment provides an earth sky light polarization basic field observation system based on atmospheric scattering. The system provided by this embodiment may implement the method for observing the earth sky light polarization basic field based on atmospheric scattering according to the first embodiment, and the system may be implemented by software, hardware, or a combination of software and hardware. For convenience of description, the present embodiment is described by dividing functions into various units and describing the units separately. Of course, the functions of the units may be implemented in the same software and/or hardware or in one or more pieces. For example, the system may comprise integrated or separate functional modules or units to perform the corresponding steps in the method of an embodiment. Since the system of the present embodiment is basically similar to the method embodiment, the description process of the present embodiment is relatively simple, and reference may be made to part of the description of the first embodiment to related points.
Specifically, the present embodiment further provides an earth sky light polarization basic field observation system based on atmospheric scattering, which includes:
the solar position information calculation unit is configured to obtain an hour angle and a declination angle of a measuring station at a certain moment based on longitude and latitude coordinates of the measuring station, and calculate a solar elevation angle and a solar azimuth angle at each longitude and latitude coordinate of the measuring station at the certain moment based on the hour angle and the declination angle;
the polarization degree calculation unit is configured to calculate the real-time polarization degree of the vertical sky at each longitude and latitude coordinate of the world at a certain moment by utilizing a Rayleigh scattering model;
the polarized light field graph drawing unit is configured to construct an earth sky light polarization field based on the sun altitude angle, the sun azimuth angle and the corresponding real-time polarization degree of each global longitude and latitude coordinate at a certain moment, calculate the earth sky light polarization field at each moment, and obtain a spatial distribution rule and a time change rule of the earth sky light polarization field. .
Example three: the present embodiment provides an electronic device corresponding to the method for observing a basic field of polarization of earth sky light based on atmospheric scattering provided in the first embodiment, where the electronic device may be an electronic device for a client, such as a mobile phone, a notebook computer, a tablet computer, a desktop computer, etc., to execute the method in the first embodiment.
The electronic equipment comprises a processor, a memory, a communication interface and a bus, wherein the processor, the memory and the communication interface are connected through the bus to complete mutual communication. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, an Extended ISA (Extended Industry Standard Component) bus, or the like. The memory stores a computer program that can be executed on the processor, and the processor executes the computer program to execute the earth sky light polarization basic field observation method based on atmospheric scattering provided by the embodiment.
In some implementations, the logic instructions in the memory may be implemented in software functional units and stored in a computer readable storage medium when sold or used as a stand-alone product. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), an optical disk, and various other media capable of storing program codes.
In other implementations, the processor may be any type of general-purpose processor such as a Central Processing Unit (CPU), a Digital Signal Processor (DSP), and the like, and is not limited herein.
Example four: the method for observing the earth sky light polarization basic field based on atmospheric scattering according to this embodiment may be embodied as a computer program product, which may include a computer readable storage medium having computer readable program instructions embodied thereon for executing the method for observing the earth sky light polarization basic field based on atmospheric scattering according to this embodiment.
The computer readable storage medium may be a tangible device that retains and stores instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but not limited to, an electronic memory device, a magnetic memory device, an optical memory device, an electromagnetic memory device, a semiconductor memory device, or any combination of the foregoing.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description herein, references to the description of "one embodiment," "some implementations," or the like, mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of an embodiment of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (9)
1. An earth sky light polarization basic field observation method based on atmospheric scattering is characterized by comprising the following steps:
obtaining a polarization degree value of each global longitude and latitude coordinate point in the vertical direction;
and drawing a global sky light polarization field obtained on a global map based on the sun altitude, the sun azimuth angle and the corresponding polarization degree value at each longitude and latitude coordinate at a certain moment, and constructing the global sky light polarization field at the certain moment so as to obtain a spatial distribution rule and a temporal change rule of the global sky light polarization field.
2. The earth sky light polarization basic field observation method based on atmospheric scattering according to claim 1, wherein constructing an earth sky light polarization field at a certain moment comprises:
acquiring a time angle and a declination angle of a measuring station at a certain moment based on longitude and latitude coordinates of the measuring station;
calculating a solar altitude angle and a solar azimuth angle at each longitude and latitude coordinate of the measuring station at a certain moment based on the hour angle and the declination angle;
calculating real-time polarization values of vertical sky at each longitude and latitude coordinate of the world at a certain moment by using a Rayleigh scattering model;
and constructing a terrestrial sky light polarization field based on the solar altitude and the solar azimuth at each longitude and latitude coordinate of the whole world at a certain moment and the corresponding real-time polarization value.
3. The method of observing an earth sky light polarization basic field based on atmospheric scattering according to claim 2, wherein obtaining a time angle and a declination angle of a measurement station at a certain time comprises:
correcting errors of the real solar time T;
based on the corrected true solar time T, the time angle ω and the declination angle δ are calculated.
4. The method of observing an earth sky light polarization basic field based on atmospheric scattering according to claim 3, wherein the formula for correcting the error of T in true sun time is as follows:
T＝t±4min*(120-L LOC )+E
wherein t is a certain time, L LOC The longitude of the measured place and the corrected time difference are E.
5. The method of observing an earth sky light polarization basic field based on atmospheric scattering according to claim 3 or 4, wherein the formula for calculating the time angle ω and the declination angle δ is:
ω＝(T-12)*15
wherein n is the age day.
6. The atmospheric scattering-based earth sky light polarization basic field observation method of claim 3, wherein the sun elevation angle and the sun azimuth angle at each latitude and longitude coordinate at a certain time of the measurement site are calculated by the following formulas:
in the formula, h s Is the altitude angle, a s Is the azimuth angle.
7. An atmospheric scattering-based earth sky light polarization elementary field observation system, comprising:
the solar position information calculation unit is configured to obtain an hour angle and a declination angle of a measuring station at a certain moment based on longitude and latitude coordinates of the measuring station, and calculate a solar elevation angle and a solar azimuth angle at each longitude and latitude coordinate of the measuring station at the certain moment based on the hour angle and the declination angle;
the polarization degree calculation unit is configured to calculate a real-time polarization degree value of a vertical sky at each longitude and latitude coordinate of the globe at a certain moment by utilizing a Rayleigh scattering model;
the polarized light field graph drawing unit is configured to construct an earth sky light polarization field based on a sun altitude angle, a sun azimuth angle and a corresponding real-time polarization value of each global longitude and latitude coordinate at a certain moment, calculate the earth sky light polarization field at each moment, and obtain a spatial distribution rule and a time change rule of the earth sky light polarization field.
8. An electronic device comprising computer program instructions, wherein the program instructions when executed by a processor are adapted to implement the atmospheric scattering based earth sky light polarization basic field observation method of any one of claims 1 to 6.
9. A computer-readable storage medium, having stored thereon computer program instructions, wherein the program instructions, when executed by a processor, are configured to implement the atmospheric scattering-based earth sky light polarization elementary field observation method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210581758.3A CN115200543B (en) | 2022-05-26 | 2022-05-26 | Earth sky light polarization basic field observation method based on atmospheric scattering |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210581758.3A CN115200543B (en) | 2022-05-26 | 2022-05-26 | Earth sky light polarization basic field observation method based on atmospheric scattering |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115200543A true CN115200543A (en) | 2022-10-18 |
| CN115200543B CN115200543B (en) | 2023-12-19 |
Family
ID=83576243
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210581758.3A Active CN115200543B (en) | 2022-05-26 | 2022-05-26 | Earth sky light polarization basic field observation method based on atmospheric scattering |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115200543B (en) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101281060A (en) * | 2008-04-30 | 2008-10-08 | 大连理工大学 | Sky polarized light radiation spectrum measuring systems |
| CN101532881A (en) * | 2009-04-03 | 2009-09-16 | 合肥工业大学 | Single factor atmospheric polarization modeling method based on Rayleigh scattering |
| CN101739504A (en) * | 2009-12-16 | 2010-06-16 | 合肥工业大学 | Two-dimensional representation method of atmospheric polarization theoretical model |
| CN101769862A (en) * | 2010-01-15 | 2010-07-07 | 合肥工业大学 | Distribution type atmospheric polarization mode detecting system and system control method |
| CN101887475A (en) * | 2010-06-25 | 2010-11-17 | 合肥工业大学 | Multi-factor atmospheric polarized modeling method based on earth surface Fresnel reflection |
| CN102538783A (en) * | 2012-02-14 | 2012-07-04 | 北京大学 | Bionic navigation method and navigation positioning system based on remote sensing sky polarization mode patterns |
| CN102564404A (en) * | 2012-02-07 | 2012-07-11 | 北京大学 | Polarized remote sensing earth-atmosphere information separation method based on atmosphere neutral point |
| HUP1100482A2 (en) * | 2011-09-05 | 2013-04-29 | Eotvos Lorand Tudomanyegyetem | Method for cloud base height measuring and device for polarization measuring |
| CN103234637A (en) * | 2013-04-04 | 2013-08-07 | 大连理工大学 | Atmospheric radiation spectrum polarization information measuring system |
| CN106767766A (en) * | 2016-11-22 | 2017-05-31 | 北京航空航天大学 | A kind of sky based on single neutrality point model polarizes localization method |
| CN108267228A (en) * | 2016-12-30 | 2018-07-10 | 南京理工大学 | A kind of atmospheric scattering light polarization analysis method based on CE318 heliographs |
| CN108759819A (en) * | 2018-06-08 | 2018-11-06 | 北京航空航天大学 | A kind of polarization navigation real-time location method based on omnimax polarization degree information |
| CN109459015A (en) * | 2018-11-09 | 2019-03-12 | 北京航空航天大学 | A kind of global autonomic positioning method of polarization navigation based on the observation of maximum degree of polarization |
| CN113280785A (en) * | 2021-07-22 | 2021-08-20 | 中国人民解放军国防科技大学 | Sky polarized light orientation method based on polarization mode consistency |
| CN113532419A (en) * | 2021-06-23 | 2021-10-22 | 合肥工业大学 | Sky polarization mode information acquisition method and device, electronic equipment and storage medium |
| CN114076594A (en) * | 2020-08-21 | 2022-02-22 | 大连理工大学 | Distributed polarized light global positioning measurement system |
-
2022
- 2022-05-26 CN CN202210581758.3A patent/CN115200543B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101281060A (en) * | 2008-04-30 | 2008-10-08 | 大连理工大学 | Sky polarized light radiation spectrum measuring systems |
| CN101532881A (en) * | 2009-04-03 | 2009-09-16 | 合肥工业大学 | Single factor atmospheric polarization modeling method based on Rayleigh scattering |
| CN101739504A (en) * | 2009-12-16 | 2010-06-16 | 合肥工业大学 | Two-dimensional representation method of atmospheric polarization theoretical model |
| CN101769862A (en) * | 2010-01-15 | 2010-07-07 | 合肥工业大学 | Distribution type atmospheric polarization mode detecting system and system control method |
| CN101887475A (en) * | 2010-06-25 | 2010-11-17 | 合肥工业大学 | Multi-factor atmospheric polarized modeling method based on earth surface Fresnel reflection |
| HUP1100482A2 (en) * | 2011-09-05 | 2013-04-29 | Eotvos Lorand Tudomanyegyetem | Method for cloud base height measuring and device for polarization measuring |
| CN102564404A (en) * | 2012-02-07 | 2012-07-11 | 北京大学 | Polarized remote sensing earth-atmosphere information separation method based on atmosphere neutral point |
| CN102538783A (en) * | 2012-02-14 | 2012-07-04 | 北京大学 | Bionic navigation method and navigation positioning system based on remote sensing sky polarization mode patterns |
| CN103234637A (en) * | 2013-04-04 | 2013-08-07 | 大连理工大学 | Atmospheric radiation spectrum polarization information measuring system |
| CN106767766A (en) * | 2016-11-22 | 2017-05-31 | 北京航空航天大学 | A kind of sky based on single neutrality point model polarizes localization method |
| CN108267228A (en) * | 2016-12-30 | 2018-07-10 | 南京理工大学 | A kind of atmospheric scattering light polarization analysis method based on CE318 heliographs |
| CN108759819A (en) * | 2018-06-08 | 2018-11-06 | 北京航空航天大学 | A kind of polarization navigation real-time location method based on omnimax polarization degree information |
| CN109459015A (en) * | 2018-11-09 | 2019-03-12 | 北京航空航天大学 | A kind of global autonomic positioning method of polarization navigation based on the observation of maximum degree of polarization |
| CN114076594A (en) * | 2020-08-21 | 2022-02-22 | 大连理工大学 | Distributed polarized light global positioning measurement system |
| CN113532419A (en) * | 2021-06-23 | 2021-10-22 | 合肥工业大学 | Sky polarization mode information acquisition method and device, electronic equipment and storage medium |
| CN113280785A (en) * | 2021-07-22 | 2021-08-20 | 中国人民解放军国防科技大学 | Sky polarized light orientation method based on polarization mode consistency |
Non-Patent Citations (1)
| Title |
|---|
| 赵开春，褚金奎，姚弘轶，张强，王体昌: "Rayleigh大气天空光偏振分布仿真与预测", 第四届中国软件工程大会, vol. 39, pages 287 - 291 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115200543B (en) | 2023-12-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104154928B (en) | Installation error calibrating method applicable to built-in star sensor of inertial platform | |
| CN104913780B (en) | The high-precision deviation of plumb line method for fast measuring of integrated GNSS and CCD zenith telescopes | |
| CN113469896B (en) | Method for improving geometric correction precision of geosynchronous orbit satellite earth observation image | |
| CN105444778A (en) | Star sensor in-orbit attitude determination error obtaining method based on imaging geometric inversion | |
| CN105203103B (en) | The method for real-time measurement of the relatively geographical position relation of ground spacecraft | |
| WO2019062166A1 (en) | Method for automatic geometric correction of cross-platform moon-based earth observation imaging | |
| CN104567869A (en) | Method and device for determining local geographic azimuth and orientation of user with sun position | |
| CN109541647B (en) | GNSS multi-path effect correction method based on semi-celestial sphere grid point model | |
| CN112066979A (en) | Polarization pose information coupling iteration autonomous navigation positioning method | |
| CN105783944B (en) | sun sensor calibration method and system | |
| CN109470237A (en) | One kind is based on polarised light and earth magnetism integrated navigation attitude measurement method | |
| CN102506894A (en) | Stationary base platform absolute space attitude reference establishing method based on precise timing | |
| Wang et al. | A bioinspired navigation system for multirotor UAV by integrating polarization compass/magnetometer/INS/GNSS | |
| CN108896038B (en) | Imaging method of lunar-based optical sensor | |
| CN114001756A (en) | Outfield ground star finding method of small-view-field star sensor | |
| CN108645410B (en) | Ship astronomical positioning method without elevation angle observation | |
| CN105043413B (en) | A kind of bearing calibration of the position of sun based on geoid feature and atmospheric parameter | |
| CN115200543A (en) | Earth sky light polarization basic field observation method based on atmospheric scattering | |
| CN111914396A (en) | High-resolution DEM data-based quick parameterization method for sub-grid terrain three-dimensional surface solar radiation forced effect | |
| CN104458653B (en) | Method and system for measuring atmospheric refraction value at large zenith distance | |
| CN115469340B (en) | Method for estimating number of visible satellites of navigation constellation | |
| KR20150053018A (en) | Visualization method of numerical weather prediction model data on six-panels grid frame and hardware device performing the same | |
| CN115218861A (en) | Astronomical azimuth measurement method based on automatic theodolite | |
| Yakubu et al. | Ramification of datum and ellipsoidal parameters on post processed differential global positioning system (DGPS) data–A case study | |
| CN108955625B (en) | Method and system for acquiring sun position information |
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 |