Satellite Cross-Domain Information Determination Method

Abstract

The present invention provides a method for determining satellite cross-domain information, which can be applied to the field of aerospace technology. The method includes: determining a plurality of target sub-satellite points related to the area to be measured from a plurality of sub-satellite points included in a search engine based on the area position information of the area to be measured; dividing the plurality of target sub-satellite points into at least one sub-satellite point sequence based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points; performing boundary point verification on the sub-satellite point sequence to obtain a verification result; and determining the satellite cross-domain information of the area to be measured based on the verification results of each of the at least one sub-satellite point sequence, thereby at least partially solving the technical problem of low calculation efficiency of satellite cross-domain information in the related art, and realizing the improvement of the speed of determining satellite cross-domain information while ensuring the accuracy of the information.

Description

Technical Field

[0001] The present invention relates to the field of aerospace, and particularly to a method, device, equipment, medium and program product for determining satellite cross-domain information. Background Art

[0002] When a satellite orbits the Earth, each position on its orbit projects a sub-satellite point on the Earth's surface, and the continuous curve formed by connecting all the sub-satellite points in sequence is the satellite trajectory. Satellite cross-domain information determination technology usually analyzes based on the satellite trajectory, and determines when and where the satellite will pass over the target area by calculating the line-plane spatial relationship between each sub-satellite point and the target area respectively.

[0003] In the process of implementing the inventive concept, it is found that there are at least the following problems in the related art: Since a large number of sub-satellite points are included in the satellite trajectory of a satellite, there is a technical problem of low calculation efficiency in determining satellite cross-domain information by calculating the line-plane spatial relationship between each sub-satellite point and the target area respectively. Summary of the Invention

[0004] In view of the above problems, the present invention provides a method, device, equipment, medium and program product for determining satellite cross-domain information.

[0005] According to one aspect of the present invention, there is provided a method for determining satellite cross-domain information, including: determining a plurality of target sub-satellite points related to the area to be measured from a plurality of sub-satellite points included in a search engine based on the area position information of the area to be measured; dividing the plurality of target sub-satellite points into at least one sub-satellite point sequence based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points; performing boundary point verification on the sub-satellite point sequence to obtain a verification result; and determining the satellite cross-domain information of the area to be measured based on the verification results of each of the at least one sub-satellite point sequence.

[0006] Another aspect of the present invention provides a device for determining satellite cross-domain information, including: a first determination module for determining a plurality of target sub-satellite points related to the area to be measured from a plurality of sub-satellite points included in a search engine based on the area position information of the area to be measured; a division module for dividing the plurality of target sub-satellite points into at least one sub-satellite point sequence based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points; a verification module for performing boundary point verification on the sub-satellite point sequence to obtain a verification result; and a second determination module for determining the satellite cross-domain information of the area to be measured based on the verification results of each of the at least one sub-satellite point sequence.

[0007] Another aspect of the present invention provides an electronic device, including: one or more processors; a memory for storing one or more computer programs, wherein the above one or more processors execute the above one or more computer programs to implement the steps of the above method.

[0008] Another aspect of the present invention also provides a computer-readable storage medium, on which a computer program or instruction is stored, and when the computer program or instruction is executed by a processor, the steps of the above method are implemented.

[0009] Another aspect of the present invention also provides a computer program product, including a computer program or instruction, and when the computer program or instruction is executed by a processor, the steps of the above method are implemented.

[0010] According to the satellite cross-domain information determination method of the present invention, by using the spatial indexing technology of a search engine, a plurality of target sub-satellite points related to the area to be measured are quickly screened from multiple sub-satellite points. And by using timestamps and satellites related to the target sub-satellite points, the sub-satellite point sequences are divided to construct at least one cross-domain trajectory. And the boundary points of the sub-satellite point sequences are respectively verified, and based on the verification results, the satellite cross-domain information of the area to be measured is determined. Since after determining a plurality of target sub-satellite points related to the area to be measured by using a search engine, the plurality of target sub-satellite points are divided into at least one sub-satellite point sequence by using timestamps and satellites related to the target sub-satellite points, and more accurate satellite cross-domain information is obtained by performing boundary point verification on the sub-satellite point sequences, thereby realizing that satellite cross-domain information can be obtained by calculating some sub-satellite points. Therefore, at least partially solving the technical problem of low calculation efficiency of satellite cross-domain information in the related art, and realizing that while improving the determination speed of satellite cross-domain information, the accuracy of the information is also guaranteed. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] Through the following description of the embodiments of the present invention with reference to the drawings, the above content and other objects, features and advantages of the present invention will become clearer.

[0012] Figure 1 The application scenario diagram of the satellite cross-domain information determination method, device, equipment, medium and program product according to the embodiments of the present invention is shown.

[0013] Figure 2 The flowchart of the satellite cross-domain information determination method according to the embodiments of the present invention is shown.

[0014] Figure 3 The flowchart of determining the sub-satellite point sequence according to the embodiments of the present invention is shown.

[0015] Figure 4(a) shows the first schematic diagram of the satellite cross-domain information according to the embodiments of the present invention.

[0016] Figure 4(b) shows the second schematic diagram of the satellite cross-domain information according to the embodiments of the present invention.

[0017] Figure 4(c) shows the third schematic diagram of the satellite cross-domain information according to the embodiments of the present invention.

[0018] Figure 4(d) shows a fourth schematic diagram of satellite cross-domain information according to an embodiment of the present invention.

[0019] Figure 5 A structural block diagram of a satellite cross-domain information determination device according to an embodiment of the present invention is shown.

[0020] Figure 6 A block diagram of an electronic device suitable for implementing a satellite cross-domain information determination method according to an embodiment of the present invention is shown. Detailed implementation manners

[0021] Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are merely exemplary and are not intended to limit the scope of the present invention. In the following detailed description, for the sake of explanation, many specific details are set forth in order to provide a comprehensive understanding of the embodiments of the present invention. However, it is obvious that one or more embodiments can be implemented without these specific details. In addition, in the following description, descriptions of well-known structures and technologies are omitted to avoid unnecessarily confusing the concepts of the present invention.

[0022] The terms used herein are merely for describing specific embodiments and are not intended to limit the present invention. The terms "including", "comprising", etc. used herein indicate the presence of the described features, steps, operations, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, or components.

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

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

[0025] During the research, it was found that satellite trajectory data is usually continuous and contains a large number of sub-satellite points, but only some of them are located in the target area. When performing cross-domain time period calculations, if the satellite trajectory is predicted, it is usually necessary to predict the sub-satellite points second by second, and when calculating the sub-satellite points passing through the target area, it is necessary to traverse all sub-satellite points and calculate the line-surface spatial relationship for each sub-satellite point, resulting in redundant calculation steps, a long algorithm time, and the inability to determine the satellite cross-domain information in a timely manner.

[0026] As well as in related technologies, when calculating satellite cross-domain, only a relatively single type of ground area is supported, ignoring the complex situations of various shapes and complex polygonal areas. For example: a fast and high-precision calculation method for the transit period of a space target in the related technology uses a large-step recursive method to calculate orbital prediction data, calculate the target sub-satellite point and analyze the transit status of the target at each moment, determine the sampling point positions before and after the target sub-satellite point passes, and combines the principle of linear interpolation to calculate the time period when the target sub-satellite point passes through a circular ground area. It is not applicable to various polygonal areas.

[0027] An embodiment of the present invention provides a method for determining satellite cross-domain information, comprising: determining multiple target sub-satellite points related to the area to be measured from multiple sub-satellite points included in a search engine based on area location information of the area to be measured; dividing the multiple target sub-satellite points into at least one sub-satellite point sequence based on timestamps of the target sub-satellite points and satellites related to the target sub-satellite points; performing boundary point verification on the sub-satellite point sequence to obtain a verification result; and determining the satellite cross-domain information of the area to be measured based on the verification result of each of the at least one sub-satellite point sequences.

[0028] Figure 1 An application scenario diagram of a satellite cross-domain information determination method, apparatus, device, medium, and program product according to an embodiment of the present invention is shown.

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

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

[0031] The first terminal device 101, the second terminal device 102, and the third terminal device 103 can be various electronic devices with a display screen and supporting web browsing, including but not limited to smartphones, tablets, laptop computers, desktop computers, and so on.

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

[0033] It should be noted that the satellite cross-domain information determination method provided by the embodiments of the present invention can generally be executed by the server 105. Correspondingly, the satellite cross-domain information determination device provided by the embodiments of the present invention can generally be set in the server 105. The satellite cross-domain information determination method provided by the embodiments of the present invention can also be executed by a server or a server cluster different from the server 105 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or the server 105. Correspondingly, the satellite cross-domain information determination device provided by the embodiments of the present invention can also be set in a server or a server cluster different from the server 105 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or the server 105.

[0034] It should be understood that Figure 1 the numbers of terminal devices, networks, and servers in

[0035] are merely illustrative. According to the implementation requirements, there can be any number of terminal devices, networks, and servers. Figure 1 described scenarios, through Figures 2 - 3 Figures 4(a) to 4(d), the satellite cross-domain information determination method of the embodiments of the invention will be described in detail.

[0036] Figure 2 shows a flowchart of the satellite cross-domain information determination method according to an embodiment of the present invention.

[0037] As Figure 2 shown, the method includes operations S210 to S240.

[0038] In operation S210, based on the regional position information of the area to be measured, a plurality of target sub-satellite points related to the area to be measured are determined from among the plurality of sub-satellite points included in the search engine.

[0039] In operation S220, based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points, the plurality of target sub-satellite points are divided into at least one sub-satellite point sequence.

[0040] In operation S230, boundary point verification is performed on the sub-satellite point sequence to obtain a verification result.

[0041] In operation S240, based on the verification results of each of the at least one sub-satellite point sequence, the satellite cross-domain information of the area to be measured is determined.

[0042] According to an embodiment of the present invention, the area to be measured may be an actual geographical area of the Earth, the sub-satellite point may be the point projected by the satellite on the Earth's surface, and the satellite trajectory of a satellite may include a plurality of sub-satellite points.

[0043] According to an embodiment of the present invention, the regional position information may be information used to characterize the regional position of the area to be measured. For example: when the regional shape of the area to be measured is circular, the regional position information may include the position of the center of the circle of the area to be measured, the radius value of the area to be measured, etc.; when the regional shape of the area to be measured is rectangular, the regional position information may include the positions of the upper left vertex and the lower right vertex of the area to be measured; when the regional shape of the area to be measured is polygonal, the regional position information may include the positions of a plurality of vertices of the area to be measured.

[0044] For example: when the shape of the area to be measured is circular, the area position information can be defined as: "central point": "E112°18'22\", N22°10'59\""; "radius": "421.61". When the shape of the area to be measured is rectangular, the area position information can be defined as: "lef up point": "E117°22'30\", N23°52'25\""; "right dow point": "E119°S3'17\", N19°34'49\"". When the shape of the area to be measured is polygonal, the area position information can be defined as: "polygon point": "E89°14'6\", N21°11'40\"; E80°51'22\", N15°28'7\"; E71°39'24\", N20°58'33\"; E74°10'35\", N26°51'50\"; E77°38'0\", N24°38'52\"; E83°11'59\", N27°7'30\"".

[0045] According to an embodiment of the present invention, the search engine can be an engine that searches for the sub-satellite point located in the area to be measured. Due to the hierarchical filtering characteristics of the search engine, the sub-satellite point in the preset area matching the area to be measured can be determined, so as to realize the filtering of the sub-satellite point, and then more quickly determine the target sub-satellite point related to the area to be measured.

[0046] According to an embodiment of the present invention, the target sub-satellite point can be the point on the satellite trajectory passing through the area to be measured.

[0047] According to an embodiment of the present invention, the method for obtaining the sub-satellite point is not limited. It can be obtained by recording the real satellite trajectory of the satellite, or can be obtained by predicting the satellite trajectory of the satellite through the orbital characteristic data of the satellite.

[0048] According to an embodiment of the present invention, if the sub-satellite point is obtained by recording the real satellite trajectory, the time stamp of the target sub-satellite point is the recording moment of the target sub-satellite point. If the sub-satellite point is obtained by predicting the satellite trajectory, the time stamp of the target sub-satellite point is the prediction moment of the target sub-satellite point.

[0049] According to an embodiment of the present invention, the cross-domain trajectory corresponding to the target sub-satellite point can be determined according to the time stamp of each target sub-satellite point and the satellite to which the target sub-satellite point belongs, so as to obtain at least one sub-satellite point sequence. The cross-domain trajectory is the trajectory segment on the satellite trajectory passing through the area to be measured. One cross-domain trajectory corresponds to one sub-satellite point sequence, and the cross-domain behaviors of the same satellite at different time periods belong to different cross-domain trajectories.

[0050] According to an embodiment of the present invention, by dividing multiple target sub-satellite points into at least one sub-satellite point sequence, the division of cross-domain trajectories can be achieved, and by respectively performing boundary point verification on each sub-satellite point sequence, each cross-domain trajectory can be verified, so as to more accurately determine the satellite cross-domain information of the area to be measured.

[0051] According to an embodiment of the present invention, the boundary point verification can be performed by verifying the boundary points located at the starting position and the ending position in the sub-satellite point sequence.

[0052] According to an embodiment of the present invention, the satellite cross-domain information of the area to be measured is not limited, and can be the specific information of the satellites passing through the area to be measured within the target time period, and can include at least one of the following: there are n satellites passing through the area to be measured within the target time period, the entry time, the exit time, the passing duration of each satellite, etc. The target time period can be a future time period or a historical time period.

[0053] According to the satellite cross-domain information determination method of the present invention, by using the spatial indexing technology of the search engine, multiple target sub-satellite points related to the area to be measured are quickly screened from multiple sub-satellite points. And by using the time stamp and the satellites related to the target sub-satellite points, the sub-satellite point sequence is divided to construct at least one cross-domain trajectory. And the boundary point verification is respectively performed on the sub-satellite point sequences, and based on the verification results, the satellite cross-domain information of the area to be measured is determined. Since after determining multiple target sub-satellite points related to the area to be measured through the search engine, multiple target sub-satellite points are divided into at least one sub-satellite point sequence by using the time stamp and the satellites related to the target sub-satellite points, and more accurate satellite cross-domain information is obtained by performing the boundary point verification of the sub-satellite point sequence, so that the satellite cross-domain information can be obtained by calculating some of the sub-satellite points, thus at least partially solving the technical problem of low calculation efficiency of satellite cross-domain information in the related art, and realizing the improvement of the determination speed of satellite cross-domain information while ensuring the accuracy of the information.

[0054] According to an embodiment of the present invention, based on the time stamp of the target sub-satellite points and the satellites related to the target sub-satellite points, dividing multiple target sub-satellite points into at least one sub-satellite point sequence may include the following operations.

[0055] Based on the satellites related to the target sub-satellite points, multiple target sub-satellite points are divided to obtain at least one candidate point group, wherein the multiple candidate sub-satellite points included in the candidate point group all belong to the same satellite; for each candidate point group, the multiple candidate sub-satellite points are sorted according to their respective time stamps to obtain a candidate point sequence; based on the time interval between adjacent candidate sub-satellite points in the candidate point sequence, the multiple candidate sub-satellite points included in the candidate point sequence are respectively divided into their respective cross-domain trajectories to obtain at least one sub-satellite point sequence.

[0056] According to an embodiment of the present invention, multiple target sub-satellite points can be grouped by satellites respectively related to the target sub-satellite points, so as to obtain candidate point groups.

[0057] According to an embodiment of the present invention, by sorting multiple candidate sub-satellite points according to timestamps, a candidate point sequence is obtained, so that the arrangement order of each candidate sub-satellite point in the candidate point sequence is the same as the order in the satellite cross-domain trajectory, facilitating the distinction of different cross-domain trajectories of the same satellite.

[0058] According to an embodiment of the present invention, there is no limitation on the way of arranging candidate sub-satellite points according to timestamps. It can be arranged in ascending order according to timestamps, or in descending order according to timestamps.

[0059] According to an embodiment of the present invention, since if two adjacent candidate sub-satellite points belong to different cross-domain trajectories, the time interval between them will be greater than the transit trajectories of two normal adjacent candidate sub-satellite points, different cross-domain trajectories can be distinguished.

[0060] According to an embodiment of the present invention, different satellite trajectories can be clearly distinguished by grouping, avoiding confusion between trajectories. And by sorting and analyzing the time interval, different cross-domain trajectories of the same satellite can be accurately distinguished, ensuring the continuity and accuracy of the trajectories, and improving the accuracy of satellite cross-domain information analysis.

[0061] According to an embodiment of the present invention, based on the time interval between adjacent candidate sub-satellite points in the candidate point sequence, multiple candidate sub-satellite points included in the candidate point sequence are respectively divided into their respective cross-domain trajectories to obtain at least one sub-satellite point sequence, which may include the following operations.

[0062] Determine the i-th candidate sub-satellite point from the candidate point sequence, where i is a positive integer greater than 1. In the case where the candidate sub-satellite point is the first candidate sub-satellite point in the candidate point sequence, construct a corresponding initial sub-satellite point sequence for the first candidate sub-satellite point; compare the timestamp of the i-th candidate sub-satellite point with the timestamp of the (i - 1)-th candidate sub-satellite point to determine the time interval between the i-th candidate sub-satellite point and the (i - 1)-th candidate sub-satellite point; in the case where it is determined that the time interval is less than the preset time threshold, add the i-th candidate sub-satellite point to the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point; repeat the above operations until multiple candidate sub-satellite points are all added to their respective corresponding initial sub-satellite point sequences to obtain at least one sub-satellite point sequence.

[0063] According to an embodiment of the present invention, each candidate sub-satellite point in the candidate point sequence can be traversed separately. If the time interval between the current i-th candidate sub-satellite point and the (i - 1)-th candidate sub-satellite point is less than a preset threshold, it can be considered that the i-th candidate sub-satellite point and the (i - 1)-th candidate sub-satellite point belong to the same cross-domain trajectory. Therefore, the i-th candidate sub-satellite point can be added to the initial sub-satellite point sequence where the (i - 1)-th candidate sub-satellite point is located.

[0064] According to an embodiment of the present invention, in the case where multiple sub-satellite points are obtained by predicting the satellite trajectory, the preset time threshold can be determined based on the time interval between adjacent sub-satellite points in the satellite trajectory when predicting the sub-satellite points. For example, if the time interval between adjacent sub-satellite points during prediction is 10s, the preset time threshold can be 11s.

[0065] According to an embodiment of the present invention, if the sub-satellite points are obtained by recording the real satellite trajectory, the preset time threshold can be determined based on the maximum time interval between adjacent sub-satellite points in the satellite trajectory. The sub-satellite points of the same satellite trajectory can be preprocessed, and the time interval between adjacent sub-satellite points can be no greater than a preset value.

[0066] According to an embodiment of the present invention, the above satellite cross-domain information determination method further includes the following operations.

[0067] In the case where it is determined that the time interval is greater than or equal to the preset time threshold, construct a new initial sub-satellite point sequence different from the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point; add the i-th candidate sub-satellite point to the new initial sub-satellite point sequence.

[0068] According to an embodiment of the present invention, in the case where it is determined that the time interval is greater than or equal to the preset time threshold, a new initial sub-satellite point sequence can be constructed to place the current i-th candidate sub-satellite point.

[0069] According to an embodiment of the present invention, in the case where it is determined that the time interval is greater than or equal to the preset time threshold, it can be considered that the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point is a sub-satellite point sequence, that is, the candidate sub-satellite points in this sub-satellite point sequence have been determined.

[0070] According to an embodiment of the present invention, by comparing the time intervals between adjacent candidate sub-satellite points, if the time interval is less than a preset time threshold, it is considered that the two belong to the same trajectory; if the time interval is greater than or equal to the preset time threshold, it is considered that the two belong to different trajectories. Thus, the sub-satellite points of candidate stars belonging to different cross-domain trajectories can be accurately divided, effectively avoiding trajectory confusion and ensuring the continuity and accuracy of each cross-domain trajectory. Moreover, by using local information, that is, by comparing the time intervals between adjacent candidate sub-satellite points in the candidate point sequence to divide the sub-satellite point sequence, the accumulation of errors can be reduced, and new candidate sub-satellite points can be processed in real time without reprocessing the entire data set, thereby improving the speed of cross-domain information determination.

[0071] According to an embodiment of the present invention, based on the regional position information of the area to be measured, determining a plurality of target sub-satellite points related to the area to be measured from among the plurality of sub-satellite points included in the search engine may include the following operations.

[0072] Based on the regional position information of the area to be measured, determine a query instruction for the area to be measured; input the query instruction and the regional position information into the search engine, so that the search engine, based on the regional position information, determines at least one preset area that matches the regional position information, and based on the coordinate values of the sub-satellite points located in each preset area, determines a plurality of target sub-satellite points located in the area to be measured; obtain the plurality of target sub-satellite points from the search engine.

[0073] According to an embodiment of the present invention, since the information included in the regional position information of the areas to be measured with different regional shapes is different, and different regional shapes correspond to different query instructions, therefore, based on the relationship between the above regional position information and the regional shape, and the relationship between the regional shape and the query instruction, a mapping relationship between the regional position information and the regional position information can be preset, so that the query instruction can be quickly determined through the regional position information.

[0074] According to an embodiment of the present invention, the search engine is not limited to a preset search engine embedded with a geographical location retrieval function, and when performing geographical location retrieval, the search engine can use a BKD tree (Block-K-Dimensional Tree), thereby dividing the space into a plurality of preset areas. When performing a query, the preset areas can be first matched with the area to be measured, and then further calculations can be performed on the sub-satellite points in the matched preset areas to determine whether the sub-satellite point is located in the area to be measured, that is, hierarchical filtering is achieved, saving a certain amount of computing resources and improving the computing speed.

[0075] According to an embodiment of the present invention, the search engine may also be Elasticsearch, etc. When the search engine is Elasticsearch, when the area shape of the area to be measured is circular, the query instruction may be geo_distance to find the sub-satellite points within a certain range from a certain center point; when the area shape is rectangular, the query instruction may be geo_bounding_box to find the sub-satellite points located within the specified rectangular box; when the area shape is polygonal, the query instruction may be geo_polygon to find the sub-satellite points located within the specified polygon.

[0076] According to an embodiment of the present invention, a query statement is automatically generated based on the area position information corresponding to the area to be measured, and the search engine is used to implement the rapid retrieval of geographical information, so as to quickly determine the target sub-satellite points passing through the area to be measured, and improve the rate of satellite cross-domain information; in addition, the target sub-satellite points can be quickly determined for different shapes of the area to be measured, at least partially solving the problem in the related art that only supports relatively single shapes of ground area types and ignores various shapes and complex polygon areas.

[0077] According to an embodiment of the present invention, the coordinate values of the sub-satellite points are determined in the following manner.

[0078] Obtain the orbital characteristic data of at least one satellite within the target time period; for each satellite, based on a preset prediction model and the orbital characteristic data, predict the satellite trajectory of the satellite within the future time period to obtain the initial coordinate values respectively corresponding to multiple sub-satellite points in the inertial coordinate system; perform an Earth coordinate system conversion on the initial coordinate values respectively corresponding to multiple sub-satellite points to obtain the coordinate values respectively corresponding to multiple sub-satellite points.

[0079] According to an embodiment of the present invention, the acquisition method of the orbital characteristic data is not limited, and it may be obtained through a TLE (Two-Line Element Set) data set. Specifically, the orbital characteristic data may include: semi-major axis, eccentricity, orbital inclination, argument of perigee, right ascension of the ascending node, and true anomaly, and these orbital characteristic data can be obtained by parsing the TLE data set or calculating using the data in the TLE data set.

[0080] According to an embodiment of the present invention, the TLE data set can be obtained from the public website Space-Track.

[0081] According to an embodiment of the present invention, the data and data locations included in the TLE data set are shown in Table 1.

[0082] Table 1

[0083]

[0084] According to an embodiment of the present invention, the preset prediction model is not limited and can be a celestial dynamics model (Simplified General Perturbations Model Version 4 / Simplified Deep Space Perturbations Model Version 4, SGP4 / SDP4).

[0085] According to an embodiment of the present invention, the celestial dynamics model mainly calculates the accurate position information of the satellite at a specified time point based on the orbital parameters of the satellite and considering the influence of time dependence and various perturbation force factors.

[0086] According to an embodiment of the present invention, the time interval of the orbit prediction, that is, the time interval between adjacent sub-satellite points, is not limited and can be 10 seconds.

[0087] According to an embodiment of the present invention, based on the orbital characteristic data and the SGP4 / SDP4 model, the initial coordinate value of the sub-satellite point of the satellite in the inertial coordinate system can be calculated, and the initial coordinate value can be converted into the coordinate value in the earth coordinate system, that is, the longitude and latitude plus altitude coordinates.

[0088] According to an embodiment of the present invention, when the search engine is Elasticsearch, the respective coordinate values of each sub-satellite point can be stored in the search engine in the geo_point type.

[0089] According to an embodiment of the present invention, the orbital characteristic data is obtained by parsing the TLE data, the sub-satellite point coordinates of the satellite are predicted using the celestial kinematics model, and the prediction time interval is customized. Since the TLE data set is rich in resources and open source, a comprehensive and accurate sub-satellite point data set can be constructed and synchronized to the search cluster, which can achieve better data support for determining the satellite over-domain information and facilitate data update.

[0090] According to an embodiment of the present invention, boundary point verification is performed on the sub-satellite point sequence to obtain a verification result, which may include the following operations.

[0091] Determine at least two boundary points located at the boundary positions from the sub-satellite point sequence; for each boundary point, based on the boundary point type determined by the boundary position of the boundary point, determine the verification point for verifying the boundary point from the satellite trajectory of the satellite related to the boundary point, where the satellite trajectory includes a plurality of other sub-satellite points except the target sub-satellite point related to the satellite, and the verification point is determined from the plurality of other sub-satellite points; verify the boundary point based on the verification rule and the verification point for the regional shape of the area to be measured to obtain a verification sub-result; determine the verification result based on the respective verification sub-results of the at least two boundary points.

[0092] According to an embodiment of the present invention, the boundary position may include at least a starting position and an ending position.

[0093] According to an embodiment of the present invention, when the sub-satellite point sequence is arranged in ascending order of time stamps, the boundary point type of the target sub-satellite point at the starting position in the sub-satellite point sequence is an in-domain point, and the boundary point type of the target sub-satellite point at the ending position in the sub-satellite point sequence is an out-domain point.

[0094] According to an embodiment of the present invention, when the sub-satellite point sequence is arranged in descending order of time stamps, the boundary point type of the target sub-satellite point at the starting position in the sub-satellite point sequence is an out-domain point, and the boundary point type of the target sub-satellite point at the ending position in the sub-satellite point sequence is an in-domain point.

[0095] According to an embodiment of the present invention, the in-domain point is the first sub-satellite point in the satellite trajectory that passes through the area to be measured, and its previous sub-satellite point in the satellite trajectory is outside the area to be measured. Therefore, the previous point of this in-domain point in the satellite trajectory can be used as the verification point of this in-domain point.

[0096] According to an embodiment of the present invention, the out-domain point is the last sub-satellite point in the satellite trajectory that passes through the area to be measured, and its next sub-satellite point is outside the area to be measured. Therefore, the next sub-satellite point of this out-domain point in the satellite trajectory can be used as the verification point of this out-domain point.

[0097] According to an embodiment of the present invention, different verification rules may be adopted for different area shapes. Through the verification rules and verification points, the verification of the boundary points in the sub-satellite point sequence can be realized, thereby the verification of the sub-satellite point sequence can be realized, and further the accurate time points when the satellite enters and leaves the area to be measured can be determined, avoiding misjudgment or omission, and improving the accuracy of satellite cross-domain information.

[0098] According to an embodiment of the present invention, based on the verification rules and verification points for the area shape of the area to be measured, verifying the boundary points to obtain a verification result may include the following operations.

[0099] When the shape of the region is circular, determining the first coordinate value of the verification point and the second coordinate value of the center of the region to be measured included in the region position information, where the first coordinate value includes a first longitude coordinate and a first latitude coordinate, and the second coordinate value includes a second longitude coordinate and a second latitude coordinate; respectively performing radian conversion on the first longitude coordinate, the first latitude coordinate, the second longitude coordinate, and the second latitude coordinate to obtain a first longitude radian, a first latitude radian, a second longitude radian, and a second latitude radian; determining the distance value between the verification point and the center of the circle based on a first radian difference between the first longitude radian and the second longitude radian, a second radian difference between the first latitude radian and the second latitude radian, the first latitude radian, and the second latitude radian; determining the verification result for the boundary point based on the comparison result between the distance value and the radius value of the region to be measured included in the region position.

[0100] According to an embodiment of the present invention, when the shape of the region is circular, it is possible to obtain a verification sub-result for the boundary point by determining whether the distance from the verification point to the center of the region to be measured is less than the radius value of the region to be measured.

[0101] According to an embodiment of the present invention, specifically, when the distance from the verification point to the center of the region to be measured is less than or equal to the radius value of the region to be measured, it can be considered that the verification sub-result indicates verification failure; when the distance from the verification point to the center of the region to be measured is greater than the radius value of the region to be measured, it can be considered that the verification sub-result indicates verification success.

[0102] According to an embodiment of the present invention, performing radian conversion on the first coordinate value and the second coordinate value and calculating the first radian difference and the second radian difference can be as shown in the following formula (1).

[0103] (1);

[0104] Wherein, is the first latitude radian, is the first latitude coordinate, is the second latitude radian, is the second latitude coordinate, is the first radian difference, is the second radian difference, is the first longitude coordinate, is the second longitude coordinate.

[0105] According to an embodiment of the present invention, determining the distance value between the verification point and the center of the circle based on a first radian difference between the first longitude radian and the second longitude radian, a second radian difference between the first latitude radian and the second latitude radian, the first latitude radian, and the second latitude radian can be as shown in the following formula (2).

[0106] (2);

[0107] Among them, d is the distance value between the check point and the center of the circle, and R is the radius of the earth.

[0108] According to an embodiment of the present invention, formula (2) is derived from the following formulas (3) to (6).

[0109] (3);

[0110] (4);

[0111] (5);

[0112] (6);

[0113] Among them, is the haversine formula, is a replaceable preset value, is an intermediate parameter.

[0114] According to an embodiment of the present invention, the actual distance between two points is calculated through the longitude and latitude coordinates of the check point and the center of the circle of the area to be measured. Since the curvature of the earth is considered, the distance between two geographical points can be calculated more accurately, thereby making the check sub-result more accurate.

[0115] According to an embodiment of the present invention, based on the check rules and check points for the area shape of the area to be measured, the boundary points are checked to obtain a check result, which may include the following operations.

[0116] In the case where the area shape is a polygon or a rectangle, based on multiple area vertices of the area to be measured in the area position information, multiple vertex groups for forming the area edges of the area to be measured are determined, where the vertex group includes at least two area vertices; for each vertex group, determine the first latitude difference between the first latitude coordinate of the check point and the third latitude coordinate of the first area vertex in the vertex group, and the second latitude difference between the first latitude coordinate and the fourth latitude coordinate of the second area vertex in the vertex group; in the case where it is determined that the first latitude difference and the second latitude difference meet the preset conditions, based on the longitude difference between the third longitude coordinate of the first area vertex and the fourth longitude coordinate of the second area vertex, the third latitude difference between the third latitude coordinate and the fourth latitude coordinate, the first latitude difference and the third longitude coordinate, determine the check longitude value; in the case where the check longitude value is greater than or equal to the first longitude coordinate of the check point, determine that there is an intersection between the ray generated with the check point as the starting point and the area edge formed by the vertex group; based on the total number of intersections between the area edges formed by multiple vertex groups and the ray generated with the check point as the starting point, determine the check result.

[0117] According to an embodiment of the present invention, the polygon can be a regular or irregular polygon.

[0118] According to an embodiment of the present invention, in the case where the area shape is a rectangle, the coordinate values of other area vertices can be determined based on the coordinate values of the upper left corner vertex and the lower right corner vertex included in the area position information, so as to obtain multiple area vertices of the area to be measured.

[0119] According to an embodiment of the present invention, in the case where the area shape is a polygon, multiple area vertices of the area to be measured can be stored in the target database in the actual connection order.

[0120] According to an embodiment of the present invention, by determining the vertex groups of the area edges that can form the area to be measured based on the storage order of multiple area vertices, and based on the coordinates of the area vertices in each vertex group, it is determined whether there is an intersection between the ray emitted from the verification point and the area edge formed by the vertex group.

[0121] According to an embodiment of the present invention, area vertices adjacent in the storage order can be formed into vertex groups, and area vertices with the first and last storage orders can be formed into vertex groups.

[0122] According to an embodiment of the present invention, by counting the total number of intersections between the verification point and the area edges respectively formed by multiple vertex groups and the ray generated with the verification point as the starting point, the verification sub-result can be determined.

[0123] According to an embodiment of the present invention, in the case where the total number of intersections is odd, it can be determined that the verification point is inside the area to be measured, that is, the verification sub-result is verification failure. In the case where the total number of intersections is even, it can be determined that the verification point is outside the area to be measured, that is, the verification sub-result is verification success.

[0124] According to an embodiment of the present invention, by respectively comparing the latitude differences of the area vertices in the vertex group of the verification point, it can be determined whether the first area vertex and the second area vertex are on both sides of the verification point. In the case where it is determined that the first area vertex and the second area vertex are on both sides of the verification point, it is determined that there may be an intersection between the ray generated by the verification point and the area edge formed by the first area vertex and the second area vertex. In the case where it is determined that the first area vertex and the second area vertex are not on both sides of the verification point, it is determined that there is no intersection between the ray generated by the verification point and the area edge formed by the first area vertex and the second area vertex. Thus, in this case, the subsequent intersection determination steps can be omitted, realizing the screening of vertex groups and improving the calculation speed.

[0125] According to an embodiment of the present invention, the situation where the first latitude difference and the second latitude difference satisfy a preset condition may be that the first latitude coordinate of the verification point determined by the first latitude difference and the second latitude difference is greater than the latitude coordinate of one of the vertices of the first region or the second region and less than the latitude coordinate of one of them.

[0126] According to an embodiment of the present invention, for example: in the case where the first latitude difference and the second latitude difference are obtained by subtracting the third latitude coordinate from the first latitude coordinate and subtracting the fourth latitude coordinate from the first latitude coordinate respectively, if one of the first latitude difference and the second latitude difference is greater than 0 and the other is less than 0, it is determined that the preset condition is satisfied.

[0127] According to an embodiment of the present invention, based on the longitude difference between the third longitude coordinate of the vertex of the first region and the fourth longitude coordinate of the vertex of the second region, the third latitude difference between the third latitude radian and the fourth latitude radian, the first latitude difference, and the third longitude coordinate, the verification longitude value can be determined as shown in the following formula (7).

[0128] (7);

[0129] Wherein, is the verification longitude value, is the first latitude coordinate of the verification point, is the third latitude coordinate, is the fourth latitude coordinate, is the third longitude coordinate, is the fourth longitude coordinate.

[0130] According to an embodiment of the present invention, through formula (7), the intersection situation between the ray in a fixed direction generated with the verification point as the starting point and the region edge can be simulated, and the fixed direction can be horizontal to the right, etc.

[0131] According to an embodiment of the present invention, in the case where it is determined that the verification longitude value is greater than or equal to the first longitude coordinate, it can be determined that the intersection point is on the right side of the ray and is valid, and it is considered that there is an intersection point between the ray generated with the verification point as the starting point and the region edge formed by the vertex group. Otherwise, there is no intersection point.

[0132] According to an embodiment of the present invention, by analyzing the latitude difference between the verification point and the region vertices in the vertex group of the region edges that can form the region to be measured, the screening of the vertex group that can have intersection points can be realized, thereby saving a certain amount of calculation, and using the coordinate values of the verification point, the vertices of the first region, and the vertices of the second region respectively to measure whether there is an intersection point between the ray formed by the verification point and the region edge, so that the intersection situation can be determined through a convenient calculation process, improving the determination speed, and accurately judging whether the verification point is within the region to be measured by counting the number of intersection points, and maintaining high accuracy even in complex polygon regions.

[0133] According to an embodiment of the present invention, for a region to be measured with different regional shapes, different query instructions and verification rules are adopted, which supports regions with complex shapes. Multiple target sub-satellite points related to the time stamp and the target sub-satellite point are divided into at least one sub-satellite point sequence, so as to effectively distinguish multiple cross-domain events, that is, cross-domain trajectories. And it can realize batch online monitoring of specific regions and specific types of satellites.

[0134] According to an embodiment of the present invention, based on the verification results of each of the at least one sub-satellite point sequence, determining the satellite cross-domain information of the region to be measured may include the following operations.

[0135] For each sub-satellite point sequence, when the verification result indicates that the verification of the target boundary point among at least two boundary points fails, based on the time stamp of the target boundary point, a new boundary point is determined from the satellite trajectory to which the target boundary point belongs; when the verification sub-result of determining the new boundary point is successful, the sub-satellite point sequence is updated based on the new boundary point to obtain an updated sub-satellite point sequence; based on the updated sub-satellite point sequence, the satellite cross-domain information is determined. According to an embodiment of the present invention, the target boundary point is determined based on the verification sub-result, and the boundary point corresponding to the verification sub-result indicating verification failure is the target boundary point.

[0136] According to an embodiment of the present invention, the new boundary point may be the verification point of the target boundary point. Or when the target boundary point is an in-domain point, the point before the in-domain point in the satellite trajectory may be used as the new boundary point; when the target boundary point is an out-domain point, the point after the out-domain point in the satellite trajectory may be used as the new boundary point.

[0137] According to an embodiment of the present invention, the new boundary point can be verified in the same verification manner as the boundary point, so as to determine the verification sub-result of the new boundary point. When the verification sub-result indicates verification success, the target boundary point in the sub-satellite point sequence is replaced with the new boundary point, so as to realize the update of the sub-satellite point sequence and obtain an updated sub-satellite point sequence. When the verification sub-result indicates verification failure, the above process is repeated for the new boundary point until a boundary point with successful verification is obtained, and all the new boundary points determined during the repeated process are added to the sub-satellite point sequence to obtain an updated sub-satellite point sequence.

[0138] According to an embodiment of the present invention, the in-domain time, out-domain time, cross-domain duration, etc. of the satellite can be determined according to the updated sub-satellite point sequence. By counting the sub-satellite point sequence and the updated sub-satellite point sequence, it can be determined how many satellites pass through the region to be measured within the target time period.

[0139] According to an embodiment of the present invention, the satellite cross-domain information can be displayed using a display interface.

[0140] According to an embodiment of the present invention, when the subsatellite point is obtained by predicting the satellite trajectory, since the satellite trajectory is dynamic, the satellite cross-domain information may change with the orbit change. In this embodiment, when predicting the subsatellite point of the satellite, the prediction result can be dynamically updated according to the real-time acquired TLE dataset. Therefore, the satellite cross-domain information can be dynamically changed to monitor the satellite in real time. According to the set backward trajectory prediction duration and backward monitoring time threshold, the specific information of the satellites passing through the area to be measured in the future can be monitored, such as how many satellites pass through the area in the future, and the entry point, exit point, and residence time in the area of a certain cross-domain satellite, etc.

[0141] Figure 3 The flowchart of determining the subsatellite point sequence according to an embodiment of the present invention is shown;

[0142] As Figure 3 shown, determining the subsatellite point sequence includes operations S301 to S311.

[0143] In operation S301, based on the satellites related to the target subsatellite point, multiple target subsatellite points are divided to obtain at least one candidate point group, where the multiple candidate subsatellite points included in the candidate point group all belong to the same satellite.

[0144] In operation S302, for each candidate point group, the multiple candidate subsatellite points are sorted according to their respective timestamps to obtain a candidate point sequence.

[0145] In operation S303, the i-th candidate subsatellite point is determined from the candidate point sequence, where i is a positive integer greater than 1. When the candidate subsatellite point is the first candidate subsatellite point in the candidate point sequence, an initial subsatellite point sequence corresponding to the first candidate subsatellite point is constructed.

[0146] In operation S304, the timestamp of the i-th candidate subsatellite point is compared with the timestamp of the (i - 1)-th candidate subsatellite point to determine the time interval between the i-th candidate subsatellite point and the (i - 1)-th candidate subsatellite point.

[0147] In operation S305, it is determined whether the time interval is less than a preset time threshold. When it is determined that the time interval is greater than or equal to the preset time threshold, operation S306 is executed. When it is determined that the time interval is less than the preset time threshold, operation S311 is executed.

[0148] In operation S306, a new initial subsatellite point sequence different from the initial subsatellite point sequence corresponding to the (i - 1)-th candidate subsatellite point is constructed.

[0149] In operation S307, the i-th candidate subsatellite point is added to the new initial subsatellite point sequence.

[0150] In operation S308, it is determined whether the i-th candidate sub-satellite point is the last candidate sub-satellite point in the candidate point sequence. In the case where it is determined that the i-th candidate sub-satellite point is the last candidate sub-satellite point in the candidate point sequence, operation S309 is executed. In the case where it is determined that the i-th candidate sub-satellite point is not the last candidate sub-satellite point in the candidate point sequence, operation S310 is executed.

[0151] In operation S309, at least one sub-satellite point sequence is determined based on at least one initial sub-satellite point sequence.

[0152] In operation S310, let i be equal to i + 1, and return to operation S303.

[0153] In operation S311, the i-th candidate sub-satellite point is added to the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point, and return to operation S308.

[0154] FIG. 4(a) shows a first schematic diagram of satellite cross-domain information according to an embodiment of the present invention.

[0155] As shown in FIG. 4(a), the display interface shows that the shape of the area to be measured is circular, and its satellite cross-domain information includes: it is expected that two satellites will pass through the area to be measured within the next 3 hours.

[0156] According to an embodiment of the present invention, the satellite cross-domain information may further include: basic information of the satellite and cross-domain trajectory information. For example: by selecting the basic information of the target satellite to be displayed from the satellites expected to pass through the area to be measured, such as the satellite identity document (ID) is ID1, the height is H1, the speed is V1, the latitude is N1, the longitude is W1, the name is S1, the current time is 2025-02-10 10:36:13, the launch time is =T1, the location it belongs to is a1, the classification is b1, and the launch location is c1. And the displayed cross-domain trajectory information of the target satellite includes: satellite ID1 is expected to enter the area to be measured at 2025-02-10 13:07:40, leave the area to be measured at 2025-02-10 13:11:50, and stay for 4.17 minutes.

