CN106599548A - Spatial and temporal scale matching method and device for land surface evapotranspiration remote sensing estimation - Google Patents
Spatial and temporal scale matching method and device for land surface evapotranspiration remote sensing estimation Download PDFInfo
- Publication number
- CN106599548A CN106599548A CN201611047244.0A CN201611047244A CN106599548A CN 106599548 A CN106599548 A CN 106599548A CN 201611047244 A CN201611047244 A CN 201611047244A CN 106599548 A CN106599548 A CN 106599548A
- Authority
- CN
- China
- Prior art keywords
- matching
- time
- remote sensing
- data
- remotely
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Image Processing (AREA)
Abstract
The invention provides a spatial and temporal scale matching method and device for land surface evapotranspiration remote sensing estimation. The method comprises the steps of: S1: acquiring a remote sensing data set with multiple spatial resolutions and multiple temporal resolutions; S2: according to the acquired remote sensing data set, carrying out spatial and temporal data fusion to obtain a remote sensing data set which simultaneously has a high temporal resolution and a high spatial resolution; S3: carrying out ET parametric inversion on the data set generated in the steps S1 and S2 by utilizing a preset data model so as to obtain a remote sensing parameter set with different temporal and spatial resolutions; S4: according to remote sensing parameters with different temporal and spatial resolutions, which are acquired in the step S3, respectively carrying out spatial and temporal scale matching judgment; and if a spatial scale matching judgment result and a temporal scale matching judgment result meet the accuracy requirement of ET estimation, using remote sensing data in the step S4, which meets the accuracy requirement of ET estimation, as a matched data set. The spatial and temporal scale matching method and device provided by the invention can solve the problem of uncertainty of ET inversion, which is caused by temporal and spatial scale matching uncertainty in existing ET estimation.
Description
Technical field
The present invention relates to remote sensing and field of ecological hydrology, and in particular to a kind of land evapotranspiration remote sensing appraising when
Empty yardstick matching process and device.
Background technology
It is method conventional at present by Remote sensing parameters inverting and Model coupling simulation estimate evapotranspiration (ET), compares tradition
Hydrological model, remote sensing ET models in precision and spatially have obvious advantage.But its temporal discreteness, underlying surface
Isomerism and the comprehensive of spectral information limit its observation to the continuous hydrology.Therefore, " time ", " space " matching are to reduce
One of uncertain and the raising suitability key issue of ET invertings.
Yet with the limitation of single-sensor, it is difficult to take into account different times, space scale matching correct.At present, ET
The yardstick matching correct of inverting is also extremely limited, and part research has done certain with reference to data fusion, parameter simulation etc. aspect
Exploratory development, but not enough system, it is to carry out the probabilistic comprehensive characterization research of yardstick even more extremely limited.
The content of the invention
For defect of the prior art, the present invention provides a kind of spatial and temporal scales match party of land evapotranspiration remote sensing appraising
Method and device, to solve prior art in because cannot take into account the time, space scale matching and caused ET invertings are probabilistic
Problem.
To solve above-mentioned technical problem, the present invention provides technical scheme below：
In a first aspect, the invention provides a kind of spatial and temporal scales matching process of land evapotranspiration remote sensing appraising, including such as
Lower step：
S1, the remotely-sensed data collection that many spaces, many temporal resolutions are obtained using sensor；
S2, the remotely-sensed data collection obtained according to step S1 carry out space-time data fusion, obtain having high time and height simultaneously
The remotely-sensed data collection of spatial resolution；
S3, the data set that step S1 and S2 are generated is carried out into ET parametric inversions using preset data model, obtain ET parameters
Collection, obtains the Remote sensing parameters collection of different time and spatial resolution；
The Remote sensing parameters of S4, the different time obtained according to step S3 and spatial resolution, carry out respectively space scale
With judgement and time scale matching judgment；
If S5, space scale matching judgment result and time scale matching judgment result meet the required precision of ET estimations,
The data set of the remotely-sensed data of ET estimation precision requirements as matching then will be met in step S4.
Further, the different time for being obtained according to step S3 and the Remote sensing parameters of spatial resolution, carry out space scale
Matching judgment, specifically includes：
S41, the land-use map for obtaining the generation of different resolution remote sensing image, by land-use map vector quantization, will be identical
The land-use map of the vector quantization that resolution remote sensing images are obtained is superimposed with the remote sensing land-use map of grid, calculates space loose
Degree index：SP=Ni/Si；
In formula, SP is space shatter value index；NiFor the Land_use change patch number generated under pre-set spatial resolution；SiFor
The area of above-mentioned Land_use change；
S42, calculating underlying surface heterogeneity index：
In formula, A be underlying surface heterogeneity index, AaFor selected remotely-sensed data pixel sum；M is mixed pixel judgement
Value, when for mixed pixel when, be worth for 1, be otherwise 0；SmIt is the ground included in mixed pixel for correspondence pixel coefficient of similarity
The species number of species type；
S43, as SP >=a, while during A≤b, being judged as space scale matched；
When only meeting SP >=a, or when only meeting A≤b, it is judged as that space scale moderate is matched；
When the two is unsatisfactory for, it is judged as that space scale is mismatched；
In formula, a and b is respectively corresponding predetermined threshold value.
Further, the different time for being obtained according to step S3 and the Remote sensing parameters of spatial resolution, carry out time scale
Matching judgment, specifically includes：
S41 ', acquisition are in t-t0The period of duration curve data of interval difference evapotranspiration parameter, and judge the curve data
And whether the dependency between the curve of satellite actual measurement meets：R2>=c, RMSE≤d, P≤0.01；In formula, c, d are default threshold
Value；t-t0For time interval；R is correlation coefficient；RMSE is root-mean-square error；P values are significance level discriminant value；
S42 ', it is segmented according to vegetation growing period or Crop growing stage, on the basis of here segmentation, is calculated different numbers
According to retrievable parametric stability index VP between the time periodz：VPz=VARtb-tc(Pztb,Pztb+1,Pztb+2,......,Pztc),
And judge VPzWhether value levels off to 0；Wherein, PzFor parameter value, tb ... tc are different time point values, VARtb-tcFor when
Between interval tb to tc variance computation model, VPzFor the variance yields of tb to tc；
S43 ' if, have one to meet condition in S41 ' and S42 ', be judged as that time scale is matched, meet the essence of ET estimations
Degree is required.
Further, if space scale and time scale matching judgment result are matching, the difference that step S3 is obtained
The Remote sensing parameters of time and spatial resolution, as the data set for meeting ET required precisions, are as a result matching, are specifically included：
If space scale matching result is matched, time scale matching result is matching, then it is assumed that meet ET estimations
Required precision, now will meet the data set of the remotely-sensed data of ET estimation precision requirements as matching in step S4；
If space scale matching result is moderate matching, time scale matching result is matching, then by spatial match and
The parameter set of time match is input in ET models, if the ET values relation for obtaining the ET estimation results and actual observation simulated meets
Condition：R2>=e, RMSE≤f, during P≤g, then it is assumed that meet ET estimation precision requirements, now will meet ET estimation essences in step S4
Data set of the remotely-sensed data that degree is required as matching；In formula, e, f, g are corresponding predetermined threshold value.
Further, the remotely-sensed data collection for being obtained according to step S1 carries out space-time data fusion, when obtaining having high simultaneously
Between and high spatial resolution remotely-sensed data collection, specifically include：
Using S-G Filtering Models and temporal-spatial fusion model, by many spaces in step S1, the remote sensing number of many temporal resolutions
Space-time data fusion is carried out according to collection, obtains that there is high time, the remotely-sensed data collection of high spatial resolution simultaneously.
Second aspect, present invention also offers a kind of spatial and temporal scales coalignment of land evapotranspiration remote sensing appraising, including：
Acquisition module, for obtaining the remotely-sensed data of many spaces, many temporal resolutions using sensor；
Fusion Module, for processing the data set that model generates the data acquisition module using preset data, is carried out
High time and the data fusion of high spatial, obtain the remotely-sensed data collection with high time and high spatial resolution；
Generation module, the data set for being obtained according to the acquisition module and Fusion Module carries out ET parametric inversions, obtains
To ET parameter sets, the Remote sensing parameters collection of different time and spatial resolution is obtained；
Matching judgment module, the remote sensing of different time and spatial resolution for being generated according to the generation module is joined
Number, carries out respectively space scale matching judgment and time scale matching judgment；
Matching result processing module, for determining space scale matching judgment result and time in the matching judgment module
When yardstick matching judgment result is satisfied by the required precision of ET estimations, will meet the remotely-sensed data of ET estimation precision requirements as
The data set matched somebody with somebody.
Further, the matching judgment module, specifically for：
S41, the land-use map for obtaining the generation of different resolution remote sensing image, by land-use map vector quantization, by vector
The land-use map of change is superimposed with the Remote sensing parameters of the grid of equal resolution, calculates space shatter value index：SP=Ni/Si；
In formula, SP is space shatter value index；NiFor the Land_use change patch number generated under pre-set spatial resolution；SiFor
The area of above-mentioned land-use map；
S42, calculating underlying surface heterogeneity index：
In formula, A be underlying surface heterogeneity index, AaFor selected remotely-sensed data pixel sum；M is mixed pixel judgement
Value, when for mixed pixel when, value is 1, is otherwise 0；SmIt is included in mixed pixel for correspondence pixel coefficient of similarity
The species number of type of ground objects；
S43, as SP >=a, while during A≤b, being judged as space scale matched；
When only meeting SP >=a, or when only meeting A≤b, it is judged as that space scale is matched；
When the two is unsatisfactory for, it is judged as that space scale is mismatched；
In formula, a and b is respectively corresponding predetermined threshold value.
Further, the matching judgment module, specifically for：
S41 ', acquisition are in t-t0The period of duration curve data of interval difference evapotranspiration parameter, and judge the curve data
And whether the dependency between the curve of satellite actual measurement meets：R2>=c, RMSE≤d, P≤0.01；In formula, c, d are default threshold
Value；t-t0For time interval；
S42 ', it is segmented according to vegetation growing period or Crop growing stage, on the basis of here segmentation, is calculated different numbers
According to retrievable parametric stability index VP between the time periodz：VPz=VARtb-tc(Pztb,Pztb+1,Pztb+2,......,Pztc),
And judge VPzWhether value levels off to 0；Wherein, PzFor parameter value, tb ... tc are different time point values, VARtb-tcFor side
Difference calculates function, VPzFor the variance yields of tb to tc；
S43 ' if, have one to meet condition in S41 ' and S42 ', be judged as that time scale is matched.
Further, the matching result processing module, specifically for：
Determine space scale matching result for matched in the matching judgment module, time scale matching result for
Timing, it is believed that meet ET estimation precision requirements, will now meet the data of the remotely-sensed data of ET estimation precision requirements as matching
Collection；
The matching judgment module determine space scale matching result for moderate matching, time scale matching result for
Timing, is input to the relevant parameter model used during spatial match and time match and illegally in ET models, and
Meet condition with the ET value relations of actual observation in the ET estimation results for obtaining simulating：R2>=e, RMSE≤f, during P≤g, it is believed that
Meet ET estimation precision requirements, will now meet the data set of the remotely-sensed data of ET estimation precision requirements as matching；In formula, e,
F, g are corresponding predetermined threshold value.
Further, the Fusion Module, specifically for：
Using S-G Filtering Models and temporal-spatial fusion model, many spaces that acquisition module is obtained, many temporal resolutions it is distant
Sense data set carries out space-time data fusion, obtains having high time, the remotely-sensed data collection of high spatial resolution simultaneously.
As shown from the above technical solution, the spatial and temporal scales matching process of the land evapotranspiration remote sensing appraising that the present invention is provided,
The remotely-sensed data of many spaces, many temporal resolutions is obtained first with sensor；Then according to many spaces, many times point for obtaining
The remotely-sensed data collection of resolution carries out space-time data fusion, obtains the remotely-sensed data with high time and high spatial resolution
Collection；Followed by preset data model by many spaces for obtaining, the remotely-sensed data collection of many temporal resolutions and while when having high
Between and the remotely-sensed data collection of high spatial resolution carry out ET parametric inversions, obtain ET parameter sets, and then obtain different time and sky
Between resolution Remote sensing parameters collection；According to the Remote sensing parameters of the different time and spatial resolution for obtaining, space chi is carried out respectively
Degree matching judgment and time scale matching judgment；If space scale matching judgment result and time scale matching judgment result are full
The required precision of sufficient ET estimations, then will meet the data set of the remotely-sensed data of ET estimation precision requirements as matching.It can be seen that, this
The spatial and temporal scales matching process of the land evapotranspiration remote sensing appraising of bright offer, when can solve the problem that in prior art because taking into account
Between, space scale matching and the probabilistic problem of caused ET invertings.
Description of the drawings
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are the present invention
Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can be with basis
These accompanying drawings obtain other accompanying drawings.
Fig. 1 is the stream of the spatial and temporal scales matching process of the land evapotranspiration remote sensing appraising that first embodiment of the invention is provided
Cheng Tu；
Fig. 2 is the knot of the spatial and temporal scales coalignment of the land evapotranspiration remote sensing appraising that second embodiment of the invention is provided
Structure schematic diagram.
Specific embodiment
To make purpose, technical scheme and the advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, clear, complete description is carried out to the technical scheme in the embodiment of the present invention, it is clear that described embodiment is
The a part of embodiment of the present invention, rather than the embodiment of whole.Based on the embodiment in the present invention, those of ordinary skill in the art
The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
First embodiment of the invention provides a kind of spatial and temporal scales matching process of land evapotranspiration remote sensing appraising, Fig. 1
Show the flow chart of the spatial and temporal scales matching process of the land evapotranspiration remote sensing appraising that first embodiment of the invention is provided.Ginseng
See Fig. 1, methods described comprises the steps：
Step 101：The remotely-sensed data of many spaces, many temporal resolutions is obtained using sensor.
Step 102：Space-time data fusion is carried out according to the remotely-sensed data collection that step 101 is obtained, when obtaining having high simultaneously
Between and high spatial resolution remotely-sensed data collection.
Step 103：Step 101 and 102 data sets for generating are carried out ET parametric inversions by profit using preset data model, are obtained
To ET parameter sets, the Remote sensing parameters collection of different time and spatial resolution is obtained.
In this step, using the remotely-sensed data of different spatial and temporal resolutions, the parameter set of evapotranspiration ET remote-sensing inversions is generated
(land cover pattern/utilization, Albedo, LAI, Coverage, Ts, soil moisture, precipitation etc.)；Compile meteorological data, flux
Data etc.；Determine region underlying surface type, natural vegetation species, natural vegetation coverage condition, crop species and crop phenology etc.
Information, and its phenological period of storage management corresponding spectral information etc., set up priori storehouse.
Step 104：The different time obtained according to step 103 and the Remote sensing parameters of spatial resolution, carry out respectively space
Yardstick matching judgment and time scale matching judgment.
Step 105：If space scale matching judgment result and time scale matching judgment result meet the precision of ET estimations
Require, then the high time for obtaining step 102 and the remotely-sensed data collection of high spatial resolution are used as the data set for matching.
In this step, the required precision for meeting ET estimations as meets matching condition set in advance, specifically can be found in
Subsequent embodiment is introduced.
Seen from the above description, the spatial and temporal scales match party of land evapotranspiration remote sensing appraising provided in an embodiment of the present invention
Method, first with sensor the remotely-sensed data of many spaces, many temporal resolutions is obtained；Then according to obtain many spaces, it is many when
Between the remotely-sensed data collection of resolution carry out space-time data fusion, obtain the remote sensing number with high time and high spatial resolution
According to collection；Followed by preset data model is by many spaces for obtaining, the remotely-sensed data collection of many temporal resolutions and while has height
The remotely-sensed data collection of time and high spatial resolution carries out ET parametric inversions, obtains ET parameter sets, so obtain different time and
The Remote sensing parameters collection of spatial resolution；According to the Remote sensing parameters of the different time and spatial resolution for obtaining, space is carried out respectively
Yardstick matching judgment and time scale matching judgment；If space scale matching judgment result and time scale matching judgment result
Meet the required precision of ET estimations, then will meet the data set of the remotely-sensed data of ET estimation precision requirements as matching.It can be seen that, this
The spatial and temporal scales matching process of the land evapotranspiration remote sensing appraising that invention is provided, when can solve the problem that in prior art because taking into account
Between, space scale matching and the probabilistic problem of caused ET invertings.
In a kind of optional embodiment, above-mentioned steps 104 are specifically included：
S41, the land-use map for obtaining the generation of different resolution remote sensing image, by land-use map vector quantization, will be identical
The land-use map of the vector quantization that resolution remote sensing images are obtained is superimposed with the remote sensing land-use map of grid, calculates space loose
Degree index：SP=Ni/(Si·Pi)；
In formula, SP is space shatter value index；NiFor the Land_use change patch number generated under pre-set spatial resolution；SiFor
The area of above-mentioned Land_use change；
S42, calculating underlying surface heterogeneity index：
In formula, A be underlying surface heterogeneity index, AaFor selected remotely-sensed data pixel sum；M is mixed pixel judgement
Value, when for mixed pixel when, be worth for 1, be otherwise 0；SmIt is the ground included in mixed pixel for correspondence pixel coefficient of similarity
The species number of species type；
S43, as SP >=a, while during A≤b, being judged as space scale matched；
When only meeting SP >=a, or when only meeting A≤b, it is judged as that space scale moderate is matched；
When the two is unsatisfactory for, it is judged as that space scale is mismatched；
In formula, a and b is respectively corresponding predetermined threshold value.
In another kind of optional embodiment, above-mentioned steps 104 are specifically included：
S41 ', acquisition are in t-t0The period of duration curve data of interval difference evapotranspiration parameter, and judge the curve data
And whether the dependency between the curve of satellite actual measurement meets：R2>=c, RMSE≤d, P≤0.01；In formula, c, d are default threshold
Value；t-t0For time interval；R is correlation coefficient；RMSE is root-mean-square error；P values are significance level discriminant value；
S42 ', it is segmented according to vegetation growing period or Crop growing stage, on the basis of here segmentation, is calculated different numbers
According to retrievable parametric stability index VP between the time periodz：VPz=VARtb-tc(Pztb,Pztb+1,Pztb+2,......,Pztc),
And judge VPzWhether value levels off to 0；Wherein, PzFor parameter value, tb ... tc are different time point values, VARtb-tcFor when
Between interval tb to tc variance computation model, VPzFor the variance yields of tb to tc；
S43 ' if, have one to meet condition in S41 ' and S42 ', be judged as that time scale is matched, meet the essence of ET estimations
Degree is required.
In a kind of optional embodiment, above-mentioned steps 105 are specifically included：
If space scale matching result is matched, time scale matching result is matching, then it is assumed that meet space-time
With requiring, the data set of the remotely-sensed data of ET estimation precision requirements as matching now will be met in step 104；
If space scale matching result is moderate matching, time scale matching result is matching, then by spatial match and
The parameter set of time match is input in ET models
If the ET value relations for obtaining the ET estimation results and actual observation simulated meet condition：R2>=e, RMSE≤f, P≤g
When, then it is assumed that meet ET estimation precision requirements, now will meet in step S4 the remotely-sensed data of ET estimation precision requirements as
The data set matched somebody with somebody；In formula, e, f, g are corresponding predetermined threshold value.
In a kind of optional embodiment, above-mentioned steps 102 are specifically included：
Using S-G Filtering Models and temporal-spatial fusion model, by many spaces in step 101, the remote sensing of many temporal resolutions
Data set carries out space-time data fusion, obtains having high time, the remotely-sensed data collection of high spatial resolution simultaneously.
Because the factors such as sensor, cloud layer air affect, inevitably produce some noises, using Savitzky and
The S-G filtering algorithms that Golay is proposed remove noise, and the time series data of its reconstruct can clearly describe the secular change of sequence
Trend and the abrupt information of local, the reconstruct to vegetation index time series data has the preferable suitability
S-G filtering restructing algorithms are as follows：
Wherein,Reconstruct time series data, Tj+1For original temporal data, CiFor filter factor, N be in sliding window when
Ordinal number data bulk.
The data of the high spatial resolution such as landsat lacked by temporal-spatial fusion modeling, are obtained the high time high
The time series data of spatial resolution.
ETARFM(Enhanced Spatial and Temporal A daptive Reflectance Fusion
Model such as landsat and high time resolution reflectivity data) are utilized in the difference of the information such as pixel distance, spectrum, acquisition time
It is different, simulate if MODIS data correspondence phase is such as landsat reflectivity datas.
During pure pixel, t0There is such as ShiShimonoseki in the data of the high spatial resolution of time and high time resolution reflectivity data
System：
Then there is simulation tkThe data of the high spatial resolution of time are：
L and M are respectively high spatial resolution and high time resolution Reflectivity for Growing Season, (xi,yj) it is pixel position, a, b are
There is the coefficient (being caused by waveband width, geometric error etc.) of linear relationship in two sensorses reflectivity data.
In view of 1. in practice pixel mostly is mixed pixel 2. 3. atural object coverage condition senses such as the change of time
Device position can change over time.
Therefore, Fusion Model establishes micro-slip window, the similar neighbouring pixel of spectrum of center pel is found, according to height
The time difference of the high time resolution of the SPECTRAL DIVERSITY of spatial resolution and high time resolution, fiducial time and simulated time
It is different, and the European geometric distance in space of center pel and neighbouring pixel, different weights are given to neighbouring pixel, simulated
Center pel reflectance.
When selecting neighbouring pixel, meet：
Wherein, L be high spatial resolution reflectance, (xi,yj) be pixel position, w be micro-slip window size, t0For
Time, σ (Bn) be the n-th wave band reflectivity data standard deviation.
The foveal reflex rate for then having simulation is：
Wherein,Centered on pixel (xw/2,yw/2) in simulated time tpHigh spatial resolution earth surface reflection
Rate,Centered on pixel (xw/2,yw/2) in fiducial time t0High spatial resolution Reflectivity for Growing Season,Respectively pixel (xi,yj) in t0, tpThe high time resolution Reflectivity for Growing Season of time, ViIt is by mixed
Close the linear coefficient that pixel analysis are obtained, WijkFor weight.
Cijk=Sijk*Tijk*Dijk
Sijk=| L (xi,yj,tk)-M(xi,yj,tk)|
Tijk=| M (xi,yj,tk)-M(xi,yj,t0)|
Sijk(x is located at for giveni,yj) high time resolution and high spatial resolution Reflectivity for Growing Season difference, the parameter
Both SPECTRAL DIVERSITY can be weighed, value is less, represent that given position adjacent to the similar height of pixel, then gives higher weight；TijkTable
Show the difference in reflectivity between twice of high time resolution data, the value is less show the time period in spectrum change more
It is little, higher weight is given in the calculation；DijkCentered on pixel point and the pixel point for participating in calculating geometric distance, the value gets over
The weight of little imparting is higher.
Can obtain further according to phase weighting：
Wherein
Wherein, m, n are two phases, and p simulates phase, and B is wave band.
The embodiment of the present invention simulates high spatial resolution data, that is, be input into Tm,TnThe high spatial resolution and Gao Shi of time
Between resolution reference images and TkThe high time resolution of time, obtains TkResolution data between the simulated altitude of time.By melting
Matched moulds type can obtain the high spatial resolution data for being limited and being lacked due to cloud pollution or revisiting period.
Second embodiment of the invention provides a kind of spatial and temporal scales coalignment of land evapotranspiration remote sensing appraising, referring to
Fig. 2, the device includes：
Acquisition module 21, for obtaining the remotely-sensed data of many spaces, many temporal resolutions using sensor；
Fusion Module 22, for processing the data set that model generates the data acquisition module using preset data, enters
Row high time and the data fusion of high spatial, obtain the remotely-sensed data collection with high time and high spatial resolution；
Generation module 23, the data set for being obtained according to the acquisition module and Fusion Module carries out ET parametric inversions,
ET parameter sets are obtained, the Remote sensing parameters collection of different time and spatial resolution is obtained；
Matching judgment module 24, the remote sensing of different time and spatial resolution for being generated according to the generation module is joined
Number, carries out respectively space scale matching judgment and time scale matching judgment；
Matching result processing module 25, for the matching judgment module determine space scale matching judgment result and when
Between yardstick matching judgment result when being satisfied by the required precision of ET estimations, will meet the remotely-sensed data of ET estimation precision requirements as
The data set of matching.
In a kind of optional embodiment, the matching judgment module 24, specifically for：
S41, the land-use map for obtaining the generation of different resolution remote sensing image, by land-use map vector quantization, by vector
The land-use map of change is superimposed with the Remote sensing parameters of the grid of equal resolution, calculates space shatter value index：SP=Ni/Si；
In formula, SP is space shatter value index；NiFor the Land_use change patch number generated under pre-set spatial resolution；SiFor
The area of above-mentioned land-use map；
S42, calculating underlying surface heterogeneity index：
In formula, A be underlying surface heterogeneity index, AaFor selected remotely-sensed data pixel sum；M is mixed pixel judgement
Value, when for mixed pixel when, value is 1, is otherwise 0；SmIt is included in mixed pixel for correspondence pixel coefficient of similarity
The species number of type of ground objects；
S43, as SP >=a, while during A≤b, being judged as space scale matched；
When only meeting SP >=a, or when only meeting A≤b, it is judged as that space scale is matched；
When the two is unsatisfactory for, it is judged as that space scale is mismatched；
In formula, a and b is respectively corresponding predetermined threshold value.
In another kind of optional embodiment, the matching judgment module 42, specifically for：
S41 ', acquisition are in t-t0The period of duration curve data of interval difference evapotranspiration parameter, and judge the curve data
And whether the dependency between the curve of satellite actual measurement meets：R2>=c, RMSE≤d, P≤0.01；In formula, c, d are default threshold
Value；t-t0For time interval；
S42 ', it is segmented according to vegetation growing period or Crop growing stage, on the basis of here segmentation, is calculated different numbers
According to retrievable parametric stability index VP between the time periodz：VPz=VARtb-tc(Pztb,Pztb+1,Pztb+2,......,Pztc),
And judge VPzWhether value levels off to 0；Wherein, PzFor parameter value, tb ... tc are different time point values, VARtb-tcFor side
Difference calculates function, VPzFor the variance yields of tb to tc；
S43 ' if, have one to meet condition in S41 ' and S42 ', be judged as that time scale is matched.
In a kind of optional embodiment, the matching result processing module 25, specifically for：
Determine space scale matching result for matched in the matching judgment module, time scale matching result for
Timing, it is believed that meet ET estimation precision requirements, will now meet the data of the remotely-sensed data of ET estimation precision requirements as matching
Collection；
The matching judgment module determine space scale matching result for moderate matching, time scale matching result for
Timing, is input to the relevant parameter model used during spatial match and time match and illegally in ET models, and
Meet condition with the ET value relations of actual observation in the ET estimation results for obtaining simulating：R2>=e, RMSE≤f, during P≤g, it is believed that
Meet ET estimation precision requirements, will now meet the data set of the remotely-sensed data of ET estimation precision requirements as matching；In formula, e,
F, g are corresponding predetermined threshold value.
In a kind of optional embodiment, the Fusion Module 22, specifically for：
Using S-G Filtering Models and temporal-spatial fusion model, many spaces that acquisition module is obtained, many temporal resolutions it is distant
Sense data set carries out space-time data fusion, obtains having high time, the remotely-sensed data collection of high spatial resolution simultaneously.
The spatial and temporal scales coalignment of the land evapotranspiration remote sensing appraising described in the present embodiment can be used for performing above-mentioned reality
The spatial and temporal scales matching process of the land evapotranspiration remote sensing appraising described in example is applied, its principle is similar with technique effect, herein no longer
Repeat.
Above example is merely to illustrate technical scheme, rather than a limitation；Although with reference to the foregoing embodiments
The present invention has been described in detail, it will be understood by those within the art that：It still can be to aforementioned each enforcement
Technical scheme described in example is modified, or carries out equivalent to which part technical characteristic；And these are changed or replace
Change, do not make the spirit and scope of the essence disengaging various embodiments of the present invention technical scheme of appropriate technical solution.
Claims (10)
1. a kind of spatial and temporal scales matching process of land evapotranspiration remote sensing appraising, it is characterised in that comprise the steps：
S1, the remotely-sensed data collection that many spaces, many temporal resolutions are obtained using sensor；
S2, the remotely-sensed data collection obtained according to step S1 carry out space-time data fusion, obtain having high time and high spatial simultaneously
The remotely-sensed data collection of resolution；
S3, the data set that step S1 and S2 are generated is carried out into ET parametric inversions using preset data model, obtain ET parameter sets, obtained
To different time and the Remote sensing parameters collection of spatial resolution；
The Remote sensing parameters of S4, the different time obtained according to step S3 and spatial resolution, carry out respectively space scale matching and sentence
Disconnected and time scale matching judgment；
If S5, space scale matching judgment result and time scale matching judgment result meet the required precision of ET estimations, will
The data set of the remotely-sensed data as matching of ET estimation precision requirements is met in step S4.
2. method according to claim 1, it is characterised in that the different time obtained according to step S3 and spatial resolution
Remote sensing parameters, carry out space scale matching judgment, specifically include：
S41, the land-use map for obtaining the generation of different resolution remote sensing image, by land-use map vector quantization, by identical resolution
The land-use map of the vector quantization that rate remote sensing images are obtained is superimposed with the remote sensing land-use map of grid, calculates space shatter value and refers to
Number：SP=Ni/Si；
In formula, SP is space shatter value index；NiFor the Land_use change patch number generated under pre-set spatial resolution；SiFor above-mentioned
The area of Land_use change；
S42, calculating underlying surface heterogeneity index：Sm≥1；
In formula, A be underlying surface heterogeneity index, AaFor selected remotely-sensed data pixel sum；M is mixed pixel decision content, when
For mixed pixel when, be worth for 1, be otherwise 0；SmIt is the ground species included in mixed pixel for correspondence pixel coefficient of similarity
The species number of type；
S43, as SP >=a, while during A≤b, being judged as space scale matched；
When only meeting SP >=a, or when only meeting A≤b, it is judged as that space scale moderate is matched；
When the two is unsatisfactory for, it is judged as that space scale is mismatched；
In formula, a and b is respectively corresponding predetermined threshold value.
3. method according to claim 2, it is characterised in that the different time obtained according to step S3 and spatial resolution
Remote sensing parameters, carry out time scale matching judgment, specifically include：
S41 ', acquisition are in t-t0The period of duration curve data of interval difference evapotranspiration parameter, and judge the curve data and satellite
Whether the dependency between the curve of actual measurement meets：R2>=c, RMSE≤d, P≤0.01；In formula, c, d are default threshold value；t-t0
For time interval；R is correlation coefficient；RMSE is root-mean-square error；P values are significance level discriminant value；
S42 ', it is segmented according to vegetation growing period or Crop growing stage, on the basis of here segmentation, when calculating different data
Between retrievable parametric stability index VP between sectionz：VPz=VARtb-tc(Pztb,Pztb+1,Pztb+2,......,Pztc), and sentence
Disconnected VPzWhether value levels off to 0；Wherein, PzFor parameter value, tb ... tc are different time point values, VARtb-tcFor time zone
Between tb to tc variance computation model, VPzFor the variance yields of tb to tc；
S43 ' if, have one to meet condition in S41 ' and S42 ', be judged as that time scale is matched, the precision for meeting ET estimations will
Ask.
4. method according to claim 3, it is characterised in that if space scale and time scale matching judgment result for
Match somebody with somebody, then using the Remote sensing parameters of the different time of step S3 acquisition and spatial resolution as the data set for meeting ET required precisions,
As a result it is matching, specifically includes：
If space scale matching result is matched, time scale matching result is matching, then it is assumed that meet ET estimation precisions
Require, the data set of the remotely-sensed data of ET estimation precision requirements as matching now will be met in step S4；
If space scale matching result is moderate matching, time scale matching result is matching, then by spatial match and time
The parameter set of matching is input in ET models, if the ET value relations for obtaining the ET estimation results and actual observation simulated meet bar
Part：R2>=e, RMSE≤f, during P≤g, then it is assumed that meet ET estimation precision requirements, now will meet ET estimation precisions in step S4
Data set of the remotely-sensed data of requirement as matching；In formula, e, f, g are corresponding predetermined threshold value.
5. the method according to any one of Claims 1 to 4, it is characterised in that according to the remotely-sensed data collection that step S1 is obtained
Space-time data fusion is carried out, the remotely-sensed data collection with high time and high spatial resolution is obtained, is specifically included：
Using S-G Filtering Models and temporal-spatial fusion model, by many spaces in step S1, the remotely-sensed data collection of many temporal resolutions
Space-time data fusion is carried out, obtains that there is high time, the remotely-sensed data collection of high spatial resolution simultaneously.
6. a kind of spatial and temporal scales coalignment of land evapotranspiration remote sensing appraising, it is characterised in that include：
Acquisition module, for obtaining the remotely-sensed data of many spaces, many temporal resolutions using sensor；
Fusion Module, for processing the data set that model generates the data acquisition module using preset data, when carrying out high
Between and high spatial data fusion, obtain the remotely-sensed data collection with high time and high spatial resolution；
Generation module, the data set for being obtained according to the acquisition module and Fusion Module carries out ET parametric inversions, obtains ET
Parameter set, obtains the Remote sensing parameters collection of different time and spatial resolution；
Matching judgment module, for the different time that generated according to the generation module and the Remote sensing parameters of spatial resolution, point
Space scale matching judgment and time scale matching judgment are not carried out；
Matching result processing module, for determining space scale matching judgment result and time scale in the matching judgment module
When matching judgment result is satisfied by the required precision of ET estimations, the remotely-sensed data of ET estimation precision requirements will be met as matching
Data set.
7. device according to claim 6, it is characterised in that the matching judgment module, specifically for：
S41, the land-use map for obtaining the generation of different resolution remote sensing image, by land-use map vector quantization, by vector quantization
Land-use map is superimposed with the Remote sensing parameters of the grid of equal resolution, calculates space shatter value index：SP=Ni/Si；
In formula, SP is space shatter value index；NiFor the Land_use change patch number generated under pre-set spatial resolution；SiFor above-mentioned
The area of land-use map；
S42, calculating underlying surface heterogeneity index：Sm≥1；
In formula, A be underlying surface heterogeneity index, AaFor selected remotely-sensed data pixel sum；M is mixed pixel decision content, when
For mixed pixel when, value is 1, is otherwise 0；SmIt is the atural object included in mixed pixel for correspondence pixel coefficient of similarity
The species number of type；
S43, as SP >=a, while during A≤b, being judged as space scale matched；
When only meeting SP >=a, or when only meeting A≤b, it is judged as that space scale is matched；
When the two is unsatisfactory for, it is judged as that space scale is mismatched；
In formula, a and b is respectively corresponding predetermined threshold value.
8. device according to claim 7, it is characterised in that the matching judgment module, specifically for：
S41 ', acquisition are in t-t0The period of duration curve data of interval difference evapotranspiration parameter, and judge the curve data and satellite
Whether the dependency between the curve of actual measurement meets：R2>=c, RMSE≤d, P≤0.01；In formula, c, d are default threshold value；t-t0
For time interval；
S42 ', it is segmented according to vegetation growing period or Crop growing stage, on the basis of here segmentation, when calculating different data
Between retrievable parametric stability index VP between sectionz：VPz=VARtb-tc(Pztb,Pztb+1,Pztb+2,......,Pztc), and sentence
Disconnected VPzWhether value levels off to 0；Wherein, PzFor parameter value, tb ... tc are different time point values, VARtb-tcFor variance meter
Calculate function, VPzFor the variance yields of tb to tc；
S43 ' if, have one to meet condition in S41 ' and S42 ', be judged as that time scale is matched.
9. device according to claim 8, it is characterised in that the matching result processing module, specifically for：
Determine that space scale matching result is matched in the matching judgment module, time scale matching result is matching
When, it is believed that meet ET estimation precision requirements, will now meet the data set of the remotely-sensed data of ET estimation precision requirements as matching；
Determine space scale matching result in the matching judgment module to match for moderate, time scale matching result is matching
When, it is input to the relevant parameter model used during spatial match and time match and illegally in ET models, and
The ET value relations for obtaining the ET estimation results and actual observation simulated meet condition：R2>=e, RMSE≤f, during P≤g, it is believed that full
Sufficient ET estimation precisions requirement, will now meet the data set of the remotely-sensed data of ET estimation precision requirements as matching；In formula, e, f,
G is corresponding predetermined threshold value.
10. the device according to any one of claim 6～9, it is characterised in that the Fusion Module, specifically for：
Using S-G Filtering Models and temporal-spatial fusion model, many spaces, the remote sensing number of many temporal resolutions that acquisition module is obtained
Space-time data fusion is carried out according to collection, obtains that there is high time, the remotely-sensed data collection of high spatial resolution simultaneously.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611047244.0A CN106599548B (en) | 2016-11-23 | 2016-11-23 | The spatial and temporal scales matching process and device of land evapotranspiration remote sensing appraising |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611047244.0A CN106599548B (en) | 2016-11-23 | 2016-11-23 | The spatial and temporal scales matching process and device of land evapotranspiration remote sensing appraising |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106599548A true CN106599548A (en) | 2017-04-26 |
| CN106599548B CN106599548B (en) | 2019-05-24 |
Family
ID=58592096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611047244.0A Active CN106599548B (en) | 2016-11-23 | 2016-11-23 | The spatial and temporal scales matching process and device of land evapotranspiration remote sensing appraising |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106599548B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112416915A (en) * | 2020-11-11 | 2021-02-26 | 北京大学 | Data correction method, device, storage medium and equipment |
| CN115964958A (en) * | 2023-02-15 | 2023-04-14 | 中国科学院、水利部成都山地灾害与环境研究所 | Matching method of flux observation footprint and remote sensing pixel spatial scale |
| CN116306265A (en) * | 2023-02-24 | 2023-06-23 | 武汉大学 | Global earth surface evapotranspiration estimation method based on water-carbon coupling method |
| CN116776651A (en) * | 2023-08-23 | 2023-09-19 | 中国科学院空天信息创新研究院 | Method and device for measuring and calculating surface evapotranspiration, electronic equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102314677A (en) * | 2010-07-07 | 2012-01-11 | 中国科学院地理科学与资源研究所 | Classification-based high time definition and high space definition remote sensing data quantitative fusing method |
| CN103177431A (en) * | 2012-12-26 | 2013-06-26 | 中国科学院遥感与数字地球研究所 | Method of spatial-temporal fusion for multi-source remote sensing data |
-
2016
- 2016-11-23 CN CN201611047244.0A patent/CN106599548B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102314677A (en) * | 2010-07-07 | 2012-01-11 | 中国科学院地理科学与资源研究所 | Classification-based high time definition and high space definition remote sensing data quantitative fusing method |
| CN103177431A (en) * | 2012-12-26 | 2013-06-26 | 中国科学院遥感与数字地球研究所 | Method of spatial-temporal fusion for multi-source remote sensing data |
Non-Patent Citations (4)
| Title |
|---|
| 位贺杰等: "基于MOD16数据的渭河流域地表实际蒸散发时空特征", 《中国沙漠》 * |
| 夏浩铭等: "遥感反演蒸散发时间尺度拓展方法研究进展", 《农业工程学报》 * |
| 左其亭等: "《现代水文学》", 30 June 2006 * |
| 张荣华等: "区域蒸散发遥感估算方法及验证综述", 《地球科学进展》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112416915A (en) * | 2020-11-11 | 2021-02-26 | 北京大学 | Data correction method, device, storage medium and equipment |
| CN115964958A (en) * | 2023-02-15 | 2023-04-14 | 中国科学院、水利部成都山地灾害与环境研究所 | Matching method of flux observation footprint and remote sensing pixel spatial scale |
| CN115964958B (en) * | 2023-02-15 | 2023-06-30 | 中国科学院、水利部成都山地灾害与环境研究所 | Matching method of flux observation footprint and remote sensing pixel space scale |
| CN116306265A (en) * | 2023-02-24 | 2023-06-23 | 武汉大学 | Global earth surface evapotranspiration estimation method based on water-carbon coupling method |
| CN116306265B (en) * | 2023-02-24 | 2024-03-08 | 武汉大学 | Global earth surface evapotranspiration estimation method based on water-carbon coupling method |
| CN116776651A (en) * | 2023-08-23 | 2023-09-19 | 中国科学院空天信息创新研究院 | Method and device for measuring and calculating surface evapotranspiration, electronic equipment and storage medium |
| CN116776651B (en) * | 2023-08-23 | 2023-11-14 | 中国科学院空天信息创新研究院 | Method and device for measuring and calculating surface evapotranspiration, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106599548B (en) | 2019-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Peleg et al. | An advanced stochastic weather generator for simulating 2‐D high‐resolution climate variables | |
| Xie et al. | Data assimilation for distributed hydrological catchment modeling via ensemble Kalman filter | |
| Bihlo | A generative adversarial network approach to (ensemble) weather prediction | |
| Montandon et al. | The impact of soil reflectance on the quantification of the green vegetation fraction from NDVI | |
| Revilla-Romero et al. | Integrating remotely sensed surface water extent into continental scale hydrology | |
| Panagos et al. | Monthly soil erosion monitoring based on remotely sensed biophysical parameters: a case study in Strymonas river basin towards a functional pan-European service | |
| Moreno et al. | Mapping daily global solar irradiation over Spain: a comparative study of selected approaches | |
| CN109344865A (en) | A kind of data fusion method of multi-data source | |
| CN106599548B (en) | The spatial and temporal scales matching process and device of land evapotranspiration remote sensing appraising | |
| Walker et al. | Hydrologic data assimilation | |
| Elgamal et al. | Impact analysis of satellite rainfall products on flow simulations in the Magdalena River Basin, Colombia | |
| Cohen et al. | A methodology for calculating the spatial distribution of the area‐slope equation and the hypsometric integral within a catchment | |
| Frances et al. | Topsoil thickness prediction at the catchment scale by integration of invasive sampling, surface geophysics, remote sensing and statistical modeling | |
| CN102034337A (en) | System and method for prairie snow disaster remote sensing monitoring and disaster situation evaluation | |
| Temme et al. | Multi‐process Late Quaternary landscape evolution modelling reveals lags in climate response over small spatial scales | |
| Crow et al. | Optimal averaging of soil moisture predictions from ensemble land surface model simulations | |
| Deschamps‐Berger et al. | Improving the spatial distribution of snow cover simulations by assimilation of satellite stereoscopic imagery | |
| Aduah et al. | Assessing suitability of the ACRU hydrological model in a rainforest catchment in Ghana, West Africa | |
| Paul et al. | Forcing a distributed glacier mass balance model with the regional climate model REMO. Part II: Downscaling strategy and results for two Swiss glaciers | |
| Casado-Rodríguez et al. | Hydrograph separation for tackling equifinality in conceptual hydrological models | |
| Yildiz et al. | Elucidating vegetation controls on the hydroclimatology of a mid-latitude basin | |
| Houser et al. | Land surface data assimilation | |
| CN107957982A (en) | Secondary Geological Hazards liability fast evaluation method and system after shake | |
| Zus et al. | Development and optimization of the IPM MM5 GPS slant path 4DVAR system | |
| Herzfeld et al. | Sea ice deformation in Fram Strait—Comparison of CICE simulations with analysis and classification of airborne remote-sensing data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |