CN107632327A - The determination method and apparatus of earth's surface solar shortwave radiation - Google Patents
The determination method and apparatus of earth's surface solar shortwave radiation Download PDFInfo
- Publication number
- CN107632327A CN107632327A CN201710865903.XA CN201710865903A CN107632327A CN 107632327 A CN107632327 A CN 107632327A CN 201710865903 A CN201710865903 A CN 201710865903A CN 107632327 A CN107632327 A CN 107632327A
- Authority
- CN
- China
- Prior art keywords
- cloud
- pixel
- ice
- earth
- atmospheric
- 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
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention discloses a kind of determination method and apparatus of earth's surface solar shortwave radiation.This method includes：Structure includes the information table of the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation；By Cloud product information, the atmospheric parameter of each pixel is obtained；And by searching the atmospheric parameter in described information table, determine the earth's surface solar shortwave radiation of each pixel.The present invention can realize following technique effect：First, earth's surface solar shortwave radiation can quickly, accurately be determined by Cloud product information and information table；Second, atmospheric conditions are divided into clear sky condition, water cloud condition and ice cloud condition, so as to the determination suitable for the earth's surface solar shortwave radiation under the conditions of these.
Description
Technical field
The present invention relates to meteorological field, more particularly to a kind of determination method and apparatus of earth's surface solar shortwave radiation.
Background technology
Earth's surface solar shortwave radiation is the important component of Radiation budget balance, and Land surface energy budget and ground
The important driving factors and land surface emissivity of atmosphere interaction and the important parameter of climate change research.The high-precision earth's surface sun is short
Wave radiation not only has realistic meaning, or human ecology ring to the forecast of the climate damage such as Global climate change and drought flood
The related decision-making in border provides help, has great importance for hazard forecasting, the estimation of solar energy etc..
Generally, earth's surface solar shortwave radiation is determined using radiation transfer equation method.This method computational accuracy is higher.So
And because its physical model is complicated, therefore it is computationally intensive, calculate time length, it is difficult to applied to extensive business.Also, the party
Method is applied to determine earth's surface solar shortwave radiation under the conditions of clear sky mostly, and can not be applied to determine earth's surface under the conditions of having cloud
Solar shortwave radiation.
The content of the invention
It is a primary object of the present invention to overcome using existing radiation transfer equation method to determine earth's surface sun shortwave spoke
The defects of calculating time is long when penetrating, and a kind of determination method and apparatus of new earth's surface solar shortwave radiation are provided.
A kind of according to an aspect of the invention, there is provided determination method of earth's surface solar shortwave radiation.This method includes：
Structure includes the information table of the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation；By Cloud product information, obtain
Take the atmospheric parameter of each pixel；And by searching the atmospheric parameter in described information table, determine each pixel
Earth's surface solar shortwave radiation.
According to an embodiment of the invention, structure includes the letter of the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation
Breath table includes：Atmospheric parameter corresponding to different atmospheric conditions is set respectively, wherein the atmospheric parameter bag corresponding to clear sky condition
Include at least one of：Solar zenith angle scope, aerosol optical depth, aerosol type；Corresponding to the air of water cloud condition
Parameter includes at least one of：Solar zenith angle scope, water efficient radius of cloud particle, water cloud optical thickness；Corresponding to ice cloud
The atmospheric parameter of condition includes at least one of：Solar zenith angle scope, ice efficient radius of cloud particle, ice cloud optical thickness；
And structure corresponds to the described information table of the different atmospheric conditions respectively.
According to an embodiment of the invention, respectively set corresponding to different atmospheric conditions atmospheric parameter include it is following at least it
One：Atmospheric parameter corresponding to clear sky condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1
Degree；Aerosol optical depth scope is 0 to 3, at intervals of 0.5, no unit；Aerosol type is arranged to urban type aerosol, township
Village's type aerosol and maritime aerosol；Atmospheric parameter corresponding to water cloud condition is set as follows：Solar zenith angle
Scope is 0 to 90 degree, at intervals of 1 degree；Water efficient radius of cloud particle is respectively 4,8,12,16,20 or 24, and unit is micron；Water
Cloud optical thickness scope is 0 to 150, at intervals of 1, no unit；The air corresponding to ice cloud condition is set to join as follows
Number：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree；Ice efficient radius of cloud particle be respectively 24,32,64,128,256 or
512, unit is micron；Ice cloud optical thickness scope is 0 to 150, at intervals of 1, no unit.
According to an embodiment of the invention, by Cloud product information, obtaining the atmospheric parameter of each pixel includes：Pass through
Cloud product information, it is determined that the atmospheric conditions of each pixel, wherein the atmospheric conditions include clear sky condition, water cloud condition and
It is at least one in ice cloud condition；And obtain the pixel corresponding atmospheric parameter under the atmospheric conditions.
According to an embodiment of the invention, obtain the pixel under the atmospheric conditions corresponding atmospheric parameter include it is following it
One：When it is clear sky that the Cloud product information, which indicates the pixel, solar zenith angle, the aerosol optical for obtaining the pixel are thick
Degree and aerosol type；When it is water cloud that the Cloud product information, which indicates the pixel, obtain the pixel solar zenith angle,
Water efficient radius of cloud particle and water cloud optical thickness；When it is ice cloud that the Cloud product information, which indicates the pixel, obtaining should
Solar zenith angle, ice efficient radius of cloud particle and the ice cloud optical thickness of pixel.
According to an embodiment of the invention, by searching the atmospheric parameter in described information table, each picture is determined
The earth's surface solar shortwave radiation of member includes one below：When the Cloud product information instruction pixel be clear sky when, by
Corresponding to the atmospheric parameter is searched in the clear sky information table of clear sky condition, the earth's surface sun shortwave spoke of each pixel is determined
Penetrate；When Cloud product information instruction pixel is water cloud, by being looked into the water cloud information table corresponding to water cloud condition
The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel；And when the Cloud product information refers to
When to show pixel be ice cloud, by searching the atmospheric parameter in the ice cloud information table corresponding to ice cloud condition, determine described every
The earth's surface solar shortwave radiation of individual pixel.
According to another aspect of the present invention, there is provided a kind of determining device of earth's surface solar shortwave radiation.The device includes：
Module is built, for building the information table for the corresponding relation for including atmospheric parameter and earth's surface solar shortwave radiation；Acquisition module,
For by Cloud product information, obtaining the atmospheric parameter of each pixel；And determining module, for by described information
The atmospheric parameter is searched in table, determines the earth's surface solar shortwave radiation of each pixel.
According to an embodiment of the invention, the structure module includes：Setting unit, correspond to difference greatly for setting respectively
The atmospheric parameter of gas bar part, wherein including at least one of corresponding to the atmospheric parameter of clear sky condition：Solar zenith angle scope,
Aerosol optical depth, aerosol type；Include at least one of corresponding to the atmospheric parameter of water cloud condition：Solar zenith angle
Scope, water efficient radius of cloud particle, water cloud optical thickness；Include at least one of corresponding to the atmospheric parameter of ice cloud condition：
Solar zenith angle scope, ice efficient radius of cloud particle, ice cloud optical thickness；And construction unit, correspond to for building respectively
The described information table of the different atmospheric conditions.
According to an embodiment of the invention, the setting unit is used at least one of：Set as follows corresponding
In the atmospheric parameter of clear sky condition：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree；Aerosol optical depth scope is 0
To 3, at intervals of 0.5, no unit；Aerosol type is arranged to urban type aerosol, rural type aerosol and maritime aerosol；
Atmospheric parameter corresponding to water cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree；Water
Efficient radius of cloud particle is respectively 4,8,12,16,20 or 24, and unit is micron；Water cloud optical thickness scope is 0 to 150, interval
For 1, no unit；Atmospheric parameter corresponding to ice cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree,
At intervals of 1 degree；Ice efficient radius of cloud particle is respectively 24,32,64,128,256 or 512, and unit is micron；Ice cloud optical thickness
Scope is 0 to 150, at intervals of 1, no unit.
According to an embodiment of the invention, the acquisition module includes：Determining unit, for by Cloud product information,
It is determined that the atmospheric conditions of each pixel, wherein the atmospheric conditions are included in clear sky condition, water cloud condition and ice cloud condition extremely
It is few one；And acquiring unit, for obtaining the pixel corresponding atmospheric parameter under the atmospheric conditions.
According to an embodiment of the invention, the acquiring unit is used for one below：When the Cloud product information indicates
When the pixel is clear sky, the solar zenith angle, aerosol optical depth and aerosol type of the pixel are obtained；When the Cloud
Product information indicates the pixel when being water cloud, obtains the solar zenith angle, water efficient radius of cloud particle and water cloud optics of the pixel
Thickness；When it is ice cloud that the Cloud product information, which indicates the pixel, the solar zenith angle of the pixel is obtained, ice cloud particle has
Imitate radius and ice cloud optical thickness.
According to an embodiment of the invention, the determining module is used for one below：When the Cloud product information indicates
When pixel is clear sky, by searching the atmospheric parameter in the clear sky information table corresponding to clear sky condition, determine described each
The earth's surface solar shortwave radiation of pixel；When Cloud product information instruction pixel is water cloud, by corresponding to water cloud
The atmospheric parameter is searched in the water cloud information table of condition, determines the earth's surface solar shortwave radiation of each pixel；And work as
When the Cloud product information instruction pixel is ice cloud, described in being searched in the ice cloud information table corresponding to ice cloud condition
Atmospheric parameter, determine the earth's surface solar shortwave radiation of each pixel.
According to above-described embodiment, there is provided a kind of determination method and apparatus of earth's surface solar shortwave radiation.First, this method
Earth's surface solar shortwave radiation can quickly, accurately be determined by Cloud product information and information table.Second, this method is by air
Condition divides into clear sky condition, water cloud condition and ice cloud condition, so as to suitable for the earth's surface sun shortwave under the conditions of these
The determination of radiation.
Brief description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, forms the part of the present invention, this hair
Bright schematic description and description is used to explain the present invention, does not form inappropriate limitation of the present invention.In the accompanying drawings：
Fig. 1 is the flow chart of the determination method of earth's surface solar shortwave radiation according to embodiments of the present invention；
Fig. 2 is the particular flow sheet of the determination method of earth's surface solar shortwave radiation according to embodiments of the present invention；
Fig. 3 is that existing radiation transfer equation method and the method progress earth's surface solar shortwave radiation of the present invention is respectively adopted
Determination consumed time contrast schematic diagram；And
Fig. 4 is the structured flowchart of the determining device of earth's surface solar shortwave radiation according to embodiments of the present invention.
Embodiment
Some illustrative embodiments shown in below with reference to the accompanying drawings come describe the present invention principle and spirit.Should
Understand, describe these embodiments just for the sake of better understood when those skilled in the art and then realize the present invention,
And the scope not limiting the invention in any way.
Fig. 1 is the flow chart of the determination method 100 of earth's surface solar shortwave radiation according to embodiments of the present invention.In method
After 100 start, in step S102, structure includes the information of the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation
Table.
The step can be realized in the following way.
First, the atmospheric parameter corresponding to different atmospheric conditions is set respectively.For example, join corresponding to the air of clear sky condition
Number can be set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree；Aerosol optical depth scope is
0 to 3, at intervals of 0.5, no unit；It is molten that aerosol type is arranged to urban type aerosol, rural type aerosol and maritime gas
Glue.For example, it can be set as follows corresponding to the atmospheric parameter of water cloud condition：Solar zenith angle scope is 0 to 90 degree,
At intervals of 1 degree；Water efficient radius of cloud particle is respectively 4,8,12,16,20 or 24, and unit is micron；Water cloud optical thickness scope
For 0 to 150, at intervals of 1, no unit.For example, it can be set as follows corresponding to the atmospheric parameter of ice cloud condition：Too
Positive zenith angular region is 0 to 90 degree, at intervals of 1 degree；Ice efficient radius of cloud particle is respectively 24,32,64,128,256 or 512,
Unit is micron；Ice cloud optical thickness scope is 0 to 150, at intervals of 1, no unit.
Second, pair for including atmospheric parameter and earth's surface solar shortwave radiation corresponding to different atmospheric conditions is built respectively
The information table that should be related to.Wherein, the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation can be by of the prior art
Various ways obtain to obtain, such as by atmospheric radiative transfer model RSTAR.
Specifically, atmospheric radiative transfer model RSTAR comes from Atmosphere and Ocean research institute of Tokyo Univ Japan Nakajima
Team is taught, it can simulate the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation under various atmospheric conditions.These are big
Gas bar part includes tropical atmosphere condition, middle latitude summer atmospheric conditions, middle latitude winter atmospheric conditions, high latitude summer big gas bar
Part, high latitude winter atmospheric conditions, United States standard atmosphere condition.The height that it is simulated can be from sea level 0km to 120km at most
It is divided into 50 layers.It additionally provides 11 kinds of aerosol types：Water, ice crystal, dust, coal smoke, volcanic ash, 75% sulfuric acid, rural type
Aerosol, urban type aerosol, sea spray, convection current stratotype aerosol, yellow sand.It is on the air ginseng under various atmospheric conditions
Number sets more, output result to enrich, computational accuracy is high.
Next, this method 100 proceeds to step S104, by Cloud product information, the air of each pixel is obtained
Parameter.
The step can be realized in the following way.First, by Cloud product information, it is determined that each pixel is big
Gas bar part, the wherein atmospheric conditions include at least one in clear sky condition, water cloud condition and ice cloud condition.Second, obtaining should
Pixel corresponding atmospheric parameter under the atmospheric conditions.For example, when it is clear sky that the Cloud product information, which indicates the pixel, obtain
Take the solar zenith angle, aerosol optical depth and aerosol type of the pixel；When the Cloud product information indicates the pixel
For water cloud when, obtain the solar zenith angle, water efficient radius of cloud particle and water cloud optical thickness of the pixel；When the Cloud product
Information indicates the pixel when being ice cloud, obtains the solar zenith angle, ice efficient radius of cloud particle and ice cloud optical thickness of the pixel.
Then, this method 100 proceeds to step S106, by searching the atmospheric parameter in the information table, determines that this is every
The earth's surface solar shortwave radiation of individual pixel.Specifically, when Cloud product information instruction pixel is clear sky, by right
The atmospheric parameter should be searched in the clear sky information table of clear sky condition, determine the earth's surface solar shortwave radiation of each pixel；When
When Cloud product information instruction pixel is water cloud, by searching the air in the water cloud information table corresponding to water cloud condition
Parameter, determine the earth's surface solar shortwave radiation of each pixel；And when Cloud product information instruction pixel is ice cloud,
By searching the atmospheric parameter in the ice cloud information table corresponding to ice cloud condition, the earth's surface sun shortwave of each pixel is determined
Radiation.
Fig. 2 is the particular flow sheet of the determination method of earth's surface solar shortwave radiation according to embodiments of the present invention.Such as Fig. 2 institutes
Show, this method includes：Step S1, based on atmospheric radiative transfer model RSTAR, structure calculates the information of earth's surface solar shortwave radiation
Table；And step S2, earth's surface solar shortwave radiation is determined based on the information table built in Cloud product information and step S1.
It is described separately below.
Step S1 can include one of following three kinds of situations.
Step S1.1 builds clear sky information table.
The construction method of clear sky information table is：Using atmospheric radiative transfer model RSTAR, under the conditions of simulating clear sky, earth's surface is too
Positive shortwave radiation.Arrange parameter is as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree, i.e., [0,1,2,3 ... 88,89,
90]；Aerosol optical depth scope is 0 to 3, at intervals of 0.5, i.e. [0.0,0.5,1.0,1.5,2.0,2.5,3.0], no list
Position；Aerosol type is arranged to urban type aerosol, rural type aerosol and maritime aerosol.Therefore constructed clear sky letter
The combination kind number of different solar zenith angles, aerosol optical depth and aerosol type is in breath table：91 × 7 × 3=1911.
Step S1.2 builds water cloud information table.
The construction method of water cloud information table is：Using atmospheric radiative transfer model RSTAR, under the conditions of simulating water cloud, earth's surface is too
Positive shortwave radiation.Arrange parameter is as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree, i.e., [0,1,2,3 ... 88,89,
90]；Water efficient radius of cloud particle is respectively 4,8,12,16,20 or 24, and unit is micron；Water cloud optical thickness scope be 0 to
150, at intervals of 1, i.e. [0,1,2,3 ... 148,149,150], no unit.Therefore different sun zenith in the information table formed
The combination kind number at angle, water efficient radius of cloud particle and water cloud optical thickness is：91 × 6 × 151=82446.
Step S1.3 builds ice cloud information table
The construction method of ice cloud information table is：Using atmospheric radiative transfer model RSTAR, under the conditions of simulating ice cloud, earth's surface is too
Positive shortwave radiation.Arrange parameter is as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree, i.e., [0,1,2,3 ... 88,89,
90]；Ice efficient radius of cloud particle is respectively 24,32,64,128,256 or 512, and unit is micron；Ice cloud optical thickness scope is
0 to 150, at intervals of 1, i.e. [0,1,2,3 ... 148,149,150], no unit.Therefore the different sun in the information table formed
The combination kind number of zenith angle, ice efficient radius of cloud particle and ice cloud optical thickness is：91 × 6 × 151=82446.
It should be noted that above-mentioned parameter set-up mode is as just an example in the embodiment of the present invention, not structure
Any restriction of the present invention in pairs.That is, setting thought of the those skilled in the art in examples detailed above, can also be selected
Other number ranges and precision are selected, these selections should all include protection scope of the present invention.
Step S2 can include successively two sub-steps：Step S2.1 and step S2.2.
Step S2.1 utilizes programming language IDL or Python, reads Cloud product information.Wherein the Cloud product is believed
Breath can include：Japanese sunflower -8, the satellite of sunflower -9, Terra/MODIS, Aqua/MODIS, NPP/VIIRS satellite in the U.S. with
And domestic wind and cloud series of satellites.
When the Cloud product information instruction pixel of reading is clear sky, is then handled according to clear sky, that is, read the picture
The aerosol product data of member, and then read the solar zenith angle, aerosol type and aerosol optical depth of the pixel.
When the Cloud product information instruction pixel of reading is cloud, then continue to read the cloud phase mark of the pixel,
When cloud phase is labeled as water cloud, then continue to read the solar zenith angle of the pixel, water efficient radius of cloud particle and water cloud optics
Thickness；When the cloud phase of reading is labeled as ice cloud, then continue to read the solar zenith angle of the pixel, ice efficient radius of cloud particle
With ice cloud optical thickness.
Step S2.2 is under programming language IDL or Python, by searching atmospheric parameter in information table, it is determined that each picture
The earth's surface solar shortwave radiation of member.
For clear sky condition, atmospheric parameter (i.e. solar zenith angle, aerosol optical depth, the gas of each pixel are being obtained
Colloidal sol type) after, the atmospheric parameter or the earth's surface sun shortwave closest to the atmospheric parameter are searched in clear sky information table
Radiation, the result as the earth's surface solar shortwave radiation under the conditions of clear sky.
For there is cloud condition, it is divided into the processing of the following two kinds method：
The first is water cloud condition, is obtaining atmospheric parameter (the i.e. effective particle of solar zenith angle, water cloud half of each pixel
Footpath, water cloud optical thickness) after, searched in water cloud information table the atmospheric parameter or closest to the atmospheric parameter earth's surface too
Positive shortwave radiation, the result as the earth's surface solar shortwave radiation under the conditions of water cloud；
Second is ice cloud condition, is obtaining atmospheric parameter (i.e. solar zenith angle, the effective particle of ice cloud half of each pixel
Footpath, ice cloud optical thickness) after, searched in ice cloud information table the atmospheric parameter or closest to the atmospheric parameter earth's surface too
Positive shortwave radiation, the result as the earth's surface solar shortwave radiation under the conditions of ice cloud.
Fig. 3 is that existing radiation transfer equation method and the method progress earth's surface solar shortwave radiation of the present invention is respectively adopted
Determination consumed time contrast schematic diagram.As shown in figure 3, using existing radiation transfer equation method (conventionally calculation side
Method), it is assumed that each pixel was fixed time as 3 seconds really, and pixel number is 2401 × 2401=5764801, then the time needed is:
5764801 × 3=17294403 seconds；Using the method for the present invention, the cloud of sunflower -8 product and aerosol data only need 55 seconds left sides
The right side, the inverse time (about 1 × 10 of average each pixel-5Second /).As can be seen here, compared to existing radiation transfer equation
Method, can significantly it be saved the consumed time using the method for the present invention.
Fig. 4 is the structured flowchart of the determining device 400 of earth's surface solar shortwave radiation according to embodiments of the present invention.Such as Fig. 4
Shown, the device 400 includes：Module 402 is built, the correspondence of atmospheric parameter and earth's surface solar shortwave radiation is included for building
The information table of relation；Acquisition module 404, for by Cloud product information, obtaining the atmospheric parameter of each pixel；And really
Cover half block 406, for by searching the atmospheric parameter in the information table, determining the earth's surface sun shortwave spoke of each pixel
Penetrate.
According to an embodiment of the invention, the structure module 402 includes：Setting unit, correspond to difference for setting respectively
The atmospheric parameter of atmospheric conditions, wherein including at least one of corresponding to the atmospheric parameter of clear sky condition：Solar zenith angle model
Enclose, aerosol optical depth, aerosol type；Include at least one of corresponding to the atmospheric parameter of water cloud condition：Sun day
Drift angle scope, water efficient radius of cloud particle, water cloud optical thickness；Corresponding to ice cloud condition atmospheric parameter include it is following at least it
One：Solar zenith angle scope, ice efficient radius of cloud particle, ice cloud optical thickness；And construction unit, for building respectively correspondingly
In the information table of the different atmospheric conditions.
According to an embodiment of the invention, the setting unit is used at least one of：Set and correspond to as follows
The atmospheric parameter of clear sky condition：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree；Aerosol optical depth scope be 0 to
3, at intervals of 0.5, no unit；Aerosol type is arranged to urban type aerosol, rural type aerosol and maritime aerosol；Press
According to following manner, the atmospheric parameter corresponding to water cloud condition is set：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 degree；Water cloud
Particle effective radius is respectively 4,8,12,16,20 or 24, and unit is micron；Water cloud optical thickness scope is 0 to 150, at intervals of
1, no unit；Atmospheric parameter corresponding to ice cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree,
It is divided into 1 degree；Ice efficient radius of cloud particle is respectively 24,32,64,128,256 or 512, and unit is micron；Ice cloud optical thickness model
Enclose for 0 to 150, at intervals of 1, no unit.
According to an embodiment of the invention, the acquisition module 404 includes：Determining unit, for by Cloud product information,
It is determined that the atmospheric conditions of each pixel, the wherein atmospheric conditions are included in clear sky condition, water cloud condition and ice cloud condition at least
One；And acquiring unit, for obtaining the pixel corresponding atmospheric parameter under the atmospheric conditions.
According to an embodiment of the invention, the acquiring unit is used for one below：When the Cloud product information indicates the picture
When member is clear sky, the solar zenith angle, aerosol optical depth and aerosol type of the pixel are obtained；When the Cloud product is believed
When the breath instruction pixel is water cloud, the solar zenith angle, water efficient radius of cloud particle and water cloud optical thickness of the pixel are obtained；When
The Cloud product information indicates the pixel when being ice cloud, obtain the solar zenith angle of the pixel, ice efficient radius of cloud particle and
Ice cloud optical thickness.
According to an embodiment of the invention, the determining module 406 is used for one below：When the Cloud product information indicates picture
When member is clear sky, by searching the atmospheric parameter in the clear sky information table corresponding to clear sky condition, each pixel is determined
Earth's surface solar shortwave radiation；When Cloud product information instruction pixel is water cloud, by the water corresponding to water cloud condition
The atmospheric parameter is searched in cloud information table, determines the earth's surface solar shortwave radiation of each pixel；And when the Cloud product
When information instruction pixel is ice cloud, by searching the atmospheric parameter in the ice cloud information table corresponding to ice cloud condition, it is determined that should
The earth's surface solar shortwave radiation of each pixel.
In summary, according to the abovementioned embodiments of the present invention, there is provided a kind of determination method of earth's surface solar shortwave radiation
And device.This method includes：Structure includes the information table of the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation；Pass through
Cloud product information, obtain the atmospheric parameter of each pixel；And by searching the atmospheric parameter in described information table,
Determine the earth's surface solar shortwave radiation of each pixel.The present invention can realize following technique effect：First, pass through Cloud
Product information and information table can quickly, accurately determine earth's surface solar shortwave radiation；Second, atmospheric conditions are divided into clear sky bar
Part, water cloud condition and ice cloud condition, so as to the determination suitable for the earth's surface solar shortwave radiation under the conditions of these.
The realization and all feature operations provided herein of the present invention with Fundamental Digital Circuit or can use computer
Software, firmware or hardware, including this specification and its structure disclosed in structural equivalents or one of those or it is more
Individual combination is realized.The realization of the present invention can be implemented as one or more computer program products, i.e., computer-readable
One or more modules of the computer program instructions encoded on medium, these instructions are performed or used by data processing equipment
With the operation of control data processing unit.The computer-readable medium can be machine readable storage device, machine readable storage
Substrate, memory devices, the composition for influenceing machine readable transmitting signal or one or more combination.Term
" data processing equipment " covers all devices, equipment and the machine for processing data, including such as programmable processor, calculating
Machine or multiple processors or computer.In addition to hardware, the device can include creating for described computer program
The code of performing environment, such as form processor firmware, protocol stack, data base management system, operating system or therein one
The code of individual or multiple combination.
Computer program (also referred to as program, software, software application, script or code) can use any type of programming language
(including compiler language or interpretative code) is sayed to write, and computer program can be disposed with any form, including conduct
Stand-alone program either as module, part, subroutine or is adapted to other units used in a computing environment.Computer program
The file not necessarily corresponded in file system.Program, which can be stored in, keeps other programs or data (such as markup language
The one or more scripts stored in document) file part in, be stored in the single text for being exclusively used in described program
In part, either it is stored in multiple coordinated files and (such as stores the text of the part of one or more modules, subprogram or code
Part) in.Computer program can be deployed on a computer to perform, or at a website or is being distributed in
At multiple websites and pass through interconnection of telecommunication network multiple computers on perform.
Process described in the present invention and logic flow can be by one or more of the one or more computer programs of execution
Individual programmable processor performs with by operating input data and generating output come perform function.The process and logic flow
It can be performed by dedicated logic circuit, and device can also be embodied as the dedicated logic circuit, the dedicated logic circuit example
Such as it is FPGA (field programmable gate array) or ASIC (application specific integrated circuit).
Being suitably executed the processor of computer program includes for example general and both special microprocessors and any species
Digital computer any one or more processors.Generally, processor is from read-only storage or random access memory
Or the two receives instruction and data.The element of computer can include for the processor of execute instruction and refer to for storing
Order and one or more memory devices of data.Generally, computer will also include one or more mass memory units so as to
Data storage, or the computer are operationally coupled to be received from mass memory unit or transmit data to mass memory unit
Or the two, the mass memory unit is, for example, disk, magneto-optic disk or CD.However, computer need not be with such
Equipment.In addition, computer can be embedded in another equipment, another equipment is, for example, mobile phone, personal digital assistant
(PDA), Mobile audio player, global positioning system (GPS) receiver etc..It is adapted to storage computer program instructions and data
Computer-readable medium includes nonvolatile memory, medium and the memory devices of form of ownership, including for example：Semiconductor is deposited
Storage device, such as EPROM, EEPROM and flash memory device；Disk, such as built-in hard disk or removable disk；Magneto-optic disk；And CD ROM
With DVD-ROM disks.The processor and memory can be supplemented or are incorporated in the dedicated logic circuit with dedicated logic circuit.
Interacted to provide with user, of the invention realize can be with setting for the display to user's display information
Standby (such as CRT (cathode-ray tube) or LCD (liquid crystal display) monitor) and keyboard and pointing device (such as mouse or with
Track ball, input can be provided to computer by its user) computer on realize.Other kinds of equipment can also be used
Interacted to provide with user；For example, the feedback provided a user can be any type of sense feedback, such as vision is anti-
Feedback, audio feedback or touch feedback；And the input from user can receive in any form, including the sense of hearing, voice
Or sense of touch.
Although the present invention includes some details, but these details should not be interpreted as to the disclosure or claimed
Content scope limitation, but be appreciated that to the present invention example implementation feature description.In the present invention
Some features described in the situation being implemented separately can also be combined to provide with single realization.On the contrary, in single realization
Each feature described in situation can also be provided in multiple realizations or carried in any suitable sub-portfolio respectively
For.Although in addition, it can describe feature as so protecting to perform and even initially require in some combination above, so
And one or more features from claimed combination can be removed from combination in some cases, and be claimed
Combination can be related to the change of sub-portfolio or sub-portfolio.
Similarly, although describing operation according to particular order in the accompanying drawings, but this is understood not to require this
The operation of sample either performs or required that all illustrated operations are all held according to shown particular order according to sequential order
OK, to realize desired result.In some circumstances, multitask and parallel processing are probably favourable.It is in addition, described above
The separation of various system units in realization is understood not to require such separation in all realizations, and should
Understanding, described program element and system generally can integrate or be packaged into single software product multiple
Software product.
Therefore, specific implementation of the invention has been described, and other are realized in the range of following claims.Example
Such as, the action described in claim can perform in a different order, and these actions still can realize expectation
Result.A large amount of realize has been described.It will be appreciated, however, that can be in situation without departing from the spirit and scope of the present invention
Under make various modifications.It is, for example, possible to use each form of flow illustrated above, wherein step can be reordered,
Addition removes.Therefore, other are realized in the range of following claims.
Claims (10)
- A kind of 1. determination method of earth's surface solar shortwave radiation, it is characterised in that including：Structure includes the information table of the corresponding relation of atmospheric parameter and earth's surface solar shortwave radiation；By Cloud product information, the atmospheric parameter of each pixel is obtained；AndBy searching the atmospheric parameter in described information table, the earth's surface solar shortwave radiation of each pixel is determined.
- 2. the method as described in claim 1, it is characterised in that by Cloud product information, obtain the air of each pixel Parameter includes：By Cloud product information, it is determined that the atmospheric conditions of each pixel, wherein the atmospheric conditions include clear sky condition, water It is at least one in cloud condition and ice cloud condition；AndObtain the pixel corresponding atmospheric parameter under the atmospheric conditions.
- 3. method as claimed in claim 2, it is characterised in that structure includes atmospheric parameter and earth's surface solar shortwave radiation The information table of corresponding relation includes：Respectively set corresponding to different atmospheric conditions atmospheric parameter, wherein corresponding to clear sky condition atmospheric parameter include it is following At least one：Solar zenith angle scope, aerosol optical depth, aerosol type；Corresponding to the atmospheric parameter bag of water cloud condition Include at least one of：Solar zenith angle scope, water efficient radius of cloud particle, water cloud optical thickness；Corresponding to ice cloud condition Atmospheric parameter includes at least one of：Solar zenith angle scope, ice efficient radius of cloud particle, ice cloud optical thickness；AndDescribed information table of the structure corresponding to the different atmospheric conditions respectively.
- 4. method as claimed in claim 3, it is characterised in that the atmospheric parameter bag corresponding to different atmospheric conditions is set respectively Include at least one of：Atmospheric parameter corresponding to clear sky condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 Degree；Aerosol optical depth scope is 0 to 3, at intervals of 0.5, no unit；Aerosol type is arranged to urban type aerosol, township Village's type aerosol and maritime aerosol；Atmospheric parameter corresponding to water cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 Degree；Water efficient radius of cloud particle is respectively 4,8,12,16,20 or 24, and unit is micron；Water cloud optical thickness scope be 0 to 150, at intervals of 1, no unit；Atmospheric parameter corresponding to ice cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 Degree；Ice efficient radius of cloud particle is respectively 24,32,64,128,256 or 512, and unit is micron；Ice cloud optical thickness scope is 0 To 150, at intervals of 1, no unit.
- 5. method as claimed in claim 3, it is characterised in that by searching the atmospheric parameter in described information table, really The earth's surface solar shortwave radiation of fixed each pixel includes one below：When Cloud product information instruction pixel is clear sky, by being looked into the clear sky information table corresponding to clear sky condition The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel；When Cloud product information instruction pixel is water cloud, by being looked into the water cloud information table corresponding to water cloud condition The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel；AndWhen Cloud product information instruction pixel is ice cloud, by being looked into the ice cloud information table corresponding to ice cloud condition The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel.
- A kind of 6. determining device of earth's surface solar shortwave radiation, it is characterised in that including：Module is built, for building the information table for the corresponding relation for including atmospheric parameter and earth's surface solar shortwave radiation；Acquisition module, for by Cloud product information, obtaining the atmospheric parameter of each pixel；AndDetermining module, for by searching the atmospheric parameter in described information table, determining the earth's surface of each pixel too Positive shortwave radiation；The acquisition module includes：Determining unit, for by Cloud product information, it is determined that the atmospheric conditions of each pixel, wherein the atmospheric conditions bag Include at least one in clear sky condition, water cloud condition and ice cloud condition；AndAcquiring unit, for obtaining the pixel corresponding atmospheric parameter under the atmospheric conditions.
- 7. device as claimed in claim 6, it is characterised in that the structure module includes：Setting unit, for setting the atmospheric parameter corresponding to different atmospheric conditions respectively, wherein corresponding to the big of clear sky condition Gas parameter includes at least one of：Solar zenith angle scope, aerosol optical depth, aerosol type；Corresponding to water cloud bar The atmospheric parameter of part includes at least one of：Solar zenith angle scope, water efficient radius of cloud particle, water cloud optical thickness；It is right The atmospheric parameter of Ying Yubing cloud conditions includes at least one of：Solar zenith angle scope, ice efficient radius of cloud particle, ice cloud light Learn thickness；AndConstruction unit, for building the described information table corresponding to the different atmospheric conditions respectively.
- 8. device as claimed in claim 7, it is characterised in that the setting unit is used at least one of：Atmospheric parameter corresponding to clear sky condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 Degree；Aerosol optical depth scope is 0 to 3, at intervals of 0.5, no unit；Aerosol type is arranged to urban type aerosol, township Village's type aerosol and maritime aerosol；Atmospheric parameter corresponding to water cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 Degree；Water efficient radius of cloud particle is respectively 4,8,12,16,20 or 24, and unit is micron；Water cloud optical thickness scope be 0 to 150, at intervals of 1, no unit；Atmospheric parameter corresponding to ice cloud condition is set as follows：Solar zenith angle scope is 0 to 90 degree, at intervals of 1 Degree；Ice efficient radius of cloud particle is respectively 24,32,64,128,256 or 512, and unit is micron；Ice cloud optical thickness scope is 0 To 150, at intervals of 1, no unit.
- 9. the device as described in claim 6~8 any one, it is characterised in that the acquiring unit is used for one below：When it is clear sky that the Cloud product information, which indicates the pixel, solar zenith angle, the aerosol optical of the pixel are obtained Thickness and aerosol type；When it is water cloud that the Cloud product information, which indicates the pixel, the solar zenith angle of the pixel is obtained, water cloud particle has Imitate radius and water cloud optical thickness；When it is ice cloud that the Cloud product information, which indicates the pixel, the solar zenith angle of the pixel is obtained, ice cloud particle has Imitate radius and ice cloud optical thickness.
- 10. the device as described in claim 6~8 any one, it is characterised in that the determining module is used for one below：When Cloud product information instruction pixel is clear sky, by being looked into the clear sky information table corresponding to clear sky condition The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel；When Cloud product information instruction pixel is water cloud, by being looked into the water cloud information table corresponding to water cloud condition The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel；AndWhen Cloud product information instruction pixel is ice cloud, by being looked into the ice cloud information table corresponding to ice cloud condition The atmospheric parameter is looked for, determines the earth's surface solar shortwave radiation of each pixel.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710865903.XA CN107632327A (en) | 2017-09-22 | 2017-09-22 | The determination method and apparatus of earth's surface solar shortwave radiation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710865903.XA CN107632327A (en) | 2017-09-22 | 2017-09-22 | The determination method and apparatus of earth's surface solar shortwave radiation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107632327A true CN107632327A (en) | 2018-01-26 |
Family
ID=61103490
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710865903.XA Pending CN107632327A (en) | 2017-09-22 | 2017-09-22 | The determination method and apparatus of earth's surface solar shortwave radiation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107632327A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109945969A (en) * | 2019-03-04 | 2019-06-28 | 张艳 | The method and device for determining earth radiation budget are observed based on meteorological satellite |
| CN110175375A (en) * | 2019-05-13 | 2019-08-27 | 中国科学院遥感与数字地球研究所 | A kind of earth's surface Calculation method for solar radiation based on deep learning |
| CN110222301A (en) * | 2019-06-03 | 2019-09-10 | 中国科学院遥感与数字地球研究所 | Earth's surface solar shortwave radiation calculation method under the conditions of a kind of haze |
Citations (2)
| 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 |
| CN103955749A (en) * | 2014-03-19 | 2014-07-30 | 内蒙古电力勘测设计院 | Sun direct radiation value prediction method and system |
-
2017
- 2017-09-22 CN CN201710865903.XA patent/CN107632327A/en active Pending
Patent Citations (2)
| 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 |
| CN103955749A (en) * | 2014-03-19 | 2014-07-30 | 内蒙古电力勘测设计院 | Sun direct radiation value prediction method and system |
Non-Patent Citations (1)
| Title |
|---|
| 赵静等: "基于SBDART模型的改进短波辐射传输模型", 《HTTP://WWW.DOC88.COM/P-7089633658212.HTML》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109945969A (en) * | 2019-03-04 | 2019-06-28 | 张艳 | The method and device for determining earth radiation budget are observed based on meteorological satellite |
| CN110175375A (en) * | 2019-05-13 | 2019-08-27 | 中国科学院遥感与数字地球研究所 | A kind of earth's surface Calculation method for solar radiation based on deep learning |
| CN110222301A (en) * | 2019-06-03 | 2019-09-10 | 中国科学院遥感与数字地球研究所 | Earth's surface solar shortwave radiation calculation method under the conditions of a kind of haze |
| CN110222301B (en) * | 2019-06-03 | 2020-09-04 | 中国科学院遥感与数字地球研究所 | Surface solar short wave radiation calculation method under haze condition |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wedi et al. | A baseline for global weather and climate simulations at 1 km resolution | |
| Themens et al. | The empirical Canadian high arctic ionospheric model (E‐CHAIM): NmF2 and hmF2 | |
| Woldemeskel et al. | Quantification of precipitation and temperature uncertainties simulated by CMIP3 and CMIP5 models | |
| Nielsen et al. | Increase in Arctic coastal erosion and its sensitivity to warming in the twenty-first century | |
| Wang et al. | Changes in global ocean wave heights as projected using multimodel CMIP5 simulations | |
| Lee et al. | Assimilation of FORMOSAT‐3/COSMIC electron density profiles into a coupled thermosphere/ionosphere model using ensemble Kalman filtering | |
| Church et al. | Evaluating the ability of process based models to project sea-level change | |
| Chau et al. | Quiet time ionospheric variability over Arecibo during sudden stratospheric warming events | |
| Zou et al. | Dynamical downscaling of historical climate over CORDEX East Asia domain: A comparison of regional ocean‐atmosphere coupled model to stand‐alone RCM simulations | |
| Hsin et al. | Seasonal fluctuations of the surface North Equatorial Countercurrent (NECC) across the Pacific basin | |
| Pedatella et al. | Multimodel comparison of the ionosphere variability during the 2009 sudden stratosphere warming | |
| Ajith et al. | On the fresh development of equatorial plasma bubbles around the midnight hours of June solstice | |
| Wang et al. | Physical processes driving the response of the F2 region ionosphere to the 21 August 2017 solar eclipse at Millstone Hill | |
| Ren et al. | Equinoctial asymmetry of ionospheric vertical plasma drifts and its effect on F‐region plasma density | |
| Lee et al. | Aura Microwave Limb Sounder observations of the polar middle atmosphere: Dynamics and transport of CO and H2O | |
| Neely III et al. | Biases in southern hemisphere climate trends induced by coarsely specifying the temporal resolution of stratospheric ozone | |
| Wang et al. | Model performance in spatiotemporal patterns of precipitation: New methods for identifying value added by a regional climate model | |
| Holmes et al. | Compensating biases and a noteworthy success in the CMIP5 representation of Antarctic sea ice processes | |
| Ciesielski et al. | Relationships between radiation, clouds, and convection during DYNAMO | |
| McIlhattan et al. | Arctic clouds and precipitation in the Community Earth System Model version 2 | |
| Balaguru et al. | Increased US coastal hurricane risk under climate change | |
| Hurwitz et al. | Assessment of the breakup of the Antarctic polar vortex in two new chemistry‐climate models | |
| Hormann et al. | Interannual North Equatorial Countercurrent variability and its relation to tropical Atlantic climate modes | |
| CN107632327A (en) | The determination method and apparatus of earth's surface solar shortwave radiation | |
| Lappen et al. | The role of tilted heating in the evolution of the MJO |
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: 20180126 |
|
| RJ01 | Rejection of invention patent application after publication |