CN104616265A - Time domain reconstruction method of remote sensing sequence data - Google Patents
Time domain reconstruction method of remote sensing sequence data Download PDFInfo
- Publication number
- CN104616265A CN104616265A CN201510074076.3A CN201510074076A CN104616265A CN 104616265 A CN104616265 A CN 104616265A CN 201510074076 A CN201510074076 A CN 201510074076A CN 104616265 A CN104616265 A CN 104616265A
- Authority
- CN
- China
- Prior art keywords
- data
- remote sensing
- time
- time series
- time domain
- 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
Landscapes
- Image Processing (AREA)
Abstract
The invention discloses a time domain reconstruction method of remote sensing sequence data. The time domain reconstruction method comprises the following steps: firstly, calculating the correlation of non-missing public parts between the remote sensing sequence data, and reconstructing reasonable distribution weights for missed parts through the correlation size; secondly, calculating a matching corresponding value of a reference image and an image to be reconstructed through moment matching; and finally, reconstructing a time sequence by a weighting harmonic analyzing method. By the aid of the time domain reconstruction method, a gray level of the reference image horizontally gets closer to the gray level of the image to be reconstructed by a moment matching principle through the correlation reasonable distribution weights between the data; time domain reconstruction of the remote sensing sequence data is carried out by the weighting harmonic analyzing method through abiding by a time dependence relationship, so that the time dependence relationship of the time sequence data is sufficiently utilized and a reconstructed result is real and available; the time domain reconstruction method is easy to realize and has high efficiency and important actual application meanings.
Description
Technical field
The invention belongs to remote sensing image processing technology field, relate to a kind of time domain method for reconstructing of remote sensing sequence data, be specifically related to the remote sensing time series data time domain method for reconstructing of a kind of associating weighting frequency analysis and match by moment.
Background technology
Remote sensing long-term sequence data are widely used in global environmental change research, but pollute and many-sided impact such as uncertainty of air owing to being subject to cloud, cause remote sensing image and may there is a large amount of noise or loss of data, the normal use of remote sensing image is hindered.Therefore, develop rational method for reconstructing, for remote sensing time series data collection carries out denoising, repair process after acquisition, the quality and the stability that improve data are far reaching problems.
At present, the reparation for remote sensing image mainly divides based on single width image with based on several images two kinds.Based on the reparation of single width image, mainly rely on the information of self from mending the angle of painting, its main method has the method for interpolation, Histogram Matching, partial differential equation and total variation method etc.When repairing single width image, if disappearance area is less, can be repaired result preferably, during if there is large area information dropout, the information of single width image self is not enough to support the reconstruction to lack part, is difficult to well be repaired result.Based on the reparation of several images, compensate for the phenomenon of spatial information deficiency, conventional method has the method such as Histogram Matching and regretional analysis, data reconstruction can be carried out by setting up funtcional relationship between different images, if different images spectral characteristic exists very large difference, it is repaired result and often occurs significantly rebuilding vestige.At present, the process for long-term sequence data carries out denoising by the mode of filtering to data, and this little filtering method does not often consider the utilization of spatial information, and the data for consecutive miss are difficult to the reconstructed results obtained.
Summary of the invention
The present invention is directed to the shortcoming of prior art to space-time complementary information underutilization in time series, the remote sensing time series data time domain method for reconstructing of a kind of associating weighting frequency analysis and match by moment is proposed, on the basis utilizing the long-term sequence data time relations of dependence, consider spatial coherence, making full use of data space information in conjunction with match by moment theory, by adopting weighting harmonic analysis method, time domain reconstruction can be carried out to missing data accurately, and counting yield is high, practical.
The technical solution adopted in the present invention is: a kind of time domain method for reconstructing of remote sensing sequence data, is characterized in that, comprise the following steps:
Step 1: remote sensing time series data is carried out geometrical registration, obtains accurate registration image;
Step 2: carry out cloud detection to remote sensing time series data, obtains the cloud mask of each time data, and calculates the related coefficient CC of the public cloud-free area of each time series data;
Step 3: with reference to data public cloudless region grey level's approximate transform to the public cloudless region grey level of data to be reconstructed, generate new time series data;
Step 4: the related coefficient according to step 2 gained distributes weights, rebuilds new time series data.
As preferably, when step 2 carries out cloud detection, be also marked as cloud sector by the exceptional value of threshold decision.
As preferably, the related coefficient CC of the public cloud-free area of each time series data of the calculating described in step 2, its computing formula is:
Wherein x, y are two remotely-sensed datas, and N is the public cloud-free area number of pixels of remotely-sensed data, x
i, y
ii-th pixel in data x, y, 1≤i≤N, m
x, m
ythe average of the corresponding public cloud-free area of remotely-sensed data respectively.
As preferably, the specific implementation process of step 3 is, by moment-matching method with reference to the horizontal approximate transform of data gray value to data gray value level to be reconstructed, calculate a new time proximity sequence data in each data moment to be reconstructed; For moment t to be reconstructed
othe gray average of data is m
to, gray variance is s
to, other moment t
fthe gray average of data is m
tf, gray variance is s
tf, gray-scale value is g
tf, then new time proximity sequence data corresponding grey scale value is:
As preferably, the specific implementation of step 4 comprises following sub-step:
Step 4.1: to reconstruction moment t
odata are recovered;
First be the new time series g of step 3 gained
newin other time data distribute weights according to the related coefficient (CC) of step 2 gained, its weight computing mode is:
Wherein, q is the number without missing data in time series, CC
iit is the related coefficient of i-th data and the public cloud-free area of data to be reconstructed;
Next converts data to the form that sine and cosine is added, and answers gray-scale value g to the new time proximity sequence pair of step 3 gained
newbe weighted and solve, its implementation is:
Wherein, nf is harmonic wave number, a
0for coefficient when harmonic wave number is 0, a
jand b
jfor the coefficient of trigonometric function each several part, f
jfor frequency, t
jfor timing node, y is time proximity sequence data;
Step 4.2: travel through all single picture point time sequences, repeating above-mentioned steps can rebuild remote sensing time series data, the reconstruction mark problem caused with the grey level's difference eliminated between different images.
The remote sensing time series data time domain method for reconstructing of the associating weighting frequency analysis that the present invention proposes and match by moment utilizes the message complementary sense of remote sensing long-term sequence data to carry out time domain reconstruction.Not only take into account the temporal relations of dependence, and by the correlativity reasonable distribution weights between image, convert image greyscale level by moment-matching method, construct new time series and missing values is repaired, the reconstruction mark problem caused with the grey level's difference eliminated between different images.In a word, the method that the present invention proposes effectively can eliminate the impact of the factors such as cloud, improves the q&r of remote sensing time series data, for subsequent applications provides reliable support.Therefore, the method that the present invention proposes has important practical application meaning.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of the embodiment of the present invention.
Embodiment
Understand for the ease of those of ordinary skill in the art and implement the present invention, below in conjunction with drawings and Examples, the present invention is described in further detail, should be appreciated that exemplifying embodiment described herein is only for instruction and explanation of the present invention, is not intended to limit the present invention.
Ask for an interview Fig. 1, the time domain method for reconstructing of a kind of remote sensing sequence data of the embodiment of the present invention, comprises the following steps:
Step 1: first need to carry out geometry correction to remote sensing time series data, obtain accurate registration image, concrete geometry correction is prior art.
Step 2: carry out cloud detection to remote sensing time series data, obtains the cloud mask of each time data, and calculates the related coefficient CC of the public cloud-free area of each time series data.
If remote sensing time series data exists quality status stamp wave band, carry out cloud Covering judgment by this wave band, if there is no this wave band then by other data or cloud detection algorithm carry out cloud cover detection set up cloud mask.Meanwhile, the exceptional value existed in data or invalid value are also labeled as cloud, are rebuild in subsequent treatment.After cloud detection, the related coefficient of sequence data computing time public cloud-free area each other, in order to reasonable distribution weights.Wherein the related coefficient of the public cloud-free area of each time series data is calculated by following formula:
Wherein x, y are two remotely-sensed datas, and N is the public cloud-free area number of pixels of remotely-sensed data, x
i, y
ii-th pixel in data x, y, 1≤i≤N, m
x, m
ythe average of the corresponding public cloud-free area of remotely-sensed data respectively.
Step 3: with reference to data public cloudless region grey level's approximate transform to the public cloudless region grey level of data to be reconstructed, generate new time series data.
It is process by picture point time series processing that this step is rebuild remote sensing time series data.Needed to carry out greyscale transformation in this pixel without missing data before execution is by picture point time series processing, generate a new time series.This process carries out match by moment by single picture point time sequence without missing data place image and this time series data place to be reconstructed image, and by the grey level of the horizontal approximate transform of nothing disappearance data gray to data place to be reconstructed image, its settling mode is:
Wherein, moment t to be reconstructed
othe cloud-free area gray average of data is m
to, gray variance is s
to, other moment t
fdata cloud-free area gray average is m
tf, gray variance is s
tf, gray-scale value is g
tf.
Step 4: distribute weights according to related coefficient, new time series is rebuild.
This step is carried out by picture point time series processing to remote sensing time series data, needs to rebuild missing data in new time series, to reconstruction moment t
odata are recovered, and utilize the method for weighting frequency analysis, convert data to the form that sine and cosine is added; Specific implementation comprises following sub-step:
Step 4.1: to reconstruction moment t
odata are recovered;
First be the new time series g of step 3 gained
newin other time data distribute weights according to the related coefficient (CC) of step 2 gained, its weight computing mode is:
Wherein, q is the number without missing data in time series, CC
iit is the related coefficient of i-th data and the public cloud-free area of data to be reconstructed;
Next converts data to the form that sine and cosine is added, and answers gray-scale value g to the new time proximity sequence pair of step 3 gained
newbe weighted and solve, its implementation is:
Wherein, nf is harmonic wave number, a
0for coefficient when harmonic wave number is 0, a
jand b
jfor the coefficient of trigonometric function each several part, f
jfor frequency, t
jfor timing node, y is time proximity sequence data;
Step 4.2: travel through all single picture point time sequences, repeating above-mentioned steps can rebuild remote sensing time series data, the reconstruction mark problem caused with the grey level's difference eliminated between different images.
Should be understood that, the part that this instructions does not elaborate all belongs to prior art.
Should be understood that; the above-mentioned description for preferred embodiment is comparatively detailed; therefore the restriction to scope of patent protection of the present invention can not be thought; those of ordinary skill in the art is under enlightenment of the present invention; do not departing under the ambit that the claims in the present invention protect; can also make and replacing or distortion, all fall within protection scope of the present invention, request protection domain of the present invention should be as the criterion with claims.
Claims (5)
1. a time domain method for reconstructing for remote sensing sequence data, is characterized in that, comprise the following steps:
Step 1: remote sensing time series data is carried out geometrical registration, obtains accurate registration image;
Step 2: carry out cloud detection to remote sensing time series data, obtains the cloud mask of each time data, and calculates the related coefficient CC of the public cloud-free area of each time series data;
Step 3: with reference to data public cloudless region grey level's approximate transform to the public cloudless region grey level of data to be reconstructed, generate new time series data;
Step 4: the related coefficient according to step 2 gained distributes weights, rebuilds new time series data.
2. the time domain method for reconstructing of remote sensing sequence data according to claim 1, is characterized in that: when step 2 carries out cloud detection, is also marked as cloud sector by the exceptional value of threshold decision.
3. the time domain method for reconstructing of remote sensing sequence data according to claim 1, is characterized in that: the related coefficient CC of the public cloud-free area of each time series data of the calculating described in step 2, and its computing formula is:
Wherein x, y are two remotely-sensed datas, and N is the public cloud-free area number of pixels of remotely-sensed data, x
i, y
ii-th pixel in data x, y, 1≤i≤N, m
x, m
ythe average of the corresponding public cloud-free area of remotely-sensed data respectively.
4. the time domain method for reconstructing of remote sensing sequence data according to claim 3, it is characterized in that: the specific implementation process of step 3 is, by moment-matching method with reference to the horizontal approximate transform of data gray value to data gray value level to be reconstructed, calculate a new time proximity sequence data in each data moment to be reconstructed; For moment t to be reconstructed
othe gray average of data is m
to, gray variance is s
to, other moment t
fthe gray average of data is m
tf, gray variance is s
tf, gray-scale value is g
tf, then new time proximity sequence data corresponding grey scale value is:
5. the time domain method for reconstructing of remote sensing sequence data according to claim 4, is characterized in that, the specific implementation of step 4 comprises following sub-step:
Step 4.1: to reconstruction moment t
odata are recovered;
First be the new time series g of step 3 gained
newin other time data distribute weights according to the related coefficient (CC) of step 2 gained, its weight computing mode is:
Wherein, q is the number without missing data in time series, CC
iit is the related coefficient of i-th data and the public cloud-free area of data to be reconstructed;
Next converts data to the form that sine and cosine is added, and answers gray-scale value g to the new time proximity sequence pair of step 3 gained
newbe weighted and solve, its implementation is:
Wherein, nf is harmonic wave number, a
0for coefficient when harmonic wave number is 0, a
jand b
jfor the coefficient of trigonometric function each several part, f
jfor frequency, t
jfor timing node, y is time proximity sequence data;
Step 4.2: travel through all single picture point time sequences, repeating above-mentioned steps can rebuild remote sensing time series data, the reconstruction mark problem caused with the grey level's difference eliminated between different images.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510074076.3A CN104616265B (en) | 2015-02-12 | 2015-02-12 | A kind of time domain method for reconstructing of remote sensing sequence data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510074076.3A CN104616265B (en) | 2015-02-12 | 2015-02-12 | A kind of time domain method for reconstructing of remote sensing sequence data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104616265A true CN104616265A (en) | 2015-05-13 |
CN104616265B CN104616265B (en) | 2018-02-09 |
Family
ID=53150698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510074076.3A Active CN104616265B (en) | 2015-02-12 | 2015-02-12 | A kind of time domain method for reconstructing of remote sensing sequence data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104616265B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106327456A (en) * | 2016-08-19 | 2017-01-11 | 中国科学院遥感与数字地球研究所 | Method and device for information loss reconstruction of remote sensing image |
CN106934783A (en) * | 2017-03-02 | 2017-07-07 | 宁波大学 | A kind of time domain method for reconstructing of high frequency time remote sensing time series data |
CN110569733A (en) * | 2019-08-09 | 2019-12-13 | 中国科学院南京地理与湖泊研究所 | Lake long time sequence continuous water area change reconstruction method based on remote sensing big data platform |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020939A (en) * | 2012-12-18 | 2013-04-03 | 武汉大学 | Method for removing large-area thick clouds for optical remote sensing images through multi-temporal data |
CN103902802A (en) * | 2012-12-29 | 2014-07-02 | 中国科学院深圳先进技术研究院 | Method for rebuilding spatial-information-considered vegetation index time series data |
CN104318270A (en) * | 2014-11-21 | 2015-01-28 | 东北林业大学 | Land cover classification method based on MODIS time series data |
-
2015
- 2015-02-12 CN CN201510074076.3A patent/CN104616265B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020939A (en) * | 2012-12-18 | 2013-04-03 | 武汉大学 | Method for removing large-area thick clouds for optical remote sensing images through multi-temporal data |
CN103902802A (en) * | 2012-12-29 | 2014-07-02 | 中国科学院深圳先进技术研究院 | Method for rebuilding spatial-information-considered vegetation index time series data |
CN104318270A (en) * | 2014-11-21 | 2015-01-28 | 东北林业大学 | Land cover classification method based on MODIS time series data |
Non-Patent Citations (2)
Title |
---|
李天祺等: "环境星NDVI时间序列重构方法研究", 《CNKI》 * |
李洪利等: "一种高保真同态滤波遥感影像薄云去除方法", 《遥感应用》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106327456A (en) * | 2016-08-19 | 2017-01-11 | 中国科学院遥感与数字地球研究所 | Method and device for information loss reconstruction of remote sensing image |
CN106327456B (en) * | 2016-08-19 | 2019-09-17 | 中国科学院遥感与数字地球研究所 | Remote sensing images loss of learning method for reconstructing and device |
CN106934783A (en) * | 2017-03-02 | 2017-07-07 | 宁波大学 | A kind of time domain method for reconstructing of high frequency time remote sensing time series data |
CN106934783B (en) * | 2017-03-02 | 2020-01-31 | 宁波大学 | Time domain reconstruction method for high-frequency time-series remote sensing data |
CN110569733A (en) * | 2019-08-09 | 2019-12-13 | 中国科学院南京地理与湖泊研究所 | Lake long time sequence continuous water area change reconstruction method based on remote sensing big data platform |
CN110569733B (en) * | 2019-08-09 | 2022-02-01 | 中国科学院南京地理与湖泊研究所 | Lake long time sequence continuous water area change reconstruction method based on remote sensing big data platform |
Also Published As
Publication number | Publication date |
---|---|
CN104616265B (en) | 2018-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104599292B (en) | A kind of anti-noise moving object detection algorithm decomposed based on low-rank matrix | |
CN107169475B (en) | A kind of face three-dimensional point cloud optimized treatment method based on kinect camera | |
CN101137003B (en) | Gray associated analysis based sub-pixel fringe extracting method | |
CN106295542A (en) | A kind of road target extracting method of based on significance in night vision infrared image | |
CN103455991A (en) | Multi-focus image fusion method | |
CN105321172A (en) | SAR, infrared and visible light image fusion method | |
CN103578092A (en) | Multi-focus image fusion method | |
CN103020939A (en) | Method for removing large-area thick clouds for optical remote sensing images through multi-temporal data | |
CN106127688A (en) | A kind of super-resolution image reconstruction method and system thereof | |
CN105023013B (en) | The object detection method converted based on Local standard deviation and Radon | |
CN104236498A (en) | Method for measuring end size of switch rail | |
CN104616265A (en) | Time domain reconstruction method of remote sensing sequence data | |
CN104091350A (en) | Object tracking method achieved through movement fuzzy information | |
CN104197933A (en) | Method for enhancing and extracting star sliding in telescope field of view | |
CN104036229A (en) | Regression-based active appearance model initialization method | |
CN104615990A (en) | Method for automatically recognizing macula based on Huairou full-disk single-color image | |
CN104766341A (en) | Detection method for construction change of remote-sensing image based on DSM and kernel density estimation | |
Zhang et al. | Methods for automatic identification and extraction of terraces from high spatial resolution satellite data (China-GF-1) | |
CN102831621B (en) | Video significance processing method based on spectral analysis | |
CN104933719A (en) | Method for detecting image edge by integral image interblock distance | |
CN105491315A (en) | Projector gamma correction method | |
CN103325104B (en) | Based on the face image super-resolution reconstruction method of iteration sparse expression | |
CN103530628B (en) | High-resolution remote sensing image ortho-rectification method based on floating control point | |
CN104899592B (en) | A kind of road semiautomatic extraction method and system based on circular shuttering | |
CN103218820A (en) | Camera calibration error compensation method based on multi-dimensional characteristics |
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 |