CN115098827A - Infrared scene numerical value calculation method in snow accumulation environment - Google Patents

Infrared scene numerical value calculation method in snow accumulation environment Download PDF

Info

Publication number
CN115098827A
CN115098827A CN202211030738.3A CN202211030738A CN115098827A CN 115098827 A CN115098827 A CN 115098827A CN 202211030738 A CN202211030738 A CN 202211030738A CN 115098827 A CN115098827 A CN 115098827A
Authority
CN
China
Prior art keywords
scene
snow
calculating
infrared
speed
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.)
Granted
Application number
CN202211030738.3A
Other languages
Chinese (zh)
Other versions
CN115098827B (en
Inventor
陈世国
黑圆圆
莫冬腊
陶克勇
薛瑞文
赖康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN202211030738.3A priority Critical patent/CN115098827B/en
Publication of CN115098827A publication Critical patent/CN115098827A/en
Application granted granted Critical
Publication of CN115098827B publication Critical patent/CN115098827B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/48Thermography; Techniques using wholly visual means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V8/00Prospecting or detecting by optical means
    • G01V8/10Detecting, e.g. by using light barriers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

The invention discloses a method for calculating an infrared scene numerical value in an accumulated snow environment, which comprises the following steps of: calculating the landing speed vector of the snowflakes; calculating the snow cover thickness in the scene based on the speed vector of the snowfall; calculating the infrared radiation brightness of the scene according to the snow cover thickness in the scene; and generating an infrared simulation image according to the infrared radiation brightness of the scene. The method aims to accurately calculate the snow cover condition of the scene surface through the ray tracing technology, so that the infrared scene radiation characteristic can be accurately calculated.

Description

Infrared scene numerical value calculation method in snow accumulation environment
Technical Field
The invention relates to the technical field of infrared scene simulation, in particular to an infrared scene numerical calculation method in an accumulated snow environment.
Background
With the wide application of infrared detection technology and infrared imaging guided weapons in informatization battlefields, modern wars increasingly pay more attention to detection, identification and attack through infrared characteristics of targets and backgrounds. In order to research the infrared radiation characteristics of snow, researchers provide several methods for calculating the infrared radiation characteristics in the snow environment, one method is to synthesize a ground target into an infrared image of the snow, but the method is only simple image superposition; one is based on a calculation model of the snow surface temperature, which in turn calculates the infrared radiation. However, these methods have two drawbacks: 1) the infrared radiation characteristic of the whole scene cannot be calculated comprehensively, so that distortion is caused; 2) the snow cover condition in the scene cannot be accurately calculated, and infrared radiation errors are large. In the article of 'infrared radiation calculation and application of snow surface' published journal paper by researchers of mechanical engineering department of armored forces engineering academy of sciences, a temperature calculation model of snow surface is established based on a thermal equilibrium equation, infrared radiation brightness is calculated according to temperature, and then infrared image simulation is carried out. The method is relatively accurate in calculating the temperature of the snow, but cannot accurately calculate the snow coverage condition in the scene, so that the calculated infrared radiation error is large.
The infrared scene numerical calculation method under the snow accumulation environment calculates the snow covered condition of the scene surface through the environment parameters, selects different material attributes, further calculates the infrared radiation brightness of the scene under the snow accumulation environment, and generates the infrared simulation image. The infrared scene calculation under the snow accumulation environment has great significance for research on the aspects of infrared target characteristic simulation, infrared camouflage/stealth effect evaluation, infrared target identification and the like. In order to accurately calculate the infrared radiation characteristic of a scene under a snow environment, researchers provide a calculation model of the surface temperature of the snow, and a basis is provided for infrared scene simulation, however, the method has certain limitations: the snow cover conditions of different parts in a scene under the influence of factors such as wind speed and wind direction cannot be accurately considered, so that the calculation error is large, and the description of a real battlefield environment is not met.
Disclosure of Invention
Aiming at the problem that the snow cover situation in a scene cannot be accurately calculated in the prior art, the invention provides the infrared scene numerical value calculation method in the snow environment.
The invention discloses a method for calculating an infrared scene numerical value in an accumulated snow environment, which comprises the following steps of:
step 1: calculating the landing speed vector of the snowflakes;
step 2: calculating the snow cover thickness in the scene based on the speed vector of the snowfall;
and 3, step 3: calculating scene infrared radiation brightness according to the snow cover thickness in the scene;
and 4, step 4: and generating an infrared simulation image according to the infrared radiation brightness of the scene.
Further, the step 1 comprises:
step 11: obtaining a wind speed vector under a scene coordinate system according to the direction and the size of the wind speed in the scene;
step 12: calculating a speed vector of the snowflakes falling to the ground at the current moment according to the gravity acceleration, the air resistance coefficient, the diameters of the snowflakes, the mass of the snowflakes and the time interval;
step 13: and (4) calculating the speed vector of snowfall at the current moment in the absence of wind, and combining the speed vector of the snowfall in the scene synthesized by the wind speed under the scene coordinate system in the step 11.
Further, the step 11 includes:
assume that the direction of the wind speed in the scene is
Figure 100002_DEST_PATH_IMAGE001
The wind speed is
Figure 230066DEST_PATH_IMAGE002
Calculating a wind speed vector under the scene coordinate system by the following formula:
Figure 100002_DEST_PATH_IMAGE003
Figure 703903DEST_PATH_IMAGE004
Figure 100002_DEST_PATH_IMAGE005
Figure 315013DEST_PATH_IMAGE006
wherein, the first and the second end of the pipe are connected with each other,
Figure 100002_DEST_PATH_IMAGE007
is the wind speed vector in the scene coordinate system,
Figure 216104DEST_PATH_IMAGE008
Figure 100002_DEST_PATH_IMAGE009
Figure 811034DEST_PATH_IMAGE010
the wind speed in the x-axis direction, the wind speed in the y-axis direction and the wind speed in the z-axis direction are respectively.
Further, the step 12 includes:
assuming a gravitational acceleration of
Figure 100002_DEST_PATH_IMAGE011
The air resistance is a number
Figure 820054DEST_PATH_IMAGE012
The diameter of the snowflake is
Figure 100002_DEST_PATH_IMAGE013
The snowflakes have the mass of
Figure 285670DEST_PATH_IMAGE014
At a time interval of
Figure 100002_DEST_PATH_IMAGE015
Calculating the speed vector of the snowflake falling to the ground at the current moment by using the following formula:
Figure 357662DEST_PATH_IMAGE016
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE017
is composed of
Figure 174309DEST_PATH_IMAGE018
The falling speed of the snow at the moment,
Figure 100002_DEST_PATH_IMAGE019
in order to be the initial speed of the vehicle,
Figure 724370DEST_PATH_IMAGE020
is subject to
Figure 100002_DEST_PATH_IMAGE021
Are evenly distributed among them.
Further, the step 13 includes:
calculating the speed of snowfall at the current moment in absence of wind by the following formula:
Figure 310072DEST_PATH_IMAGE022
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE023
the speed of snowfall at the current moment when no wind exists;
the speed vector of the snowfall in the scene is obtained through the following formula:
Figure 818545DEST_PATH_IMAGE024
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE025
the velocity vector when a snowflake lands in a scene is unitized.
Further, the step 2 comprises:
step 21: on the basis of the regular icosahedron, each edge of 30 edges of the regular icosahedron is trisected, and a trisection point of each edge is taken to intercept an angle to obtain an truncated icosahedron with 60 vertexes and 32 faces;
step 22: respectively connecting the body center point of the truncated icosahedron with 60 vertexes and 32 surface centers to obtain 92 direction vectors, unitizing the direction vectors and storing the direction vectors in a vector array with the length of 92
Figure 870290DEST_PATH_IMAGE026
Step 23: the inverse direction of snowflake landing is regarded as the direction of tracing ray, and ray tracing technology and array are utilized
Figure 942151DEST_PATH_IMAGE026
Setting a threshold value, and calculating a snow factor function in a scene;
step 24: and calculating the snow thickness of the surface unit of the intersection surface of the ray and the scene according to the average snowing thickness.
Further, the step 23 includes:
recording the negative vector of the speed direction when the snowflakes land and the second of 92 direction vectors formed by the truncated icosahedron
Figure 100002_DEST_PATH_IMAGE027
Cosine of the angle of each vector is
Figure 647939DEST_PATH_IMAGE028
Solid angle formed by stretching in the direction of speed of snow falling
Figure 100002_DEST_PATH_IMAGE029
The threshold value is
Figure 796155DEST_PATH_IMAGE030
Within this range is
Figure 100002_DEST_PATH_IMAGE031
The number of direction vectors of
Figure 587393DEST_PATH_IMAGE032
Then snow factor function in scene
Figure 100002_DEST_PATH_IMAGE033
The calculation is as follows:
Figure 479257DEST_PATH_IMAGE034
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE035
unit normal vector and scene coordinate system of snowflake landing surface
Figure 39551DEST_PATH_IMAGE036
Cosine of the angle formed by the axes, if
Figure 100002_DEST_PATH_IMAGE037
Then, then
Figure 358668DEST_PATH_IMAGE038
Further, the step 24 includes:
average snow thickness in a scene
Figure 100002_DEST_PATH_IMAGE039
The accumulated snow thickness of the surface unit of the intersection surface of the ray and the scene is given by the configuration file
Figure 674022DEST_PATH_IMAGE040
Further, the step 3 comprises:
step 31: judging the infrared property of the material of the scene surface according to the accumulated snow thickness of the surface unit;
step 32: and calculating the infrared radiation brightness of the infrared attribute of the material on the scene surface according to the Planck formula.
Further, the step 31 includes:
thickness of accumulated snow on dough units
Figure 100002_DEST_PATH_IMAGE041
If the unit is snow-free, the infrared property of the material on the surface of the scene is the original infrared property of the material of the unit;
thickness of accumulated snow on dough units
Figure 618845DEST_PATH_IMAGE042
And the infrared property of the material on the surface of the scene is the infrared property of the snow material.
Due to the adoption of the technical scheme, the invention has the following advantages: the method for calculating the snow cover condition of the scene surface comprises the steps of calculating the snow falling direction, constructing a snow factor function by using the spherical normal of 92 directions of a truncated icosahedron, calculating the snow thickness of each surface element by using a ray tracing technology to obtain more accurate snow cover condition of the infrared scene, selecting infrared material attributes of different materials according to snow coverage to calculate the infrared radiation value of a pixel, and calculating the infrared radiation value of the pixel according to the attribute of the pixel
Figure DEST_PATH_IMAGE043
The infrared simulation image is output according to the resolution ratio, and the accurate and high-resolution infrared radiation simulation image in the snow accumulation scene is obtained. The method can effectively reduce the infrared simulation error caused by the snow cover condition and improve the infrared image precision.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the description below are only some embodiments described in the embodiments of the present invention, and it is obvious for those skilled in the art that other drawings may be obtained according to the drawings.
Fig. 1 is a schematic flow chart of a method for calculating an infrared scene value in a snow accumulation environment according to an embodiment of the present invention.
Detailed Description
The invention will be further described with reference to the drawings and examples, it being understood that the examples described are only some of the examples and are not intended to be exhaustive. All other embodiments available to those of ordinary skill in the art are intended to be within the scope of the embodiments of the present invention.
Referring to fig. 1, the present invention provides an embodiment of a method for calculating an infrared scene value in a snow accumulation environment, which includes the following steps:
s1: calculating the landing speed vector of the snowflakes;
s2: calculating the snow cover thickness in the scene based on the speed vector of the snowfall;
s3: calculating scene infrared radiation brightness according to the snow cover thickness in the scene;
s4: and generating an infrared simulation image according to the infrared radiation brightness of the scene.
In this embodiment, S1 includes:
s11: obtaining a wind speed vector under a scene coordinate system according to the direction and the size of the wind speed in the scene;
s12: calculating the speed vector of the snowflakes falling to the ground at the current moment according to the gravity acceleration, the air resistance coefficient, the diameters of the snowflakes, the mass of the snowflakes and the time interval;
s13: and calculating the speed vector of snowfall at the current moment in the absence of wind, and combining the speed vector of the snowfall in the scene synthesized by the wind speed in the scene coordinate system in S11.
In this embodiment, S11 includes:
assume that the direction of the wind speed in the scene is
Figure 784378DEST_PATH_IMAGE001
The wind speed is
Figure 258084DEST_PATH_IMAGE044
Calculating a wind speed vector under a scene coordinate system by the following formula:
Figure 227177DEST_PATH_IMAGE003
Figure 178953DEST_PATH_IMAGE004
Figure 933413DEST_PATH_IMAGE005
Figure 312442DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure 768831DEST_PATH_IMAGE007
is the wind speed vector in the scene coordinate system,
Figure 524298DEST_PATH_IMAGE008
Figure 133265DEST_PATH_IMAGE009
Figure 948774DEST_PATH_IMAGE010
the wind speed in the x-axis direction, the wind speed in the y-axis direction and the wind speed in the z-axis direction are respectively.
In this embodiment, S12 includes:
assuming a gravitational acceleration of
Figure 892459DEST_PATH_IMAGE011
The air resistance is a number
Figure 186037DEST_PATH_IMAGE012
Diameter of snow flakeIs composed of
Figure 912160DEST_PATH_IMAGE013
The snowflakes have the mass of
Figure 898571DEST_PATH_IMAGE014
At a time interval of
Figure 329552DEST_PATH_IMAGE015
Calculating the speed vector of the snowflake falling to the ground at the current moment by using the following formula:
Figure 426821DEST_PATH_IMAGE016
wherein the content of the first and second substances,
Figure 744801DEST_PATH_IMAGE017
is composed of
Figure 902113DEST_PATH_IMAGE018
The speed of fall of the snow at the moment,
Figure 820391DEST_PATH_IMAGE019
in order to be the initial speed of the vehicle,
Figure 455771DEST_PATH_IMAGE020
is subject to
Figure 893837DEST_PATH_IMAGE021
Are evenly distributed in between.
In this embodiment, S13 includes:
calculating the speed of snowfall at the current moment in absence of wind by the following formula:
Figure 222050DEST_PATH_IMAGE022
wherein the content of the first and second substances,
Figure 627624DEST_PATH_IMAGE023
when there is no windThe speed of snowflake carving landing;
the speed vector of the snowfall in the scene is obtained through the following formula:
Figure 66695DEST_PATH_IMAGE024
wherein, the first and the second end of the pipe are connected with each other,
Figure 359267DEST_PATH_IMAGE025
the velocity vector when a snowflake lands in a scene is unitized.
In this embodiment, S2 includes:
s21: on the basis of the regular icosahedron, each edge of 30 edges of the regular icosahedron is trisected, and a trisection point of each edge is taken to intercept an angle to obtain an truncated icosahedron with 60 vertexes and 32 faces;
s22: respectively connecting the body center point of the truncated icosahedron with 60 vertexes and 32 surface centers to obtain 92 direction vectors, unitizing the direction vectors and storing the direction vectors in a vector array with the length of 92
Figure 858382DEST_PATH_IMAGE026
S23: the inverse direction of snowflake landing is regarded as the direction of tracing ray, and ray tracing technology and array are utilized
Figure 485672DEST_PATH_IMAGE026
Setting a threshold value, and calculating a snow factor function in a scene;
s24: and calculating the snow thickness of the surface unit of the intersection surface of the ray and the scene according to the average snowing thickness.
In this embodiment, S23 includes:
the negative vector of the speed direction when the snowflake lands and the second of 92 direction vectors formed by the truncated icosahedron
Figure 462856DEST_PATH_IMAGE027
Cosine of the included angle of the individual vectors is
Figure 138163DEST_PATH_IMAGE028
Setting a solid angle formed by stretching the snow in the direction of the speed of falling on the ground
Figure 808179DEST_PATH_IMAGE029
The threshold value is
Figure 657186DEST_PATH_IMAGE030
Within this range is
Figure 438060DEST_PATH_IMAGE031
The number of direction vectors of
Figure 705225DEST_PATH_IMAGE032
Then snow factor function in scene
Figure 811721DEST_PATH_IMAGE033
The calculation is as follows:
Figure 413604DEST_PATH_IMAGE034
wherein the content of the first and second substances,
Figure 748901DEST_PATH_IMAGE035
unit normal vector and scene coordinate system of snowflake landing surface
Figure 119840DEST_PATH_IMAGE036
Cosine of the angle formed by the axes, if
Figure 131658DEST_PATH_IMAGE037
Then, then
Figure 955257DEST_PATH_IMAGE038
In this embodiment, S24 includes:
average snow thickness in scene
Figure 94246DEST_PATH_IMAGE039
Given by the configuration file, the intersection of the ray with the sceneThe thickness of the accumulated snow of the surface unit
Figure 319691DEST_PATH_IMAGE040
In this embodiment, S3 includes:
s31: judging the infrared property of the material of the scene surface according to the accumulated snow thickness of the surface unit;
s32: and calculating the infrared radiation brightness of the infrared attribute of the material on the scene surface according to the Planck formula.
In this embodiment, S31 includes:
thickness of accumulated snow on dough units
Figure 502411DEST_PATH_IMAGE041
If the unit is snow-free, the infrared property of the material on the surface of the scene is the original infrared property of the material of the unit;
thickness of accumulated snow on dough units
Figure 78885DEST_PATH_IMAGE042
And the infrared property of the material on the surface of the scene is the infrared property of the snow material.
The method for calculating the snow cover condition of the scene surface comprises the steps of calculating the snow falling direction, constructing a snow factor function by using the spherical normal of 92 directions of a truncated icosahedron, calculating the snow thickness of each surface element by using a ray tracing technology to obtain more accurate snow cover condition of the infrared scene, selecting infrared material attributes of different materials according to snow coverage to calculate the infrared radiation value of a pixel, and calculating the infrared radiation value of the pixel according to the attribute of the pixel
Figure DEST_PATH_IMAGE045
The infrared simulation image is output according to the resolution ratio, and the accurate and high-resolution infrared radiation simulation image in the snow accumulation scene is obtained. Experimental results show that the method can effectively reduce infrared simulation errors caused by snow cover conditions and improve the infrared image precision.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims (10)

1. An infrared scene numerical calculation method under the snow accumulation environment is characterized by comprising the following steps:
step 1: calculating the landing speed vector of the snowflakes;
and 2, step: calculating the snow cover thickness in the scene based on the speed vector of the snowfall;
and step 3: calculating scene infrared radiation brightness according to the snow cover thickness in the scene;
and 4, step 4: and generating an infrared simulation image according to the infrared radiation brightness of the scene.
2. The method of claim 1, wherein step 1 comprises:
step 11: obtaining a wind speed vector under a scene coordinate system according to the direction and the size of the wind speed in the scene;
step 12: calculating the speed vector of the snowflakes falling to the ground at the current moment according to the gravity acceleration, the air resistance coefficient, the diameters of the snowflakes, the mass of the snowflakes and the time interval;
step 13: and (5) calculating the speed vector of snowfall at the current moment in the absence of wind, and combining the speed vector of snowfall in the scene synthesized by the wind speed under the scene coordinate system in the step 11.
3. The method according to claim 2, wherein the step 11 comprises:
assume that the direction of the wind speed in the scene is
Figure DEST_PATH_IMAGE001
The wind speed is
Figure 146717DEST_PATH_IMAGE002
Calculating a wind speed vector under a scene coordinate system by the following formula:
Figure DEST_PATH_IMAGE003
Figure 756821DEST_PATH_IMAGE004
Figure DEST_PATH_IMAGE005
Figure 693553DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE007
is the wind speed vector in the scene coordinate system,
Figure 680095DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE009
Figure 753093DEST_PATH_IMAGE010
the wind speed in the x-axis direction, the wind speed in the y-axis direction and the wind speed in the z-axis direction are respectively.
4. The method of claim 3, wherein step 12 comprises:
assuming gravitational acceleration of
Figure DEST_PATH_IMAGE011
The air resistance is a number
Figure 140956DEST_PATH_IMAGE012
The diameter of the snowflake is
Figure DEST_PATH_IMAGE013
The snowflakes have the mass
Figure 299404DEST_PATH_IMAGE014
At a time interval of
Figure DEST_PATH_IMAGE015
Calculating the speed vector of the snowflake falling to the ground at the current moment by using the following formula:
Figure 824058DEST_PATH_IMAGE016
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE017
is composed of
Figure 17142DEST_PATH_IMAGE018
The speed of fall of the snow at the moment,
Figure DEST_PATH_IMAGE019
in order to be the initial speed of the vehicle,
Figure 969048DEST_PATH_IMAGE020
is subject to
Figure DEST_PATH_IMAGE021
Are evenly distributed in between.
5. The method according to claim 4, wherein the step 13 comprises:
calculating the speed of snowfall at the current moment in absence of wind by the following formula:
Figure 614793DEST_PATH_IMAGE022
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE023
the speed of snowfall at the current moment in absence of wind;
the speed vector of the snowfall in the scene is obtained through the following formula:
Figure 940208DEST_PATH_IMAGE024
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE025
the velocity vector when a snowflake lands in a scene is unitized.
6. The method of claim 5, wherein the step 2 comprises:
step 21: on the basis of the regular icosahedron, each edge of 30 edges of the regular icosahedron is trisected, and a trisection point of each edge is taken to intercept an angle to obtain an truncated icosahedron with 60 vertexes and 32 faces;
step 22: respectively connecting the body center point of the truncated icosahedron with 60 vertexes and 32 surface centers to obtain 92 direction vectors, unitizing the direction vectors and storing the direction vectors in a vector array with the length of 92
Figure 987798DEST_PATH_IMAGE026
Step 23: the inverse direction of snowflake landing is regarded as the direction of tracing ray, and ray tracing technology and array are utilized
Figure 579448DEST_PATH_IMAGE026
Setting a threshold value, and calculating a snow factor function in a scene;
step 24: and calculating the snow thickness of the surface unit of the intersection surface of the ray and the scene according to the average snowing thickness.
7. The method of claim 6, wherein the step 23 comprises:
the negative vector of the speed direction when the snowflake lands and the second of 92 direction vectors formed by the truncated icosahedron
Figure DEST_PATH_IMAGE027
Cosine of the angle of each vector is
Figure 978068DEST_PATH_IMAGE028
Setting a solid angle formed by stretching the snow in the direction of the speed of falling on the ground
Figure DEST_PATH_IMAGE029
The threshold value is
Figure 578945DEST_PATH_IMAGE030
Within this range, the
Figure DEST_PATH_IMAGE031
The number of direction vectors of
Figure 746621DEST_PATH_IMAGE032
Then snow factor function in scene
Figure DEST_PATH_IMAGE033
The calculation is as follows:
Figure 40330DEST_PATH_IMAGE034
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE035
unit normal vector and scene coordinate system of snowflake landing surface
Figure 129509DEST_PATH_IMAGE036
Cosine of the angle formed by the axes, if
Figure DEST_PATH_IMAGE037
Then, then
Figure 62305DEST_PATH_IMAGE038
8. The method of claim 7, wherein step 24 comprises:
average snow thickness in scene
Figure DEST_PATH_IMAGE039
The accumulated snow thickness of the surface unit of the intersection surface of the ray and the scene is given by the configuration file
Figure 818909DEST_PATH_IMAGE040
9. The method of claim 1, wherein step 3 comprises:
step 31: judging the infrared property of the material of the scene surface according to the accumulated snow thickness of the surface unit;
step 32: and calculating the infrared radiation brightness of the infrared attribute of the material on the scene surface according to the Planck formula.
10. The method of claim 9, wherein the step 31 comprises:
thickness of accumulated snow on dough units
Figure DEST_PATH_IMAGE041
If the unit is snow-free, the infrared attribute of the material on the surface of the scene is the original infrared attribute of the material of the unit;
thickness of accumulated snow on dough units
Figure 283519DEST_PATH_IMAGE042
And the infrared property of the material on the surface of the scene is the infrared property of the snow material.
CN202211030738.3A 2022-08-26 2022-08-26 Infrared scene numerical value calculation method under snow accumulation environment Active CN115098827B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211030738.3A CN115098827B (en) 2022-08-26 2022-08-26 Infrared scene numerical value calculation method under snow accumulation environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211030738.3A CN115098827B (en) 2022-08-26 2022-08-26 Infrared scene numerical value calculation method under snow accumulation environment

