CN107247478A - A kind of quick calculation method of earth's surface solar radiation flux - Google Patents
A kind of quick calculation method of earth's surface solar radiation flux Download PDFInfo
- Publication number
- CN107247478A CN107247478A CN201710421754.8A CN201710421754A CN107247478A CN 107247478 A CN107247478 A CN 107247478A CN 201710421754 A CN201710421754 A CN 201710421754A CN 107247478 A CN107247478 A CN 107247478A
- Authority
- CN
- China
- Prior art keywords
- earth
- solar radiation
- radiation flux
- surface solar
- picture dot
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F1/00—Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
- G06F1/02—Digital function generators
- G06F1/03—Digital function generators working, at least partly, by table look-up
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/22—Indexing; Data structures therefor; Storage structures
- G06F16/2282—Tablespace storage structures; Management thereof
Abstract
This application discloses a kind of quick calculation method of earth's surface solar radiation flux, the problem of earth's surface solar radiation flux computational efficiency is low is solved.The quick calculation method of the earth's surface solar radiation flux, comprise the steps of, using atmospheric radiation transmission analog satellite data earth's surface solar radiation flux corresponding under different numerical value, earth's surface solar radiation flux data storehouse is set up, earth's surface solar radiation flux look-up table is obtained；The satellite data is downloaded, the satellite data includes solar zenith angle, cloud optical thickness, cloud Effective radius, aerosol optical depth, aerosol type, gray value；The earth's surface solar radiation flux is obtained according to the look-up table.The quick calculation method for the earth's surface solar radiation flux that the application is provided avoids the evaluation work of a large amount of repeated and redundants, substantially increases computational efficiency and computational accuracy.
Description
Technical field
The present invention relates to satellite remote sensing technology field, more particularly to a kind of computational methods of earth's surface solar radiation flux.
Background technology
At present, Surface radiation observation website is depended on solar radiation observation more, and ground observation station number is limited, and
The skewness in the whole world, it is difficult to obtain large-scale solar radiation distribution situation.
With the development of remote sensing technology, master is turned into as the data source for obtaining earth's surface solar radiation with Value of Remote Sensing Data
Stream.Moonscope is not by territorial restrictions, to estimate that the earth's surface solar radiation of arbitrary region provides new way.Wherein, it is static to defend
Star possesses very high temporal resolution, can monitor and obtain the satellite image of survey region in real time, and such as Chinese feature cloud system row are quiet
Only every 15 minutes of meteorological satellite, Japanese sunflower satellite (Himawari-8) can obtain the distant of a panel height spatial resolution in every 10 minutes
Feel image.
Satellite can obtain the situation of a wide range of air, current cloud detection and aerosol, Cloud microphysical parametric inversion technology into
Ripe, its atmospheric parameter product, such as cloud Effective radius, cloud optical thickness and aerosol optical depth have degree of precision,
Many satellites can provide comprehensive corresponding product data.Satellite atmosphere parameter product combination radiative transmission mode, can calculate arrival
The solar radiation flux of earth's surface.But radiative transmission mode needs designated centers wavelength when specifically calculating, and it is thick to input cloud optics
Need to repeat in numerous clouds and the aerosol Microphysicals such as degree, Effective radius and aerosol optical depth, calculating process
Calculate, computational efficiency is low, the earth's surface solar radiation consuming time for obtaining single picture dot is long.And weather satellite data amount is big, picture dot
It is many, only calculated with radiative transmission mode and be difficult to meet time precision requirement.
The scheme of some existing saving calculating times of tradition, is such as radiated earth's surface solar radiation by setting up atmospheric radiation top
Statistical model calculate earth's surface solar radiation, the defect of this method is its limited universality and ground observation data
It is limited.
The content of the invention
In view of this, to solve the problem of earth's surface solar radiation flux computational efficiency is low, the embodiment of the present application provides one
Plant the quick calculation method of earth's surface solar radiation flux.
The embodiment of the present application uses following technical proposals：
A kind of quick calculation method of earth's surface solar radiation flux, comprises the steps of, utilizes atmospheric radiation transmission
Analog satellite data earth's surface solar radiation flux corresponding under different numerical value, sets up earth's surface solar radiation flux data storehouse,
Obtain earth's surface solar radiation flux look-up table；The satellite data is downloaded, the satellite data includes solar zenith angle, cloud optics
Thickness, cloud Effective radius, aerosol optical depth, aerosol type, gray value；Obtain described according to the look-up table
Table solar radiation flux.
Further, the database includes cloud picture dot earth's surface solar radiation flux data storehouse and the clear sky picture dot earth's surface sun
Radiation flux data storehouse.
Further, the look-up table includes cloud picture dot earth's surface solar radiation flux look-up table and the clear sky picture dot earth's surface sun
Radiation flux look-up table.
Further, the quick calculation method of described earth's surface solar radiation flux, further comprising the steps of, judges picture dot
Type, the picture dot type includes clear sky picture dot, cloud picture dot and uncertain picture dot；Extract the institute of the clear sky picture dot and cloud picture dot
State satellite data.
Preferably for the gray value identical picture dot, by the result assignment searched first in the look-up table, no
Carry out repeating lookup again.
At least one above-mentioned technical scheme that the embodiment of the present application is used can reach following beneficial effect：
The application carries out the different observation geological information of simulation and atmospheric condition based on accurate atmospheric radiation transmission,
Each satellite geometry parameter and atmospheric condition and the solution relation of earth's surface solar radiation flux are set up using the method for setting up look-up table,
Set up the database of solar radiation flux, it is to avoid a large amount of calculating for repeating, expending time length, direct applied satellite product data are fast
Speed obtains the earth's surface solar radiation of view picture image.This method is applied to all satellites for providing input parameter, therefore, this mode
With wider range of universality, the evaluation work of a large amount of repeated and redundants is more avoided, computational efficiency is substantially increased with calculating essence
Degree.
The present invention can obtain the surface radiation flux pattern of view picture image, and input parameter is few, and computational efficiency is high, can be achieved
The acquisition earth's surface solar radiation of near real-time.
Brief description of the drawings
Accompanying drawing described herein is used for providing further understanding of the present application, constitutes the part of the application, this Shen
Schematic description and description please is used to explain the application, does not constitute the improper restriction to the application.
In the accompanying drawings：
Fig. 1 be under different solar zenith angles solar radiation flux with cloud optical thickness change curve；
Fig. 2 be under different cloud Effective radius solar radiation flux with cloud optical thickness change curve；
Fig. 3 is a kind of quick calculation method embodiment flow chart of earth's surface solar radiation flux；
Fig. 4 is the quick calculation method embodiment flow chart of another earth's surface solar radiation flux.
Embodiment
To make the purpose, technical scheme and advantage of the application clearer, below in conjunction with the application specific embodiment and
Technical scheme is clearly and completely described corresponding accompanying drawing.Obviously, described embodiment is only the application one
Section Example, rather than whole embodiments.Based on the embodiment in the application, those of ordinary skill in the art are not doing
Go out the every other embodiment obtained under the premise of creative work, belong to the scope of the application protection.
According to radiation transfer theory, scattering radiation and directly radiation can be divided into by reaching the solar radiation of earth's surface.Air spoke
It is to describe solar radiation from big pneumatic jack to penetrate mode, and earth's surface is reached by the scattering of atmospheric particles and cloud, reflection and absorption etc.
The model of process, conventional atmospheric radiation transmission has 6S, MODTRAN, LOWTRAN, RSTAR etc., and its input parameter is usual
Including atmospheric model, gas parameter setting, aerosol model, observation geological information and spectral information etc..Wherein, RSTAR exists
, it is necessary to input the satellite geometry angle of any pixel point during physical simulation, such as solar zenith angle, satellite zenith angle, the sun and defend
The wind speed of the parameters such as the difference at the azimuth of star, given calculative each band center wavelength and width, and every bit, will also
Select atmospheric model type, the secondary number of plies of air, the particle types and its volumetric mixture ratio, relative humidity for constituting aerosol model
Deng the relevant parameter of atmospheric portion.
Optical property parameter can be divided into three classes, and the first kind is Atmospheric models, including specify Atmospheric components composition, steam with
And ozone content etc.；Equations of The Second Kind is to specify aerosol optical attribute；3rd class is cloud optical properties.According to shortwave solar radiation
Feature, the rear two analogies first kind is more occupied an leading position, to reduce amount of calculation, and the first kind, which can be used, gives tacit consent to ginseng in atmospheric radiation pattern
Number.
Fig. 1 is that solar radiation flux is with the change curve of cloud optical thickness under different solar zenith angles, and Fig. 2 is that different clouds have
The change curve of solar radiation flux under particle radii with cloud optical thickness is imitated, wherein, Fig. 1 satellite zenith angle is equal to 140 °,
Relative bearing is equal to 90 °, and efficient radius of cloud particle is equal to 1 μm；Fig. 2 satellite zenith angle is equal to 140 °, relative bearing etc.
In 90 °, solar zenith angle is equal to 30 °.
The process of earth's surface solar radiation flux is being simulated with reference to the satellite data using the atmospheric radiation transmission
In, it is necessary to set observation geological information and atmospheric condition parameter information, the atmospheric condition parameter information include atmospheric gas into
Point, steam, aerosol, cloud and surface condition, the observation geological information includes solar zenith angle, view zenith angle and phase
Azimuthal.Wherein described solar zenith angle has considerable influence to the earth's surface solar radiation flux, especially described cloud optics
When thickness is less, solar radiation flux difference is larger；Cloud Effective radius has certain influence less when, works as value
When becoming big, influence weakens.
Fig. 3 is a kind of quick calculation method embodiment flow chart of earth's surface solar radiation flux, is comprised the steps of：
Step 31：Utilize atmospheric radiation transmission analog satellite data earth's surface sun spoke corresponding under different numerical value
Flux is penetrated, earth's surface solar radiation flux data storehouse is set up, earth's surface solar radiation flux look-up table is obtained.
Preferably, the atmospheric radiation transmission includes 6S, MODTRAN, LOWTRAN or RSTAR.
Further, the database includes cloud picture dot earth's surface solar radiation flux data storehouse and the clear sky picture dot earth's surface sun
Radiation flux data storehouse.
Further, the look-up table includes cloud picture dot earth's surface solar radiation flux look-up table and the clear sky picture dot earth's surface sun
Radiation flux look-up table.
Step 32：The satellite data is downloaded, it is effective that the satellite data includes solar zenith angle, cloud optical thickness, cloud
Particle radii, aerosol optical depth, aerosol type, gray value.
Further, the aerosol type includes city and rural two kinds of aerosol types.
Step 33：The earth's surface solar radiation flux is obtained according to the look-up table.
Preferably, the gray value identical picture dot, by the result assignment searched first in the look-up table, no longer enters
Row repeats to search.
Fig. 4 is the quick calculation method embodiment flow chart of another earth's surface solar radiation flux, is comprised the steps of：
Step 41：Utilize atmospheric radiation transmission analog satellite data earth's surface sun spoke corresponding under different numerical value
Flux is penetrated, cloud picture dot and clear sky picture dot earth's surface solar radiation flux data storehouse is set up, obtains cloud picture dot and clear sky picture dot earth's surface too
Positive radiation flux look-up table.
Preferably, the atmospheric radiation transmission includes 6S, MODTRAN, LOWTRAN or RSTAR.
Preferably, when setting up the clear sky picture dot earth's surface solar radiation flux data storehouse, the model of the solar zenith angle
Enclose for 0~180 °, set at intervals of 5 °, the scope of the aerosol optical depth is 0~3, is set at intervals of 0.05.
Preferably, when setting up the cloud picture dot earth's surface solar radiation flux data storehouse, the scope of the solar zenith angle
For 0~180 °, set at intervals of 5 °, the scope of the Effective radius of the cloud is 0~50 μm, at intervals of 2 μm, the cloud light
The scope for learning thickness is 0~80 μm, at intervals of 2 μm.
Preferably, when the cloud optical thickness is more than 80 μm, the cloud optical thickness change is to the earth's surface sun spoke
Penetrate variations of flux influence small, to avoid computing repeatedly, improve computational efficiency, when the cloud optical thickness is more than 80 μm, by cloud
Optical thickness is 80 μm of calculating.
Step 42：The satellite data is downloaded, it is effective that the satellite data includes solar zenith angle, cloud optical thickness, cloud
Particle radii, aerosol optical depth, aerosol type, gray value.
Further, the aerosol type includes city and rural two kinds of aerosol types.
Step 43：Judge picture dot type, the picture dot type includes clear sky picture dot, cloud picture dot and uncertain picture dot.
According to step 43 judged result, to be judged as clear sky picture dot and skies picture dot be considered as meet rapid calculating bar
Part picture dot, it is on the contrary then rapid calculating can not be carried out.
Step 44：Extract the satellite data of the clear sky picture dot and cloud picture dot.
Step 45：The earth's surface sun is obtained according to the cloud picture dot and clear sky picture dot earth's surface solar radiation flux look-up table
Radiation flux.
When step 43 judged result is clear sky picture dot, solar zenith angle, the aerosol light in the satellite data are extracted
Thickness and aerosol type are learned, according to the clear sky picture dot earth's surface solar radiation flux look-up table, the earth's surface sun of picture dot is obtained
Radiation flux.
When step 43 judged result is cloud picture dot, solar zenith angle, the cloud optical thickness in the satellite data are extracted
With cloud Effective radius, according to the cloud picture dot earth's surface solar radiation flux look-up table, the earth's surface solar radiation of picture dot is obtained
Flux.
Preferably, the gray value identical picture dot, by the result assignment searched first in the look-up table, no longer enters
Row repeats to search.
It should also be noted that, term " comprising ", "comprising" or its any other variant are intended to nonexcludability
Comprising, so that process, method or equipment including a series of key elements not only include those key elements, but also including not having
Other key elements being expressly recited, or also include for this process, method or the intrinsic key element of equipment.Not more
In the case of limitation, the key element limited by sentence "including a ...", it is not excluded that in the process including the key element, method
Or also there is other identical element in equipment
Embodiments herein is the foregoing is only, the application is not limited to.For those skilled in the art
For, the application can have various modifications and variations.It is all any modifications made within spirit herein and principle, equivalent
Replace, improve etc., it should be included within the scope of claims hereof.
Claims (10)
1. a kind of quick calculation method of earth's surface solar radiation flux, it is characterised in that comprise the steps of：
Using atmospheric radiation transmission analog satellite data earth's surface solar radiation flux corresponding under different numerical value, set up
Earth's surface solar radiation flux data storehouse, obtains earth's surface solar radiation flux look-up table；
The satellite data is downloaded, the satellite data includes solar zenith angle, cloud optical thickness, cloud Effective radius, gas
Colloidal sol optical thickness, aerosol type, gray value；
The earth's surface solar radiation flux is obtained according to the look-up table.
2. the quick calculation method of earth's surface solar radiation flux according to claim 1, it is characterised in that the database
Including cloud picture dot earth's surface solar radiation flux data storehouse and clear sky picture dot earth's surface solar radiation flux data storehouse.
3. the quick calculation method of earth's surface solar radiation flux according to claim 1, it is characterised in that the look-up table
Including cloud picture dot earth's surface solar radiation flux look-up table and clear sky picture dot earth's surface solar radiation flux look-up table.
4. the quick calculation method of earth's surface solar radiation flux according to claim 2, it is characterised in that also including following
Step：
Judge picture dot type, the picture dot type includes clear sky picture dot, cloud picture dot and uncertain picture dot；
Extract the satellite data of the clear sky picture dot and cloud picture dot.
5. the quick calculation method of earth's surface solar radiation flux according to claim 1, it is characterised in that the aerosol
Type includes city and rural two kinds of aerosol types.
6. the quick calculation method of earth's surface solar radiation flux according to claim 1, it is characterised in that for the ash
Angle value identical picture dot, by the result assignment searched in the look-up table first, no longer carries out repeating lookup.
7. the quick calculation method of the earth's surface solar radiation flux according to claim 1~6 any one, its feature exists
In the atmospheric radiation transmission is 6S, MODTRAN, LOWTRAN or RSTAR.
8. the quick calculation method of the earth's surface solar radiation flux according to claim 2~6 any one, its feature exists
In, when setting up the clear sky picture dot earth's surface solar radiation flux data storehouse, the scope of the solar zenith angle is 0~180 °,
Set at intervals of 5 °, the scope of the aerosol optical depth is 0~3, is set at intervals of 0.05.
9. the quick calculation method of earth's surface solar radiation flux according to claims 2 to 6, it is characterised in that setting up
During the cloud picture dot earth's surface solar radiation flux data storehouse, the scope of the solar zenith angle is 0~180 °, set at intervals of
5 °, the scope of the Effective radius of the cloud is 0~50 μm, and at intervals of 2 μm, the scope of the cloud optical thickness is 0~80 μ
M, at intervals of 2 μm.
10. the quick calculation method of earth's surface solar radiation flux according to claim 9, it is characterised in that when the cloud
It is 80 μm of calculating by optical thickness when optical thickness is more than 80 μm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710421754.8A CN107247478A (en) | 2017-06-07 | 2017-06-07 | A kind of quick calculation method of earth's surface solar radiation flux |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710421754.8A CN107247478A (en) | 2017-06-07 | 2017-06-07 | A kind of quick calculation method of earth's surface solar radiation flux |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107247478A true CN107247478A (en) | 2017-10-13 |
Family
ID=60018398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710421754.8A Pending CN107247478A (en) | 2017-06-07 | 2017-06-07 | A kind of quick calculation method of earth's surface solar radiation flux |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107247478A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110175375A (en) * | 2019-05-13 | 2019-08-27 | 中国科学院遥感与数字地球研究所 | A kind of earth's surface Calculation method for solar radiation based on deep learning |
| CN114003852A (en) * | 2022-01-04 | 2022-02-01 | 航天宏图信息技术股份有限公司 | Method and device for estimating down-ground solar radiation flux |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102338869A (en) * | 2011-06-20 | 2012-02-01 | 北京师范大学 | Inversion method and system of downlink shortwave radiation and photosynthetically active radiation data |
-
2017
- 2017-06-07 CN CN201710421754.8A patent/CN107247478A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102338869A (en) * | 2011-06-20 | 2012-02-01 | 北京师范大学 | Inversion method and system of downlink shortwave radiation and photosynthetically active radiation data |
Non-Patent Citations (3)
| Title |
|---|
| 孙洋；黄广辉；郝晓华: "结合极轨卫星ＭＯＤＩＳ和静止气象卫星", 《遥感技术与应用》 * |
| 张晓通: "全球陆表下行短波辐射和光合有效辐射反演算法研究", 《中国博士学位论文全文数据库基础科学辑》 * |
| 穆松宁: "《海-陆-气常规数据说明手册》", 30 September 2012, 北京理工大学出版社 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110175375A (en) * | 2019-05-13 | 2019-08-27 | 中国科学院遥感与数字地球研究所 | A kind of earth's surface Calculation method for solar radiation based on deep learning |
| CN114003852A (en) * | 2022-01-04 | 2022-02-01 | 航天宏图信息技术股份有限公司 | Method and device for estimating down-ground solar radiation flux |
| CN114003852B (en) * | 2022-01-04 | 2022-04-12 | 航天宏图信息技术股份有限公司 | Method and device for estimating down-ground solar radiation flux |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110111414B (en) | Orthographic image generation method based on three-dimensional laser point cloud | |
| Evan et al. | Derivation of an observation-based map of North African dust emission | |
| CN103942420B (en) | A kind of beam radia energy Method of fast estimating of building yardstick | |
| US11423610B2 (en) | Large-scale environment-modeling with geometric optimization | |
| US7500391B2 (en) | System and method for identifying the solar potential of rooftops | |
| CN105403664B (en) | A kind of large-scale point pollution source Atmospheric Environmental Impact Assessment method based on WRF CHEM | |
| Abreu et al. | Description of atmospheric conditions at the pierre auger observatory using the global data assimilation system (gdas) | |
| Martín et al. | Applying LIDAR datasets and GIS based model to evaluate solar potential over roofs: a review. | |
| Li et al. | The impact of observation nudging on simulated meteorology and ozone concentrations during DISCOVER-AQ 2013 Texas campaign | |
| CN106815850A (en) | The method that canopy density forest reserves very high is obtained based on laser radar technique | |
| CN102507586B (en) | Remote sensing monitoring method for carbon emission | |
| CN102854513B (en) | Cloud detection method of CCD (charge coupled device) data of environment first satellite HJ-1A/B | |
| CN103324916B (en) | Vehicle-mounted and aviation LiDAR data method for registering based on building profile | |
| An et al. | Three‐dimensional point cloud based sky view factor analysis in complex urban settings | |
| CN108874734B (en) | Global land rainfall inversion method | |
| CN109766824B (en) | Active and passive remote sensing data fusion classification method based on fuzzy evidence theory | |
| CN104318611A (en) | Three-dimensional tree re-establishing method based on laser scanning point cloud | |
| CN106780586B (en) | A kind of solar energy potential evaluation method based on ground laser point cloud | |
| CN105953921A (en) | Quick simulation method for earth observation radiation image on condition of different aerosol parameters | |
| CN114912370B (en) | Building photovoltaic potential analysis available area calculation method | |
| CN106680273A (en) | Retrieval method for high spatial resolution satellite surface reflectance | |
| CN116069882A (en) | Airspace grid diagram generating method | |
| CN107247478A (en) | A kind of quick calculation method of earth's surface solar radiation flux | |
| Zheng et al. | New insights into the correlation between lightning flash rate and size in thunderstorms | |
| Li et al. | Spatio-temporal changes of AOD in Xinjiang of China from 2000 to 2019: Which factor is more influential, natural factor or human factor? |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171013 |
|
| RJ01 | Rejection of invention patent application after publication |