Radar beam real-time distribution method and system based on geographic information system
By acquiring the scanning angle of the radar beam and satellite orbit parameters, and combining them with the GIS raster database, the allocation of space-based radar beams was optimized, solving the problem of insufficient utilization of space-based radar resources, achieving efficient detection and tracking of key targets, and saving processing resources.
Patent Information
- Application Number
- CN202511070146.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-07-31
- Publication Date
- 2025-11-07
AI Technical Summary
In existing technologies, the beam allocation method of space-based radar fails to make effective use of limited processing resources, resulting in ineffective detection of non-target areas and making it difficult to achieve efficient detection and tracking of key targets.
By acquiring the scanning angle and satellite orbit parameters of the radar beam, and combining them with the GIS raster database, the land and sea attributes of the beam center point, edge feature points, and inner circle feature points are determined, and the beam allocation strategy is optimized to avoid non-target areas.
It achieves efficient use of processing resources with low computational complexity, avoids detection of non-target areas such as land and islands, improves the detection and tracking capabilities of key targets, and saves on-board processing resources.
Smart Images
Figure CN120908758A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of radar resource scheduling, in particular to a radar beam real-time allocation method and system based on a geographic information system and a storage medium. BACKGROUND
[0002] Radar beam allocation is crucial to the efficient operation of the system, especially in space-based radars, where processing resources are limited. A reasonable beam allocation method can maximize the use of limited processing resources and improve the detection and tracking capabilities of key targets. By optimizing the beam allocation method, the radar can cover a wider detection area and save more processing resources for searching and tracking key targets.
[0003] The traditional beam allocation method is to arrange all beams in a scan period according to the radar's wave table and the predetermined beam pointing direction. This allocation method can ensure coverage of the detection area. For space-based radar sea detection, a reasonable beam allocation strategy is needed to avoid ineffective detection of non-target areas such as land, islands, and reefs under limited on-board processing resources, and to achieve detection and tracking of key targets.
[0004] Therefore, how to implement a reasonable and low-implementation-difficulty beam allocation method is the current technical difficulty in the field. SUMMARY
[0005] The purpose of the present application is to solve the problems existing in the prior art and provide a radar beam real-time allocation method based on a geographic information system, comprising the following steps: S1: Obtain the beam scanning angle according to the radar wave table, and obtain the beam edge feature point scanning angle according to the beam scanning angle, the azimuth beam width, and the elevation beam width; S2: Calculate the coordinates of the beam center point and the latitude and longitude coordinates of the beam edge feature points based on satellite orbit parameters, and obtain the inner circle feature point coordinates according to the beam center point coordinates and the latitude and longitude coordinates of the edge feature points; S3: Determine the sea-land properties of the beam center point coordinates, the latitude and longitude coordinates of the edge feature points, and the inner circle feature point coordinates through the geographic information system (GIS) grid database, and determine whether to perform current beam scheduling according to the sea-land properties.
[0006] Preferably, in step S1, the beam edge feature point scanning angle is obtained according to the beam scanning angle, the azimuth beam width, and the elevation beam width, further comprising: According to the beam scanning angle and the azimuth beam width the azimuth feature point scanning angle corresponding to the two azimuth edge feature points of the azimuth beam edge is 、 ; According to the beam scanning angle and the beam width in the elevation direction The scanning angle of the two feature points on the beam edge in the elevation direction is 、 .
[0007] Preferably, in step S2, based on the satellite orbit parameters, the coordinates of the beam center point and the latitude and longitude coordinates of the beam edge feature points are calculated, further comprising: S211: According to the distance corresponding to the beam center point and the beam scanning angle, the antenna local position coordinates of the beam center point in the antenna local coordinate system are estimated, as follows: wherein, is the distance corresponding to the beam center point, is the beam scanning angle at time t; S212: According to the antenna local position coordinates and the satellite orbit parameters, the geocentric position coordinates of the beam center point in the geocentric fixed system are calculated; S213: According to the geocentric position coordinates, the coordinates of the beam center point are obtained, and the latitude and longitude coordinates of the beam edge feature points are calculated according to the distance and scanning angle corresponding to the beam edge feature points.
[0008] Preferably, in step S212, according to the antenna local position coordinates and the satellite orbit parameters, the geocentric position coordinates of the beam center point in the geocentric fixed system are calculated, further comprising: According to the satellite orbit parameters including the satellite platform orbit height, the earth radius, the ascending node right ascension and the perigee amplitude, and the antenna local position coordinates, the geocentric position coordinates of the beam center point in the geocentric fixed system are calculated, as follows: wherein, is the current time, is the satellite platform orbit height, is the earth radius, represents the angle between the prime meridian and the equinox axis at time t, represents the ascending node right ascension, represents the perigee amplitude, denotes the inclination of the satellite orbit, denotes the argument of latitude of the satellite at the time of the perigee, i.e. the true anomaly.
[0009] Preferably, in step S213, obtaining the beam center point coordinate according to the geocentric position coordinate further comprises: obtaining the beam center point coordinate according to the geocentric position coordinate by a calculation formula as follows: wherein, denotes the beam center point longitude, denotes the beam center point latitude, arcsin(.) is the inverse sine function, cos(.) is the inverse cosine function.
[0010] Preferably, in step S2, obtaining the inner circle feature point coordinate according to the beam center point coordinate and the edge feature point longitude and latitude coordinate further comprises: obtaining the inner circle feature point coordinate according to the beam center point coordinate and the longitude and latitude coordinate of the beam edge feature point . as follows: limiting the longitude of the inner circle feature point : wherein, denotes the longitude of the inner circle feature point coordinate, denotes the latitude of the inner circle feature point coordinate, denotes the longitude of the coordinate of the beam edge feature point , denotes the latitude of the coordinate of the beam edge feature point . wherein, obtaining the longitude and latitude coordinate of the inner circle feature point other than the inner circle feature point coordinate according to the beam center point coordinate and the longitude and latitude coordinate of the edge feature point.
[0011] Preferably, in step S3, judging the sea-land attribute of the beam center point coordinate, the edge feature point longitude and latitude coordinate and the inner circle feature point coordinate by a geographic information system (GIS) raster database further comprises: S31: grid division is performed on the earth's surface, geographic coordinate conversion algorithm is used to convert the map data into raster data, and a Morton code is generated for the raster data; S32: the geographic coordinates of the beam center point coordinates, the edge feature point longitude and latitude coordinates, and the inner circle feature point coordinates are converted into coordinate row and column numbers in the raster data, and the coordinate row and column numbers are converted into Mortan codes by using the encoding rules of the Mortan code; S33: the raster data corresponding to the Mortan code is retrieved, and it is determined whether the position is a land attribute or a sea attribute.
[0012] Preferably, in step S3, it is judged whether to perform current beam scheduling according to the sea-land attribute, further comprising: If the beam center point coordinates are land, more than half of the edge feature point longitude and latitude coordinates are land, and two or more of the inner circle feature point coordinates are land, the current beam scheduling is not performed; otherwise, the current beam is normally allocated.
[0013] Based on the same concept, the application also provides a radar beam real-time allocation system based on a geographic information system, comprising: An angle calculation module obtains a beam scanning angle according to a radar beam table, and obtains a beam edge feature point scanning angle according to the beam scanning angle, an azimuth beam width, and a pitch beam width; A coordinate conversion module calculates a beam center point coordinate and a beam edge feature point longitude and latitude coordinate based on satellite orbit parameters, and obtains an inner circle feature point coordinate according to the beam center point coordinate and the edge feature point longitude and latitude coordinate; A judgment and decision module judges a sea-land attribute of the beam center point coordinate, the edge feature point longitude and latitude coordinate, and the inner circle feature point coordinate through a geographic information system (GIS) grid database, and judges whether to perform current beam scheduling according to the sea-land attribute.
[0014] Based on the same concept, the application also provides a computer readable storage medium, which stores computer codes, when the computer codes are executed, the steps of the radar beam real-time allocation method based on the geographic information system are executed.
[0015] Compared with the prior art, the application has the following beneficial effects: In the application, the beam edge feature point scanning angle is obtained according to the beam scanning angle, the azimuth beam width, and the pitch beam width, the beam center point coordinate and the beam edge feature point longitude and latitude coordinate are calculated based on the satellite orbit parameters, the inner circle feature point coordinate is obtained according to the beam center point coordinate and the edge feature point longitude and latitude coordinate, and the calculation complexity is low.
[0016] The application determines the sea-land attribute of the beam center point coordinate, the latitude and longitude coordinate of the edge feature point and the latitude and longitude coordinate of the inner circle feature point through the geographic information system (GIS) grid database, and reasonably avoids the search and detection of the non-target area by using the GIS information, so that the on-board processing resources can be effectively saved; and the application determines whether to perform the current beam scheduling according to the sea-land attribute, is convenient to realize and has high engineering application value. BRIEF DESCRIPTION OF DRAWINGS
[0017] Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments with reference made to the accompanying drawings. The drawings are for purposes of illustration only and are not intended to be limiting of the application.
[0018] Figure 1 A flow chart of a radar beam real-time distribution method based on a geographic information system according to the application; Figure 2 A schematic diagram of a beam center and feature points according to the application. DETAILED DESCRIPTION
[0019] In order to make the objects, technical solutions and advantages of the application clearer, the application will be further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the application and not to limit the application. Obviously, the described embodiments are part of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of the present application.
[0020] Those skilled in the art can understand that, unless specifically stated, the singular form "a", "an", "said" used herein can also include the plural form. It should be further understood that the phrase "comprising" used in the specification of the application means that the features, integers, steps, operations, elements and / or components exist, but does not exclude the existence or addition of one or more other features, integers, steps, operations, elements, components and / or groups thereof.
[0021] First embodiment Please refer to Figure 1 and Figure 2 As shown in the drawings, the embodiment provides a radar beam real-time distribution method based on a geographic information system, which is suitable for beam distribution of space-based radar sea detection, and includes the following steps: S1: Obtain the beam scanning angle according to the radar position table, and obtain the beam edge feature point scanning angle according to the beam scanning angle, azimuth beamwidth, and elevation beamwidth. Specifically, in this embodiment, the current beam scanning angle is obtained according to the radar position table. The beam edge feature points are A1 / A2 / E1 / E2.
[0022] Preferably, in step S1, obtaining the beam edge feature point scanning angle based on the beam scanning angle, azimuth beamwidth, and pitch beamwidth further includes: According to the beam scanning angle and azimuth beamwidth The scanning angle of the azimuth feature points corresponding to the two feature points on the edge of the azimuth beam is obtained as follows: , Specifically, in this embodiment, The scanning angle corresponding to feature point A1 on the edge of the azimuth beam. Let A2 be the scanning angle corresponding to the azimuth beam edge feature point. Given that the current azimuth beamwidth is 2° and the elevation beamwidth is 0.7°, then the scanning angles corresponding to the azimuth beam edge feature points A1 and A2 are: , ; According to the beam scanning angle and pitch beamwidth The elevation feature point scanning angle corresponding to the two feature points on the elevation beam edge is obtained as follows: , Specifically, in this embodiment, The scanning angle corresponding to the two feature points E1 at the edge of the beam is the pitch angle. The scanning angle corresponding to feature point E2 at the edge of the pitch beam is [value missing]. The scanning angles corresponding to feature points E1 and E2 at the edge of the pitch beam are [value missing]. , .
[0023] Please see Figure 2 As shown, S2: Based on the satellite orbit parameters, calculate the coordinates of the beam center point and the latitude and longitude coordinates of the beam edge feature points. Based on the coordinates of the beam center point and the latitude and longitude coordinates of the edge feature points, obtain the coordinates of the inner circle feature points. Specifically, in this embodiment, based on the satellite orbit parameters, calculate the coordinates of the beam center point O and the latitude and longitude information of the four feature points A1 / A2 / E1 / E2 on the beam edge, and calculate the latitude and longitude information of the center positions A3 / A4 / E3 / E4 (denoted as the four inner circle feature points) of the beam center point and the four feature points on the beam edge.
[0024] Preferably, in step S2, calculating the coordinates of the beam center point and the latitude and longitude coordinates of the beam edge feature points based on the satellite orbit parameters further includes: S211: Estimate the local antenna position coordinates of the beam center point in the antenna local coordinate system based on the distance corresponding to the beam center point and the beam scanning angle, as shown below: in, The distance corresponding to the beam center point. Let be the beam scanning angle at time t. Specifically, in this embodiment, ; S212: Calculate the geocentric coordinates of the beam center point under geocentric fixed connection based on the antenna local position coordinates and satellite orbit parameters; S213: Obtain the beam center point coordinates based on the geocentric position coordinates, and calculate the latitude and longitude coordinates of the beam edge feature points based on the distance and scanning angle corresponding to the beam edge feature points.
[0025] Preferably, in step S212, calculating the geocentric coordinates of the beam center point under geocentric fixed-axis configuration based on the antenna local position coordinates and satellite orbit parameters further includes: Based on satellite orbital parameters including satellite platform orbital altitude, Earth radius, right ascension of the ascending node, and argument of perigee, as well as the local position coordinates of the antenna, the geocentric coordinates of the beam center point in geocentric fixed coordinates are calculated as follows: in, For the current moment, The orbital altitude of the satellite platform. For the Earth's radius, express The angle between the prime meridian and the vernal equinox axis at any given moment. Indicates the right ascension of the ascending node. Indicates the argument of perigee. Indicates the satellite's orbital inclination. Indicates the satellite's motion from perigee. The angle after a given moment, i.e., the true perimeter angle, specifically in this embodiment... =1.9054°, =500km, =6371km, =0.1253°, = 20°, = 30°, = 44°.
[0026] Preferably, in step S213, the beam center point coordinates are obtained according to the geocentric position coordinates, further comprising: obtaining the beam center point coordinates according to the geocentric position coordinates The calculation formula is as follows: wherein, represents the beam center point longitude (east longitude is positive), represents the target latitude (north latitude is positive), arcsin(.) is the inverse sine function, cos(.) is the inverse cosine function. Specifically, in the embodiment, the beam center point coordinates are , More preferably, the coordinates of the beam edge feature point A1 in the antenna local coordinate system are estimated according to the distance corresponding to the beam edge feature point A1 and the scanning angle, as follows: wherein, is the distance corresponding to the beam center point, is the slant range difference within the main lobe of the beam, is the beam scanning angle at t, is the azimuth beam width; The coordinates of A1 in the geocentric coordinate system are calculated according to the antenna local position coordinates of the beam edge feature point A1, as follows: The longitude and latitude coordinates of the beam edge feature point A1 are obtained according to the coordinates of the beam edge feature point A1 in the geocentric coordinate system , as follows: wherein, represents the longitude of the beam edge feature point A1, represents the latitude of the beam edge feature point A1, arcsin(.) is the inverse sine function, cos(.) is the inverse cosine function. Specifically, in the embodiment, the longitude and latitude coordinates of the beam edge feature point A1 are .
[0027] Similarly, the latitude and longitude coordinates of beam edge feature points A2, E1, and E2 can be obtained. Specifically, in this embodiment, the latitude and longitude coordinates of other beam edge feature points can be obtained as A2. E1 E2 .
[0028] Preferably, in step S2, obtaining the coordinates of the inner circle feature points based on the coordinates of the beam center point and the latitude and longitude coordinates of the edge feature points further includes: Based on the beam center point coordinates and the beam edge feature points The inner circle feature points are obtained from the latitude and longitude coordinates. coordinate As shown below: Then determine the range of longitude values. limit: in, Indicates inner circle feature points Coordinates (longitude) Indicates inner circle feature points coordinates latitude, Indicates the beam edge feature point longitude coordinates Indicates the beam edge feature point The coordinates of latitude; Among them, the inner circle feature points are obtained based on the coordinates of the beam center point and the latitude and longitude coordinates of the edge feature points. The latitude and longitude coordinates of the inner circle feature points are calculated. Specifically, in this embodiment, the latitude and longitude information of the center positions of the beam center point and the four feature points along the beam edge (denoted as the four inner circle feature points A3 / A4 / E3 / E4) are calculated. The latitude and longitude of the inner circle feature point A3 is... A4 E3 E4 .
[0029] S3: Determine the land / sea attributes of the beam center point coordinates, edge feature point latitude and longitude coordinates, and inner circle feature point coordinates through the Geographic Information System (GIS) raster database. Based on the land / sea attributes, determine whether to perform current beam scheduling. Specifically, in this embodiment, GIS technology is used to obtain information on the beam center O, the four beam edge feature points A1 / A2 / E1 / E2, and the four inner circle feature points A3 / A4 / E3 / E4, which are respectively land / sea. Based on the land / sea information of the nine points, determine whether to perform current beam scheduling.
[0030] Preferably, in step S3, the sea-land attribute of the beam center point coordinate, the edge feature point longitude and latitude coordinate and the inner circle feature point coordinate is determined by the geographic information system GIS grid database, and further comprising: S31: the earth surface is meshed, the map data is converted into grid data by using the geographic coordinate conversion algorithm, and the Morton code is generated for the grid data to facilitate efficient storage and retrieval; S32: the geographic coordinates of the beam center point coordinate, the edge feature point longitude and latitude coordinate and the inner circle feature point coordinate are converted into the coordinate row and column number in the grid data, and the coordinate row and column number is converted into the Mortan code by using the encoding rule of the Morton code. Specifically, in the embodiment, the geographic coordinates of the four beam edge feature points and the four inner circle feature points are converted into the row and column number in the grid data; S33: the grid data corresponding to the Mortan code is retrieved to determine whether the position is land attribute or sea attribute.
[0031] By retrieving the grid data corresponding to the Mortan code converted from the longitude and latitude information, the sea-land situation of the beam center point and the feature points can be obtained as shown in Table 1.
[0032] Table 1: sea attribute Preferably, in step S3, whether to perform the current beam scheduling according to the sea-land attribute, further comprising: If the beam center point coordinate is land, more than half of the edge feature point longitude and latitude coordinate is land, and two or more of the inner circle feature point coordinate is land, the current beam scheduling is not performed; otherwise, the current beam is normally allocated. Specifically, in the embodiment, the beam center O is land, two of the four beam edge feature points are land, and three of the four inner circle feature points are land, so the current beam scheduling is not performed.
[0033] As can be seen from the above embodiment, for the beam with a large proportion of land in the beam coverage range, the beam allocation strategy of the present application can effectively avoid the detection of non-target regions such as land, islands and reefs, can efficiently utilize the limited processing resources, has low computational complexity, small storage requirement, and has high engineering application value.
[0034] Second embodiment Based on the same concept, the present embodiment provides a radar beam real-time allocation system based on geographic information system, comprising: an angle calculation module, which obtains the beam scanning angle according to the radar beam table, and obtains the beam edge feature point scanning angle according to the beam scanning angle, the azimuth beam width and the pitch beam width; The coordinate conversion module calculates the coordinates of the beam center point and the longitude and latitude coordinates of the edge feature points of the beam based on satellite orbit parameters, and obtains the coordinates of the inner circle feature points according to the coordinates of the beam center point and the longitude and latitude coordinates of the edge feature points. The judgment decision module judges the sea-land attribute of the beam center point coordinates, the longitude and latitude coordinates of the edge feature points and the coordinates of the inner circle feature points through a geographic information system (GIS) grid database, and judges whether to perform current beam scheduling according to the sea-land attribute.
[0035] Third embodiment In this embodiment, a computer device is provided, which includes a memory and one or more processors, and the memory stores computer code. When the computer code is executed by the one or more processors, the one or more processors execute the steps of the geographic information system-based radar beam real-time allocation method in the first embodiment.
[0036] In some embodiments of the present application, a computer readable storage medium is also provided, and the computer readable instructions are executed by one or more processors to make the one or more processors execute the steps of the geographic information system-based radar beam real-time allocation method as described in any one of the first embodiment.
[0037] It can be understood that, for the aforementioned geographic information system-based radar beam real-time allocation method, if each is implemented in the form of a software function module and sold or used as an independent product, it can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or the whole or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a number of instructions to make a computer device (which can be a personal computer server, or a network device, etc.) execute all or part of the steps of the embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various program code storage media.
[0038] The computer readable storage medium can include a tangible medium such as those previously described. Additionally, or alternatively, the computer readable storage medium can represent a transmission medium such as those described herein.
[0039] The above merely describes the preferred embodiments of the present application, and the protection scope of the present application is not limited to the above-mentioned embodiments. Any technical scheme falling within the concept of the present application shall fall within the protection scope of the present application. It should be noted that, for ordinary skilled persons in the art, some improvements and refinements without departing from the principles of the present application shall also be considered as falling within the protection scope of the present application.
Claims
1. A method for real-time allocation of radar beams based on geographic information system, characterized in that, The method comprises the following steps: S1: obtaining a beam scanning angle according to a radar wave position table, and obtaining a beam edge feature point scanning angle according to the beam scanning angle, an azimuth beam width and a pitch beam width; S2: calculating a beam center point coordinate and a longitude and latitude coordinate of a beam edge feature point based on satellite orbit parameters, and obtaining an inner circle feature point coordinate according to the beam center point coordinate and the longitude and latitude coordinate of the edge feature point; S3: judging a sea-land attribute of the beam center point coordinate, the longitude and latitude coordinate of the edge feature point and the inner circle feature point coordinate through a geographic information system (GIS) grid database, and judging whether to perform current beam scheduling according to the sea-land attribute.
2. The GIS-based real-time radar beam assignment method of claim 1, wherein, In step S1, the beam edge feature point scanning angle is obtained according to the beam scanning angle, the azimuth beam width and the pitch beam width, and the step further comprises: According to the beam scanning angle And the azimuth beam width The azimuth feature point scanning angle corresponding to the two feature points at the edge of the azimuth beam is 、 ; According to the beam scanning angle And the elevation beam width The elevation feature point scanning angle corresponding to the two feature points on the edge of the elevation beam is 、 .
3. The GIS-based real-time radar beam assignment method of claim 1, wherein, In step S2, the beam center point coordinate and the longitude and latitude coordinate of the beam edge feature point are calculated based on the satellite orbit parameters, and the step further comprises: S211: estimating an antenna local position coordinate of the beam center point in an antenna local coordinate system according to a distance corresponding to the beam center point and the beam scanning angle, as follows: wherein, is the distance corresponding to the beam center point, is the beam scanning angle at time t; S212: calculating a geocentric position coordinate of the beam center point in a geocentric fixed system according to the antenna local position coordinate and the satellite orbit parameters; S213: obtaining the beam center point coordinate according to the geocentric position coordinate, and calculating the longitude and latitude coordinate of the beam edge feature point according to a distance corresponding to the beam edge feature point and a scanning angle.
4. The GIS-based real-time radar beam assignment method of claim 3, wherein, In step S212, the geocentric position coordinate of the beam center point in the geocentric fixed system is calculated according to the antenna local position coordinate and the satellite orbit parameters, and the step further comprises: The geocentric position coordinate of the beam center point in the geocentric fixed system is calculated according to the satellite orbit parameters including a satellite platform orbit height, an earth radius, an ascending node right ascension and a perigee amplitude, and the antenna local position coordinate, as follows: wherein, is the current time, is the satellite platform orbital height, is the Earth radius, denotes the angle between the prime meridian and the equinoctial axis at the time, denotes the right ascension of the ascending node, denotes the argument of perigee, denotes the satellite orbit inclination, denotes the angle of the satellite movement from the perigee at the time, i.e. the true anomaly.
5. The GIS-based real-time radar beam assignment method of claim 4, wherein, In step S213, the beam center point coordinate is obtained according to the geocentric position coordinate, and the step further comprises: obtaining the beam center point coordinates according to the earth center position coordinates The calculation formula is as follows: wherein denotes the beam center point longitude, denotes the beam center point latitude, arcsin(.) is the inverse sine function, cos(.) is the inverse cosine function.
6. The GIS-based real-time radar beam assignment method of claim 1, wherein, In step S2, the inner circle feature point coordinate is obtained according to the beam center point coordinate and the longitude and latitude coordinate of the edge feature point, and the step further comprises: According to the beam center point coordinates and the beam edge feature points of the latitude and longitude coordinates to obtain the inner circle feature points coordinates as follows: Taking value range for longitude again Limitation: wherein, denotes the inner circle feature point coordinate longitude, denotes the inner circle feature point coordinate latitude, denotes the beam edge feature point coordinate longitude, denotes the beam edge feature point coordinate latitude; The inner circle feature point is obtained according to the beam center point coordinate and the longitude and latitude coordinate of the edge feature point The longitude and latitude coordinate of the inner circle feature point outside the coordinate.
7. The GIS-based real-time radar beam assignment method of claim 1, wherein, In step S3, the sea-land attribute of the beam center point coordinate, the longitude and latitude coordinate of the edge feature point and the inner circle feature point coordinate is judged through the geographic information system (GIS) grid database, and the step further comprises: S31: dividing the earth surface into grids, converting map data into grid data by using a geographic coordinate conversion algorithm, and generating a Morton code for the grid data; S32: converting geographic coordinates of the beam center point coordinate, the longitude and latitude coordinate of the edge feature point and the inner circle feature point coordinate into coordinate row and column numbers in the grid data, and converting the coordinate row and column numbers into the Morton code by using a coding rule of the Morton code; S33: searching the grid data corresponding to the Morton code to determine whether the position is a land attribute or a sea attribute.
8. The GIS-based real-time radar beam assignment method of claim 7, wherein, In step S3, whether to perform the current beam scheduling is judged according to the sea-land attribute, and the step further comprises: If the beam center point coordinate is land, the land proportion in the longitude and latitude coordinate of the edge feature point is more than half, and two or more feature points in the inner circle feature point coordinate are land, the current beam scheduling is not performed; otherwise, the current beam is normally distributed.
9. A real-time radar beam assignment system based on geographic information system, characterized in that, Comprise: An angle calculation module, which obtains a beam scanning angle according to a radar wave table, and obtains a beam edge feature point scanning angle according to the beam scanning angle, a beam width in an azimuth direction, and a beam width in an elevation direction; A coordinate conversion module, which calculates a beam center point coordinate and a longitude and latitude coordinate of a beam edge feature point based on satellite orbit parameters, and obtains an inner circle feature point coordinate according to the beam center point coordinate and the longitude and latitude coordinate of the edge feature point; A judgment and decision module, which judges a sea-land attribute of the beam center point coordinate, the longitude and latitude coordinate of the edge feature point, and the inner circle feature point coordinate through a geographic information system (GIS) grid database, and judges whether to perform current beam scheduling according to the sea-land attribute. 10.A computer readable storage medium, which stores computer code, when the computer code is executed, the steps of the geographic information system based radar beam real-time distribution method according to any one of claims 1-8 are executed.