Publications (2)

Publication Number Publication Date
CN115098827A true CN115098827A (en) 2022-09-23
CN115098827B CN115098827B (en) 2022-11-08

Family

ID=83300049

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211030738.3A Active CN115098827B (en) 2022-08-26 2022-08-26 Infrared scene numerical value calculation method under snow accumulation environment

Country Status (1)

Country Link
CN (1) CN115098827B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117131312A (en) * 2023-10-20 2023-11-28 西安电子科技大学 Infrared scene numerical calculation method in rainy environment

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007232652A (en) * 2006-03-02 2007-09-13 Fujitsu Ltd Device and method for determining road surface condition
RU2460968C1 (en) * 2011-03-22 2012-09-10 Федеральное государственное бюджетное учреждение "Арктический и Антарктический научно-исследовательский институт" (ФГБУ "ААНИИ") Method of determining depth of snow cover on ice over water bodies
CN104143208A (en) * 2013-05-12 2014-11-12 哈尔滨点石仿真科技有限公司 Real-time rendering method for large-scale realistic snowscape
CN109145474A (en) * 2018-08-31 2019-01-04 北京晶品镜像科技有限公司 A kind of infrared imaging device recon effect analogy method suitable for emulation of fighting
US20210065565A1 (en) * 2019-08-27 2021-03-04 International Business Machines Corporation Remote surface condition assessment
CN113534297A (en) * 2021-07-16 2021-10-22 新疆林科院森林生态研究所 Device and method capable of accurately monitoring accumulated snow depth of forest in real time
CN114492006A (en) * 2022-01-17 2022-05-13 成都众享天地网络科技有限公司 Infrared radiation brightness calculation method for simulation

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007232652A (en) * 2006-03-02 2007-09-13 Fujitsu Ltd Device and method for determining road surface condition
RU2460968C1 (en) * 2011-03-22 2012-09-10 Федеральное государственное бюджетное учреждение "Арктический и Антарктический научно-исследовательский институт" (ФГБУ "ААНИИ") Method of determining depth of snow cover on ice over water bodies
CN104143208A (en) * 2013-05-12 2014-11-12 哈尔滨点石仿真科技有限公司 Real-time rendering method for large-scale realistic snowscape
CN109145474A (en) * 2018-08-31 2019-01-04 北京晶品镜像科技有限公司 A kind of infrared imaging device recon effect analogy method suitable for emulation of fighting
US20210065565A1 (en) * 2019-08-27 2021-03-04 International Business Machines Corporation Remote surface condition assessment
CN113534297A (en) * 2021-07-16 2021-10-22 新疆林科院森林生态研究所 Device and method capable of accurately monitoring accumulated snow depth of forest in real time
CN114492006A (en) * 2022-01-17 2022-05-13 成都众享天地网络科技有限公司 Infrared radiation brightness calculation method for simulation

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHU WANG: "Real-time Simulation of Infrared Scene", 《IEEE》 *
季新杰 等: "巡航导弹排气系统红外辐射数值模拟", 《激光与红外》 *
李长春 等: "基于FY3B-MWRI数据新疆区域积雪深度反演", 《遥感技术与应用》 *
郑坤鹏: "雪地表面的红外辐射计算及应用", 《科技信息》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117131312A (en) * 2023-10-20 2023-11-28 西安电子科技大学 Infrared scene numerical calculation method in rainy environment
CN117131312B (en) * 2023-10-20 2024-01-26 西安电子科技大学 Infrared scene numerical calculation method in rainy environment

Also Published As

Publication number Publication date
CN115098827B (en) 2022-11-08

Similar Documents

Publication Publication Date Title
CN106354152B (en) A kind of reentry trajectory optimum design method of pair of radial pattern no-fly zone
US9285460B2 (en) Method and system for estimating information related to a vehicle pitch and/or roll angle
CN104050716B (en) A kind of marine multiple target SAR image visual modeling method
CN111609760B (en) Intelligent sighting telescope shooting time determination method and system
CN115098827B (en) Infrared scene numerical value calculation method under snow accumulation environment
CN101114022A (en) Navigation multiple spectrum scanner geometric approximate correction method under non gesture information condition
CN109724586B (en) Spacecraft relative pose measurement method integrating depth map and point cloud
CN108107434B (en) Regional three-dimensional wind field picture splicing method based on double-Doppler radar inversion
CN103453806A (en) Projectile nutation angle extraction method based on radar Doppler data
CN103869708A (en) Complex infrared environment modeling method supporting multispectral engagement level simulation
CN116012377B (en) Unmanned aerial vehicle virtual observation image generation and positioning method based on satellite map
CN110488871A (en) A kind of unmanned plane path dynamic programming method based on two three-dimensional integratedization
CN113902663A (en) Air small target dynamic infrared simulation method and device capable of automatically adapting to weather
CN111462182B (en) Trajectory missile three-dimensional trajectory estimation method based on infrared early warning image
CN109211960A (en) A kind of dense city building scene thermal emission directionality strength calculation method
CN117131312B (en) Infrared scene numerical calculation method in rainy environment
CN113076634B (en) Multi-machine cooperative passive positioning method, device and system
CN103743488A (en) Infrared imaging simulation method for globe limb background characteristics of remote sensing satellite
CN110992399B (en) High-precision target atmosphere disturbance detection method
CN112329627A (en) High-altitude throwing object distinguishing method
CN116187217A (en) Unmanned aerial vehicle intelligent algorithm learning training platform construction method based on parallel simulation
Zhang et al. A stereo photogrammetry system for model position and attitude measurement in hypersonic wind tunnel testing
CN115808246A (en) Space normalization method for observation temperature of remote sensing thermal infrared sensor
CN109708668A (en) Line of sight measurement error range determining method and its device for video satellite
CN112927294B (en) Satellite orbit and attitude determination method based on single sensor

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