CN109783862B - Large-range radiation transmission calculation and real-time rendering method - Google Patents
Large-range radiation transmission calculation and real-time rendering method Download PDFInfo
- Publication number
- CN109783862B CN109783862B CN201811526742.2A CN201811526742A CN109783862B CN 109783862 B CN109783862 B CN 109783862B CN 201811526742 A CN201811526742 A CN 201811526742A CN 109783862 B CN109783862 B CN 109783862B
- Authority
- CN
- China
- Prior art keywords
- absorption
- spectral transmittance
- radiation transmission
- section
- atmosphere
- 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
Abstract
The invention relates to a method for calculating and rendering in real time by atmospheric path radiation transmission, which comprises the following steps: establishing an atmosphere range radiation transmission calculation physical model; determining a mapping relation between the atmospheric path radiation transmission data and the texture coordinates; obtaining pre-calculated atmosphere path radiation transmission data according to the atmosphere path radiation transmission calculation physical model; and rendering the pre-calculated atmosphere path radiation transmission data according to the mapping relation between the atmosphere path radiation transmission data and the texture coordinates. The embodiment of the invention provides an atmosphere range radiation transmission calculation model considering atmosphere Mie scattering, Rayleigh scattering and atmospheric molecular absorption, which is used for calculating atmosphere infrared radiation and rendering the atmosphere range radiation in an infrared band and overcomes the defect that the existing model does not meet the atmosphere range radiation transmission calculation and rendering in the infrared band.
Description
Technical Field
The invention belongs to the field of atmospheric physics and computer graphics, and particularly relates to an atmospheric path radiation transmission calculation and real-time rendering method.
Background
In digital earth technology, the influence of atmospheric radiation transmission on the digital earth scene is very important. The effects of atmospheric scattering and atmospheric absorption on solar and ground radiation are not negligible.
At present, there are several methods for atmospheric scattering rendering in the visible light band: nishita et al, 1993 and 1996, have proposed the Nishita93 model. The Nishita93 model only considers single scattering, and the Nishita96 model realizes double scattering and only supports observation in the atmosphere. Therefore, the method is not suitable for atmospheric effect rendering under the space view angle.
O' Neal et al propose a rendering method ("in GPU" (Graphics Processing Unit) Gems 2.Addison Wesley,2005, pp.253-268.) that can implement rendering at any viewpoint, and it is difficult to meet the requirement of real-time rendering as the complexity of the scene and the calculation accuracy increase, since atmospheric calculation is performed during the rendering process.
The rendering method supporting any viewpoint, multiple scattering and rendering time complexity of O (1) is successively proposed by E.Bruneton et al in France and O.Elek et al in America as Bruneton models. The brunetton model reaches the same level in terms of pre-computation complexity, texture space complexity, computation accuracy and rendering time.
The above is the atmospheric scattering rendering method of the visible light band, wherein H2O、CO2And O3The absorption of the spectra by the equal components is negligible, while rayleigh scattering and mie scattering are the main factors of the visible light transmission characteristics. All the models do not satisfy the atmosphere range radiation transmission calculation and real-time rendering under the infrared band.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides an atmosphere range radiation transmission calculation and real-time rendering method. The technical problem to be solved by the invention is realized by the following technical scheme:
the embodiment of the invention provides a method for calculating atmospheric path radiation transmission and rendering in real time, which comprises the following steps:
establishing an atmosphere range radiation transmission calculation physical model;
determining a mapping relation between the atmospheric path radiation transmission data and the texture coordinates;
obtaining pre-calculated atmosphere path radiation transmission data according to the atmosphere path radiation transmission calculation physical model;
and rendering the pre-calculated atmosphere path radiation transmission data according to the mapping relation between the atmosphere path radiation transmission data and the texture coordinates.
In one embodiment of the invention, establishing an atmosphere range radiation transmission computational physical model comprises:
and establishing the atmospheric path radiation transmission calculation physical model according to the solar radiation, the reflection of the earth surface to the environment radiation, the atmospheric self-radiation, the atmospheric solar radiation scattering and the atmospheric earth surface radiation scattering.
In one embodiment of the present invention, the physical model of the atmospheric path radiation transmission calculation is:
wherein L isλ(x, v, s) represents the output of the atmospheric range radiation transmission computational physics model;
Lsun,λ(x, v, s) represents solar radiation;
Tλ(x,x0)I[Lλ](x0and s) represents the reflection of the surface of the earth from ambient radiation;
Tλ(x,x0)Latmo,λ(x0and s) represents atmospheric self-radiation;
In an embodiment of the present invention, obtaining pre-computed atmospheric path radiation transmission data according to the atmospheric path radiation transmission computation physical model includes:
acquiring a first spectral transmittance of the absorption gas by using Modtran;
presetting a first absorption section of the absorption gas, and obtaining a second spectral transmittance of the absorption gas according to the first absorption section;
obtaining a third spectral transmittance of the absorption gas according to the first spectral transmittance and the second spectral transmittance;
and obtaining the pre-calculated atmosphere path radiation transmission data according to the third spectral transmittance and the third absorption cross section.
In one embodiment of the present invention, presetting a first absorption cross section of the absorption gas from which a second spectral transmittance of the absorption gas is obtained includes:
presetting a first absorption section of the absorption gas according to a preset range;
and obtaining a second spectral transmittance of the absorption gas according to the first absorption section and the spectral transmittance calculation formula.
In one embodiment of the present invention, obtaining the third spectral transmittance of the absorption gas from the first spectral transmittance and the second spectral transmittance includes:
obtaining a second absorption cross section of the absorption gas according to the first spectral transmittance and the second spectral transmittance;
fitting the second absorption cross section and the first spectral transmittance according to an iterative method to obtain a third absorption cross section of the absorption gas;
and obtaining the third spectral transmittance according to the third absorption cross section.
In one embodiment of the present invention, obtaining the pre-calculated atmospheric range radiation transmission data according to the third spectral transmittance and the third absorption cross section comprises:
and calculating a physical model according to the third spectral transmittance, the third absorption cross section and the atmospheric path radiation transmission to obtain the pre-calculated atmospheric path radiation transmission data.
In one embodiment of the present invention, obtaining a second absorption cross section of the absorption gas from the first spectral transmittance and the second spectral transmittance includes:
and processing the first spectral transmittance and the second spectral transmittance according to a dichotomy to obtain the second absorption cross section.
In an embodiment of the present invention, fitting the second absorption cross section and the first spectral transmittance according to an iterative method to obtain a third absorption cross section of the absorption gas includes:
and fitting the second absorption cross section and the first spectral transmittance according to an iterative method to obtain a third absorption cross section meeting a preset threshold.
In one embodiment of the present invention, obtaining the third spectral transmittance according to the third absorption cross section includes:
and obtaining a third spectral transmittance of the absorption gas according to the third absorption cross section and the spectral transmittance calculation formula.
Compared with the prior art, the invention has the beneficial effects that:
the method combines three factors of atmospheric scattering, atmospheric absorption and atmospheric radiation and kirchhoff's law to establish an atmospheric physical model, considers the effects of absorption gas density distribution, atmospheric temperature and atmospheric pressure on atmospheric radiation, uses the GPU to perform accelerated calculation, adopts a pre-calculation mode to realize the calculation of atmospheric path radiation transmission in an infrared band, improves the rendering speed and realizes real-time rendering. The transmittance calculation and Modtran calculation result of the invention have high goodness of fit, and have great use value in digital earth and other infrared scene simulation.
Other aspects and features of the present invention will become apparent from the following detailed description, which proceeds with reference to the accompanying drawings. It is to be understood, however, that the drawings are designed solely for purposes of illustration and not as a definition of the limits of the invention, for which reference should be made to the appended claims. It should be further understood that the drawings are not necessarily drawn to scale and that, unless otherwise indicated, they are merely intended to conceptually illustrate the structures and procedures described herein.
Drawings
Fig. 1 is a flowchart of an atmosphere range radiation transmission calculation and real-time rendering method according to an embodiment of the present invention;
fig. 2a is a map of atmospheric path radiation texture coordinates provided in an embodiment of the present invention;
FIG. 2b is a graph of a spectral transmittance texture coordinate mapping according to an embodiment of the present invention;
FIG. 3 is a flow chart of fitting absorption cross sections of an absorption gas spectrum provided by an embodiment of the present invention;
FIGS. 4a to 4c are CO provided in the embodiments of the present invention2Spectral transmittance graph of (1), H2Spectral transmittance diagram of O, O3A spectral transmittance map of (a);
FIG. 5 is a schematic diagram of a world space derived from camera space according to an embodiment of the present invention;
fig. 6 is a diagram of an effect of atmospheric path radiation rendering provided by an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to specific examples, but the embodiments of the present invention are not limited thereto.
Referring to fig. 1, fig. 1 is a flowchart of an atmospheric path radiation transmission calculation and real-time rendering method according to an embodiment of the present invention, which includes the following steps:
Establishing an atmosphere range radiation transmission calculation physical model, which comprises the following steps: and establishing the atmosphere path radiation transmission calculation physical model according to the solar radiation, the scattering of the atmosphere on the solar radiation, the atmospheric self-radiation and the scattering of the atmospheric surface radiation.
Specifically, the physical model of atmosphere range radiation transmission calculation is combined with a Bruneton atmosphere model, and CO is introduced2、H2O、O3The absorption function and kirchhoff's law of the radiation transmission model of the atmosphere range are established.
The brunetton atmosphere model is given by:
L(x,v,s)＝(L0+R[L]+S[L])(x,v,s) (1)
wherein L is0For solar radiation, R [ L ]]For surface reflection, SL]Scattering for the atmosphere.
The spectral transmittance calculation formula is given by:
The scattering coefficients and phase functions of rayleigh scattering and mie scattering are as follows:
wherein N represents the number of scattering sites per unit volume;
HRis the Rayleigh scattering thickness, HR＝8km；
HMIs the thickness of the Mie scattering, HM＝1.2km；
mu is the cosine of the included angle between the incident light and the scattered light;
g is an asymmetry factor, and g is 0.7.
Further, the volumetric absorption coefficient of the absorption gas i is expressed as follows:
wherein in the formula (6),
showing the absorption cross section of a single particle of the i-th absorption gas, NiThe number of particles per unit volume of the i-th absorption gas is expressed.
Considering the influence of temperature and atmospheric pressure on the volume absorption coefficient of the gas molecules, the following correction is made on the volume absorption coefficient:
where n is an empirical constant, n (H)2O)＝0.9，n(CO2)＝0.75，n(O3) 0.4, P is the pressure at different heights, P0Is a standard atmospheric pressure, T0T is the temperature at different heights 273K.
In summary, under the condition of uniform temperature distribution, the calculation equation of the atmospheric path radiation transmission calculation physical model established by the method including solar radiation, reflection of the earth's surface to the environmental radiation, atmospheric self-radiation, scattering of the atmospheric solar radiation and scattering of the atmospheric earth's surface radiation is given by the following formula:
wherein, x, v and s are respectively an observation point position vector, a sight line direction vector in world space and a sun incident ray direction vector;
Lλ(x, v, s) represents the sameOutputting an atmospheric path radiation transmission calculation physical model, namely the radiation brightness received along the v direction at a detection point x when the ray direction is s;
Lsun,λ(x,v,s)＝T(x,x0)Lsun0,λor0, denotes solar radiation, in particular sunlight, from the atmospheric boundary X0The radiance propagated to viewpoint X;
Tλ(x,x0)I[Lλ](x0s) denotes the reflection of the ambient radiation from the surface, in particular from the reflection point X0The brightness of reflected light reaching the X point after atmospheric attenuation;
atmospheric scattering includes, among other things, single and multiple scattering of solar radiation by the atmosphere and multiple scattering of surface radiation.
When calculating the reflection of the surface to the ambient radiation, it is calculated in the same way as the brunetton atmosphere model, i.e. it is simply considered to be lambertian.
Further, the air conditioner is provided with a fan,
is the atmospheric radiance with a path length deltal at a temperature T, where,
is the volume absorption coefficient of the i-th gasPlease refer to formula (7).
Tλ(x,x0)Latmo,λ(x0S) represents atmospheric radiation, in particular from X0Atmospheric radiance at point-to-X;
The steps are expanded from the Bruneton atmospheric scattering model to the atmospheric path radiation transmission model, an atmospheric path radiation transmission model calculation equation suitable for the infrared band is established, and the formula (9) is used as a final equation for accelerating calculation of atmospheric path radiation by the GPU and is also the atmospheric path radiation transmission calculation physical model established by the embodiment of the invention.
And 2, determining a mapping relation between the atmospheric path radiation transmission data and the texture coordinates.
Referring to fig. 2a, fig. 2a is a map of an atmosphere path radiation texture coordinate according to an embodiment of the present invention.
According to the embodiment of the invention, the real-time rendering of the atmosphere path radiation is realized by means of atmosphere pre-calculation, and before the data is pre-calculated and stored in the texture, the mapping relation between the atmosphere path radiation transmission data and the texture coordinate needs to be determined firstly.
The embodiment of the invention adopts a precomputation mode to generate the spectrum transmittance texture and the large-range radiation texture, the size of the texture is generally variable, and the texture is normalized to the range of 0-1 by Graphics API (Application Programming Interface) such as OpenGL (Open Graphics Library). Therefore, when the OpenGL renders an object, the corresponding texture pixel points can be found through the texture coordinates. And then, carrying out adjacent sampling or linear sampling on the pixel points, and coloring the sampled color values onto corresponding fragments of the object.
Since the texture coordinates are only fractional numbers between 0-1. In calculating the spectral transmittance texture, it is first known what (θ, r) corresponds to each (x, y) coordinate, θ being the viewing zenith angle, in the range (0-180 °), r being the viewpoint height, in the range (0-90 km).
The atmospheric path radiation sampling parameters comprise: height r (0-90km) of observation point and zenith angle of sun
(0-180 degrees), a sight line zenith angle theta (0-180 degrees), and an included angle mu (0-180 degrees) between the sight line and the solar ray.
And 3, obtaining pre-calculated atmosphere path radiation transmission data according to the atmosphere path radiation transmission calculation physical model.
And 3.1, acquiring the first spectral transmittance of the absorption gas by utilizing Modtran.
Modtran (mode spectral resolution atmospheric transmittance), i.e., a radiation transmission algorithm and a calculation model.
In the embodiment of the invention, the absorption gas is mainly absorption gas CO2、、H2O、O3But the embodiments of the present invention are not limited thereto.
In an embodiment of the invention, the first spectral transmittance of the absorbing gas comprises: first CO2Transmittance, first H2O transmittance and first O3A transmittance.
Specifically, the first spectral transmittance of the absorption gas is calculated as follows:
setting Modtran calculation type as 'transmittance'; the atmosphere mode is '1976 American Standard atmosphere'; the atmospheric path type is a "bias layer path"; the temperature profile and the height profile of the gas pressure are "the us standard atmosphere in 1976"; the height profile of water vapor and ozone is "the U.S. standard atmosphere in 1976"; the wave band is 3-5 um; observation height of "0.5 km"; up to a height of "90 km"; zenith angle "60 °".
Through the arrangement, the first CO is output at the Modtran2Transmittance, first H2O transmittance and first O3The transmittances, and the three absorption gas transmittances were stored in the first texture.
And 3.2, pre-calculating the second spectral transmittance of the absorption gas.
Presetting a first absorption section of the absorption gas, and obtaining a second spectral transmittance of the absorption gas according to the first absorption section.
Specifically, a first absorption cross section of the absorption gas is preset according to a preset range, and a second spectral transmittance of the absorption gas is obtained according to the first absorption cross section and the spectral transmittance calculation formula.
In an embodiment of the present invention, the second spectral transmittance of the absorption gas is calculated in a pre-calculation manner, where the second spectral transmittance of the absorption gas includes: second CO2Transmittance, second H2O transmittance and second O3A transmittance.
Keeping the same water vapor and ozone height profile, atmospheric temperature and pressure height profile as those of Modtran in the step 3.1 as input parameters, and performing accelerated calculation by using OpenGL to obtain the second CO by using a formula (2) model2Transmittance, second H2O transmittance and second O3And the transmittance is stored in a second texture, the mapping relation of the texture coordinate (x, y) of the second spectral transmittance of the absorbed gas, the zenith angle theta and the observation point height r is shown in fig. 2b, and fig. 2b is a mapping relation chart of the texture coordinate of the spectral transmittance provided by the embodiment of the invention.
Further, if the spectral transmittance of each absorbing gas is calculated for the first time, the first absorption cross section of the absorbing gas of the absorbing particles needs to be manually set.
Preferably, the preset ranges are: the first absorption section of the absorption gas only needs to meet the requirement that the calculated spectral transmittance of each absorption gas is within the range of 0-1, and in order to improve the fitting speed, the spectral transmittance is close to 0.5 as far as possible.
And 3.3, pre-calculating a second absorption cross section of the absorption gas by adopting a dichotomy.
And obtaining a second absorption cross section of the absorption gas according to the first spectral transmittance and the second spectral transmittance.
Specifically, the first spectral transmittance and the second spectral transmittance are processed according to a dichotomy to obtain the second absorption cross section.
Further, a second absorption cross section of the absorption gas is firstly obtained by adopting a bisection method for the first absorption cross section of the absorption gas, the second absorption cross section of the absorption gas is only an intermediate parameter of the step, and the final second absorption cross section of the absorption gas can be determined only by carrying out a subsequent step 3.4 iteration method.
The concrete decomposition method of the dichotomy is as follows:
wherein the content of the first and second substances,
is calculating the first absorption cross section of the absorption gas used;
is a second absorption cross section of the absorption gas obtained by comparison;
is the upper limit of the interval corresponding to the second absorption section of the absorption gas;
is the lower limit of the interval corresponding to the second absorption cross section of the absorption gas.
Wherein, the formula (10) is used for solving a second absorption cross section of the absorption gas;
(11) the formula is used for solving the interval upper limit corresponding to the second absorption section of the absorption gas;
(12) the formula is used for obtaining the lower limit of the interval corresponding to the second absorption cross section of the absorption gas.
And 3.4, pre-calculating the third spectral transmittance of the absorption gas.
Fitting the second absorption cross section and the first spectral transmittance according to an iterative method to obtain a third absorption cross section meeting a preset threshold, and obtaining a third spectral transmittance of the absorption gas according to the third absorption cross section and the spectral transmittance calculation formula.
Preferably, the preset threshold value means that the error range of the absolute value of the first spectral transmittance minus the second spectral transmittance is 0.01-0.001.
Referring to fig. 3, fig. 3 is a flowchart illustrating fitting of a spectral absorption cross section of an absorption gas according to an embodiment of the present invention.
And (3.2) and step 3.3 are circulated, a third absorption cross section of the absorption gas is worked out by adopting an iterative method to fit, specifically, the second absorption cross section in the step 3.3 is substituted into a spectral transmittance calculation formula to obtain a second spectral transmittance of the absorption gas in the step 3.2, the second spectral transmittance of the absorption gas is compared with the first spectral transmittance of the absorption gas, if the first spectral transmittance accuracy requirement of the absorption gas is not met, the steps 3.3 and 3.2 are continuously circulated until a preset threshold value is met, and the second absorption cross section meeting the first spectral transmittance accuracy requirement of the absorption gas is taken as the third absorption cross section of the absorption gas.
Substituting the third absorption cross section of the absorption gas into a formula (2), and calculating the third spectral transmittance of the absorption gas through the atmospheric temperature and pressure, Rayleigh scattering, Mie scattering and the third absorption cross section of the absorption gas.
Specifically, please refer to fig. 4a to 4c, and fig. 4a to 4c respectively illustrate CO provided by embodiments of the present invention2Spectral transmittance graph of (1), H2Spectral transmittance diagram of O, O3Spectrum transmittance map of (1).
And 3.5, pre-calculating the atmospheric path radiation transmission data.
And obtaining a third spectral transmittance of the absorption gas according to the third absorption cross section and the spectral transmittance calculation formula, and obtaining the pre-calculated atmosphere range radiation transmission data according to the third spectral transmittance, the third absorption cross section and the atmosphere range radiation transmission calculation physical model. Specifically, the third spectral transmittance of the absorption gas is obtained by using the spectral transmittance calculation formula with the third absorption cross section as an input, the large-range radiation transmission data is calculated by using the formula (9) with the third spectral transmittance of the absorption gas as an input, and the large-range radiation transmission data is stored in the third texture.
And 4, rendering the pre-calculated atmosphere path radiation transmission data.
And rendering the pre-calculated atmosphere path radiation transmission data according to the mapping relation between the atmosphere path radiation transmission data and the texture coordinates.
Referring to fig. 5, fig. 5 is a schematic diagram of a world space derived from a camera space according to an embodiment of the present invention.
The embodiment of the invention renders the sampling of the atmospheric path radiation transmission parameter on a plane. The method comprises the following specific steps:
and 4.1, drawing a square with vertex coordinates of (-1, -1), (1, 1) and (-1, 1) by using OpenGL, wherein the square just covers the whole viewport in the OpenGL standard space.
And 4.2, inverting the coordinate P corresponding to each coloring fragment on the square plane to convert the coordinate P into a visual space. Obtaining the direction vector of the point P corresponding to the coloring fragment in the visual space
Wherein, P ═ (x, y,0,1)TAnd x and y are coordinate points corresponding to the square plane.
Step 4.3, viewing the direction vector in the space
The direction vector in the world space is obtained by inversion conversion to the world coordinate system
According to the position P of the sun under the world coordinate systemsun(x, y, z) and earth position Pearth(x, y, z) determining the direction vector of the incident solar ray
In particular, the amount of the solvent to be used,
step 4.4, according to the position P of the observation point under the world coordinate systemeye(x, y, z) and earth position Pearth(x, y, z) determining a normal vector of the observation point with respect to the earth's surface
And an observation height r, specifically,
r＝|Peye(x,y,z)-Pearth(x,y,z)|-rearthwherein r isearthRepresenting the radius of the earth.
Step 4.5, according to the direction vector under the world space
And the direction vector of the incident ray of the sun
The cosine cos mu of the angle between the line of sight and the light is found, specifically,
according to the direction vector under the world space
And the normal vector of the observation point relative to the earth's surface
The cosine of the zenith angle of the line of sight cos θ is found, specifically,
according to the direction vector of the incident ray of the sun
And the normal vector of the observation point relative to the earth's surface
Finding the zenith angle cosine of the sun
In particular, the amount of the solvent to be used,
step 4.6, for r, cos mu, cos theta,
And the four parameters are used for rendering the pre-calculated atmosphere pass radiation transmission data to a plane according to the mapping relation between the atmosphere pass radiation transmission data and the texture coordinates, so that real-time rendering is realized. Referring to fig. 6, fig. 6 is a diagram of an atmosphere range radiation rendering effect according to an embodiment of the present invention.
Compared with the prior art, the invention has the beneficial effects that:
the method combines three factors of atmospheric scattering, atmospheric absorption and atmospheric radiation and kirchhoff's law to establish an atmospheric physical model, considers the effects of absorption gas density distribution, atmospheric temperature and atmospheric pressure on atmospheric radiation, uses the GPU to perform accelerated calculation, adopts a pre-calculation mode to realize the calculation of atmospheric path radiation transmission in an infrared band, improves the rendering speed and realizes real-time rendering. The transmittance calculation and Modtran calculation result of the invention have high goodness of fit, and have great use value in digital earth and other infrared scene simulation.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For atmospheric calculations in the visible range, the model is only a subset of the entire model, so the model is equally applicable to the visible range. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention.
Claims (3)
1. A method for calculating and rendering in real time by radiation transmission in a large air range is characterized by comprising the following steps:
step 1, establishing an atmosphere range radiation transmission calculation physical model, wherein the atmosphere range radiation transmission calculation physical model is as follows:
wherein L isλ(x, v, s) represents the output of the atmospheric range radiation transmission computational physics model;
Lsun,λ(x, v, s) represents solar radiation;
Tλ(x,x0)I[Lλ](x0and s) represents the reflection of the surface of the earth from ambient radiation;
Tλ(x,x0)Latmo,λ(x0and s) represents atmospheric self-radiation;
step 2, determining a mapping relation between the atmospheric path radiation transmission data and texture coordinates;
step 3, obtaining pre-calculated atmosphere path radiation transmission data according to the atmosphere path radiation transmission calculation physical model;
the step 3 comprises the following steps:
step 3.1, acquiring a first spectral transmittance of the absorption gas by utilizing Modtran, and storing the first spectral transmittance into a first texture, wherein the first spectral transmittance comprises first CO2Transmittance, first H2O transmittance and first O3A transmittance;
step 3.2, presetting a first absorption section of the absorption gas, obtaining a second spectral transmittance of the absorption gas according to the first absorption section, and storing the second spectral transmittance into a second texture, wherein the second spectral transmittance comprises second CO2Transmittance, second H2O transmittance and second O3A transmittance;
step 3.3, calculating a second absorption cross section of the absorption gas by adopting a bisection method for the first absorption cross section of the absorption gas;
step 3.4, substituting the second absorption cross section of the absorption gas in the step 3.3 into a spectral transmittance calculation formula to obtain a second spectral transmittance of the absorption gas in the step 3.2, comparing the first spectral transmittance of the absorption gas with the second spectral transmittance of the absorption gas, if the first spectral transmittance accuracy requirement of the absorption gas is not met, continuing to circulate the step 3.2 and the step 3.3 until a preset threshold value is met, stopping circulation, and taking the second absorption cross section meeting the first spectral transmittance accuracy requirement of the absorption gas as a third absorption cross section of the absorption gas, wherein the preset threshold value means that the error range of subtracting the absolute value of the second spectral transmittance from the first spectral transmittance is 0.01-0.001;
step 3.5, obtaining a third spectral transmittance according to the third absorption cross section;
step 3.6, obtaining the pre-calculated atmosphere path radiation transmission data according to the third spectral transmittance, the third absorption cross section and the atmosphere path radiation transmission calculation physical model, and storing the pre-calculated atmosphere path radiation transmission data into a third texture;
and 4, rendering the pre-calculated atmosphere path radiation transmission data according to the mapping relation between the atmosphere path radiation transmission data and the texture coordinates.
2. The method of claim 1, wherein presetting a first absorption cross-section of the absorption gas from which a second spectral transmittance of the absorption gas is obtained comprises:
presetting a first absorption section of the absorption gas according to a preset range;
and obtaining a second spectral transmittance of the absorption gas according to the first absorption section and the spectral transmittance calculation formula.
3. The method of claim 1, wherein obtaining the third spectral transmittance from the third absorption cross-section comprises:
and obtaining a third spectral transmittance of the absorption gas according to the third absorption cross section and the spectral transmittance calculation formula.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811526742.2A CN109783862B (en) | 2018-12-13 | 2018-12-13 | Large-range radiation transmission calculation and real-time rendering method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811526742.2A CN109783862B (en) | 2018-12-13 | 2018-12-13 | Large-range radiation transmission calculation and real-time rendering method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109783862A CN109783862A (en) | 2019-05-21 |
| CN109783862B true CN109783862B (en) | 2021-04-13 |
Family
ID=66495939
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811526742.2A Active CN109783862B (en) | 2018-12-13 | 2018-12-13 | Large-range radiation transmission calculation and real-time rendering method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109783862B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113592999B (en) * | 2021-08-05 | 2022-10-28 | 广州益聚未来网络科技有限公司 | Rendering method of virtual luminous body and related equipment |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8190403B2 (en) * | 2007-06-26 | 2012-05-29 | Microsoft Corporation | Real-time rendering of light-scattering media |
| CN101915914B (en) * | 2010-07-30 | 2012-10-24 | 南京信息工程大学 | Lookup table based pixel-by-pixel atmospheric correction method of remote sensing images |
| CN103018736B (en) * | 2012-12-03 | 2014-11-26 | 北京航空航天大学 | Satellite-borne remote sensor radiation calibration method based on atmospheric parameter remote sensing retrieval |
| CN103728612B (en) * | 2013-12-23 | 2017-02-15 | 中北大学 | Passive distance measuring method based on target infrared radiation spectrum and band model |
| CN103792538B (en) * | 2014-02-24 | 2016-09-28 | 北京航空航天大学 | A kind of atmospheric profile retrieval method based on ground EO-1 hyperion microwave radiometer |
| CN104880702A (en) * | 2014-10-31 | 2015-09-02 | 航天恒星科技有限公司 | Method and device for on-orbit absolute radiation calibration |
| CN104406686B (en) * | 2014-12-10 | 2016-03-23 | 西北师范大学 | Sun shortwave incident radiation evaluation method under MODEL OVER COMPLEX TOPOGRAPHY |
| CN105158183A (en) * | 2015-09-08 | 2015-12-16 | 安徽理工大学 | Method for improving sensitivity for detecting NO3 in atmosphere by deducting vapor interference |
| CN105403201B (en) * | 2015-10-20 | 2016-09-21 | 浙江农林大学 | A kind of remote sensing images atmospheric path radiation acquisition methods based on pixel analysis |
-
2018
- 2018-12-13 CN CN201811526742.2A patent/CN109783862B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109783862A (en) | 2019-05-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111210396A (en) | Multispectral polarization image defogging method combined with sky light polarization model | |
| CN114544452B (en) | Multi-angle polarized water color remote sensor satellite atmosphere correction method | |
| CN115630256B (en) | Multi-angle polarized water color satellite atmospheric correction method based on dark pixel assumption | |
| CN109783862B (en) | Large-range radiation transmission calculation and real-time rendering method | |
| CN116228984A (en) | Volumetric cloud modeling and rendering method based on meteorological data | |
| CN114663786A (en) | Forest stand radiation flux calculation method based on point cloud data and computer graphics | |
| CN107301633B (en) | Simulation method for remote sensing imaging under cloud and fog interference | |
| CN107132545B (en) | The analysis method that particulate matter in aerosol influences the quality of remotely sensed image | |
| CN114925553B (en) | Infrared image simulation method based on theoretical/semi-empirical method | |
| Inakage | An illumination model for atmospheric environments | |
| CN108872093B (en) | O based on passive remote sensing4Absorption correction and aerosol extinction profile inversion method | |
| Richtsmeier et al. | Full spectrum cloudy scene simulation | |
| CN111859685B (en) | Rapid generation method of infrared view of ship body | |
| CN111860328B (en) | Biomass estimation method based on bidirectional reflection function and forest scene illumination effect modeling | |
| Preetham | Modeling skylight and aerial perspective | |
| CN114898040A (en) | Real-time infrared imaging simulation method and system based on hardware pipeline acceleration | |
| Wu et al. | Quantitative atmospheric rendering for real-time infrared scene simulation | |
| Morales et al. | Real-time rendering of aerial perspective effect based on turbidity estimation | |
| Yuan et al. | Near-Surface Atmospheric Scattering Rendering Method | |
| CN113656928B (en) | Rapid tail flame infrared simulation method based on single-parameter control | |
| Kahraman et al. | Influence of fog on computer vision algorithms | |
| CN114169265A (en) | Sea surface dynamic infrared simulation method | |
| CN116030179B (en) | Data processing method, device, computer equipment and storage medium | |
| Guo et al. | Efficient Simulation of Light Scattering Effects in the Atmosphere. | |
| CN116894905A (en) | Real-time water wave simulation method and device based on Niagala particle system |
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 |