Hyperspectral target height inversion method and device based on radiation transfer model

Abstract

The invention discloses a hyperspectral target height inversion method and device based on a radiation transfer model, and belongs to the field of remote sensing data analysis. The method comprises the following steps: acquiring a corresponding relation between atmospheric height and atmospheric temperature in the weather based on atmospheric data, and establishing a relation equation between target surface temperature and atmospheric temperature to determine the corresponding relation between atmospheric height and target surface temperature in the weather; sequentially simulating and calculating target spectrum entrance pupil radiation brightness and target center brightness temperature so as to transmit and determine a target center spectrum brightness temperature curve corresponding to each atmospheric height under the weather; and after wave band selection is carried out based on channel correlation of the target center spectrum brightness temperature curve corresponding to each atmospheric height, spectral angles of the real single-star image brightness temperature curve and the target center spectrum brightness temperature curve under each atmospheric height are respectively calculated, so that the atmospheric height with the highest spectral similarity is determined as the target height obtained by inversion. Therefore, the height inversion of the single-satellite imaging aerial target is realized.

Description

Technical Field

[0001] The present invention relates to the technical field of remote sensing data analysis, and particularly to a hyperspectral target height inversion method and device based on a radiative transfer model. Background Art

[0002] Due to their high-speed dynamic characteristics and potential security threats, airborne moving targets have become key objects of concern in the defense system. For such targets, it is urgent to deploy sophisticated airborne detection technologies to ensure rapid and accurate responses. Currently, the positioning technologies for airborne targets involve the comprehensive integration of multiple data sources such as satellites, ground radars, and optoelectronic sensors to achieve the comprehensive positioning and tracking of targets.

[0003] For the positioning of airborne moving targets, dual-star or multi-star joint positioning methods, or satellite-ground joint positioning are mostly used. However, although the dual-star or multi-star joint positioning method has a long enough ranging distance, it has high registration requirements and is greatly affected by the target state during shooting; satellite-ground joint positioning requires the target to be strictly within the common coverage range of the two, and the positioning area is limited.

[0004] Since the information sources that a single satellite can obtain are limited, there is currently no single-satellite positioning technology for airborne moving targets that can ensure a long enough ranging distance, is not restricted by the area, and can avoid high registration requirements.

[0005] Therefore, there is an urgent need to provide a hyperspectral target height inversion method and device based on a radiative transfer model. Summary of the Invention

[0006] In order to solve the problem that there is currently no height inversion method for airborne targets with single-satellite imaging, embodiments of the present invention provide a hyperspectral target height inversion method and device based on a radiative transfer model.

[0007] On the one hand, a hyperspectral target height inversion method based on a radiative transfer model is provided. The method includes:

[0008] Obtain the atmospheric data at this location during single-satellite imaging, and obtain the corresponding relationship between the atmospheric height and the atmospheric temperature under this weather based on the atmospheric data;

[0009] Based on the kinetic thermal equilibrium model, establish a relationship equation between the target surface temperature and the atmospheric temperature to determine the corresponding relationship between the atmospheric height and the target surface temperature under this weather;

[0010] Use the atmospheric radiative transfer model in the imaging state to sequentially simulate and calculate the target spectral pupil radiance and the target center brightness temperature, so as to transmit and determine the target center spectral brightness temperature curve corresponding to each atmospheric height under this weather;

[0011] After performing band selection based on the channel correlation of the target center spectral brightness temperature curves corresponding to each atmospheric altitude, calculate the spectral angles between the true single-star image brightness temperature curve and the target center spectral brightness temperature curves at each atmospheric altitude, respectively, to determine the atmospheric altitude with the highest spectral similarity as the retrieved target altitude; wherein, the spectral angle is used to characterize the similarity between spectra.

[0012] On the other hand, a hyperspectral target altitude retrieval device based on a radiative transfer model is provided, which is used to implement the steps described in any method embodiment of the specification. The device includes:

[0013] An acquisition unit, configured to acquire the atmospheric data at this location during single-star imaging, and based on the atmospheric data, acquire the corresponding relationship between the atmospheric altitude and the atmospheric temperature under this weather condition;

[0014] A determination unit, configured to establish a relationship equation between the target surface temperature and the atmospheric temperature based on the kinetic thermal equilibrium model, to determine the corresponding relationship between the atmospheric altitude and the target surface temperature under this weather condition;

[0015] A transfer unit, configured to use the atmospheric radiative transfer model in the imaging state to sequentially simulate and calculate the target spectral pupil radiance and the target center brightness temperature, so as to transfer and determine the target center spectral brightness temperature curves corresponding to each atmospheric altitude under this weather condition;

[0016] An inversion unit, configured to perform band selection based on the channel correlation of the target center spectral brightness temperature curves corresponding to each atmospheric altitude, and then calculate the spectral angles between the true single-star image brightness temperature curve and the target center spectral brightness temperature curves at each atmospheric altitude, respectively, to determine the atmospheric altitude with the highest spectral similarity as the retrieved target altitude; wherein, the spectral angle is used to characterize the similarity between spectra.

[0017] On the other hand, a computer device is provided. The computer device includes a memory and a processor. The memory is used to store a computer program, and the processor is used to execute the computer program stored on the memory to implement the steps of the above-mentioned method.

[0018] On the other hand, a computer-readable storage medium is provided. The storage medium stores a computer program, and when the computer program is executed by a processor, it implements the steps of the above-mentioned method.

[0019] On the other hand, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, it implements the steps of the above-mentioned method.

[0020] The technical solution provided by the present invention can at least bring the following beneficial effects:

[0021] Starting from the radiation transfer model, a simulation of the target center brightness temperature at different atmospheric heights is proposed. By separately calculating the heat balance equations at different heights on the surface of the airborne target and considering the influence of atmospheric radiation transfer, the target center spectral brightness temperature curves of the simulation image and the real single-star image are calculated. Finally, the wavelength channels with high correlation are filtered out, and the simulation image brightness temperature curve that is most similar to the real image brightness temperature curve at all atmospheric heights is selected by calculating the spectral angle to obtain the altitude inversion result, so as to realize the altitude inversion of the airborne target in single-star imaging. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0023] Figure 1 FIG. is a flowchart of a hyperspectral target altitude inversion method based on a radiation transfer model provided by an embodiment of the present invention;

[0024] Figure 2 FIG. is a structural diagram of a hyperspectral target altitude inversion device based on a radiation transfer model provided by an embodiment of the present invention;

[0025] Figure 3 FIG. is a hardware architecture diagram of a computer device provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

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

[0027] The following describes the specific implementation manners of the above concepts.

[0028] Please refer to Figure 1 , a hyperspectral target altitude inversion method based on a radiation transfer model provided by an embodiment of the present invention, the method includes:

[0029] Step 100: Obtain the atmospheric data at this location during single-star imaging, and obtain the corresponding relationship between the atmospheric height and the atmospheric temperature under this weather condition based on the atmospheric data;

[0030] Step 102: Based on the dynamic heat balance model, establish a relationship equation between the target surface temperature and the atmospheric temperature to determine the corresponding relationship between the atmospheric height and the target surface temperature under this weather;

[0031] Step 104: using the atmospheric radiation transmission model under the imaging state, the target spectral entrance pupil radiation brightness and the target center brightness temperature are simulated and calculated in sequence, so as to transfer and determine the target center spectral brightness temperature curve corresponding to each atmospheric height under this weather;

[0032] Step 106: After performing band selection based on the channel correlation of the target center spectral brightness temperature curve corresponding to each atmospheric altitude, the spectral angles of the real single star image brightness temperature curve and the target center spectral brightness temperature curve at each atmospheric altitude are calculated respectively, so as to determine the atmospheric altitude with the highest spectral similarity as the inverted target altitude; wherein the spectral angle is used to characterize the similarity between spectra.

[0033] In an embodiment of the present invention, starting from the radiation transfer model, a simulation of the target center brightness temperature at different atmospheric altitudes is proposed. By respectively calculating the heat balance equation at different altitudes of the aerial target surface and considering the influence of atmospheric radiation transfer, the target center spectral brightness temperature curves of the simulated image and the real single star image are calculated. Finally, the band channels with high correlation are screened out, and the simulated image brightness temperature curve that is most similar to the real image brightness temperature curve at all atmospheric altitudes is selected by calculating the spectral angle to obtain the height inversion result, thereby realizing the height inversion of single-star imaging aerial targets.

[0034] Described below Figure 1 How the various steps are performed.

[0035] For step 100:

[0036] For the height inversion of hyperspectral imaging targets, long-wave infrared data is selected to remove the influence of direct sunlight as much as possible, and the simulation of the target center brightness temperature under different height parameters is considered.

[0037] Based on the atmospheric data of this place during single-star imaging, after determining parameters such as atmospheric mode, atmospheric humidity, pressure data, visibility, etc., the atmospheric temperature data at different atmospheric altitudes are obtained, and the spectral transmittance, atmospheric spectral upradiation and other data are calculated in combination with MODTRAN software.

[0038] Regarding step 102:

[0039] In some embodiments, the relationship equation between the target surface temperature and the atmospheric temperature is:

[0040]

[0041] Where, T aircraft Indicates the target surface temperature, T atmrepresents the atmospheric temperature, C represents the compensation coefficient, β represents the motion coefficient, γ represents the flow field coefficient, and v aircraft represents the flight speed of the target, and T sun represents the temperature gain caused by direct solar radiation, and θ represents the solar zenith angle.

[0042] In this embodiment, the spectral radiation characteristics of an airborne target in the long-wave infrared band are closely related to the target surface temperature and emissivity. Considering the equation of dynamic thermal equilibrium, an equation relating the target surface temperature and the atmospheric temperature can be established by combining empirical parameters to utilize the atmospheric temperature to determine the corresponding relationship between the atmospheric altitude and the target surface temperature under this weather condition.

[0043] Regarding step 104:

[0044] In some embodiments, step 104 includes steps S1 - S2:

[0045] Step S1: Combining the Planck blackbody radiation formula, based on the target surface temperature, wavelength, and emissivity, calculate the target spectral radiance at the entrance pupil to obtain the target spectral radiance curve corresponding to each atmospheric altitude under this weather condition.

[0046] In this embodiment, considering data such as the spectral transmittance and upward atmospheric spectral radiation calculated by the MODTRAN software, and then combining the Planck blackbody radiation formula, the target spectral radiance at the entrance pupil can be calculated:

[0047]

[0048] In the formula, L pixel represents the target spectral radiance at the entrance pupil, λ represents the wavelength, T aircraft represents the target surface temperature, ε represents the emissivity, and c1 and c2 are the first and second radiation constants respectively. Therefore, by using the target surface temperature, the target spectral radiance curve corresponding to each atmospheric altitude under this weather condition can be determined.

[0049] Step S2: Use the Planck blackbody radiation formula for inverse derivation. Based on the target spectral radiance at the entrance pupil, Planck's constant, the speed of light, and the wavelength, calculate the target central brightness temperature to determine the target central spectral brightness temperature curve corresponding to each atmospheric altitude under this weather condition.

[0050] In this embodiment, the brightness temperature refers to the temperature of a blackbody when the spectral radiation intensity of an actual object is equal to that of a blackbody at a specific wavelength. Therefore, the target central brightness temperature can be obtained by inverse derivation using the Planck blackbody radiation formula:

[0051]

[0052] Wherein, T represents the target center brightness temperature, h represents Planck's constant, c represents the speed of light, k represents Boltzmann's constant, λ represents the wavelength, and L pixel represents the target spectral entrance pupil radiance. Thus, different altitude parameters correspond to different target center spectral brightness temperature curves.

[0053] Regarding step 106:

[0054] In some embodiments, step 106 may include B1 - B4:

[0055] B1. Calculate the correlation between every two wavelength channels in the target center spectral brightness temperature curves at all atmospheric altitudes using the Pearson coefficient, and screen out relatively independent band combinations;

[0056] B2. Calculate the true image brightness temperature curve of single - satellite imaging under this weather condition, and based on the screened band combinations, also screen the bands of the true image brightness temperature curve;

[0057] B3. Calculate the spectral angles between the screened true image brightness temperature curve and the target center spectral brightness temperature curves at each atmospheric altitude respectively; the smaller the spectral angle, the higher the spectral similarity;

[0058] B4. Determine the atmospheric altitude corresponding to the minimum spectral angle as the target altitude where the airborne target is located.

[0059] In this embodiment, in order to enhance the effectiveness of subsequent spectral angle calculations, after calculating the spectral brightness temperature curves at different atmospheric altitudes, use the Pearson coefficient to calculate the correlation between the band channels of the hyperspectral data to achieve feature band selection. It can be understood that each spectral brightness temperature curve corresponding to an atmospheric altitude has several band channels. As the abscissa points, band screening will reduce the abscissa points.

[0060] The correlation between every two wavelength channels is calculated by the following formula:

[0061]

[0062] Wherein, r is the correlation value between channels, ranging from - 1 to 1. 1 represents a perfect positive correlation, - 1 represents a perfect negative correlation, 0 represents no linear correlation, X i and Y i are the data points at the i - th atmospheric altitude of the X and Y two channels, X mean and Y mean are the averages of the X and Y two channels at all atmospheric altitudes, and n is the number of simulations of atmospheric altitudes.

[0063] Sort the correlation values between pairwise wavelength channels, and select the band combinations with correlation values close to 0, that is, relatively independent band combinations, and eliminate the remaining bands.

[0064] In step B2, obtain the real image of single - star imaging for this weather. The brightness temperature curve of the real image can be calculated using the steps in step 104. Based on the band combination selected in B1, the bands of the brightness temperature curve of the real image are also screened.

[0065] In steps B3 and B4, after band selection, calculate the spectral angle between the screened brightness temperature curve of the real image and the spectral brightness temperature curve of the target center at each atmospheric height:

[0066]

[0067] Among them, SAM is the spectral angle, representing the correlation between spectra. 0 means that the two spectra are exactly the same, 90 means that the spectra are completely different, X i is the value of the i - th band under this atmospheric height parameter, Y i is the value of the i - th band of the brightness temperature spectral curve of the real image, and n is the number of remaining bands after screening. Select the atmospheric height with the smallest SAM as the target height obtained by inversion.

[0068] In this way, the selection of the closest brightness temperature curve is achieved. Since the height parameter is in one - to - one correspondence with the brightness temperature curve of the central pixel, the hyperspectral target height inversion based on the radiative transfer principle is completed.

[0069] In summary, this solution utilizes the transitivity of the height parameter in the radiative transfer model of hyperspectral remote - sensing data. By simulating the brightness temperature curves of the target at different heights and selecting the curve closest to the brightness temperature of the real image, height inversion is achieved, which has the advantages of optimizing the target imaging process, enhancing the interpretability of remote - sensing data, and better meeting user requirements.

[0070] Please refer to Figure 2 , an embodiment of the present invention provides a hyperspectral target height inversion device based on a radiative transfer model, which is used to implement the steps of any method embodiment in the specification. The device includes:

[0071] An acquisition unit 201, configured to acquire the atmospheric data at this location during single - star imaging, and obtain the corresponding relationship between the atmospheric height and the atmospheric temperature under this weather based on the atmospheric data;

[0072] A determination unit 202, configured to establish a relationship equation between the target surface temperature and the atmospheric temperature based on the kinetic thermal equilibrium model to determine the corresponding relationship between the atmospheric height and the target surface temperature under this weather;

[0073] The transmission unit 203 is configured to use the atmospheric radiation transmission model in the imaging state to sequentially simulate and calculate the target spectral pupil radiance and the target center brightness temperature, so as to transmit and determine the target center spectral brightness temperature curve corresponding to each atmospheric altitude in this weather condition;

[0074] The inversion unit 204 is configured to perform band selection based on the channel correlation of the target center spectral brightness temperature curve corresponding to each atmospheric altitude, and then calculate the spectral angles between the true single-star image brightness temperature curve and the target center spectral brightness temperature curves at each atmospheric altitude, so as to determine that the atmospheric altitude with the highest spectral similarity is the target altitude obtained by inversion; wherein, the spectral angle is used to characterize the similarity between spectra.

[0075] In an embodiment of the present invention, the relationship equation between the target surface temperature and the atmospheric temperature in the determination unit 202 is:

[0076]

[0077] In the formula, T aircraft represents the target surface temperature, T atm represents the atmospheric temperature, C represents the compensation coefficient, β represents the motion coefficient, γ represents the flow field coefficient, v aircraft represents the flight speed of the target, T sun represents the temperature gain caused by direct solar radiation, and θ represents the solar zenith angle.

[0078] In an embodiment of the present invention, the transmission unit 203 is configured to execute:

[0079] Combined with the Planck blackbody radiation formula, based on the target surface temperature, wavelength, and emissivity, calculate the target spectral pupil radiance to obtain the target spectral pupil radiance curve corresponding to each atmospheric altitude in this weather condition;

[0080] Use the Planck blackbody radiation formula for backtracking, and based on the target spectral pupil radiance, Planck constant, speed of light, and wavelength, calculate the target center brightness temperature to determine the target center spectral brightness temperature curve corresponding to each atmospheric altitude in this weather condition.

[0081] In an embodiment of the present invention, the target spectral pupil radiance in the transmission unit 203 is calculated in the following manner:

[0082]

[0083] In the formula, L pixel represents the target spectral pupil radiance, λ represents the wavelength, T aircraft represents the target surface temperature, ε represents the emissivity, and c1 and c2 are the first and second radiation constants respectively;

[0084] The target center brightness temperature is calculated in the following manner:

[0085]

[0086] In the formula, T represents the target center brightness temperature, h represents Planck's constant, c represents the speed of light, k represents Boltzmann's constant, λ represents the wavelength, and L pixel represents the target spectral entrance pupil radiance.

[0087] In an embodiment of the present invention, the inversion unit 204 is used to perform:

[0088] Calculate the correlation between every two wavelength channels in the target center spectral brightness temperature curve at all atmospheric heights using the Pearson coefficient, and screen out relatively independent band combinations;

[0089] Calculate the true image brightness temperature curve of single-star imaging under this weather condition, and based on the screened band combinations, also screen the bands of the true image brightness temperature curve;

[0090] Calculate the spectral angles between the screened true image brightness temperature curve and the target center spectral brightness temperature curve at each atmospheric height respectively; the smaller the spectral angle, the higher the spectral similarity;

[0091] Determine the atmospheric height corresponding to the minimum spectral angle as the target height where the aerial target is located.

[0092] In an embodiment of the present invention, the correlation between every two wavelength channels in the inversion unit 204 is calculated by the following formula:

[0093]

[0094] In the formula, r is the correlation value between channels, ranging from -1 to 1. 1 represents a perfect positive correlation, -1 represents a perfect negative correlation, and 0 represents no linear correlation. X i and Y i are the data points at the i-th atmospheric height of the X and Y channels, and X mean and Y mean are the averages of the X and Y channels at all atmospheric heights, and n is the number of simulations of atmospheric heights.

