Cross-radiation calibration method, device, equipment and medium for spaceborne hyperspectral remote sensor

Abstract

The present invention relates to the field of remote sensing technology, and in particular to a method, device, equipment and medium for cross-radiation calibration of a satellite-borne hyperspectral remote sensor. The present invention reconstructs the apparent reflectance continuous spectrum of a hyperspectral remote sensor to be calibrated by referring to multispectral remote sensor data and a surface object apparent reflectance model, thereby solving the problem that in the current on-orbit cross-radiation calibration process of satellite-borne hyperspectral remote sensors, there is a lack of stable and reliable hyperspectral remote sensors as references and high-precision satellite-borne multispectral reference devices cannot be directly applied to the on-orbit cross-radiation calibration of hyperspectral remote sensors. The multi-channel data of the reference multispectral remote sensor is fully utilized to improve the calibration accuracy of the hyperspectral remote sensor, and the cross-radiation calibration of the hyperspectral remote sensor based on the reference multispectral remote sensor is realized, thereby improving the frequency and efficiency of the on-orbit absolute radiation calibration of the hyperspectral remote sensor under the condition of ensuring the calibration accuracy.

Description

Technical Field

[0001] The present invention relates to the field of remote sensing technology, and in particular to a method, device, equipment and medium for cross-radiation calibration of a satellite-borne hyperspectral remote sensor. Background Art

[0002] With the rapid development of hyperspectral remote sensing technology, its important application potential and value in industries such as agricultural and forestry classification, land and resources survey, ecological environment protection, and ground object identification are becoming increasingly apparent.

[0003] In the prior art, hyperspectral remote sensing data is usually obtained based on satellite-borne remote sensors. Therefore, absolute radiation calibration of satellite-borne hyperspectral remote sensors in orbit is performed to obtain absolute radiation calibration parameters that can accurately characterize the radiation performance of satellite-borne hyperspectral remote sensors, which is the basis for quantitative research and application of hyperspectral remote sensing data. However, after the satellite-borne remote sensor is launched into space, due to the lack of relatively stable and accurate reference hyperspectral remote sensors, parameters such as surface reflection and atmospheric characteristics data of stable targets can only be obtained through ground synchronous measurement, and the hyperspectral remote sensor is calibrated on-orbit by field substitution combined with radiation transmission calculation. This ground synchronous measurement field substitution calibration method requires manual or automated measurement equipment to measure synchronously on the ground, which is affected by weather, manpower, financial resources and other factors, limiting the calibration frequency and efficiency. Therefore, how to calibrate satellite-borne hyperspectral remote sensors is an important issue that needs to be solved urgently. Summary of the invention

[0004] The present invention provides a method, device, equipment and medium for cross-radiation calibration of a satellite-borne hyperspectral remote sensor, which are used to solve the problem of lack of high-precision satellite-borne hyperspectral remote sensor as a reference during the on-orbit cross-radiation calibration of the satellite-borne hyperspectral remote sensor, so as to improve the frequency and efficiency of the on-orbit absolute radiation calibration of the hyperspectral remote sensor while ensuring the calibration accuracy.

[0005] The present invention provides a cross-radiation calibration method for a spaceborne hyperspectral remote sensor, comprising the following steps:

[0006] Obtain reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated;

[0007] According to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated, the correction factors of the reference multispectral remote sensor and the target hyperspectral remote sensor are obtained;

[0008] Establishing the surface reflectivity model of the ground objects;

[0009] Based on the reference multispectral remote sensor data and correction factors, the surface reflectance model of the ground object is solved to obtain the hyperspectral reconstruction parameters;

[0010] Reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0011] The on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined based on the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0012] According to the cross-radiation calibration method for spaceborne hyperspectral remote sensor provided by the present invention, the correction factor of the reference multispectral remote sensor and the target hyperspectral remote sensor is obtained according to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated, and the expression is:

[0013]

[0014] in, is the apparent reflectance of the i-th band of the reference multispectral remote sensor, and its expression is:

[0015]

[0016] in, is the spectral response function of the i-th band of the reference multispectral remote sensor, R(λ) is the apparent reflectance continuous spectrum of the calibration site; is the starting wavelength of the i-th band of the reference multispectral remote sensor, is the end wavelength of the i-th band of the reference multispectral remote sensor;

[0017] is the apparent reflectance of the hyperspectral remote sensor to be calibrated corresponding to the i-th band of the reference multispectral remote sensor, and its expression is:

[0018]

[0019] Among them, R H (λ) is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated to the calibration site, which is obtained by merging the apparent reflectance continuous spectra of the hyperspectral remote sensor to be calibrated to the calibration site in multiple sub-bands. The expression is:

[0020]

[0021] in, is the starting wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the end wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the spectral response function of the jth band of the hyperspectral remote sensor to be calibrated.

[0022] According to the cross-radiation calibration method for spaceborne hyperspectral remote sensors provided by the present invention, a surface reflectivity model of ground objects is established, and its expression is:

[0023] R I =C w P I,w +C v P I,v +C s P I,s +C 4 P I,4 +r I

[0024] Among them, R I is the apparent reflectivity of the ground object in band I, P I,w is the formatted apparent reflectance spectrum of the water standard model, P I,v is the formatted apparent reflectance spectrum of the vegetation standard model, P I,s is the formatted apparent reflectance spectrum of the bare soil standard model, P I,4 The formatted apparent reflectance spectrum of the additional standard mode, C w is the water body standard model coefficient, C v is the vegetation standard model coefficient, C s is the bare soil standard model coefficient, C 4 is the additional standard mode coefficient, r I is the residual.

[0025] According to the cross-radiation calibration method for a satellite-borne hyperspectral remote sensor provided by the present invention, a surface reflectance model of a ground object is solved based on reference multispectral remote sensor data and a correction factor to obtain hyperspectral reconstruction parameters, comprising the following steps:

[0026] According to the surface reflectance model of the ground object, the apparent reflectance model of the hyperspectral remote sensor channel corresponding to the reference multispectral remote sensor band is obtained, and its expression is:

[0027] R H-M =F·(P M C+r)

[0028] Where F is the correction factor f i The obtained correction factor matrix, C is the coefficient matrix, and r is the residual matrix;

[0029] P M is the reference multispectral remote sensor channel apparent reflectance matrix, and the expression of its matrix elements is:

[0030]

[0031] Among them, P′ k (λ) is the hyperspectral apparent reflectance of the kth standard mode;

[0032] The reference multispectral remote sensor data is used to solve the hyperspectral remote sensor channel apparent reflectance model corresponding to the reference multispectral remote sensor band, and the hyperspectral reconstruction parameters are obtained, which are expressed as follows:

[0033]

[0034] Where F' is the correction factor f i The inverse of the correction factor matrix is ​​obtained.

[0035] According to the cross-radiation calibration method for spaceborne hyperspectral remote sensors provided by the present invention, the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated is reconstructed according to the hyperspectral reconstruction parameters, and its expression is:

[0036] R H =P H C+r H

[0037] Among them, r H is the reconstructed residual matrix;

[0038] P H is the standard mode apparent reflectance spectrum matrix, and the expression of its matrix elements is:

[0039]

[0040] According to the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated, the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined, and its expression is:

[0041]

[0042] Among them, D H is the quantized bit data obtained from the hyperspectral remote sensor data to be calibrated, R H is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0043] The present invention also provides a cross-radiation calibration device for a spaceborne hyperspectral remote sensor, comprising:

[0044] A data acquisition module, used to acquire reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated;

[0045] A correction factor module is used to obtain the correction factor of the reference multispectral remote sensor and the target hyperspectral remote sensor according to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated;

[0046] Model building module, used to build the surface reflectivity model of the ground object;

[0047] A model solving module is used to solve the surface reflectance model of the ground object according to the reference multispectral remote sensor data and the correction factor to obtain the hyperspectral reconstruction parameters;

[0048] A spectrum reconstruction module, used to reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0049] The calibration module is used to determine the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated based on the hyperspectral remote sensor data to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0050] The present invention also provides 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 program, the cross-radiation calibration method for a satellite-borne hyperspectral remote sensor as described above is implemented.

[0051] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the cross-radiation calibration method for a space-borne hyperspectral remote sensor as described in any one of the above is implemented.

[0052] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the cross-radiation calibration method for a spaceborne hyperspectral remote sensor as described above is implemented.

