CN109186774A - Surface temperature information acquisition method, device, computer equipment and storage medium - Google Patents
Surface temperature information acquisition method, device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN109186774A CN109186774A CN201811002215.1A CN201811002215A CN109186774A CN 109186774 A CN109186774 A CN 109186774A CN 201811002215 A CN201811002215 A CN 201811002215A CN 109186774 A CN109186774 A CN 109186774A
- Authority
- CN
- China
- Prior art keywords
- surface temperature
- temperature information
- day
- clear sky
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000004927 fusion Effects 0.000 claims abstract description 115
- 238000004590 computer program Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000007500 overflow downdraw method Methods 0.000 abstract description 2
- 238000002310 reflectometry Methods 0.000 description 24
- 230000002123 temporal effect Effects 0.000 description 18
- 238000004422 calculation algorithm Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 239000002689 soil Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 6
- 230000005855 radiation Effects 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 6
- 241000208340 Araliaceae Species 0.000 description 5
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 5
- 235000003140 Panax quinquefolius Nutrition 0.000 description 5
- 235000008434 ginseng Nutrition 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000003973 irrigation Methods 0.000 description 2
- 230000002262 irrigation Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000002844 melting Methods 0.000 description 2
- 230000008018 melting Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000000179 transient infrared spectroscopy Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
This application involves a kind of surface temperature information fusion method, device, computer equipment and storage mediums.The described method includes: obtaining the multi-source surface temperature information of target area;The multi-source surface temperature information includes at least the first surface temperature information and the second surface temperature information;The first surface temperature information is the surface temperature data of each day of clear sky;The space scale of the first surface temperature information is less than the space scale of the second surface temperature information;According to the space scale of the first surface temperature information, temporal-spatial fusion is carried out to the first surface temperature information and the second surface temperature information, calculates the fusion surface temperature information of the target area.It can be improved the accuracy of fusion surface temperature information using this method.
Description
Technical field
This application involves satellite remote sensing technology fields, more particularly to a kind of surface temperature information acquisition method, device, meter
Calculate machine equipment and storage medium.
Background technique
Surface temperature information is to study water, energy circulation and the important state parameter of environmental change, is widely used in water
The research fields such as text, meteorological, climate change and environmental monitoring.Traditional surface temperature information observation is often based upon point scale and opens
Exhibition, cannot provide time, spatially continuous observation.And the progress of remote sensing technology, this limitation is breached, so that on a large scale
The surface temperature information for obtaining different time scales and space scale (abbreviation spatial and temporal scales) is possibly realized.
The remotely-sensed data of different satellite sensors has different spatial and temporal resolution characteristics.Generally, there is higher spatial
The remotely-sensed data of resolution ratio (i.e. lower space scale), temporal resolution are lower (i.e. time scale is larger);And have higher
The remotely-sensed data of temporal resolution (i.e. lower time scale), spatial resolution are lower (i.e. space scale is larger).In order to obtain
Taking both has high spatial resolution or has the surface temperature information of high time resolution feature, and traditional technology is passed through frequently with multi-source
Multiple dimensioned Remote Sensing Data Fusion Algorithm method obtains the surface temperature information of high-spatial and temporal resolution.
However, the surface temperature information obtained according to current multi-source surface temperature information fusion method, however it remains quasi-
The problem of true property difference.
Summary of the invention
Based on this, it is necessary in view of the above technical problems, provide a kind of surface temperature information acquisition method that accuracy is high,
Device, computer equipment and storage medium.
In a first aspect, a kind of surface temperature information acquisition method, which comprises
Obtain the multi-source surface temperature information of target area;The multi-source surface temperature information includes at least the first earth's surface temperature
Spend information and the second surface temperature information;The first surface temperature information is the surface temperature data of each day of clear sky;Described
The space scale of one surface temperature information is less than the space scale of the second surface temperature information;
According to the space scale of the first surface temperature information, to the first surface temperature information and second ground
Table temperature information carries out temporal-spatial fusion, calculates the fusion surface temperature information of the target area.
Second aspect, a kind of surface temperature information acquisition device, described device include:
Data obtaining module, for obtaining the multi-source surface temperature information of target area;The multi-source surface temperature information
Including at least the first surface temperature information and the second surface temperature information;The first surface temperature information is the ground of each day of clear sky
Table temperature data;The space scale of the first surface temperature information is less than the space scale of the second surface temperature information;
Information Fusion Module, for the space scale according to the first surface temperature information, to the first earth's surface temperature
It spends information and the second surface temperature information carries out temporal-spatial fusion, calculate the fusion surface temperature information of the target area.
The third aspect, a kind of computer equipment, including memory and processor, the memory are stored with computer journey
Sequence, the processor perform the steps of when executing the computer program
Obtain the multi-source surface temperature information of target area;The multi-source surface temperature information includes at least the first earth's surface temperature
Spend information and the second surface temperature information;The first surface temperature information is the surface temperature data of each day of clear sky;Described
The space scale of one surface temperature information is less than the space scale of the second surface temperature information;
According to the space scale of the first surface temperature information, to the first surface temperature information and second ground
Table temperature information carries out temporal-spatial fusion, calculates the fusion surface temperature information of the target area.
Fourth aspect, a kind of computer readable storage medium are stored thereon with computer program, the computer program quilt
Processor performs the steps of when executing
Obtain the multi-source surface temperature information of target area;The multi-source surface temperature information includes at least the first earth's surface temperature
Spend information and the second surface temperature information;The first surface temperature information is the surface temperature data of each day of clear sky;Described
The space scale of one surface temperature information is less than the space scale of the second surface temperature information;
According to the space scale of the first surface temperature information, to the first surface temperature information and second ground
Table temperature information carries out temporal-spatial fusion, calculates the fusion surface temperature information of the target area.
Above-mentioned surface temperature information acquisition method, device, computer equipment and storage medium, by obtaining target area
First surface temperature information and the second surface temperature information, the first surface temperature information are the surface temperature data of each day of clear sky;
Then according to the space scale of the first surface temperature information, when being carried out to the first surface temperature information and the second surface temperature information
Sky fusion, calculates the fusion surface temperature information of target area;Because the space scale of the first surface temperature information is less than second
The space scale of surface temperature information, therefore the first surface temperature information spatial resolution with higher, the second surface temperature
Information temporal resolution with higher, fusion surface temperature information have the spatial resolution of the first surface temperature information simultaneously
With the temporal resolution of the second surface temperature information, i.e., with high time resolution while have high spatial resolution;It is heavier
It wants, because the first surface temperature information used corresponds to each day of clear sky, avoids because of invalid data caused by cloud pollution etc.,
Avoid the space point that ensure that fusion surface temperature information because of the extra error that the method for reconstructing such as space interpolation introduce
Resolution improves the accuracy of fusion surface temperature information.
Detailed description of the invention
Fig. 1 is the application environment schematic diagram of surface temperature information acquisition method in one embodiment;
Fig. 2 is the flow diagram of surface temperature information acquisition method in one embodiment;
Fig. 3 is the flow diagram of surface temperature information acquisition method in another embodiment;
Fig. 4 is the flow diagram that fusion surface temperature information calculates step in one embodiment;
Fig. 5 is the schematic diagram of surface temperature information acquisition method in one embodiment;
Fig. 6 is the schematic diagram one of the result obtained in one embodiment using surface temperature information acquisition method;
Fig. 7 is the schematic diagram two of the result obtained in one embodiment using surface temperature information acquisition method;
Fig. 8 is the schematic diagram three of the result obtained in one embodiment using surface temperature information acquisition method;
Fig. 9 is the schematic diagram four of the result obtained in one embodiment using surface temperature information acquisition method;
Figure 10 is the structural block diagram of surface temperature information acquisition device in one embodiment;
Figure 11 is the structural block diagram of surface temperature information acquisition device in another embodiment;
Figure 12 is the structural block diagram of surface temperature information acquisition device in further embodiment;
Figure 13 is the internal structure chart of computer equipment in one embodiment.
Specific embodiment
It is with reference to the accompanying drawings and embodiments, right in order to which the objects, technical solutions and advantages of the application are more clearly understood
The application is further elaborated.It should be appreciated that specific embodiment described herein is only used to explain the application, not
For limiting the application.
Surface temperature information acquisition method provided by the embodiments of the present application, can be applied to all kinds of terminal devices, server
Deng;Wherein, terminal device can be, but not limited to be various mainframe computers, personal computer, laptop, smart phone,
Tablet computer and portable wearable device.Illustratively, by taking terminal device as an example, surface temperature provided by the embodiments of the present application
Information acquisition method can be applied in information processing system as shown in Figure 1.Terminal device 11 can be with satellite server 12
It is interacted by network, the multi-source earth's surface for obtaining the satellite remote sensing date of multiple satellites or being obtained by multiple satellite remote sensing dates
The fusion surface temperature information obtained after temporal-spatial fusion multi-source surface temperature information can also be returned to satellite clothes by temperature information etc.
Business device 12, wherein store satellite remote sensing date or multi-source surface temperature information of multiple satellites etc. in satellite server 12;It defends
Star server 12 can be connect with satellite receiver 13, which can receive the satellite remote sensing that multiple satellites obtain
Data are simultaneously transferred to satellite server 12, wherein multiple satellites may include the first satellite 14 and the second satellite 15, the first satellite
14 can be Landsat (earth resources observation satellite), and the second satellite 15 can be TERRA (soil)/AQUA (water) satellite.It can
With understanding, terminal device 11 can also directly be connect with satellite receiver 13, obtain the satellite remote sensing that multiple satellites obtain
Data are participated in without satellite server 12.
In one embodiment, as shown in Fig. 2, providing a kind of surface temperature information acquisition method, the present embodiment is related to
Be terminal device using temporal-spatial fusion algorithm, believed according to the fusion surface temperature of multi-source surface temperature information acquisition target area
Breath, to realize the process for the surface temperature information for obtaining high spatial resolution high time resolution.This method can specifically include
Following steps:
S202 obtains the multi-source surface temperature information of target area.
Wherein, multi-source surface temperature information includes at least the first surface temperature information and the second surface temperature information;First
Surface temperature information is the surface temperature data of each day of clear sky;The space scale of first surface temperature information is less than the second earth's surface temperature
Spend the space scale of information.The multi-source surface temperature information can derive from the remotely-sensed data of multiple satellites, can pass through remote sensing
Surface temperature information is calculated in data.For example, the first surface temperature information derives from the remotely-sensed data of the first satellite;Second ground
Table temperature information derives from the remotely-sensed data of the second satellite.When the first surface temperature information derives from the first satellite, the is obtained
The surface temperature data that one satellite is obtained when target area is the day of clear sky are the first surface temperature information.
Wherein, space scale be the size of minimum unit that can be distinguished in detail on the corresponding remote sensing images of remotely-sensed data or
Size belongs to the characteristic of different satellite sensors.For example, the satellite of Landsat (earth resources observation satellite) series of satellites passes
The space scale of sensor is 30-100m, but its revisiting period (i.e. time scale) is about 16 days;MODIS (is mounted in earth observation
Sensor on system TERRA (soil)/AQUA (water) satellite) space scale be 250-1000m, but its revisiting period be 1
It.
It should be noted that the first satellite, at access target region, because of the factors such as cloud pollution, there are partial target areas
Domain can not obtain the situation of effective remotely-sensed data;And when the first satellite access target area, if target area is fine,
First satellite can completely obtain the remotely-sensed data of target area, then the data obtained at this time are referred to as the remotely-sensed data of the day of clear sky;
It is the first surface temperature information, the first surface temperature information by the surface temperature data that the remotely-sensed data of the day of clear sky obtains
It is made of the surface temperature data of each day of clear sky.It is understood that determining that the foundation of the first surface temperature information does not limit
Whether weather is sunny, but can the remotely-sensed data completely to obtain target area as foundation.
Illustratively, if target area is the city A, the first satellite is 10 to the access cycle in the city A;For March,
First satellite is March 1 to the city A access time first time, and the back-call time is March 11, and third time access time is 3
The moon 21, the 4th access time are March 31;If March 1, March 21, March 31 can daily obtain the complete distant of the city A
Feel data, then for March, the first surface temperature information of target area by March 1 the first satellite source target
The surface temperature data in region, the surface temperature data of the target area of first satellite source on March 21, March 31
The surface temperature data of the target area of one satellite source form.Wherein, the remote sensing that the access cycle of satellite and the satellite obtain
The temporal resolution of data is corresponding;Similarly, the earth's surface temperature that the access cycle of satellite obtains with the remotely-sensed data by the satellite
The temporal resolution for spending information is corresponding.
It is understood that if obtain surface temperature data cloud pollution etc. due to cause partial target region without
Method obtains valid data, then needs to rebuild the surface temperature data by the methods of space interpolation, obtain the part mesh
Mark the surface temperature data in region;Accordingly, because have passed through space reconstruction, extra error, therefore the earth's surface obtained are introduced
The accuracy relative reduction of temperature data.Therefore, it can guarantee surface temperature information using the surface temperature data of the day of clear sky
Spatial resolution improves the accuracy of surface temperature information.
In addition, for a certain space scale, it can merge that multiple space scales are identical or space scale is less than normal defends
The day of clear sky surface temperature data of star are the first surface temperature information, the day of clear sky surface temperature of the satellite less than normal to space scale
Data can be converted into the day of clear sky surface temperature data of the space scale by way of room for promotion scale;Illustratively,
The mode of room for promotion scale may include the methods of obscuring, being averaging.
Optionally, terminal device extracts target area in the earth's surface temperature of the day of clear sky from the surface temperature data of the first satellite
Degree is according to as the first surface temperature information.
S204, according to the space scale of the first surface temperature information, to the first surface temperature information and the second surface temperature
Information carries out temporal-spatial fusion, calculates the fusion surface temperature information of target area.
Wherein, the space scale of the first surface temperature information is the amount of space size used by the first surface temperature information
Degree, the space scale of remotely-sensed data generally corresponding with the first surface temperature information are identical;Similarly, corresponding with the remotely-sensed data
Satellite spatial resolution it is corresponding.
To the surface temperature information of two different time and space scales, even more than different time and space scales surface temperature information into
The algorithm of row temporal-spatial fusion may include: adaptive reflectivity Fusion Model (STARFM, Spatial and Temporal
Adaptive Reflectance Fusion Model), space time adaptive algorithm (STAARCH, Spatial
Temporal Adaptive Algorithm for Mapping Reflectance Change) and enhanced room and time
Adaptive reflectivity Fusion Model (ESTARFM, Enhanced Spatial and Temporal Adaptive
Reflectance Fusion Model) etc..
It is understood that generally, the space scale of the first surface temperature information is less than the second surface temperature information
Space scale, then the spatial resolution of the first surface temperature information compared to the second surface temperature information spatial resolution more
It is high;Correspondingly, the temporal resolution of the second surface temperature information is higher compared to the temporal resolution of the first surface temperature information.
In this way, empty blending algorithm when employed, according to the space scale of the first surface temperature information, to the first surface temperature information and the
Two surface temperature information carry out temporal-spatial fusion, or for according to the space scale and the first earth's surface temperature of the first surface temperature information
Information is spent, to the second surface temperature information NO emissions reduction;The sky of the fusion surface temperature information for the target area being finally calculated
Between scale be the first surface temperature information space scale, i.e. the spatial resolution of the fusion surface temperature information of target area is
The spatial resolution of first surface temperature information.
In above-mentioned surface temperature information acquisition method, by the first surface temperature information and the second ground that obtain target area
Table temperature information, the first surface temperature information are the surface temperature data of each day of clear sky;Then according to the first surface temperature information
Space scale, temporal-spatial fusion is carried out to the first surface temperature information and the second surface temperature information, calculates melting for target area
Close surface temperature information;Because the space scale of the first surface temperature information is less than the space scale of the second surface temperature information,
Therefore the first surface temperature information spatial resolution with higher, the second surface temperature information time resolution with higher
Rate, fusion surface temperature information simultaneously have the first surface temperature information spatial resolution and the second surface temperature information when
Between resolution ratio, i.e., with high time resolution while have high spatial resolution;It is prior, because of the first earth's surface used
Temperature information corresponds to each day of clear sky, avoids because of invalid data caused by cloud pollution etc., avoids because of the weight such as space interpolation
The extra error that construction method introduces, therefore ensure that the spatial resolution of fusion surface temperature information, improve fusion earth's surface temperature
Spend the accuracy of information.
In one embodiment, referring to shown in Fig. 3, another surface temperature information acquisition method, the present embodiment are provided
Pass through what is involved is terminal device and determine that the object time corresponding first refers to the day of clear sky with reference to the day of clear sky and second, using space-time
Blending algorithm, according to multi-source surface temperature information acquisition target area in object time corresponding fusion surface temperature information, with
Realize the process for obtaining the surface temperature information of high spatial resolution high time resolution.This method can specifically include following step
It is rapid:
S302 obtains the multi-source surface temperature information of target area.
Wherein, multi-source surface temperature information includes at least the first surface temperature information and the second surface temperature information;First
Surface temperature information is the surface temperature data of each day of clear sky;The space scale of first surface temperature information is less than the second earth's surface temperature
Spend the space scale of information.
The first surface temperature information or the second surface temperature information can directly acquire correspondence in one of the embodiments,
The data product of satellite obtains.For example, the second surface temperature information can be the land MODIS table product MOD11A1 (surface temperature)
Data, space scale 1000m.Meanwhile for the land MODIS table product, the data (earth's surface of MOD09Q1 can also be obtained
Reflectivity), the data (NDVI vegetation index) of MOD13A2.
In wherein another embodiment, the first surface temperature information or the second surface temperature information can be defended by correspondence
The remotely-sensed data of star is calculated;Remotely-sensed data includes the Reflectivity for Growing Season data of multiple wave bands.Illustratively, according to remote sensing number
It, can be specific as follows with atmospheric correction method according to the concrete processing procedure that surface temperature data are calculated:
1, remotely-sensed data is pre-processed, including radiation calibration, atmospheric correction etc., respectively reduce Satellite observation error,
Atmospheric effect bring error etc.;
2, according to pretreated remotely-sensed data, normalized site attenuation (NDVI, Normalized are calculated
Difference Vegetation Index), shown in calculation formula such as following formula (1):
NDVI=(NIR-R)/(NIR+R) (1)
Wherein, NIR represents the Reflectivity for Growing Season value near infrared band, and R represents the Reflectivity for Growing Season value at red wave band;
3, according to NDVI, vegetation coverage P is calculatedv, shown in calculation formula such as following formula (2):
Pv=[(NDVI-NDVISoil)/(NDVIVeg-NDVISoil)] (2)
Wherein, NDVISoilFor complete exposed soil or without the NDVI value in vegetative coverage region, NDVIVegTo be covered completely by vegetation
The NDVI value of the pixel of lid, generally takes empirical value NDVIVeg=0.70 and NDVISoil=0.05;In addition, when some pixel
When NDVI is greater than 0.70, PvIt can be with value for 1;When NDVI can be with value for 0 less than 0.05, Pv;
4, according to vegetation coverage Pv, radiance value of the black matrix in Thermal infrared bands when calculating surface temperature data are Ts
B (Ts), shown in calculation formula such as following formula (3):
B (Ts)=[Lλ-L↑-τ(1-ε)L↓]/(τε) (3)
Wherein, L ↑ be the upward radiance of atmosphere, L ↓ be downward terrestrial radiantion, brightness value, τ atmosphere is in Thermal infrared bands
Transmitance, L ↑, L ↓, τ be atmospheric profile parameter, public data acquisition can be inquired, such as data source can be http: //
atmcorr.gsfc.nasa.gov;LλFor the thermal infrared radiation brightness value that satellite sensor receives, pass through satellite sensor
Remotely-sensed data obtains, and ε is Land surface emissivity, is calculated by ε=0.004Pv+0.986;
5, according to B (Ts), surface temperature data Ts is calculated, shown in calculation formula such as following formula (4):
Ts=K2/ln(K1/B(Ts)+1) (4)
Wherein, K1And K2It is constant, for the remotely-sensed data of the TM satellite sensor of Landsat, K1=607.76W/
(m2μm sr), K2=1260.56K;For the remotely-sensed data of the ETM+ satellite sensor of Landsat, K1=666.09W/
(m2μm sr), K2=1282.71K;For the remotely-sensed data of the TIRS Band10 satellite sensor of Landsat, K1=
774.89W/(m2μm sr), K2=1321.08K, sr are surface of sphere.
It should be noted that above-mentioned calculating process is for a pixel, the size of pixel and above-mentioned remotely-sensed data
Space scale it is corresponding, and the surface temperature data of target area include the surface temperature number of each pixel of the target area
According to.
In wherein another embodiment, the vegetation of the surface temperature information and the first scale that can establish the first scale refers to
Corresponding relationship between number;According to the remotely-sensed data of corresponding relationship and the second scale, the surface temperature of the second scale is calculated
Information.Wherein, vegetation index can be obtained directly from satellite data product, can also be by various in satellite data product
Reflectivity for Growing Season data obtain;Vegetation index can be NDVI, can be soil adjusting vegetation index, difference environment vegetation refers to
Number, ratio vegetation index etc., such as ratio vegetation index RVI=Near/Nred;Near and Nred be respectively near infrared band and
The reflectivity of red spectral band.For example, 1 product of Landsat Level can be pre-processed (radiation calibration, atmospheric correction),
Reflectivity for Growing Season data and Thermal Infrared Data are obtained, normalized site attenuation (NDVI), surface albedo are calculated
(Albedo) and surface temperature data (LST).
Illustratively, the first surface temperature information can be by 1 data product of Landsat Level (after radiant correction
Remotely-sensed data) it is calculated, space scale 30m, then the Landsat of available each pixel corresponding to target area
1 data product of Level, the region of pixel are 30m × 30m.However, the surface temperature number of 1 data product of Landsat Level
Space scale according to (LST, Land Surface Temperature) is 60m × 60m.In practical applications, Landsat
The space scale of the Reflectivity for Growing Season value of blue, green, red, near-infrared and short-wave infrared is 30m × 30m in 1 data product of Level,
Simultaneously as the relationship between surface temperature-vegetation index is almost the same on each scale, therefore 60m × 60m space scale
LST and NDVI between relationship and 30m × 30m space scale LST and NDVI between relationship it is almost the same, therefore can be with
NO emissions reduction processing is carried out to the LST of 60m × 60m space scale of acquisition, obtains the LST of 30m × 30m space scale, it is specific to locate
Reason process can be as follows:
1, the NDVI of 30m × 30m space scale is calculated according to above formula (1), and calculates the space 30m × 30m according to formula (5)
The α of scale;
α=0.356r1+0.13r2+0.373r3+0.085r4+0.072r5-0.0018 (5)
Wherein, r1-r5Be respectively blue, red, near-infrared and short-wave infrared (1.55-1.75 and 2.09-2.35 μm) earth's surface it is anti-
Penetrate rate;
2, shown in such as formula of the relationship between LST and NDVI (6), according to LST, NDVI of 60m × 60m space scale and
α fitting, which returns, acquires fitting parameter a1~a15;
Wherein, α is surface albedo (Albedo), surface albedo refer to ground return amount of radiation and incident radiation amount it
Than characterization ground faces the absorption and albedo of solar radiation, is fitted by the Reflectivity for Growing Season of each wave band;60m × 60m's
NDVI and α can be obtained by NDVI the and α data of 30m × 30m space scale by fuzzy algorithmic approach respectively;
3, according to NDVI, α of 30m × 30m space scale and above-mentioned fitting parameter a1~a15, calculated and asked according to formula (6)
Obtain the LST of 30m × 30m space scale.
It should be noted that although above-mentioned example uses 1 data product of Landsat Level to be illustrated, originally
Application embodiment can be also used for the data product of other satellite sensors.In addition, generally also needing before temporal-spatial fusion to the
One surface temperature information or the second surface temperature information are coordinately transformed, and using unified coordinate system, are melted convenient for subsequent space-time
It closes.For example, Landsat initial data is that WGS1984 (sit by World Geodetic System world's the earth in 1984,1984
Mark system) UTM (Universal Transverse Mercator, Universal Transverse Mercator Projection coordinate system), and MODIS is former
Beginning data are sinusoidal projected coordinate system, the two coordinate system need to be carried out to unification, therefore can choose and be converted to MODIS data
With the consistent WGS1984UTM of Landsat coordinate system.
S304 determines that the object time corresponding first refers to the day of clear sky with reference to the day of clear sky and second from each day of clear sky.
Wherein, first is before the object time and with the object time apart from the nearest date with reference to the day of clear sky, second
It is after the object time and with the object time apart from the nearest date with reference to the day of clear sky.Object time can be the moment, when
It phase, can be corresponding with the time scale of the second surface temperature information.For example, if the time scale of the second surface temperature information is one
Week, then the object time can be certain week;If the time scale of the second surface temperature information is one day, the object time can be certain
It.Illustratively, if target is worked as on the exemplary basis that above-mentioned first satellite accesses the city A in the region that target area is the city A
When time is March 5, the first reference day of clear sky was March 1, and the second reference day of clear sky was March 21;It is March when the object time
At 22 days, the first reference day of clear sky was March 21, and the second reference day of clear sky was March 31.
S306 refers to the day of clear sky according in the space scale of the first surface temperature information, the first surface temperature information first
Corresponding surface temperature data and second are with reference in the day of clear sky corresponding surface temperature data, the second surface temperature information first
With reference to the day of clear sky corresponding surface temperature data and second with reference to the day of clear sky corresponding surface temperature data, the second surface temperature
Object time corresponding surface temperature data in information determine target area in object time corresponding fusion surface temperature letter
Breath.
In one embodiment, S306 can specifically include step shown in Fig. 4:
S402 refers to the day of clear sky according in the space scale of the first surface temperature information, the first surface temperature information first
Corresponding surface temperature data and second are with reference in the day of clear sky corresponding surface temperature data, the second surface temperature information first
With reference to the day of clear sky corresponding surface temperature data and second with reference to the day of clear sky corresponding surface temperature data, the second surface temperature
Object time corresponding surface temperature data in information determine that reference surface temperature of the target area in the first reference day of clear sky is believed
The reference surface temperature information that breath and target area refer to the day of clear sky second.
In the present embodiment, terminal device can space scale according to the first surface temperature information, the first surface temperature
First is corresponding with reference to the day of clear sky with reference in the day of clear sky corresponding surface temperature data, the second surface temperature information first in information
Object time corresponding surface temperature data, determine target area first in surface temperature data, the second surface temperature information
With reference to the reference surface temperature information of the day of clear sky;And the space scale according to the first surface temperature information, the first surface temperature
Second is corresponding with reference to the day of clear sky with reference in the day of clear sky corresponding surface temperature data, the second surface temperature information second in information
Object time corresponding surface temperature data, determine target area second in surface temperature data, the second surface temperature information
With reference to the reference surface temperature information of the day of clear sky.It is understood that fusion surface temperature of the target area in the object time is believed
Breath can refer to the day of clear sky second with reference to the reference surface temperature information of the day of clear sky and target area first by target area
Reference surface temperature information interpolation obtain.
In one of them embodiment, S402 can be accomplished in that
Target area is divided at least one target picture according to the space scale of the first surface temperature information by step (A)
Member.Specifically, when the space scale of the first surface temperature information is 30m × 30m, target area can be divided at least one
The subregion of a 30m × 30m, each sub-regions are respectively goal pels.
Step (B), according in the first surface temperature information first with reference to the day of clear sky corresponding surface temperature data, the second ground
First with reference to the day of clear sky corresponding surface temperature data and object time corresponding surface temperature data in table temperature information, really
The reference surface temperature information that the pixel that sets the goal refers to the day of clear sky first.
Step (C), according in the first surface temperature information second with reference to the day of clear sky corresponding surface temperature data, the second ground
Second with reference to the day of clear sky corresponding surface temperature data and object time corresponding surface temperature data in table temperature information, really
The reference surface temperature information that the pixel that sets the goal refers to the day of clear sky second.
Illustratively, reference surface temperature information of the determining goal pels in the first reference day of clear sky in step (B), or
Person, determine in step (C) goal pels second with reference to the day of clear sky reference surface temperature information, specifically can be using including pass
It is formula (7)
Or the deformation relationship formula (7.1) of the relational expression (7)
Be calculated goal pels first with reference to the day of clear sky reference surface temperature information, alternatively, goal pels are the
The reference surface temperature information of the two reference days of clear sky;Alternatively, can be directly usedThis formula obtains,
In, Fk(xw/2,yw/2,tp) it is exactly reference surface temperature of the goal pels in the first reference day of clear sky or the second reference day of clear sky
Information.
In above-mentioned formula, tpFor the object time, k=m or n, tmThe day of clear sky, t are referred to for firstnClear sky is referred to for second
Day, (xw/2,yw/2) be goal pels coordinate, w is the size of search window centered on goal pels, and N is search window
The number of the similar pixel of middle goal pels, WiFor the pixel weight of i-th in search window similar pixel, ViFor search window
In i-th of similar pixel conversion coefficient;L(xw/2,yw/2,tk) it is goal pels first reference in the first surface temperature information
The day of clear sky or second is with reference to the day of clear sky corresponding surface temperature data, M (xi,yi,tp) it is i-th of similar picture in search window
Member object time in the second surface temperature information corresponding surface temperature data, M (xi,yi,tk) it is i-th in search window
A similar pixel first reference day of clear sky or the second reference day of clear sky corresponding surface temperature in the second surface temperature information
Data.
It should be noted that in relational expression (7.1), c1For L (xw/2,yw/2,tk) weight coefficient, embody refer to clear sky
Day tkCorresponding first surface temperature information, to reference day of clear sky tkThe corresponding influence of factors with reference to surface temperature information,
c2ForWeight coefficient, that is, embody refer to day of clear sky tkIt is right respectively with the object time
The variation for the second surface temperature information answered, to reference day of clear sky tkThe corresponding influence factor with reference to surface temperature information is big
It is small.Weight coefficient c1And c2It can be obtained by regression analysis or machine learning scheduling algorithm, evaluation points can be fusion earth's surface temperature
Spend the error parameter between information and the surface temperature information of actual measurement.
Step (D) determines that target area exists according to goal pels in the reference surface temperature information of the first reference day of clear sky
The reference surface temperature information of the first reference day of clear sky;And according to goal pels second with reference to the day of clear sky reference earth's surface temperature
Information is spent, determines target area in the reference surface temperature information of the second reference day of clear sky.It is understood that target area exists
The reference surface temperature information of the first reference day of clear sky may include: each goal pels on the reference ground of the first reference day of clear sky
Table temperature information;Target area second with reference to the day of clear sky reference surface temperature information, may include: each goal pels the
The reference surface temperature information of the two reference days of clear sky.
S404 refers to the day of clear sky according in the space scale of the first surface temperature information, the second surface temperature information first
Corresponding surface temperature data and second are with reference to target in the day of clear sky corresponding surface temperature data, the second surface temperature information
Time corresponding surface temperature data determine target area in the first time weight and target area of the first reference day of clear sky
Second time weighting of the domain in the second reference day of clear sky.
In the present embodiment, terminal device can draw target area according to the space scale of the first surface temperature information
It is divided at least one goal pels;According in the second surface temperature information first with reference to the day of clear sky corresponding surface temperature data with
And second with reference to object time corresponding surface temperature number in the day of clear sky corresponding surface temperature data, the second surface temperature information
According to, determine goal pels first with reference to the day of clear sky first time weight and goal pels second with reference to the day of clear sky
Second time weighting.
In another embodiment, S404 can be accomplished in that
Step (a), according in the second surface temperature information first with reference to the day of clear sky corresponding surface temperature data and mesh
The second reference day of clear sky corresponding surface temperature data in mark time corresponding surface temperature data, the second surface temperature information,
Determine goal pels in the first time weight of the first reference day of clear sky.
Step (b), according in the second surface temperature information second with reference to the day of clear sky corresponding surface temperature data and mesh
The first reference day of clear sky corresponding surface temperature data in mark time corresponding surface temperature data, the second surface temperature information,
Determine goal pels in the second time weighting of the second reference day of clear sky.
Illustratively, goal pels are determined in step (a) in the first time weight of the first reference day of clear sky, alternatively, step
Suddenly in (b) determine goal pels second with reference to the day of clear sky the second time weighting, specifically can using comprising
Relational expression (8) or the relational expression (8) deformation relationship formula, goal pels are calculated first with reference to clear sky
The first time weight of day, alternatively, determining goal pels in the second time weighting of the second reference day of clear sky;Alternatively, can be straight
Connect useThis formula obtains, wherein TkIt is exactly target
Pixel is first with reference to the first time weight of the day of clear sky or the second time weighting of the second reference day of clear sky.
In above-mentioned formula, tpFor the object time, k=m or n, tmThe day of clear sky, t are referred to for firstnClear sky is referred to for second
Day, w is the size of the search window centered on goal pels;M(xj,yi,tp) it is that coordinate is (x in search windowj,yi) picture
Member object time corresponding surface temperature data, M (x in the second surface temperature informationj,yi,tk) it is that coordinate is in search window
(xj,yi) pixel in the second surface temperature information first with reference to the day of clear sky or second with reference to the day of clear sky corresponding earth's surface temperature
Degree evidence.
Target area is referred to the reference surface temperature information of the day of clear sky and multiplying for first time weight first by S406
Product, target area carry out summation behaviour with reference to the reference surface temperature information of the day of clear sky and the product of the second time weighting second
Make, obtains target area in object time corresponding fusion surface temperature information.
Specifically, above-mentioned steps 406 can be accomplished in that
Step I, reference surface temperature information and goal pels by goal pels in the first reference day of clear sky are joined first
Examine the product of the first time weight of the day of clear sky, the reference surface temperature information and target that goal pels refer to the day of clear sky second
Pixel carries out sum operation in the product of the second time weighting of the second reference day of clear sky, obtains goal pels in the object time pair
The fusion surface temperature information answered.
Illustratively, step I can be calculated with relational expression (9) or the deformation relationship formula of the relational expression (9):
F(xw/2,yw/2,tp)=Tm×Fm(xw/2,yw/2,tp)+Tn×Fn(xw/2,yw/2,tp) (9)
Wherein, F (xw/2,yw/2,tp) it is goal pels in object time tpCorresponding fusion surface temperature information, Fm(xw/2,
yw/2,tp) it is goal pels in object time tpThe reference surface temperature information of the corresponding first reference day of clear sky, TmFor target picture
First time weight of the member in the first reference day of clear sky, Fn(xw/2,yw/2,tp) it is goal pels in object time tpCorresponding
The reference surface temperature information of the two reference days of clear sky, TnThe second time weighting for referring to the day of clear sky second for goal pels.
Step II determines target area in mesh according to goal pels in object time corresponding fusion surface temperature information
Mark time corresponding fusion surface temperature information.It is understood that target area is in object time corresponding fusion earth's surface temperature
Degree information may include: to constitute the fusion surface temperature information of each goal pels of the target area.
In the present embodiment, similar pixel can be all pixels in search window in addition to goal pels, can also be
Similar pixel similar to the albedo of goal pels in search window.Because different albedos represents the reflection of different zones
Characteristic, while characterizing the vegetation characteristic of different regions, temperature characterisitic.Surface temperature information is more accurately merged in order to obtain, is shown
Example property, the albedo of all pixels in search window is obtained, respectively all pixels by search window in addition to goal pels
Albedo and the albedos of goal pels be compared;If the difference of the albedo of the albedo and goal pels of the pixel
Absolute value is less than preset threshold, it is determined that the pixel is the similar pixel of goal pels.The preset threshold can rule of thumb be set
It is fixed, or the percentage of the albedo for goal pels, such as 5%, 10% etc..Albedo is the reflectivity by different-waveband
Data are obtained by weighted sum or integral, are referred to above formula (5) and are calculated.
In above-mentioned surface temperature information acquisition method, by the first surface temperature information and the second ground that obtain target area
Table temperature information, the first surface temperature information are the surface temperature data of each day of clear sky;Then distance is determined from each day of clear sky
The object time rear and front end nearest first reference day of clear sky and second refers to the day of clear sky;According to the sky of the first surface temperature information
Between scale, first refer to the day of clear sky pair with reference to the day of clear sky corresponding surface temperature data and second in the first surface temperature information
First with reference to the day of clear sky corresponding surface temperature data and the second ginseng in the surface temperature data answered, the second surface temperature information
It examines the day of clear sky corresponding surface temperature data, object time corresponding surface temperature data in the second surface temperature information, determines
Target area is in object time corresponding fusion surface temperature information;Because the space scale of the first surface temperature information is less than
The space scale of two surface temperature information, therefore the first surface temperature information spatial resolution with higher, the second earth's surface temperature
Information temporal resolution with higher is spent, fusion surface temperature information has the spatial discrimination of the first surface temperature information simultaneously
The temporal resolution of rate and the second surface temperature information, i.e., with high time resolution while have high spatial resolution;More
Important, because the first surface temperature information used corresponds to each day of clear sky, avoid because of invalid number caused by cloud pollution etc.
According to avoiding because of the extra error that the method for reconstructing such as space interpolation introduce, ensure that the sky of fusion surface temperature information
Between resolution ratio, improve fusion surface temperature information accuracy;Furthermore in the present embodiment, before using the distance objective time
The surface temperature data of the both ends first nearest reference day of clear sky and the second reference day of clear sky afterwards, compared to using single-ended reference
For the surface temperature data of the day of clear sky, the accuracy for calculating object time corresponding fusion surface temperature information is improved.
Further, in the present embodiment, surface temperature information acquisition method can also include: to obtain the first earth's surface temperature
Spend corresponding first reflectivity data of information and corresponding second reflectivity data of the second surface temperature information;First reflection
Rate data and the second reflectivity data respectively include the reflectivity data of at least one wave band;For each wave band, according to first
The space scale of reflectivity data carries out temporal-spatial fusion to the first reflectivity data and the second reflectivity data, it is anti-to calculate fusion
Penetrate rate data.Specifically, the fusion reflectivity data of each pixel of target area can be calculated according to above formula (7)~(9).
Generally, the space scale of the first reflectivity data is equal to the space scale of the first surface temperature information.It is understood that
Above-mentioned fusion reflectivity data has higher temporal resolution and spatial resolution;In addition, anti-according to the fusion of each wave band
Rate data are penetrated, the albedo data of each pixel can be calculated;Therefore in the present embodiment, on the one hand according to more accurate albedo
Data can search the similar pixel of more accurate goal pels, another party in the search window centered on goal pels
Face can calculate more accurate pixel spectral weight according to more accurate reflectivity data, melt in short, can further increase
Close the accuracy of surface temperature information.
It should be noted that in search window i-th of similar pixel pixel weight Wi, pixel similar to this and target picture
Member distance relation it is related, can also pixel similar to this spectral weight it is related.Illustratively, WiIt can according to the following formula (10)
It calculates:
Wherein, the spectral weight R of similar pixeliPhysical significance be higher spatial resolution similar pixel and with the phase
Like the spectral correlation coefficient of the corresponding pixel compared with low spatial resolution of pixel;Specifically, it can be Landsat and MODIS number
According to the spectral correlation coefficient between pixel.
In one embodiment, in the space scale according to the first surface temperature information, to the first surface temperature information and
Second surface temperature information carries out temporal-spatial fusion, before the fusion surface temperature information for calculating target area, method further include: when
There are when invalid data in second surface temperature information, space reconstruction is carried out to the second surface temperature information.Optionally, according to
The space scale of one surface temperature information carries out temporal-spatial fusion, meter to the first surface temperature information and the second surface temperature information
The fusion surface temperature information for calculating target area, can specifically include: according to the space scale of the first surface temperature information, to the
The second surface temperature information after one surface temperature information and space reconstruction carries out temporal-spatial fusion, with calculating the fusion of target area
Table temperature information.
To the second surface temperature information for needing space reconstruction, neighbouring time observation and the earth's surface with complementary information are selected
Temperature data can be rebuild based on weighted linear regression algorithm as reference data, can also be based on other spatial interpolation algorithms
Rebuild, for example, gram in golden space reconstruction method, the classification of multidate surface temperature and return method for reconstructing, space-time closes on pixel
Method for reconstructing and the consideration modified multidate surface temperature space reconstruction method of energy balance etc..For example, it is directed to target area, if
The time scale of second surface temperature information is one day, and September 1st second surface temperature information is led because of situations such as cloud pollution
It causes that data lacks completely or data portion lacks, then can choose August 31st second surface temperature information or the of September 2nd
Two surface temperature information rebuild September missing data on the 1st as reference data.According to above-mentioned reference data to missing
The method that data carry out space reconstruction using weighted linear regression algorithm, specifically can be as shown in following formula (11)~(12):
Ts=aTs'+b
Di=| Ts,0'-Ts,i'+ξ|·|V0-Vi+ξ|·((xi-x0)2+(yi-y0)2)
|Ts,0-Ts,i|≤Tthd
|V0-Vi|≤Vthd (12)
Wherein, TsAnd Ts' the surface temperature data of missing data and reference data are respectively indicated, a and b are regression coefficient, a
Amount of gain between missing data and reference data, b are offset.In order to solve regression coefficient a and b, with current pixel it is
Heart windowing, available by weighted least square in window:Each pixel in the window in expression missing data
Mean value,Indicate the mean value of each pixel in the window in reference data.In addition, subscript " 0 " indicates current pixel, i.e., in window
Center pel;Subscript " i " indicates i-th of similar pixel in window, then WiFor the weighted value of i-th of similar pixel, DiIndicate the
Terrain object attribute (surface temperature data and vegetation index) variation between i similar pixels and current pixel and they and current picture
The space length relationship of member, (x0,y0) be current pixel coordinate, (xi,yi) be i-th of similar pixel coordinate, Ts,0' it is ginseng
Examine the surface temperature data of current pixel in data, Ts,0For the surface temperature data of pixel current in missing data, Ts,iIt is scarce
Lose the surface temperature data of i-th of similar pixel in data, Ts,i' for the surface temperature number of i-th of similar pixel in reference data
According to ViFor the vegetation index of i-th of similar pixel, V0For the vegetation index of current pixel, TthdAnd VthdTo respectively correspond earth's surface temperature
The similarity threshold of degree evidence and vegetation index can be calculated by Local standard deviation, and ξ indicates that one tends to 0 value (such as
0.001) to ensure D value not equal to 0, the description that vegetation index V is referred to previously mentioned vegetation index is calculated, here no longer
It repeats.It should be noted that window size can be made by oneself, illustratively, window can be M × M centered on current pixel
The window of size, then the pixel number in window in addition to current pixel is M at this time2- 1.Above-mentioned calculating process is to be directed to work as
For preceding pixel, current pixel is that the pixel of the second surface temperature information is lacked in target area, size and above-mentioned the
The space scale of two surface temperature information is corresponding.
It, can be by closing on day for missing data one day of the MODIS of 1km space scale for example, referring to shown in Fig. 5
The MODIS data of the 1km space scale with complementary information of phase are made as reference data, the NDVI data of clear sky MOD13A2
For vegetation index data, same day missing data is rebuild by weighted regression algorithm, and then obtains the earth's surface temperature of daily 1km × 1km
Degree evidence.Then according to the MODIS LST of the Landsat LST and 1km × 1km of 30m × 30m, pass through ESTARFM temporal-spatial fusion
Algorithm obtains the surface temperature data of daily 30m × 30m.It further, can be in application ESTARFM temporal-spatial fusion algorithm
The similar pixel of goal pels is determined by the Reflectivity for Growing Season data of MOD09Q1.
In the present embodiment, when in the second surface temperature information there are when invalid data, to the second surface temperature information into
Row space reconstruction ensure that the comprehensive of the second surface temperature information, after avoiding because leading to temporal-spatial fusion there are invalid data
Fusion surface temperature information accuracy reduce.
Optionally, the second surface temperature information is daily surface temperature data.That is, the second surface temperature information
Time scale be one day, then the time scale of the fusion surface temperature information after temporal-spatial fusion is also one day, ensure that and melt
Close the temporal resolution of surface temperature information.
Optionally, the space scale of the first surface temperature information is 30m-100m;Fused surface temperature information space
Scale be equal to the first surface temperature information space scale, the space scale can to the corresponding farmland of different ground products into
Row is differentiated, the surface temperature information in the farmland suitable for studying different ground products.
Referring to shown in Fig. 6,7,8,9, the present embodiment application surface temperature information acquisition method utilizes actual measurement in Hetao Irrigation District
Data is examined.Hetao Irrigation District pattern of farming is broken, and ground surface type is complicated, tests the technical side of the present embodiment in this area
Case has more demonstration effect and promotional value.Referring to shown in Fig. 6, each initial MODIS LST being classified as on the same day, space weight
MODIS LST after building and with Landsat is fused merges LST, it is seen that fused fusion LST compares MODIS LST
Surface temperature spatial texture information is finer, i.e., spatial resolution is higher.Referring to shown in Fig. 7-9, actual measurement surface temperature is
The multipoint observation data of Different Crop vegetation in 2015, are obtained by way of combining with automonitor;With observation pair
Than root-mean-square error RMSE (Root Mean Square Error) reaches 2.31 DEG C, and deviation Bias reaches 6.66%;With the same period
1km × 1km surface temperature is compared, and the fusion surface temperature information of 30m × 30m is more fine and precision is improved, root mean square
Error RMSE and deviation Bias decrease.
It should be understood that although each step in the flow chart of Fig. 2-4 is successively shown according to the instruction of arrow,
These steps are not that the inevitable sequence according to arrow instruction successively executes.Unless expressly stating otherwise herein, these steps
Execution there is no stringent sequences to limit, these steps can execute in other order.Moreover, at least one in Fig. 2-4
Part steps may include that perhaps these sub-steps of multiple stages or stage are not necessarily in synchronization to multiple sub-steps
Completion is executed, but can be executed at different times, the execution sequence in these sub-steps or stage is also not necessarily successively
It carries out, but can be at least part of the sub-step or stage of other steps or other steps in turn or alternately
It executes.
In one embodiment, as shown in Figure 10, a kind of surface temperature information acquisition device 20 is provided, may include:
Data obtaining module 201 and information Fusion Module 202, in which:
Data obtaining module 201, for obtaining the multi-source surface temperature information of target area;Multi-source surface temperature information is extremely
It less include the first surface temperature information and the second surface temperature information;First surface temperature information is the surface temperature of each day of clear sky
Data;Space scale of the space scale of first surface temperature information less than the second surface temperature information;
Information Fusion Module 202, for the space scale according to the first surface temperature information, to the first surface temperature information
Temporal-spatial fusion is carried out with the second surface temperature information, calculates the fusion surface temperature information of target area.
In one of the embodiments, as shown in figure 11, in above-mentioned surface temperature information acquisition device 20 shown in Fig. 10
On the basis of, information Fusion Module 202 may include referring to day of clear sky determination unit 2021 and first information integrated unit 2022,
Wherein:
With reference to day of clear sky determination unit 2021, for determining that the object time corresponding first refers to clear sky from each day of clear sky
Day and second refers to the day of clear sky;Wherein, first is before the object time and with the object time apart from recently with reference to the day of clear sky
Date, second with reference to the day of clear sky is after the object time and with the object time apart from the nearest date;
First information integrated unit 2022, for the space scale according to the first surface temperature information, the first surface temperature
In information first with reference to the day of clear sky corresponding surface temperature data and second with reference to the day of clear sky corresponding surface temperature data, the
First with reference to the day of clear sky corresponding surface temperature data and the second reference day of clear sky corresponding earth's surface in two surface temperature information
Object time corresponding surface temperature data, determine target area in the object time in temperature data, the second surface temperature information
Corresponding fusion surface temperature information.
Optionally, first information integrated unit 2022 can be used for the space scale according to the first surface temperature information,
First with reference to the day of clear sky corresponding surface temperature data and the second reference day of clear sky corresponding earth's surface in one surface temperature information
First refers to the day of clear sky with reference to the day of clear sky corresponding surface temperature data and second in temperature data, the second surface temperature information
Object time corresponding surface temperature data, determine target area in corresponding surface temperature data, the second surface temperature information
The reference surface temperature that the reference surface temperature information and target area for referring to the day of clear sky first refer to the day of clear sky second
Information;Correspondingly according to the first reference day of clear sky in the space scale of the first surface temperature information, the second surface temperature information
The object time is corresponding in table temperature data and the second reference day of clear sky corresponding surface temperature data, the second surface temperature information
Surface temperature data, determine target area first with reference to the day of clear sky first time weight and target area second
With reference to the second time weighting of the day of clear sky;By target area when the reference surface temperature information of the first reference day of clear sky is with first
Between the product of weight, target area second with reference to the day of clear sky reference surface temperature information and the second time weighting product into
Row sum operation obtains target area in object time corresponding fusion surface temperature information.
Further, first information integrated unit 2022 can be used for the space scale according to the first surface temperature information,
Target area is divided at least one goal pels;Correspondingly according to the first reference day of clear sky in the first surface temperature information
First is corresponding with reference to the day of clear sky corresponding surface temperature data and object time in table temperature data, the second surface temperature information
Surface temperature data, determine goal pels first with reference to the day of clear sky reference surface temperature information;According to the first earth's surface temperature
Second is corresponding with reference to the day of clear sky with reference in the day of clear sky corresponding surface temperature data, the second surface temperature information second in degree information
Surface temperature data and object time corresponding surface temperature data, determine goal pels second with reference to the day of clear sky ginseng
Examine surface temperature information;According to goal pels in the reference surface temperature information of the first reference day of clear sky, determine that target area exists
The reference surface temperature information of the first reference day of clear sky;And according to goal pels second with reference to the day of clear sky reference earth's surface temperature
Information is spent, determines target area in the reference surface temperature information of the second reference day of clear sky.
Further, first information integrated unit 2022 can be used for according to the first reference in the second surface temperature information
The day of clear sky corresponding surface temperature data and object time corresponding surface temperature data, the second surface temperature information in second
With reference to the day of clear sky corresponding surface temperature data, determine goal pels in the first time weight of the first reference day of clear sky;According to
Second with reference to the day of clear sky corresponding surface temperature data and object time corresponding surface temperature in second surface temperature information
First with reference to the day of clear sky corresponding surface temperature data in data, the second surface temperature information, determine goal pels in the second ginseng
Examine the second time weighting of the day of clear sky.
Specifically, first information integrated unit 2022 can be used for goal pels on the reference ground of the first reference day of clear sky
Table temperature information and goal pels are fine in the second reference with reference to the product of the first time weight of the day of clear sky, goal pels first
The reference surface temperature information and goal pels of empty day is summed in the product of the second time weighting of the second reference day of clear sky
Operation, obtains goal pels in object time corresponding fusion surface temperature information;It is corresponding in the object time according to goal pels
Fusion surface temperature information, determine target area in object time corresponding fusion surface temperature information.
Specifically, first information integrated unit 2022 can be used for using comprising
Relational expression, be calculated goal pels first with reference to the day of clear sky reference surface temperature information, alternatively, target
Reference surface temperature information of the pixel in the second reference day of clear sky;Wherein, tpFor the object time, k=m or n, tmFor the first ginseng
It examines the day of clear sky, tnThe day of clear sky, (x are referred to for secondw/2,yw/2) be goal pels coordinate, w be searching centered on goal pels
The size of rope window, N are the number of the similar pixel of goal pels in search window, WiFor i-th in search window similar picture
The pixel weight of member, ViFor the conversion coefficient of i-th in search window similar pixel;L(xw/2,yw/2,tk) it is goal pels the
First with reference to the day of clear sky or the second reference day of clear sky corresponding surface temperature data, M (x in one surface temperature informationi,yi,tp)
For i-th in search window similar pixel in the second surface temperature information object time corresponding surface temperature data, M
(xi,yi,tk) it is i-th of similar pixel, the first reference day of clear sky or second in the second surface temperature information in search window
With reference to the day of clear sky corresponding surface temperature data.
Specifically, first information integrated unit 2022 can be used for using comprising
Relational expression, be calculated goal pels first with reference to the day of clear sky first time weight, alternatively, determine target
Second time weighting of the pixel in the second reference day of clear sky;Wherein, tpFor the object time, k=m or n, tmFor first with reference to fine
Empty day, tnThe day of clear sky is referred to for second, w is the size of the search window centered on goal pels;M(xj,yi,tp) it is search
Coordinate is (x in windowj,yi) pixel in the second surface temperature information object time corresponding surface temperature data, M (xj,
yi,tk) it is that coordinate is (x in search windowj,yi) pixel in the second surface temperature information first with reference to the day of clear sky or the
Two with reference to the day of clear sky corresponding surface temperature data.
In one of the embodiments, as shown in figure 12, in above-mentioned surface temperature information acquisition device 200 shown in Fig. 10
On the basis of, can also include: space reconstruction module 203, for when in the second surface temperature information there are when invalid data, it is right
Second surface temperature information carries out space reconstruction;Optionally, information Fusion Module 202 may include: the second information fusion unit
2023, for the space scale according to the first surface temperature information, to second after the first surface temperature information and space reconstruction
Surface temperature information carries out temporal-spatial fusion, calculates the fusion surface temperature information of target area.
Optionally, the second surface temperature information is daily surface temperature data.
In above-mentioned surface temperature information acquisition device, by the first surface temperature information and the second ground that obtain target area
Table temperature information, the first surface temperature information are the surface temperature data of each day of clear sky;Then according to the first surface temperature information
Space scale, temporal-spatial fusion is carried out to the first surface temperature information and the second surface temperature information, calculates melting for target area
Close surface temperature information;Because the space scale of the first surface temperature information is less than the space scale of the second surface temperature information,
Therefore the first surface temperature information spatial resolution with higher, the second surface temperature information time resolution with higher
Rate, fusion surface temperature information simultaneously have the first surface temperature information spatial resolution and the second surface temperature information when
Between resolution ratio, i.e., with high time resolution while have high spatial resolution;It is prior, because of the first earth's surface used
Temperature information corresponds to each day of clear sky, avoids because of invalid data caused by cloud pollution etc., avoids because of the weight such as space interpolation
The extra error that construction method introduces, therefore ensure that the spatial resolution of fusion surface temperature information, improve fusion earth's surface temperature
Spend the accuracy of information.
Specific restriction about surface temperature information acquisition device may refer to above for surface temperature acquisition of information
The restriction of method, details are not described herein.Modules in above-mentioned surface temperature information acquisition device can be fully or partially through
Software, hardware and combinations thereof are realized.Above-mentioned each module can be embedded in the form of hardware or independently of the place in computer equipment
It manages in device, can also be stored in a software form in the memory in computer equipment, in order to which processor calls execution or more
The corresponding operation of modules.
In one embodiment, a kind of computer equipment is provided, which can be terminal, internal structure
Figure can be as shown in figure 13.The computer equipment includes the processor connected by system bus, memory, network interface, shows
Display screen and input unit.Wherein, the processor of the computer equipment is for providing calculating and control ability.The computer equipment
Memory includes non-volatile memory medium, built-in storage.The non-volatile memory medium is stored with operating system and computer
Program.The built-in storage provides environment for the operation of operating system and computer program in non-volatile memory medium.The meter
The network interface for calculating machine equipment is used to communicate with external terminal by network connection.When the computer program is executed by processor
To realize a kind of surface temperature information acquisition method.The display screen of the computer equipment can be liquid crystal display or electronic ink
Water display screen, the input unit of the computer equipment can be the touch layer covered on display screen, be also possible to computer equipment
Key, trace ball or the Trackpad being arranged on shell can also be external keyboard, Trackpad or mouse etc..
It will be understood by those skilled in the art that structure shown in Figure 13, only part relevant to application scheme
The block diagram of structure, does not constitute the restriction for the computer equipment being applied thereon to application scheme, and specific computer is set
Standby may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.
In one embodiment, a kind of computer equipment, including memory and processor are provided, is stored in memory
Computer program, the processor perform the steps of when executing computer program
Obtain the multi-source surface temperature information of target area;Multi-source surface temperature information is believed including at least the first surface temperature
Breath and the second surface temperature information;First surface temperature information is the surface temperature data of each day of clear sky;First surface temperature letter
Space scale of the space scale of breath less than the second surface temperature information;
According to the space scale of the first surface temperature information, to the first surface temperature information and the second surface temperature information into
Row temporal-spatial fusion calculates the fusion surface temperature information of target area.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated
Machine program performs the steps of when being executed by processor
Obtain the multi-source surface temperature information of target area;Multi-source surface temperature information is believed including at least the first surface temperature
Breath and the second surface temperature information;First surface temperature information is the surface temperature data of each day of clear sky;First surface temperature letter
Space scale of the space scale of breath less than the second surface temperature information;
According to the space scale of the first surface temperature information, to the first surface temperature information and the second surface temperature information into
Row temporal-spatial fusion calculates the fusion surface temperature information of target area.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with
Relevant hardware is instructed to complete by computer program, the computer program can be stored in a non-volatile computer
In read/write memory medium, the computer program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein,
To any reference of memory, storage, database or other media used in each embodiment provided herein,
Including non-volatile and/or volatile memory.Nonvolatile memory may include read-only memory (ROM), programming ROM
(PROM), electrically programmable ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include
Random access memory (RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms,
Such as static state RAM (SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhancing
Type SDRAM (ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM
(RDRAM), direct memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..
Each technical characteristic of above embodiments can be combined arbitrarily, for simplicity of description, not to above-described embodiment
In each technical characteristic it is all possible combination be all described, as long as however, the combination of these technical characteristics be not present lance
Shield all should be considered as described in this specification.
The several embodiments of the application above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously
It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art
It says, without departing from the concept of this application, various modifications and improvements can be made, these belong to the protection of the application
Range.Therefore, the scope of protection shall be subject to the appended claims for the application patent.
Claims (13)
1. a kind of surface temperature information acquisition method, which is characterized in that the described method includes:
Obtain the multi-source surface temperature information of target area;The multi-source surface temperature information is believed including at least the first surface temperature
Breath and the second surface temperature information;The first surface temperature information is the surface temperature data of each day of clear sky;First ground
The space scale of table temperature information is less than the space scale of the second surface temperature information;
According to the space scale of the first surface temperature information, to the first surface temperature information and the second earth's surface temperature
It spends information and carries out temporal-spatial fusion, calculate the fusion surface temperature information of the target area.
2. the method according to claim 1, wherein the space ruler according to the first surface temperature information
Degree carries out temporal-spatial fusion to the first surface temperature information and the second surface temperature information, calculates the target area
Fusion surface temperature information, comprising:
Determine that the object time corresponding first refers to the day of clear sky with reference to the day of clear sky and second from each day of clear sky;Wherein, institute
It is before the object time and with the object time apart from the nearest date that first, which is stated, with reference to the day of clear sky, described second
It is after the object time and with the object time apart from the nearest date with reference to the day of clear sky;
First refers to clear sky according to the space scale of the first surface temperature information, the first surface temperature information
Day corresponding surface temperature data and described second are with reference to the day of clear sky corresponding surface temperature data, second surface temperature
First with reference to the day of clear sky corresponding surface temperature data and the second reference day of clear sky corresponding earth's surface temperature described in information
Object time corresponding surface temperature data described in degree evidence, the second surface temperature information, determine the target area
In the object time corresponding fusion surface temperature information.
3. according to the method described in claim 2, it is characterized in that, the space ruler according to the first surface temperature information
First with reference to the day of clear sky corresponding surface temperature data and second reference described in degree, the first surface temperature information
The day of clear sky corresponding surface temperature data, described in the second surface temperature information first with reference to the day of clear sky corresponding earth's surface temperature
Degree is accordingly and described second refers to target described in the day of clear sky corresponding surface temperature data, the second surface temperature information
Time corresponding surface temperature data, determine the target area in the object time corresponding fusion surface temperature information,
Include:
First refers to clear sky according to the space scale of the first surface temperature information, the first surface temperature information
Day corresponding surface temperature data and described second are with reference to the day of clear sky corresponding surface temperature data, second surface temperature
First with reference to the day of clear sky corresponding surface temperature data and the second reference day of clear sky corresponding earth's surface temperature described in information
Object time corresponding surface temperature data described in degree evidence, the second surface temperature information, determine the target area
In the reference surface temperature information of the first reference day of clear sky and the target area described second with reference to the day of clear sky
With reference to surface temperature information;
First refers to clear sky according to the space scale of the first surface temperature information, the second surface temperature information
Day corresponding surface temperature data and described second are with reference to the day of clear sky corresponding surface temperature data, second surface temperature
The corresponding surface temperature data of object time described in information determine the target area described first with reference to the of the day of clear sky
The second time weighting that one time weighting and the target area refer to the day of clear sky described second;
The reference surface temperature information and the first time weight that the target area is referred into the day of clear sky described first
Product, the target area refer to the reference surface temperature information of the day of clear sky and multiplying for second time weighting described second
Product carries out sum operation, obtains the target area in the object time corresponding fusion surface temperature information.
4. according to the method described in claim 3, it is characterized in that, the space ruler according to the first surface temperature information
First with reference to the day of clear sky corresponding surface temperature data and second reference described in degree, the first surface temperature information
The day of clear sky corresponding surface temperature data, described in the second surface temperature information first with reference to the day of clear sky corresponding earth's surface temperature
Degree is accordingly and described second refers to target described in the day of clear sky corresponding surface temperature data, the second surface temperature information
Time corresponding surface temperature data determine that reference surface temperature of the target area in the first reference day of clear sky is believed
The reference surface temperature information that breath and the target area refer to the day of clear sky described second, comprising:
According to the space scale of the first surface temperature information, the target area is divided at least one goal pels;
First with reference to the day of clear sky corresponding surface temperature data, second ground according to the first surface temperature information
First refers to the day of clear sky corresponding surface temperature data and the object time corresponding earth's surface temperature described in table temperature information
Degree evidence determines the goal pels in the reference surface temperature information of the first reference day of clear sky;
Second with reference to the day of clear sky corresponding surface temperature data, second ground according to the first surface temperature information
Second refers to the day of clear sky corresponding surface temperature data and the object time corresponding earth's surface temperature described in table temperature information
Degree evidence determines the goal pels in the reference surface temperature information of the second reference day of clear sky;
According to the goal pels in the reference surface temperature information of the first reference day of clear sky, determine that the target area exists
The reference surface temperature information of the first reference day of clear sky;And the day of clear sky is referred to described second according to the goal pels
Reference surface temperature information, determine the target area it is described second with reference to the day of clear sky reference surface temperature information.
5. according to the method described in claim 4, it is characterized in that, the space ruler according to the first surface temperature information
First with reference to the day of clear sky corresponding surface temperature data and second reference described in degree, the second surface temperature information
The day of clear sky corresponding surface temperature data, object time corresponding surface temperature number described in the second surface temperature information
According to, determine the target area it is described first with reference to the day of clear sky first time weight and the target area described
The second time weighting of the second reference day of clear sky, comprising:
First refers to the day of clear sky corresponding surface temperature data and the mesh according to the second surface temperature information
Mark the second reference day of clear sky corresponding earth's surface temperature described in time corresponding surface temperature data, the second surface temperature information
Degree evidence determines the goal pels in the first time weight of the first reference day of clear sky;
Second refers to the day of clear sky corresponding surface temperature data and the mesh according to the second surface temperature information
Mark the first reference day of clear sky corresponding earth's surface temperature described in time corresponding surface temperature data, the second surface temperature information
Degree evidence determines the goal pels in the second time weighting of the second reference day of clear sky.
6. according to the method described in claim 5, it is characterized in that, described refer to clear sky described first for the target area
Day reference surface temperature information and the product of the first time weight, the target area described second refer to the day of clear sky
Reference surface temperature information and second time weighting product carry out sum operation, obtain the target area described
Object time corresponding fusion surface temperature information, comprising:
By the goal pels it is described first with reference to the day of clear sky reference surface temperature information with the goal pels described
The first reference earth's surface for referring to the day of clear sky described second with reference to the product of the first time weight of the day of clear sky, the goal pels
Temperature information and the goal pels are obtained in the product progress sum operation of the second time weighting of the second reference day of clear sky
To the goal pels in the object time corresponding fusion surface temperature information;
According to the goal pels in the object time corresponding fusion surface temperature information, determine the target area in institute
State object time corresponding fusion surface temperature information.
7. according to the method described in claim 4, it is characterized in that, the determination goal pels are described first with reference to fine
The reference surface temperature information of empty day, alternatively, determining the goal pels in the reference earth's surface temperature of the second reference day of clear sky
Spend information, comprising:
Using comprisingRelational expression, be calculated
The goal pels described first with reference to the day of clear sky reference surface temperature information, alternatively, the goal pels are described the
The reference surface temperature information of the two reference days of clear sky;
Wherein, tpFor the object time, k=m or n, tmThe day of clear sky, t are referred to for described firstnFor described second with reference to fine
Empty day, (xw/2,yw/2) be the goal pels coordinate, w is the size of search window centered on the goal pels, N
For the number of the similar pixel of goal pels described in described search window, WiIt is described similar for i-th in described search window
The pixel weight of pixel, ViFor the conversion coefficient of i-th of similar pixel in described search window;
L(xw/2,yw/2,tk) be the goal pels described in the first surface temperature information first with reference to the day of clear sky or
Described second with reference to the day of clear sky corresponding surface temperature data, M (xi,yi,tp) it is i-th of similar pixel in described search window
The corresponding surface temperature data of the object time described in the second surface temperature information, M (xi,yi,tk) it is described search
I-th of similar pixel, first reference day of clear sky described in the second surface temperature information or second reference in window
The day of clear sky corresponding surface temperature data.
8. according to the method described in claim 5, it is characterized in that, the determination goal pels are described first with reference to fine
The first time weight of empty day, alternatively, determining the goal pels in the second time weighting of the second reference day of clear sky, packet
It includes:
Using comprisingRelational expression, the target is calculated
First time weight of the pixel in the first reference day of clear sky, alternatively, determining the goal pels described second with reference to fine
The second time weighting of empty day;
Wherein, tpFor the object time, k=m or n, tmThe day of clear sky, t are referred to for described firstnFor described second with reference to fine
Empty day, w are the size of the search window centered on the goal pels;
M(xj,yi,tp) it is that coordinate is (x in described search windowj,yi) pixel described in the second surface temperature information
Object time corresponding surface temperature data, M (xj,yi,tk) it is that coordinate is (x in described search windowj,yi) pixel described
First with reference to the day of clear sky or the second reference day of clear sky corresponding surface temperature data described in second surface temperature information.
9. the method according to claim 1, wherein in the space according to the first surface temperature information
Scale carries out temporal-spatial fusion to the first surface temperature information and the second surface temperature information, calculates the target area
Before the fusion surface temperature information in domain, the method also includes:
When there are when invalid data, carrying out space weight to the second surface temperature information in the second surface temperature information
It builds;
The space scale according to the first surface temperature information, to the first surface temperature information and second ground
Table temperature information carries out temporal-spatial fusion, calculates the fusion surface temperature information of the target area, comprising:
According to the space scale of the first surface temperature information, after the first surface temperature information and the space reconstruction
The second surface temperature information carry out temporal-spatial fusion, calculate the fusion surface temperature information of the target area.
10. method according to any one of claim 1 to 9, which is characterized in that the second surface temperature information is every
The surface temperature data of day.
11. a kind of surface temperature information acquisition device, which is characterized in that described device includes:
Data obtaining module, for obtaining the multi-source surface temperature information of target area;The multi-source surface temperature information is at least
Including the first surface temperature information and the second surface temperature information;The first surface temperature information is the earth's surface temperature of each day of clear sky
Degree evidence;The space scale of the first surface temperature information is less than the space scale of the second surface temperature information;
Information Fusion Module believes first surface temperature for the space scale according to the first surface temperature information
Breath and the second surface temperature information carry out temporal-spatial fusion, calculate the fusion surface temperature information of the target area.
12. a kind of computer equipment, including memory and processor, the memory are stored with computer program, feature exists
In the step of processor realizes any one of claims 1 to 10 the method when executing the computer program.
13. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The step of method described in any one of claims 1 to 10 is realized when being executed by processor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811002215.1A CN109186774B (en) | 2018-08-30 | 2018-08-30 | Surface temperature information acquisition method, device, computer equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811002215.1A CN109186774B (en) | 2018-08-30 | 2018-08-30 | Surface temperature information acquisition method, device, computer equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109186774A true CN109186774A (en) | 2019-01-11 |
CN109186774B CN109186774B (en) | 2019-11-22 |
Family
ID=64917258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811002215.1A Active CN109186774B (en) | 2018-08-30 | 2018-08-30 | Surface temperature information acquisition method, device, computer equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109186774B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348107A (en) * | 2019-07-08 | 2019-10-18 | 中国农业科学院农业资源与农业区划研究所 | A kind of pixel real surface temperature rebuilding method under cloud |
CN111309778A (en) * | 2020-01-19 | 2020-06-19 | 中国科学院地理科学与资源研究所 | Climate mode evaluation and deviation correction method based on information entropy |
CN111487598A (en) * | 2020-03-26 | 2020-08-04 | 清华大学 | Ice layer thickness calculation method and device, computer equipment and storage medium |
CN114087731A (en) * | 2021-11-10 | 2022-02-25 | 珠海格力电器股份有限公司 | Self-adaptive adjusting system based on coupling of infrared identification module and air conditioner temperature sensing module |
CN114218756A (en) * | 2021-11-24 | 2022-03-22 | 中国农业科学院农业资源与农业区划研究所 | Subsurface surface temperature reconstruction method based on surface temperature annual change model |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101629850A (en) * | 2009-08-24 | 2010-01-20 | 中国农业科学院农业资源与农业区划研究所 | Method for inversing land surface temperature from MODIS data |
CN101908196A (en) * | 2009-06-03 | 2010-12-08 | 北京师范大学 | Method for estimating vegetation coverage based on vegetation-soil moisture response relation |
CN102103077A (en) * | 2009-12-16 | 2011-06-22 | 中国科学院沈阳应用生态研究所 | MODIS data-based agricultural drought monitoring method |
CN102915529A (en) * | 2012-10-15 | 2013-02-06 | 黄波 | Integrated fusion technique and system based on remote sensing of time, space, spectrum and angle |
CN103177431A (en) * | 2012-12-26 | 2013-06-26 | 中国科学院遥感与数字地球研究所 | Method of spatial-temporal fusion for multi-source remote sensing data |
CN104236716A (en) * | 2013-06-17 | 2014-12-24 | 中国科学院城市环境研究所 | Land surface temperature inversion method based on spatio-temporal information of paired HJ-1B images |
CN105184076A (en) * | 2015-09-02 | 2015-12-23 | 安徽大学 | Space-time integrated fusion method for remote sensing earth surface temperature data |
CN105868529A (en) * | 2016-03-18 | 2016-08-17 | 北京师范大学 | Near-surface daily mean atmospheric temperature retrieval method based on remote control |
US20160299910A1 (en) * | 2015-04-09 | 2016-10-13 | Umm Al-Qura University | Method and system for querying and visualizing satellite data |
CN106933776A (en) * | 2017-03-02 | 2017-07-07 | 宁波大学 | A kind of method that MODIS AOD products missing data is repaired |
CN107610050A (en) * | 2017-08-31 | 2018-01-19 | 北京航天绘景科技有限公司 | A kind of high-spatial and temporal resolution land table Temperature fusion method based on multi-source satellite data |
CN108268735A (en) * | 2018-01-29 | 2018-07-10 | 浙江大学 | Soil moisture NO emissions reduction method based on multi-source remote sensing satellite fused data |
CN108446999A (en) * | 2018-04-17 | 2018-08-24 | 中国水利水电科学研究院 | Irrigated area Different Crop ET evaluation methods are carried out based on canopy-air temperature difference and remote sensing information |
-
2018
- 2018-08-30 CN CN201811002215.1A patent/CN109186774B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908196A (en) * | 2009-06-03 | 2010-12-08 | 北京师范大学 | Method for estimating vegetation coverage based on vegetation-soil moisture response relation |
CN101629850A (en) * | 2009-08-24 | 2010-01-20 | 中国农业科学院农业资源与农业区划研究所 | Method for inversing land surface temperature from MODIS data |
CN102103077A (en) * | 2009-12-16 | 2011-06-22 | 中国科学院沈阳应用生态研究所 | MODIS data-based agricultural drought monitoring method |
CN102915529A (en) * | 2012-10-15 | 2013-02-06 | 黄波 | Integrated fusion technique and system based on remote sensing of time, space, spectrum and angle |
CN103177431A (en) * | 2012-12-26 | 2013-06-26 | 中国科学院遥感与数字地球研究所 | Method of spatial-temporal fusion for multi-source remote sensing data |
CN104236716A (en) * | 2013-06-17 | 2014-12-24 | 中国科学院城市环境研究所 | Land surface temperature inversion method based on spatio-temporal information of paired HJ-1B images |
US20160299910A1 (en) * | 2015-04-09 | 2016-10-13 | Umm Al-Qura University | Method and system for querying and visualizing satellite data |
CN105184076A (en) * | 2015-09-02 | 2015-12-23 | 安徽大学 | Space-time integrated fusion method for remote sensing earth surface temperature data |
CN105868529A (en) * | 2016-03-18 | 2016-08-17 | 北京师范大学 | Near-surface daily mean atmospheric temperature retrieval method based on remote control |
CN106933776A (en) * | 2017-03-02 | 2017-07-07 | 宁波大学 | A kind of method that MODIS AOD products missing data is repaired |
CN107610050A (en) * | 2017-08-31 | 2018-01-19 | 北京航天绘景科技有限公司 | A kind of high-spatial and temporal resolution land table Temperature fusion method based on multi-source satellite data |
CN108268735A (en) * | 2018-01-29 | 2018-07-10 | 浙江大学 | Soil moisture NO emissions reduction method based on multi-source remote sensing satellite fused data |
CN108446999A (en) * | 2018-04-17 | 2018-08-24 | 中国水利水电科学研究院 | Irrigated area Different Crop ET evaluation methods are carried out based on canopy-air temperature difference and remote sensing information |
Non-Patent Citations (4)
Title |
---|
CHAO ZENG: "Recovering missing pixels for Landsat ETM+ SLC-off imagery using multi-temporal regression analysis and a regularization method", 《REMOTE SENSING OF ENVIRONMENT》 * |
LIANGLIANG BAI: "Responses of field evapotranspiration to the changes of cropping pattern and groundwater depth in large irrigation district of Yellow River basin", 《AGRICULTURAL WATER MANAGEMENT》 * |
XIAOLIN ZHU: "An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions", 《REMOTE SENSING OF ENVIRONMENT》 * |
于红博: "皇甫川流域地表特征参数遥感反演研究", 《人民黄河》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348107A (en) * | 2019-07-08 | 2019-10-18 | 中国农业科学院农业资源与农业区划研究所 | A kind of pixel real surface temperature rebuilding method under cloud |
CN110348107B (en) * | 2019-07-08 | 2022-06-24 | 中国农业科学院农业资源与农业区划研究所 | Method for reconstructing real surface temperature of pixel under cloud |
CN111309778A (en) * | 2020-01-19 | 2020-06-19 | 中国科学院地理科学与资源研究所 | Climate mode evaluation and deviation correction method based on information entropy |
CN111309778B (en) * | 2020-01-19 | 2023-05-05 | 中国科学院地理科学与资源研究所 | Climate mode evaluation and deviation correction method based on information entropy |
CN111487598A (en) * | 2020-03-26 | 2020-08-04 | 清华大学 | Ice layer thickness calculation method and device, computer equipment and storage medium |
CN111487598B (en) * | 2020-03-26 | 2021-02-12 | 清华大学 | Ice layer thickness calculation method and device, computer equipment and storage medium |
CN114087731A (en) * | 2021-11-10 | 2022-02-25 | 珠海格力电器股份有限公司 | Self-adaptive adjusting system based on coupling of infrared identification module and air conditioner temperature sensing module |
CN114087731B (en) * | 2021-11-10 | 2022-10-21 | 珠海格力电器股份有限公司 | Self-adaptive adjusting system based on coupling of infrared identification module and air conditioner temperature sensing module |
CN114218756A (en) * | 2021-11-24 | 2022-03-22 | 中国农业科学院农业资源与农业区划研究所 | Subsurface surface temperature reconstruction method based on surface temperature annual change model |
Also Published As
Publication number | Publication date |
---|---|
CN109186774B (en) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109186774B (en) | Surface temperature information acquisition method, device, computer equipment and storage medium | |
CN109115995B (en) | Soil moisture information acquisition methods, device, computer equipment and storage medium | |
Xia et al. | Combining kernel-driven and fusion-based methods to generate daily high-spatial-resolution land surface temperatures | |
Hu et al. | How can we use MODIS land surface temperature to validate long‐term urban model simulations? | |
Weng | Advances in environmental remote sensing: sensors, algorithms, and applications | |
Qiu et al. | A comparison of NDVI and EVI in the DisTrad model for thermal sub-pixel mapping in densely vegetated areas: A case study in Southern China | |
Hu et al. | A framework for analyzing and designing scale invariant remote sensing algorithms | |
CN111368817B (en) | Method and system for quantitatively evaluating thermal effect based on earth surface type | |
Yao et al. | LAI retrieval and uncertainty evaluations for typical row-planted crops at different growth stages | |
Ma et al. | Generating gapless land surface temperature with a high spatio-temporal resolution by fusing multi-source satellite-observed and model-simulated data | |
Zhao et al. | Downscaling of soil moisture products using deep learning: Comparison and analysis on Tibetan Plateau | |
Rahman et al. | A comparison of four relative radiometric normalization (RRN) techniques for mosaicing H-res multi-temporal thermal infrared (TIR) flight-lines of a complex urban scene | |
Zhou et al. | A comparison of three methods for estimating the LAI of black locust (Robinia pseudoacacia L.) plantations on the Loess Plateau, China | |
Wang et al. | Spatial–spectral radial basis function-based interpolation for Landsat ETM+ SLC-off image gap filling | |
Bhattarai et al. | Validation of evaporation estimates from a modified surface energy balance algorithm for land (SEBAL) model in the south-eastern United States | |
Shen et al. | Spatial–temporal land-use/land-cover dynamics and their impacts on surface temperature in Chongming Island of Shanghai, China | |
CN110503137A (en) | Based on the determination method of the remote sensing image temporal-spatial fusion base image pair of mixing together | |
Propastin | Large-scale mapping of aboveground biomass of tropical rainforest in Sulawesi, Indonesia, using Landsat ETM+ and MODIS data | |
Landier et al. | Calibration of urban canopies albedo and 3D shortwave radiative budget using remote-sensing data and the DART model | |
Sheng et al. | Quantifying the spatial and temporal relationship between air and land surface temperatures of different land-cover types in Southeastern China | |
Arias-Rosales et al. | Shadow modeling in urban environments for solar harvesting devices with freely defined positions and orientations | |
Salehi et al. | Spatial and temporal resolution improvement of actual evapotranspiration maps using Landsat and MODIS data fusion | |
Shi | Thermal infrared inverse model for component temperatures of mixed pixels | |
Chen et al. | Definition of component effective emissivity for heterogeneous and non-isothermal surfaces and its approximate calculation | |
Du et al. | Modeling directional brightness temperature of the winter wheat canopy at the ear stage |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |