CN117851717A - Moon north-south illumination condition real-time calculation method considering high-precision topography - Google Patents
Moon north-south illumination condition real-time calculation method considering high-precision topography Download PDFInfo
- Publication number
- CN117851717A CN117851717A CN202311060000.6A CN202311060000A CN117851717A CN 117851717 A CN117851717 A CN 117851717A CN 202311060000 A CN202311060000 A CN 202311060000A CN 117851717 A CN117851717 A CN 117851717A
- Authority
- CN
- China
- Prior art keywords
- illumination
- measurement position
- sun
- ground measurement
- terrain
- 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.)
- Pending
Links
- 238000005286 illumination Methods 0.000 title claims abstract description 139
- 238000004364 calculation method Methods 0.000 title claims abstract description 49
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 title claims abstract description 39
- 238000012876 topography Methods 0.000 title claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 97
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000000873 masking effect Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to the technical field of aerospace, and provides a method for calculating the illumination condition of the north and south poles of a moon in real time by considering high-precision topography, in particular to the detection of the surface of the north and south poles, which comprises the following steps: s1: establishing a coordinate system of a ground measurement position, inputting the coordinate of the ground measurement position by an illumination calculation model, and defining the illumination state of the sun so as to illustrate the illumination condition of the moon when taking high-precision terrain into consideration; s2: calculating the maximum height angle of the ground measurement position according to the calculation parameters of the illumination condition and the appointed topographic database information; s3: and calculating sun related parameters through the high-precision ephemeris and the specified time sequence, and further calculating the illumination condition of the ground measurement position. Based on the current highest-precision DEM topographic data, the horizon is searched for and the elevation angle is calculated by azimuth angle, then the illumination calculation result of the place is obtained by combining the sun elevation angle and the disc viewing diameter, the illumination condition of the current position is calculated, and the table data and the graphic data are output through a user interface.
Description
Technical Field
The invention relates to the technical field of aerospace, in particular to a real-time calculation method for the illumination conditions of the north and south of the moon in consideration of high-precision topography.
Background
In the future lunar exploration mission planning process, especially lunar exploration mission in the north-south polar region, the illumination condition of the expected exploration area of the north-south polar region needs to be analyzed.
At present, the calculation of the lunar illumination condition often ignores the shielding of lunar topography fluctuation on the sun, and the lunar equivalent is calculated as an ellipsoid, so that the illumination result influenced by high topography space resolution is lacked. In the north-south regions of the moon, the result error obtained by the current calculation method is large, and the method is difficult to reasonably apply in engineering design and scientific research practice.
In the north-south regions of the moon, the result error obtained by the current calculation method is large, and the method is difficult to reasonably apply in engineering design and scientific research practice.
Disclosure of Invention
The invention aims to provide a real-time calculation method for the illumination conditions of the north and south of the moon in consideration of high-precision topography.
In order to achieve the above purpose, the invention provides a method for calculating the illumination condition of the north and south of the moon in real time by considering high-precision topography, which comprises the following steps:
s1: establishing a coordinate system of the ground measurement position according to the lunar ground measurement position coordinate, inputting the coordinate of the ground measurement position by an illumination calculation model, and describing the illumination condition of the ground measurement position when taking high-precision terrain into consideration by defining the illumination state of the sun by the illumination calculation model;
s2: calculating the maximum height angle of the ground measurement position according to the calculation parameters of the illumination condition of the ground measurement position and the appointed topographic database information, namely calculating the topographic shade condition of the ground measurement position;
s3: after calculating the topography shielding condition of the ground measurement position, calculating sun related parameters through high-precision ephemeris and a specified time sequence, and calculating the illumination state according to the sun related parameters, so that the illumination condition of the ground measurement position is further calculated, and table data and graphic data are output through a user interface by the topography shielding condition and the illumination condition.
Preferably, in step S1, a coordinate system of the ground measurement position is established, specifically:
using a cartesian coordinate system, inputting longitude and latitude coordinates of the ground measurement position, establishing a coordinate system of the ground measurement position, representing a position of maximum height of the sun and surrounding terrain of the ground measurement position.
Preferably, the illumination state of the sun is defined, specifically:
defining a sun, namely an illumination state, namely defining the proportion of the sun shielded by the terrain, dividing the illumination state into three states of illumination, penumbra and full shadow according to the size of the sun shielded part, and further describing the moon night, moon day or penumbra time of the ground measurement position according to the condition of the illumination state representing the moon night, moon day or penumbra of the ground measurement position, and further describing the illumination condition of the ground measurement position when the high-precision terrain is considered.
Preferably, the illumination states are divided into three states of illumination, penumbra and full-image, specifically:
illumination means that the sun is not shielded at all and is in a complete illumination state; penumbra means that the sun is partially blocked, in a penumbra state; full shadow means that the sun is completely blocked, in a fully shaded state.
Preferably, in step S2, the calculating parameters further includes: the method comprises the steps of calculating parameters such as terrain searching precision, terrain searching range, illumination calculation span, illumination calculation interval and the like.
Preferably, in step S2, calculating the terrain mask condition at the ground measurement position includes:
according to the ground measurement position, obtaining high-precision terrain data of an area around the ground measurement position, and then calculating terrain mask elevation parameters, namely terrain mask parameters, of 360-degree azimuth of the ground measurement position at intervals of 0.5 degrees or 1 degree by taking the north direction of the ground measurement position as an azimuth zero point;
and obtaining a terrain mask range of the ground measurement position according to the terrain mask parameters, simultaneously carrying out solar azimuth interpolation according to a solar azimuth, accurately obtaining a terrain elevation parameter of a corresponding azimuth, and further judging the illumination state of the sun.
Preferably, in step S3, the solar-related parameters include: searching for a terrain maximum elevation angle alpha in the illuminated state at a specified time based on the position of the maximum height of the terrain surrounding the ground measurement position t Calculating the sun altitude angle El, the sun azimuth angle Az and the sun apparent radius alpha s The direction of the sun position and the direction alpha of the maximum elevation angle of the topography of the ground measurement position t Is included angle alpha g Obtaining the illumination condition at the current moment, wherein in a Cartesian coordinate system, the maximum elevation angle alpha of the terrain can be calculated through tangent values t Said maximum elevation angle alpha t I.e. the elevation angle of the horizon.
Preferably, the calculation of the solar altitude angle El and the solar azimuth angle Az is specifically:
solar coordinates of a coordinate system position at the ground measurement position at a specified momentWherein X is s For measuring the abscissa of the position of the sun on said ground, Y s For measuring the ordinate, Z, of the position of the sun on said ground s For the vertical coordinate of the sun at the ground measurement position, T sets the sun coordinate +.>Transpose, solar coordinates->Can be expressed as:
according to the solar coordinatesThe solar altitude angle El and the solar azimuth angle Az may be calculated according to the following formulas:
preferably, the solar viewing radius alpha s Specifically, the calculation of (a) is as follows:
wherein R is o Is the solar radius.
The illumination state at the appointed moment is used for judging the illumination condition of the ground measurement position according to the sun illumination state:
preferably, in step S3, the topographic mask condition and the illumination condition output table data and graphic data through a user interface, specifically:
establishing the user interface for calculating the terrain mask condition and the illumination condition, and inputting input information including coordinates of the ground measurement position, the terrain database information and the calculation parameters on the user interface;
the user interface is hosted on a server in the form of a service, and after the server receives the input information, the server returns a result to the user interface through the topographic masking condition algorithm and the illumination condition algorithm.
Compared with the prior art, the invention provides a real-time calculation method for the illumination condition of the north and south of the moon in consideration of high-precision topography, realizes accurate calculation of illumination results, can customize the searching range and precision of topography, searches horizon and calculates elevation angles by azimuth angle based on the currently existing highest-precision DEM topography data, and combines solar altitude and circular viewing diameter to obtain illumination calculation results of the spot, and outputs form data and graphic data through a user interface.
Drawings
FIG. 1 is a schematic diagram of a horizontal coordinate system of a ground measurement position;
FIG. 2 is a schematic diagram of a local terrain mask attribute setup page;
FIG. 3 is a schematic view of an illumination calculation model of a ground measurement location;
FIG. 4 is a schematic view of solar azimuth and altitude in the horizontal coordinate system of the ground measurement position;
FIG. 5 is a schematic view of a three-dimensional graphics setup page of a terrain mask;
FIG. 6 is a graphical representation of the terrain mask output for the bottom of the moon antarctic shakelton pit;
FIG. 7 is a diagram of a results output settings page for local reports, charts, etc.;
FIG. 8 is a diagram of local terrain mask data;
FIG. 9 is a schematic illustration of the illumination period at point P7 at the top of the south pole pit of the moon;
FIG. 10 is a graph showing parameters such as P7 point illumination intensity;
fig. 11 is a schematic view of solar azimuth and altitude reporting.
Reference numerals:
1: a zenith direction ellipsoid normal;
2; a north pole;
3: a ground wire;
4: a ground measurement location;
5: a local level;
6: moon;
7: a moon surface south pole pit vertex;
8: a sun disk;
Detailed Description
For the purposes of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments, and that a skilled person may use different methods for each step to implement the described functions, but such implementation should not be considered to be beyond the scope of the present invention. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless expressly stated otherwise, as understood by those skilled in the art. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
Examples
The invention aims to provide a real-time calculation method for the illumination conditions of the north and south of the moon in consideration of high-precision topography.
In order to achieve the above purpose, the invention provides a method for calculating the illumination condition of the north and south of the moon in real time by considering high-precision topography, which comprises the following steps:
s1: establishing a coordinate system of a ground measurement position according to the coordinate of the ground measurement position of the moon, inputting the coordinate of the ground measurement position by an illumination calculation model, and defining the illumination state of the sun by the illumination calculation model so as to illustrate the illumination condition of the ground measurement position when the high-precision terrain is considered, specifically, in the embodiment, taking the vertex of a south pole pit of the moon as an example, calculating the illumination condition of the south pole, and outputting a related chart;
s2: calculating the maximum height angle of the ground measurement position according to the calculation parameters of the illumination condition of the ground measurement position and the appointed topographic database information, namely calculating the topographic shade condition of the ground measurement position;
s3: after calculating the topography shielding condition of the ground measurement position, calculating sun related parameters through a high-precision ephemeris and a specified time sequence, and calculating an illumination state according to the sun related parameters, so as to further calculate the illumination condition of the ground measurement position, wherein the topography shielding condition and the illumination condition can output form data and graphic data through a user interface, and specifically, in the embodiment, the user interface mainly relates to calculation service of the topography shielding condition and calculation service of the illumination condition.
Preferably, in step S1, a coordinate system of the ground measurement position is established, specifically:
referring to fig. 1, using a cartesian coordinate system, longitude and latitude coordinates of a ground measurement position are input, a coordinate system of the ground measurement position is established, and the maximum height position of the sun and surrounding terrain of the ground measurement position is represented, specifically, in this embodiment, when the sun is observed at the vertex of a south pole pit of a lunar surface, a local horizontal coordinate system (also called north east coordinate system) is adopted, and an LH coordinate system is defined as follows:
1) The X axis points to the local north direction;
2) The Y axis points to the local east;
3) Z-axis pointing in the direction of the local nadir
Where the XY plane is the local horizontal plane, perpendicular to the earth's ellipsoidal normal, specifically, in this embodiment, the coordinate information of the vertex of the input south pole pit is transmitted to the user interface.
Preferably, the illumination state of the sun is defined, specifically:
defining the sun as the illumination state, namely defining the proportion of the sun shielded by the terrain, dividing the illumination state into three states of illumination, penumbra and full shadow according to the size of the shielded part of the sun, and further describing the moon night, moon day or penumbra time of the ground measurement position according to the condition of the moon night, moon day or penumbra of the ground measurement position, and further describing the illumination condition of the ground measurement position when taking high-precision terrain into consideration.
Preferably, the illumination states are divided into three states of illumination, penumbra and full shadow, specifically:
illumination means that the sun is not shielded at all and is in a complete illumination state; penumbra means that the sun is partially blocked, in the penumbra state; full shadow means that the sun is completely blocked, in a fully shaded state.
Preferably, the illumination state is divided into three states of illumination, penumbra and full shadow, and the month and night, month and day or penumbra time of the ground measurement position is accurately described according to the illumination state of the ground measurement position, specifically, in this embodiment, the sun is not blocked at all during illumination, the local is in a full illumination state, and the sun illumination intensity factor is 1; penumbra means that the sun is partially blocked, and is locally in a penumbra state, and the sun illumination intensity is in the fraction between 0 and 1; full shadow means that the sun is completely blocked, and is in a completely shadow state locally, and the sun illumination intensity factor is 0.
Preferably, in step S2, calculating parameters further includes: the method comprises the steps of calculating parameters such as terrain searching precision, terrain searching range, illumination calculation span, illumination calculation interval and the like.
Preferably, in step S2, calculating the terrain mask condition for the ground measurement location includes:
according to the ground measurement position, obtaining high-precision topographic data of an area around the ground measurement position, then calculating topographic mask elevation parameters, i.e. topographic mask parameters, of 360-degree azimuth of the ground measurement position at intervals of 0.5 degrees or 1 degree by taking the north direction of the ground measurement position as an azimuth zero point, and locally storing the calculated data in a user interface;
and obtaining a terrain masking range of the ground measurement position according to the terrain masking parameter, simultaneously carrying out solar azimuth interpolation according to the solar azimuth, accurately obtaining a terrain elevation parameter of the corresponding azimuth, and further judging the illumination state of the sun.
Referring to fig. 2, the computing service of the terrain mask condition of the application interface specifically opens the object attribute window, selects "terrain data" in the "use" option drop-down frame in the "basic-terrain mask" attribute page, invokes the background algorithm, uses the terrain data to perform terrain mask computation, uses the default parameters to perform terrain mask computation when the "terrain mask parameter" page is not selected by default, uses the default parameters to perform terrain mask computation when the background computation, the step size is 5m, the maximum search distance of the terrain is 180km, the "terrain mask parameter" is selected according to the actual condition of the terrain when the application interface is actually used, then specific parameter setting is performed, and meanwhile, custom two-pole Dem data (moon_ldem_80s_20m data in the example) can be adopted, and the terrain mask condition returned to the current position after computation is stored locally.
Referring to fig. 3, in step S3, the sun-related parameters include: calculating the sun altitude angle El, the sun azimuth Az and the sun disc radius alpha by the terrain altitude angle of the ground measurement position at the designated moment s The direction of the sun position and the direction alpha of the maximum elevation angle of the ground measurement position topography g The illumination condition at the current moment is obtained, specifically, in the embodiment, the maximum elevation angle of the solar direction is alpha t When the topography is not considered (alpha t ＝0)。
Preferably, in step S3, the solar-related parameters include: searching for a maximum elevation angle alpha of the terrain in the illuminated state at a specified time based on the position of the maximum height of the terrain surrounding the ground measurement position t Calculating the sun altitude angle El, the sun azimuth angle Az and the sun apparent radius alpha s The direction of the sun position and the direction alpha of the maximum elevation angle of the terrain at the ground measurement position t Is included angle alpha g Obtaining the illumination condition at the current moment, wherein in a Cartesian coordinate system, the maximum elevation angle alpha of the terrain can be calculated through tangent values t Maximum elevation angle alpha t I.e. the elevation angle of the horizon, in particular, in the present embodiment, the position of the maximum height of the terrain surrounding the ground measurement position is acquired and the information of the height is entered into the user interface.
More preferably, the calculation of the solar altitude angle El and the solar azimuth angle Az is specifically:
solar coordinates of a coordinate system position at a ground measurement position at a specified momentWherein X is s For measuring the abscissa of the position of the sun on the ground, Y s For measuring the ordinate, Z, of the position of the sun on the ground s For the vertical coordinate of the solar measurement position on the ground, T is the solar coordinate +.>Transpose, solar coordinates->Can be expressed as:
according to the solar coordinatesThe solar altitude angle El and the solar azimuth angle Az can be calculated according to the following formula:
specifically, in the present embodiment, R o The value of the sun radius is 695700km, which corresponds to the apparent radius alpha s About 0.27 deg., the specific values vary slightly from the distance of the sun.
Referring to fig. 4, specifically, in the present embodiment, the azimuth angle is defined as: the included angle between the X axis and the projection vector of the sun direction vector in the XY plane is positive, and the +X axis is zero point and the +Y axis direction is positive; the altitude is defined as: the included angle between the solar direction vector and the XY plane, the Z axis direction is positive.
Preferably, the sun viewing radius alpha s Specifically, the calculation of (a) is as follows:
wherein R is o Is the radius of the sun.
Preferably, specifying the time illumination state includes:
and judging the illumination condition of a measurement position on the ground at the south pole according to the illumination state of the sun.
Preferably, in step S3, the topographic mask condition and the illumination condition output table data and graphic data through a user interface, specifically:
establishing the user interface for calculating the terrain mask condition and the illumination condition, and inputting input information including coordinates of the ground measurement position, the terrain database information and the calculation parameters on the user interface;
the user interface is hosted on a server in the form of a service, and after the server receives the input information, the server returns a result to the user interface through the topographic masking condition algorithm and the illumination condition algorithm.
Referring to FIG. 5, in the "three-dimensional graphics-terrain mask" attribute page, "distance settings" may be selected so that the terrain mask parameters are visually displayed in a three-dimensional window.
Referring to fig. 6, the terrain mask parameters display the results.
Referring to fig. 7, after the calculation is completed, a menu bar clicks on a "chart", a "report and statistics chart" window is opened, a chart pattern which can be generated by the present invention is listed in a report pattern on the right side, which includes terrain mask data ("Azimuth-Elevation"), illumination time period ("illumination time"), illumination Intensity parameter ("Solar Intensity"), etc., specifically, the time parameter in the illumination condition and the calculated data item are sent to the calculation service of the illumination condition, the data item includes illumination Intensity, illumination duration time, and after the calculation is completed, the result is returned, and the relevant chart can be custom generated.
As shown in fig. 8, is a local terrain mask data map obtained by an algorithm for the terrain mask conditions provided by the user interface.
As shown in fig. 9, is the illumination period of 7 points of the top of the moon south pole pit.
As shown in FIG. 10, the illumination intensity and other parameter diagrams of the south pole pit vertex 7 of the lunar surface are obtained through the illumination condition algorithm provided by the user interface.
As shown in fig. 11, is a solar azimuth, altitude report.
In summary, the real-time computing method for the illumination conditions of the north and south of the moon provided by the invention can load the custom high-precision topographic data conforming to the GIS standard, greatly expands the compatibility and flexibility of topographic data in illumination computation, and can conveniently and rapidly specify any time span and computation interval by developing a related graphical interface, thereby facilitating the user to compute the illumination conditions at any time locally.
Claims (10)
1. A method for calculating the illumination condition of the north and south of the moon in real time by considering the high-precision topography is characterized by comprising the following steps:
s1: establishing a coordinate system of the ground measurement position according to the lunar ground measurement position coordinate, inputting the coordinate of the ground measurement position by an illumination calculation model, and describing the illumination condition of the ground measurement position when taking high-precision terrain into consideration by defining the illumination state of the sun by the illumination calculation model;
s2: calculating the maximum height angle of the ground measurement position according to the calculation parameters of the illumination condition of the ground measurement position and the appointed topographic database information, namely calculating the topographic shade condition of the ground measurement position;
s3: after calculating the topography shielding condition of the ground measurement position, calculating sun related parameters through high-precision ephemeris and a specified time sequence, and calculating the illumination state according to the sun related parameters, so that the illumination condition of the ground measurement position is further calculated, and table data and graphic data are output through a user interface by the topography shielding condition and the illumination condition.
2. The method for calculating the illumination condition of the north and south of the moon in consideration of high-precision topography according to claim 1, wherein in step S1, a coordinate system of the ground measurement position is established, specifically:
using a cartesian coordinate system, inputting longitude and latitude coordinates of the ground measurement position, establishing a coordinate system of the ground measurement position, representing a position of maximum height of the sun and surrounding terrain of the ground measurement position.
3. The method for real-time calculation of the illumination condition of the north and south of the moon in consideration of high-precision topography according to claim 1, wherein in step S1, the illumination condition of the sun is defined, specifically:
defining a sun, namely an illumination state, namely defining the proportion of the sun shielded by the terrain, dividing the illumination state into three states of illumination, penumbra and full shadow according to the size of the sun shielded part, and further describing the moon night, moon day or penumbra time of the ground measurement position according to the illumination state representing the moon night, moon day or penumbra condition of the ground measurement position and further describing the illumination condition of the ground measurement position when the high-precision terrain is considered.
4. The method for calculating the illumination conditions of the north and south of the moon in real time by considering high-precision topography according to claim 3, wherein the illumination conditions are divided into three states of illumination, penumbra and full shadow, specifically:
illumination means that the sun is not shielded at all and is in a complete illumination state; penumbra means that the sun is partially blocked, in a penumbra state; full shadow means that the sun is completely blocked, in a fully shaded state.
5. The method for real-time calculation of the illumination condition of the north and south of the moon in consideration of high-precision topography according to claim 1, wherein in step S2, the calculation parameters further include: the method comprises the steps of calculating parameters such as terrain searching precision, terrain searching range, illumination calculation span, illumination calculation interval and the like.
6. The method for real-time calculation of the lunar north and south illumination condition taking into account high-precision terrain according to claim 1, wherein in step S2, calculating the terrain mask condition of the ground measurement location comprises:
according to the ground measurement position, obtaining high-precision terrain data of an area around the ground measurement position, and then calculating terrain mask elevation parameters, namely terrain mask parameters, of 360-degree azimuth of the ground measurement position at intervals of 0.5 degrees or 1 degree by taking the north direction of the ground measurement position as an azimuth zero point;
and obtaining a terrain mask range of the ground measurement position according to the terrain mask parameters, and simultaneously carrying out solar azimuth interpolation according to a solar azimuth to obtain a terrain elevation parameter of a corresponding azimuth, so as to further judge the illumination state of the sun.
7. A method for real-time calculation of the north and south lunar lighting conditions taking into account high-precision terrain according to claims 1, 2 and 6, wherein in step S3, the sun-related parameters include:
searching for a maximum elevation angle alpha of the terrain in the illuminated state at a specified time based on the position of the maximum height of the terrain surrounding the ground measurement position t Calculating the sun altitude angle El, the sun azimuth angle Az and the sun apparent radius alpha s A maximum elevation angle alpha of the topography depending on the direction of the sun position and the ground measurement position t The included angle of the direction is alpha g Obtaining the illumination condition at the current moment, wherein in a Cartesian coordinate system, the maximum elevation angle alpha of the terrain can be calculated through tangent values t Said maximum elevation angle alpha t I.e. the elevation angle of the horizon.
8. The method for calculating the illumination condition of the north and south of the moon in real time taking into account the high-precision topography as set forth in claim 7, wherein the calculation of the solar altitude angle El and the solar azimuth angle Az is specifically as follows:
solar coordinates of a coordinate system position at the ground measurement position at a specified momentWherein X is s For measuring the abscissa of the position of the sun on said ground, Y s For the ordinate of the solar measured position on the ground,Z s for the vertical coordinate of the sun at the ground measurement position, T sets the sun coordinate +.>Transpose, solar coordinates->Can be expressed as:
according to the solar coordinatesThe solar altitude angle El and the solar azimuth angle Az may be calculated according to the following formulas:
the sun viewing radius alpha s Wherein R is o The solar radius is specifically:
9. the method for calculating the illumination condition of the north and south of the moon in real time according to claim 4, 7 and 9, wherein the illumination condition judgment basis at the designated time comprises:
and judging the illumination condition of the ground measurement position according to the sun illumination state.
10. The method according to claim 1, wherein in step S3, the topographic mask condition and the illumination condition output table data and graphic data through a user interface, specifically:
establishing the user interface for calculating the terrain mask condition and the illumination condition, and inputting input information including coordinates of the ground measurement position, the terrain database information and the calculation parameters on the user interface;
the user interface is hosted on a server in the form of a service, and after the server receives the input information, the server returns a result to the user interface through the topographic masking condition algorithm and the illumination condition algorithm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311060000.6A CN117851717A (en) | 2023-08-22 | 2023-08-22 | Moon north-south illumination condition real-time calculation method considering high-precision topography |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311060000.6A CN117851717A (en) | 2023-08-22 | 2023-08-22 | Moon north-south illumination condition real-time calculation method considering high-precision topography |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117851717A true CN117851717A (en) | 2024-04-09 |
Family
ID=90535280
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311060000.6A Pending CN117851717A (en) | 2023-08-22 | 2023-08-22 | Moon north-south illumination condition real-time calculation method considering high-precision topography |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117851717A (en) |
-
2023
- 2023-08-22 CN CN202311060000.6A patent/CN117851717A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109671150A (en) | Airport Earthwork calculation method based on digital earth | |
| CN102324740A (en) | Online correction method for interference of DC (Direct Current) power transmission line on geomagnetic observation | |
| US5807113A (en) | Method and apparatus for training in the detection of nuclear, biological and chemical (NBC) contamination | |
| CN111275757A (en) | Pseudo-satellite field simulation layout method based on DEM data processing | |
| CN111121758A (en) | Rapid modeling and credible positioning method for indoor magnetic map | |
| CN110146052B (en) | Plane normal astronomical directional measurement method and system based on total station | |
| CN112859130A (en) | High-precision electronic map position matching method for field navigation patrol | |
| CN111914396B (en) | Sub-grid terrain three-dimensional earth surface solar radiation forced effect rapid parameterization method based on high-resolution DEM data | |
| CN117851717A (en) | Moon north-south illumination condition real-time calculation method considering high-precision topography | |
| Kaňuk et al. | Testing of V3. sun module prototype for solar radiation modelling on 3D objects with complex geometric structure | |
| CN115018973B (en) | Non-target evaluation method for low-altitude unmanned aerial vehicle point cloud modeling precision | |
| CN117195488A (en) | Construction method and related system of unrefined regional land-sea integrated digital elevation model | |
| CN114187399A (en) | Method and system for generating illumination simulation image on surface of extraterrestrial star | |
| CN110310370A (en) | A method of realizing that GPS is merged with SRTM point face based on integral adjustment | |
| CN108132140A (en) | A kind of star-loaded optical remote sensing device full filed column radiance determines method | |
| CN115962755A (en) | Earth and stone calculation method based on unmanned aerial vehicle oblique photography technology | |
| Maliqi et al. | Creating and analysing the Digital Terrain Model of the Slivovo area using QGIS software | |
| CN102521880B (en) | Mars-surface three-dimensional drawing method | |
| KR20110072257A (en) | Method for analyzing sunshine, analyzing system of sunshine, and computer readable media | |
| Bhattacharya et al. | A Novel 2.5 D Shadow Calculation Algorithm for Urban Environment. | |
| JP7179368B2 (en) | Haunting time providing method, haunting time providing device and haunting time providing program | |
| Vassilopoulou | A USEFUL SOFTWARE FOR TERRAIN ANALYSIS AND GEOENVIRONMENTAL APPLICATIONS–STUDY CASE ON THE GEODYNAMIC EVOLUTION OF ARGOLIS PENINSULA, GREECE. | |
| Noinak et al. | Testing Horizontal Coordinate Correction Model Used for Transformation from PPP GNSS Technique to Thai GNSS CORS Network Based on ITRF2014 | |
| Dinkov | Generation of 3D panoramic map for tourism management applications | |
| Idi | Spatio-temporal mapping of solar energy potential of Dutse, Jigawa State Nigeria using arcgis solar radiation spatial analyst tool |
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 |