CN107767448A - Consider the earth's surface solar radiation remote sensing computational methods and system of cloud orographic effect - Google Patents
Consider the earth's surface solar radiation remote sensing computational methods and system of cloud orographic effect Download PDFInfo
- Publication number
- CN107767448A CN107767448A CN201710928918.6A CN201710928918A CN107767448A CN 107767448 A CN107767448 A CN 107767448A CN 201710928918 A CN201710928918 A CN 201710928918A CN 107767448 A CN107767448 A CN 107767448A
- Authority
- CN
- China
- Prior art keywords
- radiation
- earth
- cloud
- solar radiation
- shadow region
- 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
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The present invention provides the earth's surface solar radiation remote sensing computational methods and system for considering cloud orographic effect, wherein, methods described includes：S1, using the recognizer of rolling topography cloud shadow region, cloud shade of the goal pels under MODEL OVER COMPLEX TOPOGRAPHY is obtained, distinguish earth's surface solar radiation by shadow region that cloud influences and the non-hatched area not influenceed by cloud；S2, the earth's surface solar radiation of the non-hatched area is estimated using clear sky radiance computational methods, the earth's surface solar radiation of the shadow region is estimated using skies Radiation calculation schemes；S3, Combining with terrain bearing calibration, the earth's surface solar radiation of earth's surface solar radiation and the non-hatched area to the shadow region are corrected.The present invention has considered the distortion of cloud shade and influence of the complicated landform to solar radiation, comprehensively and truly simulates situation when solar radiation reaches earth's surface through cloud layer, improves the accuracy using remote sensing technique estimation earth's surface solar radiation.
Description
Technical field
The present invention relates to geographical science field, more particularly, to the earth's surface solar radiation remote sensing for considering cloud-orographic effect
Computational methods and system.
Background technology
Solar radiation is earth's surface physics, chemistry and organismal physiological processes (snow melt, photosynthesis, evapotranspire with plant growth etc.)
Most important energy source, and in earth atmosphere various phenomenons and all physical processes basic motive.Earth's surface sun spoke
Penetrate the energy for controlling ground vapour system and flux exchange process, be cause spatial surface heterogeneous and bioprocess it is crucial because
Son, the significance all having to climatic prediction, Solar use.Because cloud layer has transmission, absorption to solar shortwave radiation
And reflex, make the solar radiation for being incident to earth's surface weakened.Cloud shade is in ground mulching situation and rolling topography
Change, add the uncertainty to earth's surface solar radiation, cloud and influence of the landform to surface radiation is turned into scientist pass
The focus of note.
During solar radiation reaches earth's surface, the decay of air and cloud layer can be passed through, wherein about 31% is reflected and dissipated
Space is emitted back towards, 24% is directly absorbed by the steam in air and ozone etc., and 45% reaches earth's surface.Reach earth's surface solar radiation by
In the masking of landform and the interaction of neighbouring landform, its spatial distribution on ground is set to have the process of reallocation.The earth's surface sun
In addition to being influenceed by cloud, the gradient, slope aspect, 3D geometric effects and cloud the shade distortion of cloud also to utilize remote sensing technique accurate for radiation
Solar simulating radiation becomes more difficult.
Many researchs or just for not considering hypsography, cloud in the solar radiation under clear sky pattern, or skies algorithm
The distortion of 3D geometric effects and cloud shade under MODEL OVER COMPLEX TOPOGRAPHY, it is impossible to which comprehensive real simulation solar radiation is arrived through cloud layer
Up to situation during earth's surface.Earth's surface cloud shade distorts and terrain masking is that the main of remote sensing solar radiation is influenceed under the conditions of rolling topography
Factor.The estimation research of some radiation is only limitted to horizontal earth's surface although it is contemplated that the 3D effects of cloud, or cloud shade is being risen and fallen
The change of earth's surface, which lacks, to be considered.Therefore, the distortion of cloud shade and terrain masking cause remote sensing solar radiation as MODEL OVER COMPLEX TOPOGRAPHY
The main source of deviation, but rarely study studied from the two problem angles while start with present.
The content of the invention
In order to overcome prior art not consider the distortion of cloud shade when estimating earth's surface solar radiation using remote sensing technique
The influence brought with terrain masking to solar radiation, the present invention provide the earth's surface solar radiation remote sensing meter for considering cloud-orographic effect
Calculate method and system.
According to an aspect of the present invention, there is provided a kind of earth's surface solar radiation remote sensing computational methods, including：
S1, using the recognizer of rolling topography cloud shadow region, obtain cloud of the goal pels under MODEL OVER COMPLEX TOPOGRAPHY
Shade, earth's surface solar radiation is distinguished by shadow region that cloud influences and the non-hatched area not influenceed by cloud；
S2, the earth's surface solar radiation of the non-hatched area is estimated using clear sky radiance computational methods, is radiated using the skies
Computational methods estimate the earth's surface solar radiation of the shadow region；
S3, Combining with terrain bearing calibration, the ground of earth's surface solar radiation and the non-hatched area to the shadow region
Table solar radiation is corrected.
Wherein, the earth's surface solar radiation of the non-hatched area is estimated described in step S2 using clear sky radiance computational methods
The step of further comprise：
The Radiation Transmittance of first band and second band is obtained, calculates the earth's surface of the non-hatched area under the conditions of fine day
Direct solar radiation；
By the earth's surface solar scattered radiation of the non-hatched area under the conditions of clear sky be divided into directly scattering radiation and due to
Multiple reflections between ground vapour and caused two parts of scattering radiation are respectively calculated.
Wherein, the earth's surface solar radiation of the shadow region is estimated using skies Radiation calculation schemes described in step S2
Step further comprises：
Using the relation between 0.63-0.69 mu m wavebands reflectivity and cloud optical thickness, high-resolution is obtained by interpolation method
The cloud optical thickness of rate remote sensing image；
Calculate the Radiation Transmittance of first band and second band cloud respectively according to the cloud optical thickness；
The earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud.
Wherein, step S3 further comprises：
According to the appraising model of domatic total solar radiation, the direct radiation from domatic upper reception, sky radiation respectively
With the ground of earth's surface solar radiation and the non-hatched area of three aspects of spurious radiation adjacent to landform to the shadow region
Table solar radiation is corrected.
Wherein, the earth's surface solar scattered radiation by the non-hatched area under the conditions of clear sky is divided into directly scattering spoke
Penetrate and caused two parts of scattering radiation the step of being respectively calculated includes due to the multiple reflections between ground vapour：
Directly scattering radiation is calculated according to following formula：
Edpi=ToiTgiTwsi[BRi(1-Tri)·Tai 0.25]+[Ba·Fi·Tri(1-Tasi 0.25)]Ioicos(Z)
In above formula, EdpiIt is direct scattering radiation, Fi is an experience constant, BRiBefore the Rayleigh scattering for representing i wave bands
To dispersion factor, Ba is the aerosol forward scattering factor, TasiFor the scattering transmitance of i wave bands, TwsiIt is saturating for the scattering of i wave bands steam
Cross rate, ToiFor the Radiation Transmittance of ozone, TgiFor the Radiation Transmittance of mixed gas, TriFor rayleigh scattering transmitance, TaiFor gas
The Radiation Transmittance of colloidal sol, IoiFor the solar constant of i wave bands, cos (Z) is solar zenith angle cosine；
Between calculating ground vapour according to following formula radiation is scattered caused by multiple reflections：
Eddi=ρgiρsi[Ii(hor)+Edpi]/(1-ρgiρsi)
Wherein, EddiRadiation, ρ are scattered caused by multiple reflections between ground vapourgiFor surface albedo, ρsiFor the sky reflection of light
Rate, Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, EdpiIt is direct scattering radiation.
Wherein, the earth's surface sun spoke that the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud
The step of penetrating further comprises：
The earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to following formula,
Ii(hor)=IoiT(cloud)i T(clear)i Ercos(Z)
In above formula, T(cloud)iRepresent the cloud transmitance of i wave bands, Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, Ioi
For the solar constant of i wave bands, T(clear)iRepresent the Radiation Transmittance of i wave bands, ErFor solar distance correction coefficient, cos (Z) is too
Positive zenith angle cosine.
According to another aspect of the present invention, there is provided a kind of earth's surface solar radiation remote sensing computing system, including：
Cloud Shadow recognition unit, for the recognizer using rolling topography cloud shadow region, goal pels are obtained multiple
Cloud shade under miscellaneous orographic condition, earth's surface solar radiation is distinguished by shadow region that cloud influences and the non-shadow not influenceed by cloud
Region；
Solar radiation unit, for estimating the earth's surface sun spoke of the non-hatched area using clear sky radiance computational methods
Penetrate, the earth's surface solar radiation of the shadow region is estimated using skies Radiation calculation schemes；
Topographical correction unit, for Combining with terrain bearing calibration, earth's surface solar radiation to the shadow region and described
The earth's surface solar radiation of non-hatched area is corrected.
Wherein, the solar radiation unit is specifically used for：
The Radiation Transmittance of first band and second band is obtained, calculates the earth's surface of the non-hatched area under the conditions of fine day
Direct solar radiation；
By the earth's surface solar scattered radiation of the non-hatched area under the conditions of clear sky be divided into directly scattering radiation and due to
Multiple reflections between ground vapour and caused two parts of scattering radiation are respectively calculated.
Wherein, the solar radiation unit also particularly useful for：
Using the relation between 0.63-0.69 mu m wavebands reflectivity and cloud optical thickness, high-resolution is obtained by interpolation method
The cloud optical thickness of rate remote sensing image；
Calculate the Radiation Transmittance of first band and second band cloud respectively according to the cloud optical thickness；
The earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud.
Wherein, the topographical correction unit is specifically used for：
According to the appraising model of domatic total solar radiation, the direct radiation from domatic upper reception, sky radiation respectively
With the ground of earth's surface solar radiation and the non-hatched area of three aspects of spurious radiation adjacent to landform to the shadow region
Table solar radiation is corrected.
The earth's surface solar radiation remote sensing computational methods and system provided by the invention for considering cloud-orographic effect, utilize fluctuating
The computational methods of orographic cloud shadow region distinguish actual shadow region and nonshaded area, with reference to earth's surface under the conditions of the skies and clear sky
The earth's surface solar radiation of solar radiation evaluation method and landform correcting method of lane to cloud and rolling topography is estimated, is considered
The distortion of cloud shade and influence of the complicated landform to solar radiation, comprehensively and truly simulate solar radiation and reach in ground through cloud layer
Situation during table, improve the accuracy using remote sensing technique estimation earth's surface solar radiation.
Brief description of the drawings
Fig. 1 is the earth's surface solar radiation remote sensing computational methods for consideration cloud-orographic effect that one embodiment of the invention provides
Schematic flow sheet；
Fig. 2 is the schematic diagram of total incident radiation that earth's surface receives under MODEL OVER COMPLEX TOPOGRAPHY；
Fig. 3 calculates system by the earth's surface solar radiation remote sensing for consideration cloud-orographic effect that another embodiment of the present invention provides
The structural representation of system；
Fig. 4 is cloud and landform condition diagram under certain simulated conditions；
Fig. 5 is by the earth's surface solar shortwave radiation distribution map that is obtained under algorithms of different；
Fig. 5 (a) is the earth's surface solar shortwave radiation distribution map for not doing any correction；
Fig. 5 (b) is to correct but do not consider the earth's surface solar shortwave radiation distribution map of orographic factor apparent place putting into racking；
Fig. 5 (c) is to correct and consider that the earth's surface solar radiation of orographic factor is distributed apparent place putting into racking；
Fig. 5 (d) be based on Fig. 5 (b) carry out topographical correction obtained by the earth's surface for not considering cloud shade landform distortion effects too
Positive shortwave radiation distribution map；
Fig. 5 (e) is the earth's surface sun that the consideration cloud shade landform distortion effects obtained by topographical correction are carried out based on Fig. 5 (c)
Shortwave radiation distribution map；
Fig. 5 (f) is Fig. 5 (d) and Fig. 5 (e) differential chart；
Fig. 6 is the radiation error standard deviation of 20 kinds of different terrain complexities and the relation of terrain complexity；
Fig. 7 is the radiation value comparison diagram for having the tailo sampling point of cloud shade under three kinds of algorithm situations；
Fig. 8 is the radiation value comparison diagram for having the Schattenseite sampling point of cloud shade under three kinds of algorithm situations；
Fig. 9 is radiation value comparison diagram of the tailo sampling point of cloudless shade under three kinds of algorithm situations；
Figure 10 is radiation value comparison diagram of the Schattenseite sampling point of cloudless shade under three kinds of algorithm situations；
Figure 11 is the cloud shade distribution map around cajaput website on July 14th, 2014；
Figure 12 is the cloud shade distribution map around cajaput website on July 24th, 2014.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, clear, complete description is carried out to the technical scheme in the embodiment of the present invention, it is clear that described embodiment is only
Only it is part of the embodiment of the present invention, rather than whole embodiments.Based on embodiments of the invention, ordinary skill people
The every other embodiment that member is obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
As described in Figure 1, the earth's surface solar radiation remote sensing of the consideration cloud-orographic effect provided for one embodiment of the invention calculates
The schematic flow sheet of method, including：
S1, using the recognizer of rolling topography cloud shadow region, obtain cloud of the goal pels under MODEL OVER COMPLEX TOPOGRAPHY
Shade, earth's surface solar radiation is distinguished by shadow region that cloud influences and the non-hatched area not influenceed by cloud；
Earth's surface cloud shade distorts and terrain masking is the principal element of influence remote sensing solar radiation under the conditions of rolling topography, because
This, the earth's surface solar radiation remote sensing computational methods that the embodiment of the present invention is provided, obtains cloud the moon under MODEL OVER COMPLEX TOPOGRAPHY first
Shadow, earth's surface solar radiation is distinguished by region that cloud influences and the region not influenceed by cloud.Wherein, the fluctuating described in step S1
The recognizer of orographic cloud shadow region, specifically includes following steps：Any cloud pixel in cloud atlas picture is chosen as target picture
Member, obtain the position of shade of the goal pels in the upright projection position, the goal pels of horizontal earth's surface in horizontal earth's surface
With the height value of the goal pels, and the linear equation of light is established；The altitude data of cloud atlas earth's surface as corresponding to is obtained, will
The altitude data is transformed under space coordinates；According to the coordinate of each point, height value and the light in the altitude data
Linear equation obtain match point coordinate, the coordinate of the match point is the position of goal pels shadow spots；Obtain cloud atlas
The position of goal pels shadow spots as in, so as to obtain the position of cloud shade.The identification of above-mentioned rolling topography cloud shadow region is calculated
Method is in patent " a kind of cloud shadow detection method and system (number of patent application：201510792310.6) " in be documented.
Behind the position for detecting cloud shade, you can using by region division of the cloud atlas as residing for medium cloud shadow positions as shadow region, cloud atlas
Then incorporate into as non-hatched area in other regions as in.
S2, the earth's surface solar radiation of the non-hatched area is estimated using clear sky radiance computational methods, is radiated using the skies
Computational methods estimate the earth's surface solar radiation of the shadow region；
1) estimated for the surface radiation of non-hatched area, using clear sky radiance computational methods, because rolling topography is to straight
It is different with the Influencing Mechanism of scattering radiation to connect radiation, thus the earth's surface solar radiation calculating under the conditions of fine day is needed respectively for straight
Connect radiation and scattering radiation individually calculates.
The step of estimating the earth's surface solar radiation of the non-hatched area using clear sky radiance computational methods specifically includes：
S21a, the Radiation Transmittance of first band and second band is obtained, calculate the non-hatched area under the conditions of fine day
Earth's surface direct solar radiation；
The directly radiation of horizontal earth's surface can be calculated with formula (1) under the conditions of fine day,
Ii(hor)=IoiT(clear)iErcos(Z) (1)
In formula (1), Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, IoiFor the solar constant of i wave bands, T(clear)i
Represent the Radiation Transmittance of i wave bands, ErFor solar distance correction coefficient, cos (Z) is solar zenith angle cosine.I value is ripple
Section 1 and wave band 2, first band refer to wave band 1, represent 0.29-0.7 μm of visible light wave range, and second band refers to wave band 2, represent
0.7-4 μm of near infrared band, typically by Io1Value is 0.038MJm-2·min-1, Io2Value is 0.043MJm-2·min-1。
Solar radiation, which reaches earth's surface, under the conditions of fine day will also pass through a series of absorption, scattering, wherein mainly including ozone
Absorption, Rayleigh scattering, mixed gas absorb, the absorption and scattering of water vapor absorption and aerosol.It is i as shown in formula (2)
The Radiation Transmittance T of wave band(clear)iCalculation formula, subscript o, r, g, w, a in formula represent ozone, Rayleigh, mixing respectively
The Radiation Transmittance of gas, steam and aerosol.
T(clear)i=ToiTriTgiTwiTai (2)
Due to ozone and the more difficult acquisition of aerosol related data, thus we using some models come simulate correlation gas it is molten
Glue and ozone parameter.It is aerosol model as shown in formula (3), (4) and (5)：
Ua=ln (1+Ma* β) (3)
Ma={ cosZ+ [0.168Z0.18/(95.318-Z)1.95]}-1 (4)
β=(0.025+0.1cos φ) exp (- 0.7alt) (5)
Wherein, UaAerosol load is represented, Ma is aerosol optical quality, and β is Angstrom Turbidity parameters, and Z is too
Positive zenith angle, φ represent latitude, and alt represents height above sea level.
As shown in formula (6) and (7), the content of ozone obtains also by parametric method,
U0=0.44-0.16 { [(φ -80)/60]2+[(d-120)/(9263-φ)2]}0.5 (6)
Wherein, U0Ozone content is represented, d is Julian date, and Jd is Julian date, and φ represents latitude.
On the right of formula (2) equation all transmitances (including：Ozone, Rayleigh, mixed gas, the spoke of steam and aerosol
Penetrate transmitance) calculating all referring to Gueymard in paper " Gueymard C A.A two-band model for the
calculation ofclear sky solar irradiance,illuminance,and photosynthetically
active radiation at the Earth’s surface[J].SolarEnergy,1989,43(5):253-265 " in
The algorithm proposed is calculated.So far, we can calculate the direct radiation obtained under the conditions of clear sky.
S21b, by the earth's surface solar scattered radiation of the non-hatched area under the conditions of clear sky be divided into directly scattering radiation and
Due to the multiple reflections between ground vapour, caused two parts of scattering radiation are respectively calculated.
Earth's surface (surface albedo 0) is being fully absorbed, directly scattering radiation is calculated according to formula (8)：
Edpi=ToiTgiTwsi[BRi(1-Tri)·Tai 0.25]+[Ba·Fi·Tri(1-Tasi 0.25)]Ioicos(Z) (8)
In above formula, EdpiIt is direct scattering radiation, Fi is an experience constant, and correcting the factor as one is introduced into, purpose
It is to make up used simple transmitance algorithm in Multiple Scattering and otherwise deficiency.BRiRepresent the auspicious of i wave bands
The forward scattering factor of profit scattering, Ba is the aerosol forward scattering factor, TasiFor the scattering transmitance of i wave bands, TwsiFor i wave bands
Steam scatters transmitance, ToiFor the Radiation Transmittance of ozone, TgiFor the Radiation Transmittance of mixed gas, TriIt is saturating for rayleigh scattering
Cross rate, TaiFor the Radiation Transmittance of aerosol, IoiFor the solar constant of i wave bands, cos (Z) is solar zenith angle cosine.
Radiation is scattered between ground vapour caused by multiple reflections, can be calculated by formula (9),
Eddi=ρgiρsi[Ii(hor)+Edpi]/(1-ρgiρsi) (9)
In above formula, EddiRadiation, ρ are scattered caused by multiple reflections between ground vapourgiFor surface albedo, ρsiFor the sky reflection of light
Rate, Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, EdpiIt is direct scattering radiation.
Therefore, the total scattering radiation of each wave band can be expressed as Edi=Edpi+Eddi.And the Globalradiation of each wave band is
For Egi=Ii(hor)+Edi.Broadband Globalradiation is then Eg=I (hor)+Ed。
2) for the earth's surface solar radiation of shadow region, estimated using skies Radiation calculation schemes.
Radiation under the conditions of the skies calculates the Radiation Transmittance that cloud is mainly multiplied by the basis of formula (1), such as formula
(10) shown in,
Ii(hor)=IoiT(cloud)i T(clear)i Ercos(Z) (10)
In above formula, T(cloud)iRepresent the cloud transmitance of i wave bands, the same formula of other specification (1).
Understood according to formula (10), the earth's surface solar radiation of the shadow region is estimated using skies Radiation calculation schemes
Step further comprises：
S22a, using the relation between 0.63-0.69 mu m wavebands reflectivity and cloud optical thickness, obtained by interpolation method high
The cloud optical thickness of resolution remote sense image；
Specifically, cloud transmitance T(cloud)iCalculating need to obtain cloud optical thickness parameter first, it is but corresponding due to lacking
High spatial resolution cloud optical thickness product, and existing MODIS (MODerate-resolution Imaging
Spectroradiometer cloud optical thickness product resolution ratio) is relatively low (1km), is differed with actually required high spatial resolution
It is larger, it is more difficult to match.Therefore, we mainly employ the mu m waveband of high spatial resolution satellite 0.63~0.69 (centre wavelength are
0.65 μm) relation between reflectivity and cloud optical thickness, cloud optical thickness product NO emissions reduction is given, then utilizes interpolation to obtain high
The cloud optical thickness of spatial resolution.0.63~0.69 mu m waveband (centre wavelength be 0.65 μm) reflectivity and cloud optical thickness it
Between relation it is as shown in table 1.
The relation of the wave band of 1 HJ stars of table the 3rd (centre wavelength is 0.65 μm) reflectivity and cloud optical thickness
| Reflectivity (0.65 μm) | Cloud optical thickness (COT) |
| 0.15-0.35 | 0-1 |
| 0.35-0.5 | 1-4 |
| 0.5-0.55 | 4-8 |
| 0.55-0.6 | 8-16 |
| 0.6-0.65 | 16-32 |
| 0.65-0.7 | 32-64 |
S22b, calculate the Radiation Transmittance of first band and second band cloud respectively according to the cloud optical thickness；
First band cloud transmitance is calculated using formula (11)：
T(cloud)1=1- [β (μ0)τn/μ0]/[1+β(μ0)τn/μ0] (11)
In above formula, T(cloud)1Represent the cloud transmitance of 1 wave band, μ0It is solar zenith angle cosine, β (μ0) represent it is corresponding after
To scattering ratio, τnFor cloud optical thickness.
The cloud transmitance of second band is then obtained using formula (12) to (15), and formula (12) is as follows to (15)：
T(cloud)2=4u/R (12)
R=(u+1)2exp(τeff)-(u-1)2exp(-τeff) (14)
During above formula is various, T(cloud)2For second band cloud transmitance,For single scattering albedo, μ0It is solar zenith angle
Cosine, β (μ0) represent corresponding Back-scattering ratio, τnFor cloud optical thickness, u, R and τeffAfter being parameterized for two-stream approximation solution
Intermediate variable, relevant parameter can be from Stephens in paper " Stephens G L.Ashortwave
parameterization revised to improve cloud absorption[J].Journal of the
Atmospheric Sciences,1984,41(4):Obtained in the look-up table provided in 687-690 ".
S23c, the earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud.
The Radiation Transmittance that the cloud of the i wave bands of acquisition will be calculated substitutes into the bottom that formula (10) can obtain the shadow region
Portion's solar radiation.
In summary, after the cloud shade under obtaining MODEL OVER COMPLEX TOPOGRAPHY, it is possible to which distinguish surface radiation is influenceed by cloud
Region and the region that is not influenceed by cloud, then use skies Radiation calculation schemes in shadow region, clear sky used in nonshaded area
Radiation calculation schemes, preresearch estimates is carried out to earth's surface solar radiation.
S3, with reference to known landform correcting method of lane, earth's surface solar radiation and the nonshaded area to the shadow region
The earth's surface solar radiation in domain is corrected.
Dubayah etc. is in paper " Dubayah R.Estimating net solar radiation using
Landsat Thematic Mapper and digital elevation data[J].Water resources
research,1992,28(9):Have ignored in 2469-2484 " between earth's surface and air on the basis of multiple reflections, by it is domatic too
Positive global radiation EtIt is divided into the direct radiation E of domatic upper receptiondir,t, sky radiation Edif,tWith the spurious radiation of neighbouring landform
Ea,t, as shown in formula (16),
Et=Edir,t+Edif,t+Ea,t (16)
As shown in Fig. 2 the schematic diagram of the total incident radiation received for earth's surface under MODEL OVER COMPLEX TOPOGRAPHY.From radiation transmission and
The angle of illumination geometrical principle is set out, and influence of the landform to surface radiation is broadly divided into following three aspects：First, directly
The influence of radiation, the difference of gradient slope aspect change domatic with sunray geometrical relationship, form domatic itself and periphery
Bridging effect, change the sunshine-duration.2nd, the influence of sky scattering, gradient slope aspect cause day a certain scattered particle, the sun in the air
Position and the triangular relative position of earth's surface goal pels change, and then change scattering light and reach earth's surface goal pels
Incidence angle, also have impact on sky shield coverage of the periphery landform to goal pels；3rd, the influence radiated to neighbouring landform, by
Geometrical relationship and its mutually masking of the goal pels with neighbouring pixel, and near reflection pixel reflected radiation BRDF are (two-way anti-
Penetrate distribution function) influence of characteristic.
Develop as GIS is calculated, the topographical correction of earth's surface solar radiation starts to rely on dem data, from three above component
Horizontal direction carry out, i.e., according to the appraising model of domatic total solar radiation, the direct radiation from domatic upper reception, sky respectively
Earth's surface solar radiation and the non-shadow of three aspects of spurious radiation of scattering radiation and neighbouring landform to the shadow region
The earth's surface solar radiation in region is corrected.
Specifically, the topographical correction model 1) directly radiated combines Dubayah in paper " Dubayah
R.Estimating net solar radiation using Landsat Thematic Mapper and digital
elevation data[J].Water resources research,1992,28(9):2469-2484 " and Zeng Yan are in paper
" Zeng Yan, Qiu Xinfa, Liu Changming, Huanghe valley direct solar radiation distributed simulation [J] Geography Journals under rolling topographies are waited,
2005,60(4):680-688 and Zeng Yan, Qiu Xinfa, Pan Ao great, the analysis that wait landform influences on Huanghe valley solar radiation are ground
Study carefully [J] Advances in Earth Sciences, 2008,23 (11):The model provided in 1185-1193 ", it is contemplated that goal pels are in itself and adjacent
The masking of both nearly pixels influences, as shown in formula (17),
In above formula, Edir,hIt is directly to radiate, I represents incident angle, if I is more than 90 expression pixels and covered by itself, L is
For judging the binary coefficient of projection, L=0 represents that goal pels are covered by neighbouring pixel, and L=1 then represents not shielded.
2) correction of scattering radiation is mainly according to Hay in paper " Hay J E, Mckay D C.Estimating solar
irradiance on inclined surfaces:a review and assessment of methodologies[J]
.International Journal of Solar Energy.1985,3(4):The anisotropy mould proposed in 203-240 "
Type, by separating for the anisotropic segment of sky scattering and isotropism part independence, with anisotropic index k
(Anisotropy Index) represents that anisotropic scattering accounts for the weight of sky total scattering, and k is received too with domatic normal direction
The ratio between positive direct projection irradiation level and the radiation of atmosphere top calculate.Scattering radiation is corrected using formula (18) and (19).
Rd=k+ (1-k) [(1+cosS)/2] (18)
In above formula, RdRadiated for scattering, S is the domatic gradient, ShIt is sun altitude, Edir,hDirectly to radiate, the sun
Constant Io=1367w/m2。
3) for neighbouring landform spurious radiation calculating Main Basiss Proy paper " Proy C, Tanre D,
Deschamps P Y.Evaluation of topographic effects in remotely sensed data[J]
.Remote Sensing of Environment,1989,30(1):The theory carried in 21-32 ".
To sum up, step S3 completes estimates the correction based on landform to earth's surface solar radiation.
The embodiment of the present invention distinguishes actual shadow region and non-the moon using the computational methods of rolling topography cloud shadow region
Shadow zone, with reference to earth's surface solar radiation evaluation method under the conditions of the skies and clear sky and landform correcting method of lane to cloud and rolling topography
Earth's surface solar radiation is estimated, has considered the distortion of cloud shade and influence of the complicated landform to solar radiation, comprehensively true
Simulation solar radiation reaches situation during earth's surface through cloud layer, improves and utilizes remote sensing technique estimation earth's surface solar radiation
Accuracy.
As shown in figure 3, the earth's surface solar radiation remote sensing of the consideration cloud-orographic effect provided by another embodiment of the present invention
The structural representation of computing system, including：Cloud Shadow recognition unit 31, solar radiation unit 32 and topographical correction unit 33, its
In,
Cloud Shadow recognition unit 31, for the recognizer using rolling topography cloud shadow region, obtain goal pels and exist
Cloud shade under MODEL OVER COMPLEX TOPOGRAPHY, earth's surface solar radiation is distinguished by shadow region that cloud influences and the non-the moon not influenceed by cloud
Shadow zone domain；
Earth's surface cloud shade distorts and terrain masking is the principal element of influence remote sensing solar radiation under the conditions of rolling topography, because
This, the earth's surface solar radiation remote sensing computing system that the embodiment of the present invention is provided includes cloud Shadow recognition unit 31, for obtaining
Cloud shade under MODEL OVER COMPLEX TOPOGRAPHY, earth's surface solar radiation is distinguished by region that cloud influences and the region not influenceed by cloud.Its
In, the recognizer of described rolling topography cloud shadow region, specifically include following steps：Choose any cloud picture in cloud atlas picture
Member, obtain position and the institute of shade of the cloud pixel in the upright projection position of horizontal earth's surface, the cloud pixel in horizontal earth's surface
The height value of cloud pixel is stated, and establishes the linear equation of light；The altitude data of cloud atlas earth's surface as corresponding to is obtained, by the height
Number of passes evidence is transformed under space coordinates；According to the straight line of the coordinate of each point, height value and the light in the altitude data
Equation obtains the coordinate of match point, and the coordinate of the match point is the position of cloud pixel shadow spots；Cloud atlas is obtained as medium cloud picture
The position of first shadow spots, so as to obtain the position of cloud shade.The recognizer of above-mentioned rolling topography cloud shadow region is in patent " one
Kind cloud shadow detection method and system (number of patent application：201510792310.6) " in be documented.Detecting cloud shade
Position after, you can using by region division of the cloud atlas as residing for medium cloud shadow positions as shadow region, other areas in cloud atlas picture
Then incorporate into as non-hatched area in domain.
Solar radiation unit 32, for estimating the earth's surface sun spoke of the non-hatched area using clear sky radiance computational methods
Penetrate, the earth's surface solar radiation of the shadow region is estimated using skies Radiation calculation schemes；
1) estimate that solar radiation unit 32 uses clear sky radiance computational methods for the surface radiation of non-hatched area, by
It is different with the Influencing Mechanism of scattering radiation to directly radiation in rolling topography, thus the earth's surface solar radiation under the conditions of fine day calculates
Need individually to calculate for directly radiation and scattering radiation respectively.
Therefore solar radiation unit 32 is specifically used for：
S21a, the Radiation Transmittance of first band and second band is obtained, calculate the non-hatched area under the conditions of fine day
Earth's surface direct solar radiation；
The directly radiation of horizontal earth's surface can be calculated with formula (1) under the conditions of fine day,
Ii(hor)=IoiT(clear)iErcos(Z) (1)
In formula (1), Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, IoiFor the solar constant of i wave bands, T(clear)i
Represent the Radiation Transmittance of i wave bands, ErFor solar distance correction coefficient, cos (Z) is solar zenith angle cosine.I value is ripple
Section 1 and wave band 2, wave band 1 represent 0.29-0.7 μm of visible light wave range, and wave band 2 represents 0.7-4 μm of near infrared band, typically by Io1
Value is 0.038MJm-2·min-1, Io2Value is 0.043MJm-2·min-1。
Solar radiation, which reaches earth's surface, under the conditions of fine day will also pass through a series of absorption, scattering, wherein mainly including ozone
Absorption, Rayleigh scattering, mixed gas absorb, the absorption and scattering of water vapor absorption and aerosol.It is i as shown in formula (2)
The Radiation Transmittance T of wave band(clear)iCalculation formula, subscript o, r, g, w, a in formula represent ozone, Rayleigh, mixing respectively
The Radiation Transmittance of gas, steam and aerosol.
T(clear)i=ToiTriTgiTwiTai (2)
So far, we can calculate the direct radiation obtained under the conditions of clear sky.
S21b, by the earth's surface solar scattered radiation of the non-hatched area under the conditions of clear sky be divided into directly scattering radiation and
Due to the multiple reflections between ground vapour, caused two parts of scattering radiation are respectively calculated.
Earth's surface (surface albedo 0) is being fully absorbed, directly scattering radiation is calculated according to formula (8)：
Edpi=ToiTgiTwsi[BRi(1-Tri)·Tai 0.25]+[Ba·Fi·Tri(1-Tasi 0.25)]Ioicos(Z) (8)
In above formula, EdpiIt is direct scattering radiation, Fi is an experience constant, and correcting the factor as one is introduced into, purpose
It is to make up used simple transmitance algorithm in Multiple Scattering and otherwise deficiency.BRiRepresent the auspicious of i wave bands
The forward scattering factor of profit scattering, Ba is the aerosol forward scattering factor, TasiFor the scattering transmitance of i wave bands, TwsiFor i wave bands
Steam scatters transmitance, ToiFor the Radiation Transmittance of ozone, TgiFor the Radiation Transmittance of mixed gas, TriIt is saturating for rayleigh scattering
Cross rate, IoiFor the solar constant of i wave bands, cos (Z) is solar zenith angle cosine.
Radiation is scattered between ground vapour caused by multiple reflections, can be calculated by formula (9),
Eddi=ρgiρsi[Ii(hor)+Edpi]/(1-ρgiρsi) (9)
In above formula, EddiRadiation, ρ are scattered caused by multiple reflections between ground vapourgiFor surface albedo, ρsiFor the sky reflection of light
Rate, Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, EdpiIt is direct scattering radiation.
Therefore, the total scattering radiation of each wave band can be expressed as Edi=Edpi+Eddi.And the Globalradiation of each wave band is
For Egi=Ii(hor)+Edi.Broadband Globalradiation is then Eg=I (hor)+Ed。
2) estimated for the earth's surface solar radiation of shadow region, solar radiation unit 32 using skies Radiation calculation schemes
Calculate.
Radiation under the conditions of the skies calculates the Radiation Transmittance that cloud is mainly multiplied by the basis of formula (1), such as formula
(10) shown in,
Ii(hor)=IoiT(cloud)i T(clear)i Ercos(Z) (10)
In above formula, T(cloud)iRepresent the cloud transmitance of i wave bands, the same formula of other specification (1).
Understand that solar radiation unit 32 is specifically used for according to formula (10)：
Using the relation between 0.63-0.69 mu m wavebands reflectivity and cloud optical thickness, high-resolution is obtained by interpolation method
The cloud optical thickness of rate remote sensing image；
Employ the mu m waveband of high spatial resolution satellite 0.63~0.69 (centre wavelength is 0.65 μm) reflectivity and cloud
Relation between optical thickness, cloud optical thickness product NO emissions reduction is given, the cloud light of high spatial resolution is then obtained using interpolation
Learn thickness.
Calculate the Radiation Transmittance of first band and second band cloud respectively according to the cloud optical thickness；
First band cloud transmitance is calculated using formula (11)：
T(cloud)1=1- [β (μ0)τn/μ0]/[1+β(μ0)τn/μ0] (11)
In above formula, T(cloud)1Represent the cloud transmitance of 1 wave band, μ0It is solar zenith angle cosine, β (μ0) represent it is corresponding after
To scattering ratio, τnFor cloud optical thickness.
The cloud transmitance of second band is then obtained using formula (12) to (15), and formula (12) is as follows to (15)：
T(cloud)2=4u/R (12)
R=(u+1)2exp(τeff)-(u-1)2exp(-τeff) (14)
During above formula is various, T(cloud)2For second band cloud transmitance,For single scattering albedo, μ0It is solar zenith angle
Cosine, β (μ0) represent corresponding Back-scattering ratio, τnFor cloud optical thickness, u, R and τeffAfter being parameterized for two-stream approximation solution
Intermediate variable.
Then, the earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud.
The Radiation Transmittance that the cloud of the i wave bands of acquisition will be calculated substitutes into the bottom that formula (10) can obtain the shadow region
Portion's solar radiation.
Topographical correction unit 33, for Combining with terrain bearing calibration, earth's surface solar radiation and institute to the shadow region
The earth's surface solar radiation for stating non-hatched area is corrected.
Wherein, the topographical correction unit is specifically used for：
According to the appraising model of domatic total solar radiation, the direct radiation from domatic upper reception, sky radiation respectively
With the ground of earth's surface solar radiation and the non-hatched area of three aspects of spurious radiation adjacent to landform to the shadow region
Table solar radiation is corrected.
The earth's surface solar radiation remote sensing computing system for consideration cloud-orographic effect that the embodiment of the present invention proposes, utilizes fluctuating
The computational methods of orographic cloud shadow region distinguish actual shadow region and nonshaded area, with reference to earth's surface under the conditions of the skies and clear sky
The earth's surface solar radiation of solar radiation evaluation method and landform correcting method of lane to cloud and rolling topography is estimated, is considered
The distortion of cloud shade and influence of the complicated landform to solar radiation, comprehensively and truly simulate solar radiation and reach in ground through cloud layer
Situation during table, improve the accuracy using remote sensing technique estimation earth's surface solar radiation.
With reference to specific simulated experiment, HJ satellite datas are based on to method proposed by the invention and carried out further
Analytic explanation.As shown in table 2, it is the domain load parameter of HJ-1B satellites.
The HJ-1B satellite domain load parameters of table 2
Earth's surface solar radiation by cloud, the gradient, slope aspect in addition to being influenceed, with the change of the relative position of the sun and cloud and its 3D
Geometric effect, the masking of distortion (including position and pixel number) and landform of the cloud shade in complicated landform etc., which turns into, utilizes remote sensing
The important factor in order of method solar simulating radiation exactly.Carried out regarding to the issue above in earth's surface solar radiation refutation process
A series of correction, analysis is compared for the result for correcting front and rear.
The global radiation change of three kinds of algorithms of Main Analysis and difference, three kinds of algorithms do not do any correction traditional algorithm,
Enter to rack 3D geometric corrections or radiation topographical correction and the algorithm for considering cloud-topographical correction proposed by the invention.
Fig. 4 is cloud and the landform situation under certain simulated conditions, and Fig. 5 is by the earth's surface sun shortwave spoke that is obtained under algorithms of different
Penetrate distribution map, wherein, Fig. 5 (a) is not do the earth's surface solar shortwave radiation distribution map of any correction, Fig. 5 (b) be into racking apparent place
Put correction but do not consider the earth's surface solar shortwave radiation distribution map of orographic factor, Fig. 5 (c) is to correct and consider apparent place putting into racking
Orographic factor earth's surface solar radiation distribution, Fig. 5 (d) be based on Fig. 5 (b) carry out topographical correction obtained by with not considering cloud shade
The earth's surface solar shortwave radiation distribution map of shape distortion effects, Fig. 5 (e) are that the consideration cloud obtained by topographical correction is carried out based on Fig. 5 (c)
The earth's surface solar shortwave radiation distribution map of shade landform distortion effects, Fig. 5 (f) are Fig. 5 (d) and Fig. 5 (e) differential chart.
By Fig. 5 (a) compared with Fig. 5 (b), it can be seen that both have a larger difference in radiation profiles, and several places are obvious
Difference be marked respectively on two figures with black circles.
By Fig. 5 (b) compared with Fig. 5 (c), if it can be found that not doing any correction or only cloud three-dimensional geometry effect
It should correct without considering that landform distorts to caused by cloud shade and can cause relatively large deviation to radiative invesion.As originally not by cloud shade
The region of covering can be by the region that is covered by cloud shade of can be regarded as of mistake, and vice versa, black circles can be marked from two figures
A few place's differences find out.
By Fig. 5 (d) compared with Fig. 5 (e), it can be seen that whether consider change of the cloud shade in complicated landform to spoke
The result for penetrating inverting has a major impact.Radiation has significant difference in the region for having cloud shadow-casting, particularly in some landform
Complexity, the larger region of unduation grade.For the ease of analyzing radiation change therebetween and difference, both are subtracted each other, tied
Shown in fruit such as Fig. 5 (f).From Fig. 5 (f) it can be seen that the area radiation value being distorted in cloud shade has different degrees of difference
It is different.The place to differ greatly is usually the tailo and the Schattenseite for having cloud shade in the cloudless shade of complicated earth surface, and difference is lessly
It is square then be usually that the position of shade is changed, and does not cover the phenomenon of superposition same domatic.
Distortion of the cloud shade in landform causes earth's surface solar radiation to produce error, the standard deviation of radiation error in distribution
The assembly average (absolute value) of the radiation error caused by the distortion of shade landform can be represented not consider in whole region.Radiation
The standard deviation of error can be expressed as：
In above formula, Dti,jTo consider point i, j of the distortion of shade landform radiation value, Di,jNot consider that shade landform distorts
Point i, j radiation value, m, n are the line number and columns of the lattice point in whole region respectively.Fig. 6 is 20 kinds of different terrain complexities
Radiation error standard deviation and terrain complexity relation.It can be seen that radiation error increases with the increase of terrain complexity,
This deviation may be bigger when solar zenith angle is larger, data resolution is higher.Therefore, provided in application high-resolution satellite
When material calculates Land surface energy budget each component in complicated landform, it is necessary to consider the effect of landform.
It is instantaneous when on June 12nd, 2012, July 4,8, August data on the 3rd, 11 are simulated to obtain with the morning 11 or so
Global radiation, sampling point is chosen with reference to the visual interpretation to actual shadow region and nonshaded area and logic judgment, being divided into has cloud shade
Tailo and Schattenseite and cloudless shade 4 kinds of tailo and Schattenseite.3 kinds of algorithms of Main Analysis, algorithm 1 are not make any correction
Traditional algorithm, algorithm 2 are the method that only considered cloud 3D effects or radiate orographic effect, and algorithm 3 is method proposed by the present invention.
As shown in fig. 7, to there is radiation value comparison diagram of the tailo sampling point of cloud shade under three kinds of algorithm situations.It can be seen that
It is obvious higher compared with inventive algorithm that the result of any correction is not done, reaches amount of radiation 1000-1100W/m during fine day2,
Average deviation about 400W/m2.And do not consider that the result of shade distortion is also higher, due to being in tailo, so after topographical correction
Radiation value is more bigger than normal, and average deviation is about 500W/m2。
As shown in figure 8, to there is radiation value comparison diagram of the Schattenseite sampling point of cloud shade under three kinds of algorithm situations.It is most of not
Consider the result of shade distortion has different degrees of higher, maximum difference 500W/m compared with inventive algorithm2, average deviation
About 300W/m2.And the result for not doing any correction is also mostly higher, some reaches amount of radiation 1000W/m during fine day2More than, by
In in Schattenseite, so the radiation value after topographical correction is more relatively low, deviation further increases, averagely about 400W/m2。
As shown in figure 9, radiation value comparison diagram of the tailo sampling point under three kinds of algorithm situations for cloudless shade.It is most of not
Consider the result of shade distortion has different degrees of relatively low, maximum difference 650W/m compared with inventive algorithm2, average deviation
About 400W/m2。
As shown in Figure 10, it is radiation value comparison diagram of the Schattenseite sampling point of cloudless shade under three kinds of algorithm situations.Do not consider
The result of shade distortion has different degrees of less than normal, maximum difference 480W/m compared with inventive algorithm2, average deviation is about
300W/m2.And the result for not doing any correction is also mostly higher, because the radiation value after Schattenseite, topographical correction is more
Relatively low, deviation further increases, averagely about 350W/m2。
Whether be corrected in the region that some clouds and shade have superposition the deviation for research at this stage be not it is very big,
Typically in 100W/m2Left and right.
Figure 11 and Figure 12 is respectively that (red point is cajaput in figure for the cloud shade distribution map on July 14th, 2014 and July 24
Stand).Due to lacking the transient data of website, while remotely-sensed data that can be illustratively is extremely limited.Its day is estimated with two
Value radiation, with station data comparative analysis.
If data on July 14 can think website position cloudless covering according to traditional algorithm, using whole day to be fine
My god, it is assumed that calculate instantaneous radiation value and reach maximum (when about 11), radiation diurnal variation meets normal distribution, estimates that the radiation of its earning in a day is
33.9MJ/m2, it is 26.2MJ/m according to day radiation obtained by inventive algorithm2, and station data is 24.5MJ/m2。
If data on July 24 can think website position have cloud covering according to traditional algorithm, with whole day (14
Hour) radiated for the skies, estimate that the radiation of its earning in a day is 26.2MJ/m2, it is 30.6MJ/m according to day radiation obtained by inventive algorithm2,
And station data is 29.5MJ/m2。
Although the estimation of earning in a day radiation has certain error, analyzed more than it can be found that inventive algorithm is in cloud
There is advantage compared to traditional algorithm under it and MODEL OVER COMPLEX TOPOGRAPHY.
Finally, the various embodiments described above of the present invention are only preferable embodiment, are not intended to limit the protection model of the present invention
Enclose.Within the spirit and principles of the invention, any modification, equivalent substitution and improvements made etc., should be included in the present invention
Protection domain within.
Claims (10)
1. consider the earth's surface solar radiation remote sensing computational methods of cloud-orographic effect, it is characterised in that including：
S1, using the recognizer of complicated landform cloud shadow region, cloud shade of the goal pels under MODEL OVER COMPLEX TOPOGRAPHY is obtained,
Earth's surface solar radiation is distinguished by shadow region that cloud influences and the non-hatched area not influenceed by cloud；
S2, the earth's surface solar radiation of the non-hatched area is estimated using clear sky radiance computational methods, radiated and calculated using the skies
Method estimates the earth's surface solar radiation of the shadow region；
S3, Combining with terrain bearing calibration, the earth's surface of earth's surface solar radiation and the non-hatched area to the shadow region is too
Sun radiation is corrected.
2. according to the method for claim 1, it is characterised in that estimated described in step S2 using clear sky radiance computational methods
The step of earth's surface solar radiation of the non-hatched area, further comprises：
The Radiation Transmittance of first band and second band is obtained, calculates the earth's surface sun of the non-hatched area under the conditions of fine day
Directly radiate；
The earth's surface solar scattered radiation of the non-hatched area under the conditions of clear sky is divided into directly scattering radiation and due to ground vapour
Between multiple reflections and it is caused scattering radiation two parts be respectively calculated.
3. according to the method for claim 1, it is characterised in that estimated described in step S2 using skies Radiation calculation schemes
The step of earth's surface solar radiation of the shadow region, further comprises：
Using the relation between 0.63-0.69 mu m wavebands reflectivity and cloud optical thickness, it is distant that high-resolution is obtained by interpolation method
Feel the cloud optical thickness of image；
Calculate the Radiation Transmittance of first band and second band cloud respectively according to the cloud optical thickness；
The earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud.
4. according to the method for claim 1, it is characterised in that step S3 further comprises：
According to the appraising model of domatic total solar radiation, respectively from direct radiation, sky radiation and the neighbour of domatic upper reception
The earth's surface of earth's surface solar radiation and the non-hatched area of three aspects of spurious radiation of near-earth shape to the shadow region is too
Sun radiation is corrected.
5. according to the method for claim 2, it is characterised in that the ground by the non-hatched area under the conditions of clear sky
Table solar scattered radiation is divided into directly scattering radiation and caused two parts of scattering radiation due to the multiple reflections between ground vapour
The step of being respectively calculated includes：
Directly scattering radiation is calculated according to following formula：
Edpi=ToiTgiTwsi[BRi(1-Tri)·Tai 0.25]+[Ba·Fi·Tri(1-Tasi 0.25)]Ioicos(Z)
In above formula, EdpiIt is direct scattering radiation, Fi is an experience constant, BRiRepresent the forward scattering of the Rayleigh scattering of i wave bands
The factor, Ba are the aerosol forward scattering factor, TasiFor the scattering transmitance of i wave bands, TwsiTransmitance is scattered for i wave bands steam,
ToiFor the Radiation Transmittance of ozone, TgiFor the Radiation Transmittance of mixed gas, TriFor rayleigh scattering transmitance, TaiFor aerosol
Radiation Transmittance, IoiFor the solar constant of i wave bands, cos (Z) is solar zenith angle cosine；
Between calculating ground vapour according to following formula radiation is scattered caused by multiple reflections：
Eddi=ρgiρsi[Ii(hor)+Edpi]/(1-ρgiρsi)
Wherein, EddiRadiation, ρ are scattered caused by multiple reflections between ground vapourgiFor surface albedo, ρsiFor sky albedo, Ii
(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, EdpiIt is direct scattering radiation.
6. according to the method for claim 3, it is characterised in that described that skies bar is calculated according to the Radiation Transmittance of the cloud
Further comprise under part the step of the earth's surface solar radiation of the shadow region：
The earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to following formula,
Ii(hor)=IoiT(cloud)i T(clear)i Ercos(Z)
In above formula, T(cloud)iRepresent the cloud transmitance of i wave bands, Ii(hor) it is the direct radiation of the horizontal earth's surface of i wave bands, IoiFor i ripples
The solar constant of section, T(clear)iRepresent the Radiation Transmittance of i wave bands, ErFor solar distance correction coefficient, cos (Z) is sun day
Drift angle cosine.
7. consider the earth's surface solar radiation remote sensing computing system of cloud-orographic effect, it is characterised in that including：
Cloud Shadow recognition unit, for the recognizer using rolling topography cloud shadow region, goal pels are obtained intricately
Cloud shade under the conditions of shape, earth's surface solar radiation is distinguished by shadow region that cloud influences and the nonshaded area not influenceed by cloud
Domain；
Solar radiation unit, for estimating the earth's surface solar radiation of the non-hatched area using clear sky radiance computational methods, adopt
The earth's surface solar radiation of the shadow region is estimated with skies Radiation calculation schemes；
Topographical correction unit, for Combining with terrain bearing calibration, earth's surface solar radiation and described non-the moon to the shadow region
The earth's surface solar radiation in shadow zone domain is corrected.
8. system according to claim 7, it is characterised in that the solar radiation unit is specifically used for：
The Radiation Transmittance of first band and second band is obtained, calculates the earth's surface sun of the non-hatched area under the conditions of fine day
Directly radiate；
The earth's surface solar scattered radiation of the non-hatched area under the conditions of clear sky is divided into directly scattering radiation and due to ground vapour
Between multiple reflections and it is caused scattering radiation two parts be respectively calculated.
9. system according to claim 7, it is characterised in that the solar radiation unit also particularly useful for：
Using the relation between 0.63-0.69 mu m wavebands reflectivity and cloud optical thickness, it is distant that high-resolution is obtained by interpolation method
Feel the cloud optical thickness of image；
Calculate the Radiation Transmittance of first band and second band cloud respectively according to the cloud optical thickness；
The earth's surface solar radiation of the shadow region under the conditions of the skies is calculated according to the Radiation Transmittance of the cloud.
10. system according to claim 7, it is characterised in that the topographical correction unit is specifically used for：
According to the appraising model of domatic total solar radiation, respectively from direct radiation, sky radiation and the neighbour of domatic upper reception
The earth's surface of earth's surface solar radiation and the non-hatched area of three aspects of spurious radiation of near-earth shape to the shadow region is too
Sun radiation is corrected.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710928918.6A CN107767448A (en) | 2017-10-09 | 2017-10-09 | Consider the earth's surface solar radiation remote sensing computational methods and system of cloud orographic effect |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710928918.6A CN107767448A (en) | 2017-10-09 | 2017-10-09 | Consider the earth's surface solar radiation remote sensing computational methods and system of cloud orographic effect |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107767448A true CN107767448A (en) | 2018-03-06 |
Family
ID=61267844
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710928918.6A Pending CN107767448A (en) | 2017-10-09 | 2017-10-09 | Consider the earth's surface solar radiation remote sensing computational methods and system of cloud orographic effect |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107767448A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108596153A (en) * | 2018-05-10 | 2018-09-28 | 四川省冶地工程勘察设计有限公司 | Remote sensing image defends piece vertical control point extraction method and data processing method |
| CN109543721A (en) * | 2018-11-05 | 2019-03-29 | 中国科学院寒区旱区环境与工程研究所 | A kind of solar irradiance ultra-short term forecasting procedure under fine with occasional clouds weather condition |
| CN110186823A (en) * | 2019-06-26 | 2019-08-30 | 中国科学院遥感与数字地球研究所 | A kind of aerosol optical depth inversion method |
| CN111199092A (en) * | 2019-11-11 | 2020-05-26 | 南京深全人工智能技术研发有限公司 | Solar radiation remote sensing estimation method and system and data processing device |
| CN111428975A (en) * | 2020-03-13 | 2020-07-17 | 中国气象局公共气象服务中心 | Solar energy resource evaluation method and device and computer readable storage medium |
| CN111488669A (en) * | 2020-02-24 | 2020-08-04 | 宿迁学院 | Calculation method for determining optimal inclination angle of fixed solar device |
| CN113029332A (en) * | 2021-03-05 | 2021-06-25 | 阳光电源股份有限公司 | Satellite cloud picture prediction method, device and processor |
| CN116188633A (en) * | 2023-04-25 | 2023-05-30 | 之江实验室 | Method, device, medium and electronic equipment for generating simulated remote sensing image |
| CN116722544A (en) * | 2023-08-02 | 2023-09-08 | 北京弘象科技有限公司 | Distributed photovoltaic short-term prediction method and device, electronic equipment and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103198314A (en) * | 2013-02-20 | 2013-07-10 | 北京农业信息技术研究中心 | Remote sensing image radiation correction method |
| CN104406686A (en) * | 2014-12-10 | 2015-03-11 | 西北师范大学 | Estimation method of incident radiation of solar shortwaves under complex terrain conditions |
| CN105242247A (en) * | 2015-09-08 | 2016-01-13 | 首都师范大学 | Improved topographic radiation correction method |
-
2017
- 2017-10-09 CN CN201710928918.6A patent/CN107767448A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103198314A (en) * | 2013-02-20 | 2013-07-10 | 北京农业信息技术研究中心 | Remote sensing image radiation correction method |
| CN104406686A (en) * | 2014-12-10 | 2015-03-11 | 西北师范大学 | Estimation method of incident radiation of solar shortwaves under complex terrain conditions |
| CN105242247A (en) * | 2015-09-08 | 2016-01-13 | 首都师范大学 | Improved topographic radiation correction method |
Non-Patent Citations (1)
| Title |
|---|
| 李彬: "起伏地形上云阴影畸变对地表太阳辐射影响研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108596153A (en) * | 2018-05-10 | 2018-09-28 | 四川省冶地工程勘察设计有限公司 | Remote sensing image defends piece vertical control point extraction method and data processing method |
| CN109543721A (en) * | 2018-11-05 | 2019-03-29 | 中国科学院寒区旱区环境与工程研究所 | A kind of solar irradiance ultra-short term forecasting procedure under fine with occasional clouds weather condition |
| CN109543721B (en) * | 2018-11-05 | 2023-07-04 | 中国科学院西北生态环境资源研究院 | Solar irradiance ultra-short-term forecasting method under sunny cloudy weather condition |
| CN110186823A (en) * | 2019-06-26 | 2019-08-30 | 中国科学院遥感与数字地球研究所 | A kind of aerosol optical depth inversion method |
| CN111199092A (en) * | 2019-11-11 | 2020-05-26 | 南京深全人工智能技术研发有限公司 | Solar radiation remote sensing estimation method and system and data processing device |
| CN111488669A (en) * | 2020-02-24 | 2020-08-04 | 宿迁学院 | Calculation method for determining optimal inclination angle of fixed solar device |
| CN111488669B (en) * | 2020-02-24 | 2023-12-22 | 宿迁学院 | Calculation method for determining optimal inclination angle of fixed solar device |
| CN111428975B (en) * | 2020-03-13 | 2023-09-29 | 中国气象局公共气象服务中心 | Solar resource evaluation method and device and computer readable storage medium |
| CN111428975A (en) * | 2020-03-13 | 2020-07-17 | 中国气象局公共气象服务中心 | Solar energy resource evaluation method and device and computer readable storage medium |
| CN113029332A (en) * | 2021-03-05 | 2021-06-25 | 阳光电源股份有限公司 | Satellite cloud picture prediction method, device and processor |
| CN116188633A (en) * | 2023-04-25 | 2023-05-30 | 之江实验室 | Method, device, medium and electronic equipment for generating simulated remote sensing image |
| CN116722544A (en) * | 2023-08-02 | 2023-09-08 | 北京弘象科技有限公司 | Distributed photovoltaic short-term prediction method and device, electronic equipment and storage medium |
| CN116722544B (en) * | 2023-08-02 | 2023-10-20 | 北京弘象科技有限公司 | Distributed photovoltaic short-term prediction method and device, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107767448A (en) | Consider the earth's surface solar radiation remote sensing computational methods and system of cloud orographic effect | |
| Hofierka et al. | Physically-based land surface temperature modeling in urban areas using a 3-D city model and multispectral satellite data | |
| Mueller et al. | The CM-SAF operational scheme for the satellite based retrieval of solar surface irradiance—A LUT based eigenvector hybrid approach | |
| Yan et al. | Topographic radiation modeling and spatial scaling of clear-sky land surface longwave radiation over rugged terrain | |
| Journée et al. | Improving the spatio-temporal distribution of surface solar radiation data by merging ground and satellite measurements | |
| Long et al. | A two-source trapezoid model for evapotranspiration (TTME) from satellite imagery | |
| Bode et al. | Subcanopy Solar Radiation model: Predicting solar radiation across a heavily vegetated landscape using LiDAR and GIS solar radiation models | |
| Dubayah et al. | Modeling topographic solar radiation using GOES data | |
| Zhang et al. | A decadal (2008–2017) daily evapotranspiration data set of 1 km spatial resolution and spatial completeness across the North China Plain using TSEB and data fusion | |
| Ma et al. | Determination of regional net radiation and soil heat flux over a heterogeneous landscape of the Tibetan Plateau | |
| CN109460532B (en) | Solar direct radiation remote sensing calculation method and device | |
| Leon et al. | Large‐scale advection of continental aerosols during INDOEX | |
| CN102338869B (en) | Inversion method and system of downlink shortwave radiation and photosynthetically active radiation data | |
| CN104535979A (en) | Remote sensing inversion method and system for land cloud optical thickness | |
| Yan et al. | Modeling surface longwave radiation over high-relief terrain | |
| CN107145848B (en) | A kind of wind erosion of soil monitoring method and system based on remotely-sensed data | |
| Isabelle et al. | Solar radiation transmittance of a boreal balsam fir canopy: Spatiotemporal variability and impacts on growing season hydrology | |
| Yoon et al. | Evaluation of hourly solar radiation on inclined surfaces at Seoul by Photographical Method | |
| CN113744249B (en) | Marine ecological environment damage investigation method | |
| CN114152350B (en) | Earth surface temperature inversion method considering urban three-dimensional geometrical structure influence | |
| Zhang et al. | Simplified SEBAL method for estimating vast areal evapotranspiration with MODIS data | |
| CN104880701A (en) | Satellite-borne sensor imaging simulation method and device | |
| Xian et al. | A uniform model for correcting shortwave downward radiation over rugged terrain at various scales | |
| Gastellu-Etchegorry et al. | Recent improvements in the dart model for atmosphere, topography, large landscape, chlorophyll fluorescence, satellite image inversion | |
| Kim et al. | Atmospheric correction of true-color RGB imagery with limb area-blending based on 6S and satellite image enhancement techniques using geo-kompsat-2A advanced meteorological imager data |
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: 20180306 |
|
| RJ01 | Rejection of invention patent application after publication |