[0095] It should be noted that: for the hyperspectral target height inversion device based on the radiation transfer model provided in the above embodiments, only the above-mentioned division of each functional unit is used for illustration. In practical applications, the above functions can be allocated to different functional units according to needs, that is, the internal structure of the device is divided into different functional units to complete all or part of the functions described above. In addition, the above device embodiments and method embodiments belong to the same concept, and the specific implementation process can be seen in the method embodiments, which will not be elaborated here.

[0096] Embodiments of the present application also provide a computer device. Please refer to Figure 3 , the computer device includes a processor and a memory. At least one instruction, at least one program, a code set or an instruction set is stored in the memory. The at least one instruction, at least one program, the code set or the instruction set is loaded and executed by the processor to implement the hyperspectral target height inversion method based on the radiative transfer model provided in the above method embodiments.

[0097] Embodiments of the present application also provide a computer-readable storage medium. At least one instruction, at least one program, a code set or an instruction set is stored on the computer-readable storage medium. The at least one instruction, at least one program, the code set or the instruction set is loaded and executed by the processor to implement the hyperspectral target height inversion method based on the radiative transfer model provided in the above method embodiments.

[0098] Embodiments of the present application also provide a computer program product. The computer program product includes a computer program. The processor of the computer device reads the computer program from the computer-readable storage medium, and the processor executes the computer program, so that the computer device executes the hyperspectral target height inversion method according to any one of the above embodiments based on the radiative transfer model.

[0099] For convenience of description, when describing the above system or device, it is divided into various modules or units according to functions for description. Of course, when implementing the present application, the functions of each unit can be implemented in one or more software and / or hardware.

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

[0101] Finally, it should also be noted that in this text, relational terms such as first, second, third, and fourth are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprising", "including" or any other variant thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device comprising a series of elements not only includes those elements, but also includes other elements not expressly listed, or further includes elements inherent to such process, method, article or device. Without further limitation, an element defined by the statement "comprising an..." does not exclude the presence of additional identical elements in the process, method, article or device comprising the element.

[0102] The above are only the preferred embodiments of the present application. It should be pointed out that for those of ordinary skill in the art, without departing from the principle of the present application, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present application.

Claims

1. A hyperspectral target height inversion method based on a radiation transfer model, characterized in that: include: Acquire the atmospheric data of the place when the single star is imaged, and acquire the corresponding relationship between the atmospheric height and the atmospheric temperature under the weather conditions based on the atmospheric data; Based on the dynamic heat balance model, the relationship equation between the target surface temperature and the atmospheric temperature is established to determine the corresponding relationship between the atmospheric height and the target surface temperature under this weather condition; Using the atmospheric radiation transfer model under imaging conditions, the target spectral entrance pupil radiation brightness and target center brightness temperature are simulated and calculated in sequence to transfer and determine the target center spectral brightness temperature curve corresponding to each atmospheric height under this weather condition; After band selection based on the channel correlation of the target center spectral brightness temperature curve corresponding to each atmospheric altitude, the spectral angles of the real single star image brightness temperature curve and the target center spectral brightness temperature curve at each atmospheric altitude are calculated respectively to determine the atmospheric altitude with the highest spectral similarity as the inverted target altitude; wherein the spectral angle is used to characterize the similarity between spectra.

2. The method according to claim 1, characterized in that The relationship equation between the target surface temperature and the atmospheric temperature is: Where, T aircraft Indicates the target surface temperature, T atm represents the atmospheric temperature, C represents the compensation coefficient, β represents the motion coefficient, γ represents the flow field coefficient, and v aircraft Indicates the target's flight speed, T sun represents the temperature gain caused by direct sunlight, and θ represents the solar zenith angle.

