CN116523859A - Coastal zone vegetation weather parameter satellite image calibration method combined with ground observation - Google Patents
Coastal zone vegetation weather parameter satellite image calibration method combined with ground observation Download PDFInfo
- Publication number
- CN116523859A CN116523859A CN202310439447.8A CN202310439447A CN116523859A CN 116523859 A CN116523859 A CN 116523859A CN 202310439447 A CN202310439447 A CN 202310439447A CN 116523859 A CN116523859 A CN 116523859A
- Authority
- CN
- China
- Prior art keywords
- ground
- satellite
- vegetation
- image
- weather
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000007781 pre-processing Methods 0.000 claims abstract description 11
- 238000011160 research Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 19
- 239000004973 liquid crystal related substance Substances 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 10
- 238000012163 sequencing technique Methods 0.000 claims description 8
- 230000015556 catabolic process Effects 0.000 claims description 7
- 238000006731 degradation reaction Methods 0.000 claims description 7
- 230000000737 periodic effect Effects 0.000 claims description 7
- 230000001174 ascending effect Effects 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 230000003203 everyday effect Effects 0.000 claims description 5
- 238000003786 synthesis reaction Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 claims description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 240000001949 Taraxacum officinale Species 0.000 description 2
- 235000005187 Taraxacum officinale ssp. officinale Nutrition 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000002354 daily effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000029553 photosynthesis Effects 0.000 description 1
- 238000010672 photosynthesis Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a satellite image calibration method for vegetation weathers parameters of a coastal zone combined with ground observation, which comprises the following steps: acquiring a ground camera photo and a satellite remote sensing image; respectively carrying out data preprocessing and fitting a GCC time sequence observation curve and an NDVI time sequence observation curve on the ground camera photo and the satellite remote sensing image; and calculating each pixel climatic parameter covering the image of the research area by applying the optimal threshold value, thereby obtaining the climatic spatial distribution. The beneficial effects of the invention are as follows: according to the invention, the object-weather curve inflection point is used for automatically acquiring object-weather parameter information of the ground object-weather photo, which is densely and continuously observed, and satellite image object-weather parameter information which is sparsely observed and is unevenly distributed is extracted through a continuously-changing threshold value so as to match the object-weather parameter information of the photo, thereby constructing a time sequence object-weather data set which is attached to the sense of real object-weather.
Description
Technical Field
The invention relates to the technical field of remote sensing time sequence processing, in particular to a satellite image calibration method for vegetation weathers parameters of a coastal zone combined with ground observation.
Background
Vegetation weather is a biological phenomenon that natural plants are influenced by genetic factors and surrounding environments to generate periodical changes, and is an important ecosystem parameter for representing the dynamic state of an ecosystem and the response of the ecological system to the environmental changes. The coastal zone is in the land and ocean system staggered transition zone, and is one of the most productive and valuable ecosystems worldwide. Climate (especially temperature) changes have a significant impact on coastal zone vegetation climate; at the same time, the vegetation climate change will affect the coastal zone ecosystem structural function and further react to the coastal climate system. Therefore, the method has important significance in improving the accuracy of the global land-sea ecological model on carbon, water and energy cycle simulation and prediction and the like for long-term observation and accurate description of vegetation waiting periods of coastal zones.
The satellite images can acquire a large-scale space continuous earth observation, so that the research on the relation between vegetation weather and climate change on the national, intercontinental and even global scale becomes possible. In order to ensure sufficient observation quantity, the high-time-resolution satellite image synthesized product is firstly used for observing the weather information. However, the above-mentioned image products have mixed pixel climate information because of a relatively coarse spatial resolution (250 m to 8 km), and are difficult to be applied to a coastal wetland region with a narrow distribution range and complex ground feature elements. The medium-fraction (less than or equal to 30 m) satellite images are gradually used for acquiring vegetation climate information of the coastal wetland. Notably, the extraction of vegetation weather parameters based on the medium-resolution satellite images is mainly determined by fitting a time-series weather curve of remote sensing observation and extracting change characteristic nodes of the weather curve. However, the time resolution of the middle satellite image is generally low, and the influence of cloud and rain weather in the coastal zone is also caused, so that the effective annual observation quantity is small, and the accuracy of extracting the vegetation weather parameters in the coastal zone is often difficult to guarantee.
The use of a terrestrial weather camera (Phenocam) to obtain continuous observation of the visible light range of vegetation becomes one of the important means for remote sensing weather observation in recent years. At present, the weather cameras are gradually distributed in a plurality of countries and ecological systems worldwide, the rapid development of the weather camera observation network promotes the research of the coastal wetland weather observed by the ground camera to be increased. While ground observation with a weatherometer can provide accurate vegetation weatherometer information at a specific location, spatial dispersion of observations presents a significant challenge in constructing a large range of weatherometer data sets.
Therefore, the method fully plays the ground-satellite observation advantages, and combines the ground observation to determine the coastal zone vegetation weather parameters of the satellite image, thereby being beneficial to realizing the coastal zone vegetation weather monitoring with large range, long period and high precision and providing scientific reference for the coastal zone resource development and ecological protection.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a satellite image calibration method for vegetation weather parameters of a coastal zone combined with ground observation.
In a first aspect, a method for calibrating satellite images of vegetation weathers parameters of a coastal zone in combination with ground observation is provided, including:
step 1, acquiring a data source, wherein the data source comprises a ground camera photo from the same year and a satellite remote sensing image;
step 2, carrying out data preprocessing on the ground camera photo, including green coordinate index GCC calculation and percentile value denoising, and constructing photo GCC time sequence observation;
step 3, fitting a photo GCC time sequence observation curve by using a harmonic function, and extracting the time corresponding to the inflection point of the curve as a real earth surface weather parameter;
step 4, combining the vision analysis and GPS correction to determine the corresponding vision range of the ground camera photo;
step 5, carrying out data preprocessing on the satellite remote sensing image, including normalized difference vegetation index NDVI calculation and inverse distance weighted synthesis, and constructing an image NDVI time sequence observation in a corresponding view range;
step 6, fitting an image NDVI time sequence observation curve by using a harmonic function, acquiring image weather parameters through different thresholds, and selecting the threshold closest to the ground surface real weather parameters as an optimal threshold;
and 7, calculating each pixel climatic parameter of the image of the coverage research area by applying an optimal threshold value, so as to obtain the climatic spatial distribution.
Preferably, step 2 includes:
step 2.1, calculating a greenness coordinate index, wherein the formula is as follows:
wherein GCC represents the green degree coordinate index, band Blue 、Band Green 、Band Red Blue light, green light and red light wave bands in the photos of the ground camera are respectively corresponding;
step 2.2, calculating GCC for all photos acquired every day, and counting the average value of the GCC of vegetation parts in the photos;
step 2.3, ascending order of GCC mean value is ordered, and a numerical value corresponding to the set quantile is selected to be output as an effective value;
and 2.4, sequencing the effective values according to the acquisition date to obtain photo GCC time sequence observation.
Preferably, step 3 includes:
step 3.1, fitting a photo GCC time sequence observation curve by using a harmonic function, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the method is characterized in that the method is used for estimating a harmonic function of GCC, t is julian day, w is annual frequency, and a constant 2 pi/365 is taken; i is the order of the harmonic function, a i And b i Is a harmonic function parameter, and is used for describing periodic changes in the year; c is the intercept, used to describe the change in annual trends; solving harmonic function parameters by adopting a general least square method;
step 3.2, calculating a second derivative of the harmonic function curve, and solving a corresponding time vector T when the second derivative is zero; the time vector T comprises two time points, which are expressed as T= (T1, T2), and the previous time point T1 is a vegetation growth period weather parameter SOS observed on the ground ground The latter time point t2 is the vegetation degradation period physical parameter EOS observed on the ground ground 。
Preferably, step 4 includes:
step 4.1, combining the ground camera height, azimuth angle, inclination angle and digital elevation model parameters, and obtaining a rough vision range by utilizing the vision analysis;
and 4.2, the observer holds the GPS to move along the boundary of the rough vision range, and the position feedback correction of the observer in the photo is carried out to determine the fine vision range.
Preferably, step 5 includes:
step 5.1, calculating a normalized difference vegetation index, wherein the formula is as follows:
wherein NDVI represents normalized difference vegetation index, band Red 、Band NIR Respectively corresponding to red light and near infrared wave bands in the satellite remote sensing image;
step 5.2, synthesizing the satellite images NDVI in the photo view range by using an inverse distance weighting mode, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,is a synthesized sanitationThe satellite images NDVI, i are the i-th satellite image pixels falling into the field of view of the weather photo, and NDVI i For pixel i, corresponding to NDVI, d i The distance from the pixel i to the ground camera;
step 5.3, willAnd sequencing according to the acquisition date to obtain the satellite image NDVI time sequence observation.
Preferably, step 6 includes:
step 6.1, fitting an image NDVI time sequence observation curve by using a harmonic function, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the method is characterized by taking a harmonic function estimated value of NDVI, wherein t is julian day, w is annual frequency, and a constant of 2 pi/365 is taken; i is the order of the harmonic function, a i And b i Is a harmonic function parameter, and is used for describing periodic changes in the year; c is the intercept, used to describe the change in annual trends; solving harmonic function parameters by adopting a general least square method;
step 6.2, acquiring vegetation growth period climate parameters SOS from satellite images satellite，x1 Defining the corresponding time of the rise of the harmonic function curve from the minimum value to the maximum value by x1%, and acquiring a vegetation degradation period weathered parameter EOS obtained by the satellite remote sensing image satellite，x2 Defining the corresponding time of the decrease of the harmonic function curve from the maximum value to the minimum value by x 2%;
and 6.3, taking x1 and x2 as undetermined thresholds, and trying with the step length of 5 from 20 to 80, wherein when the satellite image vegetation weather is closest to the ground observation vegetation weather, the x1 and x2 are the optimal thresholds, and the formula is as follows:
|SOS satellite，x1 -SOS ground |→min，x1＝20，25，…，80
|EOS satellite，x2 -EOS ground |→min，x2＝20，25，…，80。
in a second aspect, a satellite image calibration device for a vegetation weatherometer of a coastal zone in joint ground observation is provided, and the satellite image calibration method for the vegetation weatherometer of the coastal zone in joint ground observation in any aspect is executed, and includes:
the acquisition module is used for acquiring a data source, wherein the data source comprises a ground camera photo from the same year and a satellite remote sensing image;
the first construction module is used for carrying out data preprocessing on the photos of the ground camera, and comprises the steps of calculating a green coordinate index GCC and denoising a percentile value, so as to construct a photo GCC time sequence observation;
the first fitting module is used for fitting a photo GCC time sequence observation curve by using a harmonic function and extracting the time corresponding to the inflection point of the curve as a real earth surface weather parameter;
the correction module is used for combining the visibility analysis and GPS correction to determine the corresponding view range of the ground camera photo;
the second construction module is used for carrying out data preprocessing on the satellite remote sensing image, and comprises normalized difference vegetation index NDVI calculation and inverse distance weighted synthesis, and constructing an image NDVI time sequence observation in a corresponding view range;
the second fitting module is used for fitting an image NDVI time sequence observation curve by using a harmonic function, acquiring image weather parameters through different thresholds, and selecting the threshold closest to the ground surface real weather parameters as an optimal threshold;
and the calculation module is used for applying an optimal threshold value to calculate each pixel climatic parameter of the image of the coverage research area so as to obtain the climatic space distribution.
In a third aspect, a computer storage medium having a computer program stored therein is provided; the computer program when run on a computer causes the computer to execute the satellite image calibration method for vegetation parameters of the coastal zone in combination with ground observation according to any one of the first aspects.
The beneficial effects of the invention are as follows: the invention innovatively provides a satellite image coastal zone vegetation weather rating method combined with ground observation, which automatically acquires densely-observed continuous ground weather photo weather parameter information by using a weather curve inflection point, extracts satellite image weather parameter information which is sparsely observed and is unevenly distributed by continuously changing a threshold value to match the photo weather parameter information, and thereby constructs a time sequence weather data set which is fit with the sense of real weather. The data set has more definite biophysical connotation, can refract photosynthesis and carbon sink efficiency annual progress, and is further hopeful to serve the monitoring and evaluation work of the ecological environment of the coastal zone. In addition, as important vegetation biophysical information, a satellite image object data set with large range, long period and high precision is expected to expand research contents of cross subjects such as ocean, ecology, environment and the like.
Drawings
FIG. 1 is a schematic view of Landsat-8OLI images and a ground camera photograph of the natural protection area of Dandelion crane in 2020;
FIG. 2 is a schematic diagram of a daily land photograph pretreatment process;
FIG. 3 is a schematic diagram of a process for matching the vegetation ground observation weathers with the satellite image weathers in the coastal zone;
FIG. 4 is a schematic diagram of a process for matching the field of view of a terrestrial camera with satellite image pixels;
FIG. 5 is a schematic view showing the spatial distribution of the climate parameters in the natural protection zone of the Dan Ding Ji in 2020;
FIG. 6 is a flow chart of a method for satellite image calibration of the vegetation weathers parameters of the coastal zone in combination with ground observation.
Detailed Description
The invention is further described below with reference to examples. The following examples are presented only to aid in the understanding of the invention. It should be noted that it will be apparent to those skilled in the art that modifications can be made to the present invention without departing from the principles of the invention, and such modifications and adaptations are intended to be within the scope of the invention as defined in the following claims.
Example 1:
a satellite image calibration method for vegetation weathers of coastal zone combined with ground observation is shown in fig. 6, and comprises the following steps:
and step 1, acquiring a data source, wherein the data source comprises a ground camera photo from the same year and a satellite remote sensing image.
And 2, carrying out data preprocessing on the ground camera photo, including calculation of a green coordinate index GCC and denoising of a percentile value, and constructing a photo GCC time sequence observation.
The step 2 comprises the following steps:
step 2.1, calculating a greenness coordinate index, wherein the formula is as follows:
wherein GCC represents the green degree coordinate index, band Blue 、Band Green 、Band Red Blue light, green light and red light wave bands in the photos of the ground camera are respectively corresponding;
step 2.2, calculating GCC for all photos acquired every day, and counting the average value of the GCC of vegetation parts in the photos;
step 2.3, ascending order of GCC mean value is ordered, and a numerical value corresponding to the set quantile is selected to be output as an effective value; in this embodiment, the quantile may be set to 90%
And 2.4, sequencing the effective values according to the acquisition date to obtain photo GCC time sequence observation.
In the step 2, the GCC mean value is calculated for the vegetation parts in the photos collected every day, and the specific quantile is selected as the effective value output of the current day after ascending order, so that the influence of illumination and atmospheric condition differences at different moments on the GCC value can be avoided.
And step 3, fitting a photo GCC time sequence observation curve by using a harmonic function, and extracting the time corresponding to the inflection point of the curve as the real earth surface weather parameter.
The step 3 comprises the following steps:
step 3.1, fitting a photo GCC time sequence observation curve by using a harmonic function, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the method is characterized in that the method is used for estimating a harmonic function of GCC, t is julian day, w is annual frequency, and a constant 2 pi/365 is taken; i is the order of the harmonic function, a i And b i Is a harmonic function parameter, and is used for describing periodic changes in the year; c is the intercept, used to describe the change in annual trends; and solving harmonic function parameters by adopting a general least square method.
As can be seen by combining the subsequent step 6, the application adopts the harmonic function as a fitting curve of the time series observation of the ground photo and the satellite image, so that the application comprehensively considers the periodic and trend change characteristics of the vegetation in the coastal zone.
Step 3.2, calculating a second derivative of the harmonic function curve, and solving a time vector T corresponding to zero (curve inflection point) of the second derivative; the time vector T comprises two time points, which are expressed as T= (T1, T2), and the previous time point T1 is a vegetation growth period weather parameter SOS observed on the ground ground The latter time point t2 is the vegetation degradation period physical parameter EOS observed on the ground ground . According to the method, the corresponding time when the inflection point of the GCC fitting curve is solved, and the time is used as the physical parameters of vegetation growth period and decay period of ground observation, so that the actual sensory result of the vegetation change of the ground surface can be attached
And 4, combining the visibility analysis and GPS correction to determine the corresponding view range of the ground camera photo, and obtaining the view range of the ground camera from thick to thin.
Step 4 comprises:
step 4.1, combining the ground camera height, azimuth angle, inclination angle and digital elevation model parameters, and obtaining a rough vision range by utilizing the vision analysis;
and 4.2, the observer holds the GPS to move along the boundary of the rough vision range, and the position feedback correction of the observer in the photo is carried out to determine the fine vision range.
And 5, performing data preprocessing on the satellite remote sensing image, including normalized difference vegetation index NDVI calculation and inverse distance weighted synthesis, and constructing an image NDVI time sequence observation in a corresponding view range.
The step 5 comprises the following steps:
step 5.1, calculating a normalized difference vegetation index, wherein the formula is as follows:
wherein NDVI represents normalized difference vegetation index, band Red 、Band NIR Respectively corresponding to red light and near infrared wave bands in the satellite remote sensing image;
step 5.2, synthesizing the satellite images NDVI in the photo view range by using an inverse distance weighting mode, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,for the synthesized satellite image NDVI, i is the i-th satellite image pixel falling into the visible range of the weather photo, and NDVI i For pixel i, corresponding to NDVI, d i The distance from the pixel i to the ground camera;
step 5.3, willAnd sequencing according to the acquisition date to obtain the satellite image NDVI time sequence observation.
In step 5, the present application synthesizes the satellite images NDVI in the range of the photo view by using the inverse distance weighting method, and then sequences the synthesized NDVI according to the acquisition date to obtain the time series observation of the satellite images NDVI, so as to fully consider the imaging characteristics of the near-far dimension of the photo,
and 6, fitting an image NDVI time sequence observation curve by using a harmonic function, acquiring image weather parameters through different thresholds, and selecting the threshold closest to the ground surface real weather parameters as an optimal threshold.
The step 6 comprises the following steps:
step 6.1, fitting an image NDVI time sequence observation curve by using a harmonic function, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the method is characterized by taking a harmonic function estimated value of NDVI, wherein t is julian day, w is annual frequency, and a constant of 2 pi/365 is taken; i is the order of the harmonic function, a i And b i Is a harmonic function parameter, and is used for describing periodic changes in the year; c is the intercept, used to describe the change in annual trends; solving harmonic function parameters by adopting a general least square method;
step 6.2, acquiring vegetation growth period climate parameters SOS from satellite images satellite,x1 Defining the corresponding time of the rise of the harmonic function curve from the minimum value to the maximum value by x1%, and acquiring a vegetation degradation period weathered parameter EOS obtained by the satellite remote sensing image satellite,x2 Defining the corresponding time of the decrease of the harmonic function curve from the maximum value to the minimum value by x 2%;
and 6.3, taking x1 and x2 as undetermined thresholds, and trying with the step length of 5 from 20 to 80, wherein when the satellite image vegetation weather is closest to the ground observation vegetation weather, the x1 and x2 are the optimal thresholds, and the formula is as follows:
|SOS satellite，x1 -SOS ground |→min，x1＝20，25，…，80
|EOS satellite，x2 -EOS ground |→min，x2＝20，25，…，80
in step 6, the method adopts a continuous change threshold segmentation method to determine the satellite image weather parameters closest to the ground observation weather parameters, so that the accuracy of satellite image vegetation weather extraction is ensured.
And 7, calculating each pixel climatic parameter of the image of the coverage research area by applying an optimal threshold value, so as to obtain the climatic spatial distribution.
Example 2:
the natural protection area of the Dan top crane is positioned on the coast between the Jiulong port and the Xinyang port of salt city in Jiangsu province, and is considered to be one of the coastal zone wetlands which remain relatively intact in China. In order to make the content and advantages of the invention more clear, the following description uses the natural protection area of the dan crane as a case and the accompanying drawings to illustrate the specific implementation process of the invention, and the detailed steps are as follows:
step 1, as shown in fig. 1, downloading cloud-free coverage Landsat-8OLI and Landsat-7 ETM+ image 21 scenes in 2020 of a natural protection area covering a Dandelion crane; in addition, as shown in part (a) of fig. 2, a land camera photo 9862 of the year 2020 of the north bank of the natural protection area of the dan-top crane is obtained.
And 2, calculating GCC for all the ground photos collected every day, counting vegetation average values in the photos, and outputting 90% quantiles as GCC effective values of the same day after ascending order sequencing. Taking 5/13/2020 as an example, as shown in part (b) of fig. 2, 24 land photographs GCC are calculated first, and then the average ascending order of vegetation GCC is counted, wherein 90% of the digits correspond to 0.375, which is the effective value of GCC on 13/2020. Photo GCC time series observations are obtained by sequencing the GCC effective values daily in 2020 according to the dates.
Step 3, as shown in part (a) of fig. 3, solving a harmonic function fitting curve by using a least square method through GCC time series observation in 2020, wherein the specific expression is as follows:
the expression of the second derivative curve is further calculated as:
the second derivative curve expression is made to be zero, the time vector T is solved, and the vegetation growth period climate parameters SOS of ground observation are obtained ground 157 of (157)Ground-observed vegetation degradation period weathered parameter EOS ground 289, as shown in part (b) of fig. 3.
And 4, the ground camera height is 2.5m, the azimuth angle is 285 degrees, the inclination angle is 30 degrees, the digital elevation model is downloaded from USGS Earth Explorer, and relevant parameters and data are input into a Viewshed tool of ArcGIS 10.3 software to obtain a rough visual range. On this basis, the observer moves along the rough viewing area boundary, the relative position of the observer and the viewing area boundary in the photo is recorded, and the fine viewing area range is obtained by correcting the observer GPS information as shown in fig. 4.
Step 5, calculating NDVI for 21 scene Landsat images, and synthesizing pixels positioned in the field of view of the ground camera by adopting an inverse distance weighting formulaTaking Landsat-8OLI of 7.6.2020 as an example, the synthesized value of pixel NDVI +.>0.431. All images +.>And obtaining the time series observation of the image NDVI according to the date sequence.
Step 6, part (c) of fig. 3 shows the fitting of the harmonic function of the NDVI time series of the satellite image and the matching of the physical parameters thereof, and the fitting curve of the harmonic function is solved by using the least square method through the NDVI time series observation in 2020, and the specific expression is as follows:
during successive attempts from 20 to 80 in steps 5 of the pending threshold x1, x2, it was found that when x1=45, sos satellite,45 =159 and SOS ground Closest, eos when x2=40 satellite,40 =286 and EOS ground Nearest, 45 and 40 are the optimal thresholds for x1 and x2, respectively.
Step 7, let x1=45 Applying x2 = 40 to calculate each pixel climatic parameter of the image covering the natural protection region of the Dan Ding Crane, as shown in parts (a) and (b) of fig. 5, respectively obtaining vegetation period (SOS) satellite,45 ) And a recession phase (EOS) satellite,40 ) Is spatially distributed.
Claims (8)
1. The method for calibrating the satellite images of the vegetation weathers of the coastal zone by combining ground observation is characterized by comprising the following steps of:
step 1, acquiring a data source, wherein the data source comprises a ground camera photo from the same year and a satellite remote sensing image;
step 2, carrying out data preprocessing on the ground camera photo, including green coordinate index GCC calculation and percentile value denoising, and constructing photo GCC time sequence observation;
step 3, fitting a photo GCC time sequence observation curve by using a harmonic function, and extracting the time corresponding to the inflection point of the curve as a real earth surface weather parameter;
step 4, combining the vision analysis and GPS correction to determine the corresponding vision range of the ground camera photo;
step 5, carrying out data preprocessing on the satellite remote sensing image, including normalized difference vegetation index NDVI calculation and inverse distance weighted synthesis, and constructing an image NDVI time sequence observation in a corresponding view range;
step 6, fitting an image NDVI time sequence observation curve by using a harmonic function, acquiring image weather parameters through different thresholds, and selecting the threshold closest to the ground surface real weather parameters as an optimal threshold;
and 7, calculating each pixel climatic parameter of the image of the coverage research area by applying an optimal threshold value, so as to obtain the climatic spatial distribution.
2. The method for satellite image calibration of coastal zone vegetation weatherometer parameters in combination with ground observation according to claim 1, wherein step 2 comprises:
step 2.1, calculating a greenness coordinate index, wherein the formula is as follows:
wherein GCC represents the green degree coordinate index, band Blue 、Band Green 、Band Red Blue light, green light and red light wave bands in the photos of the ground camera are respectively corresponding;
step 2.2, calculating GCC for all photos acquired every day, and counting the average value of the GCC of vegetation parts in the photos;
step 2.3, ascending order of GCC mean value is ordered, and a numerical value corresponding to the set quantile is selected to be output as an effective value;
and 2.4, sequencing the effective values according to the acquisition date to obtain photo GCC time sequence observation.
3. The method for satellite image calibration of coastal zone vegetation weatherometer parameters in combination with ground observation according to claim 2, wherein step 3 comprises:
step 3.1, fitting a photo GCC time sequence observation curve by using a harmonic function, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the method is characterized in that the method is used for estimating a harmonic function of GCC, t is julian day, w is annual frequency, and a constant 2 pi/365 is taken; i is the order of the harmonic function, a i And b i Is a harmonic function parameter, and is used for describing periodic changes in the year; c is the intercept, used to describe the change in annual trends; solving harmonic function parameters by adopting a general least square method;
step 3.2, calculating a second derivative of the harmonic function curve, and solving a corresponding time vector T when the second derivative is zero; the time vector T comprises two time points, which are expressed as T= (T1, T2), and the previous time point T1 is a vegetation growth period weather parameter SOS observed on the ground ground The latter time point t2 is the vegetation degradation period physical parameter EOS observed on the ground ground 。
4. The method for satellite image calibration of coastal zone vegetation weatherometer according to claim 3, wherein step 4 comprises:
step 4.1, combining the ground camera height, azimuth angle, inclination angle and digital elevation model parameters, and obtaining a rough vision range by utilizing the vision analysis;
and 4.2, the observer holds the GPS to move along the boundary of the rough vision range, and the position feedback correction of the observer in the photo is carried out to determine the fine vision range.
5. The method for satellite image calibration of coastal zone vegetation weatherometer according to claim 4, wherein step 5 comprises:
step 5.1, calculating a normalized difference vegetation index, wherein the formula is as follows:
wherein NDVI represents normalized difference vegetation index, band Red 、Band NIR Respectively corresponding to red light and near infrared wave bands in the satellite remote sensing image;
step 5.2, synthesizing the satellite images NDVI in the photo view range by using an inverse distance weighting mode, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,for the synthesized satellite image NDVI, i is the i-th satellite image pixel falling into the visible range of the weather photo, and NDVI i Corresponds to the pixel iNDVI，d i The distance from the pixel i to the ground camera;
step 5.3, willAnd sequencing according to the acquisition date to obtain the satellite image NDVI time sequence observation.
6. The method for satellite image calibration of coastal zone vegetation weatherometer according to claim 5, wherein step 6 comprises:
step 6.1, fitting an image NDVI time sequence observation curve by using a harmonic function, wherein the calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the method is characterized by taking a harmonic function estimated value of NDVI, wherein t is julian day, w is annual frequency, and a constant of 2 pi/365 is taken; i is the order of the harmonic function, a i And b i Is a harmonic function parameter, and is used for describing periodic changes in the year; c is the intercept, used to describe the change in annual trends; solving harmonic function parameters by adopting a general least square method;
step 6.2, acquiring vegetation growth period climate parameters SOS from satellite images satellite,x1 Defining the corresponding time of the rise of the harmonic function curve from the minimum value to the maximum value by x1%, and acquiring a vegetation degradation period weathered parameter EOS obtained by the satellite remote sensing image satellite,x2 Defining the corresponding time of the decrease of the harmonic function curve from the maximum value to the minimum value by x 2%;
and 6.3, taking x1 and x2 as undetermined thresholds, and trying with the step length of 5 from 20 to 80, wherein when the satellite image vegetation weather is closest to the ground observation vegetation weather, the x1 and x2 are the optimal thresholds, and the formula is as follows:
|SOS satellite,x1 -SOS ground |→min,x1＝20,25,…,80
|EOS satellite,x2 -EOS ground |→min,x2＝20,25,…,80。
7. a satellite image calibration device for a coastal zone vegetation weather parameter of a joint ground observation, which is used for executing the satellite image calibration method for the coastal zone vegetation weather parameter of the joint ground observation according to any one of claims 1 to 6, and comprises the following steps:
the acquisition module is used for acquiring a data source, wherein the data source comprises a ground camera photo from the same year and a satellite remote sensing image;
the first construction module is used for carrying out data preprocessing on the photos of the ground camera, and comprises the steps of calculating a green coordinate index GCC and denoising a percentile value, so as to construct a photo GCC time sequence observation;
the first fitting module is used for fitting a photo GCC time sequence observation curve by using a harmonic function and extracting the time corresponding to the inflection point of the curve as a real earth surface weather parameter;
the correction module is used for combining the visibility analysis and GPS correction to determine the corresponding view range of the ground camera photo;
the second construction module is used for carrying out data preprocessing on the satellite remote sensing image, and comprises normalized difference vegetation index NDVI calculation and inverse distance weighted synthesis, and constructing an image NDVI time sequence observation in a corresponding view range;
the second fitting module is used for fitting an image NDVI time sequence observation curve by using a harmonic function, acquiring image weather parameters through different thresholds, and selecting the threshold closest to the ground surface real weather parameters as an optimal threshold;
and the calculation module is used for applying an optimal threshold value to calculate each pixel climatic parameter of the image of the coverage research area so as to obtain the climatic space distribution.
8. A computer storage medium, wherein a computer program is stored in the computer storage medium; the computer program, when run on a computer, causes the computer to perform the satellite image calibration method for the vegetation weatherometer parameters of the coastal zone in combination with the ground observation as claimed in any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310439447.8A CN116523859A (en) | 2023-04-23 | 2023-04-23 | Coastal zone vegetation weather parameter satellite image calibration method combined with ground observation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310439447.8A CN116523859A (en) | 2023-04-23 | 2023-04-23 | Coastal zone vegetation weather parameter satellite image calibration method combined with ground observation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116523859A true CN116523859A (en) | 2023-08-01 |
Family
ID=87391463
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310439447.8A Pending CN116523859A (en) | 2023-04-23 | 2023-04-23 | Coastal zone vegetation weather parameter satellite image calibration method combined with ground observation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116523859A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117457066A (en) * | 2023-12-26 | 2024-01-26 | 山东科技大学 | Winter wheat grain protein content prediction method with provincial scale |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107632967A (en) * | 2017-09-14 | 2018-01-26 | 青海省基础地理信息中心 | A kind of meadow grass yield evaluation method |
| CN111579565A (en) * | 2019-02-18 | 2020-08-25 | 深圳先进技术研究院 | Agricultural drought monitoring method, system and storage medium |
| CN112016052A (en) * | 2020-08-20 | 2020-12-01 | 广东省气象探测数据中心 | Near-surface daily maximum air temperature estimation method, system and terminal based on multi-source data |
| CN112446323A (en) * | 2020-11-24 | 2021-03-05 | 云南电网有限责任公司电力科学研究院 | HHT harmonic analysis method based on improved EMD modal aliasing and endpoint effect |
| CN113850139A (en) * | 2021-08-25 | 2021-12-28 | 南京林业大学 | Multi-source remote sensing-based forest annual phenological monitoring method |
| CN115019190A (en) * | 2022-04-11 | 2022-09-06 | 贵州师范大学 | Terrain broken region complex terrain information extraction method based on aerial remote sensing platform |
-
2023
- 2023-04-23 CN CN202310439447.8A patent/CN116523859A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107632967A (en) * | 2017-09-14 | 2018-01-26 | 青海省基础地理信息中心 | A kind of meadow grass yield evaluation method |
| CN111579565A (en) * | 2019-02-18 | 2020-08-25 | 深圳先进技术研究院 | Agricultural drought monitoring method, system and storage medium |
| CN112016052A (en) * | 2020-08-20 | 2020-12-01 | 广东省气象探测数据中心 | Near-surface daily maximum air temperature estimation method, system and terminal based on multi-source data |
| CN112446323A (en) * | 2020-11-24 | 2021-03-05 | 云南电网有限责任公司电力科学研究院 | HHT harmonic analysis method based on improved EMD modal aliasing and endpoint effect |
| CN113850139A (en) * | 2021-08-25 | 2021-12-28 | 南京林业大学 | Multi-source remote sensing-based forest annual phenological monitoring method |
| CN115019190A (en) * | 2022-04-11 | 2022-09-06 | 贵州师范大学 | Terrain broken region complex terrain information extraction method based on aerial remote sensing platform |
Non-Patent Citations (2)
| Title |
|---|
| 范德芹;赵学胜;朱文泉;郑周涛;: "植物物候遥感监测精度影响因素研究综述", 地理科学进展, no. 03, 31 March 2016 (2016-03-31), pages 304 - 319 * |
| 赵君;张晓民;: "改进的TOPKAPI模型及其在洪水预报中的应用", 河海大学学报(自然科学版), no. 02, 25 March 2011 (2011-03-25), pages 131 - 136 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117457066A (en) * | 2023-12-26 | 2024-01-26 | 山东科技大学 | Winter wheat grain protein content prediction method with provincial scale |
| CN117457066B (en) * | 2023-12-26 | 2024-03-15 | 山东科技大学 | Winter wheat grain protein content prediction method with provincial scale |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | An index and approach for water extraction using Landsat–OLI data | |
| CN111242224B (en) | Multi-source remote sensing data classification method based on unmanned aerial vehicle extraction classification sample points | |
| CN112348812B (en) | Forest stand age information measurement method and device | |
| WO2015195484A1 (en) | Automated geospatial image mosaic generation with multiple zoom level support | |
| CN113033670A (en) | Method for extracting rice planting area based on Sentinel-2A/B data | |
| CN109308688B (en) | Visible light and near-infrared band thick cloud and shadow removing method | |
| Ebrahimy et al. | Downscaling MODIS land surface temperature product using an adaptive random forest regression method and Google Earth Engine for a 19-years spatiotemporal trend analysis over Iran | |
| CN107688777B (en) | Urban green land extraction method for collaborative multi-source remote sensing image | |
| CN114563378B (en) | Method, device, medium and equipment for quantitatively describing space distribution of cyanobacterial bloom in lakes and reservoirs | |
| CN114740180A (en) | Soil organic carbon estimation method and device based on multi-source remote sensing data | |
| CN113744249B (en) | Marine ecological environment damage investigation method | |
| Chen et al. | Mapping double-cropped irrigated rice fields in Taiwan using time-series Satellite Pour I'Observation de la Terre data | |
| CN116523859A (en) | Coastal zone vegetation weather parameter satellite image calibration method combined with ground observation | |
| Kennaway et al. | Mapping land cover and estimating forest structure using satellite imagery and coarse resolution lidar in the Virgin Islands | |
| Xue et al. | Flood monitoring by integrating normalized difference flood index and probability distribution of water bodies | |
| Singh et al. | Estimation of revised capacity in Gobind Sagar reservoir using Google earth engine and GIS | |
| Wang et al. | Virtual image-based cloud removal for Landsat images | |
| CN112052720B (en) | High-space-time normalization vegetation index NDVI fusion model based on histogram clustering | |
| CN115203624A (en) | Time sequence remote sensing-based comprehensive evaluation method for earth surface environment at any moment | |
| Xie et al. | Estimation of forest stand parameters using SPOT-5 satellite images and topographic information | |
| Salman et al. | Detection of Spectral Reflective Changes for Temporal Resolution of Land Cover (LC) for Two Different Seasons in central Iraq | |
| El-Mezouar et al. | A high-resolution index for vegetation extraction in IKONOS images | |
| Wang et al. | A New Remote Sensing Change Detection Data Augmentation Method based on Mosaic Simulation and Haze Image Simulation | |
| CN104616262A (en) | H alpha cloud removing method based on quite sun chromospheres background processing | |
| Wang et al. | Virtual image patch-based cloud removal for Landsat images |
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 |