Method and system for detecting tumbling characteristics of a space object

Abstract

A method for determining the tumbling characteristics of a resident space object includes obtaining the geometry of a target object, generating a distribution of radio wave scatterers within the geometry of the target, calculating multiple radar cross sections at different angles to the distribution, determining an autocorrelation of each radar cross section versus each angle, determining a decorrelation angle using the autocorrelation, extracting decorrelation times from the radar data of the object, and determining the rotation rate of the object as a tumbling characteristic.

Description

[Technical Field]

[0001] The present application relates to a method and system for detecting the tumbling characteristics of space objects, particularly satellites. [Background technology]

[0002] In the field of space situational awareness, it is desirable to be able to characterize the tumbling characteristics of a resident space object (RSO) remotely, without the cooperation of the RSO itself. This is relevant for many applications, including, but not limited to, alerting satellite operators if a satellite begins a tumbling motion, characterizing stabilized objects that may indicate they are actively controlled satellites, and understanding the characteristics of space debris and other non-cooperative objects.

[0003] Currently available approaches include optical "light curve" measurements, which show the change in brightness of an object as a function of time. Other approaches include using radiometric data, which show the scattered or emitted signal strength of an object as a function of time. However, both of these approaches rely on some knowledge or data about the object being available. Furthermore, such approaches require the object to be observed and measured over an extended period of time, typically at least one 360-degree rotation, so that periodic changes in brightness or signal strength can be identified.

[0004] It is therefore desirable to provide an improved means of characterizing the tumbling characteristics of space objects.

[0005] The embodiments described below are not limited to implementations that address any or all of the drawbacks of the known approaches described above. Summary of the Invention

[0006] This Summary is provided to introduce a selection of concepts in a simplified form that are further described below in the Detailed Description. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used as an aid in determining the scope of the claimed subject matter.

[0007] In a first aspect, the present disclosure provides a method for determining at least one tumbling characteristic of an object, the method including: obtaining a decorrelation angle of an expected radar cross section (RCS) of the object; obtaining radar data of the object; determining a decorrelation time of the RCS of the object from the radar data of the object; and determining the at least one tumbling characteristic of the object using the obtained decorrelation angle and the determined decorrelation time.

[0008] In a second aspect, the present disclosure provides a system configured to perform the method of the first aspect.

[0009] In a third aspect, the present disclosure provides a system for determining tumbling characteristics of objects, the system comprising: a radar that acquires radar data for each object; and at least one processor configured to execute code that causes the processor to acquire a decorrelation angle of a radar cross section (RCS) of each object; determine a decorrelation time of the RCS of each object from the object's radar data; and determine at least one tumbling characteristic of each object using the obtained decorrelation angle and the determined decorrelation time for each object.

[0010] The methods described herein may be performed by software in machine-readable form on a tangible storage medium, and the computer program may be embodied on a computer-readable medium, for example, in the form of a computer program comprising computer program code means adapted to perform all the steps of any of the methods described herein when the program is run on a computer. Examples of tangible (or non-transitory) storage media include disks, thumb drives, memory cards, etc., but do not include propagated signals. The software may be suitable for execution on a parallel processor or a serial processor, such that the steps of the method may be performed in any suitable order, or simultaneously.

[0011] This application recognizes that firmware and software can have value and be separately tradable commodities. It is intended to encompass software that runs on or controls "dumb" or standard hardware to perform a desired function. It is also intended to encompass software that "describes" or defines the configuration of hardware to perform a desired function, such as HDL (Hardware Description Language) software used to design silicon chips or configure universal programmable chips.

[0012] The preferred features may be combined as appropriate and combined with any aspect of the present invention, as will be apparent to one skilled in the art.

[0013] Embodiments of the present invention will now be described, by way of example, with reference to the following drawings, in which: [Brief explanation of the drawings]

[0014] [Figure 1] 1 is an illustration of a radar system used to track a resident space object according to a first embodiment; FIG. [Figure 2] 1 is a flow diagram of an embodiment of a method for determining the rotational rate of a resident space object according to a first embodiment. [Figure 3] 4 is a flowchart of a method for determining decorrelation time in the first embodiment. [Figure 4] FIG. 1 is an explanatory diagram showing an example of a random scatterer model of an artificial satellite. [Figure 5] 1 shows a graph of autocorrelation versus angle of a resident space object for an ensemble of models. [Figure 6] 1 shows a graph of target and noise autocorrelation data versus window length. [Figure 7] 1 shows a graph of target and noise uncorrelated data versus window length. [Figure 8] Graphs related to targets are shown. [Figure 9] Graphs related to targets are shown. [Figure 10] Graphs relating to different targets are shown. [Figure 11] Graphs relating to different targets are shown. [Figure 12] 10 shows a graph relating to yet another target. [Figure 13] 10 shows a graph relating to yet another target. [Figure 14] 4 shows a flow chart of a method for determining a rotation vector of a resident space object according to a second embodiment. [Figure 15] An explanatory diagram of the method of FIG. 14 is shown. DETAILED DESCRIPTION OF THE INVENTION

[0015] Common reference numbers are used throughout the figures to denote similar features. Embodiments of the present invention are described below by way of example only. These examples represent the best ways of carrying out the invention currently known to the applicant, although they are not the only ways that this can be achieved. The specification describes the functions of the examples and the sequence of steps for constructing and operating the examples. However, the same or equivalent functions and arrangements may be achieved by different examples.