3. The method according to claim 1, characterized in that Using the atmospheric radiation transfer model under imaging conditions, the target spectral entrance pupil radiation brightness and target center brightness temperature are simulated and calculated in sequence to transfer and determine the target center spectral brightness temperature curve corresponding to each atmospheric height under this weather, including: Combined with Planck's blackbody radiation formula, based on the target surface temperature, wavelength and emissivity, the target spectrum entrance pupil radiation brightness is calculated to obtain the target spectrum entrance pupil radiation brightness curve corresponding to each atmospheric height under this weather; The target center brightness temperature is calculated based on the target spectral entrance pupil radiation brightness, Planck constant, speed of light and wavelength using Planck's blackbody radiation formula to determine the target center spectral brightness temperature curve corresponding to each atmospheric altitude under this weather condition.

4. The method according to claim 3, characterized in that The target spectral entrance pupil radiance is calculated as follows: Where, L pixel represents the target spectrum entrance pupil radiance, λ represents the wavelength, T aircraft represents the target surface temperature, ε represents the emissivity, c1 and c2 are the first and second radiation constants respectively; The target center brightness temperature is calculated as follows: Where T is the target center brightness temperature, h is Planck's constant, c is the speed of light, k is the Boltzmann constant, λ is the wavelength, and L is pixel Represents the target spectral entrance pupil radiance.

5. The method according to claim 1, characterized in that After performing band selection based on the channel correlation of the target center spectral brightness temperature curve corresponding to each atmospheric altitude, the spectral angles of the real single star image brightness temperature curve and the target center spectral brightness temperature curve at each atmospheric altitude are calculated respectively, so as to determine the atmospheric altitude with the highest spectral similarity as the inverted target altitude, including: The Pearson coefficient is used to calculate the correlation between each two wavelength channels in the target center spectral brightness temperature curve at all atmospheric altitudes, and relatively independent band combinations are selected; Calculate the real image brightness temperature curve of single star imaging under this weather, and based on the screened band combination, screen the bands of the real image brightness temperature curve as well; The spectral angles of the brightness temperature curve of the filtered real image and the brightness temperature curve of the target center spectrum at each atmospheric height are calculated respectively; the smaller the spectral angle, the higher the spectral similarity; The atmospheric height corresponding to the minimum spectral angle is determined as the target height of the aerial target.

6. The method according to claim 5, characterized in that The correlation between each two wavelength channels is calculated by the following formula: Where r is the correlation value between channels, ranging from -1 to 1, 1 represents a complete positive correlation, -1 represents a complete negative correlation, and 0 represents no linear correlation. i With Y i is the data point at the ith atmospheric height of the two channels X and Y, X mean With Y mean is the average value of all atmospheric heights for the X and Y channels, and n is the number of simulated atmospheric heights.

7. A hyperspectral target height inversion device based on a radiation transfer model, used to implement the steps of any of the methods described in claims 1-6 above, characterized in that: include: An acquisition unit, used to acquire atmospheric data of the location when a single star is imaged, and acquire a corresponding relationship between the atmospheric height and the atmospheric temperature under the weather conditions based on the atmospheric data; A determination unit is used to establish a relationship equation between the target surface temperature and the atmospheric temperature based on a dynamic heat balance model, so as to determine the corresponding relationship between the atmospheric height and the target surface temperature under this weather; The transmission unit is used to simulate and calculate the target spectral entrance pupil radiation brightness and the target center brightness temperature in sequence by using the atmospheric radiation transmission model under the imaging state, so as to transmit and determine the target center spectral brightness temperature curve corresponding to each atmospheric height under this weather; The inversion unit is used to calculate the spectral angles of the brightness temperature curve of the real single star image and the brightness temperature curve of the target center spectrum at each atmospheric height after performing band selection based on the channel correlation of the brightness temperature curve of the target center spectrum at each atmospheric height, so as to determine the atmospheric height with the highest spectral similarity as the inverted target height; wherein the spectral angle is used to characterize the similarity between the spectra.

8. A computer device, characterized in that: The computer device includes a memory and a processor, the memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory to implement the steps of any one of the methods described in claims 1-6.

9. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method described in any one of claims 1 to 6 are implemented.

10. A computer program product, characterized in that The method comprises a computer program, wherein when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images