CN104143031A - Vegetation index time series data reconstruction method based on wavelet multi-scale decomposition - Google Patents

Vegetation index time series data reconstruction method based on wavelet multi-scale decomposition Download PDF

Info

Publication number
CN104143031A
CN104143031A CN201310164200.6A CN201310164200A CN104143031A CN 104143031 A CN104143031 A CN 104143031A CN 201310164200 A CN201310164200 A CN 201310164200A CN 104143031 A CN104143031 A CN 104143031A
Authority
CN
China
Prior art keywords
time series
series data
vegetation index
model
vegetation
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
Application number
CN201310164200.6A
Other languages
Chinese (zh)
Other versions
CN104143031B (en
Inventor
邱炳文
钟鸣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN201310164200.6A priority Critical patent/CN104143031B/en
Publication of CN104143031A publication Critical patent/CN104143031A/en
Application granted granted Critical
Publication of CN104143031B publication Critical patent/CN104143031B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a vegetation index time series data reconstruction method based on wavelet multi-scale decomposition. The method is characterized by comprising the following steps: performing wavelet transformation according to the vegetation index time series data to respectively sequentially decompose the original time series data of vegetation indexes and climatic factors into corresponding half moth, month, double month, season, half a year, and time series data based on interannual scales; further selecting the climatic factor time series data under the corresponding scale according to the features of the vegetation index time series data in each scale; selecting a proper model; reconstructing time series data in different scales; finally combining the vegetation index time series data in all scales to realize the reconstruction of the original vegetation index time series data. The method has the characteristics of being high in precision and wide in applicable scope.

Description

A kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition
Technical field
The present invention relates to techniques of teime series analysis field, be specifically related to a kind of vegetation index (Vegetation index, VI) time series data reconstructing method based on Multiscale Wavelet Decomposition.
Background technology
Remote sensing vegetation index time series data is widely used in forest, crops vegetation Monitoring on Dynamic Change.In remote sensing data acquiring and image processing process, due to the interference of the various factorss such as observation angle and cloud, cause generated vegetation index time series data quality not ideal enough, therefore need further original vegetation index time series data to be carried out to denoising and reconstruct.
The method of the denoising of vegetation index time series data and reconstruct has a lot, as optimum gradient index extraction method (Best slope extraction, BISE), Fourier analysis method, polynary least square method, Geostatistics, non-linear Gaussian function method and S-G (Savizky-Golay) filter method etc.These methods all have certain rationality and practical reference value, but that its weak point is that vegetation index time series data generally all has is non-stationary, are not suitable for adopting stationarity method, therefore inevitably have some limitations.
Wavelet transformation, as a kind of mathematical tool of multiscale analysis, can be decomposed into the time series data on a series of different scales by original time series data effectively.The present invention carries out wavelet transformation to vegetation index time series data, by differentiating time series data on the different scale generating after wavelet decomposition, whether meet stationarity, and according to the changing character of himself, choose suitable method, respectively the time series data on different scale is chosen to suitable method and rebuild, will provide a kind of new thinking and method for the reconstruct of vegetation index time series data.
Summary of the invention
In view of this, the object of the present invention is to provide a kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition, for achieving the above object, technical scheme of the present invention is: a kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition, it is characterized in that: the method is based on vegetation index time series data, utilize wavelet transformation, the original time series data of vegetation index and climatic factor is decomposed into respectively to corresponding first quarter moon successively, month, bimonthly, season, time series data on half a year and year border yardstick, further according to the vegetation index time series data feature on each yardstick, select the climatic factor time series data on corresponding yardstick, choose suitable model, the time series data carrying out on different scale is rebuild, vegetation index time series data on last comprehensive all yardsticks is realized the reconstruct of original vegetation index time series data.
Method of the present invention is taken the non-stationary feature of vegetation index into account, and generalized time and space two aspect factors, according to time series data feature on the different scale generating after wavelet decomposition, adopt respectively suitable method, carries out the reconstruct of vegetation index time series data.The present invention has the features such as precision is high, applied widely.
Accompanying drawing explanation
Fig. 1 is the realization flow figure of the embodiment of the present invention.
Embodiment
Below in conjunction with drawings and the specific embodiments, the present invention is described in detail.
The present invention is based on the vegetation index time series data reconstructing method of Multiscale Wavelet Decomposition, as shown in Figure 1, utilize wavelet transformation, vegetation index original signal is decomposed into the High-frequency and low-frequency composition on a plurality of yardsticks such as corresponding first quarter moon respectively, the moon, bimonthly, season, half a year and Nian Ji, further according to the vegetation index time series data Changing Pattern on each yardstick, choose suitable model and carry out the reconstruction of vegetation index time series data, last comprehensive different scale time series data is realized the reconstruct of vegetation index time series data.
Concrete, the vegetation index time series data reconstructing method of the present embodiment based on Multiscale Wavelet Decomposition, further comprising the steps:
Step S1: set up vegetation index time series data collection.
Based on 16 days synthetic MODIS EVI (enhancement mode vegetation index) data (MOD09Q1), set up study area 2001-2011 phase vegetation index time series data collection for many years.
Step S2: set up climatic factor time series data collection.
Based on Continuous Observation climatic factor time series data collection every day, generate corresponding with the vegetation index data 16 days synthetic temperature in 2001-2011 study area, precipitation time series data collection.
Step S3: respectively vegetation index and climatic factor time series data are carried out to wavelet transformation, obtain vegetation index and climatic factor time series data collection on different scale.
Based on the female small echo of Meyer, the study area original time series data of 2001-2011 vegetation index is carried out to wavelet transformation 1 time, generate respectively low-frequency component A 1and radio-frequency component D (Approximation) 1(Detail), then to low-frequency component A 1proceed wavelet transformation, further obtain low-frequency component A 2with radio-frequency component D 2, by that analogy, until obtain the low-frequency component A after the 5th wavelet decomposition 5with radio-frequency component D 5.According to the female wavelet center frequency of Meyer and original vegetation index time series data interval, can calculate and obtain D 1, D 2, D 3, D 4, D 5, A 5the vegetation index time series data on corresponding first quarter moon, moon, bimonthly, season, half a year and year border yardstick roughly respectively.Be specially: D 1(<24 days), D 2(24-48 days), D 3(48-95 days), D 4(95-190 days), D 5(190-380 days), A 5(>380 days).
Wavelet transformation can be expressed as:
(1)
A wherein, b is respectively frequency field and time domain parameter, for original signal, for wavelet coefficient.
Because the little wave energy of Meyer is guaranteed that time series data is level and smooth and have the advantage such as to be simple and easy to use, therefore adopt Meyer small echo to carry out multi-scale wavelet transformation, be expressed as:
(2)
According to said method, climatic factor time series data is carried out to wavelet transformation, obtain the climatic factor time series data D on corresponding first quarter moon respectively, moon, bimonthly, season, half a year and year border yardstick 1, D 2, D 3, D 4, D 5, A 5.
Step S4: choose suitable model form, rebuild respectively D 1(first quarter moon), D 2(moon), D 3(bimonthly), D 4(season), D 5(half a year), A 5vegetation index time series data on different scales such as (year borders).
Step S4.1 adopts autoregression method to D 1(first quarter moon), D 2(moon) time series data is rebuild.Adopt akaike information criterion (Akaike Information Criterion) method to carry out determining rank to model, set up D 1, D 2autoregressive model.Model form is as follows:
(3)
Wherein, a 1, a 2a nfor coefficient, y (t-1), y (t-2) ... y (t-n) is respectively 1,2 ... n rank autoregression item, for residual error.
Step S4.2 adopts autoregressive model with controlled quentity controlled variable to D 3(bimonthly) and D 5(half a year), time series data was rebuild.First adopt akaike information criterion (Akaike Information Criterion) method to determine Autoregressive, then set up time series data y(D 3or D 5) with a plurality of factor time series data (D such as temperature (T) on autoregression project and corresponding yardstick and precipitation (R) 3or D 5) quantitative model.Model form is:
(4)
Wherein, a 1, a 2a n, b, c is coefficient, y (t-1), y (t-2) ... y (t-n) is respectively 1,2 ... n rank autoregression item, the D that T is temperature index 3or D 5time series data, R is the D of precipitation index 3or D 5time series data, for residual error.
Step S4.3 adopts linear regression model (LRM) method to D 4(season), time series data was reconstructed.With the time series data (D on the corresponding yardstick of the climatic factors such as temperature and precipitation 4) as independent variable, set up vegetation index time series data y(D on yardstick in season 4) answer the time series data (D on yardstick with a plurality of factor pairs such as temperature (T) and precipitation (R) 4) quantitative model.Model form is:
(5)
Wherein, a, b is coefficient, C is constant, the D that T is temperature index 4time series data, R is the D of precipitation index 4time series data, for residual error.
Step S4.4 calculates the A of each pixel 5the mean value of time series data, is called vegetation signal criterion surface layer.Spatial distribution characteristic based on vegetation signal criterion surface layer, choose respectively elevation (DEM), the gradient (slope), Land Use Degree (land use intensity, LUI), mean annual temperature (Temperature, T), mean precipitation (rainfall for many years, R) and with the distance (river) in nearest river, some factors of influence such as distance (resident) with nearest residential area, first to elevation, the gradient and the accessibility factor are (apart from river, the distance in residential area) carry out log-transformation, then all independents variable are done to normalized, finally set up the linear regression model (LRM) that obtains vegetation signal criterion surface layer (B).Its model form is as follows:
(6)
Wherein, , for coefficient, C is constant term, DEM, Slope, River, Resident, LUI, T and R represent respectively through pretreated absolute elevation, the gradient, with the distance in river recently, with distance, Land Use Degree and the temperature in residential area and the A of precipitation index recently 5the average of time series data, for residual error.
Obtaining on the basis of vegetation signal criterion surface layer (Base), by the A obtaining after wavelet decomposition by pixel 5time series data subtracts each other with corresponding vegetation signal criterion surface layer, obtains A 5the difference of level and reference field, is called the A of this pixel 5trend time series data.First by ADF (Augmented Dickey-Fuller test) check, determine A 5whether trend time series data meets time stationarity.If do not meet time stationarity, by difference method, ask calculation increment sequence, until increment sequence reaches steadily, thereby obtain A stably 5trend time series data.To A stably 5trend time series data, adopts AIC method to carry out determining rank to model, sets up A 5the ARMA model of trend time series data.Model form is as follows:
(7)
Wherein for autoregressive coefficient, for moving average coefficient, for white noise sequence.
Step 5: by the D after rebuilding 1, D 2, D 3, D 4, D 5, A 5signal plus, thus realize the reconstruct of vegetation index time series data.
D after comprehensive redevelopment 1, D 2, D 3, D 4, D 5, A 5signal, obtains the vegetation index time series data after reconstruct s ':
Wherein for the time series data after rebuilding on corresponding yardstick.
The foregoing is only preferred embodiment of the present invention, all equalizations of doing according to the present patent application the scope of the claims change and modify, and all should belong to covering scope of the present invention.

Claims (9)

1. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition, it is characterized in that: the method is based on vegetation index time series data, utilize wavelet transformation, the original time series data of vegetation index and climatic factor is decomposed into respectively to corresponding first quarter moon successively, month, bimonthly, season, time series data on half a year and year border yardstick, further according to the vegetation index time series data feature on each yardstick, select the climatic factor time series data on corresponding yardstick, choose suitable model, the time series data carrying out on different scale is rebuild, vegetation index time series data on last comprehensive all yardsticks is realized the reconstruct of original vegetation index time series data.
2. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 1, it is characterized in that: describedly based on vegetation index time series data, comprise: first, according to N days synthetic vegetation index data, generate phase vegetation index sequential image for many years; Then generate and this temperature of phase for many years, precipitation time series data that phase vegetation index time series data is corresponding for many years, N is greater than 10 natural number.
3. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 2, is characterized in that: described N is 16.
4. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 1, is characterized in that: also comprise by the time series data on described first quarter moon, the moon, bimonthly, season, half a year and year border yardstick respectively correspondence be defined as D 1, D 2, D 3, D 4,, D 5, A 5; Describedly choose suitable model and comprise: D 1and D 2adopt autoregressive model, D 3and D 5employing is with the autoregressive model of controlled quentity controlled variable, D 4adopt linear regression model (LRM), A 5first adopt trend surface model to carry out overall fit, then to residual error, partly adopt difference to move autoregressive model.
5. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 4, is characterized in that: described D 1and D 2adopt the implementation of autoregressive model to be: to adopt akaike information criterion (Akaike Information Criterion is called for short AIC) to carry out determining rank to model, set up D 1, D 2autoregressive model; Model form is as follows:
Wherein, a 1, a 2a nfor coefficient, y (t-1), y (t-2) ... y (t-n) is respectively 1,2 ... n rank autoregression item, for residual error.
6. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 4, is characterized in that: described D 3and D 5employing with the implementation of the autoregressive model of controlled quentity controlled variable is: first adopt akaike information criterion (Akaike Information Criterion, being called for short AIC) method determines Autoregressive, then sets up the quantitative model of time series data y and the time series data of temperature T on autoregression item and corresponding yardstick and precipitation R; Model form is:
Wherein, a 1, a 2a n, b, c is coefficient, y (t-1), y (t-2) ... y (t-n) is respectively 1,2 ... n rank autoregression item, the D that T is temperature index 3or D 5time series data, R is the D of precipitation index 3or D 5time series data, for residual error.
7. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 4, is characterized in that: described D 4adopt the implementation of linear regression model (LRM) to be: with the D on temperature and the corresponding yardstick of rainfall factor 4time series data, as independent variable, is set up vegetation index time series data y(D on yardstick in season 4) answer the D on yardstick with temperature T and precipitation R factor pair 4the quantitative model of time series data; This model form is:
Wherein, a, b is coefficient, C is constant, the D that T is temperature index 4time series data, R is the D of precipitation index 4time series data, for residual error.
8. the vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition according to claim 4, is characterized in that: the A that calculates each pixel 5the mean value of time series data, is called vegetation signal criterion surface layer, spatial distribution characteristic based on vegetation signal criterion face, choose respectively elevation (DEM), the gradient (slope), Land Use Degree (land use intensity, LUI), mean annual temperature (temperature, T), mean precipitation (rainfall for many years, R) and with the distance (river) in nearest river, some factors of influence such as distance (resident) with nearest residential area, first to elevation, the gradient and the accessibility factor are (apart from river, the distance in residential area) carry out log-transformation, then all independents variable are done to normalized, finally set up the linear regression model (LRM) that obtains vegetation signal criterion surface layer (B), its model form is as follows:
Wherein, , for coefficient, C is constant term, DEM, Slope, River, Resident, LUI, T and R represent respectively through pretreated absolute elevation, the gradient, with the distance in river recently, with distance, Land Use Degree and the temperature in residential area and the A of precipitation index recently 5the average of time series data, for residual error;
Obtaining on the basis of vegetation signal criterion surface layer (Base), by the A obtaining after wavelet decomposition by pixel 5time series data subtracts each other with corresponding vegetation signal criterion surface layer, obtains A 5the difference of level and reference field, is called the A of this pixel 5trend time series data; First by ADF (Augmented Dickey-Fuller test) check, determine A 5whether trend time series data meets time stationarity;
If do not meet time stationarity, by difference method, ask calculation increment sequence, until increment sequence reaches steadily, thereby obtain A stably 5trend time series data; To A stably 5trend time series data, adopts akaike information criterion (Akaike Information Criterion is called for short AIC) method to carry out determining rank to model, sets up A 5the ARMA model of trend time series data; Model form is as follows:
Wherein for autoregressive coefficient, for moving average coefficient, for white noise sequence.
9. according to claim 1 a kind of based on the multi-level vegetation index time series data reconstructing method decomposing of small echo, it is characterized in that: the method is applied in time series data denoising, reconstruct and prediction field.
CN201310164200.6A 2013-05-07 2013-05-07 A kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition Active CN104143031B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310164200.6A CN104143031B (en) 2013-05-07 2013-05-07 A kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310164200.6A CN104143031B (en) 2013-05-07 2013-05-07 A kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition

Publications (2)

Publication Number Publication Date
CN104143031A true CN104143031A (en) 2014-11-12
CN104143031B CN104143031B (en) 2017-08-11

Family

ID=51852203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310164200.6A Active CN104143031B (en) 2013-05-07 2013-05-07 A kind of vegetation index time series data reconstructing method based on Multiscale Wavelet Decomposition

Country Status (1)

Country Link
CN (1) CN104143031B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104809357A (en) * 2015-05-21 2015-07-29 北京理工大学 Analysis method and analyzer for multi-scale contact behaviors of three-dimensional rough surfaces
CN105718936A (en) * 2016-02-02 2016-06-29 福州大学 Forest dynamic change mode automatic extraction method
CN106708782A (en) * 2016-12-22 2017-05-24 贾翔 Regional pest detection diagnosis judging method based on wavelet analysis
CN107918660A (en) * 2017-11-21 2018-04-17 清华大学 High frequency time series data processing method based on multistage differential networks
CN110852585A (en) * 2019-10-30 2020-02-28 北京师范大学 Calculation method and device for vegetation growth stability
CN111753738A (en) * 2020-06-24 2020-10-09 北京师范大学 Vegetation annual change monitoring method and system based on wavelet analysis

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101604356A (en) * 2009-06-24 2009-12-16 南京大学 A kind of method for building up of uncertain mid-and-long term hydrologic forecast model
CN101650422A (en) * 2009-09-27 2010-02-17 北京师范大学 Denoising method of remote sensing vegetation index time series data
CN101718865A (en) * 2009-12-07 2010-06-02 中国科学院新疆生态与地理研究所 Method for delaying vegetation response based on remote sensing
CN102176242A (en) * 2011-02-01 2011-09-07 环境保护部卫星环境应用中心 Method for removing cloud noise effects in normalized difference vegetation index (NDVI) time sequence image
CN102323970A (en) * 2011-05-19 2012-01-18 河海大学 Hydrological time series prediction method based on multiple-factor wavelet neural network model
CN102495937A (en) * 2011-10-18 2012-06-13 南京信息工程大学 Prediction method based on time sequence

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101604356A (en) * 2009-06-24 2009-12-16 南京大学 A kind of method for building up of uncertain mid-and-long term hydrologic forecast model
CN101650422A (en) * 2009-09-27 2010-02-17 北京师范大学 Denoising method of remote sensing vegetation index time series data
CN101718865A (en) * 2009-12-07 2010-06-02 中国科学院新疆生态与地理研究所 Method for delaying vegetation response based on remote sensing
CN102176242A (en) * 2011-02-01 2011-09-07 环境保护部卫星环境应用中心 Method for removing cloud noise effects in normalized difference vegetation index (NDVI) time sequence image
CN102323970A (en) * 2011-05-19 2012-01-18 河海大学 Hydrological time series prediction method based on multiple-factor wavelet neural network model
CN102495937A (en) * 2011-10-18 2012-06-13 南京信息工程大学 Prediction method based on time sequence

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
宋涛: "基于小波变换的植被覆盖变化研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
李晓兵等: "中国典型植被类型NDVI动态变化与气温、降水变化的敏感性分析", 《植物生态学报》 *
王军邦等: "内蒙古中部MODIS植被动态监测分析", 《地球信息科学学报》 *
聂中青等: "近50a葫芦河流域气候变化特征", 《兰州大学学报(自然科学版)》 *
马明国等: "基于遥感的植被年际变化及其与气候关系研究进展", 《遥感学报》 *
黄玫等: "土地覆被的气候预测模型", 《地理学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104809357A (en) * 2015-05-21 2015-07-29 北京理工大学 Analysis method and analyzer for multi-scale contact behaviors of three-dimensional rough surfaces
CN104809357B (en) * 2015-05-21 2017-10-03 北京理工大学 The analysis method and device of a kind of multiple dimensioned touching act of three-dimensional rough surface
CN105718936A (en) * 2016-02-02 2016-06-29 福州大学 Forest dynamic change mode automatic extraction method
CN106708782A (en) * 2016-12-22 2017-05-24 贾翔 Regional pest detection diagnosis judging method based on wavelet analysis
CN106708782B (en) * 2016-12-22 2019-05-03 贾翔 Region pest and disease damage checkout and diagnosis method of discrimination based on wavelet analysis
CN107918660A (en) * 2017-11-21 2018-04-17 清华大学 High frequency time series data processing method based on multistage differential networks
CN110852585A (en) * 2019-10-30 2020-02-28 北京师范大学 Calculation method and device for vegetation growth stability
CN110852585B (en) * 2019-10-30 2022-05-10 北京师范大学 Calculation method and device for vegetation growth stability
CN111753738A (en) * 2020-06-24 2020-10-09 北京师范大学 Vegetation annual change monitoring method and system based on wavelet analysis

Also Published As

Publication number Publication date
CN104143031B (en) 2017-08-11

Similar Documents

Publication Publication Date Title
CN104143031A (en) Vegetation index time series data reconstruction method based on wavelet multi-scale decomposition
Nourani et al. Wavelet-entropy data pre-processing approach for ANN-based groundwater level modeling
Deka et al. Discrete wavelet neural network approach in significant wave height forecasting for multistep lead time
Altunkaynak et al. Prediction of daily rainfall by a hybrid wavelet-season-neuro technique
CN110361778B (en) Seismic data reconstruction method based on generation countermeasure network
Liu et al. Time-frequency analysis and simulation of the watershed suspended sediment concentration based on the Hilbert-Huang transform (HHT) and artificial neural network (ANN) methods: A case study in the Loess Plateau of China
Nguyen et al. Correcting for systematic biases in GCM simulations in the frequency domain
Hou et al. Multifractal analysis of the drought area in seven large regions of China from 1961 to 2012
CN101303764A (en) Method for self-adaption amalgamation of multi-sensor image based on non-lower sampling profile wave
Adamala Temperature based generalized wavelet-neural network models to estimate evapotranspiration in India
CN110503060B (en) Spectral signal denoising method and system
CN104809357B (en) The analysis method and device of a kind of multiple dimensioned touching act of three-dimensional rough surface
CN101930598A (en) Natural image denoising method based on non-local mean value of shearlet region
Zhang et al. VST-Net: Variance-stabilizing transformation inspired network for Poisson denoising
CN105205788A (en) Denoising method for high-throughput gene sequencing image
CN107977944A (en) One kind estimation NPP remote sensing image data generation methods
Foroumandi et al. Climate change or regional human impacts? Remote sensing tools, artificial neural networks, and wavelet approaches aim to solve the problem
CN109709585A (en) The method for removing coloured noise in GPS coordinate time series
Roushangar et al. Identifying complexity of annual precipitation variation in Iran during 1960–2010 based on information theory and discrete wavelet transform
Haddad et al. Sea level in the Mediterranean Sea: seasonal adjustment and trend extraction within the framework of SSA
CN104268408A (en) Energy consumption data macro-forecast method based on wavelet coefficient ARMA model
Zhang et al. Multi-source information fused generative adversarial network model and data assimilation based history matching for reservoir with complex geologies
Liu et al. A new framework for rainfall downscaling based on EEMD and an improved fractal interpolation algorithm
Wang et al. An image fusion algorithm based on lifting wavelet transform
Evrendilek et al. Regression model-based predictions of diel, diurnal and nocturnal dissolved oxygen dynamics after wavelet denoising of noisy time series

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant