Celestial body surface illumination simulation method and device

Abstract

The invention provides a celestial body surface illumination simulation method and device, and the method comprises the steps: converting initial topographic data under a first resolution ratio and initial topographic data under a second resolution ratio into a target coordinate system, and obtaining first resolution ratio topographic data and second resolution ratio topographic data, determining a target area in a target coordinate system based on the first resolution terrain data; according to a second target data point corresponding to the first target data point, constructing a second resolution terrain maximum elevation angle-occlusion point database; according to a preset search window and the second resolution terrain maximum elevation angle-shielding point database, constructing a first resolution terrain maximum elevation angle-shielding point database; and simulating surface illumination of the celestial body based on a planet calendar related to the celestial body, the local terrain data of each first target data point and the first resolution terrain maximum elevation angle-occlusion point database to obtain a celestial body illumination simulation image.

Description

Technical Field

[0001] The present disclosure relates to the field of image simulation technology, and more specifically, to a method and device for simulating illumination of a celestial body surface. Background Art

[0002] The unique lighting conditions and surface environmental characteristics of the lunar polar regions are ideal places for water enrichment and preservation. Therefore, accurately acquiring the lighting conditions of the lunar polar regions with high temporal and spatial resolution plays an important role in landing area selection, solar energy resource utilization, and exploration path planning for future polar exploration missions.

[0003] In the process of realizing the concept disclosed in the present invention, the inventors discovered that there are at least the following problems in the related technology: Since the terrain conditions in the lunar polar regions are relatively complex, the influence of terrain factors on the sunlight conditions is very significant, thus forming a complex lighting environment, which results in a long calculation time and low calculation accuracy for the lighting conditions in the lunar polar regions. Summary of the invention

[0004] In view of this, the present disclosure provides a method and device for simulating illumination of a celestial body surface.

[0005] One aspect of the present disclosure provides a method for simulating illumination of a celestial body surface, comprising: converting initial terrain data at a first resolution and initial terrain data at a second resolution into a target coordinate system, respectively obtaining first-resolution terrain data and second-resolution terrain data, and determining a target area in the target coordinate system based on the first-resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, where m ≥ 1 and is a positive integer; for each first target data point, based on n second target data points corresponding to the first target data point, constructing a second-resolution terrain maximum altitude angle-occlusion point database related to the first target data point, wherein the second-resolution terrain maximum altitude angle-occlusion point database includes the first target data point at any observation The second-resolution terrain maximum altitude angle and the second-resolution occlusion point on the measured azimuth, n ≥ 1 and is a positive integer; for each first target data point, according to the preset search window and the second-resolution terrain maximum altitude angle-occlusion point database, a first-resolution terrain maximum altitude angle-occlusion point database related to the first target data point is constructed, wherein the first-resolution terrain maximum altitude angle-occlusion point database includes the first-resolution terrain maximum altitude angle and the first-resolution occlusion point of the first target data point at any observation azimuth; and based on the planetary ephemeris related to the celestial body, the local terrain data of each first target data point and the first-resolution terrain maximum altitude angle-occlusion point database of each first target data point, the surface illumination of the celestial body is simulated to obtain a celestial body illumination simulation image.

[0006] Another aspect of the present disclosure provides a celestial body surface illumination simulation device, including: a determination module, used to convert initial terrain data at a first resolution and initial terrain data at a second resolution into a target coordinate system, respectively obtain first resolution terrain data and second resolution terrain data, and determine a target area in the target coordinate system based on the first resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, m ≥ 1 and is a positive integer; a first construction module, used to construct, for each first target data point, a first target data point corresponding to the first target data point according to n second target data points corresponding to the first target data point A second resolution terrain maximum altitude angle-occlusion point database related to the target data point, n≥1 and is a positive integer; a second construction module, used to construct a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point for each first target data point according to a preset search window and the second resolution terrain maximum altitude angle-occlusion point database; and a simulation module, used to simulate the illumination of the celestial body surface based on the planetary ephemeris related to the celestial body, the local terrain data of each first target data point and the first resolution terrain maximum altitude angle-occlusion point database of each first target data point, so as to obtain a celestial body illumination simulation image.

[0007] According to the embodiments of the present disclosure, by uniformly converting the initial terrain data at different resolutions into the same target coordinate system, the first resolution terrain data and the second resolution terrain data can be obtained. The second resolution terrain data with a lower resolution can be used to preliminarily construct a second resolution terrain maximum altitude angle-occlusion point database, and then all the first target data points in the same direction as each second resolution occlusion point can be quickly screened based on the preset search window, so as to further accurately calculate the first resolution terrain maximum altitude angle-occlusion point database, thereby realizing the rapid construction of a higher resolution terrain maximum altitude angle-occlusion point database. Therefore, the celestial body surface illumination simulation method provided by the embodiments of the present disclosure can realize the rapid calculation of the celestial body surface illumination conditions at high resolution under limited computing resources and tolerable time conditions, and fully consider the influence of the slope and slope direction of the local terrain on the illumination conditions, thereby improving the accuracy of the calculation, and thus can at least partially overcome the problems of long illumination condition calculation time and low calculation accuracy existing in the related art. BRIEF DESCRIPTION OF THE DRAWINGS

[0008] The above and other objects, features and advantages of the present disclosure will become more apparent through the following description of the embodiments of the present disclosure with reference to the accompanying drawings, in which:

[0009] Figure 1 A diagram schematically shows an application scenario of the method and device for simulating illumination of a celestial body surface according to an embodiment of the present disclosure;

[0010] Figure 2 A flowchart of a method for simulating illumination of a celestial body surface according to an embodiment of the present disclosure is schematically shown;

[0011] Figure 3 A schematic diagram of constructing a second-resolution terrain maximum altitude angle-occlusion point database according to an embodiment of the present disclosure is schematically shown;

[0012] Figure 4 A schematic diagram of constructing a first-resolution terrain maximum altitude angle-occlusion point database according to an embodiment of the present disclosure is schematically shown;

[0013] Figure 5 An astronomical illumination simulation image according to an embodiment of the present disclosure is exemplarily shown;

[0014] Figure 6 An astronomical illumination simulation image according to another embodiment of the present disclosure is exemplarily shown;

[0015] Figure 7 A schematic diagram of determining the solar visibility coefficient according to an embodiment of the present disclosure is schematically shown;

[0016] Figure 8 A block diagram of a celestial body surface illumination simulation device according to an embodiment of the present disclosure is schematically shown. DETAILED DESCRIPTION

[0017] Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of the present disclosure. In the following detailed description, for ease of explanation, many specific details are set forth to provide a comprehensive understanding of the embodiments of the present disclosure. However, it is apparent that one or more embodiments may also be implemented without these specific details. In addition, in the following description, descriptions of known structures and technologies are omitted to avoid unnecessary confusion of the concepts of the present disclosure.

[0018] The terms used herein are only for describing specific embodiments and are not intended to limit the present disclosure. The terms "comprise", "include", etc. used herein indicate the existence of the features, steps, operations and / or components, but do not exclude the existence or addition of one or more other features, steps, operations or components.

[0019] All terms (including technical and scientific terms) used herein have the meanings commonly understood by those skilled in the art unless otherwise defined. It should be noted that the terms used herein should be interpreted as having a meaning consistent with the context of this specification and should not be interpreted in an idealized or overly rigid manner.

[0020] When using expressions such as "at least one of A, B, and C, etc.", they should generally be interpreted according to the meaning of the expression commonly understood by those skilled in the art (for example, "a system having at least one of A, B, and C" should include but is not limited to a system having A alone, B alone, C alone, A and B, A and C, B and C, and / or A, B, C, etc.).

[0021] In the embodiments of the present disclosure, the collection, updating, analysis, processing, use, transmission, provision, disclosure, storage, etc. of the data involved (for example, including but not limited to user personal information) are in compliance with the provisions of relevant laws and regulations, are used for legitimate purposes, and do not violate public order and good morals. In particular, necessary measures are taken for user personal information to prevent illegal access to user personal information data and maintain the security of user personal information and network security.

[0022] In the embodiments of the present disclosure, the user's authorization or consent is obtained before obtaining or collecting the user's personal information.

[0023] In an embodiment of the present disclosure, the celestial body may include the moon.

[0024] The unique lighting conditions and surface environmental characteristics of the lunar polar regions are ideal places for water enrichment and preservation. Therefore, accurately acquiring the lighting conditions of the lunar polar regions with high temporal and spatial resolution plays an important role in landing area selection, solar energy resource utilization, and exploration path planning for future polar exploration missions.

[0025] In the process of realizing the concept disclosed in the present invention, the inventors discovered that there are at least the following problems in the related technology: Since the terrain conditions in the lunar polar regions are relatively complex, the influence of terrain factors on the sunlight conditions is very significant, thus forming a complex lighting environment, which results in a long calculation time and low calculation accuracy for the lighting conditions in the lunar polar regions.

[0026] According to the embodiments of the present disclosure, since the angle between the lunar equatorial plane and the ecliptic plane is about 1.54°, the sunlight radiates the lunar poles at a nearly horizontal angle. Polar impact craters are shielded by their own crater edges and the surrounding higher terrain to form seasonal or permanent shadow areas. Due to the lack of solar radiation, the temperature can be as low as ~18 K. Water brought by the impact of comets and water-containing meteorites, magma degassing, solar wind, etc. will be captured by permanent shadow craters in the form of "cold traps" and preserved for a long time. Through the detection of water ice in the permanent shadow areas of the polar regions, scientific support is provided for studying the origin and evolution of lunar water, the characteristics of volatiles in the early solar system, etc.; engineering provides an important energy foundation for the exploitation and utilization of water resources and the establishment of lunar bases. Due to the dense impact craters and large terrain undulations in the lunar polar regions, as the moon rotates, the nearly horizontal incident sunlight is blocked by the terrain, forming a complex lighting environment.

[0027] According to an embodiment of the present disclosure, the relevant method for determining the lighting conditions of the lunar polar regions may include a horizon (or maximum altitude angle) method, which calculates the lighting conditions by calculating the maximum altitude angle of all polar locations in advance once and then comparing it with the solar altitude angle information of any time period, thereby saving a lot of calculation time. However, since the calculation time increases exponentially with the increase in resolution, there is still a time-consuming problem when carrying out high-resolution polar lighting condition calculations. In addition, the relevant method for determining the lighting conditions of the lunar polar regions usually does not take into account the slope and slope aspect of the local terrain of the lunar polar regions, which directly affects the real-time lighting condition simulation of the local terrain, the effective solar irradiance, and the accurate calculation of solar energy.

[0028] In order to at least partially solve the technical problems existing in the related art, the present disclosure provides a method and device for celestial body surface illumination simulation. The celestial body surface illumination simulation method includes: converting initial terrain data at a first resolution and initial terrain data at a second resolution into a target coordinate system, respectively obtaining first resolution terrain data and second resolution terrain data, and determining a target area in the target coordinate system based on the first resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, m ≥ 1 and is a positive integer; for each first target data point, based on the n second target data points corresponding to the first target data point, constructing a second resolution terrain maximum altitude angle-occlusion point database related to the first target data point, wherein the second resolution terrain maximum altitude angle-occlusion point database includes the first target data point at any observation azimuth. two-resolution terrain maximum altitude angle and second-resolution occlusion point, n≥1 and is a positive integer; for each first target data point, according to the preset search window and the second resolution terrain maximum altitude angle-occlusion point database, a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point is constructed, wherein the first resolution terrain maximum altitude angle-occlusion point database includes the first resolution terrain maximum altitude angle and the first resolution occlusion point of the first target data point at any observation azimuth; and based on the planetary ephemeris related to the celestial body, the local terrain data of each first target data point and the first resolution terrain maximum altitude angle-occlusion point database of each first target data point, the illumination of the celestial body surface is simulated to obtain a celestial body illumination simulation image.

[0029] Figure 1 The application scenario diagram of the celestial body surface simulation method and device according to the embodiment of the present disclosure is schematically shown.

[0030] like Figure 1As shown, the application scenario 100 according to this embodiment may include a first terminal device 101, a second terminal device 102, a third terminal device 103, a network 104, and a server 105. The network 104 is used to provide a medium for a communication link between the first terminal device 101, the second terminal device 102, the third terminal device 103, and the server 105. The network 104 may include various connection types, such as wired, wireless communication links, or optical fiber cables, etc.

[0031] The user can use the first terminal device 101, the second terminal device 102, and the third terminal device 103 to interact with the server 105 through the network 104 to receive or send messages, etc. Various communication client applications can be installed on the first terminal device 101, the second terminal device 102, and the third terminal device 103, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients, social platform software, etc. (only for example).

[0032] The first terminal device 101, the second terminal device 102, and the third terminal device 103 may be various electronic devices having display screens and supporting web browsing, including but not limited to smart phones, tablet computers, laptop computers, desktop computers, and the like.

[0033] The server 105 may be a server that provides various services, such as a background management server (only as an example) that provides support for websites browsed by users using the first terminal device 101, the second terminal device 102, and the third terminal device 103. The background management server may analyze and process the received data such as user requests, and feed back the processing results (such as web pages, information, or data obtained or generated according to user requests) to the terminal device.

[0034] It should be noted that the celestial body surface illumination simulation method provided in the embodiment of the present disclosure can generally be executed by the server 105. Accordingly, the celestial body surface illumination simulation device provided in the embodiment of the present disclosure can generally be set in the server 105. The celestial body surface illumination simulation method provided in the embodiment of the present disclosure can also be executed by a server or server cluster that is different from the server 105 and can communicate with the first terminal device 101, the second terminal device 102, the third terminal device 103 and / or the server 105. Correspondingly, the celestial body surface illumination simulation device provided in the embodiment of the present disclosure can also be set in a server or server cluster that is different from the server 105 and can communicate with the first terminal device 101, the second terminal device 102, the third terminal device 103 and / or the server 105. Alternatively, the celestial body surface illumination simulation method provided in the embodiment of the present disclosure can also be executed by the first terminal device 101, the second terminal device 102 or the third terminal device 103, or by other terminal devices that are different from the first terminal device 101, the second terminal device 102 or the third terminal device 103. Correspondingly, the celestial body surface illumination simulation device provided in the embodiment of the present disclosure may also be arranged in the first terminal device 101, the second terminal device 102 or the third terminal device 103, or in other terminal devices different from the first terminal device 101, the second terminal device 102 or the third terminal device 103.

[0035] For example, the initial terrain data at the first resolution and the initial terrain data at the second resolution may be originally stored in any one of the first terminal device 101, the second terminal device 102 or the third terminal device 103 (for example, the first terminal device 101, but not limited thereto), or stored in an external storage device and may be imported into the first terminal device 101. Then, the first terminal device 101 may locally execute the celestial body surface illumination simulation method provided by the embodiment of the present disclosure, or send the initial terrain data at the first resolution and the initial terrain data at the second resolution to other terminal devices, servers, or server clusters, and the other terminal devices, servers, or server clusters that receive the initial terrain data at the first resolution and the initial terrain data at the second resolution execute the celestial body surface illumination simulation method provided by the embodiment of the present disclosure.

[0036] It should be understood that Figure 1 The number of terminal devices, networks and servers in the embodiment is only for illustration. Any number of terminal devices, networks and servers may be provided according to the implementation requirements.

[0037] Figure 2 The flowchart of the method for simulating illumination of a celestial body surface according to an embodiment of the present disclosure is schematically shown.

[0038] like Figure 2As shown, the method 200 includes operations S210 to S240.

[0039] In operation S210, initial terrain data at a first resolution and initial terrain data at a second resolution are converted into a target coordinate system to obtain first resolution terrain data and second resolution terrain data, respectively, and a target area is determined in the target coordinate system based on the first resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, where m≥1 and is a positive integer.

[0040] In operation S220, for each first target data point, a second resolution terrain maximum altitude angle-occlusion point database related to the first target data point is constructed based on the n second target data points corresponding to the first target data point, wherein the second resolution terrain maximum altitude angle-occlusion point database includes the second resolution terrain maximum altitude angle and second resolution occlusion point of the first target data point at any observation azimuth, and n≥1 is a positive integer.

[0041] In operation S230, for each first target data point, a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point is constructed according to a preset search window and a second resolution terrain maximum altitude angle-occlusion point database, wherein the first resolution terrain maximum altitude angle-occlusion point database includes the first resolution terrain maximum altitude angle and the first resolution occlusion point of the first target data point at any observation azimuth.

[0042] In operation S240, based on the planetary ephemeris associated with the celestial body, the local terrain data of each first target data point and the first resolution terrain maximum altitude angle-occlusion point database of each first target data point, the illumination of the celestial body surface is simulated to obtain a celestial body illumination simulation image.

[0043] According to an embodiment of the present disclosure, the initial terrain data may include a digital elevation model (DEM). DEM can be understood as a physical ground model that represents ground elevation in the form of a numerical array, and can be used to reflect local terrain features at a certain resolution. The resolution of DEM defines the length of the smallest cell (or grid point) in the DEM, that is, the interval between discrete data points of the surface elevation. In fact, DEM can be regarded as a grid composed of many small squares, each of which contains the elevation information of a place. The side length of these small squares is the resolution of DEM. The smaller the resolution value, the smaller the geographic space represented by each small square, that is, they more finely depict the elevation changes of the surface.

[0044] According to an embodiment of the present disclosure, the initial terrain data may include a number of grid points (which may be understood as minimum cells), each grid point having plane coordinates and elevation information corresponding to the plane coordinates, and the resolution of the initial terrain data defines the length of the grid points.

[0045] In the embodiments of the present disclosure, initial topographic data at a first resolution (hereinafter referred to as the first initial topographic data) and initial topographic data at a second resolution (hereinafter referred to as the second initial topographic data) both correspond to the same area on the surface of a celestial body (e.g., both correspond to the same local topographic area at the lunar polar region). Since the first resolution is higher than the second resolution, it can be understood that the number of data points included in the first initial topographic data is greater than the number of data points included in the second initial topographic data, and the accuracy of the first initial topographic data is higher than the accuracy of the second initial topographic data.

[0046] According to the embodiments of the present disclosure, those skilled in the art may select initial terrain data of appropriate resolution according to actual needs or application scenarios, etc., which are not specifically limited here. As an example, the first initial terrain data may be, for example, a DEM with a spatial resolution of 20 meters, and the second initial terrain data may be, for example, a DEM with a spatial resolution of 240 meters. As another example, the first initial terrain data may be, for example, a DEM with a spatial resolution of 5 meters, and the second initial terrain data may be, for example, a DEM with a spatial resolution of 20 meters.

[0047] According to an embodiment of the present disclosure, the first initial terrain data and the second initial terrain data can be read and uniformly converted to a target coordinate system (e.g., a two-dimensional Cartesian coordinate system XY) represented by rows and columns with a celestial pole (e.g., a lunar pole) as the center, to obtain first-resolution terrain data corresponding to the first initial terrain data and second-resolution terrain data corresponding to the second initial terrain data. Exemplarily, the first-resolution terrain data may include M first terrain data points, and the second-resolution terrain data may include N second terrain data points, where M and N ≥ 1 and are positive integers. Optionally, the M first terrain data points and the N second terrain data points may be in the same two-dimensional Cartesian coordinate system XY.

[0048] According to an embodiment of the present disclosure, a target area can be determined in a target coordinate system based on the first resolution terrain data, and the target area can be understood as an area of ​​interest. There are m first terrain data points falling into the target area, and the m first terrain data points can be determined as m first target data points corresponding to the target area, where m≥1 and is a positive integer.

[0049] According to an embodiment of the present disclosure, at an observation point A at a certain moment, the terrain maximum altitude angle method can determine the lighting conditions by comparing the relationship between the solar altitude angle of point A and the maximum altitude angle of the terrain in the same direction as the sun's radiation. For example, assuming that at a certain time in a certain place, approximately 32% of the solar disks and observation point A form a solar altitude angle greater than the maximum altitude angle of the terrain, it can be considered that the solar visibility coefficient of the place is 0.32. Further, the effective solar irradiance of the place can be calculated based on the solar visibility coefficient. Optionally, the maximum altitude angle of the terrain in each direction at each first target data point can be calculated and stored in the respective database of each first target data point, and the database can be called at different time periods later, thereby effectively reducing the calculation time.

[0050] According to an embodiment of the present disclosure, the maximum terrain elevation angle in each direction at the first target data point and the occlusion point corresponding to the maximum terrain elevation angle may be calculated and stored in a terrain maximum elevation angle-occlusion point database corresponding to the first target data point.

[0051] Exemplarily, for each first target data point, the second resolution terrain maximum altitude angle-occlusion point database includes the second resolution terrain maximum altitude angle and the second resolution occlusion point of the first target data point at any observation azimuth, and the first resolution terrain maximum altitude angle-occlusion point database includes the first resolution terrain maximum altitude angle and the first resolution occlusion point of the first target data point at any observation azimuth. Among them, the second resolution terrain maximum altitude angle and the second resolution occlusion point are calculated using the second resolution terrain data, and the first resolution terrain maximum altitude angle and the first resolution occlusion point are calculated using the first resolution terrain data. It can be understood that the accuracy of the first resolution terrain maximum altitude angle-occlusion point database is higher than the accuracy of the second resolution terrain maximum altitude angle-occlusion point database.

[0052] According to an embodiment of the present disclosure, each first target data point corresponds to n second target data points in the second resolution terrain data, and the first target data point is used as an observation point, and the n second target data points correspond to n groups of terrain altitude angles and occlusion points of the first target data point at the second resolution. As an example, for the first target data point B, the second terrain data points falling within the range of 0-360° in the target area can be determined as the second target data points with the point B as the center, and n second target data points are obtained.

[0053] According to an embodiment of the present disclosure, for each first target data point, n second target data points can be used to preliminarily calculate the maximum elevation angle and occlusion point of the terrain in each direction of the first target data point at a lower resolution, and obtain a second resolution terrain maximum elevation angle-occlusion point database. Based on a preset search window and the second resolution terrain maximum elevation angle-occlusion point database, all first target data points in the same direction as each occlusion point can be screened in the first resolution terrain data as candidate occlusion points, so as to further accurately calculate the maximum elevation angle and occlusion point of the terrain in each direction of the first target data point at a higher resolution, thereby realizing the rapid construction of the first resolution terrain maximum elevation angle-occlusion point database.

[0054] According to an embodiment of the present disclosure, the illumination condition of each first target data point can be determined based on the planetary ephemeris related to the celestial body, the local terrain data of each first target data point, and the first resolution terrain maximum altitude angle-occlusion point database, so that the illumination of the celestial body surface can be simulated to obtain a celestial body illumination simulation image. The local terrain data of each first target data point can include the terrain slope angle and slope azimuth at each first target data point.

[0055] According to the embodiments of the present disclosure, by uniformly converting the initial terrain data at different resolutions into the same target coordinate system, the first resolution terrain data and the second resolution terrain data can be obtained. The second resolution terrain data with a lower resolution can be used to preliminarily construct a second resolution terrain maximum altitude angle-occlusion point database, and then all the first target data points in the same direction as each second resolution occlusion point can be quickly screened based on the preset search window, so as to further accurately calculate the first resolution terrain maximum altitude angle-occlusion point database, thereby realizing the rapid construction of a higher resolution terrain maximum altitude angle-occlusion point database. Therefore, the celestial body surface illumination simulation method provided by the embodiments of the present disclosure can realize the rapid calculation of the celestial body surface illumination conditions at high resolution under limited computing resources and tolerable time conditions, and fully consider the influence of the slope and slope direction of the local terrain on the illumination conditions, thereby improving the accuracy of the calculation, and thus can at least partially overcome the problems of long illumination condition calculation time and low calculation accuracy existing in the related art.

[0056] According to an embodiment of the present disclosure, initial terrain data at a first resolution and initial terrain data at a second resolution are converted into a target coordinate system to obtain first-resolution terrain data and second-resolution terrain data, respectively, including: constructing a target coordinate system with a celestial pole as the origin; converting the initial terrain data at the first resolution into the target coordinate system to obtain first-resolution terrain data, the first-resolution terrain data including M first terrain data points, each of which includes first longitude information, first latitude information and first elevation information, M≥1 and is a positive integer; and converting the initial terrain data at the second resolution into the target coordinate system to obtain second-resolution terrain data, the second-resolution terrain data including N second terrain data points, each of which includes second longitude information, second latitude information and second elevation information, N≥1 and is a positive integer; determining a target area in the target coordinate system based on the first-resolution terrain data includes: determining a target area within a candidate area defined by the M first terrain data points, wherein there are m first target data points corresponding to the target area.

[0057] According to one embodiment of the present disclosure, a Cartesian coordinate system can be constructed with the lunar pole as the origin, the section parallel to the equatorial plane is the Xoy plain, and the axial direction is determined according to the right-hand rule. The first initial terrain data and the second initial terrain data are both DEM data. The first initial terrain data can be read, and the DEM raster data at the first resolution is converted into the first text data in a text format. The second initial terrain data can be read, and the DEM raster data at the second resolution is converted into the second text data in a text format. The first text data and the second text data can be projected onto a unified Cartesian coordinate system based on the polar cube projection method. Optionally, the first text data and the second text data can be projected into the same Cartesian coordinate system to facilitate the subsequent calculation of the terrain altitude angle at the second resolution and the screening of candidate occlusion points at the first resolution.

[0058] According to an embodiment of the present disclosure, each data point in the target coordinate system can be regarded as a grid point, and the longitude information, latitude information and elevation information corresponding to each data point can be determined according to the conversion relationship between the projection information and the geographic coordinates. For example, in the target coordinate system, the first resolution terrain data includes M first terrain data points, the second resolution terrain data includes N second terrain data points, each first terrain data point includes first longitude information, first latitude information and first elevation information, each second terrain data point includes second longitude information, second latitude information and second elevation information, N≥M≥1, and M and N are positive integers.

[0059] As an example, the (x, y) coordinates of each grid point in the target coordinate system can be obtained. Based on the polar stereographic projection method, the longitude of any two grid points can be calculated ( ) and latitude ( ) information, and establish the conversion relationship between the projection information of each raster and the geographic coordinates (longitude, latitude).

[0060] For example:

[0061] (1)

[0062] (2)

[0063] in,( , ) represents the longitude and latitude coordinates of the projection center point; (x, y) is the position coordinates of the grid point in the Cartesian coordinate system; parameters and Respectively expressed as:

[0064] (3)

[0065] (4)

[0066] in, is the average radius of a celestial body (such as the Moon).

[0067] According to an embodiment of the present disclosure, a target area may be determined within a candidate area defined by M first terrain data points to ensure data integrity of the target area, thereby ensuring validity of subsequent calculations.

[0068] According to an embodiment of the present disclosure, for each first target data point, based on the n second target data points corresponding to the first target data point, constructing a second resolution terrain maximum altitude angle-occlusion point database related to the first target data point includes: taking the first target data point as the observation point, searching for all second terrain data points falling within a preset azimuth range among N second terrain data points based on a preset azimuth step size, obtaining n second target data points, and the preset azimuth range is [0, 2π]; determining the second resolution candidate terrain altitude angles respectively between the first target data point and the n second target data points based on the first longitude information, the first latitude information, the first elevation information, the second longitude information, the second latitude information, and the second elevation information; for any selected observation azimuth, determining the second resolution candidate terrain altitude angle that is the largest at the selected observation azimuth as the second resolution terrain maximum altitude angle corresponding to the selected observation azimuth; and constructing a second resolution terrain maximum altitude angle-occlusion point database based on all second resolution terrain maximum altitude angles of the first target data point within the preset azimuth range, and the second resolution occlusion points corresponding to the second resolution terrain maximum altitude angle.

[0069] Figure 3 A schematic diagram of constructing a second-resolution terrain maximum altitude angle-occlusion point database according to an embodiment of the present disclosure is schematically shown.

[0070] According to one embodiment of the present disclosure, Figure 3 As shown, the first target data point can be used as the observation point, and with the first target data point as the center, all second terrain data points at each observation azimuth within the range of 0-360° can be searched one by one in a step of 0.5° in the clockwise direction of the Y-axis to obtain n second target data points corresponding to the aforementioned first target data point.

[0071] According to an embodiment of the present disclosure, for each first target data point, n second target data points can be used as candidate occlusion points, and the terrain altitude angle between each second target data point and the first target data point is calculated based on the first longitude information, the first latitude information, the first elevation information, the second longitude information, the second latitude information and the second elevation information, respectively, to obtain n groups of second resolution candidate terrain altitude angles and second resolution candidate occlusion points.

[0072] According to the embodiments of the present disclosure, the maximum second-resolution candidate terrain altitude angle in each observation direction can be screened out, the second-resolution candidate terrain altitude angle is determined as the second-resolution terrain maximum altitude angle, and the corresponding second-resolution candidate occlusion point is determined as the second-resolution occlusion point. According to the second-resolution terrain maximum altitude angle and the second-resolution occlusion point corresponding to each first target data point, a second-resolution terrain maximum altitude angle-occlusion point database for each first target data point can be constructed.

[0073] According to an embodiment of the present disclosure, for each second target data point, based on the first longitude information, the first latitude information, the first elevation information, the second longitude information, the second latitude information and the second elevation information, determining the second resolution candidate terrain altitude angle between the first target data point and the n second target data points includes: determining the target distance between the first target data point and the second target data point based on the first longitude information, the first latitude information, the second longitude information, the second latitude information and the celestial body radius; correcting the first elevation information and the second elevation information according to the target distance and the celestial body radius, respectively, to obtain the first elevation correction information and the second elevation correction information correspondingly; and determining the second resolution candidate terrain altitude angle between the first target data point and the second target data point according to the target distance, the celestial body radius, the first elevation correction information and the second elevation correction information.

[0074] According to an embodiment of the present disclosure, the target distance (d) between any two grid points can be calculated based on the longitude and latitude information of any grid point in the acquired target coordinate system. Exemplarily, the target distance between the first target data point and the second target data point can be determined based on the first longitude information, the first latitude information, the second longitude information, the second latitude information and the radius of the celestial body.

[0075] For example:

[0076] (5)

[0077] in, is the half-sine function of the angle between any two grid points and can be expressed as:

[0078] (6)

[0079] in, and Represents the longitude and latitude coordinates of any two grid points respectively.

[0080] According to an embodiment of the present disclosure, the elevation information of any grid point can be corrected according to the target distance and the radius of the celestial body, and the corrected elevation of any grid point in the target area can be calculated ( ). Exemplarily, the first elevation information and the second elevation information may be corrected according to the target distance and the radius of the celestial body, respectively, to obtain the first elevation correction information and the second elevation correction information accordingly.

[0081] For example:

[0082] (7)

[0083] in, is the elevation before correction.

[0084] According to one embodiment of the present disclosure, the observation point ( ) and the occlusion point ( ) and the horizontal plane ( ), that is, the terrain elevation angle. Exemplarily, the second resolution candidate terrain elevation angle between the first target data point (observation point) and the second target data point (occlusion point) can be calculated according to the target distance, the radius of the celestial body, the first elevation correction information, and the second elevation correction information.

[0085] For example:

[0086] (8)

[0087] According to an embodiment of the present disclosure, for each first target data point, according to a preset search window and a second resolution terrain maximum altitude angle-occlusion point database, a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point is constructed, including: taking the first target data point as an observation point, for any selected observation azimuth, constructing a preset search window according to the second resolution occlusion point corresponding to the selected azimuth; determining all first target data points located in the search window at the selected observation azimuth as first resolution candidate occlusion points, and obtaining p first resolution candidate occlusion points corresponding to the selected observation azimuth; ... A target data point and p first resolution candidate occlusion points corresponding to a selected observation azimuth, determining the first resolution candidate terrain altitude angles respectively between the first target data point and the p first resolution candidate occlusion points; for any selected observation azimuth, determining the largest first resolution candidate terrain altitude angle located on the selected observation azimuth as the first resolution terrain maximum altitude angle corresponding to the selected observation azimuth; and constructing a first resolution terrain maximum altitude angle-occlusion point database based on all first resolution terrain maximum altitude angles of the first target data point within a preset azimuth range and the first resolution occlusion points corresponding to the first resolution terrain maximum altitude angle.

[0088] Figure 4 A schematic diagram of constructing a first-resolution terrain maximum altitude angle-occlusion point database according to an embodiment of the present disclosure is schematically shown.

[0089] According to one embodiment of the present disclosure, Figure 4 As shown, the first target data point can be used as an observation point, and a preset search window can be constructed according to the second resolution occlusion point at the selected observation azimuth. All the first target data points in the search window at the selected observation azimuth can be determined as first resolution candidate occlusion points, and p first resolution candidate occlusion points corresponding to the selected observation azimuth can be obtained. In this way, p first resolution candidate occlusion points of each first target data point in each observation direction can be obtained.

[0090] According to an embodiment of the present disclosure, for each observation azimuth, the terrain elevation angles between the first target data point and the p first-resolution candidate occlusion points can be calculated to obtain the p first-resolution candidate terrain elevation angles for each observation azimuth. Thus, the p first-resolution candidate terrain elevation angles for each first target data point in each observation direction can be obtained. It should be noted that the calculation method of the first-resolution candidate terrain elevation angle can refer to the calculation method of the second-resolution candidate terrain elevation angle (Formula 5-Formula 8) mentioned above, which will not be repeated here.

[0091] According to the embodiments of the present disclosure, the maximum first-resolution candidate terrain altitude angle in each observation direction can be screened out, the first-resolution candidate terrain altitude angle is determined as the first-resolution terrain maximum altitude angle, and the corresponding first-resolution candidate occlusion point is determined as the first-resolution occlusion point. According to the first-resolution terrain maximum altitude angle and the first-resolution occlusion point corresponding to each first target data point, a first-resolution terrain maximum altitude angle-occlusion point database for each first target data point can be constructed.

[0092] According to an embodiment of the present disclosure, by first using the second-resolution terrain data with a lower resolution to preliminarily construct a second-resolution terrain maximum altitude angle-occlusion point database, and then quickly screening p first-resolution candidate occlusion points in each observation direction based on a preset search window, the first-resolution terrain maximum altitude angle is further accurately calculated, thereby realizing the rapid construction of a higher-resolution terrain maximum altitude angle-occlusion point database. Therefore, the celestial body surface illumination simulation method provided by the embodiment of the present disclosure can realize the rapid calculation of celestial body surface illumination conditions at high resolution under limited computing resources and tolerable time conditions. In addition, since the observation points and occlusion points involved in calculating the first-resolution terrain maximum altitude angle are all data points at the first resolution, it can be ensured that the calculation accuracy of each first-resolution terrain maximum altitude angle is the first resolution, so that the calculation accuracy of the illumination conditions at various locations in the target area can be made consistent, thereby helping to improve the calculation accuracy and calculation accuracy of the celestial body illumination simulation image.

[0093] According to an embodiment of the present disclosure, the preset search window is centered on the second resolution occlusion point and is constructed based on a preset window side length, and the preset window side length is determined based on the size of the second resolution.

[0094] According to the embodiments of the present disclosure, those skilled in the art may determine a reasonable preset window side length based on actual needs or application scenarios, etc., and no specific limitation is made here. Figure 4 As shown, a preset search window can be constructed with the second resolution occlusion point as the center and twice the second resolution as the preset window side length.

[0095] According to an embodiment of the present disclosure, local terrain data includes terrain slope angle and slope azimuth, and based on a planetary ephemeris related to the celestial body, local terrain data of each first target data point and a first resolution terrain maximum altitude angle-blocking point database of each first target data point, the illumination of the celestial body surface is simulated to obtain a celestial body illumination simulation image, including: determining the sun-celestial body distance and the longitude information and latitude information of the direct point of the sun at a selected time according to the planetary ephemeris related to the celestial body; determining the solar altitude angle and solar azimuth angle of each first target data point according to the first longitude information, the first latitude information, the longitude information of the direct point and the latitude information of the direct point; determining the solar altitude angle and solar azimuth angle of each first target data point according to each first target data point The solar visibility coefficient of each first target data point is determined based on the respective first resolution terrain maximum altitude angle-blocking point database, as well as the solar altitude angle and the solar azimuth angle, wherein the solar visibility coefficient represents the ratio of the visible area of ​​the solar disk observed at the first target data point to the area of ​​the solar disk; for each first target data point, the effective solar irradiance of each first target data point is determined based on the sun-celestial body distance, the solar visibility coefficient, the solar altitude angle, the solar azimuth angle, and the terrain slope angle and the slope azimuth angle at the first target observation point; and based on the respective effective solar irradiance of each first target data point, the illumination on the celestial body surface is simulated to obtain a celestial body illumination simulation image.

[0096] According to one embodiment of the present disclosure, the celestial body may be the moon. Exemplarily, the planetary ephemeris may be queried to obtain the distance between the sun and the moon and the coordinate information of the subsolar point (longitude information of the direct point and latitude information of the direct point) with a time resolution of 1 hour from 00:00 on January 1, 2026 to 24:00 on December 31, 2026. From this, the three component vectors of the x, y, and z axes pointing to the center of the solar disk of each first target data point in the target area may be calculated.

[0097] For example:

[0098] (9)

[0099] (10)

[0100] (11)

[0101] in, and are the latitude and longitude coordinates of the first target data point, respectively. and The coordinates of the point under the sun are the latitude and longitude coordinates. The three component vectors satisfy From this, the solar altitude angle can be calculated ( ) and the solar azimuth ( ),Right now

[0102] (12)

[0103] (13)

[0104] Figure 5 An astronomical illumination simulation image according to an embodiment of the present disclosure is exemplarily shown. Figure 6 An astronomical illumination simulation image according to another embodiment of the present disclosure is exemplarily shown.

[0105] According to an embodiment of the present disclosure, the solar visibility coefficient of each first target data point can be determined based on the first resolution terrain maximum altitude angle-blocking point database of each first target data point, as well as the solar altitude angle and the solar azimuth angle. For each first target data point, the effective solar irradiance of each first target data point can be determined based on the sun-celestial body distance, the solar visibility coefficient, the solar altitude angle, the solar azimuth angle, and the terrain slope angle and the slope azimuth angle at the first target observation point. The illumination of the celestial body surface can be simulated based on the effective solar irradiance of each first target data point to obtain a celestial body illumination simulation image.

[0106] As an example, Figure 5 The following is an example of a simulated astronomical illumination image of the target area at 03:00 on November 11, 2026. SH represents the Shackleton crater, and CR1 and SR1 represent areas with friendly illumination conditions. As another example, Figure 6 The following is an example of a simulated image of celestial body illumination at 16:00 on November 29, 2026. SH represents the Shackleton crater, and CR1 and SR1 represent areas with friendly illumination conditions.

[0107] According to an embodiment of the present disclosure, for each first target data point, according to the respective first resolution terrain maximum altitude angle-obstruction point database of each first target data point, as well as the solar altitude angle and the solar azimuth angle, determining the respective solar visibility coefficient of each first target data point includes: based on the solar azimuth angle, determining the first reference terrain maximum altitude angle and the second reference terrain maximum altitude angle in the first resolution terrain maximum altitude angle-obstruction point database, wherein the first reference terrain maximum altitude angle corresponds to the first observation azimuth angle, the second reference terrain maximum altitude angle corresponds to the second observation azimuth angle, and the first observation azimuth angle and the second observation azimuth angle are respectively located on both sides of the solar azimuth angle; dividing the visible area of ​​the solar disk based on the first adjacent terrain maximum altitude angle and the second adjacent terrain maximum altitude angle; and determining the solar visibility coefficient based on the ratio of the visible area area of ​​the solar disk to the total area of ​​the solar disk.

[0108] Figure 7 A schematic diagram of determining the solar visibility coefficient according to an embodiment of the present disclosure is schematically shown.

[0109] According to one embodiment of the present disclosure, Figure 7 As shown, the terrain maximum altitude angles a and b corresponding to the azimuth angles a' and b' can be calculated based on the first resolution terrain maximum altitude angle-blocking point database by using the interpolation method, wherein the azimuth angle a' is the first observation azimuth angle, the azimuth angle b' is the second observation azimuth angle, a is the first reference terrain maximum altitude angle, b is the second reference terrain maximum altitude angle, and the azimuth angles a' and b' are located on both sides of the solar azimuth angle. For example, it can be obtained by Figure 7 Point O in the image is moved 0.25° to the left and right respectively to obtain azimuth angles a' and b'. The straight line defined by a and b can be used to divide the solar disk into visible area and obscured area ( Figure 7 The gray area in the figure is the blocked area), thus determining the proportion of the visible area of ​​the solar disk, that is, the solar visibility coefficient Among them, when indicates that the solar disk is completely obscured; when indicates that the solar disk is fully visible; when Indicates that only part of the solar disk is visible.

[0110] For example, the solar disk area is , the visible area of ​​the solar disk ( ) can be approximately expressed as:

[0111] (14)

[0112] in, is the radius of the solar disk; is the angle between line segment oa and line segment ob (in radians), which is used to approximate the angle between the arc of the shadow section and the center of the circle. From this, the visible coefficient of the solar disk can be further calculated.

[0113] For example:

[0114] (15)

[0115] According to an embodiment of the present disclosure, determining the effective solar irradiance of each first target data point according to the sun-celestial body distance, the solar visibility coefficient, the solar altitude angle, the solar azimuth angle, and the terrain slope angle and the slope azimuth angle at the first target observation point includes: determining the solar incident angle at each first target data point according to the solar altitude angle, the solar azimuth angle, and the terrain slope angle and the slope azimuth angle at the first target observation point; and determining the effective solar irradiance of each first target data point according to the solar visibility coefficient, the solar incident angle at each first target data point and the sun-celestial body distance.

[0116] According to one embodiment of the present disclosure, the celestial body is the moon. The center of the visible area of ​​the solar disk can be obtained and the solar altitude angle ( ), combined with the terrain slope angle of each first target observation point in the target area ( ) and slope azimuth ( ), the solar incidence angle at each first target data point can be calculated.

[0117] For example:

[0118] (16)

[0119] For example, assuming that the lunar surface is approximately a Lambertian body, the effective solar irradiance of each first target observation point in the target area can be calculated according to the real-time distance between the sun and the moon ( ).

[0120] For example:

[0121] (17)

[0122] in, is the solar constant, is the distance between the sun and the moon, in AU; is the lunar albedo. On this basis, a celestial body illumination simulation image can be produced according to the effective solar irradiance of each first target data point.

[0123] According to the embodiments of the present disclosure, since the influence of the slope and aspect of the local terrain on the lighting conditions is fully considered, the real-time lighting condition simulation of the local terrain in the target area and the accurate calculation of the effective solar irradiance and solar energy can be improved.

[0124] Figure 8 The block diagram of the celestial body surface illumination simulation device 800 according to an embodiment of the present disclosure is schematically shown.

[0125] like Figure 8As shown, the celestial body surface illumination simulation device 800 includes a determination module 810 , a first construction module 820 , a second construction module 830 and a simulation module 840 .

[0126] A determination module 810 is used to convert the initial terrain data at the first resolution and the initial terrain data at the second resolution into a target coordinate system, obtain first resolution terrain data and second resolution terrain data respectively, and determine a target area in the target coordinate system based on the first resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, where m ≥ 1 and is a positive integer;

[0127] A first construction module 820 is used to construct, for each first target data point, a second resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to n second target data points corresponding to the first target data point, where n≥1 and is a positive integer;

[0128] A second construction module 830 is used to construct, for each first target data point, a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to a preset search window and a second resolution terrain maximum altitude angle-occlusion point database; and

[0129] The simulation module 840 is used to simulate the illumination of the celestial body surface based on the planetary ephemeris related to the celestial body, the local terrain data of each first target data point and the first resolution terrain maximum altitude angle-occlusion point database of each first target data point to obtain a celestial body illumination simulation image.

[0130] According to an embodiment of the present disclosure, the determination module 810 may include a coordinate system construction submodule, a first conversion submodule, a second conversion submodule, and a target area determination submodule.

[0131] The coordinate system construction submodule is used to construct the target coordinate system with the celestial pole as the origin.

[0132] The first conversion submodule is used to convert the initial terrain data at the first resolution to the target coordinate system to obtain the first resolution terrain data, the first resolution terrain data includes M first terrain data points, each of the first terrain data points includes first longitude information, first latitude information and first elevation information, M ≥ 1 and is a positive integer.

[0133] The second conversion submodule is used to convert the initial terrain data at the second resolution to the target coordinate system to obtain second resolution terrain data, the second resolution terrain data includes N second terrain data points, each second terrain data point includes second longitude information, second latitude information and second elevation information, N ≥ 1 and is a positive integer.

[0134] The target area determination submodule is used to determine the target area within the candidate area defined by the M first terrain data points, wherein the target area corresponds to m first target data points.

[0135] According to an embodiment of the present disclosure, the first construction module 820 may include a first search submodule, a first determination submodule, a second determination submodule, and a second resolution database construction submodule.

[0136] The first search submodule is used to take the first target data point as the observation point, search for all second terrain data points falling within a preset azimuth range among N second terrain data points based on a preset azimuth step size, and obtain n second target data points, and the preset azimuth range is [0, 2π].

[0137] The first determination submodule is used to determine the second resolution candidate terrain altitude angles between the first target data point and n second target data points based on the first longitude information, the first latitude information, the first elevation information, the second longitude information, the second latitude information and the second elevation information.

[0138] The second determination submodule is used to determine, for any selected observation azimuth, the largest second resolution candidate terrain altitude angle at the selected observation azimuth as the second resolution terrain maximum altitude angle corresponding to the selected observation azimuth.

[0139] The second resolution database construction submodule is used to construct a second resolution terrain maximum altitude angle-occlusion point database based on all second resolution terrain maximum altitude angles of the first target data point within a preset azimuth range, and the second resolution occlusion points corresponding to the second resolution terrain maximum altitude angle.

[0140] According to an embodiment of the present disclosure, the first determining submodule may include a target distance determining unit, a first correction unit, and a first candidate terrain altitude angle determining unit.

[0141] The target distance determination unit is used to determine the target distance between the first target data point and the second target data point based on the first longitude information, the first latitude information, the second longitude information, the second latitude information and the radius of the celestial body.

[0142] The first correction unit is used to correct the first elevation information and the second elevation information according to the target distance and the celestial body radius, respectively, to obtain the first elevation correction information and the second elevation correction information accordingly.

[0143] The first candidate terrain elevation angle determination unit is used to determine the candidate terrain elevation angle between the first target data point and the second target data point according to the target distance, the radius of the celestial body, the first elevation correction information and the second elevation correction information.

[0144] According to an embodiment of the present disclosure, the second construction module 830 may include a search window determination submodule, a second search submodule, a third determination submodule, a fourth determination submodule and a first resolution database construction submodule.

[0145] The search window determination submodule is used to take the first target data point as the observation point, and for any selected observation azimuth, to construct a preset search window according to the second resolution occlusion point corresponding to the selected azimuth.

[0146] A second search submodule is used to determine all first target data points located at the selected observation azimuth in the search window as first resolution candidate occlusion points, and obtain p first resolution candidate occlusion points corresponding to the selected observation azimuth;

[0147] The third determination submodule is used to determine, for any selected observation azimuth, the first resolution candidate terrain altitude angles respectively between the first target data point and the p first resolution candidate occlusion points corresponding to the selected observation azimuth.

[0148] The fourth determination submodule is used to determine, for any selected observation azimuth, the largest first resolution candidate terrain altitude angle at the selected observation azimuth as the first resolution terrain maximum altitude angle corresponding to the selected observation azimuth.

[0149] The first resolution database construction submodule is used to construct a first resolution terrain maximum altitude angle-occlusion point database based on all first resolution terrain maximum altitude angles of the first target data point within a preset azimuth range, and the first resolution occlusion points corresponding to the first resolution terrain maximum altitude angle.

[0150] According to an embodiment of the present disclosure, the simulation module 840 may include a fifth determination submodule, a sixth determination submodule, a seventh determination submodule, an eighth determination submodule and a simulation submodule.

[0151] The fifth determination submodule is used to determine the distance between the sun and the celestial body and the longitude information and latitude information of the direct point of the sun at a selected time according to the planetary ephemeris related to the celestial body.

[0152] The sixth determination submodule is used to determine the solar altitude angle and solar azimuth angle of each first target data point according to the first longitude information, the first latitude information, the direct point longitude information and the direct point latitude information.

[0153] The seventh determination submodule is used to determine the solar visibility coefficient of each first target data point according to the first resolution terrain maximum altitude angle-obstruction point database of each first target data point, as well as the solar altitude angle and the solar azimuth angle, wherein the solar visibility coefficient represents the ratio of the visible area of ​​the solar disk observed at the first target data point to the area of ​​the solar disk.

[0154] The eighth determination submodule is used to determine the effective solar irradiance of each first target data point for each first target data point based on the sun-celestial body distance, the solar visibility coefficient, the solar altitude angle, the solar azimuth angle, and the terrain slope angle and slope azimuth angle at the first target observation point.

[0155] The simulation submodule is used to simulate the illumination of the celestial body surface according to the effective solar irradiance of each first target data point to obtain a celestial body illumination simulation image.

[0156] According to an embodiment of the present disclosure, the seventh determination submodule may include a reference unit, a visible area determination unit, and a solar visibility coefficient determination unit.

[0157] A reference unit is used to determine a first reference terrain maximum altitude angle and a second reference terrain maximum altitude angle in a first resolution terrain maximum altitude angle-occlusion point database based on the solar azimuth, wherein the first reference terrain maximum altitude angle corresponds to a first observation azimuth, and the second reference terrain maximum altitude angle corresponds to a second observation azimuth, and the first observation azimuth and the second observation azimuth are respectively located on both sides of the solar azimuth.

[0158] The visible area determination unit is used to divide the visible area of ​​the sun disk based on the maximum altitude angle of the first adjacent terrain and the maximum altitude angle of the second adjacent terrain.

[0159] The solar visibility coefficient determination unit is used to determine the solar visibility coefficient based on the ratio of the visible region area of ​​the solar disk to the total area of ​​the solar disk.

[0160] According to an embodiment of the present disclosure, the eighth determination submodule may include a solar incident angle determination unit and a solar irradiance determination unit.

[0161] The solar incidence angle determination unit is used to determine the solar incidence angle at each first target data point according to the solar altitude angle, the solar azimuth angle, and the terrain slope angle and the slope azimuth angle at the first target observation point.

[0162] The solar irradiance determination unit is used to determine the effective solar irradiance of each first target data point according to the solar visibility coefficient, the solar incident angle at each first target data point and the distance between the sun and the celestial body.

[0163] According to the embodiments of the present invention, any one or more of the modules, submodules, units, and subunits, or at least part of the functions of any one of them can be implemented in one module. According to the embodiments of the present invention, any one or more of the modules, submodules, units, and subunits can be split into multiple modules for implementation. According to the embodiments of the present invention, any one or more of the modules, submodules, units, and subunits can be at least partially implemented as hardware circuits, such as field programmable gate arrays (FPGAs), programmable logic arrays (PLAs), systems on chips, systems on substrates, systems on packages, application specific integrated circuits (ASICs), or can be implemented by hardware or firmware in any other reasonable way of integrating or packaging circuits, or can be implemented in any one of the three implementation methods of software, hardware, and firmware, or in any appropriate combination of any of them. Alternatively, according to the embodiments of the present invention, one or more of the modules, submodules, units, and subunits can be at least partially implemented as computer program modules, and when the computer program modules are run, the corresponding functions can be performed.

[0164] For example, any multiple of the determination module 810, the first construction module 820, the second construction module 830, and the simulation module 840 can be combined in one module / unit / sub-unit for implementation, or any one of the modules / units / sub-units can be split into multiple modules / units / sub-units. Alternatively, at least part of the functions of one or more of these modules / units / sub-units can be combined with at least part of the functions of other modules / units / sub-units and implemented in one module / unit / sub-unit. According to an embodiment of the present disclosure, at least one of the determination module 810, the first construction module 820, the second construction module 830, and the simulation module 840 can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on a chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or can be implemented by hardware or firmware such as any other reasonable way of integrating or packaging the circuit, or implemented in any one of the three implementation methods of software, hardware, and firmware, or in any appropriate combination of any of them. Alternatively, at least one of the determination module 810 , the first construction module 820 , the second construction module 830 , and the simulation module 840 may be at least partially implemented as a computer program module, and when the computer program module is executed, a corresponding function may be performed.

[0165] It should be noted that the celestial body surface illumination simulation device part in the embodiment of the present disclosure corresponds to the celestial body surface illumination simulation method part in the embodiment of the present disclosure. The description of the celestial body surface illumination simulation device part specifically refers to the celestial body surface illumination simulation method part, which will not be repeated here.

[0166] The flowcharts and block diagrams in the accompanying drawings illustrate the possible implementation architecture, functions and operations of the systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each box in the flowchart or block diagram may represent a module, a program segment, or a part of a code, and the above-mentioned module, program segment, or a part of the code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box may also occur in an order different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, and the combination of boxes in the block diagram or flowchart, can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions. It can be understood by those skilled in the art that the features recorded in the various embodiments of the present disclosure can be combined and / or combined in a variety of ways, even if such a combination or combination is not explicitly recorded in the present disclosure. In particular, without departing from the spirit and teaching of the present disclosure, the features described in the various embodiments of the present disclosure may be combined and / or combined in a variety of ways. All of these combinations and / or combinations fall within the scope of the present disclosure.

[0167] The embodiments of the present disclosure are described above. However, these embodiments are only for illustrative purposes and are not intended to limit the scope of the present disclosure. Although the embodiments are described above, this does not mean that the measures in the various embodiments cannot be used in combination to advantage. Without departing from the scope of the present disclosure, those skilled in the art may make a variety of substitutions and modifications, which should all fall within the scope of the present disclosure.

Claims

1. A method for simulating illumination of a celestial body surface, include: Converting initial terrain data at a first resolution and initial terrain data at a second resolution into a target coordinate system, respectively obtaining first resolution terrain data and second resolution terrain data, and determining a target area in the target coordinate system based on the first resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, where m ≥ 1 and is a positive integer; For each of the first target data points, construct a second resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to n second target data points corresponding to the first target data point, wherein the second resolution terrain maximum altitude angle-occlusion point database includes the second resolution terrain maximum altitude angle and second resolution occlusion point of the first target data point at any observation azimuth, and n≥1 is a positive integer; For each of the first target data points, constructing a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to a preset search window and the second resolution terrain maximum altitude angle-occlusion point database, wherein the first resolution terrain maximum altitude angle-occlusion point database includes the first resolution terrain maximum altitude angle and the first resolution occlusion point of the first target data point at any observation azimuth; and Based on the planetary ephemeris associated with the celestial body, the local terrain data of each of the first target data points and the first resolution terrain maximum altitude angle-occlusion point database of each of the first target data points, the illumination of the celestial body surface is simulated to obtain a celestial body illumination simulation image.

2. The method according to claim 1, in, The converting the initial terrain data at the first resolution and the initial terrain data at the second resolution into the target coordinate system to obtain the first resolution terrain data and the second resolution terrain data respectively comprises: Taking the celestial pole as the origin, constructing the target coordinate system; Converting the initial terrain data at the first resolution to the target coordinate system to obtain terrain data at the first resolution, wherein the terrain data at the first resolution includes M first terrain data points, each of which includes first longitude information, first latitude information, and first elevation information, where M≥1 and is a positive integer; and Convert the initial terrain data at the second resolution to the target coordinate system to obtain the second resolution terrain data, where the second resolution terrain data includes N second terrain data points, each of which includes second longitude information, second latitude information, and second elevation information, where N ≥ 1 and is a positive integer; Determining the target area in the target coordinate system based on the first-resolution terrain data includes: determining the target area within a candidate area defined by the M first terrain data points, wherein the target area corresponds to the m first target data points.

3. The method according to claim 2, in, For each of the first target data points, constructing a second-resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to the n second target data points corresponding to the first target data point comprises: Taking the first target data point as an observation point, searching for all second terrain data points falling within a preset azimuth range among the N second terrain data points based on a preset azimuth step size to obtain the n second target data points, where the preset azimuth range is [0, 2π]; Determine second-resolution candidate terrain altitude angles respectively between the first target data point and the n second target data points based on the first longitude information, the first latitude information, the first elevation information, the second longitude information, the second latitude information, and the second elevation information; For any selected observation azimuth, determining the largest second-resolution candidate terrain altitude angle at the selected observation azimuth as the second-resolution terrain maximum altitude angle corresponding to the selected observation azimuth; and A second resolution terrain maximum altitude angle-occlusion point database is constructed based on all second resolution terrain maximum altitude angles of the first target data point within the preset azimuth range and second resolution occlusion points corresponding to the second resolution terrain maximum altitude angle.

4. The method according to claim 3, in, For each second target data point, the determining, based on the first longitude information, the first latitude information, the first elevation information, the second longitude information, the second latitude information, and the second elevation information, a second resolution candidate terrain altitude angle between the first target data point and the n second target data points comprises: Determine a target distance between the first target data point and the second target data point based on the first longitude information, the first latitude information, the second longitude information, the second latitude information, and a radius of a celestial body; According to the target distance and the celestial body radius, respectively correcting the first elevation information and the second elevation information to obtain first elevation correction information and second elevation correction information; and A second resolution candidate terrain elevation angle between the first target data point and the second target data point is determined according to the target distance, the celestial body radius, the first elevation correction information, and the second elevation correction information.

5. The method according to claim 1, in, For each of the first target data points, constructing a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to a preset search window and the second resolution terrain maximum altitude angle-occlusion point database includes: Taking the first target data point as an observation point, for any selected observation azimuth, constructing the preset search window according to the second resolution occlusion point corresponding to the selected azimuth; Determine all the first target data points located in the search window at the selected observation azimuth as first resolution candidate occlusion points, and obtain p first resolution candidate occlusion points corresponding to the selected observation azimuth; For any selected observation azimuth, determining first resolution candidate terrain altitude angles respectively between the first target data point and the p first resolution candidate occlusion points corresponding to the selected observation azimuth; For any selected observation azimuth, determining the largest first-resolution candidate terrain altitude angle at the selected observation azimuth as the first-resolution terrain maximum altitude angle corresponding to the selected observation azimuth; and A first resolution terrain maximum altitude angle-occlusion point database is constructed based on all first resolution terrain maximum altitude angles of the first target data point within the preset azimuth range and first resolution occlusion points corresponding to the first resolution terrain maximum altitude angle.

6. The method according to claim 5, in, The preset search window is centered on the second resolution occlusion point and is constructed based on a preset window side length, and the preset window side length is determined based on the size of the second resolution.

7. The method according to any one of claims 1 to 6, in, The local terrain data includes a terrain slope angle and a slope azimuth, and the simulating the illumination of the celestial body surface based on a planetary ephemeris related to the celestial body, the local terrain data of each of the first target data points, and a first resolution terrain maximum altitude angle-blocking point database of each of the first target data points to obtain a celestial body illumination simulation image includes: Determine the distance between the sun and the celestial body and the longitude information and latitude information of the direct point of the sun at the selected time according to the planetary ephemeris related to the celestial body; Determine the solar altitude angle and solar azimuth angle of each of the first target data points according to the first longitude information, the first latitude information, the direct point longitude information, and the direct point latitude information; Determine a solar visibility coefficient for each of the first target data points according to the first resolution terrain maximum altitude angle-blocking point database for each of the first target data points, the solar altitude angle and the solar azimuth angle, wherein the solar visibility coefficient represents the ratio of the visible area of ​​the solar disk observed at the first target data point to the area of ​​the solar disk; For each of the first target data points, determining the effective solar irradiance of each of the first target data points according to the sun-celestial body distance, the solar visibility coefficient, the solar altitude angle, the solar azimuth angle, and the terrain slope angle and slope azimuth angle at the first target observation point; and According to the effective solar irradiance of each of the first target data points, the illumination on the surface of the celestial body is simulated to obtain a celestial body illumination simulation image.

8. The method according to claim 7, in, For each of the first target data points, determining the solar visibility coefficient of each of the first target data points according to the first resolution terrain maximum altitude angle-blocking point database of each of the first target data points, and the solar altitude angle and the solar azimuth angle comprises: Based on the solar azimuth, determining a first reference terrain maximum altitude angle and a second reference terrain maximum altitude angle in the first resolution terrain maximum altitude angle-blocking point database, wherein the first reference terrain maximum altitude angle corresponds to a first observation azimuth, the second reference terrain maximum altitude angle corresponds to a second observation azimuth, and the first observation azimuth and the second observation azimuth are respectively located on both sides of the solar azimuth; Dividing the visible area of ​​the sun disk based on the maximum altitude angle of the first adjacent terrain and the maximum altitude angle of the second adjacent terrain; and The solar visibility coefficient is determined based on the ratio of the visible region area of ​​the solar disk to the total area of ​​the solar disk.

9. The method according to claim 8, in, Determining the effective solar irradiance of each of the first target data points according to the sun-celestial body distance, the solar visibility coefficient, the solar altitude angle, the solar azimuth angle, and the terrain slope angle and slope azimuth angle at the first target observation point comprises: Determine the solar incident angle at each of the first target data points according to the solar altitude angle, the solar azimuth angle, and the terrain slope angle and slope azimuth angle at the first target observation point; and The effective solar irradiance of each of the first target data points is determined according to the solar visibility coefficient, the solar incident angle at each of the first target data points and the sun-celestial body distance.

10. A celestial body surface illumination simulation device, include: A first determination module is used to convert the initial terrain data at the first resolution and the initial terrain data at the second resolution into a target coordinate system, obtain first resolution terrain data and second resolution terrain data respectively, and determine a target area in the target coordinate system based on the first resolution terrain data, wherein the first resolution is higher than the second resolution, and there are m first target data points corresponding to the target area, where m ≥ 1 and is a positive integer; A second determination module is used to construct, for each of the first target data points, a second resolution terrain maximum altitude angle-blocking point database related to the first target data point according to n second target data points corresponding to the first target data point, where n≥1 and is a positive integer; A third determination module is used to construct, for each of the first target data points, a first resolution terrain maximum altitude angle-occlusion point database related to the first target data point according to a preset search window and the second resolution terrain maximum altitude angle-occlusion point database; and A simulation module is used to simulate the illumination of the celestial body surface based on a planetary ephemeris associated with the celestial body, local terrain data of each of the first target data points, and a first-resolution terrain maximum altitude angle-occlusion point database of each of the first target data points to obtain a celestial body illumination simulation image.

Images