[0053] The satellite-borne hyperspectral remote sensor cross-radiation calibration method, device, equipment and medium provided by the present invention reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated by referring to multispectral remote sensor data and a ground object apparent reflectance model, thereby solving the problem that in the current on-orbit cross-radiation calibration process of satellite-borne hyperspectral remote sensor, there is a lack of stable and reliable hyperspectral remote sensor as a reference and the high-precision satellite-borne multispectral reference device cannot be directly applied to the on-orbit cross-radiation calibration of the hyperspectral remote sensor. The multi-channel data of the reference multispectral remote sensor is fully utilized to improve the calibration accuracy of the hyperspectral remote sensor, and the cross-radiation calibration of the hyperspectral remote sensor based on the reference multispectral remote sensor is realized, thereby improving the frequency and efficiency of the on-orbit absolute radiation calibration of the hyperspectral remote sensor while ensuring the calibration accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces 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 creative work.

[0055] Figure 1It is a schematic diagram of the flow of the cross-radiation calibration method of a space-borne hyperspectral remote sensor provided by the present invention;

[0056] Figure 2 It is a structural schematic diagram of a cross-radiation calibration device for a satellite-borne hyperspectral remote sensor provided by the present invention;

[0057] Figure 3 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

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

[0059] Combine the following Figure 1-Figure 3 The invention describes the cross-radiation calibration method, device, equipment and medium of a spaceborne hyperspectral remote sensor.

[0060] Figure 1 FIG. 1 is a flow chart of a cross-radiation calibration method for a spaceborne hyperspectral remote sensor according to an embodiment of the present invention. Figure 1 As shown, it includes steps 110 to 160, specifically:

[0061] Step 110, obtaining reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated;

[0062] Step 120, obtaining a calibration factor of the reference multispectral remote sensor and the target hyperspectral remote sensor according to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated;

[0063] Step 130, establishing a surface reflectivity model of the object;

[0064] Step 140, solving the surface reflectance model of the ground object according to the reference multispectral remote sensor data and the correction factor to obtain the hyperspectral reconstruction parameters;

[0065] Step 150, reconstructing the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0066] Step 160: Determine the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated based on the hyperspectral remote sensor data to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0067] In step 110 of the embodiment of the present invention, the reference multispectral remote sensor data is data obtained according to the reference multispectral remote sensor, and the reference multispectral remote sensor may be a currently internationally recognized high-precision satellite-borne multispectral remote sensor, such as MODIS, S-NPP / VIIRS, Landsat-8 / OLI, Sentinel-2 / MSI, etc.; the embodiment of the present invention may select the best multispectral remote sensor as a reference benchmark for cross-radiation calibration of the satellite-borne hyperspectral remote sensor from among multiple satellite-borne multispectral remote sensors based on the analysis and evaluation of the on-orbit spectral radiation performance; or may select multiple multispectral remote sensor data for joint calibration.

[0068] In step 110 of the embodiment of the present invention, the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated are remote sensing data acquired by the multispectral remote sensor and the hyperspectral remote sensor for the same calibration site. It can be understood that the calibration site should be selected as a large area of ​​uniform site, such as the Gobi land of the Dunhuang calibration site, the large area of ​​desert in Golmud, Qinghai, etc., to improve the calibration effect. Due to the difference in spatial resolution between the reference payload and the target payload, spatial scale normalization is required, usually by scaling the high-resolution image to a coarse-resolution remote sensing image, and then aligning the reference and target images to minimize the calibration error caused by the spatial difference.

[0069] In step 120 of the embodiment of the present invention, the correction factor is used to describe the influence of factors such as the sun-observation angle difference caused by the spectral difference and imaging time difference between the reference multispectral remote sensor and the target hyperspectral remote sensor, and the atmosphere. The correction factor of the reference multispectral remote sensor and the target hyperspectral remote sensor is obtained according to the reference multispectral remote sensor data and the target hyperspectral remote sensor data, and its expression is:

[0070]

[0071] in, is the apparent reflectance of the i-th band of the reference multispectral remote sensor, and its expression is:

[0072]

[0073] in, is the spectral response function of the i-th band of the reference multispectral remote sensor, R(λ) is the apparent reflectance continuous spectrum of the calibration site; is the starting wavelength of the i-th band of the reference multispectral remote sensor, is the end wavelength of the i-th band of the reference multispectral remote sensor;

[0074] is the apparent reflectance of the hyperspectral remote sensor to be calibrated corresponding to the i-th band of the reference multispectral remote sensor, and its expression is:

[0075]

[0076] Among them, R H (λ) is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated to the calibration site, which is obtained by merging the apparent reflectance continuous spectra of the hyperspectral remote sensor to be calibrated to the calibration site in multiple sub-bands. The expression is:

[0077]

[0078] in, is the starting wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the end wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the spectral response function of the jth band of the hyperspectral remote sensor to be calibrated.

[0079] In the embodiment of the present invention, the apparent reflectance continuous spectrum R(λ) of the calibration site can be obtained by simulating the atmospheric radiation transfer model based on the historical ground measured reflectance spectrum of the calibration site.

[0080] In the embodiment of the present invention, the spectral shape of any pixel extracted from the remote sensing image acquired by the remote sensor can be a linear combination of the apparent reflectance of three standard modes (water body, vegetation, bare soil) and an optional additional mode. Therefore, in step 130, a surface reflectance model of the surface object is established, and its expression is:

[0081] R I =C w P I,w +C v P I,v +C s P I,s +C 4 P I,4 +r I

[0082] Among them, R I is the apparent reflectivity of the ground object in band I, P I,w is the formatted apparent reflectance spectrum of the water standard model, P I,v is the formatted apparent reflectance spectrum of the vegetation standard model, P I,s is the formatted apparent reflectance spectrum of the bare soil standard model, P I,4 The formatted apparent reflectance spectrum of the additional standard mode, C w is the water body standard model coefficient, C v is the vegetation standard model coefficient, C s is the bare soil standard model coefficient, C 4is the additional standard mode coefficient, r I is the residual.

[0083] In step 140 of the embodiment of the present invention, solving the surface reflectance model of the ground object according to the reference multispectral remote sensor data and the correction factor to obtain the hyperspectral reconstruction parameters includes the following steps:

[0084] According to the surface reflectance model of the ground object, the surface reflectance model of the hyperspectral remote sensor channel corresponding to the reference multispectral remote sensor band is obtained;

[0085] For the reference multispectral remote sensor, based on the surface reflectance model of the ground object, the apparent reflectance matrix R of the calibration field actually observed by the reference multispectral remote sensor is M It can be expressed as:

[0086] R M =P M C+r

[0087] Based on R M and correction factor f i , we can get the hyperspectral remote sensor channel apparent reflectance model corresponding to the reference multispectral remote sensor band, and its expression is:

[0088] R H-M =F·(P M C+r)

[0089] Where F is the correction factor f i The obtained correction factor matrix, C is the coefficient matrix, and r is the residual matrix;

[0090] P M is the reference multispectral remote sensor channel apparent reflectance matrix, and the expression of its matrix elements is:

[0091]

[0092] Among them, P′ k (λ) is the hyperspectral apparent reflectance of the kth standard mode; in the embodiment of the present invention, k has three values, representing the three standard modes of water, vegetation and bare soil, respectively, and the hyperspectral apparent reflectance is 0.1cm -1 Hyperspectral apparent reflectance can be simulated by inputting the surface reflectance spectra of different standard models into the atmospheric radiation transfer model;

[0093] When the number of channels of the reference multispectral reference remote sensor is greater than or equal to the sum of the number of standard modes and additional modes, the reference multispectral remote sensor data can be used to solve the hyperspectral remote sensor channel apparent reflectance model corresponding to the reference multispectral remote sensor band to obtain the hyperspectral reconstruction parameters, which are expressed as follows:

[0094]

[0095] Where F' is the correction factor f i The inverse of the correction factor matrix is ​​obtained.

[0096] It can be understood that as the number of reference multispectral reference remote sensor channels increases, the degree of fitting of the apparent reflectance model of the hyperspectral remote sensor channels corresponding to the reference multispectral remote sensor bands can be increased, thereby improving the calibration accuracy.

[0097] In step 150 of the embodiment of the present invention, the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated is reconstructed according to the hyperspectral reconstruction parameters, and the expression thereof is:

[0098] R H =P H C+r H

[0099] Among them, r H is the reconstructed residual matrix;

[0100] P H is the standard mode apparent reflectance spectrum matrix, and the expression of its matrix elements is:

[0101]

[0102] In step 160 of the embodiment of the present invention, the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined according to the hyperspectral remote sensor data to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated, and the expression is:

[0103]

[0104] Among them, D H is the quantized bit data obtained from the hyperspectral remote sensor data to be calibrated, R H is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0105] To sum up, the cross-radiation calibration method of a satellite-borne hyperspectral remote sensor in an embodiment of the present invention reconstructs the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated by referring to multispectral remote sensor data and a ground object apparent reflectance model, and makes full use of the multi-channel data of the reference multispectral remote sensor to improve the calibration accuracy of the hyperspectral remote sensor. It solves the problem that in the current on-orbit cross-radiation calibration process of satellite-borne hyperspectral remote sensors, there is a lack of stable and reliable hyperspectral remote sensor as a reference and the high-precision satellite-borne multispectral reference device cannot be directly applied to the on-orbit cross-radiation calibration of hyperspectral remote sensors. It realizes the cross-radiation calibration of hyperspectral remote sensors based on the reference multispectral remote sensor, thereby improving the frequency and efficiency of on-orbit absolute radiation calibration of hyperspectral remote sensors while ensuring the calibration accuracy.

[0106] The satellite-borne hyperspectral remote sensor cross-radiation calibration device provided by the present invention is described below. The satellite-borne hyperspectral remote sensor cross-radiation calibration device described below and the satellite-borne hyperspectral remote sensor cross-radiation calibration method described above can correspond to each other.

[0107] Figure 2 is a schematic diagram of the structure of a cross-radiation calibration device for a spaceborne hyperspectral remote sensor according to an embodiment of the present invention. Figure 2 As shown, including:

[0108] The data acquisition module 210 is used to acquire reference multispectral remote sensor data and to-be-calibrated hyperspectral remote sensor data;

[0109] The correction factor module 220 is used to obtain the correction factor of the reference multispectral remote sensor and the target hyperspectral remote sensor according to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated;

[0110] Model building module 230, used to build a surface reflectivity model of the object;

[0111] The model solving module 240 is used to solve the surface reflectance model of the ground object according to the reference multi-spectral remote sensor data and the correction factor to obtain the hyperspectral reconstruction parameters;

[0112] A spectrum reconstruction module 250, for reconstructing a continuous spectrum of apparent reflectance of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0113] The calibration module 260 is used to determine the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated according to the hyperspectral remote sensor data to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0114] It can be seen that the satellite-borne hyperspectral remote sensor cross-radiation calibration device of the embodiment of the present invention can reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated by referring to multispectral remote sensor data and the apparent reflectance model of the ground object, and makes full use of the multi-channel data of the reference multispectral remote sensor to improve the calibration accuracy of the hyperspectral remote sensor. It solves the problem that in the current on-orbit cross-radiation calibration process of satellite-borne hyperspectral remote sensors, there is a lack of stable and reliable hyperspectral remote sensor as a reference and the high-precision satellite-borne multispectral reference device cannot be directly applied to the on-orbit cross-radiation calibration of hyperspectral remote sensors. The cross-radiation calibration of the hyperspectral remote sensor based on the reference multispectral remote sensor is realized, so that the frequency and efficiency of the on-orbit absolute radiation calibration of the hyperspectral remote sensor can be improved while ensuring the calibration accuracy.

[0115] 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 communication interface 320, a memory 330 and a communication bus 340, wherein the processor 310, the communication interface 320 and the memory 330 communicate with each other through the communication bus 340. The processor 310 may call the logic instructions in the memory 330 to execute the cross-radiation calibration method of the spaceborne hyperspectral remote sensor, the method comprising: obtaining reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated;

[0116] According to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated, the correction factors of the reference multispectral remote sensor and the target hyperspectral remote sensor are obtained;

[0117] Establishing the surface reflectivity model of the ground objects;

[0118] Based on the reference multispectral remote sensor data and correction factors, the surface reflectance model of the ground object is solved to obtain the hyperspectral reconstruction parameters;

[0119] Reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0120] The on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined based on the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0121] 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 and other media that can store program codes.

[0122] On the other hand, the present invention further 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, when the computer program is executed by a processor, the computer can execute the spaceborne hyperspectral remote sensor cross-radiation calibration method provided by the above methods, including the following steps:

[0123] Obtain reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated;

[0124] According to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated, the correction factors of the reference multispectral remote sensor and the target hyperspectral remote sensor are obtained;

[0125] Establishing the surface reflectivity model of the ground objects;

[0126] Based on the reference multispectral remote sensor data and correction factors, the surface reflectance model of the ground object is solved to obtain the hyperspectral reconstruction parameters;

[0127] Reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0128] The on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined based on the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

[0129] 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 perform the cross-radiation calibration method for a spaceborne hyperspectral remote sensor provided by the above methods, comprising the following steps:

[0130] Obtain reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated;

[0131] According to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated, the correction factors of the reference multispectral remote sensor and the target hyperspectral remote sensor are obtained;

[0132] Establishing the surface reflectivity model of the ground objects;

[0133] Based on the reference multispectral remote sensor data and correction factors, the surface reflectance model of the ground object is solved to obtain the hyperspectral reconstruction parameters;

[0134] Reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters;

[0135] The on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined based on the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

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

[0137] 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, it 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.

[0138] 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 cross-radiation calibration method for a spaceborne hyperspectral remote sensor, characterized in that: The following steps are involved: Obtain reference multispectral remote sensor data and hyperspectral remote sensor data to be calibrated; According to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated, the correction factors of the reference multispectral remote sensor and the target hyperspectral remote sensor are obtained; Establishing the surface reflectivity model of the ground objects; Based on the reference multispectral remote sensor data and correction factors, the surface reflectance model of the ground object is solved to obtain the hyperspectral reconstruction parameters; Reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters; Determine the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated based on the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated; The correction factor is expressed as: in, is the apparent reflectance of the i-th band of the reference multispectral remote sensor, and its expression is: in, is the spectral response function of the i-th band of the reference multispectral remote sensor, R(λ) is the apparent reflectance continuous spectrum of the calibration site; is the starting wavelength of the i-th band of the reference multispectral remote sensor, is the end wavelength of the i-th band of the reference multispectral remote sensor; is the apparent reflectance of the hyperspectral remote sensor to be calibrated corresponding to the i-th band of the reference multispectral remote sensor, and its expression is: Among them, R H (λ) is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated to the calibration site, which is obtained by merging the apparent reflectance continuous spectra of the hyperspectral remote sensor to be calibrated to the calibration site in multiple sub-bands. The expression is: in, is the starting wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the end wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the spectral response function of the jth band of the hyperspectral remote sensor to be calibrated; The surface reflectivity model of the ground object is expressed as follows: R I =C w P I，w +C v P I，v +C s P I，s +C4P I，4 +r I Among them, R I is the apparent reflectivity of the ground object in band I, P I，w is the formatted apparent reflectance spectrum of the water standard model, P I，v is the formatted apparent reflectance spectrum of the vegetation standard model, P I，s is the formatted apparent reflectance spectrum of the bare soil standard model, P l，4 The formatted apparent reflectance spectrum of the additional standard mode, C w is the water body standard model coefficient, C v is the vegetation standard model coefficient, C s is the bare soil standard model coefficient, C4 is the additional standard model coefficient, r I is the residual; According to the reference multispectral remote sensor data and the correction factor, the surface reflectance model of the ground object is solved to obtain the hyperspectral reconstruction parameters, including the following steps: According to the surface reflectance model of the ground object, the apparent reflectance model of the hyperspectral remote sensor channel corresponding to the reference multispectral remote sensor band is obtained, and its expression is: R H-M =F·(P M C+r) Where F is the correction factor f i The obtained correction factor matrix, C is the coefficient matrix, and r is the residual matrix; P M is the reference multispectral remote sensor channel apparent reflectance matrix, and the expression of its matrix elements is: Among them, P′ k (λ) is the hyperspectral apparent reflectance of the kth standard mode; The reference multispectral remote sensor data is used to solve the hyperspectral remote sensor channel apparent reflectance model corresponding to the reference multispectral remote sensor band, and the hyperspectral reconstruction parameters are obtained, which are expressed as follows: Where F' is the correction factor f i The inverse matrix of the correction factor is obtained by the inverse of ; The apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated is reconstructed according to the hyperspectral reconstruction parameters, and its expression is: R H =PHC+rH Among them, r H is the reconstructed residual matrix; PH is the standard mode apparent reflectance spectrum matrix, and the expression of its matrix elements is: According to the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated, the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined, and its expression is: Among them, D H is the quantized bit data obtained from the hyperspectral remote sensor data to be calibrated, R H is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

2. A cross-radiation calibration device for a satellite-borne hyperspectral remote sensor, characterized in that: include: A data acquisition module, used to acquire reference multispectral remote sensor data and to-be-calibrated hyperspectral remote sensor data; A correction factor module is used to obtain the correction factor of the reference multispectral remote sensor and the target hyperspectral remote sensor according to the reference multispectral remote sensor data and the hyperspectral remote sensor data to be calibrated; Model building module, used to build the surface reflectivity model of the ground object; A model solving module is used to solve the surface reflectance model of the ground object according to the reference multispectral remote sensor data and the correction factor to obtain the hyperspectral reconstruction parameters; A spectrum reconstruction module, used to reconstruct the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated according to the hyperspectral reconstruction parameters; A calibration module is used to determine the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated based on the hyperspectral remote sensor data to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated; The correction factor is expressed as: in, is the apparent reflectance of the i-th band of the reference multispectral remote sensor, and its expression is: in, is the spectral response function of the i-th band of the reference multispectral remote sensor, R(λ) is the apparent reflectance continuous spectrum of the calibration site; is the starting wavelength of the i-th band of the reference multispectral remote sensor, is the end wavelength of the i-th band of the reference multispectral remote sensor; is the apparent reflectance of the hyperspectral remote sensor to be calibrated corresponding to the i-th band of the reference multispectral remote sensor, and its expression is: Among them, R H (λ) is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated to the calibration site, which is obtained by merging the apparent reflectance continuous spectra of the hyperspectral remote sensor to be calibrated to the calibration site in multiple sub-bands. The expression is: in, is the starting wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the end wavelength of the jth band of the hyperspectral remote sensor to be calibrated, is the spectral response function of the jth band of the hyperspectral remote sensor to be calibrated; The surface reflectivity model of the ground object is expressed as follows: R I =C w P I，w +C v P I，v +C s P I，s +C4P I，4 +r I Among them, R I is the apparent reflectivity of the ground object in band I, P I，w is the formatted apparent reflectance spectrum of the water standard model, P I，v is the formatted apparent reflectance spectrum of the vegetation standard model, P I，s is the formatted apparent reflectance spectrum of the bare soil standard model, P I，4 The formatted apparent reflectance spectrum of the additional standard mode, C w is the water body standard model coefficient, C v is the vegetation standard model coefficient, C s is the bare soil standard model coefficient, C4 is the additional standard model coefficient, r I is the residual; The model solving module is used to solve the surface reflectance model of the ground object based on the reference multispectral remote sensor data and the correction factor to obtain the hyperspectral reconstruction parameters. It is specifically used for: According to the surface reflectance model of the ground object, the apparent reflectance model of the hyperspectral remote sensor channel corresponding to the reference multispectral remote sensor band is obtained, and its expression is: R H-M =F·(P M C+r) Where F is the correction factor f i The obtained correction factor matrix, C is the coefficient matrix, and r is the residual matrix; P M is the reference multispectral remote sensor channel apparent reflectance matrix, and the expression of its matrix elements is: Among them, P′ k (λ) is the hyperspectral apparent reflectance of the kth standard mode; The reference multispectral remote sensor data is used to solve the hyperspectral remote sensor channel apparent reflectance model corresponding to the reference multispectral remote sensor band, and the hyperspectral reconstruction parameters are obtained, which are expressed as follows: Where F' is the correction factor f i The inverse matrix of the correction factor is obtained by the inverse of ; The apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated is reconstructed according to the hyperspectral reconstruction parameters, and its expression is: R H =PHC+rH Among them, r H is the reconstructed residual matrix; PH is the standard mode apparent reflectance spectrum matrix, and the expression of its matrix elements is: According to the data of the hyperspectral remote sensor to be calibrated and the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated, the on-orbit absolute radiation calibration coefficient of the hyperspectral remote sensor to be calibrated is determined, and its expression is: Among them, D H is the quantized bit data obtained from the hyperspectral remote sensor data to be calibrated, R H is the apparent reflectance continuous spectrum of the hyperspectral remote sensor to be calibrated.

3. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the cross-radiation calibration method for a space-borne hyperspectral remote sensor as claimed in claim 1 is implemented.

4. 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 cross-radiation calibration method for a spaceborne hyperspectral remote sensor as claimed in claim 1 is implemented.

Images