[0157] FIG. 4(b) shows a second schematic diagram of satellite cross-domain information according to an embodiment of the present invention.

[0158] As shown in FIG. 4(b), the display interface shows that the shape of the area to be measured is rectangular, and its satellite cross-domain information includes: it is expected that 36 satellites will pass through the area to be measured within the next 1 hour.

[0159] According to an embodiment of the present invention, the basic information of the determined target satellite includes: ID is ID2, altitude is H2, speed is V2, latitude is N2, longitude is W2, name is S2, current time is 2025-02-10 10:41:17, launch time is T2, the location it belongs to is a2, classification is b2, and the launch location is c2. The cross-domain trajectory information includes: Satellite ID2 is expected to enter this area at 2025-02-10 11:09:40, leave this area at 2025-02-10 11:12:40, and stay for 3.00 minutes.

[0160] Figure 4 (c) shows a third schematic diagram of satellite cross-domain information according to an embodiment of the present invention.

[0161] As shown in Figure 4 (c), the display interface shows that the shape of the area to be measured is a concave polygon, and its satellite cross-domain information includes: It is expected that 4 satellites will pass through the area to be measured within the next 3 hours.

[0162] According to an embodiment of the present invention, the basic information of the determined target satellite includes: Satellite ID is ID3, altitude is H3, speed is V3, latitude is N3, longitude is W3, name is S3, current time is 2025-02-10 10:39:23, launch time is T3, the location it belongs to is a3, classification is b3, and the launch location is c3. The cross-domain trajectory information includes: Satellite ID3 is expected to enter this area at 2025-02-10 12:06:30, leave this area at 2025-02-10 12:09:50, and stay for 3.33 minutes. Enter this area at 2025-02-10 12:12:20, leave this area at 2025-02-10 12:13:50, and stay for 1.50 minutes. Satellite ID3 is expected to enter this area at 2025-02-10 11:09:40, leave this area at 2025-02-10 11:12:40, and stay for 3.00 minutes.

[0163] According to an embodiment of the present invention, as shown in Figure 4 (c), when the target satellite passes through the area to be measured, two cross-domain behaviors occur, that is, it includes two cross-domain trajectories.

[0164] Figure 4 (d) shows a fourth schematic diagram of satellite cross-domain information according to an embodiment of the present invention.

[0165] As shown in Figure 4 (d), the display interface shows that the shape of the area to be measured is a complex polygon, and its satellite cross-domain information includes: It is expected that 23 satellites will pass through the area to be measured within the next 3 hours.

[0166] According to an embodiment of the present invention, the basic information of the target satellite determined includes: ID is ID4, altitude is H4, speed is V4, latitude is N4, longitude is W4, name is S4, current time is 2025-02-10 10:38:01, launch time is T4, location is a4, classification is b4, and launch location is c4. The cross-domain trajectory information includes: Satellite ID4 is expected to enter this area at 2025-02-10 12:06:30, leave this area at 2025-02-10 12:11:00, and stay for 4.50 minutes. Enter this area at 2025-02-10 13:41:50, leave this area at 2025-02-10 13:43:30, and stay for 1.67 minutes.

[0167] Based on the above satellite cross-domain information determination method, the present invention also provides a satellite cross-domain information determination device. Figure 5 The device is described in detail.

[0168] Figure 5 A structural block diagram of a device for determining satellite cross-domain information according to an embodiment of the present invention is shown.

[0169] like Figure 5 As shown, the satellite cross-domain information determination device 500 of this embodiment includes a first determination module 510 , a division module 520 , a verification module 530 and a second determination module 540 .

[0170] The first determination module 510 is configured to determine a plurality of target sub-satellite points related to the area to be measured from a plurality of sub-satellite points included in the search engine based on the area location information of the area to be measured.

[0171] The division module 520 is used to divide the multiple target sub-satellite points into at least one sub-satellite point sequence based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points.

[0172] The verification module 530 is used to perform boundary point verification on the sub-satellite point sequence to obtain a verification result.

[0173] The second determination module 540 is configured to determine the satellite cross-domain information of the area to be measured based on the verification result of each of the at least one sub-satellite point sequences.

[0174] According to an embodiment of the present invention, the partitioning module 520 includes a first partitioning submodule, a sorting submodule and a second partitioning submodule.

[0175] The first division submodule is used to divide multiple target sub-satellite points based on satellites related to the target sub-satellite points to obtain at least one candidate point group, wherein the multiple candidate sub-satellite points included in the candidate point group all belong to the same satellite.

[0176] A sorting sub-module, for each candidate point group, sorts multiple candidate sub-satellite points according to their respective timestamps to obtain a candidate point sequence.

[0177] A second partitioning sub-module, based on the time interval between adjacent candidate sub-satellite points in the candidate point sequence, partitions the multiple candidate sub-satellite points included in the candidate point sequence into their respective cross-domain trajectories respectively to obtain at least one sub-satellite point sequence.

[0178] According to an embodiment of the present invention, the second partitioning sub-module includes a first determination unit, a comparison unit, a first addition unit, and a second addition unit.

[0179] The first determination unit is used to determine the i-th candidate sub-satellite point from the candidate point sequence, where i is a positive integer greater than 1. In the case where the candidate sub-satellite point is the first candidate sub-satellite point in the candidate point sequence, construct a corresponding initial sub-satellite point sequence for the first candidate sub-satellite point.

[0180] The comparison unit is used to compare the timestamp of the i-th candidate sub-satellite point with the timestamp of the (i - 1)-th candidate sub-satellite point to determine the time interval between the i-th candidate sub-satellite point and the (i - 1)-th candidate sub-satellite point.

[0181] The first addition unit is used to add the i-th candidate sub-satellite point to the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point when it is determined that the time interval is less than a preset time threshold.

[0182] The second addition unit is used to repeatedly execute the above operations until multiple candidate sub-satellite points are all added to their respective corresponding initial sub-satellite point sequences to obtain at least one sub-satellite point sequence.

[0183] According to an embodiment of the present invention, the satellite cross-domain information determination device 500 further includes a construction module and a sub-satellite point addition module.

[0184] The construction module is used to construct a new initial sub-satellite point sequence different from the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point when it is determined that the time interval is greater than or equal to the preset time threshold.

[0185] The sub-satellite point addition module is used to add the i-th candidate sub-satellite point to the new initial sub-satellite point sequence.

[0186] According to an embodiment of the present invention, the first determination module 510 includes an instruction determination sub-module, an input sub-module, and an acquisition sub-module.

[0187] The instruction determination sub-module is used to determine a query instruction for the area to be measured based on the area position information of the area to be measured.

[0188] An input sub-module, configured to input a query instruction and regional location information into a search engine, so that the search engine determines at least one preset region that matches the regional location information based on the regional location information, and determines a plurality of target sub-satellite points located in the region to be measured based on the coordinate values of the sub-satellite points located in each preset region.

[0189] An acquisition sub-module, configured to acquire a plurality of target sub-satellite points from the search engine.

[0190] According to an embodiment of the present invention, a satellite cross-domain information determination device 500 includes a data determination module, a prediction module, and a coordinate system conversion module.

[0191] The data determination module is configured to acquire orbital characteristic data of at least one satellite within a target time period.

[0192] The prediction module is configured to, for each satellite, predict the satellite trajectory of the satellite within a future time period based on a preset prediction model and the orbital characteristic data, and obtain initial coordinate values respectively corresponding to a plurality of sub-satellite points in an inertial coordinate system.

[0193] The coordinate system conversion module is configured to perform an Earth coordinate system conversion on the initial coordinate values respectively corresponding to a plurality of sub-satellite points to obtain coordinate values respectively corresponding to a plurality of sub-satellite points.

[0194] According to an embodiment of the present invention, a verification module 530 includes a boundary point determination sub-module, a verification point determination sub-module, a verification sub-module, and a result determination sub-module.

[0195] The boundary point determination sub-module is configured to determine at least two boundary points located at boundary positions from a sub-satellite point sequence.

[0196] The verification point determination sub-module is configured to, for each boundary point, determine a verification point for verifying the boundary point from the satellite trajectory of the satellite associated with the boundary point based on the boundary point type determined by the boundary position of the boundary point, where the satellite trajectory includes a plurality of other sub-satellite points except the target sub-satellite points associated with the satellite, and the verification point is determined from the plurality of other sub-satellite points.

[0197] The verification sub-module is configured to verify the boundary points based on a verification rule for the regional shape of the region to be measured and the verification points to obtain a verification sub-result.

[0198] The result determination sub-module is configured to determine a verification result based on the verification sub-results of at least two boundary points respectively.

[0199] According to an embodiment of the present invention, the verification sub-module includes a first coordinate determination unit, a radian determination unit, a distance value determination unit, and a first sub-result determination unit.

[0200] The first coordinate determination unit is configured to determine the first coordinate value of the verification point and the second coordinate value of the center of the area to be measured included in the area position information when the area shape is circular, where the first coordinate value includes the first longitude coordinate and the first latitude coordinate, and the second coordinate value includes the second longitude coordinate and the second latitude coordinate.

[0201] The radian determination unit is configured to perform radian conversion on the first longitude coordinate, the first latitude coordinate, the second longitude coordinate, and the second latitude coordinate respectively to obtain the first longitude radian, the first latitude radian, the second longitude radian, and the second latitude radian.

[0202] The distance value determination unit is configured to determine the distance value between the verification point and the center of the circle based on the first radian difference between the first longitude radian and the second longitude radian, the second radian difference between the first latitude radian and the second latitude radian, the first latitude radian, and the second latitude radian.

[0203] The first sub-result determination unit is configured to determine the verification sub-result for the boundary point based on the comparison result between the distance value and the radius value of the area to be measured included in the area position.

[0204] According to an embodiment of the present invention, the verification sub-module further includes a vertex group determination unit, a latitude difference determination unit, a verification value determination unit, an intersection determination unit, and a second sub-result determination unit.

[0205] The vertex group determination unit is configured to determine multiple vertex groups for forming the area edges of the area to be measured based on multiple area vertices of the area to be measured in the area position information when the area shape is a polygon or a rectangle, where the vertex group includes at least two area vertices.

[0206] The latitude difference determination unit is configured to, for each vertex group, determine the first latitude difference between the first latitude coordinate of the verification point and the third latitude coordinate of the first area vertex in the vertex group, and the second latitude difference between the first latitude coordinate and the fourth latitude coordinate of the second area vertex in the vertex group.

[0207] The verification value determination unit is configured to, when it is determined that the first latitude difference and the second latitude difference meet the preset conditions, determine the verification longitude value based on the longitude difference between the third longitude coordinate of the first area vertex and the fourth longitude coordinate of the second area vertex, the third latitude difference between the third latitude coordinate and the fourth latitude coordinate, the first latitude difference, and the third longitude coordinate.

[0208] The intersection determination unit is configured to determine that there is an intersection between the ray generated starting from the verification point and the area edge formed by the vertex group when the verification longitude value is greater than or equal to the first longitude coordinate of the verification point.

[0209] A second sub-result determining unit, configured to determine a parity check sub-result based on the total number of intersection points between the regional edges formed by multiple vertex groups and the rays generated starting from the parity check points.

[0210] According to an embodiment of the present invention, the second determining module 540 includes a new boundary point determining sub-module, an updating sub-module, and an information determining sub-module.

[0211] The new boundary point determining sub-module is configured to, for each sub-satellite point sequence, when the parity check result indicates that the target boundary point check among at least two boundary points fails, determine a new boundary point from the satellite trajectory of the satellite to which the target boundary point belongs based on the timestamp of the target boundary point.

[0212] The updating sub-module is configured to, when the parity check sub-result for the determined new boundary point is a successful parity check, update the sub-satellite point sequence based on the new boundary point to obtain an updated sub-satellite point sequence.

[0213] The information determining sub-module is configured to determine satellite cross-domain information based on the updated sub-satellite point sequence.

[0214] According to an embodiment of the present invention, any multiple of the first determining module 510, the partitioning module 520, the parity check module 530, and the second determining module 540 may be combined and implemented in one module, or any one of them may be split into multiple modules. Or, at least part of the functions of one or more of these modules may be combined with at least part of the functions of other modules and implemented in one module. According to an embodiment of the present invention, at least one of the first determining module 510, the partitioning module 520, the parity check module 530, and the second determining module 540 may be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or any other reasonable way of integrating or packaging circuits, etc., implemented by hardware or firmware, or implemented in any one of the three implementation manners of software, hardware, and firmware, or in an appropriate combination of any several of them. Or, at least one of the first determining module 510, the partitioning module 520, the parity check module 530, and the second determining module 540 may be at least partially implemented as a computer program module, and when the computer program module is run, the corresponding functions may be executed.

[0215] Figure 6 A block diagram of an electronic device suitable for implementing the satellite cross-domain information determination method according to an embodiment of the present invention is shown.

[0216] As Figure 6As shown, an electronic device 600 according to an embodiment of the present invention includes a processor 601, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 602 or a program loaded from a storage section 608 into a random access memory (RAM) 603. The processor 601 may include, for example, a general-purpose microprocessor (e.g., CPU), an instruction set processor, and / or a related chipset, and / or a dedicated microprocessor (e.g., an application-specific integrated circuit (ASIC)), etc. The processor 601 may also include on-board memory for caching purposes. The processor 601 may include a single processing unit or multiple processing units for performing different actions of a method flow according to an embodiment of the present invention.

[0217] In the RAM 603, various programs and data required for the operation of the electronic device 600 are stored. The processor 601, the ROM 602, and the RAM 603 are connected to each other via a bus 604. The processor 601 performs various operations of a method flow according to an embodiment of the present invention by executing a program in the ROM 602 and / or the RAM 603. It should be noted that the program may also be stored in one or more memories other than the ROM 602 and the RAM 603. The processor 601 may also perform various operations of a method flow according to an embodiment of the present invention by executing a program stored in the one or more memories.

[0218] According to an embodiment of the present invention, the electronic device 600 may further include an input / output (I / O) interface 605, and the input / output (I / O) interface 605 is also connected to the bus 604. The electronic device 600 may further include one or more of the following components connected to the input / output (I / O) interface 605: an input section 606 including a keyboard, a mouse, etc.; an output section 607 including, for example, a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 608 including a hard disk, etc.; and a communication section 609 including a network interface card such as a LAN card, a modem, etc. The communication section 609 performs communication processing via a network such as the Internet. A drive 610 is also connected to the input / output (I / O) interface 605 as needed. A removable medium 611, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 610 as needed so that a computer program read from it can be installed into the storage section 608 as needed.

[0219] The present invention also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiments; or may exist alone without being assembled into the device / apparatus / system. The above computer-readable storage medium carries one or more programs, and when the above one or more programs are executed, the method according to the embodiments of the present invention is implemented.

[0220] According to an embodiment of the present invention, the computer-readable storage medium may be a non-volatile computer-readable storage medium, for example, it may include but is not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the above. In the present invention, the computer-readable storage medium may be any tangible medium that contains or stores a program, and the program may be used by or in combination with an instruction execution system, apparatus, or device. For example, according to an embodiment of the present invention, the computer-readable storage medium may include the above-described ROM 602 and / or RAM 603 and / or one or more memories other than ROM 602 and RAM 603.

[0221] An embodiment of the present invention also includes a computer program product, which includes a computer program that contains program code for executing the method shown in the flowchart. When the computer program product runs in a computer system, the program code is used to cause the computer system to implement the satellite cross-domain information determination method provided by the embodiments of the present invention.

[0222] When the computer program is executed by the processor 601, the above functions defined in the system / apparatus of the embodiments of the present invention are executed. According to an embodiment of the present invention, the above-described systems, apparatuses, modules, units, etc. may be implemented by computer program modules.

[0223] In one embodiment, the computer program may rely on tangible storage media such as optical storage devices and magnetic storage devices. In another embodiment, the computer program may also be transmitted and distributed in the form of a signal on a network medium, and is downloaded and installed through the communication part 609, and / or installed from the removable medium 611. The program code included in the computer program may be transmitted by any suitable network medium, including but not limited to: wireless, wired, etc., or any suitable combination of the above.

[0224] In such an embodiment, the computer program can be downloaded and installed from a network through the communication part 609, and / or installed from the removable medium 611. When the computer program is executed by the processor 601, the above functions defined in the system of the embodiments of the present invention are executed. According to the embodiments of the present invention, the above-described systems, devices, apparatuses, modules, units, etc. can be implemented by computer program modules.

[0225] According to the embodiments of the present invention, the program code for executing the computer programs provided by the embodiments of the present invention can be written in any combination of one or more programming languages. Specifically, these computing programs can be implemented using high-level procedures and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages include, but are not limited to, such as Java, C++, Python, the "C" language, or similar programming languages. The program code can be executed entirely on the user computing device, partially on the user device, partially on a remote computing device, or entirely on a remote computing device or server. In the case of a remote computing device, the remote computing device can be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or can be connected to an external computing device (for example, by connecting through the Internet using an Internet service provider).

[0226] The flowcharts and block diagrams in the accompanying drawings illustrate the possible architectures, functions, and operations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagram may represent a module, a program segment, or a part of code, and the above module, program segment, or part of code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than marked in the accompanying drawings. For example, two consecutive blocks shown may actually be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram or flowchart, and the combination of blocks in the block diagram or flowchart, can be implemented by a dedicated hardware-based system for performing the specified functions or operations, or can be implemented by a combination of dedicated hardware and computer instructions.

[0227] Those skilled in the art can understand that the features described in the various embodiments of the present invention can be combined and / or combined in various ways, even if such combinations or combinations are not explicitly described in the present invention. In particular, without departing from the spirit and teachings of the present invention, the features described in the various embodiments of the present invention can be combined and / or combined in various ways. All such combinations and / or combinations fall within the scope of the present invention.

[0228] The embodiments of the present invention have been described above. However, these embodiments are for illustrative purposes only and are not intended to limit the scope of the present invention. Although the embodiments have been described separately above, this does not mean that the measures in each embodiment cannot be used advantageously in combination. Without departing from the scope of the present invention, those skilled in the art can make various substitutions and modifications, and all such substitutions and modifications should fall within the scope of the present invention.

Claims

1. A method for determining satellite cross-domain information, characterized in that, The method includes: Based on the regional location information of the area to be measured, determining a plurality of target sub-satellite points related to the area to be measured from among a plurality of sub-satellite points included in a search engine; Based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points, dividing the plurality of target sub-satellite points into at least one sub-satellite point sequence; Performing boundary point verification on the sub-satellite point sequence to obtain a verification result; Based on the verification results of each of the at least one sub-satellite point sequence, determining the satellite cross-domain information of the area to be measured; Wherein, the performing boundary point verification on the sub-satellite point sequence to obtain a verification result includes: Determining at least two boundary points located at the boundary positions from the sub-satellite point sequence; For each of the boundary points, based on the boundary point type determined by the boundary position of the boundary point, determining a verification point for verifying the boundary point from the satellite trajectories of the satellites related to the boundary point, wherein the satellite trajectories include a plurality of other sub-satellite points other than the target sub-satellite points related to the satellite, and the verification point is determined from among the plurality of other sub-satellite points; Based on the verification rules for the regional shape of the area to be measured and the verification point, verifying the boundary point to obtain a verification sub-result, wherein in the case where it is determined that the verification point is within the area to be measured, the verification sub-result is verification failure, and in the case where it is determined that the verification point is outside the area to be measured, the verification sub-result is verification success; Based on the verification sub-results of at least two of the boundary points, determining the verification result.

2. The method according to claim 1, characterized in that, The dividing the plurality of target sub-satellite points into at least one sub-satellite point sequence based on the timestamps of the target sub-satellite points and the satellites related to the target sub-satellite points includes: Based on the satellites related to the target sub-satellite points, dividing the plurality of target sub-satellite points to obtain at least one candidate point group, wherein the plurality of candidate sub-satellite points included in the candidate point group all belong to the same satellite; For each of the candidate point groups, sorting the plurality of candidate sub-satellite points according to their respective timestamps to obtain a candidate point sequence; Based on the time intervals between adjacent candidate sub-satellite points in the candidate point sequence, dividing the plurality of candidate sub-satellite points included in the candidate point sequence into their respective cross-domain trajectories to obtain at least one of the sub-satellite point sequences.

3. The method according to claim 2, wherein The dividing the plurality of candidate sub-satellite points included in the candidate point sequence into their respective cross-domain trajectories based on the time intervals between adjacent candidate sub-satellite points in the candidate point sequence to obtain at least one of the sub-satellite point sequences includes: Determining the i-th candidate sub-satellite point from the candidate point sequence, wherein i is a positive integer greater than 1, and in the case where the candidate sub-satellite point is the first candidate sub-satellite point in the candidate point sequence, constructing a corresponding initial sub-satellite point sequence for the first candidate sub-satellite point; Comparing the timestamp of the i-th candidate sub-satellite point with the timestamp of the (i - 1)-th candidate sub-satellite point to determine the time interval between the i-th candidate sub-satellite point and the (i - 1)-th candidate sub-satellite point; When it is determined that the time interval is less than a preset time threshold, add the i-th candidate sub-satellite point to the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point; Repeat the above operation until multiple candidate sub-satellite points are all added to their respective corresponding initial sub-satellite point sequences to obtain at least one sub-satellite point sequence.

4. The method according to claim 3, wherein The method further includes: When it is determined that the time interval is greater than or equal to the preset time threshold, construct a new initial sub-satellite point sequence different from the initial sub-satellite point sequence corresponding to the (i - 1)-th candidate sub-satellite point; Add the i-th candidate sub-satellite point to the new initial sub-satellite point sequence.

5. The method according to claim 1, wherein The determining, based on the regional position information of the area to be measured, multiple target sub-satellite points related to the area to be measured from multiple sub-satellite points included in the search engine includes: Based on the regional position information of the area to be measured, determine a query instruction for the area to be measured; Input the query instruction and the regional position information into the search engine, so that the search engine, based on the regional position information, determines at least one preset area that matches the regional position information, and based on the coordinate values of the sub-satellite points located in each preset area, determines multiple target sub-satellite points located in the area to be measured; Obtain multiple target sub-satellite points from the search engine.

6. The method according to claim 5, wherein The coordinate values of the sub-satellite points are determined by the following method: Obtain the orbital characteristic data of at least one satellite within a target time period; For each satellite, based on a preset prediction model and the orbital characteristic data, predict the satellite trajectory of the satellite within a future time period to obtain initial coordinate values respectively corresponding to multiple sub-satellite points in an inertial coordinate system; Perform an Earth coordinate system conversion on the initial coordinate values respectively corresponding to multiple sub-satellite points to obtain coordinate values respectively corresponding to multiple sub-satellite points.

7. The method according to claim 1, characterized in that, The verifying the boundary points based on the verification rule for the regional shape of the area to be measured and the verification points to obtain a verification sub-result includes: When the regional shape is circular, determine a first coordinate value of the verification point and a second coordinate value of the center of the area to be measured included in the regional position information, where the first coordinate value includes a first longitude coordinate and a first latitude coordinate, and the second coordinate value includes a second longitude coordinate and a second latitude coordinate; Respectively perform radian conversion on the first longitude coordinate, the first latitude coordinate, the second longitude coordinate, and the second latitude coordinate to obtain a first longitude radian, a first latitude radian, a second longitude radian, and a second latitude radian; Based on a first radian difference between the first longitude radian and the second longitude radian, a second radian difference between the first latitude radian and the second latitude radian, and the first latitude radian and the second latitude radian, determine a distance value between the verification point and the center of the circle; Based on a comparison result between the distance value and the radius value of the area to be measured included in the regional position, determine the verification sub-result for the boundary points.

8. The method according to claim 1, wherein Performing verification on the boundary points based on the verification rules for the regional shape of the area to be measured and the verification points to obtain a verification sub-result, including: In the case where the regional shape is a polygon or a rectangle, based on multiple regional vertices of the area to be measured in the regional position information, determining multiple vertex groups for forming the regional edges of the area to be measured, where each vertex group includes at least two of the regional vertices; For each vertex group, determining a first latitude difference between the first latitude coordinate of the verification point and the third latitude coordinate of the first regional vertex in the vertex group, and a second latitude difference between the first latitude coordinate and the fourth latitude coordinate of the second regional vertex in the vertex group; In the case where it is determined that the first latitude difference and the second latitude difference meet the preset conditions, determining a verification longitude value based on the longitude difference between the third longitude coordinate of the first regional vertex and the fourth longitude coordinate of the second regional vertex, the third latitude difference between the third latitude coordinate and the fourth latitude coordinate, the first latitude difference, and the third longitude coordinate; In the case where the verification longitude value is greater than or equal to the first longitude coordinate of the verification point, determining that there is an intersection between the ray generated with the verification point as the starting point and the regional edge formed by the vertex group; Determining the verification sub-result based on the total number of intersections between the regional edges respectively formed by multiple vertex groups and the ray generated with the verification point as the starting point.

9. The method according to claim 1, wherein Determining the satellite cross-domain information of the area to be measured based on the verification results of at least one of the sub-satellite point sequences respectively, including: For each sub-satellite point sequence, in the case where the verification result indicates that the verification of the target boundary point among at least two boundary points fails, determining a new boundary point from the satellite trajectory of the satellite to which the target boundary point belongs based on the timestamp of the target boundary point; In the case where it is determined that the verification sub-result of the new boundary point is verification success, updating the sub-satellite point sequence based on the new boundary point to obtain an updated sub-satellite point sequence; Determining the satellite cross-domain information based on the updated sub-satellite point sequence.

Images