CN111402169A - Method for repairing remote sensing vegetation index time sequence under influence of coastal tide - Google Patents
Method for repairing remote sensing vegetation index time sequence under influence of coastal tide Download PDFInfo
- Publication number
- CN111402169A CN111402169A CN202010209792.9A CN202010209792A CN111402169A CN 111402169 A CN111402169 A CN 111402169A CN 202010209792 A CN202010209792 A CN 202010209792A CN 111402169 A CN111402169 A CN 111402169A
- Authority
- CN
- China
- Prior art keywords
- vegetation index
- index
- vegetation
- water body
- time sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 66
- 238000012937 correction Methods 0.000 claims abstract description 32
- 238000002310 reflectometry Methods 0.000 claims description 12
- 230000003595 spectral effect Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000013535 sea water Substances 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 5
- 230000005855 radiation Effects 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 238000005067 remediation Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 description 9
- 238000012544 monitoring process Methods 0.000 description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 241001149258 Sporobolus alterniflorus Species 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000027066 STING-associated vasculopathy with onset in infancy Diseases 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 241001474374 Blennius Species 0.000 description 1
- 241001415959 Grus japonensis Species 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 240000002044 Rhizophora apiculata Species 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000009919 sequestration Effects 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a method for repairing remote sensing vegetation index time sequence under the influence of coastal tide, which is characterized by comprising the following steps: acquiring a vegetation index time sequence and a normalized difference water body index, acquiring sample pixels influenced by tides, identifying noise points of the time sequence, constructing a noise point value restoration model and restoring and evaluating the vegetation index time sequence, wherein each sample pixel in an experimental area is traversed, the noise point value restoration model is utilized to calculate the correction amount of the vegetation index noise points, the sum of the vegetation index noise point value and the corresponding vegetation index noise point correction amount value is used as the vegetation index correction value of the vegetation index noise points, and then the restored vegetation index time sequence is obtained; the method has the advantages of wide application range, high automation degree, strong robustness and suitability for different spatial resolutions.
Description
Technical Field
The invention relates to a method for repairing a remote sensing image vegetation index time sequence, in particular to a method for repairing a remote sensing vegetation index time sequence under the influence of tidal bodies on a coast.
Background
In coastal zone areas with closely interacting sea and land, vegetation (saline marsh, mangrove, seaweed and the like) is taken as an important carbon sink, and the 'blue carbon' process of soil is continuously promoted through spontaneous primary production and continuous shoal deposition. Remote sensing (RemoteSensing) provides a vegetation dynamic observation technical means in a large range and a long period, is beneficial to monitoring the primary productivity of the coastal vegetation, dynamically and quantitatively evaluates the carbon sequestration capacity change of the coastal zone ecosystem, and has very important significance. However, due to the periodic influence of ocean tides, the remote sensing of Vegetation on the coastal zones is particularly complicated, the spectral reflectivity of the Vegetation is obviously weakened by the submergence of seawater, noise is introduced into a Vegetation Index (VI) time sequence, and the time sequence characteristics (shape and size) of the Vegetation are further covered. Compared with a land ecosystem, the traditional vegetation index (NDVI, EVI, SAVI and the like) time sequence is low in the whole range of the coastal zone ecosystem, random noise is densely distributed, and time sequence mode distortion is obvious. Under the background, accurate identification and restoration of vegetation index noise under the influence of tides are urgently needed, and a coastal zone vegetation index time sequence is restored, so that the remote sensing technology is effectively applied to dynamic monitoring of a coastal zone ecosystem.
The Water body Index (Water Index, WI) can effectively distinguish the Water body on the remote sensing image from other earth surface coverage types, and provides an opportunity for coastal tidal inundation monitoring. At present, researchers use different image band combination modes to construct various water body indexes so as to meet the requirements of water body identification under different terrain backgrounds, and the method is successfully applied to the aspect of water body extraction in areas and even global areas. Notably, the water body index threshold selection is a prerequisite and key for identifying the water body, and is usually determined by image data of a single period and a short time span (the same month and the same quarter). However, the development of vegetation phenological periods and the difference of imaging time tide heights cause the relative difference between vegetation and seawater to change continuously, the threshold value selection mode in a single period is relatively unilateral, and the accuracy of the periodic tide submergence frequency drawn across multiple periods is poor. How to robustly determine the water body threshold in a time series and accurately identify vegetation index noise is a major current challenge. In addition, on the basis of identifying vegetation index noise, how to quantitatively depict the difference from the noise to the original vegetation index and construct a noise numerical value restoration model is also a technical difficulty to overcome when restoring the coastal vegetation index time sequence.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for repairing the remote sensing vegetation index time sequence under the influence of tide on the coast, which has the advantages of wide application range, high automation degree, strong robustness and is suitable for correcting various vegetation index time sequences of remote sensing images with different spatial resolutions.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for repairing a remote sensing vegetation index time sequence under the influence of coastal tide comprises the following steps:
(1) vegetation Index (NDVI) time series acquisition
Acquiring image pixels (DN) of a designated time and an experimental area through a European space Bureau database, converting image pixel values (DN) into surface reflectivity through radiation correction, calculating a vegetation index (NDVI) of each scene image, removing the pixels covered by clouds, cloud shadows and snow by using a roughly classified product corresponding to each scene image, sequentially arranging vegetation index images according to the sequence of imaging time, and constructing to obtain a vegetation index time sequence;
(2) normalized Differential Water Index (MNDWI) calculation
Calculating the normalized difference water body index (MNDWI) of each scene of image according to an improved normalized difference water body index formula, eliminating pixels covered by clouds, cloud shadows and snow by using a roughly classified product corresponding to each scene of image, and constructing and obtaining a water body index time sequence completely corresponding to the vegetation index time sequence;
(3) tide influenced sample pixel acquisition
Selecting image pixels which are not covered by snow in an experimental area and have cloud coverage rate lower than 10%, and acquiring typical vegetation pixels (which are mostly distributed in pioneer vegetation communities at river mouths, ports and sea sides) which are intermittently covered by seawater (ensuring observation times more than 7, exposing on the water surface at most times and submerging by the seawater at few times) as samples for identifying and repairing noise points;
(4) time series noise point identification
Each sample pixel corresponds to two time sequences of a vegetation index and a water body index, the overall variation trend of the vegetation index time sequence and the water body index time sequence of each sample pixel is observed, the abnormal minimum value of the vegetation index and the positive value of the water body index are recorded sample by sample at the same time, when the vegetation index time sequence and the water body index time sequence occur at the same time, the water body index numerical values marked by all samples are summarized, the minimum value of the water body index numerical values is taken as a water body index threshold value, and the vegetation index corresponding to the water body index numerical value in the water body index time sequence which is greater than the water body index threshold value is taken as a vegetation index noise pointnoise);
(5) Noise point numerical restoration model construction
For each sample pixel, removing vegetation index noise points, fitting the variation trend of the residual vegetation index in a vegetation index time sequence by adopting a second-order Fourier function, taking the vertical distance from the vegetation index noise points in the sample pixel to a Fourier function fitting value curve as a vegetation index noise point correction quantity, taking the vegetation index noise point correction quantity delta NDVI as a dependent variable, taking a water body index corresponding to the vegetation index noise points as an independent variable, drawing scatter diagrams of the vegetation index noise point correction quantity delta NDVI and the water body index, fitting the relation between the vegetation index noise point correction quantity delta NDVI and the water body index by adopting a polynomial equation, and quantitatively constructing a noise point numerical value restoration model as follows:
ΔNDVI=-3.84*MNDWI3+3.98*MNDWI2-0.33*MNDWI+0.14;
(6) vegetation index time series remediation assessment
Traversing each sample pixel of the experimental area, calculating the vegetation index noise point correction quantity by using a noise point numerical value restoration model, and taking the sum of the vegetation index noise point numerical value and the corresponding vegetation index noise point correction quantity numerical value as the vegetation index correction value VI of the vegetation index noise pointmAnd further obtaining the vegetation index time series after restoration.
The vegetation index (NDVI) of each scene image in the step (1) is calculated according to the following formula:wherein rho NIR is the surface spectral reflectivity of the near infrared band, rhoredThe surface spectral reflectivity of the red light wave band.
The improved normalized difference water body index (MNDWI) calculation formula in the step (2) is as follows:where ρ isgreenSpectral reflectance of the earth's surface in the green band, pswIRThe surface spectral reflectivity of the short wave infrared band.
The formula for fitting the variation trend of the residual vegetation index in the vegetation index time sequence by the second-order Fourier function in the step (5) is as follows: FF (a 0+ a1 Cos (w x) + b1 sin (w x) + a2 Cos (2w x) + b2 sin (2w x), where x is julian day; w is frequency, and for a year's time series values are typically 2 pi/365, w is 0.018, a0 is 0, a1 is 0, a2 is 0, b1 is 0, b2 is 0, and FF represents the fit of the fourier function to the index of the vegetation remaining.
Compared with the prior art, the invention has the advantages that:
firstly, a vegetation index time series noise identification method is provided. By selecting a small amount of vegetation pixels covered by seawater in a targeted manner, the water body index threshold value suitable for annual observation is automatically optimized, the noise of the vegetation index time sequence is accurately identified, and conditions are created for vegetation index time sequence restoration. The method has high identification precision and strong generalization capability, and can be applied to the vegetation index time sequence filtering and denoising process of other land and water interaction frequent areas (river banks and huff-and-puff lakes).
Secondly, a vegetation index time series numerical restoration model is constructed. And the annual variation trend of the vegetation index is fully described by a second-order Fourier function, and the noise correction quantity of the vegetation index is quantitatively measured and calculated. And constructing a polynomial numerical restoration model by depending on the correlation between the noise correction value and the water body index, thereby realizing the global batch correction of the vegetation index time sequence. The innovative method for compensating and correcting the missing data by using the fitting function is not only suitable for the field of remote sensing mapping, but also can provide beneficial reference for numerical restoration and data supplementation of sparse signal sequences in the field of communication.
In summary, the method solves the technical difficulties of identifying noise points, restoring noise values and the like of the time sequence of the vegetation index of the coastal zone under the influence of tide, overcomes the bottleneck that the time sequence of the vegetation index is low in the whole range of the coastal zone and low in applicability due to dense random noise, has high automation degree and strong robustness, is suitable for correcting the time sequences of various vegetation indexes (NDVI, EVI, SAVI and the like) of remote sensing images (MODIS, L andsat, Sentinel and the like) with different spatial resolutions, is expected to serve in the remote sensing time sequence monitoring of the coastal zone ecosystem, and lays a foundation for the effective application of the remote sensing time sequence technology to the dynamic monitoring of the coastal zone ecosystem.
Drawings
FIG. 1 is a time series restoration technique for index vegetation in a coastal zone affected by tide;
FIG. 2 shows a Sentiel-2 MSI image and sample pixel distribution thereof along the coast of Jiangsu;
FIG. 3 is an exponential time series noise point identification of the vegetation in the coastal zone before restoration;
FIG. 4 is a relational expression of a vegetation index noise correction quantity and a water body index based on a polynomial;
FIG. 5 is the sea-shore vegetation index time series noise point identification after the restoration;
FIG. 6 is a distribution of determinant coefficient differences fitted to a vegetation index time series before and after remediation.
Detailed Description
The invention is described in further detail below with reference to the accompanying examples.
The method for repairing the remote sensing vegetation index time sequence under the influence of the coastal tide is further explained by taking the natural protection area of the red-crowned crane in the middle part of Jiangsu along the sea as an experimental area and taking the Normalized Difference Vegetation Index (NDVI) constructed by the multispectral data (Sentinel-2MSI) of the Sentinel satellite No. 2 as experimental data, wherein the flow is shown in FIG. 1, and the specific steps are as follows:
the experimental data are derived from Sentinel 2 multispectral imager images (Sentinel-2MSI) provided by the European and aviation Bureau (https:// sciihub.copernius. eu /), through a designated time and an experimental area, a Sentinel-2MSI L1C product which is imaged in 68 scenes with the corresponding row code number of T51STT in 2018 and subjected to geometric correction is downloaded through experiments, the geometric precision meets the requirements of pixel-level time sequence construction and analysis, on the basis, Sen2Cor software provided by the European and aviation Bureau is adopted for radiation correction in the experiments, the image pixel value (DN) is converted into the surface reflectivity, the vegetation index (NDVI) of each scene image is calculated, and the calculation formula is as follows:where ρ isNIRSurface spectral reflectance, ρ, for the near infrared bandredThe surface spectral reflectivity of the red light wave band.
In the radiation correction process, a coarse classification product (Scene C L assessment, SC L) of each Scene image is automatically generated, pixels covered by clouds, cloud shadows and snow in each Scene NDVI image are removed through experiments, the NDVI images are sequentially arranged according to the imaging time sequence, and an NDVI time sequence is constructed.
And 2, calculating the normalized difference water body index, on the basis of obtaining the surface reflectivity, calculating the normalized difference water body index (MNDWI) of each scene of image according to an improved normalized water body index formula by an experiment, similarly, removing the pixels covered by clouds, cloud shadows and snow by using an SC L product corresponding to each scene of image, sequentially arranging the MNDWI images according to the imaging time sequence, and constructing and obtaining the MNDWI time sequence completely corresponding to the NDVI time sequence.
Wherein the Modified Normalized Difference Water Index (MNDWI) is calculated as follows:where ρ isgreenSpectral reflectance of the earth's surface in the green band, pSWIRThe surface spectral reflectivity of the short wave infrared band. MNDWI can highlight the difference of water and vegetation and restrain the confusion of buildings.
And step 3: and collecting sample pixels influenced by tide. In the 68 scenes in 2018, 17 image pixels with no snow coverage and less than 10% cloud coverage are selected and displayed by the standard false color scheme of R-4, G-3 and B-2. Through multi-temporal image contrast observation, the experiment uses the edge of Spartina alterniflora (Spartina alterniflora) salt marsh with frequent sea-land interaction as the target area selected by the vegetation sample. In order to ensure the objectivity and representativeness of sample selection, 39 sample pixels are equidistantly selected in a target area by 500m in the experiment (see the attached figure 2 for details).
And 4, step 4: and (4) identifying time series noise points. In the experiment, Matlab software is used for drawing an NDVI time sequence and an MNDWI time sequence of 39 sample pixels 2018, abnormal minimum values of vegetation indexes and positive values of water indexes, when the vegetation indexes and the water indexes appear at the same time, the water indexes are marked (see the attached figure 3 in detail), the water index values marked by all samples are summarized, and the minimum value of the water indexes is counted to be 0.023 and used as a water index threshold value. On the basis, traversing each pixel of the experimental area, and regarding the vegetation index corresponding to the water body index larger than the threshold value in the water body index time sequence as a vegetation index noise point (VI)noise)。
And 5: and constructing a noise point numerical value restoration model. According to the experiment, vegetation index noise points are removed from a sample pixel by sample pixel, and the variation trend of the residual vegetation index in the vegetation index time sequence is fitted by using a second-order Fourier function (Fourier) of a Matlab software Curve Fitting module. On the basis, the experiment leads vegetation index noise points in sample pixels to Fourier transformThe vertical distance of the function fitting value curve is used as the correction quantity of the vegetation index noise point, the correction quantity delta NDVI of the vegetation index noise point is used as a dependent variable, the water body index corresponding to the vegetation index noise point is used as an independent variable, and a scatter diagram is drawn as shown in an attached figure 4. Fitting the relation between the vegetation index correction quantity and the water body index by adopting a polynomial equation (the order is less than 4), and quantitatively constructing a noise point numerical restoration model: Δ NDVI ═ 3.84 × MNDWI3+3.98*MNDWI2-0.33 × MNDWI +0.14 wherein a second order Fourier Function (FF) is able to adequately characterize the cyclic variations in the coastal vegetation over the years due to seasonal temperature and precipitation environmental factors, which is calculated as follows:
FF (a 0+ a1 cos (w x) + b1 sin (w x) + a2 cos (2w x) + b2 sin (2w x), where x is julian day; w is frequency, and for a year's time series values are typically 2 pi/365, w is 0.018, a0 is 0, a1 is 0, a2 is 0, b1 is 0, b2 is 0, and FF represents the fit of the fourier function to the index of the vegetation remaining.
Step 6: and (4) vegetation index time series restoration evaluation. Traversing each pixel of the experimental area, calculating the vegetation index noise point correction quantity by using a noise point numerical value restoration model, and taking the sum of the vegetation index noise point numerical value and the corresponding vegetation index noise point correction quantity numerical value as a vegetation index correction value (NDVI)m) Thus, the vegetation index time series is repaired (see figure 5 in detail). Determining coefficients (R) by fitting a second order Fourier function2) The difference distribution shows that the effect of the invention is good, the fitting precision of the range of-40% experimental area is improved, especially the fitting precision of the outer edge area of the spartina alterniflora is averagely improved by 0.32, which is shown in figure 6 in detail.
The above description is not intended to limit the present invention, and the present invention is not limited to the above examples. Those skilled in the art should also realize that changes, modifications, additions and substitutions can be made without departing from the true spirit and scope of the invention.
Claims (4)
1. A method for repairing a remote sensing vegetation index time sequence under the influence of coastal tide, which is characterized by comprising the following steps:
(1) vegetation index time series acquisition
Acquiring image pixels of a designated time and an experimental area through a European space Bureau database, converting image pixel values into earth surface reflectivity through radiation correction, calculating vegetation indexes of each image, removing the pixels covered by clouds, cloud shadows and snow by using a roughly classified product corresponding to each image, sequentially arranging the vegetation index images according to the sequence of imaging time, and constructing and obtaining a vegetation index time sequence;
(2) normalized difference water body index calculation
Calculating the normalized difference water body index of each scene image according to an improved normalized difference water body index formula, eliminating pixels covered by clouds, cloud shadows and snow by using a rough classification product corresponding to each scene image, and constructing and obtaining a water body index time sequence completely corresponding to the vegetation index time sequence;
(3) tide influenced sample pixel acquisition
Selecting an image pixel which is not covered by snow in an experimental area and has a cloud coverage rate lower than 10%, and collecting a typical vegetation pixel which is intermittently covered by seawater to serve as a sample for identifying and repairing a noise point;
(4) time series noise point identification
Each sample pixel corresponds to two time sequences of a vegetation index and a water body index, the overall variation trend of the vegetation index time sequence and the water body index time sequence of each sample pixel is observed, the abnormal minimum value of the vegetation index and the positive value of the water body index are recorded sample by sample at the same time, when the vegetation index time sequence and the water body index time sequence occur at the same time, the water body index numerical values marked by all samples are summarized, the minimum value of the water body index numerical values is taken as a water body index threshold value, and the vegetation index corresponding to the water body index numerical value in the water body index time sequence which is greater than the water body index threshold value is taken as a vegetation index noise;
(5) noise point numerical restoration model construction
For each sample pixel, removing vegetation index noise points, fitting the variation trend of the residual vegetation index in a vegetation index time sequence by adopting a second-order Fourier function, taking the vertical distance from the vegetation index noise points in the sample pixel to a Fourier function fitting value curve as a vegetation index noise point correction quantity, taking the vegetation index noise point correction quantity delta NDVI as a dependent variable, taking a water body index corresponding to the vegetation index noise points as an independent variable, drawing scatter diagrams of the vegetation index noise points and the vegetation index noise point correction quantity delta NDVI and the water body index MNDWI, fitting the relation between the vegetation index noise point correction quantity delta NDVI and the water body index MNDWI by adopting a polynomial equation, and quantitatively constructing a noise point numerical value restoration:
ΔNDVI=-3.84*MNDWI3+3.98*MNDWI2-0.33*MNDWI+0.14;
(6) vegetation index time series remediation assessment
Traversing each sample pixel of the experimental area, calculating the vegetation index noise point correction quantity by using a noise point numerical value restoration model, and taking the sum of the vegetation index noise point numerical value and the corresponding vegetation index noise point correction quantity numerical value as the vegetation index correction value VI of the vegetation index noise pointmAnd further obtaining the vegetation index time series after restoration.
2. The method of repairing a time series of indexes of remote vegetation under the influence of coastal tidal conditions according to claim 1, wherein: the vegetation index calculation formula of each scene image in the step (1) is as follows:where ρ isNIRSurface spectral reflectance, ρ, for the near infrared bandredThe surface spectral reflectivity of the red light wave band.
3. The method of repairing a time series of indexes of remote vegetation under the influence of coastal tidal conditions according to claim 1, wherein: the improved normalized difference water body index calculation formula in the step (2) is as follows:where ρ isgreenSpectral reflectance of the earth's surface in the green band, pSWIRThe surface spectral reflectivity of the short wave infrared band.
4. The method of repairing a time series of indexes of remote vegetation under the influence of coastal tidal conditions according to claim 1, wherein: the formula for fitting the variation trend of the residual vegetation index in the vegetation index time sequence by the second-order Fourier function in the step (5) is as follows: FF (a 0+ a1 cos (w x) + b1 sin (w x) + a2 cos (2w x) + b2 sin (2w x), where x is julian day; w is frequency, and for a year's time series values are typically 2 pi/365, w is 0.018, a0 is 0, a1 is 0, a2 is 0, b1 is 0, b2 is 0, and FF represents the fit of the fourier function to the index of the vegetation remaining.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010209792.9A CN111402169B (en) | 2020-03-23 | 2020-03-23 | Method for repairing remote sensing vegetation index time sequence under influence of coastal tide |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010209792.9A CN111402169B (en) | 2020-03-23 | 2020-03-23 | Method for repairing remote sensing vegetation index time sequence under influence of coastal tide |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111402169A true CN111402169A (en) | 2020-07-10 |
CN111402169B CN111402169B (en) | 2023-04-11 |
Family
ID=71413997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010209792.9A Active CN111402169B (en) | 2020-03-23 | 2020-03-23 | Method for repairing remote sensing vegetation index time sequence under influence of coastal tide |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111402169B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113160359A (en) * | 2021-05-21 | 2021-07-23 | 宁波大学 | Coastal wetland remote sensing mapping method based on harmonic fitting parameters |
CN113591032A (en) * | 2021-07-15 | 2021-11-02 | 清华大学深圳国际研究生院 | Method for calculating optimal growth elevation of mangrove plant |
CN115035411A (en) * | 2022-06-17 | 2022-09-09 | 宁波大学 | Hyperspectral data-based spartina alterniflora index construction and extraction method |
CN115082809A (en) * | 2022-06-23 | 2022-09-20 | 宁波大学 | New tidal flat evolution monitoring method based on remote sensing image big data |
CN115078263A (en) * | 2022-05-27 | 2022-09-20 | 苏州科技大学 | Seaweed remote sensing information extraction method considering tidal influence |
CN116341932A (en) * | 2023-05-31 | 2023-06-27 | 自然资源部第二海洋研究所 | Tidal flat change monitoring method and system based on double remote sensing time sequence indexes |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090232349A1 (en) * | 2008-01-08 | 2009-09-17 | Robert Moses | High Volume Earth Observation Image Processing |
CN101650422A (en) * | 2009-09-27 | 2010-02-17 | 北京师范大学 | Denoising method of remote sensing vegetation index time series data |
CN103971027A (en) * | 2014-05-29 | 2014-08-06 | 福州大学 | Rice automatic mapping method based on ratio index of water body and vegetation index changes |
CN110378576A (en) * | 2019-07-01 | 2019-10-25 | 中国环境科学研究院 | The quantification detection method of urbanization vegetation effect effective distance |
-
2020
- 2020-03-23 CN CN202010209792.9A patent/CN111402169B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090232349A1 (en) * | 2008-01-08 | 2009-09-17 | Robert Moses | High Volume Earth Observation Image Processing |
CN101650422A (en) * | 2009-09-27 | 2010-02-17 | 北京师范大学 | Denoising method of remote sensing vegetation index time series data |
CN103971027A (en) * | 2014-05-29 | 2014-08-06 | 福州大学 | Rice automatic mapping method based on ratio index of water body and vegetation index changes |
CN110378576A (en) * | 2019-07-01 | 2019-10-25 | 中国环境科学研究院 | The quantification detection method of urbanization vegetation effect effective distance |
Non-Patent Citations (2)
Title |
---|
刘婷;包广道;张大伟;何怀江;罗也;张忠辉;: "基于植被指数的GF-2影像防护林快速提取研究" * |
刘永超;李加林;张宇;赵冰雪;许文轩;何改丽;刘永学;: "时间序列遥感影像支持下杭州湾南岸植被动态监测" * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113160359A (en) * | 2021-05-21 | 2021-07-23 | 宁波大学 | Coastal wetland remote sensing mapping method based on harmonic fitting parameters |
CN113591032A (en) * | 2021-07-15 | 2021-11-02 | 清华大学深圳国际研究生院 | Method for calculating optimal growth elevation of mangrove plant |
CN113591032B (en) * | 2021-07-15 | 2023-06-16 | 清华大学深圳国际研究生院 | Method for calculating optimal growth elevation of mangrove plant |
CN115078263A (en) * | 2022-05-27 | 2022-09-20 | 苏州科技大学 | Seaweed remote sensing information extraction method considering tidal influence |
CN115035411A (en) * | 2022-06-17 | 2022-09-09 | 宁波大学 | Hyperspectral data-based spartina alterniflora index construction and extraction method |
CN115082809A (en) * | 2022-06-23 | 2022-09-20 | 宁波大学 | New tidal flat evolution monitoring method based on remote sensing image big data |
CN115082809B (en) * | 2022-06-23 | 2023-02-17 | 宁波大学 | New tidal flat evolution monitoring method based on remote sensing image big data |
CN116341932A (en) * | 2023-05-31 | 2023-06-27 | 自然资源部第二海洋研究所 | Tidal flat change monitoring method and system based on double remote sensing time sequence indexes |
CN116341932B (en) * | 2023-05-31 | 2023-08-22 | 自然资源部第二海洋研究所 | Tidal flat change monitoring method and system based on double remote sensing time sequence indexes |
Also Published As
Publication number | Publication date |
---|---|
CN111402169B (en) | 2023-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111402169B (en) | Method for repairing remote sensing vegetation index time sequence under influence of coastal tide | |
Wang et al. | Tracking annual changes of coastal tidal flats in China during 1986–2016 through analyses of Landsat images with Google Earth Engine | |
Sun et al. | Classification mapping of salt marsh vegetation by flexible monthly NDVI time-series using Landsat imagery | |
CN109190538B (en) | Sediment-laden river delta coastal zone evolution analysis method based on remote sensing technology | |
Zhao et al. | Mapping large-area tidal flats without the dependence on tidal elevations: A case study of Southern China | |
Sánchez-García et al. | An efficient protocol for accurate and massive shoreline definition from mid-resolution satellite imagery | |
Lathrop et al. | A multi-scale segmentation approach to mapping seagrass habitats using airborne digital camera imagery | |
Dong et al. | Coral reef geomorphology of the Spratly Islands: A simple method based on time-series of Landsat-8 multi-band inundation maps | |
Tong et al. | An optimal waterline approach for studying tidal flat morphological changes using remote sensing data: A case of the northern coast of Vietnam | |
Ferrentino et al. | On the use of multipolarization satellite SAR data for coastline extraction in harsh coastal environments: The case of Solway Firth | |
CN107688776B (en) | Urban water body extraction method | |
Taghadosi et al. | Trend analysis of soil salinity in different land cover types using Landsat time series data (case study Bakhtegan Salt Lake) | |
CN109801306B (en) | Tidal ditch extraction method based on high resolution remote sensing image | |
CN115810155B (en) | Tidal wetland classification method | |
Jung et al. | A multi-sensor approach for detecting the different land covers of tidal flats in the German Wadden Sea—A case study at Norderney | |
Peng et al. | Accurate mapping of Chinese coastal aquaculture ponds using biophysical parameters based on Sentinel-2 time series images | |
CN114201732A (en) | Sentinel-2A image-based shallow sea water depth inversion method | |
CN114724049A (en) | Inland culture pond water surface identification method based on high-resolution remote sensing image data | |
Zhao et al. | Stability evaluation of tidal flats based on time-series satellite images: A case study of the Jiangsu central coast, China | |
Tsai | Monitoring 23-year of shoreline changes of the Zengwun Estuary in Southern Taiwan using time-series Landsat data and edge detection techniques | |
Wicaksono et al. | Semi-automatic shoreline extraction using water index transformation on Landsat 8 OLI imagery in Jepara Regency | |
Wibisana et al. | Preliminary study of total suspended solid distribution in coastal ujung Pangkah Gresik based reflectance value of Landsat satellite imagery | |
CN110274858B (en) | Remote sensing method for estimating lake suspended sediment concentration by utilizing GOCI data | |
Eom et al. | Tidal channel distribution in relation to surface sedimentary facies based on remotely sensed data | |
Toro Herrera et al. | Semi-automated production and filtering of satellite derived water quality parameters |
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 |