Method for digitizing satellite coverage space data and related device

Abstract

The present disclosure provides a kind of digital processing method of satellite coverage space data and related equipment, comprising: obtaining the space data covered by satellite, and the space data is processed to obtain the processed space data;The boundary coordinate point set of the processed space data is obtained by analyzing the processed space data, and the abnormal data in the space data is determined based on the boundary coordinate point set;The abnormal data is completed to obtain the completed abnormal data, and the space data is improved based on the completed abnormal data.The present disclosure processes the space data covered by satellite, obtains the processed space data, then analyzes the processed space data to obtain the boundary coordinate point set, and then determines the abnormal data in the space data, and finally completes the abnormal data to make the final space data more accurate and real.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of satellite communication, and in particular to a digital processing method for satellite coverage space data and related equipment. BACKGROUND

[0002] With the continuous advancement of the process of digitization, China's satellite communication industry has gradually embarked on the road of digital transformation and industrial internet, and the digitization of satellite communication related resources has become the primary problem on the road of transformation and development. Among them, the characterization of the coverage range of high-orbit communication satellites is particularly important.

[0003] The characterization of high-orbit communication satellite coverage by major satellite operators is usually represented by effective isotropic radiated power (EIRP) contours. There are two forms of static orbit communication satellite coverage data sources, one is the original antenna design and simulation results, and the data format is generally grd or pat. The antenna design simulation file is parsed to generate multi-layer EIRP contour space data. The other is to obtain different EIRP contour coverage pictures from public data sources, and the data format is generally png or jpg. In the use of satellite coverage space data, the satellite coverage model usually hopes to be characterized by the satellite visible range with EIRP greater than a certain value.

[0004] In the prior art, the data generated by the two methods has certain limitations. First, if the contour generated by parsing the antenna simulation design file exceeds the maximum visible range of the static orbit satellite, the maximum visible boundary data of the space data generated by satsoft will be missing, resulting in incorrect coverage results. Second, the coordinate range and data validity of the coverage data generated by the two methods have standard deviation problems, which will cause the coverage display, spatial operation and other functions to be down. Therefore, efficient and accurate detection of satellite coverage data validity and processing of data anomalies are very important for satellite full-link simulation. SUMMARY

[0005] In view of this, the purpose of the present disclosure is to provide a digital processing method for satellite coverage space data and related equipment.

[0006] As one aspect of the present disclosure, a digital processing method for satellite coverage space data is provided, characterized by comprising:

[0007] Obtaining satellite coverage space data and processing the space data to obtain processed space data;

[0008] Parsing the boundary coordinate point set of the processed space data, and determining the abnormal data in the space data based on the boundary coordinate point set;

[0009] The abnormal data is completed to obtain completed abnormal data, and the spatial data is perfected based on the completed abnormal data.

[0010] Optionally, the spatial data covered by the satellite is acquired, and the spatial data is processed to obtain processed spatial data, including:

[0011] The spatial data is geographically registered based on a coordinate system of a world map to obtain geographically registered spatial data;

[0012] The geographically registered spatial data is subjected to a stroke processing to obtain stroke-processed spatial data;

[0013] The stroke-processed spatial data is subjected to a format processing in a preset format to obtain format-processed spatial data.

[0014] Optionally, the processed spatial data is analyzed to obtain a boundary coordinate point set, and the abnormal data in the spatial data is determined based on the boundary coordinate point set, including:

[0015] The format-processed spatial data is split into a plurality of spatial data of coverage faces;

[0016] The spatial data of the plurality of coverage faces is analyzed to obtain a boundary coordinate point set;

[0017] A distance value, an expected value and a variance value between each two adjacent boundary coordinate points on each coverage face are calculated;

[0018] In response to an absolute value of a difference between the distance value of the adjacent boundary coordinate points and the expected value being greater than or equal to a preset multiple of the variance value, the adjacent boundary coordinate points are determined as abnormal data.

[0019] Optionally, the distance, the expected value and the variance value between each two adjacent boundary coordinate points in the boundary coordinate point set are solved, and are expressed as:

[0020] The distance value between the adjacent boundary coordinate points is expressed as:

[0021] d i =dis(p i ,p i+1 )

[0022] wherein d i represents a distance between adjacent spatial coordinate points, and (p i , p i+1 ) represents two adjacent spatial coordinate points.

[0023] The expected value of the adjacent coordinate point is represented as:

[0024]

[0025] Wherein, μ represents the expected value of the adjacent coordinate point, d i represents the distance between adjacent spatial coordinate points, and n represents the number of boundary coordinate points in the boundary coordinate point set.

[0026] The variance value of the adjacent coordinate point is represented as:

[0027]

[0028] Wherein, σ represents the variance value of the adjacent coordinate point, μ represents the expected value of the adjacent coordinate point, d i represents the distance between adjacent spatial coordinate points, and n represents the number of boundary coordinate points in the boundary coordinate point set.

[0029] Optionally, the supplement processing of the abnormal data to obtain the supplement-processed abnormal data comprises:

[0030] determining the number of coverage faces of the abnormal data;

[0031] in response to determining that the coverage face of the abnormal data is one, supplement processing the coverage face of the abnormal data according to the method of generating the elevation line;

[0032] in response to determining that the coverage face of the abnormal data is two, connecting the abnormal data points on the at least two coverage faces, and supplement processing the at least two coverage faces based on the connection line.

[0033] Optionally, the supplement processing of the abnormal data to obtain the supplement-processed abnormal data comprises:

[0034] determining the search area of the coverage face of the abnormal data according to the latitude and longitude coordinates of the abnormal data;

[0035] grid processing the search area to obtain a grid vertex coordinate set;

[0036] calculating the elevation angle of each vertex in the grid vertex coordinate set;

[0037] wherein the calculation formula of the elevation angle is represented as:

[0038]

[0039] Wherein, λ E is the longitude of the ground position, φ E is the latitude of the ground position, λ sa satellite longitude;

[0040] in response to determining that the values of each two adjacent elevation angles are in turn less than 0 and greater than 0, the vertex is a component of a first part of the elevation line;

[0041] in response to determining that the values of each two adjacent elevation angles are in turn greater than 0 and less than 0, the vertex is a component of a second part of the elevation line;

[0042] completing the coverage of the abnormal data based on the component of the first part and the component of the second part.

[0043] Optionally, the coverage of the abnormal data includes a first coverage and a second coverage.

[0044] in response to determining that the coverage of the abnormal data is at least two, connecting abnormal data points on the at least two coverages, and completing the at least two coverages based on the connection, comprising:

[0045] determining a first abnormal data point and a second abnormal data point on the first coverage, and determining a third abnormal data point and a fourth abnormal data point on the second coverage;

[0046] connecting the first abnormal data point and the third abnormal data point by a first line, and connecting the second abnormal data point and the fourth abnormal data point by a second line;

[0047] in response to determining that the first line and the second line intersect, completing the first coverage and the second coverage based on the hour order of the first abnormal data point and the second abnormal data point relative to the origin by the elevation line generation method;

[0048] in response to determining that the first line and the second line do not intersect, processing a set of boundary coordinate points on the second coverage in reverse order, and completing the first coverage and the second coverage based on the elevation line generation method.

[0049] As a second aspect of the present disclosure, the present disclosure further provides a digital processing device for satellite coverage space data, characterized in that, comprising:

[0050] a data acquisition and processing module configured to acquire satellite coverage space data and process the space data to obtain processed space data;

[0051] an abnormal data determination module configured to parse the processed space data to obtain a set of boundary coordinate points, and determine abnormal data in the space data based on the set of boundary coordinate points.

[0052] an abnormal data complementing module, configured to: perform complementing processing on the abnormal data to obtain complemented abnormal data, and perfect the spatial data based on the complemented abnormal data.

[0053] As a third aspect of the present disclosure, the present disclosure also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for digital processing of satellite coverage spatial data when executing the program.

[0054] As a fourth aspect of the present disclosure, the present disclosure also provides a non-transitory computer readable storage medium storing computer instructions for causing a computer to execute the method of any one of the above.

[0055] As described above, the present disclosure provides a method for digital processing of satellite coverage spatial data and related devices. In the present disclosure, the spatial data is first processed, then the boundary coordinate point set is obtained by analyzing the processed spatial data, and then the abnormal data in the spatial data is determined through the boundary coordinate point set, and finally the abnormal data is complemented to obtain more accurate and complete spatial data. BRIEF DESCRIPTION OF DRAWINGS

[0056] In order to more clearly illustrate the technical solutions in the present disclosure or the related art, the following will briefly introduce the drawings needed to be used in the embodiments or the related art descriptions. Obviously, the drawings in the following description are only embodiments of the present disclosure, and for those skilled in the art, other drawings can also be obtained without creative labor based on these drawings.

[0057] Figure 1A A schematic diagram of a method for digital processing of satellite coverage spatial data provided by an embodiment of the present disclosure.

[0058] Figure 1B A schematic diagram of a method for processing spatial data provided by an embodiment of the present disclosure.

[0059] Figure 1C A schematic diagram of a method for determining abnormal data provided by an embodiment of the present disclosure.

[0060] Figure 1D A schematic diagram of a method for complementing abnormal data provided by an embodiment of the present disclosure.

[0061] Figure 2 A schematic diagram of a digital processing device for satellite coverage spatial data provided by an embodiment of the present disclosure.

[0062] Figure 3 An electronic device structure schematic diagram of a satellite coverage space data digitization processing method provided by an embodiment of the present disclosure. DETAILED DESCRIPTION

[0063] In order to make the objects, technical solutions and advantages of the present disclosure clearer, the present disclosure will be further described in detail below with reference to the embodiments and the accompanying drawings.

[0064] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present disclosure should be understood as the general meaning understood by those skilled in the art to which the present disclosure belongs. The terms "first", "second", and similar terms used in the embodiments of the present disclosure do not represent any order, number, or importance, but are only used to distinguish different components. The terms "include", "contain", and similar terms mean that the elements or objects before the terms cover the elements or objects listed after the terms and their equivalents, and do not exclude other elements or objects. The terms "connect" or "connected" and similar terms are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. The terms "up", "down", "left", "right", and the like only represent relative positional relationships, and when the absolute positions of the described objects change, the relative positional relationships may also change accordingly.

[0065] In the prior art, when the space data covered by the satellite is digitized, it is mostly completed by two methods of original antenna design and simulation result operation and public data source acquisition. This may result in missing boundary data in the final digitized processing result, and thus lead to inaccurate data results.

[0066] In order to solve the above problems, the present disclosure provides a satellite coverage space data digitization processing method and related equipment. Through the above method, the present disclosure first, and then, and finally.

[0067] After introducing the basic principles of the present disclosure, the various non-limiting embodiments of the present disclosure will be specifically introduced below.

[0068] Figure 1A A satellite coverage space data digitization processing method provided by an embodiment of the present disclosure.

[0069] Figure 1A The satellite coverage space data digitization processing method further includes the following steps:

[0070] Step S10: acquiring space data covered by a satellite, and processing the space data to obtain processed space data.

[0071] In some embodiments, when it is disclosed that the satellite-covered spatial data is to be processed digitally, the satellite-covered spatial data can be acquired first and pre-processed.

[0072] In some embodiments, when the satellite-covered spatial data is acquired, different acquisition methods can be used according to different sources of the spatial data. Specifically, a computer mapping application ArcGIS can be used to extract spatial data from picture information, or an antenna simulation software SATSOFT can be used to analyze and convert spatial data in the format of pat\grd antenna simulation file. It can be understood that although only the above two methods of acquiring spatial data are described in the embodiment, the acquisition of spatial data in the present application can be completed by the above two methods.

[0073] In some embodiments, the pre-processing described above can further include vectorization processing and formatting processing of the satellite-covered spatial data. It can be understood that other processing methods can also be used when the acquired spatial data is pre-processed, and the embodiment only takes vectorization processing and formatting processing as an example for description.

[0074] Figure 1B A method for processing spatial data provided by the embodiment of the present disclosure is shown in the figure.

[0075] In some embodiments, as shown in FIG. 10, step S10 is further expanded and specifically includes the following steps. Figure 1B

[0076] S101: Geographical registration of the spatial data based on the coordinate system of the world map to obtain the spatial data after geographical registration.

[0077] In some embodiments, the spatial data can be geographically registered before being vectorized. The spatial data acquired in the above process can be in the form of a picture. Since the picture is a collection of pixels, it does not have a coordinate system itself, and therefore the boundary drawn directly on the picture does not have geographical coordinates and cannot be used directly. Therefore, the embodiment can first geographically register the spatial data on the picture.

[0078] ​In some embodiments, a world map with a coordinate system of EPSG: 4326 can be found first, and then the map can be used as a reference coordinate system to complete the geographic registration by using the geographic registration function of a computer mapping tool (ArcGIS) to match the left side of the world map and the obtained spatial data. It can be understood that the reference used for geographic registration and the tool for geographic registration can also be in other forms, and the present embodiment does not make specific limitations on this.

[0079] S102: performing a line drawing process on the spatial data after the geographic registration to obtain spatial data after the line drawing process.

[0080] In some embodiments, after the spatial data is geographically registered by the present disclosure, the spatial data after the geographic registration can be vectorized. In the present embodiment, the vectorization of the spatial data after the geographic registration can specifically be a line drawing process.

[0081] In some embodiments, the line drawing process can specifically be that a vector layer is first created, then the places where the spatial data after the geographic registration and the vector layer coincide are drawn, a plurality of face layers are obtained, and finally the face layers are digitized to obtain digitized face layers. Next, fields such as Label, BEAMID, and Name can be created to associate the satellite information of the digitized face layers. Finally, the vector layer data (spatial data after the line drawing process) is exported in the format of shp (ESRI Shapefile) or GeoJSON.

[0082] S103: performing a format processing on the spatial data after the line drawing process according to a preset format to obtain spatial data after the format processing.

[0083] In some embodiments, after the spatial data is vectorized and the vector layer data (spatial data after the line drawing process) is obtained, the present disclosure can further perform a format processing on the spatial data after the line drawing process.

[0084] In some embodiments, the format of the spatial data after the vectorization can not be unique, so the format processing is needed to standardize the format of the spatial data after the line drawing process. Specifically, the vector layer data (spatial data after the line drawing process) with the format of pat or grd can be imported into a satellite antenna design software (satsoft), the equivalent isotropically radiated power (eirp) is adjusted, and finally the shp (ESRI Shapefile) format is exported and generated, which is used as the spatial data after the format processing. It can be understood that the preset format can also be in other forms, and the present embodiment does not make specific limitations on this.

[0085] In summary, in this step, the present disclosure completes the initial processing of the spatial data, so that the spatial data is embodied into spatial data with a certain format. However, the formatted spatial data obtained by the present disclosure may have missing coverage boundary data in the maximum visible boundary part, resulting in incorrect coverage results. Therefore, the present disclosure will detect this part of data and determine it as abnormal data.

[0086] Step S20: Analyzing the processed spatial data to obtain a boundary coordinate point set, and determining abnormal data in the spatial data based on the boundary coordinate point set.

[0087] In some embodiments, when determining abnormal data in the processed spatial data, the present disclosure can first analyze the processed spatial data to obtain a boundary coordinate point set of the processed spatial data, and then determine abnormal data in the processed spatial data according to the boundary coordinate point set.

[0088] Figure 1C A schematic diagram of a method for determining abnormal data provided by an embodiment of the present disclosure.

[0089] In some embodiments, as shown in FIG. 2, the present disclosure further expands step S20, which specifically includes the following steps: Figure 1C

[0090] S201: Splitting the coverage face set of the formatted spatial data to obtain spatial data of each coverage face.

[0091] In some embodiments, since the formatted spatial data derived by the satellite antenna design software (satsoft) described above defaults to stack all the equivalent isotropically radiated power (eirp) contour coverage data under one beam, the present disclosure needs to split the coverage face set in the single formatted spatial data to obtain spatial data of each equivalent isotropically radiated power (eirp) coverage face.

[0092] S202: Analyzing the spatial data of the plurality of coverage faces to obtain a boundary coordinate point set.

[0093] In some embodiments, the formatted spatial data derived by the satellite antenna design software (satsoft) can be represented in the following form:

[0094] Coverage={Coverage_eirp1, Coverage_eirp2, …, Coverage_eirp i , …, Coverage_eirp k}​

[0095] Where Coverage represents the coverage of all equivalent isotropic radiated power (eirp) of a single satellite beam, Coverage_eirp i The equivalent isotropic radiated power (eirp) in satellite single-beam coverage is represented by the value eirp. i The coverage area.

[0096] In some embodiments, after obtaining spatial data of several coverage areas in this disclosure, these data can be parsed to obtain a set of boundary coordinate points. Specifically, this can be represented as:

[0097] Counter = {p1, p2, ..., p} i , ..., p n}

[0098] Where n is the number of boundary coordinate points in the boundary coordinate point set, p i This represents the coordinates of the i-th boundary point on the boundary. The coordinates of the i-th boundary point can be represented by (lon). i lat i ) is used to represent this.

[0099] S203: Calculate the distance, expected value, and variance between every two adjacent boundary coordinate points on each coverage surface.

[0100] In some embodiments, after obtaining the set of boundary coordinate points, the present disclosure can calculate the distance, expected value, and variance between any two adjacent boundary coordinate points on each coverage surface. The distance value can be expressed as:

[0101] d i =dis(p i p i+1 )

[0102] Where, d i Represents the distance between adjacent spatial coordinate points, (p i p i+1 () represents two adjacent spatial coordinate points.

[0103] In some embodiments, the expected value can be expressed as:

[0104]

[0105] Where μ represents the expected value of adjacent coordinate points, and d i This represents the distance between adjacent spatial coordinate points, and n represents the number of boundary coordinate points in the boundary coordinate point set.

[0106] In some embodiments, the variance value can be expressed as:

[0107]

[0108] wherein σ represents a variance value of adjacent coordinate points, μ represents an expectation value of adjacent coordinate points, d i represents a distance of adjacent spatial coordinate points, and n represents a number of boundary coordinate points in the boundary coordinate point set.

[0109] S204: In response to an absolute value of a difference between the distance value of the adjacent boundary coordinate points and the expectation value being greater than or equal to a preset multiple of the variance value, determining that the adjacent boundary coordinate points are abnormal data.

[0110] In some embodiments, after obtaining the distance value, the expectation value, and the variance value between each two adjacent boundary coordinate points on each coverage, the present disclosure can determine abnormal data in the spatial data through the relationship among the three.

[0111] In some embodiments, when the absolute value of the difference between the distance value and the expectation value is greater than or equal to three times the variance value, then the two boundary coordinate points used to calculate the distance value, the expectation value, and the variance value can be determined as abnormal data. It can be understood that the aforementioned three times is the most ideal multiple determined by the inventors through careful research and numerous experiments, but in actual operation of determining abnormal data, the multiple between the absolute value of the difference between the distance value and the expectation value and the variance value is not necessarily limited to three times.

[0112] In some embodiments, the absolute value of the difference between the distance value and the expectation value being greater than or equal to three times the variance value can be expressed as:

[0113] |d i - μ | ≥ 3σ

[0114] wherein d i represents the distance value, μ represents the expectation value, and σ represents the variance value.

[0115] In summary, in this step, the spatial data of a plurality of coverages is obtained by splitting the coverage of the formatted spatial data. Then, the boundary coordinate point set is obtained by parsing the spatial data of the plurality of coverages. Finally, the final abnormal data is determined by calculating and comparing the distance value, the expectation value, and the variance value between each two adjacent boundary coordinate points on each coverage. Next, the present disclosure will complete the abnormal data to obtain the final accurate spatial data.

[0116] Step S30: Completing the abnormal data to obtain completed abnormal data, and perfecting the spatial data based on the completed abnormal data.

[0117] In some embodiments, when the present disclosure determines the abnormal data, the abnormal data can be completed, and the original spatial data can be perfected through the completed abnormal data, and finally accurate and real spatial data can be obtained.

[0118] Figure 1D A method for completing abnormal data provided by the embodiments of the present disclosure is shown in the schematic diagram.

[0119] In some embodiments, as Figure 1D shown, it is a further expansion of step S30, specifically including the following steps:

[0120] S301: Determine the number of coverage of the abnormal data.

[0121] In some embodiments, when the present disclosure needs to complete the abnormal data, the number of coverage (Coverage_eirp i ) of the abnormal data can be determined. In normal cases, the number of coverage of the abnormal data should be one (polygon1), but due to errors in machine operation, the coverage of the abnormal data may become two ({polygon1, polygon2}) or even more. Therefore, the present disclosure needs to determine the number of coverage of the abnormal data first to determine the strategy for completing the abnormal data later.

[0122] In some embodiments, no matter whether the number of coverage of the abnormal data is one or multiple, the direction of the hour hand of the abnormal data needs to be determined first. Specifically, the area calculation method and the left-hand rule of the polygon can be used to determine the direction of the hour hand of the polygon (abnormal data), wherein the positive and negative of the area of the polygon follow the left-hand rule. If the size of the area is negative, the direction of the hour hand of the polygon is counterclockwise; if the size of the area is positive, the direction of the hour hand of the polygon is clockwise.

[0123] S302: In response to determining that the coverage of the abnormal data is one, the coverage of the abnormal data is completed by the method of generating the elevation line.

[0124] In some embodiments, after the present disclosure determines the direction of the hour hand of the abnormal data, the abnormal data can be completed. If the coverage of the abnormal data is one, the abnormal data can be completed by the method of generating the elevation line.

[0125] In some embodiments, the method of generating the elevation line to complete the abnormal data described above can be specifically:

[0126] S3021: Determine the search area of the coverage of the abnormal data according to the latitude and longitude coordinates of the abnormal data.

[0127] In some embodiments, the present disclosure can first determine the search area of the abnormal data coverage according to the latitude and longitude coordinates of the abnormal data. Specifically, the clockwise direction judgment method described above can be used to obtain the clockwise sequence clockwise of the abnormal data coverage. According to the clockwise sequence clockwise and the latitude and longitude coordinates (lon p1 , lat p1 ) and (lon p2 , lat p2 ) of the abnormal data points (P1, P2) on the abnormal data coverage, the target search area can be determined.

[0128] S3022: Perform grid processing on the search area to obtain a set of grid vertex coordinates;

[0129] In some embodiments, after the present disclosure determines the target search area, the search area can be uniformly gridded.

[0130] In some embodiments, the area in the target search area can be uniformly gridded according to latitude and longitude, and a set of gridded grid vertex coordinates can be obtained, which can be specifically represented as: {(x1, y1), (x1, y2), …, (x m , y m )}.

[0131] S3023: Calculate the elevation angle of each vertex in the set of grid vertex coordinates;

[0132] In some embodiments, after the present disclosure obtains the set of grid vertex coordinates, the elevation angle of each grid vertex in the set of grid vertex coordinates can be solved,

[0133] wherein the calculation formula of the elevation angle can be represented as:

[0134]

[0135] wherein λ E is the longitude of the ground position, φ E is the latitude of the ground position, and λ s is the satellite longitude.

[0136] S3024: In response to determining that the numerical values of each two adjacent elevation angles are in turn less than 0 and greater than 0, the vertex is a component of the first part of the elevation angle line.

[0137] In some embodiments, the elevation angle value of each grid vertex can be obtained through the aforementioned calculation formula of the elevation angle, which can be represented as E i，j . Then, the present disclosure can obtain the elevation angle value of each grid vertex from lon p1 to lonp2 traverse from lat p1 to lat p2 traverse by column, traverse by row within the target search area, if the value of each two adjacent elevation angles is sequentially less than 0 and greater than 0 (E i，j <0∩E i+1，j >0=true), the grid vertex is a component of the left half of the elevation angle line.

[0138] S3025: in response to determining that the value of each two adjacent elevation angles is sequentially greater than 0 and less than 0, the vertex is a component of the second part of the elevation angle line;

[0139] In some embodiments, if the value of each two adjacent elevation angles is sequentially greater than 0 and less than 0 (E i，j >0∩E i+1，j <0=true), the grid vertex is a component of the right half of the elevation angle line.

[0140] S3026: complete the coverage of the abnormal data based on the component of the first part and the component of the second part.

[0141] In some embodiments, after the disclosure determines the left half and the right half of the grid vertex elevation angle line, the coordinate set of the abnormal data coverage edge part to be supplemented can be obtained, and then the original abnormal coverage point set can be supplemented according to the clockwise direction obtained in the foregoing, to obtain a complete coverage.

[0142] S303: in response to determining that the coverage of the abnormal data is two, connect the abnormal data points on the at least two coverages, and complete the at least two coverages based on the connection.

[0143] In some embodiments, if the coverage of the abnormal data is the aforementioned 2 or even more, the abnormal data points on the two coverages can be connected, and then it is determined whether the connection intersects. If the connection intersects, the coverage of the abnormal data can be completed based on the method of generating the elevation angle line.

[0144] In some embodiments, when the coverage of the abnormal data is two, the coverage can be further specifically a first coverage and a second coverage. At this time, S303 can further include the following process:

[0145] S3031: determine the first abnormal data point and the second abnormal data point on the first coverage, and determine the third abnormal data point and the fourth abnormal data point on the second coverage.

[0146] S3032: connecting the first abnormal data point and the third abnormal data point; and connecting the second abnormal data point and the fourth abnormal data point.

[0147] In some embodiments, the disclosure can first determine the abnormal coverage points on the two abnormal coverage surfaces, which can be specifically represented as determining the first abnormal data point and the second abnormal data point on the first abnormal coverage surface, and determining the third abnormal data point and the fourth abnormal data point on the second abnormal coverage surface. When the abnormal data points of the two abnormal coverage surfaces are obtained, the disclosure can further determine the connection relationship between the two abnormal coverage surfaces {polygon1, polygon2}2, and then generate the connection lines l1 and l2 between the abnormal points P1→P1 and P2→P2 on the two abnormal coverage surfaces.

[0148] S3033: in response to determining that the first connection line intersects the second connection line, performing completion processing on the first coverage surface and the second coverage surface by using the elevation line generation method based on the clockwise order of the first abnormal data point and the second abnormal data point relative to the origin.

[0149] In some embodiments, if it is found that l1 and l2 intersect, the clockwise direction of P1→P2 relative to the latitude and longitude coordinates (0, 0) is determined by using the clockwise direction determination algorithm in the above steps, and the clockwise direction clockwise is obtained. Then, the coverage edge part between P2→P1 and P2→P1 is supplemented in the clockwise direction clockwise by using the elevation line generation method in the above steps, and a complete coverage surface is obtained.

[0150] S3034: in response to determining that the first connection line does not intersect the second connection line, performing reverse order processing on the boundary coordinate point set on the second coverage surface, and performing completion processing on the first coverage surface and the second coverage surface according to the elevation line generation method.

[0151] In some embodiments, if the connection lines l1 and l2 do not intersect, the boundary point set of the second coverage surface is sorted in reverse order, and then the above operation (S3033) is performed to obtain a complete coverage surface. It can be understood that the boundary point set of the first coverage surface can also be processed in reverse order, and then the above steps are continued. In this embodiment, only the boundary point set of the second coverage surface is processed in reverse order as an example.

[0152] In some embodiments, after the disclosure completes the coverage surface of the abnormal data, the disclosure can also verify the validity of the completed spatial data. First, each coverage surface polygon i Each two coordinate points P i →P i+1The formed line segments construct a polygon edge set, and then the Bentley-Ottmann algorithm is used to solve the intersection points of the edge set. If the intersection points are found by traversing the edge set, the covering surface has self-intersection, and the shorter path between the two line segment paths of the intersection needs to be removed, and one end point of one of the intersection line segments is removed, and the remaining part is connected in order, that is, the non-self-intersecting covering surface is obtained.

[0153] In some embodiments, the present disclosure can also standardize the data geographical coordinates. Since the processed spatial data sources are not the same, the coordinate ranges of the spatial data coverage also exist in two forms, that is, the coordinate range of 0° to 360° and the coordinate range of -180° to 180°.

[0154] In some embodiments, the coverage of different coordinate ranges can be standardized to between -180° and 180°. Specifically, the 180° meridian can be used to cut each beam coverage. If the cut coverage has a coordinate point greater than 180 degrees, the cut coverage coordinates are translated westward by 360 degrees to obtain the standardized coverage.

[0155] In summary, the present disclosure first processes the spatial data of satellite coverage to obtain the formatted spatial data, then analyzes the formatted spatial data to obtain the boundary coordinate point set, then calculates and compares the distance value, expectation and variance value between each two adjacent boundary coordinate points on each coverage to determine the abnormal data, and finally completes the abnormal data by the method of elevation line generation to obtain the final accurate and real spatial data.

[0156] Based on the same technical concept, the present disclosure also provides a satellite coverage spatial data digitization processing device corresponding to the method of any of the above embodiments. Through the satellite coverage spatial data digitization processing device provided by the present disclosure, the satellite coverage spatial data digitization processing method of any of the above embodiments can be realized.

[0157] Figure 2 The satellite coverage spatial data digitization processing device provided by the present disclosure has the structure shown in the schematic diagram.

[0158] Figure 2 The satellite coverage spatial data digitization processing device further includes the following modules:

[0159] The data acquisition and processing module 10, the abnormal data determination module 20 and the abnormal data completion module 30;

[0160] The data acquisition and processing module 10 is configured to acquire spatial data covered by satellites and process the spatial data to obtain processed spatial data. The following steps are specifically performed:

[0161] The spatial data is geographically registered based on a coordinate system of a world map to obtain geographically registered spatial data.

[0162] The geographically registered spatial data is subjected to a stroke processing to obtain stroke-processed spatial data.

[0163] The stroke-processed spatial data is subjected to a format processing in a preset format to obtain format-processed spatial data.

[0164] The abnormal data determination module 20 is configured to parse the processed spatial data to obtain a boundary coordinate point set and determine abnormal data in the spatial data based on the boundary coordinate point set. The following steps are specifically performed:

[0165] The coverage of the format-processed spatial data is split to obtain spatial data of a plurality of coverages.

[0166] The spatial data of the plurality of coverages is parsed to obtain a boundary coordinate point set.

[0167] A distance value, an expected value and a variance value between each two adjacent boundary coordinate points on each coverage are calculated.

[0168] In response to an absolute value of a difference between the distance value of the adjacent boundary coordinate points and the expected value being greater than or equal to a preset multiple of the variance value, the adjacent boundary coordinate points are determined as abnormal data.

[0169] The distance, the expected value and the variance value between each two adjacent boundary coordinate points in the boundary coordinate point set are solved, which is expressed as:

[0170] The distance value between the adjacent boundary coordinate points is expressed as:

[0171] d i =dis(p i ,p i+1 )

[0172] wherein d i represents a distance between adjacent spatial coordinate points, and (p i , p i+1 ) represents two adjacent spatial coordinate points.

[0173] The expected value of the adjacent coordinate points is expressed as:

[0174]

[0175] wherein μ represents the expected value of adjacent coordinate points, d i represents the distance between adjacent spatial coordinate points, and n represents the number of boundary coordinate points in the boundary coordinate point set.

[0176] The variance value of the adjacent coordinate points is represented as:

[0177]

[0178] wherein σ represents the variance value of adjacent coordinate points, μ represents the expected value of adjacent coordinate points, d i represents the distance between adjacent spatial coordinate points, and n represents the number of boundary coordinate points in the boundary coordinate point set.

[0179] The abnormal data completion module 30 is configured to perform a completion process on the abnormal data to obtain completed abnormal data, and to perfect the spatial data based on the completed abnormal data. The following steps are specifically performed:

[0180] determining the number of coverage areas of the abnormal data;

[0181] in response to determining that the coverage area of the abnormal data is one, performing a completion process on the coverage area of the abnormal data according to an elevation line generation method;

[0182] in response to determining that the coverage area of the abnormal data is two, connecting abnormal data points on the at least two coverage areas, and performing a completion process on the at least two coverage areas based on the connection.

[0183] The response to determining that the coverage area of the abnormal data is one, performing a completion process on the coverage area of the abnormal data according to an elevation line generation method, includes:

[0184] determining a search area of the coverage area of the abnormal data according to the latitude and longitude coordinates of the abnormal data;

[0185] performing a grid processing on the search area to obtain a grid vertex coordinate set;

[0186] calculating the elevation angle of each vertex in the grid vertex coordinate set;

[0187] wherein the calculation formula of the elevation angle is represented as:

[0188]

[0189] wherein λ E is the longitude of the ground position, φ E is the latitude of the ground position, and λ s is the satellite longitude.

[0190] in response to determining that the values of each two adjacent elevation angles are in turn less than 0 and greater than 0, the vertex is a component of the first part of the elevation line;

[0191] in response to determining that the values of each two adjacent elevation angles are in turn greater than 0 and less than 0, the vertex is a component of the second part of the elevation line;

[0192] complement the coverage of the abnormal data based on the component of the first part and the component of the second part.

[0193] the coverage of the abnormal data includes a first coverage and a second coverage;

[0194] in response to determining that the coverage of the abnormal data is at least two, connecting abnormal data points on the at least two coverages, and complementing the at least two coverages based on the connection, comprising:

[0195] determining a first abnormal data point and a second abnormal data point on the first coverage, and determining a third abnormal data point and a fourth abnormal data point on the second coverage;

[0196] connecting the first abnormal data point and the third abnormal data point by a first line, and connecting the second abnormal data point and the fourth abnormal data point by a second line;

[0197] in response to determining that the first line and the second line intersect, complementing the first coverage and the second coverage by the elevation line generation method based on the hour order of the first abnormal data point and the second abnormal data point relative to the origin;

[0198] in response to determining that the first line and the second line do not intersect, performing reverse order processing on a boundary coordinate point set on the second coverage, and complementing the first coverage and the second coverage according to the elevation line generation method.

[0199] Based on the same technical concept, the present disclosure also provides an electronic device corresponding to the above-mentioned any embodiment method, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to realize the satellite coverage space data digitization processing method of any one of the above embodiments.

[0200] Figure 3A more specific electronic device hardware structure schematic diagram provided by the embodiment is shown, and the device can include: a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. The processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040 are connected to each other through the bus 1050 for internal communication.

[0201] The processor 1010 can be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an ASIC (Application Specific Integrated Circuit), or one or more integrated circuits, and is used to execute related programs to implement the technical solutions provided by the embodiments of the present specification.

[0202] The memory 1020 can be implemented by a ROM (Read Only Memory), a RAM (Random Access Memory), a static storage device, a dynamic storage device, etc. The memory 1020 can store an operating system and other application programs. When the technical solutions provided by the embodiments of the present specification are implemented by software or firmware, the related program codes are stored in the memory 1020 and are called and executed by the processor 1010.

[0203] The input / output interface 1030 is used to connect input / output modules to realize information input and output. The input / output modules can be configured as components in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. The input device can include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device can include a display, a speaker, a vibrator, an indicator light, etc.

[0204] The communication interface 1040 is used to connect a communication module (not shown in the figure) to realize the communication interaction between the device and other devices. The communication module can realize communication through a wired manner (such as USB, network cable, etc.) or a wireless manner (such as mobile network, WIFI, Bluetooth, etc.).

[0205] The bus 1050 includes a channel to transmit information between various components (such as the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040) of the device.

[0206] It should be noted that although the above device only shows the processor 1010, the memory 1020, the input / output interface 1030, the communication interface 1040 and the bus 1050, but in the specific implementation process, the device can also include other components necessary for normal operation. In addition, those skilled in the art can understand that the above device can also only contain the components necessary to implement the embodiments of the present disclosure, and does not necessarily contain all the components shown in the figure.

[0207] The electronic device of the above embodiment is used to implement the digitization processing method of the satellite coverage space data in any of the preceding embodiments, and has the beneficial effects of the corresponding method embodiments, which are not described here.

[0208] Based on the same technical concept, the present disclosure also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to perform the digitization processing method of the satellite coverage space data according to any of the above embodiments.

[0209] The computer-readable medium of the present embodiment includes permanent and non-permanent, removable and non-removable media, which can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0210] The computer instructions stored in the storage medium of the above embodiment are used to cause the computer to perform the digitization processing method of the satellite coverage space data according to any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which are not described here.

[0211] Those skilled in the art should understand that the above discussion of any of the embodiments is only exemplary and is not intended to imply that the scope (including claims) of the present disclosure is limited to these examples; under the idea of the present disclosure, the above embodiments or technical features between different embodiments can also be combined, the steps can be implemented in any order, and there are many other changes of different aspects of the embodiments of the present disclosure as described above. In order to be brief, they are not provided in detail.

[0212] In addition, to simplify the description and discussion, and so as not to make the embodiments of the disclosure difficult to understand, the well-known power / ground connections to integrated circuit (IC) chips and other components can or can not be shown in the provided drawings. Further, devices can be shown in block diagram form in order to avoid making the embodiments of the disclosure difficult to understand, and this also takes into account the fact that details regarding implementation of these block diagram devices are highly dependent on the platform to which the embodiments of the disclosure are to be implemented (i.e., these details should be well within the understanding of one of skill in the art). Where specific details (e.g., circuitry) are set forth in order to describe an illustrative embodiment of the disclosure, it will be apparent to one of skill in the art that the embodiments of the disclosure can be practiced without, or with variation of, these specific details. Thus, the description is to be considered as illustrative and not restrictive, and the scope of the disclosure to be indicated by the appended claims.

[0213] Although the disclosure has been described in conjunction with specific embodiments thereof, numerous alternatives, modifications, and variations will be readily apparent to those of ordinary skill in the art in light of the foregoing descriptions. For example, other memory architectures (e.g., dynamic RAM (DRAM)) can use the embodiments discussed.

[0214] The embodiments of the disclosure are intended to cover all such alternatives, modifications, and variations as come within the scope of the broadest possible interpretation of the appended claims. Accordingly, any and all such modifications, variations, and equivalents that fall within the spirit and scope of the embodiments of the disclosure are intended to be included herein.

Claims

1. A method for digital processing of satellite coverage space data, characterized in that, The method comprises the following steps: acquiring spatial data covered by a satellite, and processing the spatial data to obtain processed spatial data; analyzing the processed spatial data to obtain a boundary coordinate point set thereof, and determining abnormal data in the spatial data based on the boundary coordinate point set; judging the number of coverage faces of the abnormal data; in response to determining that the coverage face of the abnormal data is one, determining a search area of the coverage face of the abnormal data according to the latitude and longitude coordinates of the abnormal data; performing grid processing on the search area to obtain a grid vertex coordinate set; calculating the elevation angle of each vertex in the grid vertex coordinate set; in response to determining that the numerical values of every two adjacent elevation angles are sequentially smaller than 0 and larger than 0, the vertex is a component of a first part of an elevation angle line; in response to determining that the numerical values of every two adjacent elevation angles are sequentially larger than 0 and smaller than 0, the vertex is a component of a second part of the elevation angle line; completing the coverage face of the abnormal data based on the components of the first part and the components of the second part; in response to determining that the coverage face of the abnormal data is two, the coverage face of the abnormal data further comprises a first coverage face and a second coverage face; determining a first abnormal data point and a second abnormal data point on the first coverage face, and determining a third abnormal data point and a fourth abnormal data point on the second coverage face; performing first connection on the first abnormal data point and the third abnormal data point, and performing second connection on the second abnormal data point and the fourth abnormal data point; in response to determining that the first connection and the second connection intersect, completing the first coverage face and the second coverage face based on the hour order of the first abnormal data point and the second abnormal data point relative to the origin by using the elevation angle line generation method; in response to determining that the first connection and the second connection do not intersect, performing reverse processing on a boundary coordinate point set on the second coverage face, and completing the first coverage face and the second coverage face based on the elevation angle line generation method; completing the abnormal data to obtain completed abnormal data, and perfecting the spatial data based on the completed abnormal data.

2. The method of claim 1, wherein, The method comprises the following steps: geographically registering the spatial data based on a coordinate system of a world map to obtain geographically registered spatial data; performing edge processing on the geographically registered spatial data to obtain edge-processed spatial data; performing format processing on the edge-processed spatial data in a preset format to obtain format-processed spatial data.

3. The method of claim 2, wherein, The method comprises the following steps: splitting a coverage face set of the format-processed spatial data to obtain spatial data of a plurality of coverage faces; analyzing the spatial data of the plurality of coverage faces to obtain a boundary coordinate point set; Calculate the distance value, the expected value and the variance value between each two adjacent boundary coordinate points on each coverage; In response to the absolute value of the difference between the distance value of the adjacent boundary coordinate points and the expected value being greater than or equal to the preset multiple of the variance value, determine that the adjacent boundary coordinate points are abnormal data.

4. The method of claim 3, wherein, The distance, the expected value and the variance value between each two adjacent boundary coordinate points in the boundary coordinate point set are calculated, and are represented as: The distance value between the adjacent boundary coordinate points is represented as: wherein, denotes the distance between adjacent boundary coordinate points, denotes two adjacent boundary coordinate points; The expected value of the adjacent boundary coordinate points is represented as: wherein, represents an expected value of the adjacent boundary coordinate point, represents a distance of the adjacent boundary coordinate point, represents a number of the boundary coordinate points in the set of boundary coordinate points; The variance value of the adjacent boundary coordinate points is represented as: wherein, denotes a variance value of the adjacent boundary coordinate points, denotes an expectation value of the adjacent boundary coordinate points, denotes a distance of the adjacent boundary coordinate points, denotes a number of the boundary coordinate points in the set of the boundary coordinate points.

5. The method of claim 1, wherein, The calculation formula of the elevation angle is represented as: wherein, is the longitude of the ground position, is the latitude of the ground position, is the satellite longitude.

6. A device for digital processing of satellite coverage space data, characterized in that, Comprise: A data acquisition and processing module configured to acquire spatial data covered by a satellite and process the spatial data to obtain processed spatial data; An abnormal data determination module configured to parse the processed spatial data to obtain a boundary coordinate point set thereof, and determine abnormal data in the spatial data based on the boundary coordinate point set; Determine the number of coverage faces of the abnormal data; In response to determining that the coverage face of the abnormal data is one, determine a search area of the coverage face of the abnormal data according to the latitude and longitude coordinates of the abnormal data; Perform grid processing on the search area to obtain a grid vertex coordinate set; Calculate the elevation angle of each vertex in the grid vertex coordinate set; In response to determining that the numerical values of each two adjacent elevation angles are sequentially less than 0 and greater than 0, the vertex is a component of a first part of an elevation angle line; In response to determining that the numerical values of each two adjacent elevation angles are sequentially greater than 0 and less than 0, the vertex is a component of a second part of the elevation angle line; Perform completion processing on the coverage face of the abnormal data based on the components of the first part and the second part; In response to determining that the coverage face of the abnormal data is two, the coverage face of the abnormal data further comprises a first coverage face and a second coverage face; Determine a first abnormal data point and a second abnormal data point on the first coverage face, and determine a third abnormal data point and a fourth abnormal data point on the second coverage face; Perform first connection on the first abnormal data point and the third abnormal data point, and perform second connection on the second abnormal data point and the fourth abnormal data point; In response to determining that the first connection and the second connection intersect, perform completion processing on the first coverage face and the second coverage face by the elevation angle line generation method based on the hour clock order of the first abnormal data point and the second abnormal data point relative to the origin; In response to determining that the first connection and the second connection do not intersect, perform reverse order processing on a boundary coordinate point set on the second coverage face, and perform completion processing on the first coverage face and the second coverage face according to the elevation angle line generation method; An abnormal data completion module configured to perform completion processing on the abnormal data to obtain completed abnormal data, and perfect the spatial data based on the completed abnormal data.

7. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the method of any one of claims 1 to 5 when executing the program. The processor implements the method of any one of claims 1 to 5 when executing the program.

8. A non-transitory computer-readable storage medium storing computer instructions for causing a computer to perform the method of any one of claims 1 to 5.

Images