[0016] FIG. 1 illustrates a system 1 used to gather information about resident space objects (RSOs). The system includes a radar array 10 residing at a fixed location on Earth, transmitting signals toward RSOs, such as RSO 16, and receiving return signals reflected from the RSOs, and a signal analysis system 18 that determines data about the RSOs from the return signals. The RSOs may be, for example, debris, such as fragmented debris, or operational satellites. In FIG. 1, RSO 16 is an operational satellite. As the RSOs move across the sky, radar array 10 gathers data about them.

[0017] In addition to characterizing the tumbling characteristics of RSOs, the data collected by radar array 10 may be used for other purposes. In an example, the data collected by radar array 10 is also used by a system that tracks RSOs and determines their orbital paths. However, this is not required. In some examples, system 1 may be directed to data collection by a party, such as a satellite operator interested in "sniffing" satellites and determining whether they are tumbling, or by a general program that simply characterizes all RSOs, or all RSOs with particular characteristics. In the embodiment described herein, radar array 10 collects data from RSOs, such as RSO 16, that are typically in low Earth orbit. Low Earth orbit is typically defined as an altitude between 160 and 2,000 kilometers (99 and 1,200 miles) above the Earth's surface. In other examples, radar array 10 may alternatively or additionally collect data from RSOs in other orbits.

[0018] In embodiments herein, the data collected by the radar array 10 enables the determination of an object's tumbling characteristics. In most of the embodiments herein, the determined tumbling characteristics take the form of the object's rotational velocity, although other characteristics are possible.

[0019] In the embodiment illustrated in FIG. 1 , the determination of the rotation rate of the RSO 16 is performed by one or more processors 12 of the signal analysis system 18. In FIG. 1 , the signal analysis system 18 is shown external to the radar system 10. However, in other examples, the processor 12 may comprise a processor resident in the radar system 10, one or more processors resident in the radar element itself. The one or more processors 12 execute code that enables the determination of the tumbling characteristics of a space object, such as the RSO 16. The signal analysis system 18 may further comprise a database 14. The database 14 may store various rotation profiles used to determine the rotation of a target, and / or possibly the identification of the object by its geometric shape, and the resulting tumbling characteristics. Like the processor 12, the database 14, if used, may alternatively reside in the radar system 10 or may be an external database. In some examples, the signal analysis system 18 may reside in a radar control module of the radar system 10. In an embodiment, the processes of Figures 2, 3 and 14 may be performed by processor 12 executing suitable code.

[0020] FIG. 2 shows a flow chart of one embodiment of an overall process for determining the tumbling characteristics of resident space objects (RSOs), also referred to herein as objects. In the illustrated embodiment of FIG. 1, radar array 10 tracks most, if not all, of the objects in low Earth orbit. The collected radar data for each of the objects is then analyzed, and the object undergoing analysis is referred to as a target. The process illustrated in FIG. 2 determines the tumbling characteristics of a single target object 16. It will be understood that the process of FIG. 2 may be performed separately for multiple different target objects, either sequentially or simultaneously, by signal analysis system 18.

[0021] The process begins by selecting a target object to analyze in collected radar data and acquiring or obtaining the geometry of the target object 16 at block 30. In the example, the process determines the geometry of the target 16 by estimation based on radar data collected by the radar array 10 about the target 16, commonly referred to as the radar signature of the target 16. As used herein, the geometry of the target 16 includes at least one physical dimension of the target 16, e.g., the one-, two-, or three-dimensional extent of the target 16, or the detailed three-dimensional shape of the target 16.

[0022] In other examples, the geometry of the target object may be obtained in other ways. In some examples, additional sensors may be used in addition to the radar array 10 to determine the geometry of the target object. In some examples, the identity of the target object may be confirmed, and stored information regarding the geometry of the target object may be obtained. For example, a satellite operator may provide details of satellites of interest, including their ephemeris and geometry details, so that the system can identify the satellites of interest in the radar data and determine their tumbling characteristics.

[0023] At block 32, the process generates a model of the radar reflectivity characteristics of the target object by synthesizing a random distribution of radio wave scatterers within the target volume based on the target's acquired geometry. An example of a generated radar reflectivity model is shown in Figure 4, where 100 radio wave scatterers are randomly distributed throughout the target object's acquired geometry. In other examples, a different number of radar scatterers may be used.

[0024] Next, in block 34, the process uses a physical modeling process to calculate multiple radar cross sections (RCSs) from at least one generated model instance, but typically from an ensemble of multiple generated model instances, to average out model uncertainties. Because the actual number and location of radio / radar scattering points, such as edges, corners, defects, and screw heads, generally cannot be accurately determined remotely, it is understood that an ensemble of multiple different generated model instances with varying random distributions of radio scatterers is used accordingly to collectively provide a closer approximation of the actual target object. Detailed information about the scattering geometry of the target object is known. In a known example, for example, if the target object is known to be a particular satellite, this information can be incorporated into the scattering model.

[0025] Next, in block 36, multiple RCS autocorrelations are determined for each generated model instance by rotating each model relative to the virtual radar observer and calculating the model's RCS at different angles. It is understood that if the model is rotated beyond the autocorrelation angle, the RCS will appear effectively random. The RCS calculated at different angles for each model is then used to calculate an RCS autocorrelation function for each model. The RCS angular autocorrelation functions are calculated for various parameter values ​​for each model. The different angles may include, for example, different axes of rotation of the model and different angular positions about these axes. Other parameters may also be used. The resulting group of angular autocorrelation functions is referred to herein as an ensemble of functions, and this ensemble of functions can be used to interpret measured RCS signatures, as described in more detail below.

[0026] In block 36, the ensemble of RCS angular autocorrelation functions is analyzed to determine a best-fit curve for the ensemble. The best-fit curve may be an averaged curve; however, other types of best-fit curves may be used. The determined best-fit curve is used as an estimated angular autocorrelation function for the modeled geometry of the target object. FIG. 5 shows a graphical representation of a set of normalized angular autocorrelation curves corresponding to an exemplary group or ensemble of angular autocorrelation functions calculated for the target object, with curve 60 being the best fit for the set of curves. The best-fit curve is an estimate of the expected angular autocorrelation function given the modeled geometry.

[0027] Next, at block 38, the angle at which the estimated angular autocorrelation function corresponding to the best fit curve drops to half its peak value is determined. This angle is identified as the expected decorrelation angle. Advantageously, the identified expected decorrelation angle may be stored in database 14.

[0028] The method for determining the best fit curve of the RCS angular autocorrelation function according to blocks 30-36 can also be considered a Monte Carlo method.

[0029] At block 40, decorrelation times are determined from radar data measurements of the target object. These radar data measurements may be referred to as a radar data set. The determined decorrelation times are then used in conjunction with the expected decorrelation angles determined at block 38 to determine the rotation rate of the target, as described in more detail below. Figure 3 shows a flow diagram of one embodiment of a process for determining decorrelation times that may be used at block 40. Conveniently, the expected decorrelation angles may also be obtained from database 14.

[0030] At block 42, the process calculates the rotation rate of the target. This is done by dividing the determined decorrelation angle by the determined decorrelation time to calculate the rotation rate. It will be appreciated that the RCS decorrelation angle of most objects is relatively small, typically only a few degrees. As a result, in order for the present approach to detect target rotation and determine rotation rate, it is only necessary to track the target with the radar array 10 as it rotates through this small RCS decorrelation angle of only a few degrees during a single pass. In contrast, conventional approaches using light curves or radiometry techniques require the target to be observed over one or more 360-degree rotations.

[0031] The output of block 42 is a determined rotation rate, i.e., angle / time, for the target based on radar data from a single radar data set from a single pass of the target across the field of view of the radar array 10. The determined rotation rate may be stored in database 14.

[0032] For example, if the decorrelation time cannot be determined in block 40 because the radar data measurements of the target object are not changing, it may be concluded that the target object is not rotating. Optionally, in this case, the process may identify the rotation rate as zero, skip block 42, and proceed directly to block 44. This may avoid wasting computing resources in performing block 42 unnecessarily.

[0033] This process is then repeated for multiple independent radar data sets, as indicated by block 44. These independent radar sets may be, for example, radar data from different passes of the target across the radar array 10, or radar data from different portions of the same pass of the target across the radar array 10, or radar data for targets observed by other radar systems. It will be appreciated that target objects may be identified as the same target object for different passes over the same or different radar arrays, for example, by using the collected radar data to determine the orbital path and timing, or ephemeris, of the target object, and assuming that target objects having matching orbital paths and timing, or ephemeris, must be the same object.

[0034] Next, at block 46, the process compares the multiple determined rotation rates of the target output by block 42 based on each of the multiple independent radar data sets and identifies the highest determined rotation rate from all of the multiple independent radar data sets. At block 48, this identified highest determined rotation rate is output as the best fit rotation rate of the target object. The best fit rotation rate may be stored in database 14.

[0035] As described above, for the present approach to detect target rotation and determine a best-fit rotation rate over many passes, the target only needs to be tracked by the radar array 10 as it rotates through a small RCS decorrelation angle of only a few degrees during each pass. In contrast, conventional approaches using light curves or radiometry techniques require observing the target through one or more 360-degree rotations in a single pass to determine the rotation rate. Therefore, the process described herein can make it possible to determine the rotation of objects and their rotation rates in situations where conventional approaches cannot be used. An example of such a situation is an object that is not within the radar's field of view for a sufficient time to complete a 360-degree rotation.

[0036] This best-fit rotational speed of the target object may be used as the tumbling characteristic of the target object. In another example, the best-fit rotational speed of the target object may be compared to one or more thresholds, and whether the best-fit rotational speed of the target object is above or below a predetermined threshold or within a predetermined band may be used as the tumbling characteristic of the target object instead of the actual rotational speed. In another example, a determination of whether the target object is non-rotating (i.e., has no detectable rotation) may be used as the tumbling characteristic of the target object.

[0037] As explained above, the goal of the process is to determine the rotational velocity of the target object. In practice, for a rotating target object, this rotation is about an axis of rotation that, at a given time, forms an angle θ with respect to the line of sight between the radar array 10 and the target object. The direction of the axis of rotation corresponds to the direction of the target object's spin angular momentum vector. The orientation of the target object's axis of rotation / angular momentum vector is generally unknown. Furthermore, because the orientation of the target object's axis of rotation / angular momentum vector is fixed to the object's own coordinate system and not fixed relative to the Earth's surface, the angle θ changes over time as the radar array 10 moves with the Earth's rotation. Therefore, the angle θ may generally have any value for any particular radar data set.

[0038] For geometric reasons, the rotation rate determined from a single radar data set will be the actual "true" rotation rate of the target object for that radar data set multiplied by sin θ, where θ is the angle between the line of sight from the radar array 10 to the target object and the angular momentum vector of the target object. This is because the rate of change of relative range of different radar scattering elements from the radar array 10 as a result of target object rotation, where a change in relative range will cause a change in the target object's RCS, is geometrically multiplied by sin θ for a radar array 10 with a line of sight of angle θ. As a result, the determined decorrelation time is divided by sin θ, and therefore the determined rotation rate is multiplied by sin θ.

[0039] Thus, when the multiple determined rotation rates are compared in block 46, the data set with the highest determined rotation rate provides the determined rotation rate that is closest to the actual true rotation rate of the target object, because this is the data set for which the angle θ between the line of sight to the target and its angular momentum vector is closest to 90 degrees, and therefore the value of sin(θ) is closest to 1. As a result, the determined maximum determined rotation rate determined in block 46 and output in block 48 is the best estimate of the actual rotation rate of the target object.

[0040] Repeating the process for multiple independent radar data sets collected on different passes or by other radar systems, and selecting the highest determined rotation rate in block 46, avoids a potential problem that, for a particular pass of a target object, if the angle θ between the field of view for the target and its angular momentum vector is 0 degrees or very close to 0 degrees, rotation may not produce a measurable change in the RCS of the target object, resulting in a rotating target object being mistaken for a non-rotating object. As noted above, the angle θ changes over time as the radar array 10 moves with the rotating Earth's surface, and the angle θ will, of course, change for radars at different locations, so repeating the process in this manner overcomes this potential problem.

[0041] As mentioned above, Figure 3 shows a flow diagram of one embodiment of a process for determining decorrelation times from radar data measurements of target objects that may be used in block 40 of Figure 2. The flow diagram of Figure 3 shows a process used to determine decorrelation times from radar data measurements of target objects in a single radar data set. If the process of Figure 2 is repeated for multiple independent radar data sets, the process of Figure 3 is repeated for each of the multiple radar data sets. In some examples, alternative methods may be used to determine decorrelation times of target objects.

[0042] 3, at block 50, the process extracts radar cross section (RCS) time series data at regular intervals from the radar data provided by the radar array 10. The RCS time series data is calculated from the signal-to-noise (SNR) ratio of the received radar data.

[0043] Optionally, the extracted RCS time series data may undergo noise suppression at block 52. In one embodiment, noise suppression may involve applying a low-pass filter to the RCS time series data at block 52. However, in some instances, noise suppression may not be necessary.

[0044] Next, in block 54, the RCS time series data, if noise suppressed / filtered in optional block 52, is multiplied by a window and the mean subtracted to form windowed data. By way of example, the window may be, but is not limited to, Gaussian, Hamming, Hann, Sine, or Tukey, although other types of windows may be used in some examples.

[0045] A window is a number of samples. Because the samples that make up the RCS time series data are taken at regular intervals over a period of time, the window represents the period over which the samples were taken. In one embodiment, the autocorrelation function is calculated over a window spanning 11,000 samples of a predetermined length of time, referred to herein as the window length. In one embodiment, the window length is 18 milliseconds, so the samples range from 18 milliseconds to 18 seconds.

[0046] Next, at block 56, the process calculates an autocorrelation function versus time using the windowed RCS time series data from block 54. This autocorrelation function calculation is repeated for an ensemble of different windows, such as various combinations of window lengths and positions within the time series of the RCS time series data, as shown in block 57. The resulting autocorrelation function versus time data is then used to determine the decorrelation time of the target object. In one embodiment, similar to finding the autocorrelation angle, the determined decorrelation time is the consensus time delay at which the autocorrelation falls to a value that is half of its peak value for the ensemble of autocorrelation functions.

[0047] The process of FIG. 2 then uses the determined decorrelation time and the expected decorrelation angle of block 42 to determine the rotation rate.

[0048] The purpose of windowing in block 54 of the process of FIG. 3 is to reduce the analysis' susceptibility to edge effects and slowly varying systematic errors. FIG. 6 shows an example graph of calculated target autocorrelation function values ​​versus delay time. The graph in FIG. 6 shows a grayscale plot of various calculated target autocorrelation function values ​​versus delay time for various window lengths, with darkness in the graph corresponding to more lines plotting the calculated autocorrelation function values ​​versus delay time. The graph in FIG. 6 shows that the autocorrelation function depends on the window size, but at certain values ​​of various window lengths, there is a convergence 60 of the autocorrelation curves, indicating a robust signal corresponding to the actual signal. The process of FIG. 3 interprets this convergence 60 of values ​​as a consensus time delay caused by the rotation of the target object, as described above.

[0049] In the embodiment of Figure 2, blocks 30-38 calculate an expected decorrelation angle for target object 16. This is not required, and in other examples, various methods may be used to determine the expected decorrelation angle for the target object. In such examples, blocks 40-48 may use the expected decorrelation angle to determine rotation rate in a manner similar to the expected decorrelation angle calculated by blocks 30-38.

[0050] In some instances where the identity of the target object is known, it may be possible to calculate the expected decorrelation angle from the target object's structure if the target object's structure and resulting radar reflectivity are known in sufficient detail. The radar reflectivity of the structure may be determined, for example, from structural details provided by the satellite operator or from a satellite design catalog. In other instances where the target object or the same object can be accessed before it is launched, it may be possible to measure the object's radar cross section at various angles and use these measurements to calculate the expected decorrelation angle. Other methods of determining an object's expected decorrelation angle may also be used.

[0051] FIG. 7 shows a graph of calculated decorrelation time as a function of window length used for a single radar data set of a target object. The decorrelation time is calculated from the autocorrelation analysis for each different window length. At short window lengths, there is a steep slope, such as at 64, that leads to noise decorrelation. However, there is typically a broad plateau spanning several seconds, such as at 66, where the actual rotation signal is likely to be present. At longer lengths, i.e., larger window sizes, there is a discontinuity at the end of the time series and a rise due to slowly varying systematic errors. One solution to avoid this rise at larger window sizes is to calculate the rotational autocorrelation over a time window with a window length within the wide plateau, such as a window of less than 12 seconds in the illustrated embodiment. Thus, in some examples, calculating the autocorrelation function versus time using the windowed RCS time series data at block 56 may include determining a range of window lengths that corresponds to this plateau and then calculating the autocorrelation function versus time using window lengths within this range.

[0052] Returning to FIG. 2, the process yields a best-fit rotational velocity for the target object, which can be used as the target object's tumbling characteristic. Having described the process in general terms, the discussion now turns to examples of analysis applied to real data. FIGS. 8-15 illustrate various examples of determining the tumbling characteristics of a space object. FIGS. 8 and 9 illustrate one embodiment of time series processing of signals indicative of a tumbling object and its associated velocity. A time series consists of a series of data points indexed in time. As implemented here, a time series is an array considered to be successive, equally spaced time points. Time series can be subject to edge effects that can skew the results.

[0053] In the top graph of Figure 8, the outer line of the top line is the raw SNR sample data, shown as solid line 70. Dashed line 72 is the raw SNR sample data after low-pass filtering. In this example, low-pass filtering was applied to the autocorrelation of a single object data set. The noise autocorrelation revealed a very narrow but significant spike at zero delay. Low-pass filtering is appropriate in this example because its primary effect is to suppress the zero-delay noise spike; the general autocorrelation shape was unaffected. The bottom dotted line 74 represents windowed data with a window length of approximately 18 seconds. Each graph in Figure 8 represents a single pass of an object over the radar, i.e., across the radar's field of view. The same convention for lines is used in Figures 8, 10, 12, and 14.

[0054] Figure 9 shows the autocorrelation results for the data in Figure 8. Again, each graph represents a single pass of an object over the radar. The convergence of the curves at long decorrelation times is due to edge effects in the time series. Each graph represents a different sample. In the graph below, curve 80 is the autocorrelation curve. Curve 82 represents the theoretical model curve for a rotation rate of 5 degrees / second. Curve 84 shows the model curve at 7 degrees / second, and curve 86 represents the model at 9 degrees / second. The fit of autocorrelation curve 80 to curves 82-86 corresponding to different rotation rates provides an illustrative illustration of how the rotation rate is determined in block 42 of Figure 2. The same convention for lines is used in Figures 9, 11, 13, and 15.

[0055] The top graph in Figure 9 shows a set of samples with the strongest rotation signal consistent with the theoretical model at 7-8 deg / s rotation. The middle graph shows a set of samples with a rotation signal consistent with the theoretical model at 5 deg / s. The bottom graph shows a set of samples with no strong rotation signal below the window envelope. Thus, the strongest indication is that the satellite is rotating at 7-8 deg / s.

[0056] Figure 10 shows the time series for a second object. This set of samples consists of four passes of the object, with each of the four graphs showing samples corresponding to a single pass of the object. Again, the raw sample data is shown as a solid line 70, the low-pass filtered data as a dashed line 72, and the windowed data as a dotted line 74. Figure 11 shows the autocorrelation curves for the time series from Figure 10, with each of the four graphs corresponding to each of the four graphs in Figure 10. Again, the convergence of the curves at long decorrelation times is due to edge effects in the time series. In the graphs in Figure 11, the top graph shows the set of curves containing the 7-9 deg / s signal, the second graph below it shows the set of curves containing the 8-9 deg / s signal, the third graph below it shows the set of curves without a strong rotation signal below the window envelope, and the bottom graph shows the set of curves containing the 7-9 deg / s signal. Therefore, the strongest constraint is the rotation rate of 8-9 deg / s.

[0057] In the above embodiment, the data may optionally be filtered at block 52. In the described embodiment, this optional filtering is low-pass filtering. In other examples, the data may be subjected to low-pass filtering in addition to, or as an alternative to, high-pass filtering. This process appears to isolate frequencies of interest relative to noise, so which type of filtering is appropriate in any particular example may depend on the data set and the equipment used to collect the data.

[0058] Another problem that may occur and be identified may be due to pointing errors seen at high altitudes. Figure 12 shows a time series for a fourth object. This object has data from two passes, and each graph in Figure 12 shows the SNR data from each pass. Figure 13 shows two graphs of the autocorrelation curve for the data in Figure 12. Neither data set provides strong rotation constraints. The hint of signal in the bottom graph is likely due to pointing errors seen at high altitudes. In some instances where the radar array is affected by pointing errors at high altitudes, the system can define the signal level so that the subtle signal in the bottom panel can be used to identify pointing errors. The system can also use additional diagnostics from the radar equipment to identify data that may be affected by pointing errors. In some instances, the identified data may not then be used to calculate rotation rates. This may eliminate irrelevant data sets and inaccuracies due to pointing errors at high altitudes.

[0059] The above-described embodiments allow the rotational velocity of an orbiting object to be determined and used as a tumbling characteristic of the object.

[0060] In further embodiments, the process can be extended to determine the orientation of the angular momentum vector of the rotating object, i.e., the orientation of the axis of rotation. In some examples, the orientation of the angular momentum vector may then be used as the tumbling property of the object. In other examples, both the rotation rate and the orientation of the angular momentum vector may be used as the tumbling property of the object.

[0061] Figure 14 shows a flow chart of a further embodiment of an overall process for determining the tumbling characteristics of a resident space object (RSO) that can determine the orientation of the rotating object's angular momentum vector. The process of the embodiment of Figure 14 may be implemented, for example, by the process of the embodiment of Figures 2 and 3 by the system of Figure 1. Figure 15 is an exemplary diagram illustrating some of the concepts of the process of Figure 14.

[0062] The process of FIG. 14 begins by selecting a target object 16, analyzing the collected radar data, and detecting the target object 16 at time t n , the process begins by determining the position of the target object 16 at . Next, at block 92, the process determines the line of sight vector between the radar array 10 and the target object 16.

[0063]

number

[0064] Determine the time t n The line of sight vector between the radar array 10 and the target object 16 at

[0065]

number

[0066] It can be specified as follows. This process is repeated at least three times, as indicated by block 94. The number of times that blocks 92 and 94 are repeated may vary in different examples, so long as the number is at least three. An example of this is illustrated in Figure 15, where the position of target object 16 at first, second, third, and fourth different times t1, t2, and t3 are represented by the line of sight vectors calculated at times t1, t2, and t3, respectively.

[0067]

number

[0068] It is shown together with. Next, in block 96, at least three different times t n, the rotational speeds of the target object 16 at each of the points 16a, 16b, 16c, 16d, 16e, 16f, 16g, 16h, 16i, 16j ...

[0069] Next, at block 98, the process calculates the angular momentum vector

[0070]

number

[0071] A description of one way this can be done is as follows: Line of sight vector between radar array 10 and target object 16

[0072]

number

[0073] can be rewritten in terms of Cartesian unit vectors as follows:

[0074]

number

[0075] Equation (1) can be rewritten to define a unit vector as follows:

[0076]

number

[0077] During the ceremony,

[0078]

number

[0079] Each calculated rotation rate w(t) corresponds to the line of sight vector between the radar array 10 and the target object 16 at time t.

[0080]

number

[0081] Angular rotation vector into a plane perpendicular to

[0082]

number

[0083] The calculated angular velocity is therefore:

[0084]

number

[0085] In the formula, x represents the cross product. As explained above, the calculated angular rotation rate can also be expressed as:

[0086]

number

[0087] where θ(t) is the line of sight vector

[0088]

number

[0089] and the angular rotation vector

[0090]

number

[0091] is the angle between From the above formula, it can be shown as follows:

[0092]

number

[0093] During the ceremony,

[0094]

number

[0095]

number

[0096] The ingredients are:

[0097]

number

[0098] is. Therefore, the result is as follows:

[0099]

number

[0100] The values ​​of x(t), y(t), and z(t) are known from determining the position of the target object 16 in block 90, and w(t) is the measured and calculated projected rotation velocity obtained in block 96. Thus, equation (7) contains three unknowns, w x , w y , and w z where, as shown in equation (6),

[0101]

number

[0102] is. Thus, if at least three measurements were taken at different times in block 90, then in block 98 the process solves equation (7) to find w x , w y , and w z ,therefore,

[0103]

number

[0104] can be determined. At block 98, the solution can be found by nonlinear least squares fitting. In other examples, equation (7) can be solved in other ways.

[0105] Block 98 decided

[0106]

number

[0107] is then output in block 100. Fitted w x , w y , and w z The uncertainty in the value, therefore,

[0108]

number

[0109] The uncertainty in is the uncertainty of the line of sight vector measured at different times.

[0110]

number

[0111] The accuracy of the accuracy depends on the angular spread of the line of sight vectors, the greater the angular spread of the line of sight vectors, the greater the certainty. In practice, for a single pass of the RSO 16 over the radar array 10, there are several different times t n Different gaze vectors measured at

[0112]

number

[0113] are generally nearly coplanar, so the angular rotation vector

[0114]

number

[0115] The estimate of t may be highly uncertain and may be inaccurate. Therefore, the estimate of t may be highly uncertain and may be inaccurate at different times t during different passes of the RSO 16 over the radar array 10. n Gaze vector in

[0116]

number

[0117] It may be desirable to measure the n gaze vector of

[0118]

number

[0119] are non-coplanar, so the angular rotation vector

[0120]

number

[0121] can be accurately estimated. Generally, the angular rotation vector

[0122]

number

[0123] The accuracy with which the rotation speed estimate w(t n ) accuracy, gaze vector

[0124]

number

[0125] Accuracy and gaze vector

[0126]

number

[0127] It depends on the dispersion of the pointing direction. In some instances, the process of FIG.

[0128]

number

[0129] The value of may be used to refine the best fit rotational velocity of the target object identified in block 46 of FIG.

[0130] In the embodiment of the process of FIG. 14 described above, time t n The position of the target object 16 at time t is determined from the radar data, and this position is then used to n The line of sight vector between the radar array 10 and the target object 16 at

[0131]

number

[0132] It is not necessary to use such a two-step process. In some instances, the line-of-sight vectors may be determined directly from radar data.

[0133]

number

[0134] It may also be possible to determine In the manner described above, systems and methods can determine and track the tumbling characteristics of resident space objects, such as satellites. As explained above, the determined tumbling characteristics may include either or both of the object's rotational rate and the orientation of the object's angular momentum vector. Identifying the tumbling characteristics of a space object may make it possible to determine the correct functioning of the satellite by determining whether those tumbling characteristics match their intended or desired tumbling characteristics. For example, if a satellite is intended not to rotate and is determined to be rotating at any measurable rate, this may indicate some malfunction or failure of the satellite. Furthermore, identifying the tumbling characteristics of a space object may make it possible to determine the life cycle stage of the satellite. Furthermore, identifying the tumbling characteristics may enable system operators to notify object owners that their objects are tumbling or to adjust the paths of other objects to allow them to pass between them.

[0135] The above methods may also further enable the development of an object's tumbling characteristic profile that can be used as an object signature to enable identification and / or long-term tracking of the object. The determined tumbling characteristic of a resident space object may be determined and recorded as an object fingerprint. Identification of a detected resident space object having tumbling characteristics that match the fingerprint may then indicate that the detected resident space object is the same resident space object.

[0136] In some examples, measured or determined parameters of a resident space object may be used instead of or in addition to the determined tumbling characteristics themselves, such as RCS characteristics or metrics, or calculated decorrelation times, or any number of these parameters may be recorded as at least part of the Object Fingerprint.

[0137] Fingerprint matching may be performed as a probabilistic assessment, whereby a number of different parameter values ​​of a resident space object are determined and compared with corresponding parameter values ​​of the recorded object fingerprint, and based on how similar each of the determined parameters is to the recorded values ​​of the fingerprint, the probability that the resident space object is the fingerprinted object can be assessed.

[0138] In some examples, the Fingerprint of Things may include a tumbling characteristic profile and other parameters of the resident space object. In one example, the Fingerprint of Things may include a tumbling characteristic profile and an orbital path of the resident space object.

[0139] In one example, fingerprinting may include comparing the determined at least one tumbling characteristic of the object to a stored record of previously determined at least one tumbling characteristic of a previous object, and determining whether the object and the previous object are the same object based on the comparison.

[0140] The above-described embodiments use windowing of the RCS time series data to remove edge effects. In some instances where the RCS time series data has suitable properties, this windowing may not be necessary.

[0141] While the above-described embodiments comprise one or more fixed radar arrays on Earth, in other examples some or all of the radars may be mounted on one or more moving platforms such as satellites or on the surface of other bodies.

[0142] The embodiments described above use one or more radar arrays. In other examples, other types of radar that are not arrays may be used.

[0143] While the above-described embodiments characterize RSOs in Earth orbit, other examples may characterize objects in orbit around other celestial bodies.

[0144] The above-described embodiments include several examples in which radar data acquired at different passes of an object and / or from different radar devices is used. It will be understood that the system includes suitable storage and communication devices to make this possible. Many configurations for storing and transmitting data are well known to those skilled in the art and need not be described herein.

[0145] The above-described embodiments, for simplicity and clarity, describe how the tumbling characteristics of a single object may be determined. It will be appreciated that the described process may be repeated for multiple different objects in order to determine the tumbling characteristics of each of the multiple different objects.

[0146] In the above embodiments, some functions may be provided by software, whereas in other examples, the functions may be provided wholly or partly in hardware, for example by dedicated electronic circuitry.

[0147] In the above embodiments, the system may be implemented as any form of computing and / or electronic device. Such devices may include one or more processors, which may be microprocessors, controllers, or any other suitable type of processor, for processing computer-executable instructions that control the operation of the device to collect and record routing information. In some examples, for example, when a system-on-chip architecture is used, the processor may include one or more fixed function blocks (also referred to as accelerators) that implement portions of the method in hardware (rather than software or firmware). Platform software, including an operating system or any other suitable platform software, may be provided with the computing-based device to enable application software to be executed on the device.

[0148] Computer programs and computer-executable instructions may be provided using any computer-readable medium accessible by a computing-based device. Computer-readable media may include, for example, computer storage media such as memory and communication media. Computer storage media, such as memory, may include volatile and non-volatile, removable or non-removable media implemented in any method or technology for storage of information such as computer-readable instructions, data structures, program modules, or other data. Computer storage media include, but are not limited to, RAM, ROM, EPROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVDs) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that may be used to store information for access by a computing device. In contrast, communication media may embodi computer-readable instructions, data structures, program modules, or other data in a modulated data signal such as a carrier wave or other transport mechanism. As defined herein, computer storage media does not include communication media.

[0149] Although the system is shown as a single device, the system may be distributed or located remotely and accessed via a network or other communications link (eg, using a communications interface).

[0150] The term "computer" is used herein to refer to any device with processing capability such that it is capable of executing instructions. Those skilled in the art will understand that such processing capability is incorporated into many different devices, and thus the term "computer" includes PCs, servers, mobile phones, personal digital assistants, and many other devices.

[0151] Those skilled in the art will recognize that storage devices utilized to store program instructions can be distributed across a network. For example, a remote computer can store an example of a process described as software. A local or terminal computer can access a remote computer and download some or all of the software to execute the program. Alternatively, a local computer can download portions of the software as needed, or execute some software instructions at a local terminal and some at a remote computer (or computer network). Those skilled in the art will also recognize that, utilizing conventional techniques known to those skilled in the art, all or a portion of the software instructions may be executed by dedicated circuitry, such as a DSP, programmable logic array, or the like.

[0152] It will be understood that the benefits and advantages described above may pertain to one embodiment or to several embodiments, and the embodiments are not limited to those that solve any or all of the stated problems or have any or all of the stated benefits and advantages.

[0153] Any reference to "an" item refers to one or more of those items. The term "comprising" is used herein to mean including specified method steps or elements, but such steps or elements do not comprise an exclusive list and a method or apparatus may include additional steps or elements.

[0154] The ordering of steps in the methods described herein is exemplary, but these steps may be performed in any suitable order, or simultaneously as appropriate. Additionally, steps may be added or substituted, or individual steps may be deleted, from any of these methods without departing from the scope of the subject matter described herein. Aspects of any of the above-described examples may be combined with aspects of any of the other examples described to form further examples without losing the desired effect.

[0155] It will be understood that the above description of preferred embodiments is given by way of example only, and that various modifications may be made by those skilled in the art. While various embodiments have been described above with a certain degree of particularity, or with reference to one or more individual embodiments, those skilled in the art may make numerous modifications to the disclosed embodiments without departing from the spirit or scope of the invention.

Claims

1. 1. A method for determining at least one tumbling characteristic of an object, comprising: Obtaining an expected radar cross section (RCS) decorrelation angle of the object; acquiring radar data of the object; determining a decorrelation time of the RCS of the object from the radar data of the object; and determining at least one tumbling characteristic of the object using the obtained decorrelation angle and the determined decorrelation time.

2. Determining the at least one tumbling characteristic of the object using the obtained decorrelation angle and the determined decorrelation time includes: The method of claim 1 , comprising dividing the obtained decorrelation angle by the determined decorrelation time to determine the at least one tumbling characteristic of the object.

3. The method of claim 1 or 2, wherein the at least one tumbling characteristic of the object comprises a rotational speed of the object.

4. repeating the acquiring of radar data, determining a decorrelation time, and dividing the acquired decorrelation angle by the determined decorrelation time for a plurality of radar data sets of the object to determine a rotation rate of the object for each of the plurality of radar data sets; identifying a maximum rotation rate from the determined rotation rates for each of the plurality of radar data sets; The method of claim 3 , further comprising: using the identified maximum rotational speed as the rotational speed of the object.

5. The method of claim 4 , wherein the plurality of radar data sets comprises radar data sets acquired at different orbital passes of the object.

6. The method of claim 4 or 5, wherein the plurality of radar data sets comprises radar data sets obtained from different radars.

7. Obtaining the RCS decorrelation angle of the object comprises: obtaining a geometric shape of the object; generating a distribution of radio wave scatterers within the geometric shape of the object; calculating a plurality of radar cross sections at different angles for the distribution; determining an autocorrelation as a function of angle for each of the plurality of radar cross sections; and determining a decorrelation angle using the autocorrelation.

8. repeating the steps of generating a distribution of radio wave scatterers, calculating a plurality of radar cross sections, and determining the autocorrelation as a function of angle for each of a plurality of different distributions of scatterers; The method of claim 7 further comprising: determining a decorrelation angle using the autocorrelation.

9. The method of claim 7 or 8, wherein obtaining the geometry of a target object comprises obtaining at least one physical dimension of the object or a three-dimensional shape of the object.

10. The method according to any one of claims 1 to 6, wherein the RCS decorrelation angle of the object is obtained by simulation or measurement.

11. Determining the decorrelation time includes: extracting time series data from the radar data of the object; defining a time window for the time series data to generate windowed data; determining an autocorrelation function versus time from a combination of a plurality of positions in the time series of the windowed data and a window length of the time window; and using the autocorrelation function versus time to determine the decorrelation time as the time delay at which the autocorrelation falls to half of its peak value.

12. defining a plurality of different time windows for the time series data to generate a plurality of windowed data; determining an autocorrelation function versus time for each of the plurality of windowed data to generate a plurality of autocorrelation functions versus time; analyzing the plurality of autocorrelation functions versus time to identify a time delay at which the autocorrelation has half the peak value of the autocorrelation; The method of claim 10 , further comprising: identifying the time delay as the determined decorrelation time.

13. applying noise suppression to the time series data to generate filtered time series data; The method of claim 10 or 11, further comprising: using the filtered time series data to generate the windowed data.

14. The method of claim 12 , wherein applying noise suppression comprises applying at least one of a low-pass filter, a high-pass filter.

15. Defining a time window on the time series data to generate windowed data includes: defining a period for the time series data; multiplying the time series data within the time period by a window; and subtracting an average of a product of the window and the time series data within the period from the product to generate the windowed data.

16. repeating the acquiring of radar data, determining a decorrelation time, and dividing the acquired decorrelation angle by the determined decorrelation time at different times for a plurality of radar data sets of the object to determine a rotational velocity of the object as the tumbling characteristic for each of at least three radar data sets; determining a line-of-sight vector between the object and a radar collecting the radar data at each of the different times; 16. A method according to claim 3, or any one of claims 4 to 15 when directly or indirectly dependent on claim 3, further comprising: solving a set of linear equations for each of the different times to determine an angular rotation vector of the object.

17. The linear equation for each of the different times is: [Equation 1] In the formula, w(t ｎ ) is the time t ｎ is the measured rotational velocity of the target object at [Equation 2] is the time t ｎ the line of sight vector between the radar and the target object at [Equation 3] The method of claim 16 , wherein: is the angular rotation vector of the body.

18. 18. The method of claim 16 or 17, wherein the at least one tumbling characteristic comprises the angular rotation vector of the object.

19. The method of any preceding claim, wherein the at least one tumbling characteristic of the object comprises whether the object is rotating.

20. The method of any one of claims 1 to 19, wherein the object is a Resident Space Object (RSO).

21. comparing the determined at least one tumbling characteristic of the object with a stored record of previously determined at least one tumbling characteristic of a previous object; The method of any one of claims 1 to 20, further comprising: determining whether the object and the previous object are the same object based on the comparison.

22. A system comprising a processor configured to perform the method of any one of claims 1 to 20.

23. 1. A system for determining a tumbling characteristic of an object, comprising: a radar for acquiring radar data of each object; and at least one processor configured to execute code, the code causing the processor to: Obtaining the decorrelation angle of the radar cross section (RCS) of each object; determining a decorrelation time of the RCS of each object from the radar data of the object; The system uses the obtained decorrelation angle and the determined decorrelation time for each object to determine at least one tumbling characteristic of the object.

24. The processor is further configured to execute code that causes the processor to:

24. The system of claim 23, wherein the at least one tumbling characteristic of the objects is determined by using the obtained decorrelation angle and the determined decorrelation time for each object and dividing the obtained decorrelation angle by the determined decorrelation time to determine the at least one tumbling characteristic of the object.

25. 24. The system of claim 23, wherein the at least one tumbling characteristic of each object comprises a rotational speed of the object.

26. The system of any one of claims 23 to 25, further comprising a database.

27. 26. The system of claim 25, wherein the processor is further configured to execute code to access the database to obtain decorrelation angles of the RCS of each object.

28. 27. The system of claim 26, wherein the processor further executes code to store the rotational speed of the object associated with the object in the database.

Images