Method and device for processing transmission result between multi-path radiation reference satellites

Abstract

The invention provides a processing method and device for a multi-path radiation reference inter-satellite transmission result, and is applied to the technical field of space radiation measurement. The method comprises the steps that # imgabs 0 # atmospheric top spectral radiance correction and # imgabs 1 # uncertainty of a target spaceborne instrument are obtained and calculated through a multi-path radiation reference inter-satellite transmission method, and one reference satellite corresponds to one atmospheric top spectral radiance correction and one uncertainty; the outliers in the # imgabs 2 # atmospheric top spectrum radiance correction values are removed, and # imgabs 3 # atmospheric top spectrum radiance correction values and corresponding # imgabs 4 # uncertainties are obtained; under the condition that the # imgabs5 uncertainty does not meet a first condition, calculating a target transmission result based on a weighted average method; if the # imgabs6 uncertainty satisfies the first condition, a target delivery result is calculated on the basis of an arithmetic average method.

Description

Technical Field

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

[0002] Space radiation measurement reference satellites are the key to achieving the unification of remote sensing constellation radiation measurement benchmarks and overall improvement of high-precision quantitative remote sensing capabilities. They are also the high ground of quantitative remote sensing technology that major aerospace remote sensing countries are competing for. Using reference satellites to transfer radiation values ​​is currently the most accurate on-orbit calibration method.

[0003] In the prior art, the radiation benchmark intersatellite transfer method usually selects the atmospheric top spectral radiation brightness data measured by the reference star sensor with complete calibration equipment and high calibration accuracy as a reference, and performs spatial matching and spectral matching on the near-orbit overlapping data of the reference instrument and the instrument to be calibrated, which are relatively close in time and space, to obtain the atmospheric top spectral radiation brightness correction or calibration coefficient 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 benchmark transfer refers to the use of different reference stars and different transfer methods to transfer the intersatellite radiation benchmark values ​​to the target star, and obtain a series of atmospheric top spectral radiation brightness corrections and their uncertainties for the target satellite-borne instrument.

[0004] Due to the different reference satellites and transfer methods, the correlation between the series of transferred top-of-atmosphere spectral radiance corrections is weak. How to scientifically and rationally analyze and process the series of radiation top-of-atmosphere spectral radiance corrections and their uncertainties obtained by multi-path radiation reference inter-satellite transfer to obtain the final radiation top-of-atmosphere spectral radiance corrections and uncertainties of the target satellite-borne instrument is a common problem encountered by meteorological, oceanographic, high-resolution land and other series of satellites when using radiation reference value transfer. Summary of the invention

[0005] The present invention provides a method and device for processing the results of multi-path radiation benchmark intersatellite transfer, which is used to solve the problem of how to scientifically and reasonably analyze and process the series of radiation atmospheric top spectrum radiation brightness corrections and their uncertainties obtained by multi-path radiation benchmark intersatellite transfer.

[0006] The present invention provides a method for processing a multi-path radiation benchmark inter-satellite transfer result, comprising: obtaining and calculating the target satellite-borne instrument through a 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; 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.

[0007] According to a method for processing multi-path radiation benchmark inter-satellite transfer results provided by the present invention, the elimination of the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding Uncertainties include: 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.

[0008] According to a method for processing multi-path radiation benchmark inter-satellite transfer results provided by the present invention, the 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 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 according to 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 weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount are determined as the target transfer result.

[0009] According to a method for processing multi-path radiation benchmark inter-satellite transfer results provided by the present invention, 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 The corresponding cutoff uncertainty is determined by the size of the median of each uncertainty, 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 according to 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 according to the uncertainty of each top-of-atmosphere spectral radiance correction amount; the weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount are determined as the target transfer result.

[0013] 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 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.

[0014] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for processing the multi-path radiation benchmark intersatellite transfer result as described in any one of the above is implemented.

[0015] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a method for processing a 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, comprising a computer program, wherein when the computer program is executed by a processor, the computer program implements any of the above-mentioned methods for processing the multi-path radiation reference inter-satellite transfer results.

[0017] The method and device for processing the multi-path radiation benchmark intersatellite transfer results provided by the present invention can first eliminate The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty, and then according to The method of calculating the target transfer result can be selected based on the difference of the uncertainty. In this way, the series of radiation top-of-atmosphere spectrum radiance corrections and their uncertainties obtained by transferring between multi-path radiation benchmark stars can be scientifically and reasonably analyzed and processed in a unified way, so as to obtain the final radiation top-of-atmosphere spectrum radiance corrections and uncertainties of the target satellite-borne instrument, thereby improving the scientific rationality of the multi-path transfer fusion results. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the technical solutions in the present invention or the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0019] Figure 1 It is a flow chart of a method for processing multi-path radiation benchmark inter-satellite transfer results provided by the present invention; Figure 2 It is a structural schematic diagram of a device for processing multi-path radiation benchmark inter-satellite transfer results provided by the present invention; Figure 3 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in this application will be clearly and completely described below in conjunction with the drawings in this application. Obviously, the described embodiments are part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0021] 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 descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0022] It should be noted that, in this article, the terms "comprise", "include" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises one..." does not exclude the presence of other identical elements in the process, method, article or device including the element. In addition, it should be noted that the scope of the method and device in the embodiment 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 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 also be added, omitted, or combined. In addition, the features described with reference to certain examples may be combined in other examples.

[0023] In order to facilitate the clear description of the technical solutions of the embodiments of the present application, in the embodiments of the present application, the words "first", "second", etc. are used to distinguish the same items or similar items with basically the same functions and effects. Those skilled in the art can understand that the words "first", "second", etc. are not limiting the quantity and execution order.

[0024] The embodiments of the present application describe some exemplary embodiments for the purpose of explanation. It should be understood that the present application can be implemented in other ways that are not specifically shown in the drawings.

[0025] like Figure 1 As shown, the embodiment of the present application provides a method for processing a multi-path radiation benchmark inter-satellite transfer result, and the method for processing a multi-path radiation benchmark inter-satellite transfer result can be applied to a processing device for a multi-path radiation benchmark inter-satellite transfer result. The method for processing a multi-path radiation benchmark inter-satellite transfer result can include S101-S104: S101, a processing device for multi-path radiation reference inter-satellite transfer results obtains and calculates the target satellite-borne instrument through a multi-path radiation reference inter-satellite transfer method. The top-of-atmosphere spectral radiance correction and An uncertainty.

[0026] Among them, one reference star corresponds to one atmospheric top spectral radiation brightness correction and one uncertainty.

[0027] In the field of remote sensing, in order to accurately measure the radiation value of the target star, it is necessary to use a 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 multiple ways of radiation reference transfer results. For example, when satellite remote sensing observes meteorological, oceanic, and land targets, radiation reference transfer may be performed through instruments carried by satellites in different orbital positions, using different optical band observation methods or different calibration models.

[0028] There will be errors when the onboard instrument of an uncalibrated target star measures the radiation value. By transferring the accurate radiation value of the reference star to the target star, the The top-of-atmosphere spectral radiance correction is used to correct the target star measurement value. Each top-of-atmosphere spectral radiance correction corresponds to an uncertainty. These uncertainties come from many aspects, including the measurement uncertainty of the reference star radiation value itself, the instrument error due to the differences in time domain, space domain, spectrum domain, and angular domain during the transmission process, the influence of environmental factors, and the uncertainty of the model. The uncertainty reflects the credibility of the top-of-atmosphere spectral radiance correction. Smaller uncertainty means that the top-of-atmosphere spectral radiance correction is more reliable.

[0029] S102, the processing device for multi-path radiation benchmark inter-satellite transmission results removes the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding An uncertainty.

[0030] Optionally, the processing device for the multi-path radiation benchmark inter-satellite transfer result can be 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.

[0031] Specifically, due to some sudden and non-repeatable factors, Outliers may appear in the top-of-atmosphere spectral radiance correction. For example, when performing remote sensing observations to obtain radiation benchmarks to transfer top-of-atmosphere spectral radiance corrections, it may be affected by factors such as cosmic ray bursts, instantaneous loss of satellite attitude control, and temporary electromagnetic interference. These factors occur by chance and are difficult to repeat under normal experimental conditions. They will cause some top-of-atmosphere spectral radiance corrections to deviate from the normal range, so they need to be eliminated, which includes the following steps: Step 1: The processing device for the multi-path radiation benchmark intersatellite transfer result can first calculate The average value of the top-of-atmosphere spectral radiance correction .

[0032] calculate The average value of the atmospheric top spectral radiation brightness correction is: .

[0033] Step 2: The processing device for the multi-path radiation benchmark intersatellite transfer result calculates the residual of each atmospheric top spectrum radiation brightness correction value based on the average value. .

