Method and device for processing multi-path radiation reference star - to - star transfer results

Abstract

The present invention provides a method and device for processing the results of multi-path radiative reference inter-satellite transfer, which are applied to the technical field of space radiation measurement. The method includes: obtaining and calculating the atmospheric top spectral radiance correction amount and the uncertainty of a target spaceborne instrument through the multi-path radiative reference inter-satellite transfer method, where one reference satellite corresponds to one atmospheric top spectral radiance correction amount and one uncertainty; removing the outliers in the atmospheric top spectral radiance correction amounts to obtain the atmospheric top spectral radiance correction amounts and the corresponding uncertainties; calculating the target transfer result based on the weighted average method when the uncertainties do not meet the first condition; and calculating the target transfer result based on the arithmetic average method when the uncertainties meet the first condition.

Description

Technical Field

[0001] The present invention relates to the technical field of space radiation measurement, and particularly to a method and device for processing the results of multi-path radiation reference inter-satellite transfer. Background Art

[0002] The space radiation measurement reference satellite is the key to realizing the unification of the radiation measurement reference for remote sensing constellations and overall improving the high-precision quantitative remote sensing ability, and is also the high ground of quantitative remote sensing technology that major space remote sensing countries compete for. Using the reference satellite for radiation quantity transfer is currently the most accurate on-orbit calibration method.

[0003] In the prior art, the radiation reference inter-satellite transfer method usually selects the spectral radiance data of the top of the atmosphere measured by the reference star sensor with perfect calibration equipment and relatively high calibration accuracy as a reference, and performs spatial matching, spectral matching, etc. on the near-orbit overlapping data pairs of the reference instrument and the instrument to be calibrated that are relatively close in time and space, so as to obtain the spectral radiance correction amount or calibration coefficient of the top of the atmosphere of the instrument to be calibrated, thereby completing the on-orbit radiation calibration of the target instrument and the evaluation of the original calibration deviation. Multi-path radiation reference transfer refers to using different reference stars and different transfer methods to perform inter-satellite radiation reference quantity transfer on the target star, and obtaining a series of spectral radiance correction amounts of the top of the atmosphere of the target satellite-borne instrument and their uncertainties.

[0004] Due to the different reference stars and transfer methods used, the correlation between the series of transferred spectral radiance correction amounts of the top of the atmosphere is weak. How to scientifically and reasonably analyze and process the series of spectral radiance correction amounts of the top of the atmosphere and their uncertainties obtained from multi-path radiation reference inter-satellite transfer to obtain the final spectral radiance correction amount and uncertainty of the top of the atmosphere of the target satellite-borne instrument is a common problem encountered by series satellites such as meteorological, oceanographic, and high-resolution land satellites when using radiation reference quantity transfer. Summary of the Invention

[0005] The present invention provides a method and device for processing the results of multi-path radiation reference inter-satellite transfer, so as to solve the problem of how to scientifically and reasonably analyze and process the series of spectral radiance correction amounts of the top of the atmosphere and their uncertainties obtained from multi-path radiation reference inter-satellite transfer.

[0006] The present invention provides a method for processing the results of multi-path radiation reference inter-satellite transfer, including: through the multi-path radiation reference inter-satellite transfer method, obtaining and calculating spectral radiance correction amounts of the top of the atmosphere of the target satellite-borne instrument and uncertainties, one reference star corresponding to one spectral radiance correction amount of the top of the atmosphere and one uncertainty; removing the outliers in the spectral radiance correction amounts of the top of the atmosphere to obtain the spectral radiance correction amount at the top of the atmosphere and the corresponding uncertainties; when the uncertainties do not meet the first condition, calculating the target transfer result based on the weighted average method; when the uncertainties meet the first condition, calculating the target transfer result based on the arithmetic average method; wherein, the first condition is that the gap between any uncertainty and the average uncertainty of the uncertainties is less than the first threshold, greater than or equal to and and are both positive integers.

[0007] According to a method for processing the interstellar transfer result of a multi-path radiation reference star provided by the present invention, removing the outliers from the spectral radiance correction amounts at the top of the atmosphere to obtain spectral radiance correction amounts at the top of the atmosphere and the corresponding uncertainties, including: removing outliers from the judgment criteria from the spectral radiance correction amounts at the top of the atmosphere to obtain spectral radiance correction amounts at the top of the atmosphere and the corresponding uncertainties.

[0008] According to a method for processing the interstellar transfer result of a multi-path radiation reference star provided by the present invention, when the uncertainties do not meet the first condition, calculating the target transfer result based on the weighted average method, including: determining the corresponding cut-off uncertainty by comparing the size of each uncertainty with the median of the uncertainties, and adjusting the original uncertainty based on the cut-off uncertainty to obtain the adjusted target uncertainty; determining the weight coefficient of each spectral radiance correction amount at the top of the atmosphere according to the adjusted target uncertainty, and the weight coefficient is inversely proportional to the square of the target uncertainty; determining the weighted average spectral radiance correction amount at the top of the atmosphere based on the weight coefficient and the spectral radiance correction amount at the top of the atmosphere; determining the uncertainty of the weighted average spectral radiance correction amount at the top of the atmosphere according to the uncertainty of each spectral radiance correction amount at the top of the atmosphere; determining the weighted average spectral radiance correction amount at the top of the atmosphere and the uncertainty of the weighted average spectral radiance correction amount at the top of the atmosphere as the target transfer result.

[0009] According to a method for processing the interstellar transfer result of a multi-path radiation reference star provided by the present invention, when the When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method, including: calculating the The average top-of-atmosphere spectral radiance correction of the top-of-atmosphere spectral radiance correction and the The standard uncertainty of an uncertainty; the average top-of-atmosphere spectral radiance correction amount and the standard uncertainty are determined as the target transfer result.

[0010] The present invention also provides a processing device for multi-path radiation benchmark inter-satellite transfer results, comprising the following modules: an acquisition module and a processing module; the acquisition module is used to acquire and calculate the target satellite-borne instrument through the multi-path radiation benchmark inter-satellite transfer method; The top-of-atmosphere spectral radiance correction and uncertainty, one reference star corresponds to one atmospheric top spectral radiation brightness correction and one uncertainty; the processing module is used to eliminate the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty; in the When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method; When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method; wherein the first condition is that any uncertainty and the The difference between the average uncertainties of the uncertainties is smaller than the first threshold. Greater than or equal to ,and , All are positive integers.

[0011] According to a processing device for multi-path radiation benchmark inter-satellite transfer results provided by the present invention, the processing module is used to The judgment criteria are as follows Remove outliers from the atmospheric top spectral radiance correction and obtain The top-of-atmosphere spectral radiance correction and the corresponding An uncertainty.

[0012] According to a processing device for multi-path radiation benchmark inter-satellite transfer results provided by the present invention, the processing module is used to compare each uncertainty with the Determine the corresponding cut-off uncertainty based on the magnitude of the median of the uncertainties, and adjust the original uncertainty based on the cut-off uncertainty to obtain the adjusted target uncertainty; determine the weight coefficient of each top-of-atmosphere spectral radiance correction amount according to the adjusted target uncertainty, and the weight coefficient is inversely proportional to the square of the target uncertainty; determine the weighted average top-of-atmosphere spectral radiance correction amount based on the weight coefficient and the top-of-atmosphere spectral radiance correction amount; determine the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount according to the uncertainty of each top-of-atmosphere spectral radiance correction amount; determine the weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount as the target transfer result.

[0013] According to a processing device for multi-path radiation reference inter-satellite transfer results provided by the present invention, the processing module is used to calculate the average top-of-atmosphere spectral radiance correction amount of the top-of-atmosphere spectral radiance correction amounts and the standard uncertainty of the uncertainties; determine the average top-of-atmosphere spectral radiance correction amount and the standard uncertainty as the target transfer result.

[0014] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the processing method of the multi-path radiation reference inter-satellite transfer result as described in any one of the above.

[0015] The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it implements the processing method of the multi-path radiation reference inter-satellite transfer result as described in any one of the above.

[0016] The present invention also provides a computer program product, including a computer program. When the computer program is executed by a processor, it implements the processing method of the multi-path radiation reference inter-satellite transfer result as described in any one of the above.

[0017] The processing method and device for multi-path radiation reference inter-satellite transfer results provided by the present invention can first eliminate the outliers in the top-of-atmosphere spectral radiance correction amounts to obtain the top-of-atmosphere spectral radiance correction amounts and the corresponding uncertainties, and then according to The way to calculate the target transfer result is selected according to the difference degree of uncertainties. In this way, a scientific and reasonable analysis and processing can be carried out on the series of radiation exoatmospheric spectral radiance correction amounts and their uncertainties obtained from the multi-path radiometric reference star inter-satellite transfer based on a unified method, so as to obtain the final radiation exoatmospheric spectral radiance correction amount and uncertainty of the target spaceborne instrument, and further improve the scientific rationality of the multi-path transfer fusion result. Brief Description of the Drawings

[0018] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0019] Figure 1 It is a schematic flowchart of the method for processing the multi-path radiometric reference star inter-satellite transfer result provided by the present invention;

[0020] Figure 2 It is a schematic structural diagram of the device for processing the multi-path radiometric reference star inter-satellite transfer result provided by the present invention;

[0021] Figure 3 It is a schematic structural diagram of the electronic device provided by the present invention. Detailed Embodiments

[0022] To make the objectives, technical solutions, and advantages of the present application clearer, the following will clearly and completely describe the technical solutions in the present application with reference to the drawings in the present application. Obviously, the described embodiments are some, but not all, of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present application.

[0023] It should be noted that in the embodiments of the present application, words such as "exemplary" or "for example" are used to indicate examples, illustrations, or explanations. Any embodiment or design solution described as "exemplary" or "for example" in the embodiments of the present application should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Exactly speaking, using words such as "exemplary" or "for example" is intended to present relevant concepts in a specific manner.

[0024] It should be noted that in this text, the term "including", "comprising", or any other variant thereof is intended to cover non-exclusive inclusion, such that a process, method, article, or device comprising a series of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or device. Without further limitation, an element defined by the phrase "including a..." does not exclude the presence of additional identical elements in the process, method, article, or device that includes such element. In addition, it should be pointed out that the scope of the methods and devices in the embodiments of the present application is not limited to performing functions in the order shown or discussed, and may also include performing functions in a substantially simultaneous manner or in a reverse order according to the functions involved. For example, the described method may be performed in an order different from that described, and various steps may be added, omitted, or combined. Additionally, features described with reference to certain examples may be combined in other examples.

[0025] For the convenience of clearly describing the technical solutions of the embodiments of the present application, in the embodiments of the present application, terms such as "first" and "second" are used to distinguish identical or similar items with substantially the same functions and roles. Those skilled in the art can understand that terms such as "first" and "second" do not limit the quantity and execution order.

[0026] Some exemplary embodiments are described for the purpose of illustration in the embodiments of the present application. It should be understood that the present application can be implemented in other ways not specifically shown in the drawings.

[0027] As Figure 1 shown, the embodiments of the present application provide a method for processing the results of multi-path radiative reference inter-satellite transfer. The method for processing the results of multi-path radiative reference inter-satellite transfer can be applied to a device for processing the results of multi-path radiative reference inter-satellite transfer. The method for processing the results of multi-path radiative reference inter-satellite transfer may include S101 - S104:

[0028] S101. The device for processing the results of multi-path radiative reference inter-satellite transfer obtains and calculates atmospheric top spectral radiance correction amounts and uncertainties of the target spaceborne instrument through the multi-path radiative reference inter-satellite transfer method.

[0029] Wherein, one reference satellite corresponds to one atmospheric top spectral radiance correction amount and one uncertainty.

[0030] In the field of remote sensing, in order to accurately measure the radiation value of the target star, it is necessary to use the reference star for radiation reference transfer. For the reference star, its radiation value is a standard value that has been accurately measured. Due to factors such as different reference stars and transfer methods, there will be radiation reference transfer results in multiple ways. For example, when observing meteorological, ocean, and land targets by satellite remote sensing, radiation reference transfer may be carried out through instruments carried by satellites at different orbital positions, using different optical band observation methods, or different calibration models.

[0031] When measuring the radiation value with an uncalibrated target star-borne instrument, there will be errors. By transferring the accurate radiation value of the reference star to the target star, an atmospheric top spectral radiance correction amount for correcting the measured value of the target star can be obtained. Each atmospheric top spectral radiance correction amount corresponds to an uncertainty. These uncertainties come from multiple aspects, including the measurement uncertainty of the radiation value of the reference star itself, instrument errors due to differences in time domain, spatial domain, spectral domain, and angular domain during the transfer process, environmental factor impacts, and model uncertainties. The uncertainty reflects the credibility of the atmospheric top spectral radiance correction amount. A smaller uncertainty means that the atmospheric top spectral radiance correction amount is more reliable.

[0032] S102. The processing device for the transfer results between multi-path radiation reference stars eliminates the outliers among the atmospheric top spectral radiance correction amounts to obtain atmospheric top spectral radiance correction amounts and the corresponding uncertainties.

[0033] Optionally, the processing device for the transfer results between multi-path radiation reference stars can eliminate outliers from the judgment criterion from the atmospheric top spectral radiance correction amounts to obtain atmospheric top spectral radiance correction amounts and the corresponding uncertainties.

[0034] Specifically, due to some sudden and non-repeatable factors, there may be outlier abnormal values among the

[0035] atmospheric top spectral radiance correction amounts. For example, when conducting remote sensing observations to obtain the atmospheric top spectral radiance correction amount for radiation reference transfer, it may be affected by factors such as cosmic ray bursts, sudden loss of satellite attitude control, and temporary electromagnetic interference. These factors occur accidentally and are difficult to repeat under normal experimental conditions. They will cause some atmospheric top spectral radiance correction amounts to deviate from the normal range, so they need to be eliminated. The specific steps are as follows:

[0035] Step 1. The processing device for the transfer results between multi-path radiation reference stars can first calculate Average value of the spectral radiance correction amount at the top of the atmosphere 。

[0036] Calculate the average value of the spectral radiance correction amount at the top of the atmosphere, and the formula is: 。

[0037] Step 2: The processing device for the multi-path radiation inter-satellite transfer result calculates the residual of each spectral radiance correction amount at the top of the atmosphere based on this average value 。

[0038] The residual represents the difference between each spectral radiance correction amount at the top of the atmosphere and the average value. Calculate the residual of each spectral radiance correction amount at the top of the atmosphere, and the formula is: 。

[0039] Step 3: The processing device for the multi-path radiation inter-satellite transfer result calculates the experimental standard deviation of each spectral radiance correction amount at the top of the atmosphere based on this average value 。

[0040] The experimental standard deviation measures the degree of dispersion of the spectral radiance correction amount at the top of the atmosphere. It comprehensively considers the difference between each spectral radiance correction amount at the top of the atmosphere and the average value. The larger the standard deviation, the greater the degree of dispersion between the spectral radiance correction amounts at the top of the atmosphere. Calculate the experimental standard deviation of each spectral radiance correction amount at the top of the atmosphere, and the formula is: 。

[0041] Step 4: The processing device for the multi-path radiation inter-satellite transfer result eliminates the outliers based on the residuals and the experimental standard deviation

[0042] In the case of a normal distribution, as long as the number of measurements is not too small, the probability that the absolute value of the residual exceeds three times the experimental standard deviation is very small, so it can be considered an impossible event. If is the one with the largest absolute value among the residuals, when it satisfies the condition , then it is considered that this spectral radiance correction amount at the top of the atmosphere belongs to an outlier and should be eliminated. This is based on the " " principle of the normal distribution. For example, if the absolute value of the residual of a spectral radiance correction amount at the top of the atmosphere is much larger than three times the experimental standard deviation, then this spectral radiance correction amount at the top of the atmosphere is very likely to be affected by the sudden abnormal factors mentioned above, so it is eliminated as an outlier to ensure the accuracy of subsequent data processing and analysis.

[0043] After outliers are removed, the above procedure should be repeatedly used until no outliers are included in the spectral radiance correction at the top of the atmosphere. This is because after removing one outlier, the statistics such as the mean, residual, and standard deviation of the remaining spectral radiance corrections at the top of the atmosphere will change, which may cause the spectral radiance corrections at the top of the atmosphere that were not originally judged as outliers to become outliers in the new statistical environment. Therefore, this process needs to be repeated to ensure that the finally remaining spectral radiance corrections at the top of the atmosphere are all relatively reasonable and reliable data, so as to provide an accurate data basis for subsequent operations such as radiometric reference transfer and data fusion.

[0044] S103. When the uncertainties of the multi-path radiometric reference inter-satellite transfer results do not meet the first condition, calculate the target transfer result based on the weighted average method. Among them, the first condition is that the gap between any uncertainty and the average uncertainty of the

[0045] uncertainties is less than the first threshold, greater than or equal to and , and and are both positive integers.

[0046] Optionally, the processing device for the multi-path radiometric reference inter-satellite transfer results calculates the target transfer result based on the weighted average method, including: determining the corresponding cut-off uncertainty by comparing each uncertainty with the median of the uncertainties, and adjusting the original uncertainty based on the cut-off uncertainty to obtain the adjusted target uncertainty; determining the weight coefficient of each spectral radiance correction at the top of the atmosphere according to the adjusted target uncertainty, and the weight coefficient is inversely proportional to the square of the adjusted target uncertainty; determining the weighted average spectral radiance correction at the top of the atmosphere based on the weight coefficient and the spectral radiance correction at the top of the atmosphere; determining the uncertainty of the weighted average according to the uncertainty of each spectral radiance correction at the top of the atmosphere; and determining the weighted average spectral radiance correction at the top of the atmosphere and the uncertainty of the weighted average spectral radiance correction at the top of the atmosphere as the target transfer result.

[0047] Specifically, after multi-path radiometric reference transfer, for the spectral radiance correction at the top of the atmosphere (TOA) of the target star at the wavelength position, the corresponding standard uncertainties have significant differences. This difference makes it impossible to simply perform arithmetic averaging on the spectral radiance corrections at the top of the atmosphere to fuse the results.

[0048] The weighted average method is considered the optimal linear estimation method for obtaining the fused average value. It can assign different weights according to the reliability of each top-of-atmosphere spectral radiance correction amount (reflected by the uncertainty), so as to more reasonably integrate the information of each top-of-atmosphere spectral radiance correction amount.

[0049] To avoid the inappropriate influence of some too low and unreasonable uncertainties on the final fusion result, the concept of cut-off uncertainty is proposed in this application.

[0050] When When . That is, when the uncertainty is less than or equal to the median of the uncertainties , its cut-off uncertainty is set to the average value of all uncertainties. This can prevent the top-of-atmosphere spectral radiance correction amount with too small uncertainty from occupying an unreasonable weight in the fusion process.

[0051] Adjust the original uncertainty according to the cut-off uncertainty to obtain the adjusted target uncertainty. If the uncertainty , then the target uncertainty ; if the uncertainty , then the target uncertainty .

[0052] Calculate the weight coefficient of each top-of-atmosphere spectral radiance correction amount according to the target uncertainty. The weight coefficient is inversely proportional to the square of the adjusted target uncertainty, which means that the smaller the uncertainty of the top-of-atmosphere spectral radiance correction amount, the larger its weight coefficient and the more important it is in the fusion process.

[0053] The weighted average top-of-atmosphere spectral radiance correction amount of the target TOA radiance can be obtained by multiplying each top-of-atmosphere spectral radiance correction amount by its weight coefficient and summing, which integrates the information of each top-of-atmosphere spectral radiance correction amount.

[0054] The target uncertainty corresponding to the weighted average top-of-atmosphere spectral radiance correction amount of the target TOA radiance. This uncertainty reflects the reliability of the weighted average correction value.

[0055] Express the difference between each target top-of-atmosphere spectral radiance correction amount and the weighted average top-of-atmosphere spectral radiance correction amount as the equivalence degree . The smaller the equivalence degree, the closer the correction amount of the spectral radiance at the top of the target atmosphere is to the weighted average correction value, and the more accurate the result is.

[0056] When the th correction amount of the spectral radiance at the top of the atmosphere participates in the calculation of the weighted average correction value, the uncertainty of ; when the th correction amount of the spectral radiance at the top of the atmosphere does not participate in the calculation of the weighted average correction value, . These uncertainties can help evaluate the reliability of the difference between the correction amount of the spectral radiance at the top of the atmosphere and the weighted average correction value.

[0057] The weighted average method needs to use the chi-square statistic to conduct a compatibility test on the correction amount of the spectral radiance at the top of the atmosphere. This is to ensure that the correction amounts of the spectral radiance at the top of the atmosphere are compatible with each other during the fusion process, that is, the differences between them are within a reasonable range, ensuring the rationality and reliability of the fusion result. If the correction amounts of the spectral radiance at the top of the atmosphere are incompatible, it may lead to deviations in the fusion result, so the compatibility test is a crucial step.

[0058] S104. When the uncertainties of the multi-path radiation reference inter-star transfer results meet the first condition, the processing device calculates the target transfer result based on the arithmetic mean method. Specifically, when the uncertainties of the correction amounts of the spectral radiance at the top of the atmosphere of the target star obtained through multi-path reference transfer at wavelength

[0059] are similar, the arithmetic mean method can be considered for fusion. The definition of "similar uncertainties" here is that the maximum deviation between the uncertainty and the average uncertainty does not exceed the first threshold, and the first threshold can be the average correction amount of the spectral radiance at the top of the atmosphere of the correction amounts of the spectral radiance at the top of the atmosphere and the combined standard uncertainty of the

[0060] ; determine the average correction amount of the spectral radiance at the top of the atmosphere and the standard uncertainty as the target transfer result.

[0061] . .

[0061] For example, in the process of performing remote sensing observations to obtain the atmospheric top spectral radiance correction amount of the target star's radiance, if the uncertainties corresponding to the atmospheric top spectral radiance correction amounts obtained through different channels are not significantly different, this may mean that the reliabilities and accuracies of each transmission channel are roughly the same, and the uncertainty of no single channel is significantly prominent. In this case, each atmospheric top spectral radiance correction amount can be regarded as equally important during the fusion process.

[0062] Calculate using the arithmetic mean method The average transmitted atmospheric top spectral radiance correction amount of multiple-channel radiation reference transfer , and the formula is . The meaning of this formula is to add up all the atmospheric top spectral radiance correction amounts and then divide by the number of atmospheric top spectral radiance correction amounts. What is obtained is the arithmetic mean of these atmospheric top spectral radiance correction amounts.

[0063] Calculate The standard uncertainty of the average correction value of multiple-channel radiation reference transfer , and the formula is . This formula takes into account the uncertainty of each atmospheric top spectral radiance correction amount. By taking the square root of the sum of the squares of all uncertainties and then dividing by the number of atmospheric top spectral radiance correction amounts, the standard uncertainty of the average correction value is obtained. It reflects the reliability degree of the atmospheric top spectral radiance correction amount after arithmetic averaging.

[0064] If the uncertainties of all atmospheric top spectral radiance correction amounts are the same, all being , then the formula for the standard uncertainty can be simplified to .

[0065] Like the weighted average method, the arithmetic average method also needs to use the chi-square statistic to perform a compatibility test on the atmospheric top spectral radiance correction amounts. This is because even if the uncertainties are similar, there may be some potential inconsistencies among the atmospheric top spectral radiance correction amounts. Through the compatibility test, it can be ensured that these atmospheric top spectral radiance correction amounts are compatible during the fusion process, that is, the differences between them are within a reasonable range, thereby ensuring the rationality and reliability of the arithmetic average fusion result. If the atmospheric top spectral radiance correction amounts are incompatible, the result of the arithmetic average may be affected, resulting in inaccurate average transmitted atmospheric top spectral radiance correction amount and the corresponding uncertainty.

[0066] In the embodiments of this application, the outliers in atmospheric top spectral radiance correction amounts can be removed first to obtain The spectral radiance correction amount at the top of the atmosphere and the corresponding uncertainties, and then select the method for calculating the target transfer result according to the difference degree of the uncertainties. In this way, a series of spectral radiance correction amounts at the top of the radiation atmosphere obtained by multi-path radiation reference star transfer and their respective uncertainties can be scientifically and reasonably analyzed and processed based on a unified method, so as to obtain the final spectral radiance correction amount and uncertainty at the top of the radiation atmosphere of the target spaceborne instrument, and further improve the scientific rationality of the multi-path transfer fusion result.

[0067] The above mainly introduces the solution provided by the embodiments of the present application from the perspective of the method. To implement the above functions, it includes the corresponding hardware structure and / or software module for executing each function. Those skilled in the art should easily realize that, combining the units and algorithm steps of each example described in the embodiments disclosed herein, the embodiments of the present application can be implemented in the form of hardware or a combination of hardware and computer software. Whether a certain function is executed in the form of hardware or computer software driving hardware depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present application.

[0068] It should be noted that the device in the embodiments of the present application includes a virtual device and a physical device. The virtual device can be a processing device for the transfer result between multi-path radiation reference stars, and the physical device can include an electronic device, a computer storage medium, and a computer program product.

[0069] The method for processing the transfer result between multi-path radiation reference stars provided by the embodiments of the present application may have an execution subject as a processing device for the transfer result between multi-path radiation reference stars, or a control module for processing the transfer result between multi-path radiation reference stars in the processing device for the transfer result between multi-path radiation reference stars. In the embodiments of the present application, taking the processing device for the transfer result between multi-path radiation reference stars to execute the method for processing the transfer result between multi-path radiation reference stars as an example, the processing device for the transfer result between multi-path radiation reference stars provided by the embodiments of the present application is described.

[0070] It should be noted that the embodiments of the present application can divide the functional modules of the processing device for the transfer result between multi-path radiation reference stars according to the above method examples. For example, each functional module can be divided corresponding to each function, or two or more functions can be integrated into one processing module. The above integrated module can be implemented in the form of hardware or in the form of a software functional module. Optionally, the division of modules in the embodiments of the present application is illustrative, only a logical function division, and there may be other division methods in actual implementation.

[0071] like Figure 2 As shown, the embodiment of the present application provides a processing device 200 for multi-path radiation benchmark inter-satellite transfer results. The processing device 200 for multi-path radiation benchmark inter-satellite transfer results includes: an acquisition module 201 and a processing module 202. The acquisition module 201 can be used to acquire and calculate the target satellite instrument through the multi-path radiation benchmark inter-satellite transfer method. The top-of-atmosphere spectral radiance correction and uncertainty, one reference star corresponds to one atmospheric top spectral radiance correction and one uncertainty; the processing module 202 is used to eliminate the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty; in the When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method; When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method; wherein the first condition is that any uncertainty and the The difference between the average uncertainties of the uncertainties is smaller than the first threshold. Greater than or equal to ,and , All are positive integers.

[0072] Optionally, the processing module 202 is used based on The judgment criteria are as follows Remove outliers from the atmospheric top spectral radiance correction and obtain The top-of-atmosphere spectral radiance correction and the corresponding An uncertainty.

[0073] Optionally, the processing module 202 is used to compare each uncertainty with the Determine the corresponding cut-off uncertainty according to the magnitude of the median of the uncertainties, and adjust the original uncertainty based on the cut-off uncertainty to obtain the adjusted target uncertainty; determine the weight coefficient of each top-of-atmosphere spectral radiance correction amount according to the adjusted target uncertainty, and the weight coefficient is inversely proportional to the square of the target uncertainty; determine the weighted average top-of-atmosphere spectral radiance correction amount based on the weight coefficient and the top-of-atmosphere spectral radiance correction amount; determine the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount according to the uncertainty of each top-of-atmosphere spectral radiance correction amount; determine the weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount as the target transfer result.

[0074] Optionally, the processing module 202 is configured to calculate the average top-of-atmosphere spectral radiance correction amount of the top-of-atmosphere spectral radiance correction amounts and the standard uncertainty of the uncertainties; determine the average top-of-atmosphere spectral radiance correction amount and the standard uncertainty as the target transfer result.

[0075] In the embodiments of the present application, outliers in the top-of-atmosphere spectral radiance correction amounts can be removed first to obtain top-of-atmosphere spectral radiance correction amounts and corresponding uncertainties, and then a method for calculating the target transfer result is selected according to the difference degree of the uncertainties. In this way, a series of radiance top-of-atmosphere spectral radiance correction amounts and their respective uncertainties obtained by multi-path radiation reference star transfer can be scientifically and reasonably analyzed and processed based on a unified method, so as to obtain the final radiance top-of-atmosphere spectral radiance correction amount and uncertainty of the target spaceborne instrument, thereby improving the scientific rationality of the multi-path transfer fusion result.

[0076] Figure 3 An example of a schematic physical structure diagram of an electronic device is shown in Figure 3 as shown. The electronic device may include: a processor 310, a communication interface 320, a memory 330, and a communication bus 340. Among them, the processor 310, the communication interface 320, and the memory 330 communicate with each other through the communication bus 340. The processor 310 can call the logical instructions in the memory 330 to execute the processing method for the transfer result between multi-path radiation reference stars, and the method includes: through the multi-path radiation reference star transfer method, obtain and calculate the top-of-atmosphere spectral radiance correction amounts of the target spaceborne instrument and An uncertainty, one reference star corresponding to one correction amount of spectral radiance at the top of the atmosphere and one uncertainty; excluding the outliers in the correction amounts of spectral radiance at the top of the atmosphere to obtain the correction amounts of spectral radiance at the top of the atmosphere and the corresponding uncertainties; when the uncertainties do not meet the first condition, calculating the target transfer result based on the weighted average method; when the uncertainties meet the first condition, calculating the target transfer result based on the arithmetic average method; wherein, the first condition is that the gap between any one uncertainty and the average uncertainty of the uncertainties is less than the first threshold, greater than or equal to , and are both positive integers.

[0077] In addition, when the logical instructions in the above-mentioned memory 330 can be implemented in the form of software functional units and sold or used as an independent product, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0078] On the other hand, the present invention also provides a computer program product. The computer program product includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the processing method for the transfer result between multi-path radiation reference stars provided by the above-mentioned various methods. The method includes: through the multi-path radiation reference star transfer method, obtaining and calculating the correction amounts of spectral radiance at the top of the atmosphere and uncertainties of the target spaceborne instrument. One reference star corresponds to one correction amount of spectral radiance at the top of the atmosphere and one uncertainty; excluding the outliers in the The top-of-atmosphere spectral radiance correction and the corresponding uncertainty; in the When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method; When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method; wherein the first condition is that any uncertainty and the The difference between the average uncertainties of the uncertainties is smaller than the first threshold. Greater than or equal to ,and , All are positive integers.

[0079] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which is implemented when the computer program is executed by a processor to execute the processing method for the multi-path radiation reference inter-satellite transfer result provided by the above methods, the method comprising: obtaining and calculating the target satellite-borne instrument through the multi-path radiation reference inter-satellite transfer method; The top-of-atmosphere spectral radiance correction and uncertainty, one reference star corresponds to one atmospheric top spectral radiance correction and one uncertainty; excluding the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty; in the When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method; When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method; wherein the first condition is that any uncertainty and the The difference between the average uncertainties of the uncertainties is smaller than the first threshold. Greater than or equal to ,and , All are positive integers.

[0080] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0081] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0082] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for processing multi-path radiation benchmark intersatellite transfer results, characterized in that: include: The target satellite instrument’s inter-satellite transfer method is used to obtain and calculate the target satellite instrument’s The top-of-atmosphere spectral radiance correction and One reference star corresponds to one atmospheric top spectral radiance correction and one uncertainty. Eliminate the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty; In the When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method; In the When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method; The first condition is any uncertainty and the The difference between the average uncertainties of the uncertainties is smaller than the first threshold. Greater than or equal to ,and , All are positive integers; mentioned in When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method, including: By comparing each uncertainty with the The corresponding cutoff uncertainty is determined based on the median of the uncertainties, and the original uncertainty is adjusted based on the cutoff uncertainty to obtain the adjusted target uncertainty; Determine a weight coefficient of each top-of-atmosphere spectral radiance correction amount according to the adjusted target uncertainty, wherein the weight coefficient is inversely proportional to the square of the target uncertainty; Determining a weighted average top-of-atmosphere spectral radiance correction amount based on the weight coefficient and the top-of-atmosphere spectral radiance correction amount; Determine the uncertainty of the weighted average top-of-atmosphere spectral radiance correction value based on the uncertainty of each top-of-atmosphere spectral radiance correction value; Determining the weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount as the target transfer result; mentioned in When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method, including: Calculate the The average top-of-atmosphere spectral radiance correction of the top-of-atmosphere spectral radiance correction and the The standard uncertainty of an uncertainty; The average top-of-atmosphere spectral radiance correction value and the standard uncertainty are determined as the target transfer result.

2. The method for processing the multi-path radiation benchmark inter-satellite transfer result according to claim 1, characterized in that: The elimination of The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainties, including: based on The judgment criteria are as follows Remove outliers from the atmospheric top spectral radiance correction and obtain The top-of-atmosphere spectral radiance correction and the corresponding An uncertainty.

3. A processing device for multi-path radiation benchmark intersatellite transfer results, characterized in that: include: Acquisition module and processing module; The acquisition module is used to acquire and calculate the target satellite-borne instrument through a multi-path radiation reference inter-satellite transfer method. The top-of-atmosphere spectral radiance correction and One reference star corresponds to one atmospheric top spectral radiance correction and one uncertainty. The processing module is used to remove the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty; in the When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method; In the When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method; The first condition is any uncertainty and the The difference between the average uncertainties of the uncertainties is smaller than the first threshold. Greater than or equal to ,and , All are positive integers; The processing module is used to compare each uncertainty with the The corresponding cutoff uncertainty is determined based on the size of the median of the uncertainties, and the original uncertainty is adjusted based on the cutoff uncertainty to obtain the adjusted target uncertainty; the weight coefficient of each top-of-atmosphere spectral radiance correction amount is determined based on the adjusted target uncertainty, and the weight coefficient is inversely proportional to the square of the target uncertainty; the weighted average top-of-atmosphere spectral radiance correction amount is determined based on the weight coefficient and the top-of-atmosphere spectral radiance correction amount; the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount is determined based on the uncertainty of each top-of-atmosphere spectral radiance correction amount; the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount and the weighted average top-of-atmosphere spectral radiance correction amount is determined as the target transfer result; or, Calculate the The average top-of-atmosphere spectral radiance correction of the top-of-atmosphere spectral radiance correction and the The standard uncertainty of an uncertainty; the average top-of-atmosphere spectral radiance correction amount and the standard uncertainty are determined as the target transfer result.

4. The processing device for multi-path radiation benchmark inter-satellite transfer results according to claim 3, characterized in that: The processing module is used to The judgment criteria are as follows Remove outliers from the atmospheric top spectral radiance correction and obtain The top-of-atmosphere spectral radiance correction and the corresponding An uncertainty.

5. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method for processing the multi-path radiation reference inter-satellite transfer result as claimed in any one of claims 1 to 2 is implemented.

6. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for processing the multi-path radiation reference inter-satellite transfer result as claimed in any one of claims 1 to 2 is implemented.

7. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method for processing the multi-path radiation reference inter-satellite transfer result as claimed in any one of claims 1 to 2 is implemented.

Images