Radar observation candidate body screening method based on Yarkovsky effect

Abstract

The invention relates to the technical field of radar observation candidate body screening, in particular to a radar observation candidate body screening method based on a Yarkovsky effect, and the method comprises the steps: determining a cataloged asteroid meeting a target screening condition as a first asteroid in all cataloged asteroids; wherein the target screening condition is that a Yarkovsky observation value exists, the semi-major axis drift rate is in a first threshold range, and the asteroid diameter is in a second threshold range; and according to radar parameter information of the target radar and ephemeris data of each first asteroid, determining the first asteroid with the visible segmental arc duration as a second asteroid, and determining the second asteroid with the visible segmental arc duration not less than a third threshold value as an observation candidate body of the target radar. According to the method, candidate bodies with high cost performance can be screened for asteroid Yarkovsky effect radar observation in the future, and technical support is provided for improving asteroid Yarkovsky effect measuring capacity.

Description

Background Art

[0002] The classical dynamic model of asteroids only considers gravitational effects and collisions. The number of near-Earth asteroids and cosmic ray exposure time deduced by this model are different from actual observations. After introducing the Yarkovsky effect caused by non-uniform thermal radiation into the dynamic model, the results are consistent with the observations, proving that the role of the Yarkovsky effect on the long-term orbital evolution of asteroids cannot be ignored. The accurate intensity of the Yarkovsky effect can enhance the long-term orbital deduction and collision warning capabilities.

[0003] At present, the main detection methods of the Yarkovsky effect on asteroids include optical detection, radar detection and deep space probe detection. Among them, optical observation is the most important observation method because of its low cost, long range and wide viewing angle. However, due to the influence of weather and sunlight, the optical observation arc is short and the determination period of the Yarkovsky effect is long. The Yarkovsky effect of asteroids measured by deep space probes has the highest accuracy, but the cost of this method is too high and the number of asteroids detected is limited. Radar observation of asteroids has the advantages of all-day and all-weather, high ranging accuracy and high imaging resolution. With the development of distributed ultra-long-range radar technology, the ability of radar to observe asteroids has been continuously enhanced, and the demand for screening candidates for radar observation of Yarkovsky effect asteroids is increasing.

[0004] Therefore, it is urgent to provide a technical solution to solve the above problems. Summary of the invention

[0005] In order to solve the above technical problems, the present invention provides a method for screening radar observation candidates based on the Yarkovsky effect.

[0006] In a first aspect, the present invention provides a method for screening radar observation candidates based on the Yarkovsky effect, and the technical solution of the method is as follows:

[0007] Among all catalogued asteroids, a catalogued asteroid that meets a target screening condition is determined as a first asteroid; wherein the target screening condition is: having a Yarkovsky observation value, a semi-major axis drift rate within a first threshold range, and an asteroid diameter within a second threshold range;

[0008] According to the radar parameter information of the target radar and the ephemeris data of each first asteroid, it is determined whether each first asteroid has a visible arc duration, the first asteroid with a visible arc duration is determined as the second asteroid, and the second asteroid with a visible arc duration not less than a third threshold is determined as the observation candidate of the target radar.

[0009] The beneficial effects of the radar observation candidate screening method based on the Yarkovsky effect of the present invention are as follows:

[0010] The method of the present invention can screen cost-effective candidates for future Yarkovsky effect radar observations of asteroids, and provide technical support for improving the ability to measure the Yarkovsky effect of asteroids.

[0011] On the basis of the above scheme, the radar observation candidate screening method based on the Yarkovsky effect of the present invention can also be improved as follows.

[0012] In an optional manner, the method further includes:

[0013] Based on the IMCCE SsODNet database, obtain the name, astronomical provisional name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid;

[0014] According to the name of any catalogued asteroid and the temporary astronomical name, the SPKID name of any catalogued asteroid is determined, and according to the SPKID name of any catalogued asteroid, the ephemeris data of any catalogued asteroid is determined until the ephemeris data of each catalogued asteroid is determined.

[0015] In an optional manner, the radar parameter information includes: a radar power threshold and a radar elevation angle threshold.

[0016] In an optional manner, the step of determining whether any first asteroid has a visible arc duration according to the radar parameter information of the target radar and the ephemeris data of any first asteroid includes:

[0017] According to the ephemeris data of any first asteroid, obtaining a first coordinate time series of the orbit of any first asteroid in the Earth-centered Earth-fixed coordinate system during a planned observation period;

[0018] According to the latitude and longitude information of the ground observation station where the target radar is located, the first coordinate time series is converted into a second coordinate time series in the northeast celestial coordinate system of the station center;

[0019] Determine the trajectory point distance modulus according to the second coordinate time series, and convert the second coordinate time series into a third coordinate time series in the radar AER coordinate system according to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast sky coordinate system;

[0020] It is determined whether any of the first asteroids has a visible arc duration according to the third coordinate time series, the radar power threshold and the radar elevation angle threshold.

[0021] In a second aspect, the present invention provides a radar observation candidate screening system based on the Yarkovsky effect, and the technical solution of the system is as follows:

[0022] It includes: a first screening module and a second screening module;

[0023] The first screening module is used to: determine, among all catalogued asteroids, a catalogued asteroid that meets a target screening condition as a first asteroid; wherein the target screening condition is: having a Yarkovsky observation value, a semi-major axis drift rate within a first threshold range, and an asteroid diameter within a second threshold range;

[0024] The second screening module is used to determine whether each first asteroid has a visible arc duration based on the radar parameter information of the target radar and the ephemeris data of each first asteroid, determine the first asteroid with a visible arc duration as the second asteroid, and determine the second asteroid with a visible arc duration not less than a third threshold as the observation candidate of the target radar.

[0025] The beneficial effects of the radar observation candidate screening system based on the Yarkovsky effect of the present invention are as follows:

[0026] The system of the present invention can screen cost-effective candidates for future Yarkovsky effect radar observations of asteroids, and provide technical support for improving the ability to measure the Yarkovsky effect of asteroids.

[0027] On the basis of the above scheme, the radar observation candidate screening system based on Yarkovsky effect of the present invention can also be improved as follows.

[0028] In an optional manner, the method further includes: an acquisition module; the acquisition module is used to:

[0029] Based on the IMCCE SsODNet database, obtain the name, astronomical provisional name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid;

[0030] According to the name of any catalogued asteroid and the temporary astronomical name, the SPKID name of any catalogued asteroid is determined, and according to the SPKID name of any catalogued asteroid, the ephemeris data of any catalogued asteroid is determined until the ephemeris data of each catalogued asteroid is determined.

[0031] In an optional manner, the radar parameter information includes: a radar power threshold and a radar elevation angle threshold.

[0032] In an optional manner, the second screening module is specifically used to:

[0033] According to the ephemeris data of any first asteroid, obtaining a first coordinate time series of the orbit of any first asteroid in the Earth-centered Earth-fixed coordinate system during a planned observation period;

[0034] According to the latitude and longitude information of the ground observation station where the target radar is located, the first coordinate time series is converted into a second coordinate time series in the northeast celestial coordinate system of the station center;

[0035] Determine the trajectory point distance modulus according to the second coordinate time series, and convert the second coordinate time series into a third coordinate time series in the radar AER coordinate system according to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast sky coordinate system;

[0036] It is determined whether any of the first asteroids has a visible arc duration according to the third coordinate time series, the radar power threshold and the radar elevation angle threshold.

[0037] In a third aspect, a technical solution of an electronic device of the present invention is as follows:

[0038] The invention comprises a memory, a processor and a program stored in the memory and running on the processor. When the processor executes the program, the steps of the radar observation candidate screening method based on the Yarkovsky effect of the present invention are implemented.

[0039] In a fourth aspect, the present invention provides a computer-readable storage medium having the following technical solution:

[0040] The computer-readable storage medium stores instructions. When the computer-readable storage medium reads the instructions, the computer-readable storage medium executes the steps of the method for screening radar observation candidate objects based on the Yarkovsky effect of the present invention.

[0041] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented according to the contents of the specification. In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] The accompanying drawings are only used to illustrate the embodiments and are not to be considered as limiting the present invention. In addition, the same reference symbols are used to represent the same components throughout the accompanying drawings. In the accompanying drawings:

[0043] Figure 1 It is a flow chart of an embodiment of a method for screening radar observation candidate objects based on the Yarkovsky effect of the present invention;

[0044] Figure 2Schematic diagram of the orbital distribution of asteroids for which Yarkovsky observations have been published;

[0045] Figure 3 Schematic diagram of the semi-major axis drift rate and diameter distribution of asteroids with published Yarkovsky observation values;

[0046] Figure 4 Schematic diagram of radar power and pitch angle screening for asteroid candidates;

[0047] Figure 5 It is a structural schematic diagram of an embodiment of a radar observation candidate screening system based on the Yarkovsky effect of the present invention;

[0048] Figure 6 The figure is a schematic structural diagram of an embodiment of an electronic device of the present invention. DETAILED DESCRIPTION

[0049] Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0050] Figure 1 A flow chart of an embodiment of a method for screening radar observation candidates based on the Yarkovsky effect provided by the present invention is shown. The method for screening radar observation candidates based on the Yarkovsky effect can be executed by an electronic device such as a terminal device or a server. The terminal device can be any fixed or mobile terminal such as a user equipment (UE), a mobile device, a user terminal, a terminal, a cellular phone, a cordless phone, a personal digital assistant (PDA), a handheld device, a computing device, a vehicle-mounted device, a wearable device, etc. The server can be a single server or a server cluster composed of multiple servers. Any electronic device can implement the method for screening radar observation candidates based on the Yarkovsky effect by calling computer-readable instructions stored in a memory by a processor. As Figure 1 As shown, the following steps are included:

[0051] S1. Among all catalogued asteroids, the catalogued asteroid that meets the target screening criteria is determined as the first asteroid.

[0052] Each cataloged asteroid includes: name, provisional astronomical name, asteroid diameter, semi-major axis drift rate and ephemeris data. ① The names of asteroids usually come from various cultural and historical backgrounds, which help to better understand and remember these celestial bodies. Asteroids may be named after mythological figures, place names, personal names, etc., or they may be named according to the characteristics of the asteroid or the process of discovery. ② Provisional astronomical name (Provisional designation inastronomy) is the name given to a celestial body immediately after it is discovered. When reliable orbital data is calculated, the provisional name will be replaced by a formal name. However, due to the large number of asteroids discovered, most of them will not have their orbits calculated within a short period of time (several years to decades) after their discovery, so there will be a long period of time when provisional names are used instead of formal names. The first element in the provisional name of an asteroid is the year of discovery, followed by two letters, and then subscript numbers added when necessary. The first letter represents the period of discovery, with one letter for each half month: "A" represents the first half of January, "D" represents the second half of February, "J" represents the first half of May ("I" is not used), and so on until "Y" represents the second half of December. Every month, whether it is a long month or a short month, has the first half of the month from the 1st to the 15th, and the second half from the 16th. The second letter and the numbers added after it represent the order of discovery within that half month, and the letter "I" is also omitted. For example, the 8th asteroid discovered in the second half of March 1950 would be recorded as 1950FH. However, with modern technology, more than 25 asteroids will be discovered in a half month, so subscript numbers are added to indicate how many times the letters have been cycled. In this way, the 28th asteroid discovered in the second half of March 1950 is 1950FC1. However, due to the limitation of computer ASCII code, the numbers that should be subscripted are sometimes directly arranged at the same height, so it may also be directly marked as 1950FC1. 1926GA1 is the first asteroid to use the subscript number marking system. An interesting phenomenon derived from the subscript system is that although the second letter "advances" faster than the following numbers, it is placed before. This is different from other naming systems in the world. ③ The diameter of an asteroid refers to its physical size, usually measured and expressed in kilometers. The diameter of an asteroid varies with its size and type, ranging from a few meters or tens of meters to hundreds of kilometers or even thousands of kilometers. The diameter of an asteroid is one of its basic physical properties and is of great significance for understanding the structure, composition and evolution of asteroids. By observing and measuring the diameter of an asteroid, its physical parameters such as mass and density can be inferred, and its internal structure and composition can be studied. ④ The semi-major axis drift rate refers to the rate at which the semi-major axis of a celestial body's orbit changes over time. In astronomy, the semi-major axis is one of the important parameters that describe the shape and size of a celestial body's orbit.For celestial bodies such as asteroids, their semi-major axis drift rate may be affected by many factors, such as gravitational perturbations, non-gravitational perturbations (such as solar radiation pressure, planetary atmospheric drag, etc.) and the physical properties of the celestial body itself (such as shape, rotation, etc.). ⑤ Ephemeris data is important information that describes the position and motion state of an asteroid in the solar system, including: the asteroid number, name, reduced Julian day (epoch), semi-major axis (SemiAxis), eccentricity (Eccentricity), orbital inclination (Inclination), perihelion angular distance (W), ascending node right ascension (Omega), mean anomaly angle (Mean) and other parameters.

[0053] The target screening conditions are: having Yarkovsky observation values, the semi-major axis drift rate is within the first threshold range, and the asteroid diameter is within the second threshold range. The first threshold range is set by default to: 5×10 -5 -5×10 -2 AU / Myr (Astronomical Unit per million years, astronomical unit is equal to the average distance between the earth and the sun, defined as 149597870.7km), the second threshold range is set by default to: 1m-30km, which can also be set according to actual conditions, and there is no limit here. Yarkovsky observation values ​​refer to the Yarkovsky observation values ​​published by asteroids that have been preliminarily observed by optical detection, other radar detection, and deep space probe detection (such as Figure 2 shown).

[0054] It should be noted that Figure 3 The diagram shows the distribution of the semi-major axis drift rate and diameter of asteroids with published Yarkovsky observations, i.e. Figure 3 The corresponding asteroid is the first asteroid in this embodiment.

[0055] S2. Determine whether each first asteroid has a visible arc duration according to the radar parameter information of the target radar and the ephemeris data of each first asteroid, determine the first asteroid with a visible arc duration as the second asteroid, and determine the second asteroid with a visible arc duration not less than a third threshold as an observation candidate of the target radar.

[0056] Among them, the first asteroid is: a cataloged asteroid that meets the target screening conditions, that is, an asteroid that has been screened by the observation effect. The second asteroid is: the first asteroid with a visible arc segment duration, and the observation candidate is: the second asteroid with a visible arc segment duration not less than the third threshold. The third threshold is set by default to: 30s, and can also be adjusted according to actual conditions. The target radar is: a radar for observing asteroids affected by the Yarkovsky effect in this embodiment. The target radar in this embodiment defaults to the next generation Arecibo radar (NGAT-130 solution), which can also be adjusted according to actual conditions. Radar parameter information includes but is not limited to: radar power threshold and radar elevation angle threshold.

[0057] In an optional manner, the method further includes:

[0058] Based on the IMCCE SsODNet database, the name, astronomical provisional name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid are obtained.

[0059] The full name of the IMCCE SsODNet database is: IMCCE (Institut de mécanique céleste et de calcul deséphémérides) SsODNet (Solar system Open Database Network) database. The IMCCE SsODNet database stores information such as the name, temporary astronomical name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid.

[0060] According to the name of any catalogued asteroid and the temporary astronomical name, the SPKID name of any catalogued asteroid is determined, and according to the SPKID name of any catalogued asteroid, the ephemeris data of any catalogued asteroid is determined until the ephemeris data of each catalogued asteroid is determined.

[0061] Among them, in the NASA JPL (Jet Propulsion Laboratory) SDB (Small-Body DataBase) database, the SPKID name of any cataloged asteroid is obtained according to the name and temporary astronomical name of the cataloged asteroid; in the Horizon System, the ephemeris data of the cataloged asteroid is obtained according to the SPKID name of the cataloged asteroid. Repeat the above process until the ephemeris data of each cataloged asteroid is obtained.

[0062] It should be noted that the Horizon System refers to the online solar system data and ephemeris calculation service system opened by NASA JPL (Jet Propulsion Laboratory of National Aeronautics and Space Administration), through which the ephemeris data of asteroids can be obtained.

[0063] In an optional manner, the step of determining whether any first asteroid has a visible arc duration according to the radar parameter information of the target radar and the ephemeris data of any first asteroid includes:

[0064] According to the ephemeris data of any first asteroid, a first coordinate time series of the orbit of any first asteroid in the Earth-centered Earth-fixed coordinate system within a planned observation period is obtained.

[0065] Among them, the first coordinate time series is a coordinate time series in the geocentric earth-fixed coordinate system, including the coordinates of the orbit of the first asteroid in the geocentric earth-fixed coordinate system corresponding to multiple time points in the planned observation period. The planned observation period is a pre-planned observation period, such as the time period corresponding to 0:00 on January 1, 2025 to 0:00 on January 31, 2025. The ephemeris data includes the orbital parameters of the asteroid (such as semi-major axis, eccentricity, inclination, etc.) and time-related position information. Specifically: ① Determine the specific parameters of the geocentric earth-fixed coordinate system used, such as WGS84. ② Make the time system used consistent with the time system of the ephemeris data, usually using Julian day or UTC time. ③ Using the orbital parameters of any of the first asteroids, calculate the position of any of the first asteroids at each time point in the planned observation period through a mathematical model (such as Kepler's equation), and these positions are usually expressed in coordinates under the celestial equatorial coordinate system (such as the equatorial inertial coordinate system). ④ Use mathematical transformation (rotation and translation operations) to convert the coordinates in the celestial equatorial coordinate system into the coordinates in the Earth-centered Earth-fixed coordinate system. ⑤ For each time point in the planned observation period, repeat the steps in ④ until the first coordinate time series of the orbit of any first asteroid in the planned observation period in the Earth-centered Earth-fixed coordinate system is obtained.

[0066] According to the latitude and longitude information of the ground observation station where the target radar is located, the first coordinate time series is converted into a second coordinate time series in the station-centered northeast celestial coordinate system.

[0067] Among them, the second coordinate time series is a coordinate time series in the station-centered northeast sky coordinate system, which contains the coordinates of the first asteroid's orbit in the station-centered northeast sky coordinate system corresponding to multiple time points in the planned observation period. The longitude and latitude height information includes: latitude, longitude and altitude. The longitude and latitude height information of the ground observation station where the target radar is located in this embodiment is (40, 90, 0). Specifically: ① According to the longitude and latitude height information, calculate the rotation matrix from the geocentric earth-fixed coordinate system to the station-centered northeast sky coordinate system. Involving three rotation steps: first rotate around the Z axis (longitude), then rotate around the new Y axis (latitude), and fine-tune in combination with the curvature and elevation of the earth's surface to obtain a rotation matrix. ② Use the rotation matrix to convert the coordinates of the asteroid in the geocentric earth-fixed coordinate system (each coordinate in the first coordinate time series) into coordinates in the station-centered northeast sky coordinate system. ③ Since the U axis of the station-centered northeast celestial coordinate system is perpendicular to the ground, it is also necessary to adjust the coordinates in the converted station-centered northeast celestial coordinate system according to the altitude of the ground observation station, and finally obtain the second coordinate time series of the orbit of any first asteroid in the station-centered northeast celestial coordinate system during the planned observation period.

[0068] According to the second coordinate time series, the trajectory point distance modulus is determined, and according to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast celestial coordinate system, the second coordinate time series is converted into a third coordinate time series in the radar AER coordinate system.

[0069] Among them, the third coordinate time series is the coordinate time series in the radar AER coordinate system, which contains the coordinates of the first asteroid's orbit in the radar AER coordinate system corresponding to multiple time points in the planned observation period. Using the Euclidean distance formula and based on the second coordinate time series, the trajectory point distance modulus is calculated. According to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast sky coordinate system, the process of converting the second coordinate time series into the third coordinate time series in the radar AER coordinate system is as follows: ① Record the trajectory point P corresponding to the second coordinate time series in the station center northeast sky coordinate system O enu x enu y enu z enu The Cartesian coordinates (trajectory point coordinates) under enu =[x enu ,y enu ,z enu ], then the distance modulus of the trajectory point is ② The angle between the coordinates of the trajectory point and the coordinate plane of the northeast celestial coordinate system of the station center is the trajectory point vector, that is, P enu =[x enu ,y enu ,z enu ] and the coordinate plane of the northeast celestial coordinate system, denoted by Penu =[x enu ,y enu ,z enu ] and x enu O enu z enu The angle between the coordinate planes is θ enu , P enu =[x enu ,y enu ,z enu ] and x enu O enu y enu The angle between the coordinate planes is Then the radar AER coordinates of the trajectory point are That is: the third coordinate time series of the first asteroid in the radar AER coordinate system.

[0070] It is determined whether any of the first asteroids has a visible arc duration according to the third coordinate time series, the radar power threshold and the radar elevation angle threshold.

[0071] in, Figure 4 A schematic diagram of radar power and elevation angle screening of asteroid candidates is shown. The radar power threshold and the radar elevation angle threshold are determined according to the parameters of the target radar, and different target radars have different radar power thresholds and radar elevation angle thresholds. In this embodiment, the radar elevation angle threshold defaults to 10°, and the radar power defaults to the radar power corresponding to the next generation Arecibo radar (NGAT-130 solution). According to the third coordinate time series, the radar power threshold and the radar elevation angle threshold, the specific process of determining whether any of the first asteroids has a visible arc duration is: ① Calculate the observation distance threshold (radar power threshold) according to the following radar power equation P t G t G r is the radar power aperture product of the target radar, S imin is the minimum detectable signal power, λ is the radar signal wavelength of the target radar, and σ is the asteroid radar cross-section RCS. ② In the third coordinate time series of any first asteroid, R≤R max and The arc segment is the visible arc segment, and the starting and ending times are t begin With t end , then the length of the visible arc of the first asteroid is t visible =t end -t begin ; If there is no third coordinate time series of any first asteroid that satisfies R≤R max and arc segment, it is determined that the first asteroid does not have a visible arc segment duration.

[0072] The technical solution of this embodiment can screen cost-effective candidates for future asteroid Yarkovsky effect radar observations and provide technical support for improving the asteroid Yarkovsky effect measurement capability.

[0073] Figure 5 FIG. 2 is a schematic diagram showing a structure of an embodiment of a radar observation candidate screening system 200 based on the Yarkovsky effect provided by the present invention. Figure 5 As shown, the system 200 includes: a first screening module 210 and a second screening module 220;

[0074] The first screening module 210 is used to: determine, among all catalogued asteroids, a catalogued asteroid that meets a target screening condition as a first asteroid; wherein the target screening condition is: having a Yarkovsky observation value, a semi-major axis drift rate within a first threshold range, and an asteroid diameter within a second threshold range;

[0075] The second screening module 220 is used to determine whether each first asteroid has a visible arc duration according to the radar parameter information of the target radar and the ephemeris data of each first asteroid, determine the first asteroid with a visible arc duration as the second asteroid, and determine the second asteroid with a visible arc duration not less than a third threshold as the observation candidate of the target radar.

[0076] In an optional manner, the method further includes: an acquisition module; the acquisition module is used to:

[0077] Based on the IMCCE SsODNet database, obtain the name, astronomical provisional name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid;

[0078] According to the name of any catalogued asteroid and the temporary astronomical name, the SPKID name of any catalogued asteroid is determined, and according to the SPKID name of any catalogued asteroid, the ephemeris data of any catalogued asteroid is determined until the ephemeris data of each catalogued asteroid is determined.

[0079] In an optional manner, the radar parameter information includes: a radar power threshold and a radar elevation angle threshold.

[0080] In an optional manner, the second screening module 220 is specifically used for:

[0081] According to the ephemeris data of any first asteroid, obtaining a first coordinate time series of the orbit of any first asteroid in the Earth-centered Earth-fixed coordinate system during a planned observation period;

[0082] According to the latitude and longitude information of the ground observation station where the target radar is located, the first coordinate time series is converted into a second coordinate time series in the northeast celestial coordinate system of the station center;

[0083] Determine the trajectory point distance modulus according to the second coordinate time series, and convert the second coordinate time series into a third coordinate time series in the radar AER coordinate system according to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast sky coordinate system;

[0084] It is determined whether any of the first asteroids has a visible arc duration according to the third coordinate time series, the radar power threshold and the radar elevation angle threshold.

[0085] It should be noted that the beneficial effects of the radar observation candidate screening system based on the Yarkovsky effect provided in the above embodiment are the same as the beneficial effects of the radar observation candidate screening method based on the Yarkovsky effect, which will not be described in detail here. In addition, when the system provided in the above embodiment realizes its functions, it only takes the division of the above functional modules as an example. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the system can be divided into different functional modules according to actual conditions to complete all or part of the functions described above. In addition, the system and method embodiments provided in the above embodiment belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be described in detail here.

[0086] Among them, the radar observation candidate screening system based on Yarkovsky effect of the present invention can be a computer program (including program code) running in a computer device. For example, the radar observation candidate screening system based on Yarkovsky effect of the present invention is an application software that can be used to execute the corresponding steps in the radar observation candidate screening method based on Yarkovsky effect of the present invention.

[0087] In some embodiments, the radar observation candidate screening system based on the Yarkovsky effect of the present invention can be implemented in a combination of software and hardware. As an example, the radar observation candidate screening system based on the Yarkovsky effect of the present invention can be a processor in the form of a hardware decoding processor, which is programmed to execute the radar observation candidate screening method based on the Yarkovsky effect of the present invention. For example, the processor in the form of a hardware decoding processor can adopt one or more application specific integrated circuits (ASICs), DSPs, programmable logic devices (PLDs), complex programmable logic devices (CPLDs), field programmable gate arrays (FPGAs) or other electronic components.

[0088] The modules involved in the embodiments of the present invention may be implemented in software or hardware. The name of a module does not limit the module itself in some cases.

[0089] An electronic device according to an embodiment of the present invention includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, any of the above-mentioned radar observation candidate screening methods based on the Yarkovsky effect is implemented. That is, an electronic device according to an embodiment of the present invention may include but is not limited to: a processor and a memory; the memory is used to store the computer program; the processor is used to execute the radar observation candidate screening method based on the Yarkovsky effect shown in any embodiment of the present invention by calling the computer program.

[0090] In an alternative embodiment, an electronic device is provided, such as Figure 6 As shown, Figure 6 The electronic device 4000 shown includes: a processor 4001 and a memory 4003. The processor 4001 and the memory 4003 are connected, such as through a bus 4002. Optionally, the electronic device 4000 may also include a transceiver 4004, which may be used for data interaction between the electronic device and other electronic devices, such as data transmission and / or data reception. It should be noted that in actual applications, the transceiver 4004 is not limited to one, and the structure of the electronic device 4000 does not constitute a limitation on the embodiments of the present invention.

[0091] Processor 4001 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array) or other programmable logic devices, transistor logic devices, hardware components or any combination thereof. It may implement or execute various exemplary logic blocks, modules and circuits described in conjunction with the disclosure of the present invention. Processor 4001 may also be a combination that implements computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, etc.

[0092] The bus 4002 may include a path to transmit information between the above components. The bus 4002 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus. The bus 4002 may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 6 In the figure, only one thick line is used to represent the bus 4002, but this does not mean that there is only one bus or one type of bus.

[0093] The memory 4003 can be a ROM (Read Only Memory) or other types of static storage devices that can store static information and instructions, a RAM (Random Access Memory) or other types of dynamic storage devices that can store information and instructions, or an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read Only Memory) or other optical disk storage, optical disk storage (including compressed optical disk, laser disk, optical disk, digital versatile disk, Blu-ray disk, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited to this.

[0094] The memory 4003 is used to store the application code (computer program) for executing the solution of the present invention, and the execution is controlled by the processor 4001. The processor 4001 is used to execute the application code stored in the memory 4003 to implement the content shown in the above method embodiment.

[0095] Among them, the electronic device can also be a terminal device, and the terminal device can be any terminal device that can install applications and access web pages through applications, including at least one of a smart phone, a tablet computer, a laptop computer, a desktop computer, a smart speaker, a smart watch, a smart TV, and a smart car device.

[0096] It should be noted that Figure 6 The electronic device shown is only an example and should not bring any limitation to the functions and scope of use of the embodiments of the present invention.

[0097] A computer-readable storage medium according to an embodiment of the present invention stores a computer program, and when the computer program is executed by a processor, any of the above-mentioned radar observation candidate screening methods based on the Yarkovsky effect is implemented.

[0098] Optionally, the computer-readable storage medium may be a read-only memory (ROM), a random access memory (RAM), a compact disc (CD-ROM), a magnetic tape, a floppy disk, an optical data storage device, and the like.

[0099] In an exemplary embodiment, a computer program product or a computer program is also provided, the computer program product or the computer program including computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the electronic device performs the above-mentioned radar observation candidate screening method based on the Yarkovsky effect.

[0100] Computer program code for performing the operations of the present invention may be written in one or more programming languages, or a combination thereof, including object-oriented programming languages, such as Java, Smalltalk, C++, and conventional procedural programming languages, such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., through the Internet using an Internet service provider).

[0101] It should be understood that the flow charts and block diagrams in the accompanying drawings illustrate the possible implementation architecture, functions and operations of the method and computer program product according to various embodiments of the present invention. In this regard, each box in the flow chart or block diagram can represent a module, a program segment, or a part of a code, and the module, the program segment, or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a sequence different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flow chart, and the combination of the boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0102] The computer-readable storage medium provided by the embodiment of the present invention may be, but is not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or component, or any combination thereof. More specific examples of computer-readable storage media may include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, device, or component.

[0103] The computer-readable storage medium carries one or more programs. When the one or more programs are executed by the electronic device, the electronic device executes the method shown in the above embodiment.

[0104] The above description is only a preferred embodiment of the present invention and an explanation of the technical principles used. Those skilled in the art should understand that the disclosure scope involved in the present invention is not limited to the technical solution formed by a specific combination of the above technical features, but should also cover other technical solutions formed by any combination of the above technical features or their equivalent features without departing from the above disclosed concept. For example, the above features are replaced with the technical features with similar functions disclosed in the present invention (but not limited to) by each other.

[0105] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application are used to distinguish similar objects and represent the definition of a specific order or sequence. The order of use of similar objects can be interchanged where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than the order shown or described.

[0106] Those skilled in the art know that the present invention can be implemented as a system, method or computer program product. Therefore, the present invention can be specifically implemented in the following forms, namely: it can be complete hardware, it can be complete software (including firmware, resident software, microcode, etc.), or it can be a combination of hardware and software, which is generally referred to as "circuit", "module" or "system" herein. In addition, in some embodiments, the present invention can also be implemented in the form of a computer program product in one or more computer-readable media, and the computer-readable medium contains computer-readable program code.

[0107] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A method for screening radar observation candidates based on the Yarkovsky effect, characterized in that: include: Among all catalogued asteroids, a catalogued asteroid that meets a target screening condition is determined as a first asteroid; wherein the target screening condition is: having a Yarkovsky observation value, a semi-major axis drift rate within a first threshold range, and an asteroid diameter within a second threshold range; According to the radar parameter information of the target radar and the ephemeris data of each first asteroid, it is determined whether each first asteroid has a visible arc duration, the first asteroid with a visible arc duration is determined as the second asteroid, and the second asteroid with a visible arc duration not less than a third threshold is determined as the observation candidate of the target radar.

2. The method for screening radar observation candidates based on the Yarkovsky effect according to claim 1, characterized in that: Also includes: Based on the IMCCE SsODNet database, obtain the name, astronomical provisional name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid; According to the name of any catalogued asteroid and the temporary astronomical name, the SPKID name of any catalogued asteroid is determined, and according to the SPKID name of any catalogued asteroid, the ephemeris data of any catalogued asteroid is determined until the ephemeris data of each catalogued asteroid is determined.

3. The method for screening radar observation candidates based on the Yarkovsky effect according to claim 1, characterized in that: The radar parameter information includes: a radar power threshold and a radar elevation angle threshold.

4. The method for screening radar observation candidates based on the Yarkovsky effect according to claim 3, characterized in that: The step of determining whether any of the first asteroids has a visible arc duration according to the radar parameter information of the target radar and the ephemeris data of any of the first asteroids comprises: According to the ephemeris data of any first asteroid, obtaining a first coordinate time series of the orbit of any first asteroid in the Earth-centered Earth-fixed coordinate system during a planned observation period; According to the latitude and longitude information of the ground observation station where the target radar is located, the first coordinate time series is converted into a second coordinate time series in the northeast celestial coordinate system of the station center; Determine the trajectory point distance modulus according to the second coordinate time series, and convert the second coordinate time series into a third coordinate time series in the radar AER coordinate system according to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast sky coordinate system; It is determined whether any of the first asteroids has a visible arc duration according to the third coordinate time series, the radar power threshold and the radar elevation angle threshold.

5. A radar observation candidate screening system based on the Yarkovsky effect, characterized in that: include: a first screening module and a second screening module; The first screening module is used to: determine, among all catalogued asteroids, a catalogued asteroid that meets a target screening condition as a first asteroid; wherein the target screening condition is: having a Yarkovsky observation value, a semi-major axis drift rate within a first threshold range, and an asteroid diameter within a second threshold range; The second screening module is used to determine whether each first asteroid has a visible arc duration based on the radar parameter information of the target radar and the ephemeris data of each first asteroid, determine the first asteroid with a visible arc duration as the second asteroid, and determine the second asteroid with a visible arc duration not less than a third threshold as the observation candidate of the target radar.

6. The radar observation candidate screening system based on Yarkovsky effect according to claim 5, characterized in that: Also includes: Acquisition module; the acquisition module is used to: Based on the IMCCE SsODNet database, obtain the name, astronomical provisional name, semi-major axis drift rate and asteroid diameter of each catalogued asteroid; According to the name of any catalogued asteroid and the temporary astronomical name, the SPKID name of any catalogued asteroid is determined, and according to the SPKID name of any catalogued asteroid, the ephemeris data of any catalogued asteroid is determined until the ephemeris data of each catalogued asteroid is determined.

7. The radar observation candidate screening system based on Yarkovsky effect according to claim 6, characterized in that: The radar parameter information includes: a radar power threshold and a radar elevation angle threshold.

8. The radar observation candidate screening system based on Yarkovsky effect according to claim 7, characterized in that: The second screening module is specifically used for: According to the ephemeris data of any first asteroid, obtaining a first coordinate time series of the orbit of any first asteroid in the Earth-centered Earth-fixed coordinate system during a planned observation period; According to the latitude and longitude information of the ground observation station where the target radar is located, the first coordinate time series is converted into a second coordinate time series in the northeast celestial coordinate system of the station center; Determine the trajectory point distance modulus according to the second coordinate time series, and convert the second coordinate time series into a third coordinate time series in the radar AER coordinate system according to the trajectory point distance modulus and the angle between the trajectory point coordinates and the coordinate plane of the station center northeast sky coordinate system; It is determined whether any of the first asteroids has a visible arc duration according to the third coordinate time series, the radar power threshold and the radar elevation angle threshold.

9. An electronic device, characterized in that: The electronic device includes a processor, the processor is coupled to a memory, at least one computer program is stored in the memory, and the at least one computer program is loaded and executed by the processor so that the electronic device implements the radar observation candidate screening method based on the Yarkovsky effect as described in any one of claims 1 to 4.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores at least one computer program, and the at least one computer program is loaded and executed by the processor so that the computer-readable storage medium implements the radar observation candidate screening method based on the Yarkovsky effect as described in any one of claims 1 to 4.

Images