CN111105402A - SEVI (sequence independent variable) adjustment factor optimization method based on information entropy - Google Patents
SEVI (sequence independent variable) adjustment factor optimization method based on information entropy Download PDFInfo
- Publication number
- CN111105402A CN111105402A CN201911341317.0A CN201911341317A CN111105402A CN 111105402 A CN111105402 A CN 111105402A CN 201911341317 A CN201911341317 A CN 201911341317A CN 111105402 A CN111105402 A CN 111105402A
- Authority
- CN
- China
- Prior art keywords
- sevi
- adjustment factor
- information entropy
- calculating
- delta
- 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 30
- 238000005457 optimization Methods 0.000 title claims abstract description 14
- 238000011160 research Methods 0.000 claims abstract description 12
- 238000002310 reflectometry Methods 0.000 claims abstract description 10
- 230000008030 elimination Effects 0.000 claims abstract description 8
- 238000003379 elimination reaction Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000001105 regulatory effect Effects 0.000 claims description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
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
- 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)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention relates to an SEVI (sequence independent variable) adjustment factor optimization method based on information entropy, which comprises the steps of firstly selecting a research area according to a research target and obtaining a corresponding long-time sequence remote sensing image; calculating a vegetation index for eliminating the shade of the research area SEVI by using the reflectivity data: then selecting a pure vegetation complex terrain area, and calculating a Shannon information entropy H (v) by utilizing a shadow elimination vegetation index SEVI; starting the adjustment factor from 0, sequentially increasing the adjustment factor at preset intervals T, and respectively calculating SEVI and H (v); and when the value of H (v) is maximum, taking the adjustment factor corresponding to the maximum H (v) as the optimal adjustment factor. The method does not need DEM data and remote sensing image classification, and does not depend on ground survey data.
Description
Technical Field
The invention relates to the field of vegetation detection, in particular to an SEVI (sequence-independent vegetation index) regulating factor optimization method based on information entropy.
Background
The existing optimization methods for regulating factors in topographic shadow elimination vegetation index SEVI and early achievements TAVI and TCVI thereof mainly comprise 4 types: "matching method (national patent No. 200910111688X)", "extreme value (Max) method (national patent No. 201010180895.3)", "correlation coefficient (r) method (national patent No. 2015108077580)", and solar altitude method (national patent No. 201611127461.0).
The 4 optimization algorithms can effectively reduce the influence of terrain shadow on the vegetation information in the mountainous area without the support of data such as DEM and the like. However, the matching method and the extreme method need to classify the remote sensing images first, wherein the matching method also needs the support of ground data and the like; the correlation coefficient method is an empirical algorithm, and the theoretical basis is weak; the s factor in the solar elevation angle method is not easy to determine, and the s factor influences the large-scale popularization and application of the SEVI technology in long-time sequence monitoring of vegetation in mountainous areas and the like.
Disclosure of Invention
In view of this, the invention aims to provide an SEVI adjustment factor optimization method based on information entropy, which does not need DEM data and remote sensing image classification and does not depend on ground survey data.
The invention is realized by adopting the following scheme: an SEVI adjustment factor optimization method based on information entropy comprises the following steps:
selecting a research area according to a research target, and acquiring a corresponding long-time sequence remote sensing image;
calculating a vegetation index for eliminating the shade of the research area SEVI by using the reflectivity data:
wherein f (Delta) is a regulatory factor, BrFor red band reflectivity data of remote-sensing images, BnirThe data is the reflectivity data of the remote sensing image near infrared band;
selecting a pure vegetation complex terrain area (such as a mountainous area natural protection area with good vegetation protection), and calculating a Shannon information entropy H (v) by utilizing a shadow elimination vegetation index SEVI;
starting the adjustment factor f (delta) from 0, sequentially increasing the adjustment factor f (delta) at preset intervals T, and respectively calculating SEVI and H (v);
when the value of H (v) is maximum, taking the adjustment factor f (delta) corresponding to the maximum H (v) as the optimal adjustment factor fopt:
fopt=argmax(H(v))。
Further, selecting a pure vegetation complex terrain area, and calculating a shannon information entropy h (v) by using a shadow elimination vegetation index SEVI specifically comprises:
wherein,
in the formula, piIs the proportion of the i pixel information; n is the number of pixels in the region of interest,and f (delta) of the i pixel is equal to the SEVI calculation result when the value is equal to f.
Further, the preset interval T is 0.001.
compared with the prior art, the invention has the following beneficial effects:
1. the method uses the information entropy to calculate the adjustment factor, has stronger theoretical basis, gets rid of the defect that the original empirical algorithm is difficult to understand, and can be applied to monitoring and analyzing the vegetation change in the mountainous area for a long time sequence;
2. the optimal solution determined by the invention ensures that SEVI can effectively eliminate the interference of various terrain shadows on the vegetation information in the mountainous area;
3. the method only needs the waveband data carried by the remote sensing image, does not need the support of ground survey data or field investigation data and the like, and realizes the minimization of data cost and time cost;
4. the method has simple process and strong operability.
Drawings
FIG. 1 is a schematic flow chart of a method according to an embodiment of the present invention.
FIG. 2 is a diagram illustrating an f (Δ) optimization calculation curve according to an embodiment of the present invention.
Fig. 3 shows a long-term sequence topographic shadow map according to an embodiment of the present invention, (a) shows a map in 2011, (b) shows a map in 2012, (c) shows a map in 2013, (d) shows a map in 2014, (e) shows a map in 2015, (f) shows a map in 2016, (g) shows a map in 2017, and (h) shows a map in 2018.
Detailed Description
The invention is further explained below with reference to the drawings and the embodiments.
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
As shown in fig. 1, the present embodiment provides an SEVI adjustment factor optimization method based on information entropy, including the following steps:
selecting a research area according to a research target, and acquiring a corresponding long-time sequence remote sensing image;
calculating a vegetation index for eliminating the shade of the research area SEVI by using the reflectivity data:
wherein f (Delta) is a regulatory factor, BrFor red band reflectivity data of remote-sensing images, BnirThe data is the reflectivity data of the remote sensing image near infrared band;
selecting a pure vegetation complex terrain area (such as a mountainous area natural protection area with good vegetation protection), and calculating a Shannon information entropy H (v) by utilizing a shadow elimination vegetation index SEVI;
starting the adjustment factor f (delta) from 0, sequentially increasing the adjustment factor f (delta) at preset intervals T, and respectively calculating SEVI and H (v);
when the value of H (v) is maximum, taking the adjustment factor f (delta) corresponding to the maximum H (v) as the optimal adjustment factor fopt(as shown in fig. 2, the dot in fig. 2 is f (Δ) corresponding to the maximum h (v)):
fopt=argmax(H(v))。
in this embodiment, the selecting a pure vegetation complex terrain area, and calculating a shannon information entropy h (v) by using a shadow elimination vegetation index SEVI specifically includes:
wherein,
in the formula, piIs the proportion of the i pixel information; n is the number of pixels in the region of interest,and f (delta) of the i pixel is equal to the SEVI calculation result when the value is equal to f.
In this embodiment, the preset interval T is 0.001.
the SEVI calculation result of the embodiment is shown in fig. 3, and application and verification of the method to the multispectral remote sensing image in the 8 th period of 2018 in the natural protection area 2011-2018 in the wuyishan mountain of china show that the average errors of the SEVI in the 8 th period calculated by the method in the terrain ghost and the shadow relative to the non-shadow sunny slope are only 0.93% and-1.94% respectively, the average value of the determination coefficient of the SEVI and the cosine value of the sun incident angle (cosi) correlation analysis is only 0.005, the sample variation coefficient is less than 9.69%, and a good terrain shadow correction effect of the remote sensing image of the long-time sequence is obtained.
The foregoing is directed to preferred embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow. However, any simple modification, equivalent change and modification of the above embodiments according to the technical essence of the present invention are within the protection scope of the technical solution of the present invention.
Claims (4)
1. An SEVI regulation factor optimization method based on information entropy is characterized by comprising the following steps:
selecting a research area according to a research target, and acquiring a corresponding long-time sequence remote sensing image;
calculating a vegetation index for eliminating the shade of the research area SEVI by using the reflectivity data:
wherein f (Delta) is a regulatory factor, BrFor red band reflectivity data of remote-sensing images, BnirThe data is the reflectivity data of the remote sensing image near infrared band;
selecting a pure vegetation complex terrain area, and calculating a Shannon information entropy H (v) by utilizing a shadow elimination vegetation index SEVI;
starting the adjustment factor f (delta) from 0, sequentially increasing the adjustment factor f (delta) at preset intervals T, and respectively calculating SEVI and H (v);
and when the value of H (v) is maximum, taking the adjustment factor f (delta) corresponding to the maximum H (v) as the optimal adjustment factor.
2. The SEVI regulation factor optimization method based on information entropy according to claim 1, wherein the method for selecting the pure vegetation complex terrain area and calculating the Shannon information entropy H (v) by using the shadow elimination vegetation index SEVI specifically comprises the following steps:
wherein,
3. An information entropy based SEVI adjustment factor optimization method according to claim 1, wherein the preset interval T is 0.001.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911341317.0A CN111105402A (en) | 2019-12-24 | 2019-12-24 | SEVI (sequence independent variable) adjustment factor optimization method based on information entropy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911341317.0A CN111105402A (en) | 2019-12-24 | 2019-12-24 | SEVI (sequence independent variable) adjustment factor optimization method based on information entropy |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111105402A true CN111105402A (en) | 2020-05-05 |
Family
ID=70423937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911341317.0A Pending CN111105402A (en) | 2019-12-24 | 2019-12-24 | SEVI (sequence independent variable) adjustment factor optimization method based on information entropy |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111105402A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112964643A (en) * | 2021-02-03 | 2021-06-15 | 福州大学 | Method for correcting landform falling shadow of visible light wave band of remote sensing image |
CN114332645A (en) * | 2021-12-31 | 2022-04-12 | 福州大学 | SEVI (gradient information entropy) regulation factor optimization method based on steep slope block information entropy |
CN114778483A (en) * | 2022-04-25 | 2022-07-22 | 福州大学 | Method for correcting terrain shadow of remote sensing image near-infrared wave band for monitoring mountainous region |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006116769A1 (en) * | 2005-04-28 | 2006-11-02 | The General Hospital Corporation | Systems, processes and software arrangements for evaluating information associated with an anatomical structure by an optical coherence ranging technique |
CN108304798A (en) * | 2018-01-30 | 2018-07-20 | 北京同方软件股份有限公司 | The event video detecting method of order in the street based on deep learning and Movement consistency |
CN109471125A (en) * | 2018-10-17 | 2019-03-15 | 福州大学 | A kind of SEVI regulatory factor method of global optimization |
CN109638810A (en) * | 2018-11-02 | 2019-04-16 | 中国电力科学研究院有限公司 | A kind of energy storage method and system for planning based on electric power system transient stability |
-
2019
- 2019-12-24 CN CN201911341317.0A patent/CN111105402A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006116769A1 (en) * | 2005-04-28 | 2006-11-02 | The General Hospital Corporation | Systems, processes and software arrangements for evaluating information associated with an anatomical structure by an optical coherence ranging technique |
CN108304798A (en) * | 2018-01-30 | 2018-07-20 | 北京同方软件股份有限公司 | The event video detecting method of order in the street based on deep learning and Movement consistency |
CN109471125A (en) * | 2018-10-17 | 2019-03-15 | 福州大学 | A kind of SEVI regulatory factor method of global optimization |
CN109638810A (en) * | 2018-11-02 | 2019-04-16 | 中国电力科学研究院有限公司 | A kind of energy storage method and system for planning based on electric power system transient stability |
Non-Patent Citations (1)
Title |
---|
江洪: "基于森林碳汇信息的福建省低碳发展研究", 《中国优秀博硕士学位论文全文数据库(博士)经济与管理科学辑 J145-4》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112964643A (en) * | 2021-02-03 | 2021-06-15 | 福州大学 | Method for correcting landform falling shadow of visible light wave band of remote sensing image |
CN114332645A (en) * | 2021-12-31 | 2022-04-12 | 福州大学 | SEVI (gradient information entropy) regulation factor optimization method based on steep slope block information entropy |
CN114332645B (en) * | 2021-12-31 | 2024-06-07 | 福州大学 | SEVI (sea-level-difference-like elevation) adjusting factor optimization method based on abrupt slope block information entropy |
CN114778483A (en) * | 2022-04-25 | 2022-07-22 | 福州大学 | Method for correcting terrain shadow of remote sensing image near-infrared wave band for monitoring mountainous region |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108662991B (en) | Land block scale winter wheat leaf area index estimation method based on remote sensing satellite data | |
US20120155714A1 (en) | Vegetation indices for measuring multilayer microcrop density and growth | |
CN111105402A (en) | SEVI (sequence independent variable) adjustment factor optimization method based on information entropy | |
Gitelson | Wide dynamic range vegetation index for remote quantification of biophysical characteristics of vegetation | |
CN106372592B (en) | A kind of winter wheat planting area calculation method based on winter wheat area index | |
AU2020101709A4 (en) | Crop yield prediction method and system based on low-altitude remote sensing information from unmanned aerial vehicle | |
US10585210B2 (en) | Apparatus for radiometric correction and orthorectification of aerial imagery | |
WO2014103181A1 (en) | Image measuring method, system, device, and program | |
CN102313526A (en) | Method for obtaining leaf area index based on quantitative fusion and inversion of multi-angle and multi-spectral remote sensing data | |
CN103278503B (en) | Multi-sensor technology-based grape water stress diagnosis method and system therefor | |
CN109974854B (en) | Radiation correction method for frame-type FPI (field programmable Gate array) hyperspectral image | |
CN116519557B (en) | Aerosol optical thickness inversion method | |
CN109784300A (en) | A kind of crops science survey production method and system | |
Yamasaki et al. | Development of a high-accuracy autonomous sensing system for a field scouting robot | |
CN109031343B (en) | Automatic optimization algorithm for SEVI (sequence independent optimization) adjustment factors of window traversal | |
Shibayama et al. | Continuous monitoring of visible and near-infrared band reflectance from a rice paddy for determining nitrogen uptake using digital cameras | |
CN107529726A (en) | Check device, inspection method and program | |
CN107657206A (en) | A kind of method based on remote sensing technology estimation forest coverage rate | |
Liu et al. | Application of UAV-retrieved canopy spectra for remote evaluation of rice full heading date | |
CN114332645B (en) | SEVI (sea-level-difference-like elevation) adjusting factor optimization method based on abrupt slope block information entropy | |
CN113474635A (en) | Image processing apparatus, image processing method, and program | |
Shang et al. | Estimation of crop ground cover and leaf area index (LAI) of wheat using RapidEye satellite data: Prelimary study | |
Meng et al. | Hot dark spot index method based on multi-angular remote sensing for leaf area index retrieval | |
Li et al. | Optimal color composition method for generating high-quality daily photographic time series from PhenoCam | |
CN109471125B (en) | Globally optimized SEVI (sequence independent variable) adjustment factor method |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200505 |
|
RJ01 | Rejection of invention patent application after publication |