Lunar surface illumination condition parallel analysis method

Abstract

The embodiment of the invention discloses a lunar surface illumination condition parallel analysis method. The method comprises the steps that lunar surface illumination calculation correlation parameters of a lunar surface set area are configured; the lunar surface illumination calculation correlation parameters comprise lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters; parallel computing resources are distributed to lunar surface illumination computing correlation parameters of a lunar surface set area; the parallel computing resources comprise parallel thread resources and memory resources; according to the parallel computing resources and the lunar surface illumination computing correlation parameters, converting the pixel coordinates and the latitude and longitude coordinates of the set area of the lunar surface to obtain converted latitude and longitude coordinates; according to the parallel computing resources, the lunar surface illumination computing correlation parameters and the converted longitude and latitude coordinates, illumination computing results of all points in the set area of the lunar surface are computed in parallel. According to the technical scheme of the embodiment of the invention, the calculation efficiency of the lunar surface illumination condition can be improved, and the efficiency of future lunar exploration task planning is improved.

Description

Technical Field

[0001] Embodiments of the present invention relate to the field of aerospace technology, and in particular to a method, device, electronic device, storage medium, and program for parallel analysis of lunar illumination conditions. Background Art

[0002] In lunar exploration missions, especially for soft landing missions on the lunar surface (referred to as the lunar surface) and the construction of future lunar scientific research stations, the analysis of lunar lighting conditions is of great significance for site selection and overall mission planning.

[0003] At present, when analyzing the lighting conditions on the lunar surface, a serial algorithm is usually used, that is, the lighting conditions of each point on the lunar surface are calculated in turn.

[0004] During the development of this invention, the inventors discovered that illumination conditions at various points on the lunar surface vary over time due to factors such as the complex and rugged lunar topography and the changing position of the sun relative to the moon. Traditional methods for analyzing lunar illumination conditions rely on serial calculations, which are inefficient when considering large-scale and long-term illumination variations on the lunar surface, and therefore cannot meet the needs of real-time lunar mission planning. Summary of the Invention

[0005] Embodiments of the present invention provide a method, device, electronic device, storage medium, and program for parallel analysis of lunar illumination conditions, which can efficiently analyze large-scale and long-period lunar illumination conditions, improve the efficiency of calculating lunar surface illumination conditions, and thereby improve the efficiency of future lunar exploration mission planning.

[0006] According to one aspect of the present invention, a method for parallel analysis of lunar illumination conditions is provided, comprising:

[0007] Configuring lunar surface illumination calculation parameters for a set area on the lunar surface; wherein the lunar surface illumination calculation parameters include lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters;

[0008] Allocating parallel computing resources to lunar illumination calculation-related parameters for the set area on the lunar surface; wherein the parallel computing resources include parallel thread resources and memory resources;

[0009] Converting pixel coordinates and longitude and latitude coordinates of a set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation-related parameters to obtain converted longitude and latitude coordinates;

[0010] The illumination calculation results of each point in the set area on the lunar surface are calculated in parallel based on the parallel computing resources, the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates.

[0011] According to another aspect of the present invention, there is provided a device for parallel analysis of lunar illumination conditions, comprising:

[0012] A lunar illumination calculation parameter configuration module is used to configure lunar illumination calculation parameters for a set area on the lunar surface; wherein the lunar illumination calculation parameters include lunar illumination calculation condition parameters and lunar terrain digital elevation parameters;

[0013] a parallel computing resource allocation module, configured to allocate parallel computing resources to the lunar illumination calculation-related parameters of the set area on the lunar surface; wherein the parallel computing resources include parallel thread resources and memory resources;

[0014] a pixel longitude and latitude coordinate conversion module, configured to convert the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation-related parameters to obtain converted longitude and latitude coordinates;

[0015] The illumination calculation result calculation module is used to calculate the illumination calculation results of each point in the set area on the lunar surface according to the parallel computing resources, the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates.

[0016] According to another aspect of the present invention, an electronic device is provided, comprising:

[0017] at least one processor; and

[0018] a memory communicatively connected to the at least one processor; wherein,

[0019] The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the parallel analysis method of lunar illumination conditions described in any embodiment of the present invention.

[0020] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the parallel analysis method for lunar illumination conditions described in any embodiment of the present invention when executed.

[0021] According to another aspect of the present invention, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the parallel analysis method for lunar illumination conditions according to any embodiment of the present invention.

[0022] The embodiment of the present invention configures lunar surface illumination calculation condition parameters and lunar surface topography digital elevation parameters and other lunar surface illumination calculation-related parameters of a set lunar surface area to allocate parallel computing resources such as parallel thread resources and memory resources to the lunar surface illumination calculation-related parameters of the set lunar surface area, thereby converting the pixel coordinates and longitude and latitude coordinates of the set lunar surface area according to the parallel computing resources and the lunar surface illumination calculation-related parameters to obtain converted longitude and latitude coordinates. Finally, the illumination calculation results of each point in the set lunar surface area are calculated in parallel according to the parallel computing resources, the lunar surface illumination calculation-related parameters and the converted longitude and latitude coordinates. This can solve the problem of low computational efficiency in the existing serial calculation of lunar surface illumination conditions, can efficiently analyze the illumination conditions of the lunar surface over a large area and over a long period, improve the computational efficiency of the lunar surface illumination conditions, and thereby improve the efficiency of future lunar exploration mission planning.

