CN109710718A - A kind of Method of fast estimating of the vegetative coverage factor - Google Patents
A kind of Method of fast estimating of the vegetative coverage factor Download PDFInfo
- Publication number
- CN109710718A CN109710718A CN201811598780.9A CN201811598780A CN109710718A CN 109710718 A CN109710718 A CN 109710718A CN 201811598780 A CN201811598780 A CN 201811598780A CN 109710718 A CN109710718 A CN 109710718A
- Authority
- CN
- China
- Prior art keywords
- use pattern
- land
- vegetative coverage
- coverage factor
- land use
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Image Processing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of Method of fast estimating of vegetative coverage factor.Method first determines land use pattern grid stability, obtains stable land use pattern grid and non-stable land use pattern grid.Vegetative coverage factor grid value range on stable land use pattern grid is determined secondly, analyzing using meta.Again, the dynamic analog of the vegetative coverage factor on non-stable land use pattern grid is carried out using linear model.Finally, merging to obtain complete vegetative coverage factor space distribution map by space.Vegetative coverage factor Characteristics of spatio-temporal in the quantitative analysis of the present invention land use pattern of same position different time, vegetative coverage factor variations in dynamic analog land use pattern conversion process, this drawback of the dynamic change of the vegetative coverage factor in process of land use change cannot be estimated by overcoming conventional method, provide new approaches for soil erosion simulation in big regional scale.
Description
Technical field
The present invention relates to soil erosion method and technology fields, and in particular to a kind of quick estimation side of the vegetative coverage factor
Method.
Background technique
In soil erosion appraising model, the vegetative coverage factor refers in the certain time of the same terms, is there is vegetation
Soil loss amount and clear the ratio between the soil loss amount ploughed on the soil continuously to lie fallow on the soil of covering.The vegetative coverage factor is anti-
Vegetative coverage has been answered to act on (Wischmeier etc., 1965) to the inhibition of the soil erosion.Under normal conditions, the vegetative coverage factor
Estimation needs the information such as plant canopy, ground mulching degree, roughness of ground surface, soil moisture content, and the vegetative coverage factor is estimated
Calculation value can only also represent the analog result on point scale or normal cell scale, and the estimation for big regional scale can not be real
It is existing, it is therefore desirable to seek new method and establish the simple algorithm for being suitable for the vegetative coverage factor of big regional scale.
The method for studying the big regional scale vegetative coverage factor both at home and abroad at present is mainly land use pattern indirect assignment
Method.Empirical value assignment is carried out to the vegetative coverage factor in different land use type.Land use pattern indirect assignment method
Think that the vegetative coverage factor in same land use pattern is identical and is a constant constant.However, in land use pattern
During changing, vegetative coverage factor values are gradually changed.Such as during forest land is changed into meadow,
Vegetation coverage constantly successively decreases from high to lower until reaching a stable state, during this transformation, vegetative coverage because
Subvalue gradually increases with the reduction of vegetation coverage.Traditional assignment method has ignored same place different time land use pattern
Vegetative coverage factor dynamic variation feature in change procedure, therefore there is errors.
Land use pattern and vegetation coverage information are comprehensively utilized, the mould between the two and the vegetative coverage factor is established
Type relationship facilitates the fast quantification estimation and drawing of the big regional scale vegetative coverage factor.
Summary of the invention
It is an object of the invention to solve the problems, such as existing method, and provide a kind of the quick of the vegetative coverage factor
Evaluation method.
The step of the technical solution adopted by the present invention, is as follows:
A kind of Method of fast estimating of the vegetative coverage factor comprising following steps:
Step (1) data acquisition: according to time span range to be evaluated, the soil benefit within the scope of the time span is obtained
With type raster data and MODIS vegetation raster data；
Step (2) land use pattern grid stability preliminary judgement: the land use pattern grid that step (1) is obtained
Data carry out multidate land use pattern grid stability preliminary judgement, obtain land use pattern grid steady in a long-term and
Long-term non-stable land use pattern grid；
Step (3) land use pattern grid stability further determines: the long-term non-stable soil that step (2) is obtained
Ground use pattern grid further determines, obtain short-term stability land use pattern grid and short-term non-stable land use
Type grid；
Step (4) land use pattern grid conversion identification: short-term non-stable soil benefit obtained in analytical procedure (3)
With type grid, the transition form of land use pattern grid in such land use pattern grid is identified and tracked；
The vegetative coverage factor assignment in step (5) arable land: to land use pattern steady in a long-term obtained in step (2)
In arable land grid and step (3) obtained in short-term stability land use pattern in arable land grid directly assigned
Value；
Step (6) carries out assignment to the vegetative coverage factor value range of the land use pattern other than arable land: according to step
Suddenly the land use pattern raster data that (1) obtains, determine the vegetative coverage of different land use type in addition to arable land because
The grid value range of son；
The vegetative coverage factor values that step (7) stablizes land use pattern grid determine: to long-term obtained in step (2)
The land use pattern grid these two types of short-term stability obtained in stable land use pattern grid and step (3) are stablized
Land use pattern grid, the vegetative coverage factor grid threshold value obtained by establishment step (5) and step (1) acquisition
Linear model between MODIS vegetation raster data obtains two classes and stablizes on land use pattern raster data in addition to arable land
Vegetative coverage factor values, and carry out grid assignment；
The vegetative coverage factor values of the short-term unstable land use pattern grid of step (8) determine: to step (4) identification
The transition form of land use pattern grid, the vegetative coverage factor grid threshold value obtained by establishment step (5) and step (1)
Linear model between the MODIS vegetation raster data of acquisition obtains on short-term unstable land use pattern raster data
Vegetative coverage factor values；
Step (9) the vegetative coverage factor merges: the arable land vegetative coverage obtained using the method that space merges to step (6)
The vegetative coverage factor values and step (8) for the stabilization land use pattern grid that factor values, step (7) obtain obtain short
The phase vegetative coverage factor values of unstable land use pattern grid merge, and obtain complete vegetative coverage factor space point
Butut.
Preferably, land use pattern grid stability preliminary judgement described in step (2), is to pass through land use
What the first-level class of type was determined；Land use pattern includes arable land, 4 forest land, meadow and Unutilized Energy fractions
Class and 16 secondary classifications；The MODIS vegetation raster data is MOD44B data, and temporal resolution is year, space point
Resolution is 250 meters, and each grid includes tree coverage rate, three layer data of non-tree coverage rate and exposed soil coverage rate.
Preferably, land use pattern grid steady in a long-term described in step (2) referred in the time to be evaluated
In span scope, there is no the grids of variation for land use pattern；Long-term non-stable land use pattern grid refers to
Within the scope of time span to be evaluated, the changed grid of land use pattern.
Preferably, the land use pattern grid of the short-term stability referred to the time to be evaluated in step (3)
Span scope is divided into several subspan ranges, and land use pattern is not sent out in any two adjacent sub- span scopes
The grid for changing；The short-term non-stable land use pattern grid refers in all two adjacent sub- span scopes
The interior changed grid of land use pattern.
Preferably, the transition form of land use pattern grid described in step (4) includes that forest land and meadow are mutual
Conversion and meadow are mutually converted with arable land.
Preferably, arable land described in the vegetative coverage factor values and step (6) in arable land described in step (5)
The vegetative coverage factor value range of land use pattern in addition includes maximum value, minimum value, these values pass through meta points
Analysis obtains.
Preferably, the vegetative coverage factor values of stabilization land use pattern grid described in step (7) pass through three kinds
Linear model is simulated, the vegetative coverage factor in the simulation of the vegetative coverage factor, stable meadow on respectively stable forest land
Simulation, and stablize the vegetative coverage factor simulation on unused land, simulation formula is specific as follows:
CWoods=MinWoods+(MaxWoods-MinWoods)×(100-TC)×0.01
Wherein, CWoodsIt is the vegetative coverage factor values stablized on forest land, MaxWoodsIt is to stablize on forest land the vegetative coverage factor most
Big value, MinWoodsIt is the minimum value for stablizing the vegetative coverage factor on forest land, TC is the coverage rate containing tree of MODIS vegetation raster data
Value；
CGrass=MinGrass+(MaxGrass-MinGrass)×NVS×0.01
Wherein, CGrassIt is the vegetative coverage factor values stablized on meadow, MaxGrassIt is to stablize on meadow the vegetative coverage factor most
Big value, MinGrassIt is the minimum value for stablizing the vegetative coverage factor on meadow, NVS is the exposed soil coverage rate of MODIS vegetation raster data
Value；
CNot=MinNot+ (Max not-MinNot)×NVS×0.01
Wherein, CNotIt is the vegetative coverage factor values stablized on unused land, MaxNotIt is to stablize vegetative coverage on unused land
The maximum value of the factor, MinNotIt is the minimum value for stablizing the vegetative coverage factor on unused land, NVS is MODIS vegetation raster data
Exposed soil coverage value.
Preferably, the vegetative coverage factor values of short-term unstable land use pattern grid described in step (8) are true
It surely is simulated by two kinds of linear models, the vegetative coverage factor that respectively short-term unstable forest land and meadow mutually convert
The vegetative coverage factor that simulation and short-term unstable meadow and arable land mutually convert is simulated, specific formula is as follows:
CWoods grass=MinWoods grass+(MaxWoods grass-MinWoods grass)×NVS×0.01
Wherein, CWoods grassIt is the vegetative coverage factor values that short-term unstable forest land and meadow mutually convert, MaxWoods grassIt is short-term non-
Stablize the vegetative coverage factor maximum value that forest land and meadow mutually convert, MaxWoods grassThe minimum value of the vegetative coverage factor on=meadow；
MinWoods grassIt is the vegetative coverage factor minimum value that short-term unstable forest land and meadow mutually convert, MinWoods grassVegetative coverage on=forest land
The maximum value of the factor；NVS is the exposed soil coverage value of MODIS vegetation raster data；
CGrass is ploughed=MinGrass is ploughed+(MaxGrass is ploughed-MinGrass is ploughed)×NVS×0.01
Wherein, CGrass is ploughedIt is the vegetative coverage factor values that short-term unstable meadow and arable land mutually convert, MaxGrass is ploughedIt is short-term non-
Stablize the vegetative coverage factor maximum value that meadow and arable land mutually convert, MaxGrass is ploughedThe minimum value of the vegetative coverage factor on=arable land；
MinGrass is ploughedIt is the vegetative coverage factor minimum value that short-term unstable meadow and arable land mutually convert, MinGrass is ploughedVegetative coverage on=meadow
The maximum value of the factor；NVS is the exposed soil coverage value of MODIS vegetation raster data.
Compared with prior art, the invention has the advantages that:
The present invention estimates the vegetative coverage factor in extensive area, can be to same place different time
The variability feature of the vegetative coverage factor carries out quantitative analysis in land use pattern, while can also be to forest land-meadow, grass
Vegetative coverage factor variations caused by during mutually converting between ground-arable land carry out dynamic analog.This method overcomes biography
It is special that system assignment method ignores vegetative coverage factor dynamic variation in the land use pattern change procedure of same place different time
The drawbacks of sign, provides new method for big regional scale vegetative coverage factor dynamic estimation, and the water and soil of big regional scale is protected
Policy making is held with reference significance.
Detailed description of the invention
Fig. 1 is 2015 Tibetan Plateau vegetative coverage factor space distribution map of present embodiment.
Specific embodiment
The invention will be further described with example with reference to the accompanying drawing.
The Method of fast estimating of the vegetative coverage factor of the invention the following steps are included:
Step (1) data acquisition: according to time span range to be evaluated, the soil benefit within the scope of the time span is obtained
With type raster data and MODIS vegetation raster data.
Step (2) land use pattern grid stability preliminary judgement: the land use pattern grid that step (1) is obtained
Data carry out multidate land use pattern grid stability preliminary judgement, obtain land use pattern grid steady in a long-term and
Long-term non-stable land use pattern grid；
Wherein, land use pattern grid stability preliminary judgement is carried out by the first-level class of land use pattern
Determine；Land use pattern includes arable land, forest land, 4 first-level class in meadow and Unutilized Energy and 16 secondary classifications.
For the MODIS vegetation raster data used in the present invention for MOD44B data, temporal resolution is year, and spatial resolution is 250 meters,
Each grid includes tree coverage rate, three layer data of non-tree coverage rate and exposed soil coverage rate.
The specific differentiation side of land use pattern grid and long-term non-stable land use pattern grid steady in a long-term
Method is as follows: land use pattern grid steady in a long-term refers within the scope of time span to be evaluated, land use pattern one
Directly there is no the grids of variation；Long-term non-stable land use pattern grid refers in time span range to be evaluated
Interior, the grid of variation occurred for land use pattern.
Step (3) land use pattern grid stability further determines: the long-term non-stable soil that step (2) is obtained
Ground use pattern grid further determines, obtain short-term stability land use pattern grid and short-term non-stable land use
Type grid, specific judgment method are as follows:
The land use pattern grid of short-term stability refer to by time span range to be evaluated be divided into several height across
Range is spent, there is no the grids of variation for land use pattern in any two adjacent sub- span scopes；And it is ungratified
It is then divided into short-term non-stable land use pattern grid, that is, refers to the land use in all two adjacent sub- span scopes
The changed grid of type.
Step (4) land use pattern grid conversion identification: short-term non-stable soil benefit obtained in analytical procedure (3)
With type grid, the transition form of land use pattern grid in such land use pattern grid is identified and tracked.Soil benefit
Transition form with type grid includes that forest land mutually converts and meadow and arable land mutually convert two classes with meadow.
The vegetative coverage factor assignment in step (5) arable land: to land use pattern steady in a long-term obtained in step (2)
In arable land grid and step (3) obtained in short-term stability land use pattern in arable land grid directly assigned
Value.The vegetative coverage factor values in arable land are analyzed to obtain by meta.
Step (6) carries out assignment to the vegetative coverage factor value range of the land use pattern other than arable land: according to step
Suddenly the land use pattern raster data that (1) obtains, determine the vegetative coverage of different land use type in addition to arable land because
The grid value range of son.The vegetative coverage factor value range of land use pattern other than arable land includes the maximum of grid
Value, minimum value, these values pass through meta and analyze to obtain.
The vegetative coverage factor values that step (7) stablizes land use pattern grid determine: to long-term obtained in step (2)
The land use pattern grid these two types of short-term stability obtained in stable land use pattern grid and step (3) are stablized
Land use pattern grid, the vegetative coverage factor grid threshold value obtained by establishment step (6) and step (1) acquisition
Linear model between MODIS vegetation raster data obtains two classes and stablizes on land use pattern raster data in addition to arable land
Vegetative coverage factor values, and carry out grid assignment.Its specific calculation method is as follows:
The vegetative coverage factor values for stablizing land use pattern grid are simulated by three kinds of linear models, respectively surely
Determine the vegetative coverage factor simulation on forest land, stablize the vegetative coverage factor simulation on meadow, and stablizes on unused land
The simulation of the vegetative coverage factor, simulation formula are specific as follows:
CWoods=MinWoods+(MaxWoods-MinWoods)×(100-TC)×0.01
Wherein, CWoodsIt is the vegetative coverage factor values stablized on forest land, MaxWoodsIt is to stablize on forest land the vegetative coverage factor most
Big value, MinWoodsIt is the minimum value for stablizing the vegetative coverage factor on forest land, TC is the coverage rate containing tree of MODIS vegetation raster data
Value；
CGrass=MinGrass+(MaxGrass-MinGrass)×NVS×0.01
Wherein, CGrassIt is the vegetative coverage factor values stablized on meadow, MaxGrassIt is to stablize on meadow the vegetative coverage factor most
Big value, MinGrassIt is the minimum value for stablizing the vegetative coverage factor on meadow, NVS is the exposed soil coverage rate of MODIS vegetation raster data
Value；
CNot=MinNot+(MaxNot-MinNot)×NVS×0.01
Wherein, CNotIt is the vegetative coverage factor values stablized on unused land, MaxNotIt is to stablize vegetative coverage on unused land
The maximum value of the factor, MinNotIt is the minimum value for stablizing the vegetative coverage factor on unused land, NVS is MODIS vegetation raster data
Exposed soil coverage value.
The vegetative coverage factor values of the short-term unstable land use pattern grid of step (8) determine: to step (4) identification
The transition form of land use pattern grid, the vegetative coverage factor grid threshold value obtained by establishment step (6) and step (1)
Linear model between the MODIS vegetation raster data of acquisition obtains on short-term unstable land use pattern raster data
Vegetative coverage factor values；
The vegetative coverage factor values determination of short-term unstable land use pattern grid is carried out by two kinds of linear models
Simulation, the vegetative coverage factor modulus that respectively short-term unstable forest land and meadow mutually convert fit short-term unstable meadow and plough
The vegetative coverage factor simulation that ground mutually converts, specific formula are respectively as follows:
CWoods grass=MinWoods grass+(MaxWoods grass-MinWoods grass)×NVS×0.01
Wherein, CWoods grassIt is the vegetative coverage factor values that short-term unstable forest land and meadow mutually convert, MaxWoods grassIt is short-term non-
Stablize the vegetative coverage factor maximum value that forest land and meadow mutually convert, MaxWoods grassThe minimum value of the vegetative coverage factor on=meadow；
MinWoods grassIt is the vegetative coverage factor minimum value that short-term unstable forest land and meadow mutually convert, MinWoods grassVegetative coverage on=forest land
The maximum value of the factor；NVS is the exposed soil coverage value of MODIS vegetation raster data；
CGrass is ploughed=MinGrass is ploughed+(MaxGrass is ploughed-MinGrass is ploughed)×NVS×0.01
Wherein, CGrass is ploughedIt is the vegetative coverage factor values that short-term unstable meadow and arable land mutually convert, MaxGrass is ploughedIt is short-term non-
Stablize the vegetative coverage factor maximum value that meadow and arable land mutually convert, MaxGrass is ploughedThe minimum value of the vegetative coverage factor on=arable land；
MinGrass is ploughedIt is the vegetative coverage factor minimum value that short-term unstable meadow and arable land mutually convert, MinGrass is ploughedVegetative coverage on=meadow
The maximum value of the factor；NVS is the exposed soil coverage value of MODIS vegetation raster data.
Step (9) the vegetative coverage factor merges: the arable land vegetative coverage obtained using the method that space merges to step (5)
The vegetative coverage factor values and step (8) for the stabilization land use pattern grid that factor values, step (7) obtain obtain short
The phase vegetative coverage factor values of unstable land use pattern grid merge, and obtain complete vegetative coverage factor space point
Butut.
Below based on this method, technical effect of the invention is described further in conjunction with specific embodiments.
Embodiment
Choose Qinghai-Tibet Platean be survey region, to the vegetative coverage in the different land use type of 2000-2015 because
Son is quickly estimated, the Qinghai-Tibet vegetative coverage factor distribution characteristics figure of 1km spatial resolution is finally obtained:
(1) it obtains in Qinghai-Tibet range between 2000-2015, the land use pattern grid of 1km spatial resolution
Data, the data acquisition is from the more period land use windy and sandy soils of China in Chinese Academy of Sciences's Resources and environmental data cloud platform
Remote sensing monitoring data set (CNLUCC) (http://www.resdc.cn/doi/doi.aspx? doiid=54), the data are with 5
Nian Weiyi period, i.e. data in a land use pattern figure represent 5 years results.Land use pattern include arable land,
Forest land, 4 first-level class in meadow and Unutilized Energy and 16 secondary classifications.Obtain 2000-2015 in Qinghai-Tibet range
Between year, temporal resolution is year, MODIS (the Moderate Resolution Imaging that spatial resolution is 250 meters
Spectroradiometer, Moderate Imaging Spectroradiomete) vegetation raster data product MOD44B data, each grid includes
Set coverage rate, three layer data of non-tree coverage rate and exposed soil coverage rate.The influence of climate variation and mankind's activity, 2000-
Qinghai-Tibet land use pattern is there is certain variation between 2015, including between forest land-meadow mutually convert with
And mutually converting between meadow-arable land.These land use pattern change procedures will be to earth's surface vegetation coverage (as tree is covered
Lid rate, non-tree coverage rate and exposed soil coverage rate) it has an impact, to change vegetative coverage of the same grid not between the same year border
Factor values.Therefore need that the stability of land use pattern and variation are determined and tracked, and simulate vegetative coverage because
The change procedure of son.
(2) land use pattern grid stability preliminary judgement: the soil the 2000-2015 benefit that step (1) will be obtained
Land use pattern determination of stability is carried out with type raster data, obtains 2000-2015 land use class steady in a long-term
Type grid and the long-term non-stable land use pattern grid of 2000-2015.Land use pattern grid steady in a long-term is
Refer between 2000-2015, vegetative coverage factor values of the land use pattern there is no the grid of variation, on these grids
Variation does not surpass the threshold range of vegetative coverage factor values in this land use pattern；Long-term non-stable land use pattern
Grid refers between 2000-2015, the changed grid of land use pattern, the vegetative coverage factor on these grids
Value changes in the different vegetative coverage factor values threshold ranges in different land use type.
(3) land use pattern grid stability further determines: 2000-2015 is divided into 2000-2005,2005-
2010 and 2010-2015 three phases, and the long-term non-stable land use pattern of 2000-2015 that step (2) is obtained
Grid further determines, respectively obtains the land use pattern grid and short-term non-stable land use pattern grid of short-term stability
Lattice.The land use pattern grid of short-term stability referred in these three stages of 2000-2005,2005-2010 and 2010-2015
In, there are the land use classes of an adjacent two stages (2000-2005 and 2005-2010,2005-2010 and 2010-2015)
There is no the grids of variation for type；Short-term non-stable land use pattern grid refer in 2000-2005,2005-2010 and
In these three stages of 2010-2015, the changed grid of the land use pattern in two neighboring stage, on these grids
Land use pattern variation acutely, and different vegetation of the vegetative coverage factor values on grid in different land use type
Significant change occurs in blanketing fctor value threshold range.
(4) land use pattern grid conversion identification: the short-term non-stable land use pattern grid of analytical procedure (3),
It identifies and tracks land use pattern grid transition form.Determine land use pattern grid transition form between 2000-2015
It is mainly mutually converted including forest land with meadow and meadow is mutually converted with arable land.
(5) the vegetative coverage factor assignment ploughed: the land use pattern raster data obtained according to step (1), to step
Suddenly short-term stability obtained in the arable land grid in land use pattern steady in a long-term obtained in (2) and step (3)
Arable land grid in land use pattern carries out indirect assignment.The vegetative coverage factor values in arable land are analyzed to obtain by meta.
(6) assignment is carried out to the vegetative coverage factor value range of the land use pattern other than arable land: according to step (1)
The land use pattern raster data of acquisition determines the grid value range of the vegetative coverage factor of different land use type.
The vegetative coverage factor value range of land use pattern other than arable land includes the maximum value of grid, minimum value, these values are equal
It analyzes to obtain by meta.
(7) the vegetative coverage factor values for stablizing land use pattern grid determine: steady in a long-term to mentioning in step (2)
Land use pattern grid and step (3) in the land use pattern grid of short-term stability mentioned, pass through establishment step
(6) linear model between MODIS vegetation raster data that the vegetative coverage factor grid threshold value mentioned and step (1) are mentioned,
Obtain stablizing the vegetative coverage factor values on land use pattern raster data.The corresponding vegetative coverage of different land use type
Factor values range is as shown in table 1.Linear model includes three kinds, and the vegetative coverage factor on respectively stable forest land is simulated, stablized
Vegetative coverage factor simulation on meadow, and stablize the vegetative coverage factor simulation on unused land, simulate formula specifically such as
Under:
CWoods=MinWoods+(MaxWoods-MinWoods)×(100-TC)×0.01
Wherein, CWoodsIt is the vegetative coverage factor values stablized on forest land, MaxWoodsIt is to stablize on forest land the vegetative coverage factor most
Big value, MinWoodsIt is the minimum value for stablizing the vegetative coverage factor on forest land, TC is the coverage rate containing tree of MODIS vegetation raster data
Value；
CGrass=MinGrass+(MaxGrass-MinGrass)×NVS×0.01
Wherein, CGrassIt is the vegetative coverage factor values stablized on meadow, MaxGrassIt is to stablize on meadow the vegetative coverage factor most
Big value, MinGrassIt is the minimum value for stablizing the vegetative coverage factor on meadow, NVS is the exposed soil coverage rate of MODIS vegetation raster data
Value；
CNot=MinNot+(MaxNot-MinNot)×NVS×0.01
Wherein, CNotIt is the vegetative coverage factor values stablized on unused land, MaxNotIt is to stablize vegetative coverage on unused land
The maximum value of the factor, MinNotIt is the minimum value for stablizing the vegetative coverage factor on unused land, NVS is MODIS vegetation raster data
Exposed soil coverage value.
(8) the short-term unstable land use pattern vegetative coverage factor determines: the land use pattern mentioned to step (4)
Transition form obtains the vegetative coverage factor threshold that step (6) is mentioned and the MOD44B number that step (1) is mentioned by linear model
The relationship between obtains the vegetative coverage factor values of each grid in short-term unstable land use pattern.Linear model includes
Two kinds, the vegetative coverage factor modulus that respectively short-term unstable forest land and meadow mutually convert fits short-term unstable meadow and ploughs
The vegetative coverage factor simulation that ground mutually converts, specific formula is as follows:
CWoods grass=MinWoods grass+(MaxWoods grass-MinWoods grass)×NVS×0.01
Wherein, CWoods grassIt is the vegetative coverage factor values that short-term unstable forest land and meadow mutually convert, MaxWoods grassIt is short-term non-
Stablize the vegetative coverage factor maximum value that forest land and meadow mutually convert, MinWoods grassIt is that short-term unstable forest land mutually turns with meadow
The vegetative coverage factor minimum value of change, NVS are the exposed soil coverage values of MODIS vegetation raster data；
CGrass is ploughed=MinGrass is ploughed+(MaxGrass is ploughed-MinGrass is ploughed)×NVS×0.01
Wherein, CGrass is ploughedIt is the vegetative coverage factor values that short-term unstable meadow and arable land mutually convert, MaxGrass is ploughedIt is short-term non-
Stablize the vegetative coverage factor maximum value that meadow and arable land mutually convert, MinGrass is ploughedIt is that short-term unstable meadow and arable land mutually turn
The vegetative coverage factor minimum value of change, NVS are the exposed soil coverage values of MODIS vegetation raster data.
(9) the vegetative coverage factor merges: the arable land vegetative coverage factor obtained using the method that space merges to step (5)
The vegetative coverage factor values and step (8) for the stabilization land use pattern grid that value, step (7) obtain obtain short-term non-
The vegetative coverage factor values for stablizing land use pattern grid merge, and are planted based on 2015 Tibetan Plateaus that this method calculates
The spatial distribution characteristic of the capped factor is as shown in Figure 1.
The corresponding vegetative coverage factor values range of the Qinghai-Tibet land use pattern of table 1
The present embodiment being capable of vegetative coverage factor space-time in the land use pattern of quantitative analysis same position different time
Variation features, the vegetative coverage factor variations in dynamic analog land use pattern conversion process, overcoming conventional method cannot
It estimates this drawback of the dynamic change of the vegetative coverage factor in process of land use change, is soil erosion mould in big regional scale
It is quasi- to provide new approaches.
Claims (8)
1. a kind of Method of fast estimating of the vegetative coverage factor, which comprises the following steps:
Step (1) data acquisition: according to time span range to be evaluated, the land use class within the scope of the time span is obtained
Type raster data and MODIS vegetation raster data；
Step (2) land use pattern grid stability preliminary judgement: the land use pattern raster data that step (1) is obtained
Multidate land use pattern grid stability preliminary judgement is carried out, land use pattern grid steady in a long-term and long-term is obtained
Non-stable land use pattern grid；
Step (3) land use pattern grid stability further determines: the long-term non-stable soil benefit that step (2) is obtained
Further determined with type grid, obtain short-term stability land use pattern grid and short-term non-stable land use pattern
Grid；
Step (4) land use pattern grid conversion identification: short-term non-stable land use class obtained in analytical procedure (3)
Type grid identifies and tracks the transition form of land use pattern grid in such land use pattern grid；
The vegetative coverage factor assignment in step (5) arable land: in land use pattern steady in a long-term obtained in step (2)
Arable land grid in the land use pattern of short-term stability obtained in grid and step (3) of ploughing carries out indirect assignment；
Step (6) carries out assignment to the vegetative coverage factor value range of the land use pattern other than arable land: according to step (1)
The land use pattern raster data of acquisition determines the grid of the vegetative coverage factor of the different land use type in addition to arable land
Value range；
The vegetative coverage factor values that step (7) stablizes land use pattern grid determine: to steady in a long-term obtained in step (2)
Land use pattern grid and step (3) obtained in short-term stability land use pattern grid these two types stablize soil
The MODIS of use pattern grid, the vegetative coverage factor grid threshold value obtained by establishment step (5) and step (1) acquisition plants
By the linear model between raster data, the vegetative coverage on the stable land use pattern raster data of two classes in addition to arable land is obtained
Factor values, and carry out grid assignment；
The vegetative coverage factor values of the short-term unstable land use pattern grid of step (8) determine: to the soil of step (4) identification
The transition form of use pattern grid, the vegetative coverage factor grid threshold value obtained by establishment step (6) and step (1) obtain
MODIS vegetation raster data between linear model, obtain the vegetation on short-term unstable land use pattern raster data
Blanketing fctor value；
Step (9) the vegetative coverage factor merges: the arable land vegetative coverage factor obtained using the method that space merges to step (5)
The vegetative coverage factor values and step (8) for the stabilization land use pattern grid that value, step (7) obtain obtain short-term non-
The vegetative coverage factor values for stablizing land use pattern grid merge, and obtain complete vegetative coverage factor space distribution
Figure.
2. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, which is characterized in that in step (2)
The land use pattern grid stability preliminary judgement, is determined by the first-level class of land use pattern；
Land use pattern includes arable land, forest land, 4 first-level class in meadow and Unutilized Energy and 16 secondary classifications；Described
MODIS vegetation raster data is MOD44B data, and temporal resolution is year, and spatial resolution is 250 meters, and each grid includes tree
Three layer data of coverage rate, non-tree coverage rate and exposed soil coverage rate.
3. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, which is characterized in that in step (2)
The land use pattern grid steady in a long-term refers to that within the scope of time span to be evaluated, land use pattern does not have
Changed grid；Non-stable land use pattern grid refers within the scope of time span to be evaluated for a long time, soil
The changed grid of use pattern.
4. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, it is characterised in that in step (3),
The land use pattern grid of the short-term stability, which refers to, is divided into several subspans for time span range to be evaluated
Range, there is no the grids of variation for land use pattern in any two adjacent sub- span scopes；Described is short-term non-
Stable land use pattern grid refers to that land use pattern is changed in all two adjacent sub- span scopes
Grid.
5. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, which is characterized in that in step (4)
The transition form of the land use pattern grid include forest land mutually converted with meadow and meadow and arable land mutually turn
Change.
6. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, which is characterized in that in step (5)
The vegetation of land use pattern other than arable land described in the vegetative coverage factor values and step (6) in the arable land is covered
Lid factor value range includes maximum value, minimum value, these values pass through meta and analyze to obtain.
7. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, which is characterized in that in step (7)
The vegetative coverage factor values of the stabilization land use pattern grid are simulated by three kinds of linear models, respectively stable
The vegetative coverage factor simulation in the simulation of the vegetative coverage factor, stable meadow on forest land, and stablize the plant on unused land
Capped factor simulation, simulation formula are specific as follows:
CWoods=MinWoods+(MaxWoods-MinWoods)×(100-TC)×0.01
Wherein, CWoodsIt is the vegetative coverage factor values stablized on forest land, MaxWoodsIt is the maximum value for stablizing the vegetative coverage factor on forest land,
MinWoodsIt is the minimum value for stablizing the vegetative coverage factor on forest land, TC is the coverage value containing tree of MODIS vegetation raster data；
CGrass=MinGrass+(MaxGrass-MinGrass)×NVS×0.01
Wherein, CGrassIt is the vegetative coverage factor values stablized on meadow, MaxGrassIt is the maximum value for stablizing the vegetative coverage factor on meadow,
MinGrassIt is the minimum value for stablizing the vegetative coverage factor on meadow, NVS is the exposed soil coverage value of MODIS vegetation raster data；
CNot=MinNot+(MaxNot-MinNot)×NVS×0.01
Wherein, CNotIt is the vegetative coverage factor values stablized on unused land, MaxNotIt is to stablize the vegetative coverage factor on unused land
Maximum value, MinNotIt is the minimum value for stablizing the vegetative coverage factor on unused land, NVS is the naked of MODIS vegetation raster data
Native coverage value.
8. a kind of Method of fast estimating of vegetative coverage factor according to claim 1, which is characterized in that in step (8)
The vegetative coverage factor values determination of the short-term unstable land use pattern grid is to carry out mould by two kinds of linear models
Quasi-, the vegetative coverage factor modulus that respectively short-term unstable forest land and meadow mutually convert fits short-term unstable meadow and arable land
The vegetative coverage factor simulation mutually converted, specific formula is as follows:
CWoods grass=MinWoods grass+(MaxWoods grass-MinWoods grass)×NVS×0.01
Wherein, CWoods grassIt is the vegetative coverage factor values that short-term unstable forest land and meadow mutually convert, MaxWoods grassIt is short-term unstable
The vegetative coverage factor maximum value that forest land and meadow mutually convert, MaxWoods grassThe minimum value of the vegetative coverage factor on=meadow；
MinWoods grassIt is the vegetative coverage factor minimum value that short-term unstable forest land and meadow mutually convert, MinWoods grassVegetative coverage on=forest land
The maximum value of the factor；NVS is the exposed soil coverage value of MODIS vegetation raster data；
CGrass is ploughed=MinGrass is ploughed+(MaxGrass is ploughed-MinGrass is ploughed)×NVS×0.01
Wherein, CGrass is ploughedIt is the vegetative coverage factor values that short-term unstable meadow and arable land mutually convert, MaxGrass is ploughedIt is short-term unstable
The vegetative coverage factor maximum value that meadow and arable land mutually convert, MaxGrass is ploughedThe minimum value of the vegetative coverage factor on=arable land；
MinGrass is ploughedIt is the vegetative coverage factor minimum value that short-term unstable meadow and arable land mutually convert, MinGrass is ploughedVegetative coverage on=meadow
The maximum value of the factor；NVS is the exposed soil coverage value of MODIS vegetation raster data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811598780.9A CN109710718B (en) | 2018-12-26 | 2018-12-26 | Method for quickly estimating plant coverage factor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811598780.9A CN109710718B (en) | 2018-12-26 | 2018-12-26 | Method for quickly estimating plant coverage factor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109710718A true CN109710718A (en) | 2019-05-03 |
| CN109710718B CN109710718B (en) | 2020-09-04 |
Family
ID=66257707
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811598780.9A Active CN109710718B (en) | 2018-12-26 | 2018-12-26 | Method for quickly estimating plant coverage factor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109710718B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112801487A (en) * | 2021-01-21 | 2021-05-14 | 中国科学院新疆生态与地理研究所 | Method for converting land use type into plant function type, terminal and storage medium |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102288954A (en) * | 2011-08-01 | 2011-12-21 | 高吉喜 | Remote sensing estimation method for vegetation coverage of grassland |
| CN102592056A (en) * | 2012-01-16 | 2012-07-18 | 浙江大学 | Method for estimating vegetation covering-managing factors of soil erosion by remote sensing |
| US20140201667A1 (en) * | 2011-03-02 | 2014-07-17 | Barbara Schoeberl | System and Method for Generating and Displaying Climate System Models |
| KR20140131704A (en) * | 2013-05-06 | 2014-11-14 | (주)웹솔루스 | Terrain management method for analyzing hydrologic model and hydraulic model, and system for the same |
| CN104268560A (en) * | 2014-09-15 | 2015-01-07 | 复凌科技（上海）有限公司 | Land utilization identification method based on remote sensing interpretation |
| CN105004725A (en) * | 2015-08-04 | 2015-10-28 | 珠江水利委员会珠江水利科学研究院 | Method for quantitatively monitoring soil erosion change amount in real time for water and soil conservation comprehensive treatment |
| CN106165571A (en) * | 2016-07-12 | 2016-11-30 | 交通运输部科学研究院 | A kind of vegetation restorability evaluation methodology based on investigation with factor correction |
| CN106548017A (en) * | 2016-10-25 | 2017-03-29 | 中国科学院地理科学与资源研究所 | A kind of ecological construction data processing method based on LU data and NDVI data |
| CN108334662A (en) * | 2018-01-03 | 2018-07-27 | 浙江大学 | A kind of soil erosion Method of fast estimating |
| CN108846352A (en) * | 2018-06-08 | 2018-11-20 | 广东电网有限责任公司 | A kind of vegetation classification and recognition methods |
-
2018
- 2018-12-26 CN CN201811598780.9A patent/CN109710718B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140201667A1 (en) * | 2011-03-02 | 2014-07-17 | Barbara Schoeberl | System and Method for Generating and Displaying Climate System Models |
| CN102288954A (en) * | 2011-08-01 | 2011-12-21 | 高吉喜 | Remote sensing estimation method for vegetation coverage of grassland |
| CN102592056A (en) * | 2012-01-16 | 2012-07-18 | 浙江大学 | Method for estimating vegetation covering-managing factors of soil erosion by remote sensing |
| KR20140131704A (en) * | 2013-05-06 | 2014-11-14 | (주)웹솔루스 | Terrain management method for analyzing hydrologic model and hydraulic model, and system for the same |
| CN104268560A (en) * | 2014-09-15 | 2015-01-07 | 复凌科技（上海）有限公司 | Land utilization identification method based on remote sensing interpretation |
| CN105004725A (en) * | 2015-08-04 | 2015-10-28 | 珠江水利委员会珠江水利科学研究院 | Method for quantitatively monitoring soil erosion change amount in real time for water and soil conservation comprehensive treatment |
| CN106165571A (en) * | 2016-07-12 | 2016-11-30 | 交通运输部科学研究院 | A kind of vegetation restorability evaluation methodology based on investigation with factor correction |
| CN106548017A (en) * | 2016-10-25 | 2017-03-29 | 中国科学院地理科学与资源研究所 | A kind of ecological construction data processing method based on LU data and NDVI data |
| CN108334662A (en) * | 2018-01-03 | 2018-07-27 | 浙江大学 | A kind of soil erosion Method of fast estimating |
| CN108846352A (en) * | 2018-06-08 | 2018-11-20 | 广东电网有限责任公司 | A kind of vegetation classification and recognition methods |
Non-Patent Citations (2)
| Title |
|---|
| ZIQIANG MA 等: "Improving TMPA 3B43 V7 Data Sets Using Land-Surface Characteristics and Ground Observations on the Qinghai-Tibet Plateau", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
| 滕洪芬: "基于多源信息的潜在土壤侵蚀估算与数字制图研究", 《中国博士学位论文全文数据库 农业科技辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112801487A (en) * | 2021-01-21 | 2021-05-14 | 中国科学院新疆生态与地理研究所 | Method for converting land use type into plant function type, terminal and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109710718B (en) | 2020-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108921885B (en) | Method for jointly inverting forest aboveground biomass by integrating three types of data sources | |
| Clawges et al. | The use of airborne lidar to assess avian species diversity, density, and occurrence in a pine/aspen forest | |
| CN103824077B (en) | Urban impervious layer rate information extraction method based on multi-source remote sensing data | |
| CN111091079B (en) | TLS-based method for measuring vegetation advantage single plant structural parameters in friable region | |
| CN114781011A (en) | High-precision calculation method and system for pixel-level global forest carbon reserves | |
| Ørka et al. | Subalpine zone delineation using LiDAR and Landsat imagery | |
| CN103398957A (en) | Hyperspectrum and laser radar-based method for extracting vertical distribution of leaf area | |
| Du et al. | Agricultural land abandonment and retirement mapping in the Northern China crop-pasture band using temporal consistency check and trajectory-based change detection approach | |
| Hu et al. | Integrating CART algorithm and multi-source remote sensing data to estimate sub-pixel impervious surface coverage: a case study from Beijing Municipality, China | |
| Mustafa | Spatiotemporal analysis of vegetation cover in Kurdistan region-Iraq using MODIS image data | |
| Ahongshangbam et al. | Drone‐based photogrammetry‐derived crown metrics for predicting tree and oil palm water use | |
| Sudhakar Reddy et al. | Development of national database on long-term deforestation in Sri Lanka | |
| Zhang et al. | Vegetation fractional coverage change in a typical oasis region in Tarim River Watershed based on remote sensing | |
| CN101788664A (en) | Land degradation response unit-oriented high-spectrum land degradation information extracting method | |
| CN112014542A (en) | Method and device for making soil moisture map of vegetation coverage area, storage medium and equipment | |
| Borak et al. | Time series vegetation aerodynamic roughness fields estimated from MODIS observations | |
| Lam et al. | Double sampling for stratification for the monitoring of sparse tree populations: the example of Populus euphratica Oliv. forests at the lower reaches of Tarim River, Southern Xinjiang, China | |
| CN109710718A (en) | A kind of Method of fast estimating of the vegetative coverage factor | |
| CN107067384B (en) | A kind of loess plateau terraced fields extracting method based on ray method | |
| Gu et al. | Mapping forest structure and uncertainty in an urban area using leaf-off lidar data | |
| Du et al. | Characterizing spatiotemporal variations of forest canopy gaps using aerial laser scanning data | |
| Michez et al. | Characterization of forests with LiDAR technology | |
| CN111913185B (en) | TLS (TLS (visual inspection) measuring method for low shrub pattern investigation of severe cold fragile region | |
| Zhao et al. | Spatio-temporal analysis of land use changes using remote sensing in Horqin sandy land, China | |
| Zhang et al. | Spatial variability and temporal stability of actual evapotranspiration on a hillslope of the Chinese Loess Plateau |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |