CN101825702B - Method for adjusting and optimizing terrain by utilizing terrain adjusting vegetation index - Google Patents
Method for adjusting and optimizing terrain by utilizing terrain adjusting vegetation index Download PDFInfo
- Publication number
- CN101825702B CN101825702B CN2010101808953A CN201010180895A CN101825702B CN 101825702 B CN101825702 B CN 101825702B CN 2010101808953 A CN2010101808953 A CN 2010101808953A CN 201010180895 A CN201010180895 A CN 201010180895A CN 101825702 B CN101825702 B CN 101825702B
- Authority
- CN
- China
- Prior art keywords
- tavi
- vegetation index
- landform
- terrain
- vegetation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a method for adjusting and optimizing a terrain by utilizing a terrain adjusting vegetation index. The method is characterized in that waveband information of an optical remote sensing image and derived information thereof are used for optimizing; and the method comprises the steps of constructing a terrain adjusting vegetation index (TAVI), classifying images and optimizing matching, can accurately and quickly invert the vegetation index of complicated terrain mountain areas, provide an important means for drawing complicated terrain mountain area vegetation (forest), monitoring vegetation coverage, estimating forest coverage rate and inverting the leaf area index, can also provide important calculation parameters for estimating biophysical and biochemical parameters such as forest biomass, productivity, photosynthetic active radiation and absorption and the like, and has favorable practical value.
Description
Technical field
The present invention relates to a kind of landform and regulate the landform adjusting optimization method of vegetation index.
Background technology
The complex-terrain mountain area; Because topographic relief has changed the distribution of solar radiation on the face of land, for same wave band, the radiance value of massif Schattenseite part is less than normal in the multispectral remote sensing image; And the radiance value of tailo part is bigger than normal; Cause image information to distort, the vegetation index information distortion by corresponding wave band generates has had a strong impact on the accurate inverting of atural object vegetation information.The influence of topography becomes one of major obstacle of the accurate inverting of mountain area vegetation information remote sensing.
At present, conventional disposal route is before calculating vegetation index, to carry out topographic correction earlier.The method of topographic correction can be divided into two big types, and one type of topographic correction that is based on dem data mainly contains: (1) empirical statistics model, classic algorithm have cosine correction, C to proofread and correct, SCS proofreaies and correct, SCS+C proofreaies and correct, Proy proofreaies and correct etc.; (2) based on the physical model of radiation transfer theory, carry out topographic correction through the physical process of research light and face of land effect; (3) based on the landform radiant correction of statistical study, like the square matching algorithm.These methods all produce effect, but degree varies is also different to the improvement degree of vegetation information; Even more serious problem is that these conventional landform correcting method of lane all need the support of high accuracy DEM data.Because the introducing of dem data might be introduced new error; The renewal of dem data often lags behind the variation that actual landform brings because of nature, human factor again; Particularly the high accuracy DEM Data Acquisition is restricted, and has restricted the effect and the application of this type of landform correcting method of lane to a great extent.Another kind of landform correcting method of lane need not the support of dem data, and as adopting the method for wave band than model, because that landform changes the influence of different-waveband is different, so this model can only play the effect of constraining the influence of topography, and do not reach desirable topographic correction effect; The linear matching process of gray scale be because the sample of selecting can not be contained all images, and the correction of therefore only carrying out wave band with a linear function of this derivation has bigger limitation; Do not see correlative study and can improve mountain area vegetation index information inversion accuracy to the method for cloud, buildings shadow removal in the high-resolution remote sensing image of city.
To the deficiency of above-mentioned the whole bag of tricks, I invent multi landform regulate the building method (national patent application number 200910111688.X) of vegetation index (having another name called topographic correction vegetation index).The vegetation index that makes up according to this method need not the support of data such as DEM, just can effectively eliminate the influence of landform to vegetation information.For the key factor in this vegetation index---the optimized Algorithm of landform regulatory factor (having another name called the topographic correction factor) has proposed a kind of optimized Algorithm based on the ground investigation data among " building method of topographic correction vegetation index " (national patent application number 200910111688.X); Because this optimized Algorithm needs the support of ground investigation data or field study data, receives certain restriction in the practical application.For this reason; Design a kind of landform that does not rely on ground investigation data or field study data and regulate optimized Algorithm, landform is regulated vegetation index have important scientific meaning and economic worth with the businessization popularization in the successful Application of the accurate inverting of complex-terrain mountain area vegetation information.
Summary of the invention
The purpose of this invention is to provide a kind of landform and regulate the landform adjusting optimization method of vegetation index; It is inverting complex-terrain mountain area vegetation information accurately and rapidly; For complex-terrain mountain area vegetation (forest) drawing, vegetation coverage monitoring, afforestation rate assessment, leaf area index inverting provide important means, and can the important computations parameter be provided for the estimation of biophysics such as forest biomass, yield-power, photosynthetically active radiation absorption and biochemical parameters.
Invention adopts following technical scheme to realize: a kind of landform is regulated the landform of vegetation index and is regulated optimization method, it is characterized in that: utilize the band class information of optical remote sensing image self and derived information thereof to be optimized, may further comprise the steps:
(1) make up landform and regulate vegetation index TAVI, its computing formula is: TAVI=CVI+TAC * SVI; Wherein, CVI representes vegetation index commonly used, and SVI representes the shade vegetation index; TAC representes the landform regulatory factor, and the corresponding calculated formula is:
SVI=[MAX (B
r) – B
r]/B
rWherein, B
rRepresent red wave band data, MAX (B
r) maximal value of expression study area red wave band data.
(2) select the appearance district: inspection remote sensing image quality, select a certain size planar zone in the complex-terrain mountain area, guarantee appearance district's image " noise " disturb minimum, to have the strong influence of topography, an area enough big;
(3) image classification: use conventional not supervised classification and be divided into Schattenseite and tailo two big classes to the vegetation in the appearance district remote sensing image;
(4) optimization of matching: the design cycle program, make TAC since 0,0.001 to be the interval, increase progressively successively, investigate the maximal value M of Schattenseite part TAVI simultaneously
TAVI is cloudyMaximal value M with tailo part TAVI
The TAVI sun, work as M
TAVI is cloudyWith M
The TAVI sunEquate or during approximately equal, withdraw from circulation, finally confirm the TAC Optimization result; If when TAC tires out when increasing to 5 M
TAVI is cloudyWith M
The TAVI sunAlso do not satisfy condition, then return step (2), reselect the appearance district, repeating step (2) is to step (4), until M
TAVI is cloudyWith M
The TAVI sunSatisfy and equate or close condition.
The present invention has following advantage:
The first, the topographic correction effect is remarkable.The landform regulatory factor that this optimization method is confirmed, the assurance landform is regulated vegetation index can effectively eliminate the interference of the influence of topography to vegetation information.Regulate the linear regression analysis of vegetation index and solar incident angle cosine value through the landform of this method being confirmed in experiment appearance district; Show that landform that this method is confirmed regulates the related coefficient and the equation of linear regression slope of vegetation index and solar incident angle cosine value and can reduce to minimum; Obviously be superior to other vegetation index commonly used, and more excellent than common atmospheric correction or topographic correction effect.
The second, demand data is few, and cost is low.This optimization method only needs the self-contained wave band data of remote sensing image, need not the support of ground investigation data or on-the-spot investigation data etc., and data cost and time cost are realized minimizing.
The 3rd, flow process is simple, and is workable.This optimization method mainly just can be confirmed the Optimization result of landform regulatory factor by " selecting the appearance district ", " image classification " and " optimization of matching " three steps, flow process is simple, and operation easily helps landform to regulate the business application and the popularization of vegetation index.
Description of drawings
Fig. 1 is a techniqueflow synoptic diagram of the present invention.
Embodiment
As shown in Figure 1, the landform that the present invention provides a kind of landform to regulate vegetation index is regulated optimization method, it is characterized in that: utilize the band class information of optical remote sensing image self and derived information thereof to be optimized, may further comprise the steps:
(1) make up landform and regulate vegetation index TAVI: the basic calculating formula that landform is regulated vegetation index is following
TAVI?=?CVI?+?TAC?*?SVI
Wherein, CVI representes vegetation index commonly used, like normalized differential vegetation index NDVI, ratio vegetation index RVI etc.; SVI representes the shade vegetation index; TAC representes the landform regulatory factor, and corresponding concrete computing formula is following:
NDVI?=?(B
nir?-?B
r?)?/?(B
nir?+?B
r?)
RVI?=?B
nir?/?B
r
SVI?=?[MAX(B
r)?–?B
r]?/?B
r
Wherein, B
NirExpression near-infrared band data, B
rRepresent red wave band data, MAX (B
r) represent that the maximal value of the red wave band data of study area (is noted that: in present embodiment here; The remote sensing image wave band data all with apparent reflectance as experimental data; Do not specify that if having it all is the basis with the apparent reflectance that vegetation index hereinafter calculates).
(2) select the appearance district: inspection remote sensing image quality; Select a certain size planar zone in the complex-terrain mountain area, guarantee that appearance district's image " noise " disturbs minimum, has the strong influence of topography, area enough big (preferable area is greater than 10 square kilometres) here;
(3) image classification: use conventional not supervised classification and be divided into Schattenseite and tailo two big classes to the vegetation in the appearance district remote sensing image, like maximum likelihood method or fuzzy clustering method etc.;
(4) optimization of matching: according to formula TAVI=CVI+TAC * SVI, design TAVI optimizes loop program: make TAC since 0,0.001 to be the interval, increase progressively successively, investigate the maximal value M of Schattenseite part TAVI simultaneously
TAVI is cloudyMaximal value M with tailo part TAVI
The TAVI sun, work as M
TAVI is cloudyWith M
The TAVI sunEquate or during approximately equal, withdraw from circulation, finally confirm the TAC Optimization result; If when TAC tires out when increasing to 5 M
TAVI is cloudyWith M
The TAVI sunAlso do not satisfy condition, then return step (2), reselect the appearance district, repeating step (2) is to step (4), until M
TAVI is cloudyWith M
The TAVI sunSatisfy and equate or close condition.
The above is merely preferred embodiment of the present invention, and all equalizations of doing according to claim of the present invention change and modify, and all should belong to covering scope of the present invention.
Claims (3)
1. the landform of a landform adjusting vegetation index is regulated optimization method, it is characterized in that: utilize the band class information of optical remote sensing image self and derived information thereof to be optimized, may further comprise the steps:
(1) make up landform and regulate vegetation index TAVI, its computing formula is: TAVI=CVI+TAC * SVI; Wherein, CVI representes vegetation index commonly used, and SVI representes the shade vegetation index; TAC representes the landform regulatory factor, and the corresponding calculated formula is:
SVI=[MAX (B
r) – B
r]/B
rWherein, B
rRepresent red wave band data, MAX (B
r) maximal value of expression study area red wave band data;
(2) select the appearance district: inspection remote sensing image quality, select a certain size planar zone in the complex-terrain mountain area, guarantee appearance district's image " noise " disturb minimum, to have the strong influence of topography, an area enough big;
(3) image classification: use conventional not supervised classification and be divided into Schattenseite and tailo two big classes to the vegetation in the appearance district remote sensing image;
(4) optimization of matching: the design cycle program, make TAC since 0,0.001 to be the interval, increase progressively successively, investigate the maximal value M of Schattenseite part TAVI simultaneously
TAVI is cloudyMaximal value M with tailo part TAVI
The TAVI sun, work as M
TAVI is cloudyWith M
The TAVI sunWhen equating, withdraw from circulation, finally confirm the TAC Optimization result; If when TAC tires out when increasing to 5 M
TAVI is cloudyWith M
The TAVI sunAlso do not satisfy condition, then return step (2), reselect the appearance district, repeating step (2) is to step (4), until M
TAVI is cloudyWith M
The TAVI sunSatisfy the condition that equates.
2. landform according to claim 1 is regulated the landform of vegetation index and is regulated optimization method, and it is characterized in that: the vegetation index CVI commonly used in the said step (1) comprises normalized differential vegetation index NDVI, ratio vegetation index RVI; The corresponding calculated formula is:
NDVI=(B
Nir-B
r)/(B
Nir+ B
r); RVI=B
Nir/ B
rWherein, B
NirExpression near-infrared band data.
3. landform according to claim 1 is regulated the landform of vegetation index and regulated optimization method, it is characterized in that: the appearance district area of said step (2) is greater than 10 square kilometres.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010101808953A CN101825702B (en) | 2010-05-24 | 2010-05-24 | Method for adjusting and optimizing terrain by utilizing terrain adjusting vegetation index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010101808953A CN101825702B (en) | 2010-05-24 | 2010-05-24 | Method for adjusting and optimizing terrain by utilizing terrain adjusting vegetation index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101825702A CN101825702A (en) | 2010-09-08 |
| CN101825702B true CN101825702B (en) | 2012-11-21 |
Family
ID=42689730
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010101808953A Expired - Fee Related CN101825702B (en) | 2010-05-24 | 2010-05-24 | Method for adjusting and optimizing terrain by utilizing terrain adjusting vegetation index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101825702B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105389538B (en) * | 2015-10-09 | 2018-07-13 | 南京大学 | A method of based on a cloud hemisphere slice estimation Forest Leaf Area Index |
| CN105487066B (en) * | 2015-11-20 | 2018-02-09 | 福州大学 | A kind of new TAVI landform regulatory factor optimized algorithm |
| CN106324614B (en) * | 2016-08-10 | 2018-08-17 | 福州大学 | A kind of new TAVI combinational algorithms |
| WO2018028191A1 (en) | 2016-08-10 | 2018-02-15 | 福州大学 | Tavi calculation method based on waveband ration model and solar elevation angle |
| CN107730502A (en) * | 2017-11-15 | 2018-02-23 | 福州大学 | A kind of algorithm on the s factors in new TAVI models |
| CN109086661B (en) * | 2018-06-14 | 2019-05-03 | 中科禾信遥感科技（苏州）有限公司 | A kind of crops relative radiometric normalization method and device |
| CN113553549B (en) * | 2021-07-26 | 2023-04-14 | 中国科学院西北生态环境资源研究院 | Method and device for inversion of coverage degree of planting, electronic equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101561502A (en) * | 2009-05-07 | 2009-10-21 | 福州大学 | Constructing method for topographic correction vegetation index |
| CN101699315A (en) * | 2009-10-23 | 2010-04-28 | 北京农业信息技术研究中心 | Monitoring device and method for crop growth uniformity |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008079549A (en) * | 2006-09-28 | 2008-04-10 | Mitsubishi Paper Mills Ltd | Method for evaluating tree growth |
-
2010
- 2010-05-24 CN CN2010101808953A patent/CN101825702B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101561502A (en) * | 2009-05-07 | 2009-10-21 | 福州大学 | Constructing method for topographic correction vegetation index |
| CN101699315A (en) * | 2009-10-23 | 2010-04-28 | 北京农业信息技术研究中心 | Monitoring device and method for crop growth uniformity |
Non-Patent Citations (7)
| Title |
|---|
| 汪小钦.福建省长汀县植被覆盖度遥感动态监测研究.《自然资源学报》.2006,第21卷(第1期), |
| 温仲明.基于TM遥感影像的陕北黄土区结构化植被因子指数提取.《应用生态学报》.2009,第20卷(第11期), |
| 王钦敏 |
| 江洪 |
| 江洪; 王钦敏; 汪小钦.福建省长汀县植被覆盖度遥感动态监测研究.《自然资源学报》.2006,第21卷(第1期), * |
| 雷婉宁 |
| 雷婉宁; 温仲明.基于TM遥感影像的陕北黄土区结构化植被因子指数提取.《应用生态学报》.2009,第20卷(第11期), * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101825702A (en) | 2010-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101825702B (en) | Method for adjusting and optimizing terrain by utilizing terrain adjusting vegetation index | |
| CN101561502B (en) | Constructing method for topographic correction vegetation index | |
| Chirici et al. | Non-parametric and parametric methods using satellite images for estimating growing stock volume in alpine and Mediterranean forest ecosystems | |
| Sehgal et al. | Inversion of radiative transfer model for retrieval of wheat biophysical parameters from broadband reflectance measurements | |
| Tonolli et al. | Fusion of airborne LiDAR and satellite multispectral data for the estimation of timber volume in the Southern Alps | |
| Xu et al. | Estimating soil total nitrogen in smallholder farm settings using remote sensing spectral indices and regression kriging | |
| CN103294792B (en) | Based on the polarization SAR terrain classification method of semantic information and polarization decomposing | |
| CN114021348A (en) | Mining area vegetation carbon sink remote sensing inversion method of fine land utilization type | |
| CN103345707A (en) | Crop maturation stage remote sensing prediction method based on multi-source remote sensing data | |
| Browning et al. | Digital soil mapping in the absence of field training data: A case study using terrain attributes and semiautomated soil signature derivation to distinguish ecological potential | |
| CN111582575B (en) | Method for identifying urban thermal environment formation development leading factors under multiple space-time scales | |
| Torbick et al. | Mapping rice greenhouse gas emissions in the Red River Delta, Vietnam | |
| CN102636779A (en) | Extraction method for coverage rate of sub-pixel accumulated snow based on resampling regression analysis | |
| CN106324614B (en) | A kind of new TAVI combinational algorithms | |
| Xiao et al. | Bias correction framework for satellite precipitation products using a rain/no rain discriminative model | |
| Zhou et al. | Spatiotemporal fusion of multi-source remote sensing data for estimating aboveground biomass of grassland | |
| CN108205718A (en) | Production method and system are surveyed in a kind of cereal crops sampling | |
| Villoslada et al. | Reindeer control over shrubification in subarctic wetlands: spatial analysis based on unoccupied aerial vehicle imagery | |
| Zhu et al. | UAV flight height impacts on wheat biomass estimation via machine and deep learning | |
| Tian et al. | Retrieving leaf area index for coniferous forest in Xingguo County, China with Landsat ETM+ images | |
| Zou et al. | Combining spectral and texture feature of UAV image with plant height to improve LAI estimation of winter wheat at jointing stage | |
| CN105487066B (en) | A kind of new TAVI landform regulatory factor optimized algorithm | |
| CN109086661A (en) | A kind of crops relative radiometric normalization method and device | |
| Dahms et al. | Modelling biophysical parameters of maize using Landsat 8 time series | |
| Wijaya et al. | Fusion of ALOS Palsar and Landsat ETM data for land cover classification and biomass modeling using non-linear methods |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121121 Termination date: 20150524 |
|
| EXPY | Termination of patent right or utility model |