[0023] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it intended to limit the scope of the present invention. Other features of the present invention will become readily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0025] Figure 1 This is a flow chart of a parallel analysis method for lunar illumination conditions provided by an embodiment of the present invention;

[0026] Figure 2 This is a flow chart of another method for parallel analysis of lunar illumination conditions provided by an embodiment of the present invention;

[0027] Figure 3 2 is a schematic diagram of a parallel analysis device for lunar illumination conditions provided by an embodiment of the present invention;

[0028] Figure 4 A schematic structural diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0029] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of the present invention.

[0030] It should be noted that the terms "first", "second", etc. in the description and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that the numbers used in this way can be interchanged where appropriate so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0031] Figure 1 This is a flow chart of a method for parallel analysis of lunar illumination conditions provided by an embodiment of the present invention. This embodiment is applicable to the case of parallel calculation of illumination conditions for each point in a set area on the lunar surface. The method can be executed by a device for parallel analysis of lunar illumination conditions. The device can be implemented by software and / or hardware and can generally be integrated into an electronic device. The electronic device can be a terminal device or a server device. As long as it can execute the method for parallel analysis of lunar illumination conditions, the embodiment of the present invention does not limit the specific device type of the electronic device. Accordingly, if Figure 1 As shown, the method includes the following operations:

[0032] S110. Configure lunar surface illumination calculation associated parameters for a set area on the lunar surface; wherein the lunar surface illumination calculation associated parameters include lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters.

[0033] The lunar surface set area may be a region of the lunar surface for which illumination conditions need to be calculated. Embodiments of the present invention do not limit the specific size of the lunar surface set area. Lunar surface illumination calculation-related parameters may be relevant parameters required for reference when calculating lunar surface illumination conditions. Lunar surface illumination calculation condition parameters may be relevant parameters for illumination calculation conditions. Lunar surface topography digital elevation parameters may be preprocessing parameters for lunar surface digital elevation model (DEM) data.

[0034] Before calculating lunar surface illumination conditions, you can first determine the lunar surface area for which illumination conditions need to be calculated, which serves as the lunar surface setting area. Optionally, the lunar surface setting area can encompass a larger area of ​​the lunar surface. After determining the lunar surface setting area, you can configure lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters for the lunar surface setting area as parameters associated with lunar surface illumination calculation for the lunar surface setting area.

[0035] S120. Allocate parallel computing resources to the lunar illumination calculation-related parameters of the set area on the lunar surface; wherein the parallel computing resources include parallel thread resources and memory resources.

[0036] The parallel computing resources may be computing resources capable of computing the illumination conditions of the lunar surface in parallel, and the parallel thread resources may be thread resources capable of computing the illumination conditions of the lunar surface in parallel.

[0037] To enable parallel calculation of lunar illumination conditions, after configuring the lunar illumination calculation parameters for a specific area on the lunar surface, you can allocate parallel computing resources, such as parallel thread resources and memory resources, to these parameters. Parallel computing resources allow for simultaneous calculation of lunar illumination conditions at each point in the specified area.

[0038] S130. Convert the pixel coordinates and the longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation-related parameters to obtain converted longitude and latitude coordinates.

[0039] The converted longitude and latitude coordinates may be longitude and latitude coordinates calculated by converting pixel coordinates of a point in a set area on the lunar surface.

[0040] Correspondingly, after allocating corresponding parallel computing resources to the lunar surface illumination calculation-related parameters of the set area on the lunar surface, a mapping relationship between pixel coordinates and parallel computing resources can be established. Based on the mapping relationship between pixel coordinates and parallel computing resources, the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface are converted in combination with the parallel computing resources and the lunar surface illumination calculation-related parameters to obtain the converted longitude and latitude coordinates of each point in the set area on the lunar surface.

[0041] S140 , parallelly calculating illumination calculation results of each point in the set area on the lunar surface according to the parallel computing resources, the lunar surface illumination calculation associated parameters, and the converted longitude and latitude coordinates.

[0042] Correspondingly, after calculating the converted longitude and latitude coordinates of each point in the set area on the lunar surface, we can further use parallel computing methods based on parallel computing resources, lunar surface illumination calculation related parameters and converted longitude and latitude coordinates to synchronously calculate the illumination calculation results of each point in the set area on the lunar surface.

[0043] It can be seen that the above technical solution allocates parallel computing resources to the lunar surface illumination calculation-related parameters according to the specific configuration of the lunar surface illumination calculation-related parameters in the set area of ​​the lunar surface. After calculating the converted longitude and latitude coordinates of each point in the set area of ​​the lunar surface, the parallel computing resources are used to parallelly calculate the illumination calculation results of each point in the set area of ​​the lunar surface. This can greatly improve the calculation efficiency of the illumination calculation results of each point in the set area of ​​the lunar surface, thereby efficiently analyzing the illumination conditions of the lunar surface over a large range and long period, which will help to improve the efficiency of future lunar exploration mission planning.

[0044] The embodiment of the present invention configures lunar surface illumination calculation condition parameters and lunar surface topography digital elevation parameters and other lunar surface illumination calculation-related parameters of a set lunar surface area to allocate parallel computing resources such as parallel thread resources and memory resources to the lunar surface illumination calculation-related parameters of the set lunar surface area, thereby converting the pixel coordinates and longitude and latitude coordinates of the set lunar surface area according to the parallel computing resources and the lunar surface illumination calculation-related parameters to obtain converted longitude and latitude coordinates. Finally, the illumination calculation results of each point in the set lunar surface area are calculated in parallel according to the parallel computing resources, the lunar surface illumination calculation-related parameters and the converted longitude and latitude coordinates. This can solve the problem of low computational efficiency in the existing serial calculation of lunar surface illumination conditions, can efficiently analyze the illumination conditions of the lunar surface over a large area and over a long period, improve the computational efficiency of the lunar surface illumination conditions, and thereby improve the efficiency of future lunar exploration mission planning.

[0045] Figure 2 This is a flow chart of another parallel analysis method for lunar illumination conditions provided by an embodiment of the present invention. This embodiment is specific based on the above embodiment. In this embodiment, multiple specific optional implementation methods are provided for configuring lunar illumination calculation parameters associated with a set area on the lunar surface, allocating parallel computing resources, converting pixel coordinates and latitude and longitude coordinates of the set area on the lunar surface, and parallel computing illumination calculation results for each point in the set area on the lunar surface. Figure 2 As shown, the method of this embodiment may include:

[0046] S210. Configure lunar surface illumination calculation associated parameters for a set area on the lunar surface; wherein the lunar surface illumination calculation associated parameters include lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters.

[0047] In an optional embodiment of the present invention, the configuration of lunar surface illumination calculation associated parameters for a set lunar surface area may include: configuring the lunar surface longitude range, lunar surface latitude range, terrain search range, time span, longitude resolution, latitude resolution and time resolution of the set lunar surface area as the lunar surface illumination calculation condition parameters; configuring the first lunar surface terrain digital elevation parameters of the set lunar surface area for the equidistant cylindrical projection mode; wherein the first lunar surface terrain digital elevation parameters include a first resolution, a scale factor, an offset of the projection system origin relative to the target point, the origin longitude, the origin latitude, a first number of rows and a first number of columns; and / or, configuring the second lunar surface terrain digital elevation parameters of the set lunar surface area for the polar stereographic projection mode; wherein the second lunar surface terrain digital elevation parameters include a second resolution, a scale factor, a second number of rows and a pole identifier.

[0048] The lunar longitude range may be the longitude range of a set area on the lunar surface. The lunar latitude range may be the latitude range of a set area on the lunar surface. The terrain search range may be the distance searched outward from a certain point on the lunar surface. The first lunar terrain digital elevation parameter may be a lunar terrain digital elevation parameter configured for an equidistant cylindrical projection. The first resolution may be a map resolution. The first number of rows may be a row number parameter involved in the first lunar terrain digital elevation parameter. The second lunar terrain digital elevation parameter may be a lunar terrain digital elevation parameter configured for a polar stereographic projection. The second resolution may be a map scale. The second number of rows may be a row number parameter involved in the second lunar terrain digital elevation parameter.

[0049] Specifically, you can first set the lighting calculation conditions, and configure the lunar surface longitude range, lunar surface latitude range, terrain search range, time span, longitude resolution, latitude resolution and time resolution of the set area on the lunar surface as lunar surface lighting calculation condition parameters.

[0050] In a specific example, the lunar longitude range can be set to [ ], the lunar latitude range can be set to [ ], the terrain search range is set to R, the time span is set to D, the longitude resolution is set to N_lon, the latitude resolution is set to N_lat, and the time resolution is set to N_t.

[0051] At the same time, different types of lunar topography digital elevation parameters can be set according to different projection methods. Optionally, for the equidistant cylindrical projection method, the first resolution, scale factor, offset of the projection system origin relative to the target point, origin longitude, origin latitude, first number of rows, and number of columns can be configured as the first lunar topography digital elevation parameters corresponding to the equidistant cylindrical projection method. For the polar stereographic projection method, the second resolution, scale factor, second number of rows, and pole identifier can be configured as the second lunar topography digital elevation parameters corresponding to the polar stereographic projection method.

[0052] In a specific example, continuing with the above example, after configuring the lunar illumination calculation condition parameters, you can continue to set the lunar terrain digital elevation parameters according to different projection methods. Optionally, for the equidistant cylindrical projection method, you can configure the first resolution MAP_RESOLUTION, the scaling factor SCALING_FACTOR, the origin longitude CENTER_LONGITUDE, the origin latitude CENTER_LATITUDE, the first number of rows LINES, and the number of columns LINE_SAMPLES. The offset of the projection system origin relative to the target point can be the offset of the projection system origin relative to the first grid point in the upper left corner, which can include two offsets, namely the horizontal offset SAMPLE_PROJECTION_OFFSET and the vertical offset LINE_PROJECTION_OFFSET.

[0053] Optionally, for the polar stereographic projection mode, the second resolution MAP_SCALE, the scaling factor SCALING_FACTOR, the second number of rows N and the pole identifier NS can be configured, where the South Pole NS=-1 and the North Pole NS=1.

[0054] S220: Determine resource allocation reference parameters according to the target projection method used for the set area on the lunar surface.

[0055] The target projection method may be the projection method used for lunar illumination analysis. The resource allocation reference parameter may be a parameter used to allocate parallel computing resources. It is understood that different target projection methods correspond to different types of resource allocation reference parameters.

[0056] When allocating parallel computing resources, we can first determine the target projection method used when performing lunar illumination analysis on a set area on the lunar surface, and then determine the resource allocation reference parameters for allocating parallel computing resources based on the target projection method.

[0057] S230 , using a two-dimensional thread block, allocating the memory resources for the lunar topography digital elevation parameters in a processor video memory according to the resource allocation reference parameters.

