SO2 product background error correction method and device and electronic equipment

Abstract

The invention provides an SO2 product background error correction method and device and electronic equipment, and the method comprises the steps: obtaining a whole layer of SO2 vertical column total amount product; dividing the data blocks into data blocks with different integration times, screening effective pixels, and estimating a mean value of each rail-crossing position in the rail-crossing direction; separating the SO2 vertical column total amount product from the mean value of the rail penetrating position to obtain a first correction result; based on the first correction result, the concentration threshold value and the sliding window, eliminating the pixels exceeding the concentration threshold value; and estimating the background error matrix and continuously separating to obtain a secondary correction result, and carrying out loop iteration until the two continuous correction results meet a preset condition. According to the SO2 product background error correction method, the SO2 product background error correction device and the electronic equipment provided by the embodiment of the invention, the data blocks are divided into different data blocks according to the integration time, and two background error schemes are used for iteratively processing each data block with different integration time, so that an obvious seam phenomenon of a satellite-borne atmospheric component SO2 inversion result at an integration time change part can be effectively solved.

Description

Technical Field

[0001] The present invention relates to the field of atmospheric remote sensing technology, and in particular to a method, device, electronic equipment and computer-readable storage medium for correcting background errors of SO2 products. Background Art

[0002] The trace gas SO2 column concentration obtained by direct inversion based on satellite-observed irradiance and radiance is always affected by factors such as calibration error and inversion model error, which leads to systematic errors in the inversion results. Since the residual calibration error in the satellite-observed irradiance and radiance data is difficult to monitor and correct, the satellite-borne SO2 column concentration product often corrects the initial SO2 column concentration inversion results based on the background error correction method. The existing background error correction method has achieved good results for payloads such as GOME2 (Global Ozone Monitoring Experiment-2) and OMI (Ozone Monitoring Instrument), but for payloads with discontinuous orbital observation integration time (such as the FY-3F ultraviolet hyperspectral ozone detector (Ozone Monitoring Suite, OMS)), the current background error correction method cannot effectively correct the systematic deviations of the SO2 inversion results in the through-track and along-track directions, resulting in obvious seams in the trace gas SO2 inversion results at the integration time change. Summary of the invention

[0003] In order to solve the existing technical problems, the embodiments of the present invention provide a method, device, electronic device and computer-readable storage medium for correcting background error of SO2 products.

[0004] In a first aspect, an embodiment of the present invention provides a method for correcting background errors of SO2 products, comprising: obtaining a whole layer of SO2 vertical column total product, wherein the whole layer of SO2 vertical column total product includes SO2 vertical column total products of multiple tracks, each of which includes multiple scanning lines; dividing the SO2 vertical column total product of each track into data blocks with different integration times; screening effective pixels for each data block with the same integration time, and estimating the mean of each track crossing position in the track crossing direction; cyclically processing each scanning line in the data block with the same integration time, separating the SO2 vertical column total product from the mean of the track crossing position. , obtain the total amount of SO2 vertical column after the first correction; based on the total amount of SO2 vertical column after the first correction and the preset SO2 concentration threshold, exclude the pixels exceeding the SO2 concentration threshold; use the sliding window to estimate the background error matrix corresponding to each of the scanning lines, perform difference separation on the total amount of SO2 vertical column after the first correction and the background error, and obtain the total amount of SO2 vertical column after the second correction; based on the corrected total amount of SO2 vertical column product, the SO2 concentration threshold and the sliding window, perform a cyclic iteration until the total amount of SO2 vertical column products obtained by two consecutive corrections meet the preset conditions.

[0005] Optionally, obtaining the total amount product of the vertical column of SO2 for the entire layer includes: taking the observation data of the irradiance and radiance of the ultraviolet hyperspectral satellite payload as input, and obtaining the total amount product of the vertical column of SO2 for the entire layer with the help of a forward model and a classical differential absorption theory inversion model.

[0006] Optionally, the SO2 vertical column total product of each of the tracks is divided into data blocks with different integration times, including: based on the observation integration time of each observation pixel, the SO2 vertical column total product of each of the tracks is divided into data blocks with different integration times.

[0007] Optionally, screening valid pixels for each data block with the same integration time includes: screening valid pixels for each data block with the same integration time based on the same load cloud product and quality identification code.

[0008] Optionally, the preset SO2 concentration threshold is 2DU; the size of the sliding window changes with the actual spatial resolution of the satellite.

[0009] Optionally, the preset condition includes: a relative deviation between the total vertical column amount products of SO2 obtained from two consecutive calibrations is less than or equal to a preset percentage.

[0010] In a second aspect, an embodiment of the present invention further provides a SO2 product background error correction device, comprising: a data acquisition module, a first correction module, a second correction module and a loop iteration module; the data acquisition module is used to obtain the total amount product of the entire layer of SO2 vertical column, and the total amount product of the entire layer of SO2 vertical column includes the total amount product of the SO2 vertical column of multiple tracks, and each track contains multiple scanning lines; the first correction module is used to divide the total amount product of the SO2 vertical column of each track into data blocks with different integration times; for each data block with the same integration time, effective pixels are selected, and the mean value of each track crossing position in the track crossing direction is estimated; each scanning line in the data block with the same integration time is looped to convert the total amount product of the SO2 vertical column into a data block with different integration times; the second correction module is used to divide the total amount product of the SO2 vertical column of each track into a data block with different integration times; the second correction module is used to divide the total amount product of the SO2 vertical column ... The first correction module is used to separate the total amount of SO2 vertical column product from the mean of the crossing position to obtain the total amount of SO2 vertical column after the first correction; the second correction module is used to exclude pixels exceeding the SO2 concentration threshold based on the total amount of SO2 vertical column product after the first correction and the preset SO2 concentration threshold; the background error matrix corresponding to each of the scanning lines is estimated by using a sliding window, and the total amount of SO2 vertical column product after the first correction is separated from the background error to obtain the total amount of SO2 vertical column product after the second correction; the loop iteration module is used to perform loop iteration based on the corrected total amount of SO2 vertical column product, the SO2 concentration threshold and the sliding window until the total amount of SO2 vertical column products obtained by two consecutive corrections meet the preset conditions.

[0011] In a third aspect, an embodiment of the present invention provides an electronic device, including a processor and a memory, wherein the memory stores a computer program, and the processor executes the computer program stored in the memory, and when the computer program is executed by the processor, the SO2 product background error correction method described in the first aspect is implemented.

[0012] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the SO2 product background error correction method described in the first aspect is implemented.

[0013] In a fifth aspect, the present application further provides a computer program product, including a computer program. When the computer program is executed, it can implement the SO2 product background error correction method described in the first aspect or any possible design method of the first aspect.

[0014] The SO2 product background error correction method, device, electronic device and computer-readable storage medium provided by the embodiments of the present invention are compared with the existing background error correction method. The scheme is divided into different data blocks according to the integration time, and two background error schemes are used for iterative processing for each different integration time data block. It can effectively solve the obvious seam phenomenon of the trace gas SO2 inversion result at the integration time change, and can be more suitable for orbital observation of loads with discontinuous integration time, and improve the accuracy and application range of the satellite-borne SO2 vertical column total inversion product. In addition, this method has a certain reference value for other trace gas inversion algorithms. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the background technology, the drawings required for use in the embodiments of the present invention or the background technology will be described below.

[0016] Figure 1 The results of the correction using the existing background error correction method are shown in FIG. Figure 1 ;

[0017] Figure 2 The results of the correction using the existing background error correction method are shown in FIG. Figure 2 ;

[0018] Figure 3 A flow chart of a method for correcting background error of a SO2 product provided by an embodiment of the present invention is shown;

[0019] Figure 4 A schematic diagram showing the result of correction using the SO2 product background error correction method provided by an embodiment of the present invention;

[0020] Figure 5 A schematic structural diagram of a SO2 product background error correction device provided by an embodiment of the present invention is shown;

[0021] Figure 6 A schematic structural diagram of an electronic device for executing a method for correcting a background error of an SO2 product provided by an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0022] SO2 (i.e. sulfur dioxide) is the most common sulfur oxide. SO2 is one of the trace gas components in the atmosphere and is also an important indicator for evaluating urban air quality. Changes in the concentration of SO2 in the atmosphere and changes in its temporal and spatial distribution will have a significant impact on the environment and climate. The natural source of SO2 emissions is mainly volcanic eruptions, but human activities will also cause the content of SO2 in the atmosphere to increase. The main sources of human emissions are the combustion of coal, oil, and natural gas, mineral smelting, industrial waste gas from the production of sulfuric acid and phosphate fertilizers, and exhaust from motor vehicles.

[0023] Since SO2 has obvious absorption characteristics in the ultraviolet band, especially in the wavelength range of 300-400nm, it is suitable for the detection of low-concentration SO2, such as volcanic gas or atmospheric background value. Therefore, since the large amount of SO2 produced by the eruption of El Chichón volcano was observed for the first time in the ultraviolet band in 1982, ultraviolet hyperspectral satellite payloads have been used for quantitative global SO2 concentration monitoring. Various ultraviolet hyperspectral satellite payloads launched later can realize global SO2 emission monitoring.

