Method for on-orbit real-time judgment of sunny area and shadow area by remote sensing satellite

Abstract

The embodiment of the invention discloses a method for on-orbit real-time judgment of a sunny area and a shadow area by a remote sensing satellite, and relates to the technical field of remote sensing satellites, and the method comprises the steps: obtaining a first position vector of the remote sensing satellite in an earth-fixed system, and obtaining a sunlight direction vector of the earth-fixed system; acquiring a second position vector of the projection point of the remote sensing satellite in the earth-fixed system; constructing a judgment vector based on the model of the first position vector relative to the second position vector, the reverse vector of the sunlight direction vector of the earth-fixed system and the position coordinate of the starting point; the height of the tail end point of each judgment vector is calculated through continuous iteration, and whether the remote sensing satellite is located in a shadow area or a sunshine area is judged according to the comparison result of the height and a preset threshold value. According to the embodiment of the invention, only few iterations are needed, the robustness of calculation can be ensured through iterative calculation, the area where the remote sensing satellite is located can be stably and efficiently determined, real-time on-orbit judgment can be realized, and reliable data support can be provided for pose adjustment of the satellite.

Description

Technical Field

[0001] The present application relates to the field of remote sensing satellite technology, and in particular to a method for real-time determination of sunlight areas and shadow areas by remote sensing satellites in orbit. Background Art

[0002] For high-precision optical remote sensing satellites, the workload is a high-resolution camera, and its optical and electronic components have certain requirements for operating temperature. Therefore, the temperature control of the camera is the key to ensure the normal operation of the camera.

[0003] During the satellite's in-orbit flight phase, the satellite's exposure to sunlight changes due to changes in the geometric position relationship between the satellite, the earth and the sun. When the satellite maintains a heliocentric attitude, the camera's optical axis generally points in the opposite direction of the sun. When the satellite is in the sunlit area, the satellite's heliocentric attitude will cause the camera's optical axis to point to the earth; and when the satellite is in the shadow area, the satellite's heliocentric attitude will cause the camera's optical axis to point to deep space. It should be noted that when the camera's optical axis points to deep space, it will accelerate the dissipation of heat, and the camera's temperature control system cannot control the camera's temperature within the operating temperature range. Therefore, when the satellite is in the shadow area, the satellite needs to adjust its heliocentric attitude to an earth-centric attitude, pointing the camera's optical axis toward the earth to reduce heat loss.

[0004] Therefore, there is a need for a method that can efficiently and simply determine in real time on orbit whether a remote sensing satellite is located in a sunny area or a shadow area, so as to ensure that the satellite can adjust its attitude in a timely manner so that related components can work normally. Summary of the invention

[0005] The present application provides a method for determining the sunshine area and the shadow area of ​​a remote sensing satellite in real time on orbit, so as to solve the defects in the above-mentioned related technologies. The technical solution is as follows:

[0006] In a first aspect, an embodiment of the present application provides a method for a remote sensing satellite to determine a sunlit area and a shadowed area in real time on orbit, comprising:

[0007] Obtaining the first position vector of the remote sensing satellite in the ground-fixed system and obtaining the sunlight direction vector of the ground-fixed system;

[0008] Taking the position of the remote sensing satellite as a starting point, projecting the starting point onto the ground to obtain a projection point, and obtaining a second position vector of the projection point on a ground fixed system;

[0009] Constructing a judgment vector: obtaining the modulus of the first position vector corresponding to the starting point relative to the second position vector, obtaining the reverse vector of the sunlight direction vector of the ground-fixed system, and constructing a judgment vector based on the position coordinates of the starting point, the modulus and the reverse vector;

[0010] Based on the position vector of the end point of the judgment vector at the ground fixed system, a new first position vector is updated; based on the new first position vector, a height of the end point at the ground fixed system is calculated; and the height is compared with a preset threshold;

[0011] If the height is less than the preset threshold, it is determined that the remote sensing satellite is located in a shadow area;

[0012] If the height is greater than or equal to the preset threshold, the number of projections is accumulated. If the number of projections is less than the threshold, the end point is updated to the next starting point, and the process jumps to the step of projecting the starting point to the ground to obtain the projection point. If the number of projections is greater than or equal to the threshold, it is determined that the remote sensing satellite is located in a sunny area.

[0013] In an optional solution of the first aspect, before obtaining the first position vector of the remote sensing satellite in the ground-fixed system and obtaining the sunlight direction vector of the ground-fixed system, the method further includes:

[0014] Obtain the current position information of the remote sensing satellite and the direction information of sunlight;

[0015] Acquire a position vector of the remote sensing satellite in an inertial system based on the position information, and acquire a sunlight direction vector of the inertial system based on the sunlight direction information;

[0016] The position vector of the remote sensing satellite in the inertial system and the sunlight direction vector in the inertial system are converted to the earth-fixed system, and the first position vector and the earth-fixed system sunlight direction vector are output.

[0017] In an optional solution of the first aspect, converting the position vector of the remote sensing satellite in the inertial system and the sunlight direction vector of the inertial system to an earth-fixed system to obtain the first position vector and the sunlight direction vector of the earth-fixed system includes:

[0018] Obtaining a coordinate transformation matrix from the inertial system to the earth-fixed system;

[0019] Based on the coordinate transformation matrix, the position vector of the remote sensing satellite in the inertial system is transformed into the earth-fixed system to obtain the first position vector, and the formula is applied:

[0020] R sat_ECF =M ECI2ECF ·R sat_ECI ;

[0021] Based on the coordinate conversion matrix, the inertial system sunlight direction vector is converted to the ground-fixed system to obtain the ground-fixed system sunlight direction vector, and the formula is applied:

[0022] r sun_ECF =M ECI2ECF ·rsun_ECI ;

[0023] Among them, ECI is the inertial system, ECF is the earth-fixed system, R sat_ECI is the position vector of the remote sensing satellite in the inertial system, M ECI2ECF is the coordinate transformation matrix, R sat_ECF is the first position vector, r sun_ECI is the sunlight direction vector of the inertial system, r sun_ECF is the direction vector of the sunlight on the ground.

[0024] In an optional solution of the first aspect, taking the position of the remote sensing satellite as a starting point, projecting the starting point to the ground to obtain a projection point, and obtaining a second position vector of the projection point on a ground fixed system, comprises:

[0025] Obtaining the major semi-axis and the minor semi-axis of the earth reference ellipsoid, converting the rectangular coordinates of the first position vector in the earth-fixed system into longitude and latitude coordinates, and obtaining the longitude and latitude coordinates of the starting point;

[0026] The initial value of the altitude of the latitude and longitude coordinates of the starting point is:

[0027]

[0028] The initial latitude value of the latitude and longitude coordinates of the starting point is:

[0029]

[0030] Calculate variable N1:

[0031]

[0032] Calculate the altitude, latitude and longitude of the starting point based on the variable N1 and the initial altitude value:

[0033]

[0034] L1=arctan(Y sat / X sat );

[0035] Output the latitude and longitude coordinates of the starting point:

[0036]

[0037] The longitude and latitude coordinates of the projection point are determined based on the longitude and latitude coordinates of the starting point, and the following is obtained:

[0038]

[0039] Calculate the variable N:

[0040]

[0041] Based on the variable N, the longitude and latitude coordinates of the projection point are converted into rectangular coordinates in the ground fixed system to obtain the second position vector, and the formula is applied:

[0042]

[0043] Among them, a e is the semi-major axis of the Earth reference ellipsoid, b e is the minor semi-axis of the earth reference ellipsoid, e is the eccentricity of the earth reference ellipsoid, H0 is the initial value of height, B0 is the initial value of latitude, H1 is the height, B1 is the latitude, L1 is the longitude, R satn_ECF is the second position vector.

[0044] In an optional solution of the first aspect, after calculating the altitude, latitude and longitude of the starting point, the method further includes:

[0045] Calculate the latitude difference between the latitude and the initial latitude value:

[0046] dB = |B1-B0|;

[0047] Calculate the height difference between the height and the initial height value:

[0048] dH=|H1-H0|;

[0049] Comparing the latitude difference with a preset latitude difference threshold, and comparing the height difference with a preset height difference threshold;

[0050] If the latitude difference and the altitude difference are both less than the corresponding thresholds, then executing the step of outputting the latitude and longitude coordinates of the starting point;

[0051] Otherwise, the step of calculating the altitude, latitude and longitude of the starting point is performed, wherein dB is the latitude difference and dH is the altitude difference.

[0052] In a second aspect, an embodiment of the present application further provides a device for determining a sunlit area and a shadowed area in real time on a remote sensing satellite, comprising:

[0053] A data input module, used to obtain a first position vector of the remote sensing satellite in the ground-fixed system, and also used to obtain a sunlight direction vector of the ground-fixed system;

[0054] A projection module is used to take the position of the remote sensing satellite as a starting point, project the starting point to the ground to obtain a projection point, and obtain a second position vector of the projection point on the ground fixed system;

[0055] A judgment module, used to construct a judgment vector: obtain the modulus of the second position vector corresponding to the starting point relative to the third position vector, obtain a reverse vector to the second sunlight direction vector, and construct a judgment vector based on the position coordinates of the starting point, the modulus and the reverse vector;

[0056] The judgment module is further used to update the second position vector based on the end point of the judgment vector, calculate the height of the end point on the ground fixed system based on the second position vector, and compare the height with a preset threshold;

[0057] If the height is less than the preset threshold, it is determined that the remote sensing satellite is located in a shadow area;

[0058] If the height is greater than or equal to the preset threshold, the number of projections is accumulated. If the number of projections is less than the threshold, the end point is updated to the next starting point, and the process jumps to the step of projecting the starting point to the ground to obtain the projection point. If the number of projections is greater than or equal to the threshold, it is determined that the remote sensing satellite is located in a sunny area.

[0059] In a third aspect, an embodiment of the present application further provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method provided in the first aspect of the embodiment of the present application or any one of the implementation methods of the first aspect is implemented.

[0060] In a fourth aspect, the present application also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method provided by the first aspect of the embodiment of the present application or any one of the implementations of the first aspect.

[0061] The beneficial effects brought about by the technical solutions provided by some embodiments of the present application include at least:

[0062] The embodiment of the present application provides a method for real-time judgment of the sunlit area and the shadowed area by a remote sensing satellite in orbit. The method only needs to input the orbital information of the remote sensing satellite at the current moment, and does not require other information uploaded by ground assistance. Multiple judgment vectors are calculated cyclically, and the height of the end point is determined each time based on the latitude and longitude coordinates of the end point of the judgment vector, so as to judge whether there is an intersection between the earth and a vector starting from the position coordinates of the remote sensing satellite and in the opposite direction of the vector of the sunlight direction of the earth-fixed system. If there is an intersection, it can be determined that the remote sensing satellite is located in the shadow area, otherwise it is located in the sunlit area. The method provided in the embodiment of the present application only requires a small number of iterations. The robustness of the calculation can be ensured by iteratively calculating multiple judgment vectors, and the area where the remote sensing satellite is located can be determined stably and efficiently, and real-time on-orbit judgment can be realized, so as to provide reliable data support for the satellite's posture adjustment based on the calculation results. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] In order to more clearly illustrate the technical solutions in the present application or related technologies, the drawings required for use in the embodiments or related technical descriptions are briefly introduced below. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0064] Figure 1 It is a flow chart of a method for determining a sunlit area and a shadowed area in real time on a remote sensing satellite in orbit provided by an embodiment of the present application;

[0065] Figure 2 It is a schematic diagram of the relative position relationship of a method for a remote sensing satellite to judge the sunshine area and the shadow area in real time on orbit provided by an embodiment of the present application;

[0066] Figure 3 This is a schematic diagram of the structure of a device for real-time determination of a sunlit area and a shadowed area of ​​a remote sensing satellite in orbit provided by an embodiment of the present application;

[0067] Figure 4 It is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0068] In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in this application will be clearly and completely described below in conjunction with the drawings in this application. Obviously, the described embodiments are part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0069] The terms "including" and "having" and any variations thereof in the specification and claims of the present application and the above drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device including a series of steps or modules is not limited to the listed steps or modules, but may optionally include steps or modules that are not listed, or may optionally include other steps or modules that are inherent to these processes, methods, products or devices.

[0070] It should be noted that the terms "first\second" involved in the present application are only used to distinguish similar objects, and do not represent a specific order for the objects. It is understandable that "first\second" can be interchanged with a specific order or sequence where permitted. It should be understood that the objects distinguished by "first\second" can be interchanged where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those described or illustrated herein.

[0071] The present application is described in detail below with reference to specific embodiments.

[0072] Next, combine Figure 1 , introduces a method for a remote sensing satellite to determine the sunlit area and shadowed area in real time on orbit provided by an embodiment of the present application. For details, please refer to Figure 2 , Figure 2 The flowchart of a method for determining the sunshine area and shadow area of ​​a remote sensing satellite in real time on orbit provided by an embodiment of the present application is shown. Figure 1 As shown, the method comprises the following steps:

[0073] S101, obtaining a first position vector of a remote sensing satellite in a ground-fixed system, and obtaining a sunlight direction vector of the ground-fixed system.

[0074] Taking the position of the remote sensing satellite as the starting point, perform the following steps:

[0075] S102, projecting the starting point onto the ground to obtain a projection point, and obtaining a second position vector of the projection point fixed to the ground.

[0076] S103, constructing a judgment vector: obtaining the modulus of the first position vector corresponding to the starting point relative to the second position vector, obtaining the reverse vector of the ground-fixed system sunlight direction vector, and constructing a judgment vector based on the position coordinates of the starting point, the modulus and the reverse vector.

[0077] S104: A new first position vector is obtained based on the position vector of the end point of the judgment vector on the ground fixed system, and a height of the end point on the ground fixed system is calculated based on the new first position vector.

[0078] S105: Compare the height with a preset threshold.

[0079] Further, if the calculated height of the end point in the ground fixed system is less than the preset threshold, it indicates that a vector starting from the position coordinates of the remote sensing satellite and in the opposite direction of the sunlight direction vector of the ground fixed system has an intersection with the earth, and the following steps are performed:

[0080] S106, determining that the remote sensing satellite is located in the shadow area at the current moment.

[0081] Furthermore, if the calculated height of the end point on the ground fixed system is greater than or equal to the preset threshold, the following steps are further performed:

[0082] S107, accumulating the number of projections and comparing the number of projections with a preset number threshold. If the number of projections is less than the number threshold, executing the following steps:

[0083] S108, updating the end point to the next start point, jumping to S102, and executing S102 and subsequent steps;

[0084] If the projection times are greater than or equal to the times threshold, perform the following steps:

[0085] S109, determining whether the remote sensing satellite is located in a sunny area at the current moment.

[0086] In some embodiments, before S101, the following steps are also included:

[0087] Obtain the current position information of the remote sensing satellite and the direction information of sunlight;

[0088] Acquire a position vector of the remote sensing satellite in an inertial system based on the position information, and acquire a sunlight direction vector of the inertial system based on the sunlight direction information;

[0089] The position vector of the remote sensing satellite in the inertial system and the sunlight direction vector in the inertial system are converted to the earth-fixed system, and the first position vector and the earth-fixed system sunlight direction vector are output.

[0090] It can be understood that the position information of the remote sensing satellite at the corresponding time can be determined by the satellite orbit information obtained in real time. Specifically, the vector of the inertial system (Earth-Centered Inertial Frame, ECI) can be converted to the earth-fixed system (Earth-Centered Fixed Frame, ECF) through a preset coordinate conversion matrix, which specifically includes:

[0091] Obtaining a coordinate transformation matrix from the inertial system to the earth-fixed system;

[0092] Based on the coordinate transformation matrix, the position vector of the remote sensing satellite in the inertial system is transformed into the earth-fixed system to obtain the first position vector, and the formula is applied:

[0093] R sat_ECF =M ECI2ECF ·R sat_ECI ;

[0094] Based on the coordinate conversion matrix, the inertial system sunlight direction vector is converted to the ground-fixed system to obtain the ground-fixed system sunlight direction vector, and the formula is applied:

[0095] r sun_ECF =M ECI2ECF ·r sun_ECI ;

[0096] Among them, ECI is the inertial system, ECF is the earth-fixed system, R sat_ECI is the position vector of the remote sensing satellite in the inertial system, M ECI2ECF is the coordinate transformation matrix, R sat_ECF is the first position vector, r sun_ECI is the sunlight direction vector of the inertial system, r sun_ECF is the direction vector of the sunlight on the ground.

[0097] In some embodiments, in S102, firstly, a position vector of a sub-satellite point of the remote sensing satellite on a ground fixed system, that is, a second position vector, is determined;

[0098] Among them, the sub-satellite point of the remote sensing satellite can be understood as taking the position of the remote sensing satellite as the starting point, generating a line connecting the starting point and the center of the earth, and the intersection of the line and the earth's surface is the sub-satellite point, that is, the projection point of the remote sensing satellite's position on the ground.

[0099] S102 specifically includes:

[0100] S1021, obtaining the major semi-axis and the minor semi-axis of the earth reference ellipsoid, converting the rectangular coordinates of the first position vector in the earth-fixed system into longitude and latitude coordinates, and obtaining the longitude and latitude coordinates of the starting point.

[0101] The rectangular coordinates of the first position vector of the initial point in the ground-fixed system can be expressed as:

[0102]

[0103] The initial value of the altitude of the latitude and longitude coordinates of the starting point is:

[0104]

[0105] The initial latitude value of the latitude and longitude coordinates of the starting point is:

[0106]

[0107] S1022, calculating variable N1, calculating the altitude, latitude and longitude of the starting point based on variable N1 and the initial altitude value, and applying the formula:

[0108]

[0109] L1=arctan(Y sat / X sat ).

[0110] S1023, output the latitude and longitude coordinates of the starting point:

[0111]

[0112] S1024, determining the latitude and longitude coordinates of the projection point based on the latitude and longitude coordinates of the starting point. Since the projection point is located on the ground, the height H needs to be set to 0.

[0113] Based on this, the altitude value of the starting point is set to 0, and the longitude and latitude coordinates of the projection point LLA can be obtained. satn_ECF :

[0114]

[0115] S1025, calculating a variable N, and based on the variable N, converting the longitude and latitude coordinates of the projection point into rectangular coordinates in the ground fixed system to obtain the second position vector, and applying the formula:

[0116]

[0117] For any geographic longitude and latitude coordinates The formula can be applied:

[0118]

[0119] Since the height H of the projection point is set to 0, the second position vector obtained is as follows:

[0120]

[0121] Among them, a e is the semi-major axis of the Earth reference ellipsoid, b e is the minor semi-axis of the earth reference ellipsoid, e is the eccentricity of the earth reference ellipsoid, H0 is the initial value of height, B0 is the initial value of latitude, H1 is the height, B1 is the latitude, L1 is the longitude, R satn_ECF is the second position vector.

[0122] In some embodiments, after calculating the altitude, latitude and longitude of the starting point in S1022, the latitude difference dB between the latitude and the initial latitude value is further calculated, and the height difference dH between the altitude and the initial altitude value is calculated, using the following formula:

[0123] dB = |B1-B0|;

[0124] dH=|H1-H0|;

[0125] Comparing the latitude difference with a preset latitude difference threshold, and comparing the height difference with a preset height difference threshold;

[0126] If the latitude difference and the altitude difference are both less than the corresponding thresholds, executing step S1023 to output the longitude and latitude coordinates of the starting point;

[0127] Otherwise, re-execute step S1022.

[0128] In one embodiment, next combine Figure 2 , introduces the specific calculation process of a method for real-time determination of sunlit areas and shadowed areas by a remote sensing satellite in orbit provided by an embodiment of the present application, such as Figure 2 The relative position diagram of the satellite, the earth and the sunlight shown in the figure, S is the position of the satellite, T0 is the position of the satellite sub-satellite point, specifically includes the following steps:

[0129] S201, obtain the position vector R of the remote sensing satellite in the inertial system sat_ECI , get the sunlight direction vector r in the inertial system sun_ECI .

[0130] S202, applying the coordinate transformation matrix from the inertial system to the earth-fixed system, transforming the vector of the inertial system to the earth-fixed system:

[0131] R sat_ECF =M ECI2ECF ·R sat_ECI ;

[0132] r sun_EcF =M ECI2ECF ·r sun_ECI ;

[0133] Output the first position vector R of the remote sensing satellite on the ground sat_ECF and the ground-fixed sunlight direction vector r sun_ECF .

[0134] S203, obtain the major semi-axis a of the earth reference ellipsoid e and the minor axis b e , the first position vector R under the ground fixed system sat_ECFThe rectangular coordinates are converted into longitude and latitude coordinates, and the position S of the remote sensing satellite is taken as the starting point to obtain the longitude and latitude coordinates of the starting point. For details, please refer to the description of S1021-S1023, which will not be repeated here.

[0135] Furthermore, the projection point of the starting point on the ground is determined based on the position vector of the starting point. When the position S of the remote sensing satellite is the starting point, the obtained projection point, that is, the sub-satellite point is Figure 2 The second position vector R of the point T0 on the ground fixed system can be determined based on S1024-S1025. T0_ECF ,Right now:

[0136] R T0_ECF =R satn_ECF ;

[0137] Among them, R satn_ECF That is, the position vector of the satellite's sub-satellite point on the ground fixed system.

[0138] It can be understood that the first position vector of the ground-fixed system can be regarded as Figure 2 The vector shown The intersection of the first position vector and the earth's surface is T0, and the second position vector is

[0139] S204, constructing a judgment vector, specifically including:

[0140] First determine the first position vector R sat_ECF Relative to the second position vector R T0_ECF Model R ST , apply the formula:

[0141] R ST =|R T0_ECF -R sat_ECF |;

[0142] Furthermore, taking the satellite position in the ground-fixed system as the starting point, the direction is selected as -r sun_ECF , the construction model is R ST Vector The details are as follows:

[0143]

[0144] Among them, T1 is the judgment vector The end point of T1 is located on a vector that takes the position coordinates of the remote sensing satellite as its starting point and is in the opposite direction of the sunlight direction vector of the earth-fixed system as its direction.

[0145] S205, based on the judgment vector The end point T1 is updated to obtain the new first position vector Application formula:

[0146]

[0147] Based on the new first position vector Execute steps S1021-S1023 to determine the longitude and latitude coordinates of point T1

[0148]

[0149] Then according to the longitude and latitude coordinates Get the height of point T1.

[0150] S206, compare the height of point T1 with a preset threshold, for example, the height threshold may be 10 -6 .

[0151] Specifically, if the height of point T1 is less than a preset threshold, it indicates that a vector starting from the position coordinates of the remote sensing satellite and in the opposite direction of the direction of the sunlight vector of the earth-fixed system intersects the earth, and the following steps are performed:

[0152] S207, determining that the remote sensing satellite is located in the shadow area at the current moment.

[0153] Furthermore, if the height of point T1 is greater than or equal to the preset threshold, the following steps are further performed:

[0154] S208, cumulative projection times, let i=i+1

[0155] It should be noted that, before generating the satellite sub-satellite point T0 of the satellite position S in S203, the initial projection number i=0 may be set.

[0156] The number of projections is compared with the preset number threshold. If the number of projections is less than the number threshold, for example, the number of projections i=1, the number threshold can be set to 20, and i is less than 20, then the following steps are performed:

[0157] S209, update the end point T1 to the next starting point, jump to S203, and execute S203 and subsequent steps.

[0158] It can be understood that after T1 is updated as the next starting point, the projection point T′1 of T1 on the ground is further generated based on the first position vector corresponding to T1, and the second position vector obtained in step S203 is updated, that is, the rectangular coordinate of the position vector of T1′ in the ground fixed system is

[0159] Based on the rectangular coordinates of the position vector and A new judgment vector is constructed, and then it is iterated in a loop to continuously calculate the virtual point T on the vector with the position coordinates of the remote sensing satellite as the starting point and the direction opposite to the direction of the sunlight vector of the fixed system on the ground i The longitude and latitude coordinates of the virtual point T i The coordinates of the remote sensing satellite are compared with the preset threshold to determine whether the remote sensing satellite is in the sunlit area or the shadow area.

[0160] It should be noted that in S207, as long as the calculated height is less than the preset threshold, it means that a vector starting from the position coordinates of the remote sensing satellite and with a direction in the opposite direction of the direction of the sunlight vector of the earth-fixed system intersects the earth, and the judgment result that the satellite is located in the shadow area can be directly output.

[0161] In S209, if the number of projections is less than the given number threshold, it is necessary to repeat S203 and subsequent steps, and perform iterative calculations in a loop. When the number of projections is greater than or equal to the given number threshold, step S210 is performed:

[0162] S210, determining that the remote sensing satellite is located in a sunny area at the current moment.

[0163] Through continuous iteration in this way, it is possible to quickly and accurately determine whether the satellite is in the sunlit area or the shadowed area, and to achieve real-time on-orbit judgment without the need for additional input information. Only the on-orbit position and time of the satellite are required. It has high robustness and real-time performance. Based on the judgment results, the satellite's position and posture can be adjusted in time to avoid defects such as excessive heat dissipation when the camera's optical axis points to deep space, which affects the working state of the satellite's components, reduce heat consumption, and ensure that the original components can work normally.

[0164] The following are device embodiments of the present application, which can be used to execute the method embodiments of the present application. For details not disclosed in the device embodiments of the present application, please refer to the method embodiments of the present application.

[0165] Next, please refer to Figure 30, which is a schematic diagram of the structure of a device for real-time determination of sunlit areas and shadowed areas of a remote sensing satellite in orbit provided by an exemplary embodiment of the present application. The device can be implemented as all or part of a terminal through software, hardware, or a combination of both, and can also be integrated on a server as an independent module. A device for real-time determination of sunlit areas and shadowed areas of a remote sensing satellite in orbit in an embodiment of the present application can be applied to a terminal or a cloud. The device 30 includes a data input module 301, a projection module 302, and a determination module 303, wherein:

[0166] The data input module 301 is used to obtain the first position vector of the remote sensing satellite in the ground-fixed system, and is also used to obtain the sunlight direction vector of the ground-fixed system;

[0167] The projection module 302 is used to take the position of the remote sensing satellite as a starting point, project the starting point to the ground to obtain a projection point, and obtain a second position vector of the projection point on the ground fixed system;

[0168] The judgment module 303 is used to construct a judgment vector: obtain the modulus of the first position vector corresponding to the starting point relative to the second position vector, obtain the reverse vector of the sunlight direction vector of the ground-fixed system, and construct a judgment vector based on the position coordinates of the starting point, the modulus and the reverse vector;

[0169] The judgment module 303 is further configured to obtain a new first position vector based on the position vector of the end point of the judgment vector at the ground fixed system, calculate the height of the end point at the ground fixed system based on the new first position vector, and compare the height with a preset threshold;

[0170] If the height is less than the preset threshold, it is determined that the remote sensing satellite is located in a shadow area;

[0171] If the height is greater than or equal to the preset threshold, the number of projections is accumulated. If the number of projections is less than the threshold, the end point is updated to the next starting point, and the process jumps to the step of projecting the starting point to the ground to obtain the projection point. If the number of projections is greater than or equal to the threshold, it is determined that the remote sensing satellite is located in a sunny area.

[0172] It should be noted that when the device 30 provided in the above embodiment executes a method for real-time determination of a sunlit area and a shadowed area by a remote sensing satellite in orbit, only the division of the above functional modules is used as an example. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the device provided in the above embodiment and the embodiment of the method for real-time determination of a sunlit area and a shadowed area by a remote sensing satellite in orbit belong to the same concept. The implementation process thereof is detailed in the method embodiment and will not be repeated here.

[0173] An embodiment of the present application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of any of the above-mentioned method embodiments when executing the program.

[0174] See also Figure 4 , which is a structural block diagram of an electronic device provided in an embodiment of the present application.

[0175] like Figure 4 As shown, the electronic device 400 includes: a processor 401 and a memory 402 .

[0176] In the embodiment of the present application, the processor 401 is the control center of the computer system, which can be a processor of a physical machine or a processor of a virtual machine. The processor 401 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 401 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array).

[0177] The processor 401 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also called a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state.

[0178] The memory 402 may include one or more computer-readable storage media, which may be non-transitory. The memory 402 may also include a high-speed random access memory, and a non-volatile memory, such as one or more disk storage devices, flash memory storage devices. In some embodiments of the present application, the non-transitory computer-readable storage medium in the memory 402 is used to store at least one instruction, which is used to be executed by the processor 401 to implement the method in the embodiment of the present application.

[0179] In some embodiments, the electronic device 400 further includes: a peripheral device interface 403 and at least one peripheral device 404. The processor 401, the memory 402 and the peripheral device interface 403 can be connected via a bus or a signal line. Each peripheral device 404 can be connected to the peripheral device interface 403 via a bus, a signal line or a circuit board. Specifically, the peripheral device 404 includes: a display screen, a camera and an audio circuit. The peripheral device interface 403 can be used to connect at least one peripheral device related to I / O (Input / Output) to the processor 401 and the memory 402.

[0180] In some embodiments of the present application, the processor 401, the memory 402, and the peripheral device interface 403 are integrated on the same chip or circuit board; in some other embodiments of the present application, any one or two of the processor 401, the memory 402, and the peripheral device interface 403 can be implemented on a separate chip or circuit board. This embodiment of the present application does not specifically limit this.

[0181] The electronic device structure block diagram shown in the embodiment of the present application does not constitute a limitation on the electronic device 400. The electronic device 400 may include more or fewer components than shown in the figure, or combine certain components, or adopt a different component arrangement.

[0182] The embodiment of the present application also provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the method of any of the above embodiments are implemented. The computer-readable storage medium may include, but is not limited to, any type of disk, including a floppy disk, an optical disk, a DVD, a CD-ROM, a micro drive, and a magneto-optical disk, a ROM, a RAM, an EPROM, an EEPROM, a DRAM, a VRAM, a flash memory device, a magnetic card or an optical card, a nanosystem (including a molecular memory IC), or any type of medium or device suitable for storing instructions and / or data.

[0183] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, by hardware. Based on this understanding, the above technical solution can be essentially or in other words, the part that contributes to the relevant technology can be embodied in the form of a software product, which can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiment.

[0184] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit it. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for determining the sunshine area and shadow area in real time by a remote sensing satellite on orbit, characterized in that: include: Obtaining the first position vector of the remote sensing satellite in the ground-fixed system and obtaining the sunlight direction vector of the ground-fixed system; Taking the position of the remote sensing satellite as a starting point, projecting the starting point onto the ground to obtain a projection point, and obtaining a second position vector of the projection point on a ground fixed system; Constructing a judgment vector: obtaining the modulus of the first position vector corresponding to the starting point relative to the second position vector, obtaining the reverse vector of the sunlight direction vector of the ground-fixed system, and constructing a judgment vector based on the position coordinates of the starting point, the modulus and the reverse vector; Based on the position vector of the end point of the judgment vector at the ground fixed system, a new first position vector is updated; based on the new first position vector, a height of the end point at the ground fixed system is calculated; and the height is compared with a preset threshold; If the height is less than the preset threshold, it is determined that the remote sensing satellite is located in a shadow area; If the height is greater than or equal to the preset threshold, the number of projections is accumulated. If the number of projections is less than the threshold, the end point is updated to the next starting point, and the process jumps to the step of projecting the starting point to the ground to obtain the projection point. If the number of projections is greater than or equal to the threshold, it is determined that the remote sensing satellite is located in a sunny area.

2. The method for determining the sunshine area and shadow area of ​​a remote sensing satellite in real time on orbit according to claim 1, characterized in that: Before obtaining the first position vector of the remote sensing satellite in the ground-fixed system and obtaining the sunlight direction vector of the ground-fixed system, the method further includes: Obtain the current position information of the remote sensing satellite and the direction information of sunlight; Acquire a position vector of the remote sensing satellite in an inertial system based on the position information, and acquire a sunlight direction vector of the inertial system based on the sunlight direction information; The position vector of the remote sensing satellite in the inertial system and the sunlight direction vector in the inertial system are converted to the earth-fixed system, and the first position vector and the earth-fixed system sunlight direction vector are output.

3. The method for determining the sunshine area and shadow area of ​​a remote sensing satellite in real time on orbit according to claim 2, characterized in that: The converting the position vector of the remote sensing satellite in the inertial system and the sunlight direction vector of the inertial system to the earth-fixed system to obtain the first position vector and the sunlight direction vector of the earth-fixed system includes: Obtaining a coordinate transformation matrix from the inertial system to the earth-fixed system; Based on the coordinate transformation matrix, the position vector of the remote sensing satellite in the inertial system is transformed into the earth-fixed system to obtain the first position vector, and the formula is applied: R sat_ECF =M ECI2ECF ·R sat_ECI ; Based on the coordinate conversion matrix, the inertial system sunlight direction vector is converted to the ground-fixed system to obtain the ground-fixed system sunlight direction vector, and the formula is applied: r sun_ECF ＝M ECI2ECF ·r sun_ECI ； Among them, ECI is the inertial system, ECF is the earth-fixed system, R sat_ECI is the position vector of the remote sensing satellite in the inertial system, M ECI2ECF is the coordinate transformation matrix, R sat_ECF is the first position vector, r sun_ECI is the sunlight direction vector of the inertial system, r sun_ECF is the direction vector of the sunlight on the ground.

4. The method for determining the sunshine area and shadow area of ​​a remote sensing satellite in real time on orbit according to claim 3, characterized in that: The method of taking the position of the remote sensing satellite as a starting point, projecting the starting point onto the ground to obtain a projection point, and obtaining a second position vector of the projection point on a ground fixed system includes: Obtaining the major semi-axis and the minor semi-axis of the earth reference ellipsoid, converting the rectangular coordinates of the first position vector in the earth-fixed system into longitude and latitude coordinates, and obtaining the longitude and latitude coordinates of the starting point; The initial value of the altitude of the latitude and longitude coordinates of the starting point is: The initial latitude value of the latitude and longitude coordinates of the starting point is: Calculate variable N1: Calculate the altitude, latitude and longitude of the starting point based on the variable N1 and the initial altitude value: L1=arctan(Y sat / X sat ); Output the latitude and longitude coordinates of the starting point: The longitude and latitude coordinates of the projection point are determined based on the longitude and latitude coordinates of the starting point, and the following is obtained: Calculate the variable N: Based on the variable N, the longitude and latitude coordinates of the projection point are converted into rectangular coordinates in the ground fixed system to obtain the second position vector, and the formula is applied: Among them, a e is the semi-major axis of the Earth reference ellipsoid, b e is the minor semi-axis of the earth reference ellipsoid, e is the eccentricity of the earth reference ellipsoid, H0 is the initial value of height, B0 is the initial value of latitude, H1 is the height, B1 is the latitude, L1 is the longitude, R satn_ECF is the second position vector.

5. The method for determining the sunshine area and shadow area of ​​a remote sensing satellite in real time on orbit according to claim 4, characterized in that: After calculating the altitude, latitude and longitude of the starting point, the method further includes: Calculate the latitude difference between the latitude and the initial latitude value: dB = |B1-B0|; Calculate the height difference between the height and the initial height value: dH=|H1-H0|; Comparing the latitude difference with a preset latitude difference threshold, and comparing the height difference with a preset height difference threshold; If the latitude difference and the altitude difference are both less than the corresponding thresholds, then executing the step of outputting the latitude and longitude coordinates of the starting point; Otherwise, executing the step of calculating the altitude, latitude and longitude of the starting point; Among them, dB is the latitude difference, and dH is the height difference.

6. A device for determining the sunshine area and shadow area of ​​a remote sensing satellite in real time on orbit, characterized in that: include: A data input module is used to obtain a first position vector of the remote sensing satellite in the ground-fixed system, and is also used to obtain a sunlight direction vector of the ground-fixed system; A projection module is used to take the position of the remote sensing satellite as a starting point, project the starting point to the ground to obtain a projection point, and obtain a second position vector of the projection point on the ground fixed system; A judgment module, used to construct a judgment vector: obtain the modulus of the first position vector corresponding to the starting point relative to the second position vector, obtain the reverse vector of the sunlight direction vector of the ground-fixed system, and construct a judgment vector based on the position coordinates of the starting point, the modulus and the reverse vector; The judgment module is further used to update a new first position vector based on the position vector of the end point of the judgment vector on the ground fixed system, calculate the height of the end point on the ground fixed system based on the new first position vector, and compare the height with a preset threshold; If the height is less than the preset threshold, it is determined that the remote sensing satellite is located in a shadow area; If the height is greater than or equal to the preset threshold, the number of projections is accumulated. If the number of projections is less than the threshold, the end point is updated to the next starting point, and the process jumps to the step of projecting the starting point to the ground to obtain the projection point. If the number of projections is greater than or equal to the threshold, it is determined that the remote sensing satellite is located in a sunny area.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the method according to any one of claims 1 to 5 are implemented.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented.

Images