[0058] A two-dimensional thread block refers to a thread block with multiple threads in both the x and y directions, forming a two-dimensional array. Each two-dimensional thread block contains multiple threads, which are distributed to different stream processors on a processor such as a GPU (Graphics Processing Unit) during execution. Each thread executes the same instructions but operates on different data, thus achieving parallel computing. Processor memory refers to the processor's video memory. Exemplary processor types include, but are not limited to, CPUs (Central Processing Units) and GPUs.

[0059] After determining the resource allocation reference parameters according to the target projection method used in the set area on the lunar surface, a two-dimensional thread block can be used to allocate memory resources for the lunar surface terrain digital elevation parameters in the processor memory according to the resource allocation reference parameters through a function function that dynamically allocates memory.

[0060] For example, if the target projection method is equidistant cylindrical projection, the resource allocation reference parameters can be the first number of rows LINES and columns LINE_SAMPLES in the first lunar surface topography digital elevation parameter. Accordingly, a two-dimensional thread block can be used to allocate memory resources for the first lunar surface topography digital elevation parameter in the GPU memory through a dynamic memory allocation function. For example, for the equidistant cylindrical projection method, the memory size that can be allocated is: ,in, It is a memory capacity measurement function used to calculate the size of a variable.

[0061] For example, if the target projection mode is polar stereographic, the resource allocation reference parameter can be the second row number N in the second lunar topography digital elevation parameter. Accordingly, a two-dimensional thread block can be used to allocate memory resources for the second lunar topography digital elevation parameter in the GPU memory through a dynamic memory allocation function. For example, for polar stereographic projection, the memory size that can be allocated is: .

[0062] S240 . Allocate the parallel thread resources to the set area on the lunar surface according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block.

[0063] In a specific example, assuming the pixel coordinates of a set area on the lunar surface are (SAMPLE, LINE) and the thread index of a two-dimensional thread block is (I, J), then the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index is I = SAMPLE, J = LINE. Here, I = threadIdx.x + blockIdx.x * blockDim.x, and J = threadIdx.y + blockIdx.y * blockDim.y. threadIdx.x and threadIdx.y are the thread index values ​​of a specified thread in the two-dimensional thread block, blockIdx.x and blockIdx.y are the indices of the two-dimensional thread block, and blockDim.x and blockDim.y are the dimensions of the threads in the two-dimensional thread block. Since the number of pixels in the set area on the lunar surface is fixed, its pixel coordinates are also fixed values. Accordingly, based on the known pixel coordinates (SAMPLE, LINE), the thread index of the two-dimensional thread block corresponding to each pixel coordinate can be calculated sequentially. Therefore, the thread indices of all two-dimensional thread blocks corresponding to the pixel coordinates of all points in the entire set area on the lunar surface are the total amount of parallel thread resources that need to be allocated to the set area on the lunar surface. In other words, each point in the set area on the lunar surface corresponds to one thread resource.

[0064] S250: Determine the target projection method used for the set area on the lunar surface.

[0065] S260: Convert the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block, and the lunar topography digital elevation parameters corresponding to the target projection mode, to obtain the converted longitude and latitude coordinates.

[0066] It is understandable that different target projection methods are used for setting areas on the lunar surface, and the methods of converting their pixel coordinates and longitude and latitude coordinates are also different.

[0067] In a specific example, continuing with the above example, for the equidistant cylindrical projection method, according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block, and the first lunar surface topography digital elevation parameter corresponding to the equidistant cylindrical projection method, the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface are converted to obtain the converted longitude and convert latitude They are:

[0068]

[0069] I=threadIdx.x+blockIdx.x*blockDim.x

[0070] J=threadIdx.y+blockIdx.y*blockDim.y

[0071] For the polar stereographic projection, the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface are converted according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block, as well as the second lunar surface topography digital elevation parameters corresponding to the polar stereographic projection. The converted longitude is and convert latitude They are:

[0072]

[0073] I=threadIdx.x+blockIdx.x*blockDim.x

[0074] J=threadIdx.y+blockIdx.y*blockDim.y

[0075] Where R represents the terrain search range. If the pole is marked as the North Pole, NS = 1; if the pole is marked as the South Pole, NS = -1.

[0076] S270: Calculate the maximum altitude angle between points in the set area on the lunar surface based on the lunar illumination calculation associated parameters and the converted longitude and latitude coordinates.

[0077] Furthermore, according to the illumination calculation conditions set in the lunar illumination calculation associated parameters, the lunar longitude range is calculated in a processor such as a GPU. ] and lunar latitude range[ ] the altitude angles between all points on the lunar surface , and count the maximum altitude angles of each point in different directions according to the azimuth, and form an altitude angle database based on the maximum altitude angles between each point in a set area on the lunar surface.

[0078] In an optional embodiment of the present invention, the calculation of the maximum altitude angle between each point in the set area of ​​the lunar surface based on the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates may include: based on the time resolution in the lunar surface illumination calculation associated parameters, with each point where each thread in the parallel thread resource is located as the center, searching outward the distance of the terrain search range according to the azimuth calculation interval and the distance interval, and calculating the corresponding geodetic coordinates of each point at equal distances by interpolation; wherein the geodetic coordinates include the converted longitude and latitude coordinates and elevation; converting the geodetic coordinates of each point in the set area of ​​the lunar surface into a three-dimensional position vector of lunar fixed coordinates; parallel calculation of the altitude angles of other points in each azimuth direction relative to the points where the threads are located within the terrain search range centered on the points where the threads are located according to the thread index of each thread in the parallel thread resource; performing reduction calculations in parallel according to the azimuth angles of each point in the set area of ​​the lunar surface, and selecting the maximum altitude angle of the azimuth angles of each point in the set area of ​​the lunar surface.

[0079] The azimuth angle calculation interval may be an interval for calculating the azimuth angle. For example, assuming that the azimuth angle calculation interval is 5°, it may mean that the azimuth angle is calculated once every 5°.

[0080] Specifically, we can first calculate the time resolution of the parameters associated with the lunar illumination calculation, take the points where the threads in the processor (e.g., GPU) are located as the center, and search the terrain search range R from 0-360° in the direction of the azimuth angle based on the calculation interval and the distance interval, and calculate the corresponding longitude, latitude and elevation of each point in the set area on the lunar surface by equidistant interpolation. The azimuth angle calculation interval can be ,in is the lunar rotation speed, and N_t is the time resolution. The distance interval can be determined according to the DEM resolution or a specified distance resolution. Optionally, the DEM resolution can be the first resolution MAP_RESOLUTION or the second resolution MAP_SCALE. For example, the specific process of calculating the corresponding latitude, longitude, and elevation of each point in a set area on the lunar surface by equidistant interpolation can be:

[0081] S1. Use a dynamic memory allocation function to allocate memory for the altitude angles of various azimuths in the GPU.

[0082] S2. In the GPU, use the threads corresponding to each point to parallelly calculate the latitude, longitude and elevation of each point within the terrain search range R. Assume that the center grid point of the projection coordinate system corresponding to the target projection method is marked as O, and its geodetic coordinates are ,in, , To directly read the elevation value of DEM, SCALING_FACTOR is the scaling factor. The grid point at the distance H from the center grid point in the azimuth α direction is denoted as M, and its geodetic coordinates are Optionally, the above-mentioned geodetic coordinates can be converted longitude coordinates and converted latitude coordinates. Then search for n known points near point M, and the elevation of point M is: .

[0083] in, is the distance from point M to n nearby known points, is the elevation value of n nearby known points, and m is the degree of fit. Optionally, the lunar surface digital elevation data, i.e., DEM data, can be read in the GPU by pixel coordinates and thread index to obtain the corresponding elevation value of each point on the lunar surface. .

[0084] Furthermore, according to the requirements of longitude and latitude resolution, the geodetic coordinates of each point within the specified longitude and longitude range corresponding to the set area on the lunar surface are converted into three-dimensional position vectors in the lunar fixed system. Assume that the geodetic coordinates of each point in the set area on the lunar surface are , are the converted longitude, converted latitude and elevation respectively, then the position vector of each point in the set area on the lunar surface in the lunar fixed system can be:

[0085]

[0086] in, , is the equatorial radius of the lunar reference ellipsoid, is the flattening of the reference ellipsoid.

[0087] Furthermore, each thread calculates in parallel the altitude angles of other points relative to the point in each azimuth direction, centered on the thread's location and within the terrain search range R, according to the thread index. The position vector of point i in the lunar fixed system is , the position vector of point j is First, the position vector of point j relative to point i is The conversion method from the lunar solid system to the horizontal coordinate system with a point on the lunar surface as the origin is: .in, and Represents the rotation matrices around the y-axis and around the z-axis respectively, and are the converted longitude and latitude of point i respectively. , then the altitude angle of point j relative to point i is: .

[0088] After obtaining the altitude angle between each point, further reduction calculations can be performed in the GPU according to the azimuth, that is, the maximum altitude angle in each azimuth direction is selected to form a database of the maximum altitude angle on the lunar surface. Specifically, all points on the same azimuth can be divided into blocks, and the altitude angles of the points where each thread is located can be compared two by two. The larger of the two values ​​can be taken and stored in a temporary variable, and all threads in the thread block are synchronized through thread synchronization commands used to ensure data consistency and collaborative calculations within the thread block. Then, the step size is halved, and the above reduction process is repeated for the temporary variables to obtain the maximum altitude angle corresponding to the azimuth as the maximum altitude angle in that direction. Similarly, the above calculations are performed in parallel on the GPU for different azimuths to obtain the maximum altitude angle in each azimuth direction and form a database of altitude angles.

[0089] S280 . Divide the parallel thread resources into blocks at different times according to the time resolution in the lunar illumination calculation-related parameters.

[0090] Temporal resolution refers to the minimum time interval between two illumination calculations performed on the same area. Based on the temporal resolution and time span D in the lunar illumination calculation parameters, parallel thread resources are divided into blocks at different times. Specifically, within time span D, parallel thread resources are divided into blocks based on temporal resolution, with each moment corresponding to a group of thread resources.

[0091] S290: The threads at each moment obtain the solar altitude angle and the solar azimuth angle in parallel, and calculate the maximum altitude angle corresponding to the solar azimuth angle by interpolation in parallel based on the maximum altitude angle between points in the set area on the lunar surface.

[0092] Optionally, at each time t, the thread can call the precise ephemeris to obtain the sun's position at time t, and convert the sun's position at time t to obtain the sun's altitude angle. and the solar azimuth Furthermore, the maximum altitude angle between each point in the set area on the lunar surface calculated by the above scheme is used, and the threads at each moment perform parallel interpolation to calculate the solar azimuth angle. Corresponding maximum altitude angle .

[0093] S2110. The threads at each moment compare in parallel the solar altitude angle and the maximum altitude angle corresponding to the solar altitude angle to obtain illumination calculation results for each point in the set area on the lunar surface.

[0094] Furthermore, the threads at each moment compare the solar altitude angle and the maximum horizontal altitude angle corresponding to the solar altitude angle in parallel. Suppose the solar altitude angle of a point in the set area on the lunar surface is The maximum altitude angle at this point is ,like , then the point is illuminated; if , then the point is shaded; if , then the point is the penumbra. Among them, is the apparent radius of the solar disk.

[0095] This allows the processor to parallelize the calculation of the solar altitude and azimuth for all lunar points within a specified time span, D, and a specified lunar longitude and latitude range using parallel thread resources. By comparing these calculations with the corresponding maximum altitude angle for each point, the lunar illumination conditions for each point within the specified time span and lunar range can be obtained. Furthermore, the calculated illumination conditions for a set area of ​​the lunar surface can be output, for example, by copying the illumination condition calculation results from the GPU to the CPU output.

[0096] The above technical solution provides a high-performance parallel analysis method for lunar illumination conditions. It uses parallel computing resources to parallelize the lunar illumination calculation process, which can efficiently analyze large-scale and long-period lunar illumination conditions, improve the calculation efficiency of lunar surface illumination conditions, and thus improve the efficiency of future lunar exploration mission planning.

[0097] In the technical solutions disclosed herein, the collection, storage, use, processing, transmission, provision and disclosure of user personal information are in compliance with the relevant laws and regulations and do not violate public order and good morals.

[0098] It should be noted that the relevant information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for display, data for analysis, etc.) involved in this disclosure are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data comply with the relevant laws, regulations and standards of the relevant regions.

[0099] It should be noted that any arrangement and combination of the technical features in the above embodiments also falls within the protection scope of the present invention.

[0100] Figure 3 FIG. 1 is a schematic diagram of a parallel analysis device for lunar illumination conditions provided by an embodiment of the present invention. Figure 3 As shown, the apparatus includes: a lunar illumination calculation associated parameter configuration module 310, a parallel computing resource allocation module 320, a pixel latitude and longitude coordinate conversion module 330, and an illumination calculation result calculation module 340, wherein:

[0101] A lunar illumination calculation parameter configuration module 310 is used to configure lunar illumination calculation parameters for a set area on the lunar surface; wherein the lunar illumination calculation parameters include lunar illumination calculation condition parameters and lunar terrain digital elevation parameters;

[0102] A parallel computing resource allocation module 320 is configured to allocate parallel computing resources to the lunar illumination calculation parameters of the set lunar surface area; wherein the parallel computing resources include parallel thread resources and memory resources;

[0103] A pixel longitude and latitude coordinate conversion module 330 is configured to convert the pixel coordinates and the longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation-related parameters to obtain converted longitude and latitude coordinates;

[0104] The illumination calculation result calculation module 340 is used to calculate the illumination calculation result of each point in the set area on the lunar surface according to the parallel computing resources, the lunar surface illumination calculation related parameters and the converted longitude and latitude coordinates.

[0105] The embodiment of the present invention configures lunar surface illumination calculation condition parameters and lunar surface topography digital elevation parameters and other lunar surface illumination calculation-related parameters of a set lunar surface area to allocate parallel computing resources such as parallel thread resources and memory resources to the lunar surface illumination calculation-related parameters of the set lunar surface area, thereby converting the pixel coordinates and longitude and latitude coordinates of the set lunar surface area according to the parallel computing resources and the lunar surface illumination calculation-related parameters to obtain converted longitude and latitude coordinates. Finally, the illumination calculation results of each point in the set lunar surface area are calculated in parallel according to the parallel computing resources, the lunar surface illumination calculation-related parameters and the converted longitude and latitude coordinates. This can solve the problem of low computational efficiency in the existing serial calculation of lunar surface illumination conditions, can efficiently analyze the illumination conditions of the lunar surface over a large area and over a long period, improve the computational efficiency of the lunar surface illumination conditions, and thereby improve the efficiency of future lunar exploration mission planning.

[0106] Optionally, the lunar illumination calculation associated parameter configuration module 310 is further configured to: configure the lunar surface longitude range, lunar surface latitude range, terrain search range, time span, longitude resolution, latitude resolution, and time resolution of the set lunar surface area as the lunar surface illumination calculation condition parameters;

[0107] For the equidistant cylindrical projection method, configuring first lunar surface topographic digital elevation parameters of the set area on the lunar surface; wherein the first lunar surface topographic digital elevation parameters include a first resolution, a scale factor, an offset of the projection system origin relative to the target point, the origin longitude, the origin latitude, a first number of rows and columns; and / or

[0108] For the polar stereographic projection mode, second lunar surface topographic digital elevation parameters of the set area on the lunar surface are configured; wherein the second lunar surface topographic digital elevation parameters include a second resolution, a scale factor, a second number of rows, and a pole identifier.

[0109] Optionally, the parallel computing resource allocation module 320 is further used to: determine resource allocation reference parameters based on the target projection method adopted by the set area on the lunar surface; use a two-dimensional thread block to allocate the memory resources for the lunar surface terrain digital elevation parameters in the processor video memory according to the resource allocation reference parameters; and allocate the parallel thread resources to the set area on the lunar surface according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block.

[0110] Optionally, the pixel longitude and latitude coordinate conversion module 330 is also used to: determine the target projection method adopted for the set area on the lunar surface; convert the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block, and the lunar topography digital elevation parameters corresponding to the target projection method to obtain the converted longitude and latitude coordinates.

[0111] Optionally, the illumination calculation result calculation module 340 is also used to: calculate the maximum altitude angle between each point in the set area of ​​the lunar surface according to the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates; divide the parallel thread resources into blocks according to different times according to the time resolution in the lunar surface illumination calculation associated parameters; the threads at each moment obtain the solar altitude angle and the solar azimuth angle in parallel, and according to the maximum altitude angle between each point in the set area of ​​the lunar surface, parallelly interpolate and calculate the maximum altitude angle corresponding to the solar azimuth angle; the threads at each moment compare the solar altitude angle and the maximum altitude angle corresponding to the solar altitude angle in parallel to obtain the illumination calculation results of each point in the set area of ​​the lunar surface.

[0112] Optionally, the illumination calculation result calculation module 340 is also used to: according to the time resolution in the lunar illumination calculation associated parameters, with each point where each thread in the parallel thread resource is located as the center, search the distance of the terrain search range outward according to the azimuth calculation interval and the distance interval, and calculate the corresponding geodetic coordinates of each point at equal distances by interpolation; wherein the geodetic coordinates include the converted latitude and longitude coordinates and elevation; convert the geodetic coordinates of each point in the set area of ​​the lunar surface into a three-dimensional position vector of the lunar fixed coordinates; parallelly calculate the altitude angles of other points in each azimuth direction relative to the points where the threads are located within the terrain search range centered on the points where each thread is located according to the thread index of each thread in the parallel thread resource; perform reduction calculations in parallel according to the azimuth angles of each point in the set area of ​​the lunar surface, and select the maximum altitude angle of the azimuth angles of each point in the set area of ​​the lunar surface.

[0113] The aforementioned parallel lunar illumination condition analysis device can execute the parallel lunar illumination condition analysis method provided in any embodiment of the present invention, and possesses the corresponding functional modules and beneficial effects of the execution method. For technical details not fully described in this embodiment, please refer to the parallel lunar illumination condition analysis method provided in any embodiment of the present invention.

[0114] Since the aforementioned parallel lunar illumination condition analysis device is capable of executing the parallel lunar illumination condition analysis method according to the embodiments of the present invention, those skilled in the art will be able to understand the specific implementation and various variations of the parallel lunar illumination condition analysis device according to the embodiments of the present invention based on the parallel lunar illumination condition analysis method according to the embodiments of the present invention. Therefore, how the parallel lunar illumination condition analysis device implements the parallel lunar illumination condition analysis method according to the embodiments of the present invention will not be described in detail here. As long as those skilled in the art can implement the device used in the parallel lunar illumination condition analysis method according to the embodiments of the present invention, it falls within the scope of protection of this application.

[0115] Figure 4 A schematic diagram of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital assistants, cellular phones, smartphones, wearable devices (such as helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0116] like Figure 4 As shown, electronic device 10 includes at least one processor 11 and memory, such as read-only memory (ROM) 12 and random access memory (RAM) 13, communicatively connected to at least one processor 11. The memory stores computer programs executable by the at least one processor. Processor 11 can perform various appropriate actions and processes based on the computer programs stored in ROM 12 or loaded from storage unit 18 into RAM 13. RAM 13 can also store various programs and data required for the operation of electronic device 10. Processor 11, ROM 12, and RAM 13 are interconnected via bus 14. An input / output (I / O) interface 15 is also connected to bus 14.

[0117] Multiple components in the electronic device 10 are connected to the I / O interface 15, including an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0118] Processor 11 can be any general-purpose and / or specialized processing component with processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any other suitable processor, controller, microcontroller, etc. Processor 11 executes the various methods and processes described above, such as the parallel analysis method for lunar illumination conditions.

[0119] Optionally, the parallel analysis method for lunar surface illumination conditions may include: configuring lunar surface illumination calculation associated parameters for a set area on the lunar surface; wherein the lunar surface illumination calculation associated parameters include lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters; allocating parallel computing resources to the lunar surface illumination calculation associated parameters for the set area on the lunar surface; wherein the parallel computing resources include parallel thread resources and memory resources; converting the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar surface illumination calculation associated parameters to obtain converted longitude and latitude coordinates; and parallelly calculating the illumination calculation results of each point in the set area on the lunar surface according to the parallel computing resources, the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates.

[0120] In some embodiments, the parallel lunar illumination condition analysis method can be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the parallel lunar illumination condition analysis method described above can be performed. Alternatively, in other embodiments, processor 11 can be configured to execute the parallel lunar illumination condition analysis method via any other appropriate means (e.g., via firmware).

[0121] Various embodiments of the systems and techniques described above can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system that includes at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0122] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program may be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0123] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, device, or apparatus. A computer-readable storage medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or apparatus, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media may include an electrical connection based on one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0124] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device that has: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0125] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0126] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0127] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in this disclosure can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solutions disclosed in this disclosure can be achieved. This is not a limitation herein.

[0128] The above specific embodiments do not constitute a limitation on the scope of protection of this disclosure. Those skilled in the art will appreciate that various modifications, combinations, sub-combinations, and substitutions may be made based on design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this disclosure shall be included within the scope of protection of this disclosure.

Claims

1. A parallel analysis method for lunar illumination conditions, characterized in that: include: Configuring lunar surface illumination calculation associated parameters for a set area on the lunar surface; wherein the lunar surface illumination calculation associated parameters include lunar surface illumination calculation condition parameters and lunar surface terrain digital elevation parameters; Allocating parallel computing resources to the lunar illumination calculation-related parameters of the set area on the lunar surface; wherein the parallel computing resources include parallel thread resources and memory resources; Convert the pixel coordinates and the longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation associated parameters to obtain converted longitude and latitude coordinates; The illumination calculation results of each point in the set area on the lunar surface are calculated in parallel according to the parallel computing resources, the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates.

2. The method according to claim 1, characterized in that The configuration of the lunar surface illumination calculation associated parameters for the set area on the lunar surface includes: Configuring the lunar surface longitude range, lunar surface latitude range, terrain search range, time span, longitude resolution, latitude resolution and time resolution of the set area on the lunar surface as the lunar surface illumination calculation condition parameters; For the equidistant cylindrical projection mode, configuring the first lunar surface topographic digital elevation parameters of the set area on the lunar surface; wherein the first lunar surface topographic digital elevation parameters include a first resolution, a scale factor, an offset of the projection system origin relative to the target point, the origin longitude, the origin latitude, a first number of rows and columns; and / or For the polar azimuth stereographic projection mode, second lunar surface topographic digital elevation parameters of the set area on the lunar surface are configured; wherein the second lunar surface topographic digital elevation parameters include a second resolution, a scale factor, a second number of rows and a pole mark.

3. The method according to claim 1, characterized in that The step of allocating parallel computing resources to the parameters associated with the lunar illumination calculation for the set area on the lunar surface includes: Determining resource allocation reference parameters according to the target projection method adopted by the set area on the lunar surface; Using a two-dimensional thread block, allocating the memory resources for the lunar surface topography digital elevation parameters in a processor video memory according to the resource allocation reference parameters; The parallel thread resources are allocated to the set area on the lunar surface according to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block.

4. The method according to claim 3, characterized in that The step of converting the pixel coordinates and the longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation associated parameters to obtain converted longitude and latitude coordinates includes: Determining a target projection method to be used for the set area on the lunar surface; According to the correspondence between the pixel coordinates of the set area on the lunar surface and the thread index of the two-dimensional thread block, and the digital elevation parameters of the lunar topography corresponding to the target projection method, the pixel coordinates and the longitude and latitude coordinates of the set area on the lunar surface are converted to obtain the converted longitude and latitude coordinates.

5. The method according to claim 1, characterized in that The step of calculating the illumination calculation results of each point in the set area on the lunar surface according to the parallel computing resources, the lunar illumination calculation associated parameters and the converted longitude and latitude coordinates includes: Calculate the maximum altitude angle between points in the set area on the lunar surface according to the lunar illumination calculation associated parameter and the converted longitude and latitude coordinates; According to the time resolution in the lunar illumination calculation associated parameter, the parallel thread resources are divided into blocks at different times; The threads at each moment obtain the solar altitude angle and the solar azimuth angle in parallel, and according to the maximum altitude angle between the points in the set area on the lunar surface, interpolate and calculate the maximum altitude angle corresponding to the solar azimuth angle in parallel; The threads at each moment compare in parallel the solar altitude angle and the maximum altitude angle corresponding to the solar altitude angle to obtain the illumination calculation results of each point in the set area on the lunar surface.

6. The method according to claim 5, characterized in that The step of calculating the maximum altitude angle between points in the set area on the lunar surface according to the lunar illumination calculation associated parameter and the converted longitude and latitude coordinates includes: According to the time resolution in the moon illumination calculation associated parameter, with each point where each thread in the parallel thread resource is located as the center, search the distance of the terrain search range outward according to the azimuth calculation interval and the distance interval, and interpolate and calculate the corresponding geodetic coordinates of each point at equal distances; wherein the geodetic coordinates include the converted latitude and longitude coordinates and elevation; Converting the geodetic coordinates of each point in the set area on the lunar surface into a three-dimensional position vector of the lunar solid system coordinates; Parallel calculation of the ground elevation angles of other points in various azimuth directions relative to the points where the threads are located within a terrain search range centered on the points where the threads are located according to the thread index of each thread in the parallel thread resource; The reduction calculation is performed in parallel according to the azimuth angles of each point in the set area on the lunar surface, and the maximum altitude angle of the azimuth angles of each point in the set area on the lunar surface is selected.

7. A parallel analysis device for lunar illumination conditions, characterized in that: include: A lunar illumination calculation-related parameter configuration module is used to configure lunar illumination calculation-related parameters of a set area on the lunar surface; wherein the lunar illumination calculation-related parameters include lunar illumination calculation condition parameters and lunar terrain digital elevation parameters; A parallel computing resource allocation module, used to allocate parallel computing resources to the lunar surface illumination calculation associated parameters of the set area on the lunar surface; wherein the parallel computing resources include parallel thread resources and memory resources; A pixel longitude and latitude coordinate conversion module, used to convert the pixel coordinates and longitude and latitude coordinates of the set area on the lunar surface according to the parallel computing resources and the lunar illumination calculation associated parameters to obtain converted longitude and latitude coordinates; The illumination calculation result calculation module is used to calculate the illumination calculation result of each point in the set area on the lunar surface according to the parallel computing resources, the lunar surface illumination calculation associated parameters and the converted longitude and latitude coordinates.

8. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the parallel analysis method of lunar illumination conditions as described in any one of claims 1-6.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the parallel analysis method for lunar illumination conditions described in any one of claims 1-6 when executed.

10. A computer program product comprising a computer program / instructions, wherein: When the computer program / instructions are executed by a processor, the parallel analysis method for lunar illumination conditions described in any one of claims 1-6 is implemented.

Images