CN102564574B - Method for measuring radiant illumination of earth albedo - Google Patents
Method for measuring radiant illumination of earth albedo Download PDFInfo
- Publication number
- CN102564574B CN102564574B CN201110460940.5A CN201110460940A CN102564574B CN 102564574 B CN102564574 B CN 102564574B CN 201110460940 A CN201110460940 A CN 201110460940A CN 102564574 B CN102564574 B CN 102564574B
- Authority
- CN
- China
- Prior art keywords
- earth
- sun
- star sensor
- bin
- coordinate system
- 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
Landscapes
- Navigation (AREA)
Abstract
The invention discloses a method for measuring the radiant illumination of earth albedo, and the method comprises the following steps of: establishing an earth coordinate system by using an earth core as a coordinate original point, and determining the longitude and latitude positions of the earth coordinate system of a satellite and the sun at a special moment according to the satellite orbit parameters and the mounting positions of a star sensor on the satellite; dividing the surface of an earth in a longitude and latitude manner, selecting a small surface element on the earth surface, and judging whether the small surface element can be irradiated by the sun and further be viewed by a light screening cover of the star sensor simultaneously or not; and if the small surface element can be irradiated by the sun and further be viewed by the light screening cover of the star sensor simultaneously, calculating stray light influence of sunlight scattered by the small surface element on the star sensor according to a radiation transmission relation of the sun, the surface elements and the star sensor. Comprehensive influences of the earth surface on the star sensor can be obtained by traversing each surface element of the earth surface and overlapping the stray light influences of each surface element on the star sensor. The method disclosed by the invention has the advantages of fast calculation speed and high measurement accuracy.
Description
Technical field
The present invention relates to a kind of method of measuring radiant illumination of earth albedo, belong to Star-Sensor Design technical field.
background technology
Star sensor has that precision is high, integrated level high, and three-axis attitude information can be provided, and this makes it become gradually important attitude measurement equipment in the spacecrafts such as satellite or guided missile.Star sensor is measured satellite with respect to the three-axis attitude of celestial coordinate system by responsive stellar radiation.Stellar radiation generally a little less than, the parasitic light in space environment can form and disturb the faint starlight of star sensor sensitivity, reduces the image quality of system, when serious, can flood fixed star light intensity, and star sensor cannot normally be worked.Therefore determine star sensor veiling glare impact in-orbit, guarantee that its normal work becomes the problem that must consider in Star-Sensor Design.
Star sensor source of stray light in-orbit mainly contains sunshine, moon reflection of light light, earth light and satellite surface reflected light etc.Sunshine subtended angle is less, can be by the quick shadow region that is arranged on of star is hidden; The relative sunshine brightness of moonlight reflection of light light is much lower, and subtended angle is less, therefore less on star sensor impact; Celestial body surface scattering light can reduce celestial body surface scattering by the method such as turned black in celestial body surface; Earth light have scope wide, affect the features such as large, though can carry out part by the method such as star sensor optical axis is faced upward, evade, the impact in-orbit time is still very difficult on star sensor to avoid earth light completely.
Existing Stray light analysis software is as many in simulation analysis such as TracePro, Lighttools, ASAP based on Monte-Carlo technology.Monte-Carlo technology is a kind of random device of calculating based on random number.Monte-Carlo technology in Stray light analysis needs just can obtain the veiling glare impact that earth light produces after a large amount of earth light light of trace.For the star sensor in space flight, earth light region area is very big, obtain earth light impact accurately, needs a large amount of earth light light of trace, and this has all proposed great requirement to existing Stray light analysis software and computer hardware; When Geometrical Optics number is less, the correctness of its result of calculation and precision are difficult to be guaranteed.And adopt existing software simulation, the sun while being difficult to star sensor in-orbit and the relative position of satellite, the surface albedo in the corresponding geographic position of earth surface combine.This has brought great difficulty to the veiling glare impact of star sensor to analysis earth light.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiency of the existing analytical approach based on Monte-Carlo technology, a kind of method of measuring earth light irradiation illumination is provided, have computing velocity fast, the advantage that accuracy of measurement is high.
Technical solution of the present invention is: a kind of method of measuring radiant illumination of earth albedo, and step is as follows:
(1) take earth the earth's core as true origin, set up terrestrial coordinate system, the longitude and latitude position alpha of the sun in terrestrial coordinate system
sun, β
sun, according to satellite orbit parameter, determine the longitude and latitude position alpha of satellite in terrestrial coordinate system
sat, β
satand the star sensor optical axis in terrestrial coordinate system points to v
st; α wherein
sunfor the longitude station of the sun in terrestrial coordinate system, β
sunfor the Position Latitude of the sun in terrestrial coordinate system, α
satfor the longitude station of satellite in terrestrial coordinate system, β
satfor the Position Latitude of satellite in terrestrial coordinate system;
(2) earth surface is carried out to longitude and latitude and divide M earth's surface bin ds of formation, M >=10000;
(3) an optional earth's surface bin ds, preserves the residing longitude and latitude α of this earth's surface bin ds
ds, β
ds, calculate the angle of local normal N and the sun incident vector S of described earth's surface bin ds
connect any one earth's surface bin ds and star sensor light shield entrance center and form the quick connecting line of earth's surface-star, calculate the local normal N of each earth's surface bin ds and the angle of the quick connecting line of earth's surface-star
the length of the quick connecting line of earth's surface-star is l
i,, i=1,2...M-1, M, l
ilongitude and latitude position alpha by satellite in terrestrial coordinate system
sat, β
satwith each earth's surface bin ds location positioning; Calculate the quick connecting line of any one earth's surface-star and described star sensor optical axis and point to v
stangle
(4) judge whether each earth's surface bin ds can be shined upon, the angle of the local normal N of local surface element ds and sun incident vector S
while being less than or equal to 90 °, this earth's surface bin ds can be shined upon, otherwise this earth's surface bin ds do not shined upon, wherein sun incident vector S by the sun longitude and latitude position alpha in terrestrial coordinate system
sun, β
sundetermine;
(5) judge each earth's surface bin ds whether can by star sensor light shield depending on seeing, the local normal N of local surface element ds and the angle of the quick connecting line of earth's surface-star
while being less than or equal to 90 ° this earth's surface bin ds can by star sensor light shield depending on seeing, otherwise this earth's surface bin ds not by star sensor light shield depending on seeing;
(6) judge whether each earth's surface bin ds can bring light disturbance to star sensor, if earth's surface bin ds was both shined upon, again by star sensor light shield depending on seeing, the earth light of this earth's surface bin ds can bring light disturbance to star sensor, otherwise the earth light of this earth's surface bin ds can not bring light disturbance to star sensor;
(7) according to the Reflectivity for Growing Season ρ of each earth's surface bin ds, calculate the radiance L of this earth's surface bin ds
i, suppose that the earth's surface reflection of light is Lambertian radiator, the radiance of each earth's surface bin ds is
e wherein
sunfor the irradiance of the sun outside earth atmosphere, C is the ratio of emittance in solar radiation energy in star sensor work spectral coverage and the full spectral coverage of the sun;
(8) according to the radiance L of each earth's surface bin ds
iobtain each earth's surface bin ds at the radiance dE of light shield porch
i:
wherein dA is the area of earth's surface bin ds;
(9) the radiance dE in light shield porch according to each earth's surface bin ds
iobtain each earth's surface bin ds radiant illumination dE in star sensor image planes
o:
the point-source transmisson ratio that wherein PST is star sensor;
(10) to all earth's surface bin ds, the radiant illumination in star sensor image planes is sued for peace, and can measure the radiant illumination of earth albedo E at star sensor image planes place
o.
The present invention's beneficial effect is compared with prior art:
(1) the method is with respect to the existing Stray light analysis method based on Monte-Carlo technology, and computing velocity is fast, and accuracy is high; The abundant light of Monte-Carlo Technology Need trace just can reach certain accuracy, for spacecraft, earth light district generally have area large, affect the features such as extensive, it is carried out to Stray light analysis, Stray light analysis method based on Monte-Carlo technology needs a large amount of light of emulation, makes its calculated amount be how much multiples and increases progressively, and computing is more time-consuming, the present invention can reach higher accuracy of measurement by longitude and latitude grid is suitably divided in earth's surface, and computing velocity is very fast.
(2) in existing analysis, hypothesis surface albedo is even mostly, and this and actual conditions differ greatly sometimes.The present invention can analyze complicated ground surface environment according to satellite orbit analysis software and Reflectivity for Growing Season database interface, and analysis result is closer to truth; The present invention is not only applicable to calculate the veiling glare impact of earth light on star sensor, is applicable to calculate the analytical calculation of earth light to other kinds space camera yet.
Accompanying drawing explanation
Fig. 1 is measuring principle figure of the present invention.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is further described in detail.
The principle that realizes of the present invention is: take earth the earth's core as true origin, set up terrestrial coordinate system, the installation site on star according to satellite orbit parameter and star sensor, determines satellite and the sun longitude and latitude position in terrestrial coordinate system in particular moment; Earth surface is carried out to longitude and latitude division, choose the little bin in earth's surface, judge that whether this little bin can be shined upon simultaneously but can by star sensor light shield depending on seeing, if the sunshine that can be calculates this little bin scattering by the radiation transmission relation of the sun, bin and star sensor affects the veiling glare of star sensor.Each bin of traversal earth's surface, and each bin is superposeed on the veiling glare impact of star sensor, the combined influence of earth's surface bin to star sensor can be obtained.
As shown in Figure 1, specific implementation process of the present invention is:
(1) take earth the earth's core as true origin, set up terrestrial coordinate system, the longitude and latitude position alpha of the sun in terrestrial coordinate system
sun, β
sun, according to satellite orbit parameter, determine the longitude and latitude position alpha of satellite in terrestrial coordinate system
sat, β
satand the star sensor optical axis in terrestrial coordinate system points to v
st; α wherein
sunfor the longitude station of the sun in terrestrial coordinate system, β
sunfor the Position Latitude of the sun in terrestrial coordinate system, α
satfor the longitude station of satellite in terrestrial coordinate system, β
satfor the Position Latitude of satellite in terrestrial coordinate system;
(2) earth surface is carried out to longitude and latitude and divide M earth's surface bin ds of formation, M >=10000;
(3) an optional earth's surface bin ds, preserves the residing longitude and latitude α of this earth's surface bin ds
ds, β
ds, calculate the angle of local normal N and the sun incident vector S of described earth's surface bin ds
connect any one earth's surface bin ds and star sensor light shield entrance center and form the quick connecting line of earth's surface-star, calculate the local normal N of each earth's surface bin ds and the angle of the quick connecting line of earth's surface-star
the length of the quick connecting line of earth's surface-star is l
i,, i=1,2...M-1, M, l
ilongitude and latitude position alpha by satellite in terrestrial coordinate system
sat, β
satwith each earth's surface bin ds location positioning; Calculate the quick connecting line of any one earth's surface-star and described star sensor optical axis and point to v
stangle
(4) judge whether each earth's surface bin ds can be shined upon, the angle of the local normal N of local surface element ds and sun incident vector S
while being less than or equal to 90 °, this earth's surface bin ds can be shined upon, otherwise this earth's surface bin ds do not shined upon, wherein sun incident vector S by the sun longitude and latitude position alpha in terrestrial coordinate system
sun, β
sundetermine;
(5) judge each earth's surface bin ds whether can by star sensor light shield depending on seeing, the local normal N of local surface element ds and the angle of the quick connecting line of earth's surface-star
while being less than or equal to 90 ° this earth's surface bin ds can by star sensor light shield depending on seeing, otherwise this earth's surface bin ds not by star sensor light shield depending on seeing;
(6) judge whether each earth's surface bin ds can bring light disturbance to star sensor, if earth's surface bin ds was both shined upon, again by star sensor light shield depending on seeing, the earth light of this earth's surface bin ds can bring light disturbance to star sensor, otherwise the earth light of this earth's surface bin ds can not bring light disturbance to star sensor;
(7) according to the Reflectivity for Growing Season ρ of each earth's surface bin ds, calculate the radiance L of this earth's surface bin ds
i, suppose that the earth's surface reflection of light is Lambertian radiator, the radiance of each earth's surface bin ds is
e wherein
sunfor the irradiance of the sun outside earth atmosphere, C is the ratio of emittance in solar radiation energy in star sensor work spectral coverage and the full spectral coverage of the sun;
(8) according to the radiance L of each earth's surface bin ds
iobtain each earth's surface bin ds at the radiance dE of light shield porch
i:
wherein dA is the area of earth's surface bin ds;
(9) the radiance dE in light shield porch according to each earth's surface bin ds
iobtain each earth's surface bin ds radiant illumination dE in star sensor image planes
o:
the point-source transmisson ratio that wherein PST is star sensor;
(10) to all earth's surface bin ds, the radiant illumination in star sensor image planes is sued for peace, and can measure the radiant illumination of earth albedo E at star sensor image planes place
o.
Wherein the point-source transmisson of star sensor is star sensor illuminance of image plane and the ratio of illumination perpendicular on the incident light direction of propagation than (PST (θ)), and this incident light and optical axis included angle are θ; The PST curve of star sensor can obtain by star sensor is carried out to veiling glare trace.The present invention can quantitatively calculate the veiling glare impact of the earth light of star sensor when multiple orbital attitudes and installation position on star sensor, for the use in-orbit of star sensor provides foundation, has good using value in practical engineering application.
The content that the present invention is not described in detail belongs to professional and technical personnel in the field's known technology.
Claims (1)
1. measure a method for radiant illumination of earth albedo, it is characterized in that step is as follows:
(1) take earth the earth's core as true origin, set up terrestrial coordinate system, the longitude and latitude position alpha of the sun in terrestrial coordinate system
sun, β
sun, according to satellite orbit parameter, determine the longitude and latitude position alpha of satellite in terrestrial coordinate system
sat, β
satand the star sensor optical axis in terrestrial coordinate system points to v
st; α wherein
sunfor the longitude station of the sun in terrestrial coordinate system, β
sunfor the Position Latitude of the sun in terrestrial coordinate system, α
satfor the longitude station of satellite in terrestrial coordinate system, β
satfor the Position Latitude of satellite in terrestrial coordinate system;
(2) earth surface is carried out to longitude and latitude and divide M earth's surface bin ds of formation, M >=10000;
(3) an optional earth's surface bin ds, preserves the residing longitude and latitude α of this earth's surface bin ds
ds, β
ds, calculate the angle of local normal N and the sun incident vector S of described earth's surface bin ds
connect any one earth's surface bin ds and star sensor light shield entrance center and form the quick connecting line of earth's surface-star, calculate the local normal N of each earth's surface bin ds and the angle of the quick connecting line of earth's surface-star
the length of the quick connecting line of earth's surface-star is l
i,, i=1,2...M-1, M, l
ilongitude and latitude position alpha by satellite in terrestrial coordinate system
sat, β
satwith each earth's surface bin ds location positioning; Calculate the quick connecting line of any one earth's surface-star and described star sensor optical axis and point to v
stangle
(4) judge whether each earth's surface bin ds can be shined upon, the angle of the local normal N of local surface element ds and sun incident vector S
while being less than or equal to 90 °, this earth's surface bin ds can be shined upon, otherwise this earth's surface bin ds do not shined upon, wherein sun incident vector S by the sun longitude and latitude position alpha in terrestrial coordinate system
sun, β
sundetermine;
(5) judge each earth's surface bin ds whether can by star sensor light shield depending on seeing, the local normal N of local surface element ds and the angle of the quick connecting line of earth's surface-star
while being less than or equal to 90 ° this earth's surface bin ds can by star sensor light shield depending on seeing, otherwise this earth's surface bin ds not by star sensor light shield depending on seeing;
(6) judge whether each earth's surface bin ds can bring light disturbance to star sensor, if earth's surface bin ds was both shined upon, again by star sensor light shield depending on seeing, the earth light of this earth's surface bin ds can bring light disturbance to star sensor, otherwise the earth light of this earth's surface bin ds can not bring light disturbance to star sensor;
(7) according to the Reflectivity for Growing Season ρ of each earth's surface bin ds, calculate the radiance L of this earth's surface bin ds
i, suppose that the earth's surface reflection of light is Lambertian radiator, the radiance of each earth's surface bin ds is
e wherein
sunfor the irradiance of the sun outside earth atmosphere, C is the ratio of emittance in solar radiation energy in star sensor work spectral coverage and the full spectral coverage of the sun;
(8) according to the radiance L of each earth's surface bin ds
iobtain each earth's surface bin ds at the radiance dE of light shield porch
i:
wherein dA is the area of earth's surface bin ds;
(9) the radiance dE in light shield porch according to each earth's surface bin ds
iobtain each earth's surface bin ds radiant illumination dE in star sensor image planes
o:
the point-source transmisson ratio that wherein PST is star sensor;
(10) to all earth's surface bin ds, the radiant illumination in star sensor image planes is sued for peace, and can measure the radiant illumination of earth albedo E at star sensor image planes place
o.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110460940.5A CN102564574B (en) | 2011-12-29 | 2011-12-29 | Method for measuring radiant illumination of earth albedo |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110460940.5A CN102564574B (en) | 2011-12-29 | 2011-12-29 | Method for measuring radiant illumination of earth albedo |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102564574A CN102564574A (en) | 2012-07-11 |
| CN102564574B true CN102564574B (en) | 2014-04-02 |
Family
ID=46410603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110460940.5A Active CN102564574B (en) | 2011-12-29 | 2011-12-29 | Method for measuring radiant illumination of earth albedo |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102564574B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928075B (en) * | 2012-10-25 | 2015-04-08 | 西安工业大学 | Method for analyzing illumination of stray light of low-rail space camera by earth-atmosphere stray light |
| CN107543603B (en) * | 2016-06-27 | 2019-10-18 | 广州供电局有限公司 | Real-time solar radiation quantity measuring method and system |
| CN107238435B (en) * | 2017-06-07 | 2018-06-12 | 江苏省无线电科学研究所有限公司 | Solar radiation flux Vertical Profile survey meter |
| CN108827279B (en) * | 2018-04-11 | 2021-10-26 | 西北工业大学 | Simulated sun sensor anti-interference method based on earth albedo |
| CN113361163B (en) * | 2021-06-03 | 2022-09-30 | 中国人民解放军战略支援部队航天工程大学 | Satellite attitude estimation method for correcting earth reflected light |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01196509A (en) * | 1988-02-01 | 1989-08-08 | Nec Corp | Polyhedral type satellite attitude detector |
| JPH1196509A (en) * | 1997-09-26 | 1999-04-09 | Hitachi Ltd | Magnetic disk device and semiconductor device used for it |
| CN100492050C (en) * | 2007-02-09 | 2009-05-27 | 上海大学 | Ground reflectivity distribution map based remote sensing inversion explanation of underground latency construction |
| CN101089553A (en) * | 2007-07-06 | 2007-12-19 | 中国科学院上海技术物理研究所 | Method of preventing earth light for digital solar sensor |
-
2011
- 2011-12-29 CN CN201110460940.5A patent/CN102564574B/en active Active
Non-Patent Citations (3)
| Title |
|---|
| 李道勇,王云强,宫彦军等.空间目标的光散射特性研究.《烟台大学学报（自然科学与工程版）》》.2004,第17卷(第3期),第183-187页. * |
| 王磊, 张涛, 崔文楠等.空间目标和背景的辐照特性仿真技术研究.《科学技术与工程》.2010,第10卷(第15期),3610-3614. |
| 空间目标和背景的辐照特性仿真技术研究;王磊, 张涛, 崔文楠等;《科学技术与工程》;20100531;第10卷(第15期);第3610-3614页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102564574A (en) | 2012-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Calcabrini et al. | A simplified skyline-based method for estimating the annual solar energy potential in urban environments | |
| Adams Jr et al. | An evaluation of the exposure in nadir observation of the JEM-EUSO mission | |
| CN102564574B (en) | Method for measuring radiant illumination of earth albedo | |
| CN104880702A (en) | Method and device for on-orbit absolute radiation calibration | |
| CN104573251A (en) | Method for determining full-field-of-view apparent spectral radiance of satellite-borne optical remote sensor | |
| Cognion | Observations and modeling of GEO satellites at large phase angles | |
| Rezk et al. | Experimental implementation of meteorological data and photovoltaic solar radiation monitoring system | |
| Melo et al. | Development of a closed and open loop solar tracker technology | |
| CN103852078A (en) | Device and method for measuring stray light protection angle of space optical attitude sensor | |
| Alam et al. | Shadow effect on photovoltaic potentiality analysis using 3D city models | |
| Kaynak et al. | A software tool development study for solar energy potential analysis | |
| CN104101297A (en) | Space object dimension acquisition method based on photoelectric observation | |
| CN109375650A (en) | Different roofing situation barrier shadow occlusion calculation methods | |
| Wong et al. | A feasibility study of the 1.5-axis tracking model in utility-scale solar PV plants | |
| Zhang et al. | Influence of incidence angle effects on the performance of bifacial photovoltaic modules considering rear-side reflection | |
| Cellura et al. | A photographic method to estimate the shading effect of obstructions | |
| CN106452301A (en) | Photovoltaic component installation design method for non-due-south slope | |
| CN103411896B (en) | Low rail space junk Extra-low Dispersion spectral signature acquisition methods | |
| Xu et al. | Stellar radiation modeling and image simulation for airborne daytime star sensor | |
| Yousefian et al. | Fabrication, characterization, and error mitigation of non-flat sun sensor | |
| Wei et al. | Stray light from galactic sky and zodiacal light for JWST | |
| Ivanova | Use of fisheye orthographic projection for cumulative estimation of direct and diffuse irradiation under partially obstructed sky | |
| ur Rehman et al. | A novel methodology for determining sky blocking by obstacles viewed virtually from any location on site | |
| CN112417670B (en) | GEO target photometric characteristic calculation model considering windsurfing board offset effect | |
| Wu et al. | Light scattering and visibility condition of spatial objects |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |