CN-121977694-A - Method, apparatus and storage medium for correcting relative radiation of satellite-borne hyperspectral camera

CN121977694ACN 121977694 ACN121977694 ACN 121977694ACN-121977694-A

Abstract

The application provides a relative radiation correction method, relative radiation correction equipment and relative radiation correction storage medium of a satellite-borne hyperspectral camera, wherein the relative radiation correction method comprises the steps of obtaining hyperspectral original image data and dark current observation data for statistical analysis; determining dark current average response values of each probe element on continuous spectrum channels according to dark current observation data, obtaining hyperspectral image data of dark current, calculating to obtain spectrum radiation response reference values of each probe element, calculating to obtain relative radiation correction coefficients of each probe element on all spectrum channels according to the spectrum radiation response reference values of each probe element, and carrying out relative radiation correction on all hyperspectral original image data obtained by a satellite-borne hyperspectral camera within a preset time range according to the relative radiation correction coefficients. The embodiment of the application can realize the effect of improving the calculation efficiency on the premise of ensuring the correction precision, thereby realizing the effective balance of the calculation efficiency and the correction precision.

Inventors

WANG XINMENG
Lan Kongkong
HAN QIJIN
HUANG XIAOYU
WANG HAIYAN

Assignees

中国资源卫星应用中心

Dates

Publication Date: 20260505
Application Date: 20260206

Claims (10)

1. A method of relative radiation correction for an on-board hyperspectral camera, the on-board hyperspectral camera comprising a plurality of probe elements, the method comprising: Acquiring hyperspectral original image data for statistical analysis and dark current observation data for evaluating dark current influence in a preset time range, wherein the hyperspectral original image data comprises original response values of each probe element on different spectrum channels; Determining a dark current average response value of each probe element on a continuous spectrum channel according to the dark current observation data, wherein the dark current average response value is used for reflecting the influence degree of the dark current level on the probe element under the corresponding spectrum channel; Obtaining hyperspectral image data of dark current according to an original response value and the dark current average response value in the hyperspectral original image data, wherein the hyperspectral image data of dark current comprises a plurality of dark current response values; carrying out statistical sorting on the hyperspectral image data of the dark current removal, and calculating to obtain a spectral radiation response reference value of each probe element, wherein the spectral radiation response reference value is used for representing the radiation response characteristic of the corresponding probe element under a continuous spectrum channel; According to the spectrum radiation response reference value of each probe element, calculating to obtain the relative radiation correction coefficient of each probe element on all spectrum channels; And carrying out relative radiation correction on all hyperspectral original image data acquired by the satellite-borne hyperspectral camera within a preset time range according to the relative radiation correction coefficient.
2. A method of relative radiation correction according to claim 1, wherein said determining dark current average response values for each probe element over a continuous spectrum channel from said dark current observations comprises: Determining the number of pixels of each probe element under a continuous spectrum channel; And calculating to obtain the dark current average response value of each probe element according to the dark current observation data under the corresponding spectrum channel and the pixel quantity.
3. The method according to claim 1, wherein obtaining the hyperspectral image data with dark current removed according to the original response value and the dark current average response value in the hyperspectral original image data comprises: subtracting the average response value of the dark current from each original response value of the probe element corresponding to each spectrum channel to obtain hyperspectral image data of dark current.
4. The method according to claim 1, wherein the statistically ordering the hyperspectral image data of the dark current and calculating the spectral radiation response reference value of each probe cell comprises: Sorting the dark current response values along the track direction according to the magnitude relation to obtain a sorted dark current response value set; removing partial dark current response values with large influence factors in the dark current response value set to obtain a screened dark current response value set; and calculating the average value of all the dark current response values in the screened dark current response value set to obtain a spectral radiation response reference value corresponding to each probe element.
5. The method according to claim 4, wherein removing the partial dark current response values with the larger influencing factors from the dark current response value set to obtain the filtered dark current response value set comprises: Determining the total data range corresponding to the sequenced dark current response value set; Removing dark current response values in a first preset range ranked in front and dark current response values in a second preset range ranked behind in the ranked dark current response value set; and removing the first preset range and the second preset range from the total data range according to the filtered dark current response value set in the middle range of the total data range.
6. The method according to claim 1, wherein the calculating the relative radiation correction coefficient of each probe element on all the spectral channels according to the spectral radiation response reference value of each probe element comprises: Obtaining the average value of the spectrum radiation response reference values of all the probe elements according to the spectrum radiation response reference value corresponding to each probe element; And calculating the relative radiation correction coefficient of the corresponding probe element on each continuous spectrum channel based on the spectrum radiation response reference value of each probe element and the spectrum radiation response reference value average value of all the probe elements by adopting a linear method.
7. The method according to claim 1, wherein after said statistically ordering the hyperspectral image data of the dark current, and calculating a spectral radiation response reference value for each probe cell, the method further comprises: And carrying out filtering processing on the spectral radiation response reference value of each probe element by adopting least square fitting processing in the space dimension and adopting Gaussian filtering in the spectral dimension.
8. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the relative radiation correction method according to any one of claims 1 to 7 when the computer program is executed.
9. A computer-readable storage medium storing computer-executable instructions for causing a computer to perform the relative radiation correction method as claimed in any one of claims 1 to 7.
10. A computer program product comprising a computer program, characterized in that the computer program, when executed by a processor, implements the relative radiation correction method of any one of claims 1 to 7.

Description

Method, apparatus and storage medium for correcting relative radiation of satellite-borne hyperspectral camera Technical Field The invention relates to the technical field of remote sensing, in particular to a relative radiation correction method, equipment and a storage medium of a satellite-borne hyperspectral camera. Background The push-broom hyperspectral camera is currently loaded on a remote sensing satellite produced in each country, is applied to the fields of military and civil use, and is applied to the fields of aerospace, homeland resource investigation, environment monitoring, atmosphere monitoring and the like. However, due to the limitations of manufacturing process and materials, the problem of inconsistent output response exists among the probe elements of the hyperspectral camera, and the inconsistency is usually marked as vertical stripes on the image, so that the imaging quality is seriously affected, and the application and development of the hyperspectral camera are limited. To solve the problem of inconsistent output responses, it is often necessary to perform relative radiation correction on the hyperspectral camera, and adjust the output of each probe element to a consistent reference by establishing a correction model, so as to improve the quality and usability of the image data. Relative radiation correction methods in the related art, such as laboratory relative radiation correction methods, ground relative radiation correction methods, and statistical relative radiation correction methods. The laboratory relative radiation correction method relies on integrating sphere to construct standard signals, and uses integrating sphere as input and sensor as output to calculate relative radiation correction coefficient. But the device state changes during orbit after satellite transmission, and thus this method is not applicable to the in-orbit data. The ground relative radiation correction method mainly relies on outdoor calibration, and under the sun illumination, an on-orbit satellite obtains a calibration image by observing uniform targets distributed on the ground, but the method needs to lay full-view targets with different reflectivities, and has the advantages of great workload and low processing efficiency. The statistical relative radiation correction method is divided into conventional statistics and yaw data statistics, but the former has long statistical time and low timeliness, and the latter obtains a nonlinear relative radiation correction coefficient and has poor effect in a middle-low end response interval. Therefore, how to improve the correction effect of the relative radiation and increase the calculation efficiency is a problem to be solved. Disclosure of Invention The invention aims to at least solve one of the technical problems in the prior art, and provides a relative radiation correction method, equipment and a storage medium of a satellite-borne hyperspectral camera, which can achieve the effect of improving the calculation efficiency on the premise of ensuring the correction precision, thereby achieving the effective balance of the calculation efficiency and the correction precision. In a first aspect, an embodiment of the present invention provides a method for correcting relative radiation, where the satellite-borne hyperspectral camera includes a plurality of probe elements, and the method for correcting relative radiation includes: Acquiring hyperspectral original image data for statistical analysis and dark current observation data for evaluating dark current influence in a preset time range, wherein the hyperspectral original image data comprises original response values of each probe element on different spectrum channels; Determining a dark current average response value of each probe element on a continuous spectrum channel according to the dark current observation data, wherein the dark current average response value is used for reflecting the influence degree of the dark current level on the probe element under the corresponding spectrum channel; Obtaining hyperspectral image data of dark current according to an original response value and the dark current average response value in the hyperspectral original image data, wherein the hyperspectral image data of dark current comprises a plurality of dark current response values; carrying out statistical sorting on the hyperspectral image data of the dark current removal, and calculating to obtain a spectral radiation response reference value of each probe element, wherein the spectral radiation response reference value is used for representing the radiation response characteristic of the corresponding probe element under a continuous spectrum channel; According to the spectrum radiation response reference value of each probe element, calculating to obtain the relative radiation correction coefficient of each probe element on all spectrum channels; And carrying out relative radiation correction on all