CN113378419B - Infrared imaging simulation method based on MODTRAN optimization - Google Patents
Infrared imaging simulation method based on MODTRAN optimization Download PDFInfo
- Publication number
- CN113378419B CN113378419B CN202110934162.2A CN202110934162A CN113378419B CN 113378419 B CN113378419 B CN 113378419B CN 202110934162 A CN202110934162 A CN 202110934162A CN 113378419 B CN113378419 B CN 113378419B
- Authority
- CN
- China
- Prior art keywords
- degrees
- zenith angle
- range
- observation
- taking
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention discloses an infrared imaging simulation method based on MODTRAN optimization, which comprises the following steps: s1: collecting atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance and sky background irradiance data, and determining a zenith angle needing to be optimized according to the collected atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance and sky background irradiance data; s2: observing zenith angle optimization and solar zenith angle optimization simultaneously through the optimized zenith angle, and determining a specific zenith angle by inquiring an optimization method table; s3: and calling MODTRAN software to calculate a generated result through a target infrared irradiance formula received at the entrance pupil of the detector, and outputting the result.
Description
Technical Field
The invention relates to an infrared imaging simulation method, in particular to an infrared imaging simulation method based on MODTRAN optimization.
Background
In order to enable the sky background radiation intensity spectrum distribution to finally generate a background radiation power distribution image on the infrared detector focal plane, all influencing factors of sky background radiation need to be analyzed, and an accurate radiation energy transmission model of background radiation to the detector optical system focal plane is established. The infrared simulation is carried out by using MODTRAN software, and zenith angles are optimized, because when an observation geometric model is set, if a general observation range is large, the calculated amount is large, and the repeated work is high. The precision is slightly reduced, a radiation energy transmission model of background radiation to a focal plane of an optical system of the detector is established, a radiation power calculation method of the sky background on the focal plane is provided, and an infrared image of the sky background is generated by using the radiation power distribution of the sky background. And in the infrared simulation modeling process, when an observation geometric condition is established, the existing zenith angle range is optimized under the condition of large observation zenith angle or large solar zenith angle range to obtain one or more reasonable zenith angle in the range.
Disclosure of Invention
The invention aims to solve the technical problems that in the process of using MODTRAN atmospheric radiation transmission simulation software, when an observation geometric model is set, if the observation range is large, the calculated amount is large, and the repeated work is high. The accuracy is slightly reduced, and the purpose is to provide an infrared imaging simulation method based on MODTRAN optimization, so as to solve the problems.
The invention is realized by the following technical scheme: s1: determining a zenith angle type to be optimized according to five types of atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance and sky background irradiance, wherein the optimization type comprises an observation zenith angle and a solar zenith angle; s2: optimizing the observation zenith angle or optimizing the range of the solar zenith angle according to the optimized zenith angle type, and determining a specific zenith angle by inquiring an optimization method table; s3: and (4) calling MODTRAN software to set observation geometric conditions according to the optimized zenith angles, calculating by using a target infrared irradiance formula received at the entrance pupil of the detector to generate a result, and outputting the result.
Further, the optimization method table in step S2 includes three categories, which are calculation types, optimization methods, and remark information, where the calculation types include atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance, and sky background irradiance.
Furthermore, the optimization mode of the radiation brightness of the atmospheric path is that 10 degrees are taken when the zenith angle of observation is within the range of 0-19 degrees; when the observation zenith angle is within the range of 20-34 degrees, taking 27 degrees; when the observation zenith angle is within the range of 35-44 degrees, taking 40 degrees; when the observation zenith angle is within the range of 45-54 degrees, taking 50 degrees; when the observation zenith angle is within the range of 55-64 degrees, 60 degrees are taken; when the observation zenith angle is within the range of 65-79 degrees, taking down the odd observation zenith angles; when the zenith angle is observed to exceed 90 degrees, the test piece is symmetrical about the 90-degree center.
Further, the optimization mode of the atmospheric transmittance is that when the observation zenith angle is in the range of 0-29 degrees, 15 degrees are taken; when the observation zenith angle is within the range of 30-44 degrees, taking 37 degrees; when the observation zenith angle is within the range of 45-54 degrees, taking 50 degrees; when the observation zenith angle is within the range of 55-59 degrees, taking 57 degrees; when the observation zenith angle is within the range of 60-69 degrees, downwards taking an even number of observation zenith angles; when the zenith angle is observed to exceed 90 degrees, the test piece is symmetrical about the 90-degree center.
Further, the optimization mode of the solar radiation illumination is that when the sun zenith angle is in the range of 0-30 degrees, 15 degrees are taken; taking Esun _30- (Esun30-Esun40) × (X-30) when the zenith angle of the sun is within the range of 30-39 degrees; taking Esun _40- (Esun _40-Esun _ 50) × (X-40) when the zenith angle of the sun is in the range of 40-49 degrees; taking Esun _50- (Esun _50-Esun _ 60) × (X-50) when the zenith angle of the sun is within the range of 50-59 degrees; taking Esun _60- (Esun _60-Esun _ 70) (X-60) when the zenith angle of the sun is in the range of 60-69 degrees; taking Esun _70- (Esun _70-Esun _ 80) × (X-70) when the zenith angle of the sun is in the range of 70-79 degrees; wherein Esun _ x is the solar irradiance at x degrees.
Furthermore, the optimization mode of the earth radiation illumination is that when the observation zenith angle is within the range of 0-29 degrees, 15 degrees are taken; when the observation zenith angle is within the range of 30-49 degrees, taking 40 degrees; when the observation zenith angle is within the range of 50-64 degrees, 60 degrees are taken; when the observation zenith angle is within the range of 65-79 degrees, taking down odd number zenith angles; when the zenith angle is observed to exceed 90-180 degrees, the observation is symmetrical about the 90-degree center.
Furthermore, the optimization mode of the sky background irradiance is that 15 degrees are taken when the observed zenith angle is within the range of 0-19 degrees; when the observation zenith angle is within the range of 20-29 degrees, 25 degrees are taken; when the observation zenith angle is within the range of 30-39 degrees, taking 35 degrees; when the observation zenith angle is within the range of 40-49 degrees, taking 45 degrees; when the observation zenith angle is within the range of 50-69 degrees, downwards taking even number zenith angles; when the zenith angle is observed to exceed 90-180 degrees, the observation is symmetrical about the 90-degree center.
Further, the target infrared irradiance formula in step S3 is:
wherein L ispathRadiating luminance for an atmospheric path; l isselfIs the intrinsic radiance of the target; esunThe irradiance of the sun on the surface of the target skin; eearthThe irradiance of the earth on the surface of the target skin; eskyThe radiation illumination of sky background radiation on the surface of a target skin;is the atmospheric average transmittance of the path, dimensionless;is the reflectivity of the skin to ambient radiation. Here the BRDF coefficient, sr-1。
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. according to the infrared imaging simulation method based on MODTRAN optimization, the following method can be performed on the observation zenith angle, the solar zenith angle or the observation wave band when the observation geometric model is set when calculating different types of radiometric degree or transmittance in the process of using MODTRAN atmospheric radiation transmission simulation software. When the observation range is large, the method can effectively reduce the calculation amount without influencing the calculation result precision.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
FIG. 1 is a simulation flow chart of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
Examples
As shown in fig. 1, the infrared imaging simulation method based on MODTRAN optimization of the present invention includes the following steps: s1: determining a zenith angle type to be optimized according to five types of atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance and sky background irradiance, wherein the optimization type comprises an observation zenith angle and a solar zenith angle; s2: optimizing the observation zenith angle or optimizing the range of the solar zenith angle according to the optimized zenith angle type, and determining a specific zenith angle by inquiring an optimization method table; s3: and (4) calling MODTRAN software to set observation geometric conditions according to the optimized zenith angles, calculating by using a target infrared irradiance formula received at the entrance pupil of the detector to generate a result, and outputting the result.
The target infrared irradiance formula in step S3 is:
wherein L ispathRadiating luminance for an atmospheric path; l isselfIs the intrinsic radiance of the target; esunThe irradiance of the sun on the surface of the target skin; eearthThe irradiance of the earth on the surface of the target skin; eskyThe radiation illumination of sky background radiation on the surface of a target skin;mean atmospheric permeability for pathwayThe excess rate is dimensionless;is the reflectivity of the skin to ambient radiation. Here the BRDF coefficient, sr-1。
The optimization method table in step S2 is shown in table 1:
TABLE 1 optimization method Table
In actual use, the optimal optimized observation zenith angle and the optimal optimized observation zenith angle are selected through the table 1, large-angle and large-range calculation can be rapidly carried out, the calculation speed is increased, the calculation precision is guaranteed, the output speed in MODTRAN software can be effectively improved and optimized, and the working efficiency is improved.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (2)
1. An infrared imaging simulation method based on MODTRAN optimization is characterized by comprising the following steps:
s1: determining a zenith angle type to be optimized according to five types of atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance and sky background irradiance, wherein the optimization type comprises an observation zenith angle and a solar zenith angle;
s2: determining a specific zenith angle by inquiring an optimization method table;
s3: setting an observation geometric condition according to the optimized zenith angle by calling MODTRAN software, calculating by using a target infrared irradiance formula received at an entrance pupil of the detector to generate a result, and outputting the result;
the optimization method table in the step S2 includes three major categories, i.e., calculation type, optimization mode and remark information, wherein the calculation type includes atmospheric path radiance, atmospheric transmittance, solar radiance, earth radiance and sky background irradiance;
the optimization mode of the atmospheric path radiation brightness is that when the observation zenith angle is in the range of 0-19 degrees, 10 degrees are taken; when the observation zenith angle is within the range of 20-34 degrees, taking 27 degrees; when the observation zenith angle is within the range of 35-44 degrees, taking 40 degrees; when the observation zenith angle is within the range of 45-54 degrees, taking 50 degrees; when the observation zenith angle is within the range of 55-64 degrees, 60 degrees are taken; when the observation zenith angle is within the range of 65-79 degrees, taking down the odd observation zenith angles; when the observation zenith angle exceeds 90 degrees, the observation zenith angle is symmetrical about the 90-degree center;
the optimization mode of the atmospheric transmittance is that when the zenith angle is observed within the range of 0-29 degrees, 15 degrees are taken; when the observation zenith angle is within the range of 30-44 degrees, taking 37 degrees; when the observation zenith angle is within the range of 45-54 degrees, taking 50 degrees; when the observation zenith angle is within the range of 55-59 degrees, taking 57 degrees; when the observation zenith angle is within the range of 60-69 degrees, downwards taking an even number of observation zenith angles; when the observation zenith angle exceeds 90 degrees, the observation zenith angle is symmetrical about the 90-degree center;
the optimization mode of the solar radiation illumination is that when the zenith angle of the sun is within the range of 0-30 degrees, 15 degrees are taken; taking Esun _30- (Esun30-Esun40) × (X-30) when the zenith angle of the sun is within the range of 30-39 degrees; taking Esun _40- (Esun _40-Esun _ 50) × (X-40) when the zenith angle of the sun is in the range of 40-49 degrees; taking Esun _50- (Esun _50-Esun _ 60) × (X-50) when the zenith angle of the sun is within the range of 50-59 degrees; taking Esun _60- (Esun _60-Esun _ 70) (X-60) when the zenith angle of the sun is in the range of 60-69 degrees; taking Esun _70- (Esun _70-Esun _ 80) × (X-70) when the zenith angle of the sun is in the range of 70-79 degrees; wherein Esun _ x is the solar irradiance at x degrees;
the optimization mode of the earth radiation illumination is that 15 degrees are taken when the observation zenith angle is in the range of 0-29 degrees; when the observation zenith angle is within the range of 30-49 degrees, taking 40 degrees; when the observation zenith angle is within the range of 50-64 degrees, 60 degrees are taken; when the observation zenith angle is within the range of 65-79 degrees, taking down odd number zenith angles; when the observation zenith angle exceeds 90-180 degrees, the observation zenith angle is symmetrical about a 90-degree center;
the optimization mode of the sky background irradiance is that 15 degrees are taken when the observed zenith angle is in the range of 0-19 degrees; when the observation zenith angle is within the range of 20-29 degrees, 25 degrees are taken; when the observation zenith angle is within the range of 30-39 degrees, taking 35 degrees; when the observation zenith angle is within the range of 40-49 degrees, taking 45 degrees; when the observation zenith angle is within the range of 50-69 degrees, downwards taking even number zenith angles; when the zenith angle is observed to exceed 90-180 degrees, the observation is symmetrical about the 90-degree center.
2. The MODTRAN optimization-based infrared imaging simulation method according to claim 1, wherein the target infrared irradiance formula in step S3 is:
wherein L ispathRadiating luminance for an atmospheric path; l isselfIs the intrinsic radiance of the target; esunThe irradiance of the sun on the surface of the target skin; eearthThe irradiance of the earth on the surface of the target skin; eskyThe radiation illumination of sky background radiation on the surface of a target skin;is the atmospheric average transmittance of the path, dimensionless;is the reflectivity of the skin to ambient radiation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110934162.2A CN113378419B (en) | 2021-08-16 | 2021-08-16 | Infrared imaging simulation method based on MODTRAN optimization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110934162.2A CN113378419B (en) | 2021-08-16 | 2021-08-16 | Infrared imaging simulation method based on MODTRAN optimization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113378419A CN113378419A (en) | 2021-09-10 |
CN113378419B true CN113378419B (en) | 2021-11-23 |
Family
ID=77577188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110934162.2A Active CN113378419B (en) | 2021-08-16 | 2021-08-16 | Infrared imaging simulation method based on MODTRAN optimization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113378419B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080067402A (en) * | 2007-01-16 | 2008-07-21 | 연세대학교 산학협력단 | Retrieval method of aerosol optical depth using a visible channel of geostationary satellite |
CN101976275A (en) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | Airplane infrared radiation and atmospheric transmittance modeling method |
CN102628940A (en) * | 2012-04-20 | 2012-08-08 | 中国科学院遥感应用研究所 | Remote sensing image atmospheric correction method |
CN108896188A (en) * | 2018-06-06 | 2018-11-27 | 北京空间飞行器总体设计部 | A kind of high rail remote sensing satellite infrared channel atmospheric correction method of Shared aperture |
CN111191380A (en) * | 2020-01-08 | 2020-05-22 | 北京大学 | Atmospheric aerosol optical thickness estimation method and device based on measurement data of foundation spectrometer |
CN113189562A (en) * | 2021-07-02 | 2021-07-30 | 成都众享天地网络科技有限公司 | Terrain detection algorithm based on elevation |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103198314B (en) * | 2013-02-20 | 2015-11-25 | 北京农业信息技术研究中心 | Remote sensing images radiation correction method |
CN106772682B (en) * | 2016-12-31 | 2017-10-31 | 华中科技大学 | A kind of infrared radiation spectrum Simulation Analysis method of moving-target |
-
2021
- 2021-08-16 CN CN202110934162.2A patent/CN113378419B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080067402A (en) * | 2007-01-16 | 2008-07-21 | 연세대학교 산학협력단 | Retrieval method of aerosol optical depth using a visible channel of geostationary satellite |
CN101976275A (en) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | Airplane infrared radiation and atmospheric transmittance modeling method |
CN102628940A (en) * | 2012-04-20 | 2012-08-08 | 中国科学院遥感应用研究所 | Remote sensing image atmospheric correction method |
CN108896188A (en) * | 2018-06-06 | 2018-11-27 | 北京空间飞行器总体设计部 | A kind of high rail remote sensing satellite infrared channel atmospheric correction method of Shared aperture |
CN111191380A (en) * | 2020-01-08 | 2020-05-22 | 北京大学 | Atmospheric aerosol optical thickness estimation method and device based on measurement data of foundation spectrometer |
CN113189562A (en) * | 2021-07-02 | 2021-07-30 | 成都众享天地网络科技有限公司 | Terrain detection algorithm based on elevation |
Non-Patent Citations (3)
Title |
---|
"Gradient-Based Automatic Lookup Table Generator for Radiative Transfer Models";Jorge Vicent Servera等;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20180910;第57卷(第2期);第1040-1047页 * |
中红外大气辐射传输解析模型及遥感成像模拟;杨贵军等;《光谱学与光谱分析》;20090331;第29卷(第03期);第629-634页 * |
综合利用环境星CCD和红外数据反演大气气溶胶光学厚度;盛莉等;《红外与毫米波学报》;20130615;第32卷(第03期);第277-282页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113378419A (en) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bell | Exploiting a transmission grating spectrometer | |
Shanks et al. | Theoretical investigation considering manufacturing errors of a high concentrating photovoltaic of cassegrain design and its experimental validation | |
CN107993281B (en) | Method for simulating optical characteristics of space target visible light | |
CN109253976B (en) | High-spectrum real-time radiometric calibration method based on light sensing module | |
CN111077093A (en) | Method and device for quickly detecting coal gangue based on multispectral technology | |
CN108731542A (en) | Auxiliary sighting device, gun sight with debugging functions and auxiliary aim at modification method | |
CN108225552B (en) | Method for measuring concentration energy flux density distribution of heliostat field in tower-type power station | |
CN105973850B (en) | Visible light wave band atmospheric transmittance measurement method based on single-frame color image | |
CN105425394A (en) | Optical system of high-energy and high-collimated angle solar simulator | |
Cole | Modelling CPV. | |
CN113378419B (en) | Infrared imaging simulation method based on MODTRAN optimization | |
Fatehi et al. | Modeling the incidence angle dependence of photovoltaic modules in PVsyst | |
CN111770284A (en) | Backlight compensation shooting method and related device for transmission tower | |
CN104391466A (en) | Quick design method of self-adaptive optical controller | |
CN113790798B (en) | Seamless spectral imaging device, system and method for dynamic point target tracking measurement | |
CN107830928A (en) | A kind of heliograph for cirrus optical characteristic measurement | |
Miener et al. | IACT event analysis with the MAGIC telescopes using deep convolutional neural networks with CTLearn | |
CN106704898B (en) | A kind of light channel structure of space structure formula solar simulator | |
CN107271142A (en) | Real-time optical fiber positioning device and localization method based on center drilling type 4 quadrant detector | |
CN116404983A (en) | Photovoltaic module monitoring method and device, computer equipment and storage medium | |
CN110765641B (en) | Method for quantitatively calculating optical observation light pollution | |
CN103076154A (en) | Optical efficiency analysis method for light condensation and heat collection system of solar thermal power generation | |
CN106154542A (en) | A kind of quantum imaging method and quantum imaging system | |
CN211262484U (en) | Portable full-field lightning spectral imager | |
Schöttl et al. | Raytrace3D by Fraunhofer ISE |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |