Radio frequency interference data fusion method for one-dimensional and two-dimensional synthetic aperture radiometers

Abstract

The application discloses a radio frequency interference data fusion method and device of one-dimensional and two-dimensional synthetic aperture radiometers, and the method comprises the following steps: obtaining measurement data of the one-dimensional synthetic aperture radiometer, and determining first parameters of each radio frequency interference (RFI) source; obtaining measurement data of the two-dimensional synthetic aperture radiometer, and determining second parameters of each radio frequency interference (RFI) source; constructing a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model, wherein the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model comprises a plurality of objective functions; determining third parameters corresponding to all radio frequency interference sources, forming fusion data of radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometers; and comparing the performance of the measurement data and the fusion data. The method of the application is suitable for simultaneously carrying the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer to carry out RFI detection, RFI positioning and RFI suppression data processing.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of satellite microwave ocean remote sensing, and particularly relates to a radio frequency interference data fusion method of one-dimensional and two-dimensional synthetic aperture radiometers. BACKGROUND

[0002] Ocean salinity detection technology is a technology for detecting ocean salinity radiance brightness temperature and inversely obtaining ocean salinity. The best frequency band for seawater salinity detection is the L frequency band. In order to achieve high-precision and high-spatial-resolution detection of seawater salinity radiance brightness temperature, the synthetic aperture radiometer system is the first choice. Synthetic aperture radiometers are usually divided into one-dimensional synthetic aperture radiometers and two-dimensional synthetic aperture radiometers.

[0003] Radio frequency interference (RFI) is a radio frequency interference signal generated by human beings on the ground. In the process of ocean salinity detection by the SMOS satellite represented by the first two-dimensional synthetic aperture radiometer (MARIS) carried by the satellite, it is found that RFI is the main error source. Therefore, scholars from all over the world have carried out in-depth research and application of RFI on the data of the two-dimensional synthetic aperture radiometer system. The research and application directions include RFI detection, RFI positioning and RFI suppression. Due to the complexity of the signal type, energy intensity and spatial position distribution of RFI, the effects of RFI detection, RFI positioning and RFI suppression have long been unable to meet the error requirements.

[0004] In view of the defects and deficiencies of low detection accuracy of the two-dimensional synthetic aperture radiometer and low spatial resolution of the one-dimensional synthetic aperture radiometer under the same antenna aperture condition, it is necessary to use the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer for joint detection. The two synthetic aperture radiometer systems complement each other, thereby realizing high-precision and high-resolution ocean salinity detection. Through joint detection and fusion data processing of the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer, the RFI detection accuracy, the RFI positioning accuracy and the RFI suppression residual error reduction ability are improved. SUMMARY

[0005] Therefore, the present application provides a radio frequency interference data fusion method of one-dimensional and two-dimensional synthetic aperture radiometers, which can improve the poor quality of RFI detection, RFI positioning and RFI suppression of the prior art one-dimensional synthetic aperture radiometer and two-dimensional synthetic aperture radiometer.

[0006] In order to solve the above technical problems, the present application is implemented as follows.

[0007] A radio frequency interference data fusion method of one-dimensional and two-dimensional synthetic aperture radiometers, the method comprising the following steps:

[0008] Step S1: the satellite simultaneously carries a one-dimensional synthetic aperture radiometer and a two-dimensional synthetic aperture radiometer to measure the same region; obtaining measurement data of the one-dimensional synthetic aperture radiometer to determine first parameters of each radio frequency interference (RFI) source; obtaining measurement data of the two-dimensional synthetic aperture radiometer to determine second parameters of each radio frequency interference (RFI) source;

[0009] Step S2: based on the first parameters and the second parameters, a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model is constructed, the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model includes a plurality of objective functions; third parameters corresponding to all radio frequency interference sources are determined to form fusion data of radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometers;

[0010] Step S3: comparing the performance of the measurement data and the fusion data.

[0011] Preferably, in the step S1, obtaining the measurement data of the one-dimensional synthetic aperture radiometer to determine the first parameters of each radio frequency interference (RFI) source includes:

[0012] Step S111: obtaining measurement data of the one-dimensional synthetic aperture radiometer to form n groups of data;

[0013] Step S112: determining the first parameters of each radio frequency interference (RFI) source, the first parameters including a first detection position and first energy intensity data of each radio frequency interference source, wherein each radio frequency interference source has n groups of first parameters, and in the i-th group of first parameters of the k-th radio frequency interference source, the first detection position is for the value of the i-th group of data ξ in the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, for the value of the i-th group of data η in the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, and the first energy intensity data is for the i-th group of energy intensity data of the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, the first energy intensity data is used to represent the energy intensity of the radio frequency interference source.

[0014] Preferably, in the step S1, obtaining the measurement data of the two-dimensional synthetic aperture radiometer to determine the second parameters of each radio frequency interference (RFI) source includes:

[0015] Step S121: obtaining measurement data of the two-dimensional synthetic aperture radiometer to form n groups of data;

[0016] Step S122: determining second parameters of each radio frequency interference (RFI) source, the second parameters including a second detection position, second energy intensity data of each radio frequency interference source, wherein each radio frequency interference source has n sets of second parameters, and the second detection position in the i-th set of second parameters of the k-th radio frequency interference source is is a value of the i-th set of data ξ of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, is a value of the i-th set of data η of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, is the i-th set of energy intensity data of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, and the second energy intensity data is used to characterize the energy intensity of the radio frequency interference source.

[0017] Preferably, the step S2, wherein the fusion processing model of the one-dimensional and two-dimensional synthetic aperture radiometers includes a fusion processing position objective function and a fusion processing energy intensity objective function; the fusion processing position objective function includes a fusion processing ξ dimension objective function and a fusion processing η dimension objective function.

[0018] The fusion processing ξ dimension objective function is

[0019] When the fusion processing ξ dimension objective function takes the minimum value, the parameter ξ corresponding thereto is the optimal estimation value of the position of the RFI source in the ξ dimension, respectively, the variance of the ξ dimension measurement data of the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer, and n is the number of data sets.

[0020] The fusion processing η dimension objective function is

[0021] When the fusion processing η dimension objective function takes the minimum value, the parameter η corresponding thereto is the optimal estimation value of the position of the RFI source in the η dimension, is the measurement data of the k-th RFI source in the η dimension by the two-dimensional synthetic aperture radiometer, is the variance of the measurement data of the two-dimensional synthetic aperture radiometer;

[0022] The fusion processing energy intensity objective function is as follows:

[0023]

[0024] Wherein, for each radio frequency interference source, the function value of the fusion processing energy intensity objective function is calculated based on the corresponding n sets of first parameters and n sets of second parameters, and when the fusion processing energy intensity objective function takes the minimum value, the parameter T corresponding thereto is the optimal estimation value of the energy intensity of the radio frequency interference source; measurement data of the radio frequency interference source about energy intensity by the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer respectively, variances of the measurement data of the radio frequency interference source about energy intensity by the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer respectively;

[0025] the third parameter corresponding to all radio frequency interference sources is determined, forming fusion data of radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometers, including:

[0026] based on the optimal estimation value of the position of each radio frequency interference source and the optimal estimation value of the energy intensity, the third parameter corresponding to all radio frequency interference sources is determined, including fusion data first optimal point ξ opt , second optimal point η opt and third optimal point T opt ;

[0027]

[0028]

[0029] wherein, is the value of the i-th group of data ξ is the value of the i-th group of data η is the i-th group of energy intensity data of the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, is the value of the i-th group of data ξ is the value of the i-th group of data η is the i-th group of energy intensity data of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, variances of the measurement data of the radio frequency interference source about energy intensity by the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer respectively, is the variance of the η measurement data of the two-dimensional synthetic aperture radiometer,

[0030] measurement data of the radio frequency interference source about energy intensity by the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer respectively.

[0031] The first root mean square error T_RMSE:

[0032]

[0033] wherein, N represents the total number of measurement data, is the energy intensity of the radio frequency interference source detected based on the i-th group of data;

[0034] the second root mean square error P_RMSE,

[0035]

[0036] wherein, is the positioning position of the radio frequency interference source detected based on the i-th group of data, x represents the position of the radio frequency interference source determined based on the optimal fusion data, and the value is the optimal estimation value of the position of the RFI source in the ξ dimension or the optimal estimation value of the position of the RFI source in the η dimension;

[0037] the third root mean square error RE,

[0038] RE=T mitigated -T raw

[0039] wherein, T mitigated represents the suppressed brightness temperature determined based on the optimal fusion data, and T raw represents the original scene brightness temperature.

[0040] The application provides a radio frequency interference data fusion device of a one-dimensional and two-dimensional synthetic aperture radiometer, and the device comprises:

[0041] A data acquisition module is configured to simultaneously carry the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer on a satellite, measure the same region, acquire the measurement data of the one-dimensional synthetic aperture radiometer, determine the first parameters of each radio frequency interference (RFI) source, acquire the measurement data of the two-dimensional synthetic aperture radiometer, and determine the second parameters of each radio frequency interference (RFI) source.

[0042] A fusion module is configured to construct a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model based on the first parameters and the second parameters, wherein the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model comprises a plurality of objective functions, determine the third parameters corresponding to all radio frequency interference sources, and form the fusion data of the radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometers.

[0043] A comparison module is configured to compare the performance of the measurement data and the fusion data.

[0044] The application provides a computer readable storage medium, and the storage medium stores a plurality of instructions; the plurality of instructions are used for loading and executing the method by the processor.

[0045] The electronic device provided by the application has the characteristics that the electronic device comprises:

[0046] A processor is configured to execute a plurality of instructions.

[0047] A memory is configured to store the plurality of instructions.

[0048] The plurality of instructions are configured to be stored in the memory and loaded and executed by the processor according to the method described above.

[0049] The application has the following beneficial technical effects:

[0050] (1) Compared with the SMOS satellite launched internationally, the application adds the one-dimensional synthetic aperture radiometer RFI detection, RFI positioning and RFI suppression capabilities, and compared with the Aquarius, the application adopts the one-dimensional and two-dimensional synthetic aperture radiometer system to improve the RFI detection and RFI positioning capabilities.

[0051] (2) The application is suitable for simultaneously carrying the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer to carry out the RFI detection, RFI positioning and RFI suppression data processing.

[0052] (3) The application can improve the RFI detection and RFI positioning capabilities and reduce the RFI suppression residual through the fusion processing of the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer. BRIEF DESCRIPTION OF DRAWINGS

[0053] Figure 1 A one-dimensional and two-dimensional synthetic aperture radiometer radio frequency interference data fusion method flowchart is provided for the application.

[0054] Figure 2 Another fusion method of the one-dimensional and two-dimensional synthetic aperture radiometer radio frequency interference data is provided for the application.

[0055] Figure 3 A Y-shaped array original scene brightness temperature image and a simulation scene brightness temperature image after adding an RFI source are provided for the application.

[0056] Figure 4 A one-dimensional array original scene brightness temperature image and a simulation scene brightness temperature image after adding an RFI source are provided for the application.

[0057] Figure 5 A Y-shaped array RFI suppression image and an RFI suppression image using data fusion results are provided for the application.

[0058] Figure 6The Y-shaped array brightness temperature estimation residual error map provided by the application and the brightness temperature estimation residual error map using the data fusion result;

[0059] Figure 7 The one-dimensional array RFI suppression image provided by the application and the RFI suppression image using the data fusion result;

[0060] Figure 8 The one-dimensional array brightness temperature estimation residual error map provided by the application and the brightness temperature estimation residual error map using the data fusion result;

[0061] Figure 9 The structure schematic diagram of the radio frequency interference data fusion device of the one-dimensional and two-dimensional synthetic aperture radiometer provided by the application. DETAILED DESCRIPTION

[0062] The application will be described in detail below in combination with the drawings and embodiments.

[0063] As Figures 1-2 shown, the application provides a radio frequency interference data fusion method of one-dimensional and two-dimensional synthetic aperture radiometer, which comprises the following steps:

[0064] Step S1: simultaneously carrying one-dimensional synthetic aperture radiometer and two-dimensional synthetic aperture radiometer on a satellite to measure the same region; obtaining the measurement data of the one-dimensional synthetic aperture radiometer to determine the first parameters of each radio frequency interference (RFI) source; obtaining the measurement data of the two-dimensional synthetic aperture radiometer to determine the second parameters of each radio frequency interference (RFI) source;

[0065] Step S2: based on the first parameters and the second parameters, constructing a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model, the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model comprising a plurality of objective functions; determining the third parameters corresponding to all radio frequency interference sources to form the fusion data of the radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometer;

[0066] Step S3: comparing the performance of the measurement data and the fusion data.

[0067] The step S1, obtaining the measurement data of the one-dimensional synthetic aperture radiometer to determine the first parameters of each radio frequency interference (RFI) source, comprises:

[0068] Step S111: obtaining the measurement data of the one-dimensional synthetic aperture radiometer to form n groups of data;

[0069] Step S112: determining the first parameters of each radio frequency interference (RFI) source, the first parameters comprising the first detection position and the first energy intensity data of each radio frequency interference source, wherein each radio frequency interference source has n groups of first parameters, and in the i-th group of first parameters of the k-th radio frequency interference source, the first detection position is Xi(k) is the value of the i-th dimension of the data of the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, Xi(k) is the value of the i-th dimension of the data of the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, Xi(k) is the i-th energy intensity data of the k-th radio frequency interference source in the field of view of the one-dimensional synthetic aperture radiometer, and the first energy intensity data is used to represent the energy intensity of the radio frequency interference source.

[0070] In this embodiment, the method for determining the first parameters of each radio frequency interference (RFI) source from the measurement data of the one-dimensional synthetic aperture radiometer is to add the position data of the RFI source and the energy intensity data of the RFI source in the earth environment through simulation analysis, or to obtain the actual position data of the RFI source and the energy intensity data of the RFI source by adding the RFI source in the ground experimental environment, or to obtain the position data of the RFI source and the energy intensity data of the RFI source in the earth environment through the on-orbit measurement of the one-dimensional synthetic aperture radiometer of the ocean salinity remote sensing satellite.

[0071] The step S1 of obtaining the measurement data of the two-dimensional synthetic aperture radiometer and determining the second parameters of each radio frequency interference (RFI) source comprises:

[0072] The step S121 of obtaining the single snapshot measurement data and the multi-snapshot measurement data of the two-dimensional synthetic aperture radiometer and forming n groups of data;

[0073] The step S122 of determining the second parameters of each radio frequency interference (RFI) source, wherein the second parameters comprise the second detection position and the second energy intensity data of each radio frequency interference source, and each radio frequency interference source has n groups of second parameters, in the i-th group of second parameters of the k-th radio frequency interference source, the second detection position is Xi(k) is the value of the i-th dimension of the data of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, Xi(k) is the value of the i-th dimension of the data of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, Xi(k) is the i-th energy intensity data of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, and the second energy intensity data is used to represent the energy intensity of the radio frequency interference source.

[0074] In this embodiment, the two-dimensional synthetic aperture radiometer is in the form of a two-dimensional array, specifically a Y-shaped array. Each group of data obtains a group of parameters of each radio frequency interference source.

[0075] In the embodiment, the method for determining the second parameters of each RFI source from the measurement data of the two-dimensional synthetic aperture radiometer is to add the position data of the RFI source and the energy intensity data of the RFI source in the earth environment through simulation analysis, or to obtain the position data of the RFI source and the energy intensity data of the RFI source by adding the RFI source in the ground experimental environment, or to obtain the position data of the RFI source and the energy intensity data of the RFI source by the two-dimensional synthetic aperture radiometer of the ocean salinity remote sensing satellite in the on-orbit measurement.

[0076] The step S2, wherein the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model comprises a fusion processing position target function and a fusion processing energy intensity target function; the fusion processing position target function comprises a fusion processing ξ-dimensional target function and a fusion processing η-dimensional target function.

[0077] The fusion processing ξ-dimensional target function is

[0078] When the fusion processing ξ-dimensional target function takes the minimum value, the parameter ξ corresponding to the minimum value is the optimal estimation value of the position of the RFI source in the ξ dimension, respectively, are the variances of the measurement data of the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer, and n is the number of data groups.

[0079] The one-dimensional synthetic aperture radiometer can only perform synthetic aperture measurement in the η dimension, so when η is estimated, only the measurement data of the two-dimensional synthetic aperture radiometer is available, and the η-dimensional target function becomes as follows, that is,

[0080] The fusion processing η-dimensional target function is

[0081] When the fusion processing η-dimensional target function takes the minimum value, the parameter η corresponding to the minimum value is the optimal estimation value of the position of the RFI source in the η dimension, is the measurement data of the two-dimensional synthetic aperture radiometer on the position of the kth RFI source in the η dimension, is the variance of the measurement data of the two-dimensional synthetic aperture radiometer, and n is the number of data groups.

[0082] The fusion processing energy intensity target function is as follows:

[0083]

[0084] For each radio frequency interference source, the function value of the fusion processing energy intensity objective function is calculated based on its corresponding n sets of first parameters and n sets of second parameters. The parameter T corresponding to the minimum value of the fusion processing energy intensity objective function is the optimal estimate of the energy intensity of the radio frequency interference source. The energy intensity value corresponding to the minimum position of ξ and η is the position of the minimum energy intensity. That is, the energy intensity value corresponding to the position determined by the optimal estimate of the RFI source position in the η dimension and the optimal estimate of the RFI source position in the ξ dimension is the position when the fusion processing energy intensity objective function takes the minimum value. The parameter T corresponding to the minimum value of the fusion processing energy intensity objective function is determined based on the position when the fusion processing energy intensity objective function takes the minimum value. The data represent the energy intensity measurements of radio frequency interference sources taken by one-dimensional and two-dimensional synthetic aperture radiometers, respectively. , respectively, represent the variances of the measurement data of the radio frequency interference source with respect to energy intensity by the one-dimensional and two-dimensional synthetic aperture radiometers, where n is the number of data sets.

[0085] Step S2 involves determining the third parameter corresponding to all radio frequency interference sources, forming fused data of radio frequency interference data from one-dimensional and two-dimensional synthetic aperture radiometers, including:

[0086] Based on the optimal estimates of the location and energy intensity of each radio frequency interference source, a third parameter corresponding to all radio frequency interference sources is determined. This third parameter includes the first optimal value ξ of the fused data. opt The second optimal point η opt And the third best point T opt ;

[0087]

[0088]

[0089] in, For the i-th data set ξ of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, For the η-th data set of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, For the i-th set of energy intensity data of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, Let ξ be the value of the i-th data set of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer. Let η be the value of the i-th data set of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer. This represents the energy intensity data of the i-th group of the k-th radio frequency interference source in the field of view of a two-dimensional synthetic aperture radiometer. variance of the measurement data of the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer, respectively, variance of the measurement data of the two-dimensional synthetic aperture radiometer, variance of the measurement data of the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer on the energy intensity of the radio frequency interference source, respectively.

[0090] combining f(ξ), f(η) and f(T) as a fusion target function of the radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometers, and the optimal solution of the combination of ξ opt , η opt , T opt as the optimal fusion data of the radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometers.

[0091] The step S3: comparing the performance of the measurement data and the fusion data, including calculating one or more of the first root mean square error, the second root mean square error and the third root mean square error, wherein,

[0092] The first root mean square error T_RMSE:

[0093]

[0094] wherein N represents the total number of measurement data, is the energy intensity of the radio frequency interference source detected based on the i-th group of data;

[0095] The second root mean square error P_RMSE,

[0096]

[0097] wherein, is the positioning position of the radio frequency interference source detected based on the i-th group of data, x represents the position of the radio frequency interference source determined based on the optimal fusion data, and takes the optimal estimation value of the position of the RFI source in the ξ dimension or the optimal estimation value of the position of the RFI source in the η dimension;

[0098] The third root mean square error RE,

[0099] RE = T mitigated -T raw

[0100] wherein T mitigated represents the suppressed brightness temperature determined based on the optimal fusion data, and T raw represents the original scene brightness temperature.

[0101] In this embodiment, the first root mean square error of the measured energy intensity before and after data fusion and the true energy intensity is taken as the performance index, the second root mean square error of the positioning position before and after data fusion and the true position, and the third root mean square error represented by the standard deviation of the residual graph between the residual results after calculation and the ideal inversion results before adding the RFI source are taken as the representation of the good and bad of the suppression results to quantitatively describe the good and bad of the fusion results.

[0102] The application provides the fusion method of the RFI data of the one-dimensional and two-dimensional synthetic aperture radiometers.

[0103] (1) The ideal original scene visibility data without RFI source is generated by using the darkroom background measurement data of the two-dimensional synthetic aperture radiometer and the one-dimensional synthetic aperture radiometer in the RFI processing experiment, then the brightness temperature data is obtained according to the relationship between the visibility and the brightness temperature data, an RFI source with a specified energy intensity is added at a specified position (the coordinates are (0.2, 0.2) in the (ξ, η) coordinate system) in the ideal original scene to generate a simulated RFI scene. The inversion images of the simulated RFI scene of the two-dimensional synthetic aperture radiometer and the one-dimensional synthetic aperture radiometer are shown in Figure 3 and Figure 4 .

[0104] (2) 50 groups of different RFI-containing single-snapshot positioning and energy intensity result data of the one-dimensional array / two-dimensional array are obtained. 50 snapshots in the darkroom background measurement data of the one-dimensional array and the two-dimensional array are selected for processing to obtain 50 groups of simulated scenes used as different snapshots (one snapshot is a group of data at a moment).

[0105] (3) 50 groups of snapshot positioning and energy intensity result data of the one-dimensional array / two-dimensional array are obtained. The 50 groups of snapshot suppression results are obtained, and the 50 groups of RFI source positioning and energy intensity results of the Y-shaped array and the one-dimensional array are obtained, the corresponding parameters when the minimization is obtained, and the optimal estimation value of the RFI source position and the energy intensity is obtained.

[0106] Table 1 optimal estimation value

[0107] ξ η RFI_position RFI_T(K) 0.19994 0.19985 0.19994+0.19985i 997.7111

[0108] (4) Root mean square error analysis is performed on the positioning and energy intensity results before and after data fusion

[0109] The root mean square error between the positioning and energy intensity results before and after data fusion and the true data is calculated, and the calculation result is shown in Table 2:

[0110] Table 2 positioning and energy intensity RMSE

[0111] Status err or]]> ​ err or]]> ​ T err or(K)]]> Before fusion 1.4862e-04 2.8995e-04 5.5368 After fusion 5.8399e-05 1.5177e-04 2.2889

[0112] As can be seen from the table data, the positioning accuracy of RFI in the ξ dimension has been improved by an order of magnitude in both the ξ and η dimensions, and the accuracy of the energy intensity measurement has been improved from 5.5368 to 2.2889.

[0113] (5) Residual suppression analysis before and after data fusion

[0114] A new simulation scene is regenerated, and the brightness temperature residual image and its standard deviation between the suppressed visibility and the original visibility are calculated. The Y-shaped array result is as follows. Figure 5 and Figure 6 As shown.

[0115] Figure 5 and Figure 6 The standard deviations are 0.5148 and 0.3661, respectively. It can be seen that the standard deviation of the brightness temperature residual map after using the fusion result is smaller, which proves that data fusion reduces the suppression residual of the RFI source of the two-dimensional synthetic aperture radiometer.

[0116] (6) The result of the one-dimensional matrix is ​​as follows Figure 7 and Figure 8 As shown. Figure 7 and Figure 8 The standard deviations were 1.3008 and 0.5581, respectively. The standard deviation of the brightness temperature residual map after using the fusion result was smaller, which proves that data fusion reduced the suppression residual of the one-dimensional synthetic aperture radiometer RFI source.

[0117] like Figure 9 As shown, the present invention also provides a device for fusing radio frequency interference data from one-dimensional and two-dimensional synthetic aperture radiometers, the device comprising:

[0118] Data acquisition module: The satellite is configured to carry both a one-dimensional synthetic aperture radiometer and a two-dimensional synthetic aperture radiometer to measure the same area; acquire the measurement data from the one-dimensional synthetic aperture radiometer to determine the first parameter of each radio frequency interference (RFI) source; acquire the measurement data from the two-dimensional synthetic aperture radiometer to determine the second parameter of each radio frequency interference (RFI) source;

[0119] Fusion module: configured to construct a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model based on the first parameter and the second parameter, wherein the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model includes multiple objective functions; determine the third parameter corresponding to all radio frequency interference sources, and form fused data of radio frequency interference data of one-dimensional and two-dimensional synthetic aperture radiometers;

[0120] Comparison module: configured to compare the performance of the measurement data with the fused data.

[0121] The specific embodiments described above only illustrate the design principles of the present invention. The shapes and names of the components in this description may differ and are not limited. Therefore, those skilled in the art can modify or make equivalent substitutions to the technical solutions described in the foregoing embodiments; and these modifications and substitutions do not depart from the inventive spirit and technical solutions of the present invention, and should all fall within the protection scope of the present invention.

Claims

1. A method for fusing radio frequency interference data from one-dimensional and two-dimensional synthetic aperture radiometers, characterized in that, include: Step S1: The satellite carries both a one-dimensional synthetic aperture radiometer and a two-dimensional synthetic aperture radiometer to measure the same area. Acquire measurement data from a one-dimensional synthetic aperture radiometer to determine the first parameter of each radio frequency interference (RFI) source; acquire measurement data from a two-dimensional synthetic aperture radiometer to determine the second parameter of each radio frequency interference (RFI) source; Step S2: Based on the first parameter and the second parameter, construct a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model, which includes multiple objective functions; determine the third parameter corresponding to all radio frequency interference sources to form fused data of radio frequency interference data of the one-dimensional and two-dimensional synthetic aperture radiometer; Step S3: Compare the performance of the measured data with the fused data; Step S1 involves acquiring measurement data from a one-dimensional synthetic aperture radiometer and determining the first parameters of each radio frequency interference (RFI) source, including: Step S111: Obtain measurement data from a one-dimensional synthetic aperture radiometer and form n sets of data; Step S112: Determine the first parameters of each radio frequency interference (RFI) source. The first parameters include the first detection position and first energy intensity data of each RFI source. Each RFI source has n sets of first parameters. In the i-th set of first parameters of the k-th RFI source, the first detection position is ( , ), For the i-th data set ξ of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, For the value of the ηth dimension of the i-th data set of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, the first energy intensity data is: , The first set of energy intensity data for the kth radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer is used to characterize the energy intensity of the radio frequency interference source. In step S2, the one-dimensional and two-dimensional integrated aperture radiometer fusion processing model includes a fusion processing position objective function and a fusion processing energy intensity objective function; the fusion processing position objective function includes a fusion processing ξ-dimensional objective function and a fusion processing η-dimensional objective function. The ξ-dimensional objective function for the fusion processing is: The parameter ξ corresponding to the minimum value of the ξ-dimensional objective function of the fusion processing is the optimal estimate of the position of the RFI source in the ξ-dimensional dimension. , ξ represents the variance of the measurement data from the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer, respectively, where n is the number of data sets; The η-dimensional objective function for the fusion processing is: The parameters corresponding to the minimum value of the η-dimensional objective function of the fusion processing That is, the RFI source is The optimal estimate of the position in the dimension. For two-dimensional synthetic aperture radiometers to measure RFI sources in Measurement data of the location of the k-th RFI source in dimension 1. The variance of the η-dimensional measurement data of the two-dimensional synthetic aperture radiometer; The objective function for the energy intensity of the fusion process is as follows: Specifically, for each radio frequency interference source, the function value of the fusion processing energy intensity objective function is calculated based on its corresponding n sets of first parameters and n sets of second parameters. The parameter corresponding to the minimum value of the fusion processing energy intensity objective function is... T This is the optimal estimate of the energy intensity of the radio frequency interference source; This provides measurement data on the energy intensity of radio frequency interference sources using a two-dimensional synthetic aperture radiometer. , The variances of the measurement data of radio frequency interference sources with respect to energy intensity by the one-dimensional and two-dimensional synthetic aperture radiometers are respectively. The determination of the third parameter corresponding to all radio frequency interference sources, forming fused data of radio frequency interference data from one-dimensional and two-dimensional synthetic aperture radiometers, includes: Based on the optimal estimates of the location and energy intensity of each radio frequency interference source, a third parameter corresponding to all radio frequency interference sources is determined. This third parameter includes the first optimal value of the fused data. ξ opt Second best advantage η opt And the third best advantage T opt ; in, To the first in the field of view of a one-dimensional synthetic aperture radiometer k The first radio frequency interference source i The values ​​of the ξ dimension of the data set. Let ξ be the value of the i-th data set of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer. To the first in the field of view of a two-dimensional synthetic aperture radiometer k The first radio frequency interference source i The values ​​of the n-dimensional set of data, To the first in the field of view of a two-dimensional synthetic aperture radiometer k The first radio frequency interference source i Group energy intensity data, , The variances of the ξ-dimensional measurement data from the one-dimensional and two-dimensional synthetic aperture radiometers are respectively. Two-dimensional synthetic aperture radiometer η Variance of the measurement data , These represent the variances of the energy intensity measurement data of the radio frequency interference source by the one-dimensional and two-dimensional synthetic aperture radiometers, respectively.

2. The method as described in claim 1, characterized in that, Step S1 involves acquiring measurement data from a two-dimensional synthetic aperture radiometer and determining the second parameters of each radio frequency interference (RFI) source, including: Step S121: Obtain measurement data from the two-dimensional integrated aperture radiometer and form n sets of data; Step S122: Determine the second parameters of each radio frequency interference (RFI) source. The second parameters include the second detection position and second energy intensity data of each RFI source. Each RFI source has n sets of second parameters. In the i-th set of second parameters for the k-th RFI source, the second detection position is ( , ), Let ξ be the value of the i-th data set of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer. Let η be the value of the i-th data set of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer. The second energy intensity data is the i-th set of energy intensity data for the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer, and is used to characterize the energy intensity of the radio frequency interference source.

3. The method as described in claim 2, characterized in that, Step S3: Comparing the performance of the measured data and the fused data includes calculating one or more of a first root mean square error, a second root mean square error, and a third root mean square error, wherein... First root mean square error ： in, N Represents the total number of measurement data. For the first i The energy intensity of radio frequency interference sources detected by the group data; Second root mean square error , in, For the first i The location of the radio frequency interference source detected by the group data. This indicates the location of the radio frequency interference source determined based on the optimal fusion data, and the value is taken as the location of the RFI source in... ξ The optimal estimate of the position in the dimension or the RFI source in The optimal estimate of the position in the dimension; The third root mean square error , RE=T mitigated - T raw in, T mitigated This indicates the brightness temperature after suppression, determined based on the optimal fusion data. T raw This indicates the original scene brightness temperature.

4. A radio frequency interference data fusion device for one-dimensional and two-dimensional synthetic aperture radiometers, characterized in that, include: Data acquisition module: configured to carry both a one-dimensional and a two-dimensional synthetic aperture radiometer on the satellite to measure the same area; Acquire measurement data from a one-dimensional synthetic aperture radiometer to determine the first parameter of each radio frequency interference (RFI) source; acquire measurement data from a two-dimensional synthetic aperture radiometer to determine the second parameter of each radio frequency interference (RFI) source; Fusion module: configured to construct a one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model based on the first parameter and the second parameter, wherein the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model includes multiple objective functions; determine the third parameter corresponding to all radio frequency interference sources, and form fused data of radio frequency interference data of one-dimensional and two-dimensional synthetic aperture radiometers; Comparison module: configured to compare the performance of the measurement data with the fused data; The data acquisition module includes: The first submodule is configured to acquire measurement data from a one-dimensional synthetic aperture radiometer, forming n sets of data. The second submodule is configured to determine the first parameters of each radio frequency interference (RFI) source. These first parameters include the first detection location and first energy intensity data for each RFI source. Each RFI source has n sets of first parameters. In the i-th set of first parameters for the k-th RFI source, the first detection location is (…). , ), For the i-th data set ξ of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, For the value of the ηth dimension of the i-th data set of the k-th radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer, the first energy intensity data is: , The first set of energy intensity data for the kth radio frequency interference source in the field of view of a one-dimensional synthetic aperture radiometer is used to characterize the energy intensity of the radio frequency interference source. The fusion module, wherein the one-dimensional and two-dimensional synthetic aperture radiometer fusion processing model includes a fusion processing position objective function and a fusion processing energy intensity objective function; the fusion processing position objective function includes a fusion processing ξ-dimensional objective function and a fusion processing η-dimensional objective function; The ξ-dimensional objective function for the fusion processing is: The parameter ξ corresponding to the minimum value of the ξ-dimensional objective function of the fusion processing is the optimal estimate of the position of the RFI source in the ξ-dimensional dimension. , ξ represents the variance of the measurement data from the one-dimensional synthetic aperture radiometer and the two-dimensional synthetic aperture radiometer, respectively, where n is the number of data sets; The η-dimensional objective function for the fusion processing is: The parameters corresponding to the minimum value of the η-dimensional objective function of the fusion processing That is, the RFI source is The optimal estimate of the position in the dimension. For two-dimensional synthetic aperture radiometers to measure RFI sources in Measurement data of the location of the k-th RFI source in dimension 1. The variance of the η-dimensional measurement data of the two-dimensional synthetic aperture radiometer; The objective function for the energy intensity of the fusion process is as follows: Specifically, for each radio frequency interference source, the function value of the fusion processing energy intensity objective function is calculated based on its corresponding n sets of first parameters and n sets of second parameters. The parameter corresponding to the minimum value of the fusion processing energy intensity objective function is... T This is the optimal estimate of the energy intensity of the radio frequency interference source; This provides measurement data on the energy intensity of radio frequency interference sources using a two-dimensional synthetic aperture radiometer. , The variances of the measurement data of radio frequency interference sources with respect to energy intensity by the one-dimensional and two-dimensional synthetic aperture radiometers are respectively. The determination of the third parameter corresponding to all radio frequency interference sources, forming fused data of radio frequency interference data from one-dimensional and two-dimensional synthetic aperture radiometers, includes: Based on the optimal estimates of the location and energy intensity of each radio frequency interference source, a third parameter corresponding to all radio frequency interference sources is determined. This third parameter includes the first optimal value of the fused data. ξ opt Second best advantage η opt And the third best advantage T opt ; in, To the first in the field of view of a one-dimensional synthetic aperture radiometer k The first radio frequency interference source i The values ​​of the ξ dimension of the data set. Let ξ be the value of the i-th data set of the k-th radio frequency interference source in the field of view of the two-dimensional synthetic aperture radiometer. To the first in the field of view of a two-dimensional synthetic aperture radiometer k The first radio frequency interference source i The values ​​of the n-dimensional set of data, To the first in the field of view of a two-dimensional synthetic aperture radiometer k The first radio frequency interference source i Group energy intensity data, , The variances of the ξ-dimensional measurement data from the one-dimensional and two-dimensional synthetic aperture radiometers are respectively. Two-dimensional synthetic aperture radiometer η Variance of the measurement data , These represent the variances of the energy intensity measurement data of the radio frequency interference source by the one-dimensional and two-dimensional synthetic aperture radiometers, respectively.

5. A computer-readable storage medium storing a plurality of instructions; the plurality of instructions being loaded by a processor and executing the method as claimed in any one of claims 1-3.

6. An electronic device, characterized in that, The electronic device includes: A processor is used to execute multiple instructions; Memory, used to store multiple instructions; The plurality of instructions are to be stored in the memory and loaded by the processor and executed as described in any one of claims 1-3.

Images