[0034] The residual represents the difference between the spectral radiance correction value of each top of atmosphere and the average value. The residual of the spectral radiance correction value of each top of atmosphere is calculated using the formula: .

[0035] Step 3: The processing device for the intersatellite transfer results of the multi-path radiation benchmark calculates the experimental standard deviation of the spectral radiation brightness correction value of each top of atmosphere based on the average value. .

[0036] The experimental standard deviation measures the degree of dispersion of the top-of-atmosphere spectral radiance corrections, which comprehensively considers the difference between each top-of-atmosphere spectral radiance correction and the average value. The larger the standard deviation, the greater the dispersion between the top-of-atmosphere spectral radiance corrections. The experimental standard deviation of each top-of-atmosphere spectral radiance correction is calculated using the formula: .

[0037] Step 4: The processing device of the multi-path radiation benchmark intersatellite transfer result removes outliers based on the residual and experimental standard deviation.

[0038] Under 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. It is the residual with the largest absolute value among all residuals. When it satisfies the condition , then the atmospheric top spectral radiation brightness correction is considered to be an outlier and should be removed. This is based on the normal distribution of “ " principle. For example, if the residual absolute value of a top-of-atmospheric spectral radiance correction is much larger than three times the experimental standard deviation, then this top-of-atmospheric spectral radiance correction is likely to be affected by the sudden abnormal factors mentioned above, so it will be removed as an outlier to ensure the accuracy of subsequent data processing and analysis.

[0039] After removing outliers, the above procedure should be used repeatedly until the top-of-atmosphere spectral radiance correction value no longer contains outliers. This is because after removing an outlier, the remaining statistics such as the mean, residual, and standard deviation of the top-of-atmosphere spectral radiance correction value will change, which may cause the top-of-atmosphere spectral radiance correction value that was not originally judged as an outlier to become an outlier under the new statistical environment. Therefore, this process needs to be repeated to ensure that the top-of-atmosphere spectral radiance correction values ​​that are finally left are relatively reasonable and reliable data, thereby providing an accurate data basis for subsequent radiation benchmark transfer and data fusion operations.

[0040] S103, a processing device for multi-path radiation benchmark inter-satellite transmission results is described in When the uncertainty does not meet the first condition, the target transfer result is calculated based on the weighted average method.

[0041] 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.

[0042] Optionally, the processing device for the multi-path radiation reference intersatellite transfer result calculates the target transfer result based on the weighted average method, including: by comparing 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 brightness correction amount is determined based on the adjusted target uncertainty, and the weight coefficient is inversely proportional to the square of the adjusted target uncertainty; the weighted average top-of-atmosphere spectral radiance brightness correction amount is determined based on the weight coefficient and the top-of-atmosphere spectral radiance brightness correction amount; the uncertainty of the weighted average amount is determined based on the uncertainty of each top-of-atmosphere spectral radiance brightness correction amount; the uncertainty of the weighted average top-of-atmosphere spectral radiance brightness correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance brightness correction amount are determined as the target transfer result.

[0043] Specifically, after the multi-path radiation benchmark is transferred, for the target star at wavelength The TOA spectral radiance at the location is transferred to the TOA spectral radiance correction , and their corresponding standard uncertainties are significantly different. This difference makes it impossible to simply take the arithmetic average of the top-of-atmosphere spectral radiance corrections to fuse the results.

[0044] The weighted average method is considered to be 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 value (reflected by uncertainty), so as to more reasonably integrate the information of each top-of-atmosphere spectral radiance correction value.

[0045] In order to avoid some too low and unreasonable uncertainties from having an inappropriate impact on the final fusion result, this application proposes the concept of cutoff uncertainty.

[0046] when hour, That is, when the uncertainty Less than or equal to the median uncertainty When the cutoff uncertainty is Set to the average of all uncertainties. This can prevent the atmospheric top spectral radiance correction with too small uncertainty from occupying unreasonable weight in the fusion process.

[0047] The original uncertainty is adjusted according to the cut-off uncertainty to obtain the adjusted target uncertainty. , then the target uncertainty If the uncertainty , then the target uncertainty .

[0048] Calculate the spectral radiance correction value of each top of atmosphere according to the target uncertainty The weight coefficient 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, the larger its weight coefficient and the more important it is in the fusion process.

[0049] Weighted average top-of-atmosphere spectral radiance correction for target TOA radiance By multiplying each top-of-atmosphere spectral radiance correction by its weight coefficient and summing them up, the weighted average top-of-atmosphere spectral radiance correction can be obtained, which integrates the information of each top-of-atmosphere spectral radiance correction.

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

[0051] The spectral radiance correction value of each target top of atmosphere is Correction to the weighted average top-of-atmosphere spectral radiance The difference is expressed as equivalence The smaller the equivalence, the closer the target atmospheric top spectral radiance correction is to the weighted average correction value, and the more accurate the result is.

[0052] When When the top-of-atmosphere spectral radiance correction value participates in the calculation of the weighted average correction value, Uncertainty ; When the When the spectral radiance correction value of the top of atmosphere does not participate in the calculation of the weighted average correction value, These uncertainties can help assess the reliability of the difference between the top-of-atmosphere spectral radiance correction and the weighted average correction.

[0053] The weighted average method needs to use the chi-square statistic to perform a compatibility test on the top-of-atmosphere spectral radiance correction. This is to ensure that the top-of-atmosphere spectral radiance corrections 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 results. If the top-of-atmosphere spectral radiance corrections are incompatible, it may cause deviations in the fusion results, so compatibility testing is a critical step.

[0054] S104, a processing device for multi-path radiation benchmark inter-satellite transmission results is described in When the uncertainty satisfies the first condition, the target transfer result is calculated based on the arithmetic mean method.

[0055] Optionally, the processing device for multi-path radiation reference intersatellite transfer results calculates the target transfer results 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 is synthesized from the uncertainties; and the average top-of-atmosphere spectral radiance correction amount and the standard uncertainty are determined as the target transfer result.

[0056] Specifically, the wavelength of the target star is obtained by multi-path reference transfer. The top of atmosphere (TOA) spectral radiance transfers the top of atmosphere spectral radiance correction and its corresponding standard uncertainty Finally, when the uncertainties of these atmospheric top spectral radiance corrections are similar, the arithmetic mean method can be considered for fusion. Here, “similar uncertainties” are defined as the maximum deviation between the uncertainty and the average uncertainty does not exceed the first threshold, which can be .

[0057] For example, in the process of obtaining the top-of-atmosphere spectral radiance correction of the target star radiance through remote sensing observation, if the uncertainty of the top-of-atmosphere spectral radiance correction obtained by different channels is not much different, this may mean that the reliability and accuracy of each transmission channel are roughly the same, and no one channel has a significant uncertainty. In this case, each top-of-atmosphere spectral radiance correction can be regarded as equally important in the fusion process.

[0058] Calculated using arithmetic mean method Average transferred top-of-atmosphere spectral radiance correction for multipath radiance benchmark transfer , the formula is The meaning of this formula is to add up all the top-of-atmosphere spectral radiation brightness corrections, and then divide them by the number of top-of-atmosphere spectral radiation brightness corrections, and the result is the arithmetic mean of these top-of-atmosphere spectral radiation brightness corrections.

[0059] calculate Standard uncertainty of the mean correction value transferred by the multipath radiation reference , the formula is This formula takes into account the uncertainty of each top-of-atmosphere spectral radiance correction. By taking the square root of the sum of all uncertainties and dividing it by the number of top-of-atmosphere spectral radiance corrections, the standard uncertainty of the average correction value is obtained. It reflects the reliability of the top-of-atmosphere spectral radiance correction after arithmetic averaging.

[0060] If the uncertainty of all top-of-atmosphere spectral radiance corrections is the same, they are , then the standard uncertainty The formula can be simplified to .

[0061] The arithmetic mean method, like the weighted mean method, also requires the use of the chi-square statistic to perform a compatibility test on the top-of-atmosphere spectral radiance corrections. This is because even if the uncertainties are similar, there may be some potential inconsistencies between the top-of-atmosphere spectral radiance corrections. The compatibility test can ensure that these top-of-atmosphere spectral radiance corrections are compatible with each other during the fusion process, that is, the differences between them are within a reasonable range, thereby ensuring the rationality and reliability of the arithmetic mean fusion results. If the top-of-atmosphere spectral radiance corrections are incompatible, the results of the arithmetic mean may be affected, resulting in inaccurate average transferred top-of-atmosphere spectral radiance corrections and corresponding uncertainties.

[0062] In the embodiment of the present application, the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty, and then according to The method of calculating the target transfer result can be selected based on the difference of the uncertainty. In this way, the series of radiation top-of-atmosphere spectrum radiance corrections and their respective uncertainties obtained by multi-path radiation reference star transfer can be scientifically and reasonably analyzed and processed in a unified way, so as to obtain the final radiation top-of-atmosphere spectrum radiance correction and uncertainty of the target satellite-borne instrument, thereby improving the scientific rationality of the multi-path transfer fusion results.

[0063] The above mainly introduces the solution provided by the embodiment of the present application from the perspective of the method. In order to realize the above functions, it includes hardware structures and / or software modules corresponding to the execution of each function. Those skilled in the art should easily realize that, in combination with the units and algorithm steps of each example described in the embodiment disclosed herein, the embodiment of the present application can be implemented in the form of hardware or a combination of hardware and computer software. Whether a 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 and technical personnel 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.

[0064] It should be noted that the device in the embodiment of the present application includes a virtual device and a physical device. The virtual device can be a processing device for multi-path radiation benchmark inter-satellite transfer results, and the physical device can include an electronic device, a computer storage medium, and a computer program product.

[0065] The method for processing the multi-path radiation benchmark inter-satellite transfer result provided in the embodiment of the present application can be executed by a processing device for the multi-path radiation benchmark inter-satellite transfer result, or a control module for processing the multi-path radiation benchmark inter-satellite transfer result in the processing device for the multi-path radiation benchmark inter-satellite transfer result. In the embodiment of the present application, the processing method for processing the multi-path radiation benchmark inter-satellite transfer result is executed by the processing device for the multi-path radiation benchmark inter-satellite transfer result as an example to illustrate the processing device for processing the multi-path radiation benchmark inter-satellite transfer result provided in the embodiment of the present application.

[0066] It should be noted that the embodiment of the present application can divide the functional modules of the processing device for the multi-path radiation benchmark intersatellite transfer results according to the above method example. 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 software functional modules. Optionally, the division of modules in the embodiment of the present application is schematic and is only a logical functional division. There may be other division methods in actual implementation.

[0067] 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.

[0068] Optionally, the processing module 202 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.

[0069] Optionally, the processing module 202 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 according to 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 weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount are determined as the target transfer result.

[0070] Optionally, the processing module 202 is used to 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.

[0071] In the embodiment of the present application, the The outliers in the top-of-atmosphere spectral radiance correction are obtained The top-of-atmosphere spectral radiance correction and the corresponding uncertainty, and then according to The method of calculating the target transfer result can be selected based on the difference of the uncertainty. In this way, the series of radiation top-of-atmosphere spectrum radiance corrections and their respective uncertainties obtained by multi-path radiation reference star transfer can be scientifically and reasonably analyzed and processed in a unified way, so as to obtain the final radiation top-of-atmosphere spectrum radiance correction and uncertainty of the target satellite-borne instrument, thereby improving the scientific rationality of the multi-path transfer fusion results.

[0072] Figure 3 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 3 As shown, the electronic device may include: a processor 310, a communications interface 320, a memory 330 and a communications bus 340, wherein the processor 310, the communications interface 320 and the memory 330 communicate with each other through the communications bus 340. The processor 310 may call the logic instructions in the memory 330 to execute the method for processing the multi-path radiation benchmark inter-satellite transfer result, the method comprising: obtaining and calculating 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 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.

[0073] In addition, the logic instructions in the above-mentioned memory 330 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.

[0074] 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 in a non-transitory computer-readable storage medium, and when the computer program is executed by a processor, the computer can execute the processing method of 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.

[0075] 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.

[0076] 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.

[0077] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for 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.

[0078] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions 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.

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. The method for processing the multi-path radiation benchmark inter-satellite transfer result according to claim 1, characterized in that: 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; The weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount are determined as the target transfer result.

4. The method for processing the multi-path radiation benchmark inter-satellite transfer result according to claim 1, characterized in that: 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.

5. 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.

6. The processing device for multi-path radiation benchmark inter-satellite transfer results according to claim 5, 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.

7. The processing device for multi-path radiation benchmark inter-satellite transfer results according to claim 5, characterized in that: The processing module is used to compare each uncertainty with the The size of the median of the uncertainties determines the corresponding cutoff uncertainty, 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; A 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 according to the uncertainty of each top-of-atmosphere spectral radiance correction amount; the weighted average top-of-atmosphere spectral radiance correction amount and the uncertainty of the weighted average top-of-atmosphere spectral radiance correction amount are determined as the target transfer result.

8. 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 4 is implemented.

9. 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 4 is implemented.

10. 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 4 is implemented.

Images