[0024] However, for most of the ultraviolet hyperspectral trace gas satellite payload observation data, the residual systematic measurement errors (including the primary solar and earth observation radiation data) after L1 radiation and spectral calibration processing lead to the inhomogeneity of satellite payload solar irradiance and earth observation radiance observation values ​​in the spectral dimension and spatial dimension. This residual calibration error is difficult to monitor and correct, but because the signal-to-noise ratio of satellite-borne observations in the ultraviolet band is low and the trace gas SO2 absorption effect is only a weak signal in the satellite detection signal, these errors are transmitted to the atmospheric component SO2 inversion value through the inversion algorithm, which will further amplify the errors in the along-track and through-track directions, such as along-track stripes, asymmetric track through-track directions, and regional systematic high or low errors. In particular, when the integration time of the same track changes, it leads to obvious seams. This background error in SO2 products is different from the influence of random noise. The trace gas content information is often mixed with the systematic error, which seriously affects the application scope of trace gas remote sensing products. For example, most of the systematic stripe deviations exist along the track, and the high values ​​of the along-track stripes contain the atmospheric trace gas concentration absorption information. Therefore, in the process of SO2 satellite remote sensing inversion, it is necessary to correct the background error of the inversion results to reduce the error of SO2 concentration inversion results.

[0025] At present, the reference area method and the sliding median residual correction method are often used to correct the cross-track and latitude errors. However, the inventors found that although the existing error correction method has achieved good results for payloads such as GOME-2 (Global Ozone Monitoring Experiment-2, the second generation of the global ozone monitor) and OMI ozone monitor (Ozone Monitoring Instrument), for payloads with discontinuous orbital observation integration time (taking the FY-3F ultraviolet hyperspectral ozone detector (Ozone Monitoring Suite, OMS) as an example), after applying the above sliding background error correction method, the trace gas SO2 inversion results have obvious seams at the integration time changes, such as Figure 1 and Figure 2As shown. To address this problem, the present invention proposes a SO2 product background error correction method, which is a trace gas inversion product background error correction algorithm for orbital integration time changes, thereby obtaining high-quality satellite remote sensing trace gas products. At the same time, the present invention also has a certain reference value for other trace gas inversion algorithms.

[0026] The embodiments of the present invention are described below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0027] Figure 3 The flowchart of a method for correcting the background error of a SO2 product provided by an embodiment of the present invention is shown. The method is a method for correcting the background error of a SO2 product for the variation of the onboard integration time. Figure 3 As shown, the method includes the following steps 301-304.

[0028] Step 301: Obtain the total amount product of the entire layer of SO2 vertical column, which includes the total amount product of the SO2 vertical column of multiple tracks, and each track includes multiple scanning lines.

[0029] Among them, the remote sensing product of the total vertical column of SO2 in the whole layer refers to the vertical column total data (Vertical Column Density, VCD) of SO2 in the atmosphere obtained by remote sensing technology, which can be understood as the total content of SO2 from the surface to the top of the atmosphere per unit area, for example, it can be recorded as VCD0. Furthermore, the total vertical column product of SO2 in the whole layer includes the total vertical column products of SO2 in multiple orbits. The orbit mentioned here can actually be understood as an observation orbit synchronized with the sun observed by a polar-orbiting satellite, with about 14-15 orbits a day. Furthermore, each orbit has multiple scanning lines from the South Pole to the North Pole.

[0030] Step 302: Divide the total vertical column product of SO2 of each track into data blocks with different integration times; select valid pixels for each data block with the same integration time, and estimate the mean of each track-crossing position in the track-crossing direction; cyclically process each scanning line in the data block with the same integration time, separate the total vertical column product of SO2 from the mean of the track-crossing position, and obtain the total vertical column product of SO2 after the first correction. Among them, the integration time refers to the length of time that the satellite sensor accumulates and receives signals (such as photons and radiation energy) within a specific time. It can be understood as the time for the sensor to "expose" a target area. In the embodiment of the present invention, the total vertical column product VCD0 of SO2 obtained by inverting each orbit of polar-orbiting satellite observation is divided into data blocks with different integration times according to the integration time. For each data block with the same integration time, firstly select the valid pixels of all scanning lines. In the embodiment of the present invention, the pixels with cloud cover less than 0.3 and normal data quality are called valid pixels, and then the mean of each track-crossing position in the track-crossing direction is estimated. Secondly, each scanning line in the data block is processed cyclically, and the total vertical column of SO2 VCD0 is separated from the above-mentioned mean value of the crossing position by difference, so as to obtain the total vertical column of SO2 product after the first correction, namely VCD1.

[0031] Step 303: Based on the SO2 vertical column total product after the first correction and the preset SO2 concentration threshold, pixels exceeding the SO2 concentration threshold are excluded; using a sliding window, the background error matrix corresponding to each scanning line is estimated, and the SO2 vertical column total product after the first correction is difference-separated from the background error to obtain the SO2 vertical column total product after the second correction.

[0032] The embodiment of the present invention can further screen out unqualified pixels that need to be eliminated based on the total amount of the SO2 vertical column VCD1 after the first correction and the preset SO2 concentration threshold. For example, pixels that exceed the SO2 concentration threshold can be excluded. The size of the sliding window is used to estimate the background error matrix corresponding to each scanning line. For example, the mean of the crossing position is used as the background error to obtain the background error matrix for each scanning line. The total amount of the SO2 vertical column VCD1 after the first correction is differentially separated from the background error (such as the mean of the crossing position) to obtain the total amount of the SO2 vertical column VCD2 after the second correction. It should be noted that the differential separation removes the background error from the total amount of the SO2 vertical column. Simply put, it is to subtract these errors contained in the total amount product.

[0033] Step 304: loop iteration is performed based on the corrected SO2 vertical column total amount product, the SO2 concentration threshold and the sliding window until the SO2 vertical column total amount products obtained by two consecutive corrections meet the preset conditions.

[0034] Continue to use the corrected SO2 vertical column total product, the preset concentration threshold and the sliding window to loop through step 303 to obtain each corrected SO2 vertical column total product, such as the three-corrected SO2 vertical column total product VCD3, the four-corrected SO2 vertical column total product VCD4, etc. When the SO2 vertical column total products obtained by two consecutive corrections meet the preset conditions, such as the error (such as the difference) between the three-corrected SO2 vertical column total product VCD3 and the four-corrected SO2 vertical column total product VCD4 is less than the preset error, it is considered that the background error correction is completed. Figure 4 As shown, Figure 4 The schematic diagram of the result obtained by the method provided by the embodiment of the present invention is shown, and it can be clearly seen that Figure 4 There are no obvious seams.

[0035] Compared with the existing background error correction method, the SO2 product background error correction method provided by the embodiment of the present invention can effectively solve the obvious seam phenomenon of trace gas SO2 inversion results at the integration time change point, because the technical solution proposed by the present invention is divided into different data blocks according to the integration time, and two background error schemes are used for iterative processing for each different integration time data block. It is more suitable for orbital observation of loads with discontinuous integration time (such as FY3F OMS), and improves the accuracy and application range of satellite-borne SO2 vertical column total inversion products. At the same time, the present invention has a certain reference value for other trace gas inversion algorithms.

[0036] Optionally, obtaining the total amount product of the vertical column of SO2 for the entire layer includes: using the observation data of the ultraviolet hyperspectral satellite payload irradiance and radiance as input, using a forward model and a classical differential absorption theory inversion model to obtain the total amount product of the vertical column of SO2 for the entire layer.

[0037] Specifically, L1 observation data usually refers to the first-level radiation data product in satellite remote sensing data. These data products contain the original observation data directly obtained from the satellite sensor, which is processed preliminarily to facilitate subsequent quantitative analysis and application. The embodiment of the present invention can use the L1 observation data of ultraviolet hyperspectral satellite payload irradiance and radiance as input, and with the help of the forward model SCIATRAN, obtain the remote sensing product VCD0 of the total vertical column of SO2 in the whole layer through the classical differential absorption theory DOAS (Differential Optical Absorption Spectroscopy, differential optical absorption spectroscopy technology) inversion model. Among them, the classical DOAS algorithm is based on the Beer-Lambert law, using the ratio of the solar spectrum measured at the top of the atmosphere and the backscattered spectrum of the earth, taking the logarithm, and passing through a high-pass filter, and the molecular differential absorption cross section and the Ring effect (referring to the phenomenon caused by the rotational Raman scattering of sunlight by O2 and N2 molecules in the atmosphere) differential pseudo-absorption cross section are least squares fitted to obtain the total oblique column of the target gas. The atmospheric mass factor (AMF) is obtained by calculation using the radiation transfer model software, and the total amount of the oblique column is converted into the total amount of the vertical column.

[0038] Optionally, the above step 302 "divides the total vertical column product of SO2 of each orbit into data blocks with different integration times" may include: based on the observation integration time of each observation pixel, dividing the total vertical column product of SO2 of each orbit into data blocks with different integration times. That is, the observation integration time of each observation pixel may be obtained according to the Integration time data set in the onboard satellite SO2 product, and the total vertical column product VCD0 of SO2 obtained by inverting each orbit of polar-orbiting satellite observation may be divided into data blocks with different integration times.

[0039] Optionally, in the above step 302, "filtering valid pixels for each data block with the same integration time" may include: filtering valid pixels for each data block with the same integration time based on the same load cloud product and quality identification code. Specifically, the valid pixels of all scan lines may be filtered based on the same load cloud product and identification code. The QA quality identification code generally refers to a code used to identify the quality status of a product in the process of quality management and quality control.

[0040] Optionally, in the embodiment of the present invention, the preset SO2 concentration threshold is 2DU (DU = Dobson Units or milli atm cm, 1DU = 2.69 × 10 16 Number of molecules / cm 2). The size of the sliding window changes with the actual spatial resolution of the satellite. For example, taking the spatial resolution of the subsatellite point as 7km×7km, the sliding window can be set to 200 scanning lines. The subsatellite point refers to the subsatellite point of the satellite, specifically the intersection of the line connecting the instantaneous position of the satellite and the center of the earth and the surface of the earth, expressed in geographic longitude and latitude. It can be understood that when the satellite takes pictures at the subsatellite point, the geometric distortion of the image is minimal.

[0041] Optionally, the preset condition includes: a relative deviation between the total vertical column amount products of SO2 obtained from two consecutive calibrations is less than or equal to a preset percentage.

[0042] In the embodiment of the present invention, the preset condition can be set as the relative deviation between the total amount of SO2 vertical column products obtained by two consecutive calibrations is less than or equal to a preset percentage, for example, the relative deviation between the total amount of SO2 vertical column products obtained by two consecutive calibrations is less than or equal to 5%. Wherein, relative deviation = (ab) / b, a is the latter calibration result of the total amount of SO2 vertical column products obtained by two consecutive calibrations, and b is the former calibration result of the total amount of SO2 vertical column products obtained by two consecutive calibrations. When the relative deviation between the total amount of SO2 vertical column products obtained by two consecutive calibrations is less than or equal to the preset percentage, it is considered that the correction of the background error is completed.

[0043] The above describes in detail the SO2 product background error correction method provided by the embodiment of the present invention. The method can also be implemented by a corresponding device. The following describes in detail the SO2 product background error correction device provided by the embodiment of the present invention.

[0044] Figure 5 FIG. 2 shows a schematic diagram of the structure of a SO2 product background error correction device provided by an embodiment of the present invention. Figure 5 As shown, the SO2 product background error correction device includes a processor, which includes: a data acquisition module 61 , a first correction module 62 , a second correction module 63 and a loop iteration module 64 .

[0045] The data acquisition module 61 is used to acquire the total amount product of the vertical column of the entire layer of SO2, wherein the total amount product of the vertical column of the entire layer of SO2 includes the total amount product of the vertical column of SO2 of multiple tracks, and each track includes multiple scanning lines.

[0046] The first correction module 62 is used to divide the SO2 vertical column total product of each track into data blocks with different integration times; for each data block with the same integration time, valid pixels are selected to estimate the mean of each cross-track position in the cross-track direction; each scanning line in the data block with the same integration time is cyclically processed to separate the SO2 vertical column total product from the mean of the cross-track position to obtain the SO2 vertical column total product after the first correction.

[0047] The second correction module 63 is used to exclude pixels exceeding the SO2 concentration threshold based on the SO2 vertical column total amount product after the first correction and the preset SO2 concentration threshold; use a sliding window to estimate the background error matrix corresponding to each of the scanning lines, perform difference separation on the SO2 vertical column total amount product after the first correction and the background error, and obtain the SO2 vertical column total amount product after the second correction.

[0048] The loop iteration module 64 is used to perform loop iteration based on the corrected SO2 vertical column total volume product, the SO2 concentration threshold and the sliding window until the SO2 vertical column total volume products obtained by two consecutive corrections meet a preset condition.

[0049] Optionally, the data acquisition module 61 includes: taking the observation data of ultraviolet hyperspectral satellite payload irradiance and radiance as input, and obtaining the total vertical column product of the entire layer of SO2 with the help of a forward model and a classical differential absorption theory inversion model.

[0050] Optionally, the first correction module 62 includes: dividing the total vertical column product of SO2 of each track into data blocks with different integration times based on the observation integration time of each observation pixel.

[0051] Optionally, the first correction module 62 includes: based on the same load cloud cover product and quality identification code, screening valid pixels for each data block with the same integration time.

[0052] Optionally, the preset SO2 concentration threshold is 2DU; the size of the sliding window changes with the actual spatial resolution of the satellite.

[0053] Optionally, the preset condition includes: a relative deviation between the total vertical column amount products of SO2 obtained from two consecutive calibrations is less than or equal to a preset percentage.

[0054] Compared with the existing background error correction device, the device provided by the embodiment of the present invention can effectively solve the obvious joint phenomenon of the trace gas SO2 inversion results at the change of integration time, because the device proposed by the present invention is divided into different data blocks according to the integration time, and two background error modules are used for iterative processing for each data block with different integration time. It is more suitable for payloads with discontinuous integration time for orbital observation (such as OMS), and improves the accuracy and application range of the satellite-borne SO2 vertical column total amount inversion product.

[0055] It should be noted that the SO2 product background error correction device provided in the above embodiment only uses the division of the above functional modules as an example to implement the corresponding functions. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the SO2 product background error correction device provided in the above embodiment and the SO2 product background error correction method embodiment belong to the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.

[0056] According to one aspect of the present application, an embodiment of the present invention further provides a computer program product, the computer program product comprising a computer program, the computer program comprising a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through a communication part. When the computer program is executed by a processor, the SO2 product background error correction method provided in the embodiment of the present application is executed.

[0057] In addition, an embodiment of the present invention further provides an electronic device, including a bus, a transceiver, a memory, a processor, and a computer program stored in the memory and executable on the processor. The transceiver, the memory, and the processor are respectively connected via a bus. When the computer program is executed by the processor, each process of the above-mentioned SO2 product background error correction method embodiment is implemented, and the same technical effect can be achieved. To avoid repetition, it will not be described here.

[0058] For details, see Figure 6 As shown, the electronic device includes a bus 1110 , a processor 1120 , a transceiver 1130 , a bus interface 1140 , a memory 1150 , and a user interface 1160 .

[0059] In the embodiment of the present invention, the electronic device further includes: a computer program stored in the memory 1150 and executable on the processor 1120, and when the computer program is executed by the processor 1120, each process of the above-mentioned SO2 product background error correction method embodiment is implemented.

[0060] The transceiver 1130 is configured to receive and send data under the control of the processor 1120 .

[0061] In an embodiment of the present invention, the bus architecture (represented by bus 1110), bus 1110 may include any number of interconnected buses and bridges, and bus 1110 connects various circuits including one or more processors represented by processor 1120 and a memory represented by memory 1150.

[0062] Bus 1110 represents one or more of any of several types of bus structures, including a memory bus and memory controller, a peripheral bus, an Accelerate Graphical Port (AGP), a processor, or a local bus using any of a variety of bus architectures. By way of example and not limitation, such architectures include: Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA), Peripheral Component Interconnect (PCI) bus.

[0063] The processor 1120 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method embodiment may be completed by an integrated logic circuit of hardware in the processor or an instruction in the form of software. The above processors include: a general-purpose processor, a central processing unit (CPU), a network processor (NP), a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA), a complex programmable logic device (CPLD), a programmable logic array (PLA), a microcontroller unit (MCU) or other programmable logic devices, discrete gates, transistor logic devices, discrete hardware components. The methods, steps and logic block diagrams disclosed in the embodiments of the present invention may be implemented or executed. For example, the processor may be a single-core processor or a multi-core processor, and the processor may be integrated into a single chip or located in multiple different chips.

[0064] Processor 1120 can be a microprocessor or any conventional processor. The method steps disclosed in conjunction with the embodiments of the present invention can be directly executed by a hardware decoding processor, or by a combination of hardware and software modules in the decoding processor. The software module can be located in a readable storage medium known in the art, such as a random access memory (RAM), a flash memory (FlashMemory), a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), a register, etc. The readable storage medium is located in a memory, and the processor reads the information in the memory and completes the steps of the above method in conjunction with its hardware.

[0065] The bus 1110 may also connect various other circuits such as peripheral devices, voltage regulators or power management circuits, and the bus interface 1140 provides an interface between the bus 1110 and the transceiver 1130, which are well known in the art. Therefore, the embodiments of the present invention will not be further described.

[0066] The transceiver 1130 may be one element or multiple elements, such as multiple receivers and transmitters, and provides a unit for communicating with various other devices on a transmission medium. For example, the transceiver 1130 receives external data from other devices, and the transceiver 1130 is used to send data processed by the processor 1120 to other devices. Depending on the nature of the computer system, a user interface 1160 may also be provided, such as a touch screen, a physical keyboard, a display, a mouse, a speaker, a microphone, a trackball, a joystick, and a stylus.

[0067] It should be understood that in an embodiment of the present invention, the memory 1150 may further include a memory remotely arranged relative to the processor 1120, and these remotely arranged memories may be connected to the server through a network. One or more parts of the above-mentioned network may be an ad hoc network, an intranet, an extranet, a virtual private network (VPN), a local area network (LAN), a wireless local area network (WLAN), a wide area network (WAN), a wireless wide area network (WWAN), a metropolitan area network (MAN), the Internet, a public switched telephone network (PSTN), a plain old telephone service network (POTS), a cellular telephone network, a wireless network, a wireless fidelity (Wi-Fi) network, and a combination of two or more of the above-mentioned networks. For example, the cellular telephone network and the wireless network can be a Global System for Mobile Communications (GSM) system, a Code Division Multiple Access (CDMA) system, a Worldwide Interoperability for Microwave Access (WiMAX) system, a General Packet Radio Service (GPRS) system, a Wideband Code Division Multiple Access (WCDMA) system, a Long Term Evolution (LTE) system, a LTE Frequency Division Duplex (FDD) system, a LTE Time Division Duplex (TDD) system, an Advanced Long Term Evolution (LTE-A) system, a Universal Mobile Telecommunications (UMTS) system, an Enhanced Mobile Broadband (Enhance Mobile Broadband, eMBB) system, a massive Machine Type of Communication (mMTC) system, an Ultra-Reliable Low Latency Communications (UltraReliable Low Latency Communications, uRLLC) system, etc.

[0068] It should be understood that the memory 1150 in the embodiment of the present invention may be a volatile memory or a non-volatile memory, or may include both a volatile memory and a non-volatile memory. Among them, the non-volatile memory includes: a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory.

[0069] Volatile memory includes: Random Access Memory (RAM), which is used as an external cache. By way of example but not limitation, many forms of RAM are available, such as: Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate Synchronous Dynamic Random Access Memory (DDRSDRAM), Enhanced Synchronous Dynamic Random Access Memory (ESDRAM), Synchronous Link Dynamic Random Access Memory (SLDRAM) and Direct Rambus RAM (DRRAM). The memory 1150 of the electronic device described in the embodiment of the present invention includes but is not limited to the above and any other suitable types of memory.

[0070] In the embodiment of the present invention, the memory 1150 stores the following elements of the operating system 1151 and the application program 1152: executable modules, data structures, or subsets thereof, or extended sets thereof.

[0071] Specifically, the operating system 1151 includes various system programs, such as a framework layer, a core library layer, a driver layer, etc., which are used to implement various basic services and process hardware-based tasks. The application 1152 includes various applications, such as a media player (Media Player) and a browser (Browser), which are used to implement various application services. The program for implementing the method of the embodiment of the present invention may be included in the application 1152. The application 1152 includes applets, objects, components, logic, data structures, and other computer system executable instructions that perform specific tasks or implement specific abstract data types.

[0072] In addition, an embodiment of the present invention further provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the various processes of the above-mentioned SO2 product background error correction method embodiment are implemented, and the same technical effect can be achieved. To avoid repetition, it will not be repeated here.

[0073] Computer readable storage media include: permanent and non-permanent, removable and non-removable media, which are tangible devices that can retain and store instructions for use by instruction execution devices. Computer readable storage media include: electronic storage devices, magnetic storage devices, optical storage devices, electromagnetic storage devices, semiconductor storage devices, and any suitable combination of the above. Computer readable storage media include: phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), non-volatile random access memory (NVRAM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassette storage, magnetic tape disk storage or other magnetic storage devices, memory sticks, mechanical encoding devices (such as punched cards or raised structures in grooves with instructions recorded thereon) or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined in the embodiments of the present invention, computer-readable storage media do not include temporary signals themselves, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through waveguides or other transmission media (such as light pulses passing through fiber optic cables), or electrical signals transmitted through wires.

[0074] In the several embodiments provided in the present application, it should be understood that the disclosed devices, electronic devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the modules or units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or it can be an electrical, mechanical or other form of connection.

[0075] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, and may be located in one location or distributed on multiple network units. Some or all of the units may be selected according to actual needs to solve the problem to be solved by the embodiments of the present invention.

[0076] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0077] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiment of the present invention is essentially or part of the contribution to the prior art, or all or part of the technical solution can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (including: a personal computer, a server, a data center or other network device) to perform all or part of the steps of the method described in each embodiment of the present invention. The above-mentioned storage medium includes various media that can store program codes as listed above.

[0078] In the description of the embodiments of the present invention, those skilled in the art should know that the embodiments of the present invention can be implemented as methods, devices, electronic devices and computer-readable storage media. Therefore, the embodiments of the present invention can be specifically implemented in the following forms: complete hardware, complete software (including firmware, resident software, microcode, etc.), a combination of hardware and software. In addition, in some embodiments, the embodiments of the present invention can also be implemented in the form of a computer program product in one or more computer-readable storage media, and the computer-readable storage medium contains computer program code.

[0079] The above-mentioned computer-readable storage medium may adopt any combination of one or more computer-readable storage media. Computer-readable storage media include: electrical, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or devices, or any combination of the above. More specific examples of computer-readable storage media include: portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM), flash memories, optical fibers, compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices or any combination of the above. In an embodiment of the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program, which may be used by or in combination with an instruction execution system, device, or device.

[0080] The computer program code contained in the above-mentioned computer-readable storage medium can be transmitted using any appropriate medium, including: wireless, wire, optical cable, radio frequency (RF) or any suitable combination thereof.

[0081] The computer program code for performing the operation of the embodiments of the present invention can be written in assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, state setting data, integrated circuit configuration data, or in one or more programming languages ​​or a combination thereof, wherein the programming language includes an object-oriented programming language, such as Java, Smalltalk, C++, and also includes a conventional procedural programming language, such as C language or a similar programming language. The computer program code can be executed entirely on the user's computer, partially on the user's computer, as an independent software package, partially on the user's computer, partially on a remote computer, and completely on a remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer or to an external computer through any type of network, including a local area network (LAN) or a wide area network (WAN).

[0082] The embodiments of the present invention describe the provided methods, devices, and electronic devices through flowcharts and / or block diagrams.

[0083] It should be understood that each box in the flowchart and / or block diagram and the combination of boxes in the flowchart and / or block diagram can be implemented by computer-readable program instructions. These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer or other programmable data processing device to produce a machine, and these computer-readable program instructions are executed by a computer or other programmable data processing device to produce a device that implements the functions / operations specified by the boxes in the flowchart and / or block diagram.

[0084] These computer-readable program instructions may also be stored in a computer-readable storage medium that enables a computer or other programmable data processing device to work in a specific manner. In this way, the instructions stored in the computer-readable storage medium produce an instruction device product including functions / operations specified in the blocks in the flowchart and / or block diagram.

[0085] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device, so that a series of operational steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby enabling the instructions executed on the computer or other programmable data processing apparatus to provide a process for implementing the functions / operations specified in the blocks in the flowchart and / or block diagram.

[0086] The above is only a specific implementation of the embodiment of the present invention, but the protection scope of the embodiment of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the embodiment of the present invention, which should be included in the protection scope of the embodiment of the present invention. Therefore, the protection scope of the embodiment of the present invention should be based on the protection scope of the claims.

Claims

1. A method for correcting background error of SO2 product, characterized in that: include: Obtaining a total amount product of a vertical column of SO2 for a whole layer, wherein the total amount product of a vertical column of SO2 for a whole layer includes total amount products of vertical columns of SO2 for a plurality of tracks, and each track includes a plurality of scanning lines; The total amount product of the vertical column of SO2 of each track is divided into data blocks with different integration times; for each data block with the same integration time, effective pixels are selected to estimate the mean value of each track crossing position in the track crossing direction; each scanning line in the data block with the same integration time is processed cyclically, and the total amount product of the vertical column of SO2 is separated from the mean value of the track crossing position to obtain the total amount product of the vertical column of SO2 after the first correction; Based on the first corrected SO2 vertical column total product and the preset SO2 concentration threshold, pixels exceeding the SO2 concentration threshold are excluded; the background error matrix corresponding to each of the scanning lines is estimated using a sliding window, and the first corrected SO2 vertical column total product and the background error are difference-separated to obtain a second corrected SO2 vertical column total product; A loop iteration is performed based on the corrected total SO2 vertical column product, the SO2 concentration threshold and the sliding window until the total SO2 vertical column products obtained by two consecutive corrections meet a preset condition.

2. The method according to claim 1, characterized in that The method of obtaining the total amount of the vertical column of SO2 of the entire layer includes: The observation data of ultraviolet hyperspectral satellite payload irradiance and radiance are used as input, and the total vertical column product of the entire layer of SO2 is obtained with the help of a forward model and through a classical differential absorption theory inversion model.

3. The method according to claim 1, characterized in that The step of dividing the total amount of SO2 vertical column of each track into data blocks with different integration times includes: dividing the total amount of SO2 vertical column of each track into data blocks with different integration times based on the observation integration time of each observation pixel.

4. The method according to claim 1, characterized in that: The screening of valid pixels for each data block with the same integration time includes: screening of valid pixels for each data block with the same integration time based on the same load cloud product and quality identification code.

5. The method according to claim 1, characterized in that The preset SO2 concentration threshold is 2DU; the size of the sliding window changes with the actual spatial resolution of the satellite.

6. The method according to claim 1, characterized in that The preset condition includes: the relative deviation between the total amount of SO2 vertical column products obtained from two consecutive calibrations is less than or equal to a preset percentage.

7. A SO2 product background error correction device, characterized in that: include: A data acquisition module, a first correction module, a second correction module and a loop iteration module; The data acquisition module is used to acquire the total amount product of the entire layer of SO2 vertical column, and the total amount product of the entire layer of SO2 vertical column includes the total amount product of the SO2 vertical column of multiple tracks, and each track includes multiple scanning lines; The first correction module is used to divide the total amount of SO2 vertical column products of each track into data blocks with different integration times; for each data block with the same integration time, effective pixels are selected to estimate the mean value of each track crossing position in the track crossing direction; each scanning line in the data block with the same integration time is cyclically processed to separate the total amount of SO2 vertical column products from the mean value of the track crossing position to obtain the total amount of SO2 vertical column products after the first correction; The second correction module is used to exclude pixels exceeding the SO2 concentration threshold based on the first corrected SO2 vertical column total product and the preset SO2 concentration threshold; use a sliding window to estimate the background error matrix corresponding to each of the scanning lines, perform difference separation on the first corrected SO2 vertical column total product and the background error, and obtain a second corrected SO2 vertical column total product; The loop iteration module is used to perform loop iteration based on the corrected SO2 vertical column total volume product, the SO2 concentration threshold and the sliding window until the SO2 vertical column total volume products obtained by two consecutive corrections meet a preset condition.

8. An electronic device comprising a processor and a memory, wherein the memory stores a computer program, characterized in that: The processor executes the computer program stored in the memory to implement the steps in the SO2 product background error correction method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps in the SO2 product background error correction method according to any one of claims 1 to 6 are implemented.

10. A computer program product, characterized in that The invention comprises a computer program, and when the computer program is executed, the steps in the SO2 product background error correction method according to any one of claims 1 to 6 can be implemented.

Images