CN109657994A - A kind of the effectiveness appraisal procedure and system of ecology class protected natural area - Google Patents
A kind of the effectiveness appraisal procedure and system of ecology class protected natural area Download PDFInfo
- Publication number
- CN109657994A CN109657994A CN201811581119.7A CN201811581119A CN109657994A CN 109657994 A CN109657994 A CN 109657994A CN 201811581119 A CN201811581119 A CN 201811581119A CN 109657994 A CN109657994 A CN 109657994A
- Authority
- CN
- China
- Prior art keywords
- index
- ecosystem
- effectiveness
- moment
- area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000008859 change Effects 0.000 claims abstract description 34
- 238000012544 monitoring process Methods 0.000 claims abstract description 23
- 238000011156 evaluation Methods 0.000 claims abstract description 13
- 241000894007 species Species 0.000 claims description 96
- 230000012010 growth Effects 0.000 claims description 48
- 230000007306 turnover Effects 0.000 claims description 39
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 31
- 230000000694 effects Effects 0.000 claims description 22
- 238000009826 distribution Methods 0.000 claims description 17
- 238000009825 accumulation Methods 0.000 claims description 15
- 238000013467 fragmentation Methods 0.000 claims description 12
- 238000006062 fragmentation reaction Methods 0.000 claims description 12
- 238000006467 substitution reaction Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 8
- 238000005065 mining Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 4
- 239000002689 soil Substances 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000010200 validation analysis Methods 0.000 claims 1
- 238000012512 characterization method Methods 0.000 abstract 1
- 238000011835 investigation Methods 0.000 description 4
- 230000006641 stabilisation Effects 0.000 description 4
- 238000011105 stabilization Methods 0.000 description 4
- 241000124873 Grus nigricollis Species 0.000 description 3
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 3
- 235000011613 Pinus brutia Nutrition 0.000 description 3
- 241000018646 Pinus brutia Species 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000009191 jumping Effects 0.000 description 3
- 241000272471 Circus aeruginosus Species 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 241000288140 Gruiformes Species 0.000 description 2
- 241000897717 Mergus squamatus Species 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 235000011611 Pinus yunnanensis Nutrition 0.000 description 2
- 241000018652 Pinus yunnanensis Species 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- OLEWMKVPSUCNLG-UHFFFAOYSA-N 4-methylthioamphetamine Chemical compound CSC1=CC=C(CC(C)N)C=C1 OLEWMKVPSUCNLG-UHFFFAOYSA-N 0.000 description 1
- 241000272473 Aquila chrysaetos Species 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 235000005018 Pinus echinata Nutrition 0.000 description 1
- 241001236219 Pinus echinata Species 0.000 description 1
- 235000017339 Pinus palustris Nutrition 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008635 plant growth Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses the effectiveness appraisal procedures and system of a kind of ecological class protected natural area.This method comprises: obtaining the state and the state at the second moment at the first moment of protected natural area；First moment was the initial stage in monitoring and evaluation period, and the second moment was the latter stage in monitoring and evaluation period；Obtain the state at the first moment of characterization and the index of the state at the second moment；According to the variable of index, the change rate of many kinds of parameters of protected natural area effectiveness assessment is calculated；The change rate of every kind of parameter is substituted into the corresponding exponential model constructed in advance, the effectiveness index of various parameters index is obtained；Each effectiveness index is summed to obtain integrated protection into valid value.The method of the present invention and system can reduce effectiveness and calculate error, improve effectiveness and assess accuracy.
Description
Technical field
The present invention relates to conservation of nature fields, more particularly to a kind of effectiveness assessment side of ecological class protected natural area
Method and system.
Background technique
Protecting field (nature reserve area, national park, wetland park, plasm resource protection, Forest Park, water head site, wind
Scape Scenic Area etc.) effectiveness, refer to the survival condition of protected object, the protection and management of ecological environmental condition etc.
Effect.The determination of effectiveness gradually becomes an area or even nationwide bio-diversity conservation, the ecosystem, ecology
An element task in environment supervision supervises function pipe related to supervising for analyzing and finding the problem of protecting
Reason department carries out effective protection work, has important practice significance and application value.
The determination method of traditional effectiveness is the determination that by stages carries out effectiveness, and this method is many industry portions
Door is observed and predicted, is assessed and the basic skills of early warning, has the characteristics that intuitive, succinct, readily understood and application.But according to mathematics
For, the determination that by stages carries out effectiveness can regard a kind of " step " function as, have discontinuity and jumping characteristic, can band
To calculate error.Under specific step point boundary condition, this error be it is very big, for the model with polytomy variable
For, this error is to lead to the definitive result inaccuracy of effectiveness very greatly.
Summary of the invention
The object of the present invention is to provide the effectiveness appraisal procedures and system of a kind of ecological class protected natural area, reduce and protect
It protects effect and calculates error, improve the accuracy of effectiveness assessment.
To achieve the above object, the present invention provides following schemes:
A kind of effectiveness appraisal procedure of ecology class protected natural area, comprising:
Obtain the state and the state at the second moment at the first moment of the protected natural area；First moment is monitoring
The initial stage of assessment cycle, second moment are the latter stage in monitoring and evaluation period；The monitoring cycle be to protected natural area into
The period of the state of protected natural area is monitored in row effectiveness determination process；
Obtain the index of the state and the state at second moment that characterize first moment；
According to the variable of the index, the change rate of many kinds of parameters of the protected natural area effectiveness assessment is calculated；
The change rate of every kind of parameter is substituted into the corresponding exponential model constructed in advance, the guarantor of various parameters index is obtained
Protect effect index；
Each effectiveness index is summed to obtain integrated protection into valid value.
Optionally, it is described each effectiveness index is summed to obtain integrated protection into valid value after, further includes:
The integrated protection is compared at valid value and effectiveness semantic criteria, the language for obtaining effectiveness is retouched
It states.
Optionally, the variable according to the index calculates a variety of ginsengs of the protected natural area effectiveness assessment
Several change rates, specifically includes:
Utilize formulaLandscape pattern's Effective exponent is calculated, wherein x1For landscape lattice
Office's Effective exponent；EP is landscape indices,
T was the first moment；T+1 was the second moment；m1For protectiveness landscape mosaic sum；i1, j1It is protectiveness landscape mosaic
Serial number；Si1It is i-th1The area of a protectiveness landscape mosaic, TS are the gross area of conservation of nature；FI is broken for landscape cell
Change index, between 0-1, the value of FI is bigger, and protectiveness landscape mosaic fragmentation degree is bigger, otherwise smaller；
Utilize formulaTurnover ecosystem Effective exponent is calculated,
Middle x2For turnover ecosystem Effective exponent；Assessment initial stage i-th2Class turnover ecosystem is assessment latter stage jth2Class is raw
The area of state system；m2、n2The quantity at respectively the first moment and the second moment ecosystem-type；Q is Different ecosystems class
The importance assignment of type；
When the ecosystem is non-wetlands ecosystems, formula is utilizedCalculate ecosystem quality validity
Index, wherein x31Ecosystem quality Effective exponent when for the ecosystem being non-wetlands ecosystems, n3For pixel quantity,
For pixel i3With remote sensing obtain normalized differential vegetation index change rate,NDVIst
It is the first moment pixel value by the normalized differential vegetation index value of Growing season, NDVIs(t+1)For the normalization of the second moment pixel vegetation growing season
Vegetation index value；When the ecosystem is wetlands ecosystems, formula is utilized
Ecosystem quality Effective exponent is calculated, wherein x32Ecosystem quality when for the ecosystem being wetlands ecosystems is effective
Sex index,For the first wetland wet season at moment area or water quality indicator the i-th '3A variable,When being second
Carve the wetland wet season area or water quality indicator the i-th '3A variable；i′3Value be 1,2 ..., m '3, m '3For the wetland water surface,
The variable number of water quality；
Utilize formulaForest Growth Effective exponent is calculated, wherein x4Refer to for Forest Growth validity
Number；For Forest Growth Effective exponent i-th4The secondary index change rate of item, secondary index includes canopy density EyWith accumulation Ex,
Wherein accumulation is high according to the diameter of a cross-section of a tree trunk 1.3 meters above the ground and tree, consults standard volume table and determines；For secondary index weights, w1Take 0.4, w2Take 0.6；
Utilize formulaIt calculates vegetation and stablizes Effective exponent, wherein x5Have for forest cover stabilization
Sex index is imitated, H is the measurement index of species diversity, For species i5Number of individuals Zhan it is always a
The ratio of body number, i5=1,2 ..., n5；n5For total Number of Species；
Utilize formulaProtected object species quantity Effective exponent is calculated,
Middle x6For protected object species quantity Effective exponent；For protected natural area i-th6The weight of kind protection species, according to protection
The relative importance degree of species determines；n6The species number of species is protected for protected natural area；For the first moment i-th6Kind is protected
Protect the population quantity or distribution area of species；For the second moment i-th6The population quantity of kind protection species or distribution face
Product；
Utilize formulaCalculate interference variations Effective exponent, wherein x7For interference variations validity
Index；I is mankind's moving obstacle index, For functional areas coefficient, i.e. Human impact
Land type is in protected natural area jth7The interference assignment of functional areas；It is i-th7The interference power of kind Human impact land type
Weight, protected natural area Human impact land type includes mining lease, industrial land, communications and transportation land used, Tourism land, house
Land used and Farm Land, weighted value take 0.42,0.22,0.16,0.12,0.05 and 0.03 respectively；It is i-th7Kind people
Class interferes land type in functional areas j7In area.
Optionally, the change rate by every kind of parameter substitutes into the corresponding exponential model constructed in advance, obtains various
The effectiveness index of parameter index, specifically includes:
Landscape pattern's Effective exponent is substituted into landscape pattern's Effective exponent model F (x1)=15/ (1+EXP (-
1.635-0.757*x1)) in, obtain landscape pattern's effectiveness index；
The turnover ecosystem Effective exponent is substituted into turnover ecosystem Effective exponent model F (x2)=15/
(1+EXP(-1.635-1.557*x2)) in, obtain turnover ecosystem effectiveness index；
When the ecosystem is non-wetlands ecosystems, it is raw that the ecosystem quality Effective exponent is substituted into non-wetland
State mass of system Effective exponent model F (x31(the 1+EXP (- 1.21-0.69*x of)=18/31)) in, obtain non-wetland ecosystem
System quality protection effect index；When the ecosystem is wetlands ecosystems, by the ecosystem quality Effective exponent generation
Enter Wetland Area variability index model F (x32(the 1+EXP (- 1.35-0.65*x of)=10/32)) in, obtain Wetland Area protection at
Imitate index；
The Forest Growth Effective exponent is substituted into Forest Growth Effective exponent model F (x4)=10/ (1+EXP (-
0.73-0.54*x4)) in, obtain Forest Growth effectiveness index；
The vegetation is stablized into Effective exponent and substitutes into the stable Effective exponent model F (x of vegetation5)=10/ (1+EXP (-
1.83-0.81*x5)) in, it obtains vegetation and stablizes effectiveness index；
The protected object species quantity Effective exponent is substituted into protected object species quantity Effective exponent model F
(x6(the 1+EXP (- 1.23-0.69*x of)=18/6)) in, obtain protected object species quantity effectiveness index；
The interference variations Effective exponent is substituted into interference variations Effective exponent model F (x7)=14/ (1+EXP (-
1.135+0.71*x7)) in, obtain mankind's activity Annoyance Index.
Invention additionally discloses a kind of effectiveness assessment systems of ecological class protected natural area, comprising:
Whole story state acquisition module, for obtaining the state and the shape at the second moment at the first moment of the protected natural area
State；First moment is the initial stage in monitoring and evaluation period, and second moment is the latter stage in monitoring and evaluation period；The period
For the period for carrying out the state for monitoring protected natural area in effectiveness determination process to protected natural area；
Index selection module, for obtaining the finger of the state and the state at second moment that characterize first moment
Mark；
Parameters variation computing module calculates the protected natural area effectiveness and comments for the variable according to the index
The change rate for many kinds of parameters estimated；
Effectiveness index computing module, for the change rate of every kind of parameter to be substituted into the corresponding index mould constructed in advance
In type, the effectiveness index of various parameters index is obtained；
Summation module, for summing to obtain integrated protection into valid value for the effectiveness index of various parameters.
Optionally, the effectiveness assessment system further include:
Semantic contrast module is obtained for comparing the integrated protection at valid value and effectiveness standard semantic
The language description of effectiveness.
Optionally, the Parameters variation computing module includes:
Landscape pattern's validity computing unit, for utilizing formulaCalculate landscape lattice
Office's Effective exponent, wherein x1For landscape pattern's Effective exponent；EP is landscape indices,T was the first moment；T+1 was the second moment；m1
For protectiveness landscape mosaic sum；i1, j1It is the serial number of protectiveness landscape mosaic；Si1It is i-th1A protectiveness landscape
The area of mosaics, TS are the gross area of conservation of nature；FI is landscape cell fragmentation index, and between 0-1, the value of FI is got over
Greatly, protectiveness landscape mosaic fragmentation degree is bigger, otherwise smaller；
Turnover ecosystem validity computing unit, for utilizing formula
Turnover ecosystem Effective exponent is calculated, wherein x2For turnover ecosystem Effective exponent；Assessment initial stage i-th2Class
Turnover ecosystem is assessment latter stage jth2The area of the class ecosystem；m2、n2Respectively the first moment and the second moment ecosystem
The quantity for type of uniting；Q is the importance assignment of Different ecosystems type；
Ecosystem quality validity computing unit, for utilizing formula when the ecosystem is non-wetlands ecosystemsEcosystem quality Effective exponent is calculated, wherein x31Life when for the ecosystem being non-wetlands ecosystems
State mass of system Effective exponent, n3For pixel quantity,For pixel i3With remote sensing obtain normalized differential vegetation index
Change rate,NDVIstIt is planted for the first moment pixel value by the normalization of Growing season
By index value, NDVIs(t+1)For the normalized differential vegetation index value of the second moment pixel vegetation growing season；When the ecosystem is wetland
When the ecosystem, formula is utilizedEcosystem quality validity is calculated to refer to
It counts, wherein x32Ecosystem quality Effective exponent when for the ecosystem being wetlands ecosystems,For the first moment
The i-th ' of the area of wetland wet season or water quality indicator3A variable,For the area or water of the second wetland wet season at moment
The i-th ' of matter index3A variable；i′3Value be 1,2 ..., m '3, m '3For the wetland water surface, the variable number of water quality；
Forest Growth validity computing unit, for utilizing formulaForest Growth validity is calculated to refer to
It counts, wherein x4For Forest Growth Effective exponent；For Forest Growth Effective exponent i-th4The secondary index change rate of item, it is secondary
Index includes canopy density EyWith accumulation Ex, wherein accumulation is high according to the diameter of a cross-section of a tree trunk 1.3 meters above the ground and tree, consults standard volume table and determines；It is secondary
Grade index weights, w1Take 0.4, w2Take 0.6；
Vegetation stablizes validity computing unit, for utilizing formulaThe stable validity of vegetation is calculated to refer to
It counts, wherein x5Stablizing Effective exponent for forest cover, H is the measurement index of species diversity, For species i5Number of individuals account for the ratio of total individual number, i5=1,2 ..., n5；n5For total Number of Species；
Protected object species quantity validity computing unit, for utilizing formula
Protected object species quantity Effective exponent is calculated, wherein x6For protected object species quantity Effective exponent；For natural guarantor
Shield ground i-th6The weight of kind protection species, determines according to the relative importance degree of protection species；n6For protected natural area protection
The species number of kind；For the first moment i-th6The population quantity or distribution area of kind protection species；For the second moment
I-th6The population quantity or distribution area of kind protection species；
Interference variations validity computing unit, for utilizing formulaInterference variations validity is calculated to refer to
Number, wherein x7For interference variations Effective exponent；I is mankind's moving obstacle index, For functional areas coefficient, i.e. Human impact land type is in protected natural area jth7The interference assignment of functional areas；It is
i7The interference weight of kind of Human impact land type, protected natural area Human impact land type includes mining lease, industrial
Ground, communications and transportation land used, Tourism land, residential land and Farm Land, weighted value takes 0.42 respectively, 0.22,0.16,
0.12,0.05 and 0.03；It is i-th7Kind Human impact land type is in functional areas j7In area.
Optionally, the effectiveness index computing module includes:
Landscape pattern's efficacy data substitutes into unit, for landscape pattern's Effective exponent substitution landscape pattern to have
Imitate sex index model F (x1(the 1+EXP (- 1.635-0.757*x of)=15/1)) in, obtain landscape pattern's effectiveness index；
Turnover ecosystem efficacy data substitutes into unit, gives birth to for substituting into the turnover ecosystem Effective exponent
The system converting Effective exponent model F (x of state2(the 1+EXP (- 1.635-1.557*x of)=15/2)) in, obtain turnover ecosystem
Effectiveness index；
Ecosystem quality efficacy data substitutes into unit, is used for when the ecosystem is non-wetlands ecosystems, by institute
It states ecosystem quality Effective exponent and substitutes into non-wetland ecosystem quality Effective exponent model F (x31(the 1+EXP of)=18/
(-1.21-0.69*x31)) in, obtain non-wetland ecosystem quality effectiveness index；When the ecosystem is Wetland ecological system
When system, the ecosystem quality Effective exponent is substituted into Wetland Area variability index model F (x32)=10/ (1+EXP (-
1.35-0.65*x32)) in, obtain Wetland Area effectiveness index；
Forest Growth efficacy data substitutes into unit, for Forest Growth Effective exponent substitution Forest Growth to have
Imitate sex index model F (x4(the 1+EXP (- 0.73-0.54*x of)=10/4)) in, obtain Forest Growth effectiveness index；
Vegetation stablizes efficacy data and substitutes into unit, has for the vegetation to be stablized Effective exponent substitution vegetation stabilization
Imitate sex index model F (x5(the 1+EXP (- 1.83-0.81*x of)=10/5)) in, it obtains vegetation and stablizes effectiveness index；
Protected object species quantity efficacy data substitutes into unit, for referring to the protected object species quantity validity
Number substitutes into protected object species quantity Effective exponent model F (x6(the 1+EXP (- 1.23-0.69*x of)=18/6)) in, it is protected
Protect object species quantity effectiveness index；
Interference variations efficacy data substitutes into unit, for there is interference variations Effective exponent substitution interference variations
Imitate sex index model F (x7(the 1+EXP (- 1.135+0.71*x of)=14/7)) in, obtain mankind's activity Annoyance Index.
The specific embodiment provided according to the present invention, the invention discloses following technical effects: life disclosed in this invention
The effectiveness appraisal procedure and system of state class protected natural area, with the change of each element of protected object in a time cycle,
Effectiveness is portrayed, each variable is all the index variable quantity in assessment cycle, thus using the consecutive variations of variable come really
Determine effectiveness, avoid error brought by discontinuity and jumping characteristic, effectively improves the accuracy of effectiveness assessment.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is the method flow diagram of the effectiveness appraisal procedure of the ecological class protected natural area of the embodiment of the present invention 1；
Fig. 2 is the system construction drawing of the effectiveness assessment system of the ecological class protected natural area of the embodiment of the present invention 3.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide the effectiveness appraisal procedures and system of a kind of ecological class protected natural area, reduce and protect
Shield effect determines error, improves the accuracy that effectiveness determines.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Embodiment 1:
Fig. 1 is the method flow diagram of the effectiveness appraisal procedure of the ecological class protected natural area of the embodiment of the present invention 1.
Referring to Fig. 1, the effectiveness appraisal procedure of the ecology class protected natural area, comprising:
Step 101: obtaining the state and the state at the second moment at the first moment of the protected natural area；When described first
It carves the initial stage for the monitoring and evaluation period, second moment is the latter stage in monitoring and evaluation period；The monitoring cycle is to nature
Protecting field carries out the period that the state of protected natural area is monitored in effectiveness determination process.
The data source of the step 101 of the invention includes:
(1) data collection.By searching for pertinent literature and the personnel's interview of protection zone management board, collect in protected natural area
Protection zone function division figure, in the recent scientific investigation material in protection zone in protection zone the Main Types of Vegetation spatial distribution data, protect
Shield area mainly protect it is dynamic plant be distributed, the basic datas number such as admin site distribution, protection zone Mankind action situation in protection zone
According to.
(2) remotely-sensed data.Landsat remotely-sensed data is downloaded from the NASA website (NASA), land is defended
Star remotely-sensed data can interpret and provide different times vegetative coverage/land use (ecosystem-type) as basic data source
Spatial information, the spaces such as Mankind action objective information, ecosystem quality information.
(3) secondary data.Pass through relevant industrial department's internal data, administrative department's publicity data, NGO survey of organization number
According to the collection of, scientific and technical literature data, acquisition and integration section secondary data.These data contents specifically include that tourism development number
According to, protection species data, pollution situation data etc..
(4) ground investigation data.The key data of on-site inspection acquisition includes: protection zone management condition data, study plot
Data, landscape photo, tourism development photo etc..
Step 102: obtaining the index of the state and the state at second moment that characterize first moment；
Step 103: according to the variable of the index, calculating many kinds of parameters of the protected natural area effectiveness assessment
Change rate.
Step 104: the change rate of every kind of parameter being substituted into the corresponding exponential model constructed in advance, various parameters are obtained
The effectiveness index of index.
Step 105: each effectiveness index is summed to obtain integrated protection into valid value.
Step 106: the integrated protection being compared at valid value and effectiveness semantic criteria, obtains effectiveness
Language description.
Table 1 is effectiveness standard scale.
Table 1
| Integrated protection is at valid value section | [85,100] | (70,85] | (55,70] | (40,55] | [0,40] |
| Effectiveness degree | It is outstanding | Well | Generally | Difference | It is very poor |
Referring to table 1, effectiveness degree can be determined according to table 1.
As an alternative embodiment, step 103 specifically includes:
Utilize formulaLandscape pattern's Effective exponent is calculated, wherein x1It is effective for landscape pattern
Sex index；EP is landscape indices, T is
One moment；T+1 was the second moment；m1For protectiveness landscape mosaic sum；i1, j1It is the sequence of protectiveness landscape mosaic
Number；Si1It is i-th1The area of a protectiveness landscape mosaic, TS are the gross area of protected natural area；FI is landscape cell fragmentation
Index, between 0-1, the value of FI is bigger, and protectiveness landscape mosaic fragmentation degree is bigger, otherwise smaller；
Utilize formulaTurnover ecosystem validity is calculated to refer to
It counts, wherein x2For turnover ecosystem Effective exponent；Assessment initial stage i-th2Class turnover ecosystem is assessment latter stage the
j2The area of the class ecosystem；m2、n2The quantity at respectively the first moment and the second moment ecosystem-type；Q is different ecological
The importance assignment of system type；Ecosystem-type is divided into the most important ecosystem (I), the secondary important ecosystem
(II), third important ecosystem (III), slight disturbance ecology system (IV) and severe disturbance ecology system (V) 5 classes are specific
Assignment.Table 2 is ecosystem-type importance assignment reference table.
Table 2
| Ecosystem-type | Importance assignment |
| The most important ecosystem | 1.0 |
| The secondary important ecosystem | 0.6 |
| The important ecosystem of third | 0.4 |
| The ecosystem slightly interfered | -0.4 |
| The ecosystem of severe interference | -1.0 |
It, can be according to table 2 to Q assignment referring to table 2.
Ecosystem-type mainly includes forest, spinney, meadow, wetland, glacier and firn, bare area, agricultural life
Produce the classes such as land used and building site (including residential land, communications and transportation land used, Tourism land, industrial land and mining lease)
Type.According to each protected natural area feature, classify to the ecosystem being affected.Firstly, determining that its master is to be protected
The ecosystem, is usually no more than 3 classes, it is respectively most important, secondary it is important, third is important.Except part Farm Land is as guarantor
Species habitat is protected, other Farm Lands are classified as slight interference type.Building site is classified as severe interference type.
When the ecosystem is non-wetlands ecosystems, formula is utilizedCalculate ecosystem quality validity
Index, wherein x31Ecosystem quality Effective exponent when for the ecosystem being non-wetlands ecosystems, n3For pixel quantity,
For pixel i3With remote sensing obtain normalized differential vegetation index (NDVI) change rate,
NDVIstIt is the first moment pixel value by the normalized differential vegetation index value of Growing season, NDVIs(t+1)It is raw for the second moment pixel vegetation
The normalized differential vegetation index value in long season；When the ecosystem is wetlands ecosystems, ecosystem quality Effective exponent abundance of water
The wet floor space variation of phase, water quality indicator variable quantity calculate, and utilize formula
Ecosystem quality Effective exponent is calculated, wherein x32Ecosystem quality when for the ecosystem being wetlands ecosystems is effective
Sex index,For the first wetland wet season at moment area or water quality indicator the i-th '3A variable,When being second
Carve the wetland wet season area or water quality indicator the i-th '3A variable；i′3Value be 1,2 ..., m '3, m '3For the wetland water surface,
The variable number of water quality；
Utilize formulaForest Growth Effective exponent is calculated, wherein x4Refer to for Forest Growth validity
Number；For Forest Growth Effective exponent i-th4The secondary index change rate of item, secondary index includes canopy density EyWith accumulation Ex,
Wherein accumulation is high according to the diameter of a cross-section of a tree trunk 1.3 meters above the ground and tree, consults standard volume table and determines；For secondary index weights, w1Take 0.4, w2Take 0.6；
C1It is accumulation variable quantity for canopy density variable quantity.Wherein the accumulation of bottom class or block can be according to bottom class, block
Mean DBH increment and tree are high, consult standard volume table and obtain.If investigation object area very little, it is also possible to every wood dipping data meter
It calculates.
Utilize formulaIt calculates vegetation and stablizes Effective exponent, wherein x5Have for forest cover stabilization
Sex index is imitated, H is the measurement index of species diversity, using Shannon-prestige Nasdaq index number (Shannon-Wiener index), according to Qiao
The species composition situation of timber layer, shrub layer, herbaceous layer, by 0.5,0.3,0.2 weighted calculation, For species i5Number of individuals account for the ratio of total individual number, i5=1,2 ..., n5；n5For total Number of Species；
Utilize formulaProtected object species quantity Effective exponent is calculated,
Middle x6For protected object species quantity Effective exponent；For protected natural area i-th6The weight of kind protection species, according to protection
The relative importance degree of species determines；n6The species number that species are protected for protected natural area, when there are many main for protected natural area
When protecting species, 5 kinds are at most chosen, i.e. n6≤5；For the first moment i-th6The population quantity or distribution of kind protection species
Area；Animal population quantity absolute quantity or relative populations (such as meeting rate)；Plant population indicates that quantity is dilute with distribution area
Few Precious, Rare, Endangered, endemic species population quantity are indicated with strain (clump).For the second moment i-th6The population of kind protection species
Quantity or distribution area；
Utilize formulaCalculate interference variations Effective exponent, wherein x7For interference variations validity
Index；I is mankind's moving obstacle index, For functional areas coefficient, i.e. Human impact
Land type is in protected natural area jth7The interference assignment of functional areas；Work as j7When=1 (core space),Work as j7=2
When (buffer area),Work as j7When=3 (test block),It is i-th7Kind Human impact land type
Weight is interfered, protected natural area Human impact land type includes mining lease, industrial land, communications and transportation land used, travelling
Ground, residential land and Farm Land, weighted value take 0.42,0.22,0.16,0.12,0.05 and 0.03 respectively；It is
i7Kind Human impact land type is in functional areas j7In area.Influence and diffusion in view of communications and transportation land used to biology are imitated
It answers, communications and transportation land area is to influence width measuring and calculating multiplied by it according to traffic route length to get.Wherein, railway and high speed are public
The influence width on road takes 3 times of its width of subgrade, and Class I highway takes 2 times of its width of subgrade, other roads take its width of subgrade
1 times.The area of road does buffer area (buffer distance influences width for it) to both sides of the road by ArcGIS software and acquires.
As an alternative embodiment, step 104 specifically includes:
Landscape pattern's Effective exponent is substituted into landscape pattern's Effective exponent model F (x1)=15/ (1+EXP (-
1.635-0.757*x1)) in, obtain landscape pattern's effectiveness index；
The turnover ecosystem Effective exponent is substituted into turnover ecosystem Effective exponent model F (x2)=15/
(1+EXP(-1.635-1.557*x2)) in, obtain turnover ecosystem effectiveness index；
When the ecosystem is non-wetlands ecosystems, it is raw that the ecosystem quality Effective exponent is substituted into non-wetland
State mass of system Effective exponent model F (x31(the 1+EXP (- 1.21-0.69*x of)=18/31)) in, obtain non-wetland ecosystem
System quality protection effect index；When the ecosystem is wetlands ecosystems, by the ecosystem quality Effective exponent generation
Enter Wetland Area variability index model F (x32(the 1+EXP (- 1.35-0.65*x of)=10/32)) in, obtain Wetland Area protection at
Imitate index；
The Forest Growth Effective exponent is substituted into Forest Growth Effective exponent model F (x4)=10/ (1+EXP (-
0.73-0.54*x4)) in, obtain Forest Growth effectiveness index；
The vegetation is stablized into Effective exponent and substitutes into the stable Effective exponent model F (x of vegetation5)=10/ (1+EXP (-
1.83-0.81*x5)) in, it obtains vegetation and stablizes effectiveness index；
The protected object species quantity Effective exponent is substituted into protected object species quantity Effective exponent model F
(x6(the 1+EXP (- 1.23-0.69*x of)=18/6)) in, obtain protected object species quantity effectiveness index；
The interference variations Effective exponent is substituted into interference variations Effective exponent model F (x7)=14/ (1+EXP (-
1.135+0.71*x7)) in, obtain mankind's activity Annoyance Index.
Wherein above-mentioned each model is in statistical product and service solution (SPSS, Statistical Product
And Service Solutions) in building.
Embodiment 2:
The embodiment 2 carries out specifically for determining region as effectiveness using Lijiang, yunnan Lake Lashihai protected natural area
It is bright.With the initial stage for monitoring cycle in 2010, the latter stage for monitoring cycle in 2016 collected, Landsat in 2016 in 2010
Other data of remote sensing images, protection zone and protected object are calculated according to above-mentioned steps.
The first step, the calculating of independent variable
Landscape pattern's Effective exponent is calculated according to the calculation formula of each Validity Index, turnover ecosystem validity refers to
Number, ecosystem quality Effective exponent, Forest Growth Effective exponent, vegetation stablize Effective exponent, protected object species
Quantity Effective exponent, interference variations Effective exponent.
(1) value of landscape pattern's Effective exponent is calculated
Satellite image map is interpreted, the area of all kinds of landscape mosaics is calculated and obtain, table 3 is that Lake Lashihai Plateau Wetland is provincial
Protected natural area is all kinds of to inlay bulk area.
Table 3
| Mosaics title | (hm in 20102) | (hm in 20162) |
| Water body | 1164.02 | 1209.21 |
| Lakeside zone | 238.69 | 70.22 |
| Pinus yunnanensis | 2689.31 | 2684.46 |
| Huashan pine | 213.16 | 223.74 |
| Agricultural land | 1831.83 | 1946.13 |
| Tourism land | 38.83 | 39.98 |
| Residential land | 223.39 | 226.10 |
| Golf land used | 34.03 | 34.03 |
According to table 3, the gross area TS=6433.26 of protected natural area, landscape mosaic area summation in 2010It is calculated by Fragstats3.3
The ecosystem mosaics fragmentation index FI in initial stage and latter stage is respectively 0.582 and 0.541.EP is calculated according to each variablet、
EPt+1。
Landscape mosaic area in 2016
Then
Calculate x1,
(2) value of turnover ecosystem Effective exponent is calculated
It is wetlands ecosystems that Lake Lashihai protected natural area, which carrys out the most important ecosystem,；The secondary important ecosystem is cloud
The southern pine ecosystem and the Huashan pine ecosystem；Important third is the lakeside zone ecosystem.The soil of interpretation and two phases of production benefit
With/vegetation distribution map, and ecosystem-type transformed matrix is calculated, table 4 is the provincial protected natural area of Lake Lashihai Plateau Wetland
Each ecosystem-type transformed matrix table
Table 4
Turnover ecosystem Effective exponent is calculated in conjunction with table 4.
(3) value of ecosystem quality Effective exponent
The provincial protected natural area of Lake Lashihai Plateau Wetland belongs to wetland protected natural area, so its ecosystem quality is effective
Property mainly by Wetland Area variation and wet land water quality variation determine.Table 5 is Lake Lashihai Plateau Wetland area table.
Table 5
| Computational item | (hm in 20102) | (hm in 20162) |
| Wetland water body area | 1164.02 | 1209.21 |
This example is wetland class protecting field, SwUsing coverage of water index, then:
(4) value of Forest Growth Effective exponent
Forest Growth validity mainly stablizes validity by vegetation growth validity and vegetation to determine.Acquisition is due to lacking
The comparison sample of two phases, so plant growth validity can not be calculated using comparison sample.Table 6 is the variation of accumulation, canopy density
Tables of data.
Table 6
By area specific gravity method of weighting, the canopy density changing value of 3 protection zone totality of protection zone is calculated
C1=((2205.52 × 0.55+478.94 × 0.95+223.74 × 0.70)-(2193.34+491.05+
213.16) × 0.65)/(2193.34+491.05+213.16) × 0.65 × 100=-3.12
Accumulation changing value
C2=((2205.52 × 75.56+478.94 × 177.09+223.74 × 220.37)-(2193.34+
491.05+213.16) × 100.65 × 100=3.10
Vegetation growth Effective exponent x4
(5) vegetation stablizes the value of Effective exponent
Main vegetation is distributed as pinus yunnanensis and Huashan pine, compares two phase satellite images, calculates to obtain x according to distribution area5
(6) value of protected object species quantity Effective exponent
I grades of protection boneblack cranes of country, Mergus squamatus, black-necked crane, Circus aeruginosus, woods carving are selected to protect Species validity
Computing object.Choose I grades of protection animal black-necked cranes of country, woods carving, black crane, Mergus squamatus, Circus aeruginosus, golden eagle, Dong Fangbaihe
It is calculated as five kinds of main protection species.Since above seven kinds be all national I grades of protection animal, therefore assign equal weight, respectively for
1/7.According to investigations, black-necked crane, woods carving population disappears from 25,12 is 0, behind 5 kinds of population quantities and monitoring preliminary phase
Together, belong to stable state.
Main protected object population quantity Effective exponent
(7) value of interference variations Effective exponent
Table 7 is each mankind's moving obstacle area statistics table of Lake Lashihai wetland national level protected natural area
Table 7
According to table 7, monitoring initial stage, mankind's activity Annoyance Index in 2010 are calculated
It=0+0+0.16 × 29.05+0.12 × 38.83+0.05 × 223.38+0.03 × 1831.83=75.46
Assess latter stage, i.e., 2016 mankind's activity Annoyance Indexes
It+1=0+0+0.16 × 63.21+0.12 × 39.98+0.05 × 226.09+0.03 × 1946.13=84.59
Thus interference variations Effective exponent is calculated
Second step, the calculating of effectiveness index value
By x1=3.35 substitute into landscape pattern's Effective exponent model F (x1(the 1+EXP (- 1.635-0.757*x of)=15/1))
In obtain:
F(x1)=15/ (1+exp (- 1.635-0.757 × 3.35))=14.8
By x2=0.667 substitutes into turnover ecosystem Effective exponent model F (x2(the 1+EXP (- 1.635- of)=15/
1.557*x2)) in, it obtains:
F(x2)=15/ (1+exp (- 1.635-1.557 × 0.667))=14.0
By x32=3.9 substitute into Wetland Area variability index model F (x32(the 1+EXP (- 1.35-0.65*x of)=10/32)) in,
It obtains
F(x32)=10/ (1+exp (- 1.35-0.65 × 3.9))=17.6
By x4=0.62 substitutes into Forest Growth Effective exponent model F (x4(the 1+EXP (- 0.73-0.54*x of)=10/4))
In, it obtains
F(x4)=10/ (1+exp (- 0.73-0.54 × 0.62))=7.4
By x5=-14.28, which substitute into vegetation, stablizes Effective exponent model F (x5(the 1+EXP (- 1.83-0.81*x of)=10/5))
In, it obtains
F(x5)=10/ (1+exp (- 1.83-0.81 × (- 14.28)))=0.0006
By x6=0.0 substitutes into protected object species quantity Effective exponent model F (x6(the 1+EXP (- 1.23- of)=18/
0.69*x6)) in, it obtains
F(x6)=18/ (1+exp (- 1.23-0.69 × 0.0))=13.9
By x7=12.1 substitute into interference variations Effective exponent model F (x7(the 1+EXP (- 1.135+0.71*x of)=14/7))
In, it obtains
F(x7)=14/ (1+exp (- 1.135+0.71 × 12.1))=0.008
Third step adds up
By F (x1)、F(x2)、……、F(x7) add up
F (x)=F (x1)+F(x2)+F(x3)+F(x4)+F(x5)+F(x6)+F(x7)
=14.8+14.0+17.6+7.4+0.0006+13.9+0.008=67.7
4th step, effectiveness degree determine
The table of comparisons 1, F (x)=67.7 fall into " general " section.The effectiveness of the protection zone is general.
Embodiment 3:
Fig. 2 is the system construction drawing of 3 protected natural area effectiveness assessment system of the embodiment of the present invention.
Referring to fig. 2, the protected natural area effectiveness assessment system, comprising:
Whole story state acquisition module 201, for obtaining state and second moment at the first moment of the protected natural area
State；First moment is the initial stage in monitoring and evaluation period, and second moment is the latter stage in monitoring and evaluation period；It is described
Period is the period for carrying out monitoring the state of protected natural area to protected natural area in effectiveness determination process.
Index selection module 202, for obtaining the state and the state at second moment that characterize first moment
Index；
Parameters variation computing module 203 calculates the protected natural area effectiveness for the variable according to the index
The change rate of many kinds of parameters of assessment.
Effectiveness index computing module 204, for the change rate of every kind of parameter to be substituted into the corresponding finger constructed in advance
In exponential model, the effectiveness index of various parameters index is obtained.
Summation module 205, for summing to obtain integrated protection into valid value for the effectiveness index of various parameters.
Semantic contrast module 206 is obtained for comparing the integrated protection at valid value and effectiveness standard semantic
To the language description of effectiveness.
As an alternative embodiment, the Parameters variation computing module 203 includes:
Landscape pattern's validity computing unit, for utilizing formulaCalculate landscape lattice
Office's Effective exponent, wherein x1For landscape pattern's Effective exponent；EP is landscape indices,T was the first moment；T+1 was the second moment；m1
For protectiveness landscape mosaic sum；i1, j1It is the serial number of protectiveness landscape mosaic；Si1It is i-th1A protectiveness landscape
The area of mosaics, TS are the gross area of conservation of nature；FI is landscape cell fragmentation index, and between 0-1, the value of FI is got over
Greatly, protectiveness landscape mosaic fragmentation degree is bigger, otherwise smaller；
Turnover ecosystem validity computing unit, for utilizing formula
Turnover ecosystem Effective exponent is calculated, wherein x2For turnover ecosystem Effective exponent；Assessment initial stage i-th2Class
Turnover ecosystem is assessment latter stage jth2The area of the class ecosystem；m2、n2Respectively the first moment and the second moment ecosystem
The quantity for type of uniting；Q is the importance assignment of Different ecosystems type；
Ecosystem quality validity computing unit, for utilizing formula when the ecosystem is non-wetlands ecosystemsEcosystem quality Effective exponent is calculated, wherein x31Life when for the ecosystem being non-wetlands ecosystems
State mass of system Effective exponent, n3For pixel quantity,For pixel i3With remote sensing obtain normalized differential vegetation index
Change rate,NDVIstIt is planted for the first moment pixel value by the normalization of Growing season
By index value, NDVIs(t+1)For the normalized differential vegetation index value of the second moment pixel vegetation growing season；When the ecosystem is wetland
When the ecosystem, formula is utilizedEcosystem quality validity is calculated to refer to
It counts, wherein x32Ecosystem quality Effective exponent when for the ecosystem being wetlands ecosystems,For the first moment
The i-th ' of the area of wetland wet season or water quality indicator3A variable,For the area or water of the second wetland wet season at moment
The i-th ' of matter index3A variable；i′3Value be 1,2 ..., m '3, m '3For the wetland water surface, the variable number of water quality；
Forest Growth validity computing unit, for utilizing formulaForest Growth validity is calculated to refer to
It counts, wherein x4For Forest Growth Effective exponent；For Forest Growth Effective exponent i-th4The secondary index change rate of item, it is secondary
Grade index includes canopy density EyWith accumulation Ex, wherein accumulation is high according to the diameter of a cross-section of a tree trunk 1.3 meters above the ground and tree, consults standard volume table and determines；For
Secondary index weights, w1Take 0.4, w2Take 0.6；
Vegetation stablizes validity computing unit, for utilizing formulaThe stable validity of vegetation is calculated to refer to
It counts, wherein x5Stablizing Effective exponent for forest cover, H is the measurement index of species diversity, For species i5Number of individuals account for the ratio of total individual number, i5=1,2 ..., n5；n5For total Number of Species；
Protected object species quantity validity computing unit, for utilizing formula
Protected object species quantity Effective exponent is calculated, wherein x6For protected object species quantity Effective exponent；For natural guarantor
Shield ground i-th6The weight of kind protection species, determines according to the relative importance degree of protection species；n6For protected natural area protection
The species number of kind；For the first moment i-th6The population quantity or distribution area of kind protection species；For the second moment
i6The population quantity or distribution area of kind protection species；
Interference variations validity computing unit, for utilizing formulaInterference variations validity is calculated to refer to
Number, wherein x7For interference variations Effective exponent；I is mankind's moving obstacle index, For functional areas coefficient, i.e. Human impact land type is in protected natural area jth7The interference assignment of functional areas；It is
i7The interference weight of kind of Human impact land type, protected natural area Human impact land type includes mining lease, industrial
Ground, communications and transportation land used, Tourism land, residential land and Farm Land, weighted value takes 0.42 respectively, 0.22,0.16,
0.12,0.05 and 0.03；It is i-th7Kind Human impact land type is in functional areas j7In area.
As an alternative embodiment, the effectiveness index computing module 204 includes:
Landscape pattern's efficacy data substitutes into unit, for landscape pattern's Effective exponent substitution landscape pattern to have
Imitate sex index model F (x1(the 1+EXP (- 1.635-0.757*x of)=15/1)) in, obtain landscape pattern's effectiveness index；
Turnover ecosystem efficacy data substitutes into unit, gives birth to for substituting into the turnover ecosystem Effective exponent
The system converting Effective exponent model F (x of state2(the 1+EXP (- 1.635-1.557*x of)=15/2)) in, obtain turnover ecosystem
Effectiveness index；
Ecosystem quality efficacy data substitutes into unit, is used for when the ecosystem is non-wetlands ecosystems, by institute
It states ecosystem quality Effective exponent and substitutes into non-wetland ecosystem quality Effective exponent model F (x31(the 1+EXP of)=18/
(-1.21-0.69*x31)) in, obtain non-wetland ecosystem quality effectiveness index；When the ecosystem is Wetland ecological system
When system, the ecosystem quality Effective exponent is substituted into Wetland Area variability index model F (x32)=10/ (1+EXP (-
1.35-0.65*x32)) in, obtain Wetland Area effectiveness index；
Forest Growth efficacy data substitutes into unit, for Forest Growth Effective exponent substitution Forest Growth to have
Imitate sex index model F (x4(the 1+EXP (- 0.73-0.54*x of)=10/4)) in, obtain Forest Growth effectiveness index；
Vegetation stablizes efficacy data and substitutes into unit, has for the vegetation to be stablized Effective exponent substitution vegetation stabilization
Imitate sex index model F (x5(the 1+EXP (- 1.83-0.81*x of)=10/5)) in, it obtains vegetation and stablizes effectiveness index；
Protected object species quantity efficacy data substitutes into unit, for referring to the protected object species quantity validity
Number substitutes into protected object species quantity Effective exponent model F (x6(the 1+EXP (- 1.23-0.69*x of)=18/6)) in, it is protected
Protect object species quantity effectiveness index；
Interference variations efficacy data substitutes into unit, for there is interference variations Effective exponent substitution interference variations
Imitate sex index model F (x7(the 1+EXP (- 1.135+0.71*x of)=14/7)) in, obtain mankind's activity Annoyance Index.
The specific embodiment provided according to the present invention, the invention discloses following technical effects: it is disclosed in this invention from
Right protecting field effectiveness appraisal procedure and system are protected with the change of each element of protected object in a time cycle to portray
Protect effect, each variable is all the index variable quantity in assessment cycle, thus assessed using the consecutive variations of variable protection at
Effect, avoids error brought by discontinuity and jumping characteristic, effectively improves the accuracy of effectiveness assessment.With parameter
Change rate is assessed, and influence of the natural conditions to Evaluated effect of the protecting field of natural endowment difference can be eliminated, with a kind of unification
Scale laterally compares the effectiveness of national natural ecology class protecting field.
For the system disclosed in the embodiment, since it is corresponded to the methods disclosed in the examples, so the ratio of description
Relatively simple, reference may be made to the description of the method.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said
It is bright to be merely used to help understand method and its core concept of the invention；At the same time, for those skilled in the art, foundation
Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.
Claims (8)
1. a kind of effectiveness appraisal procedure of ecology class protected natural area characterized by comprising
Obtain the state and the state at the second moment at the first moment of the protected natural area；First moment is monitoring and evaluation
The initial stage in period, second moment are the latter stage in monitoring and evaluation period；The monitoring cycle is to protect to protected natural area
It protects in effect determination process and monitors the period of the state of protected natural area；
Obtain the index of the state and the state at second moment that characterize first moment；
According to the variable of the index, the change rate of many kinds of parameters of the protected natural area effectiveness assessment is calculated；
The change rate of every kind of parameter is substituted into the corresponding exponential model constructed in advance, obtain the protection of various parameters index at
Imitate index；
Each effectiveness index is summed to obtain integrated protection into valid value.
2. a kind of effectiveness appraisal procedure of ecological class protected natural area according to claim 1, which is characterized in that
It is described each effectiveness index is summed to obtain integrated protection into valid value after, further includes:
The integrated protection is compared at valid value and effectiveness semantic criteria, obtains the language description of effectiveness.
3. a kind of effectiveness appraisal procedure of ecological class protected natural area according to claim 1, which is characterized in that institute
The variable according to the index is stated, the change rate of many kinds of parameters of the protected natural area effectiveness assessment is calculated, it is specific to wrap
It includes:
Utilize formulaLandscape pattern's Effective exponent is calculated, wherein x1It is effective for landscape pattern
Sex index；EP is landscape indices, T is
One moment；T+1 was the second moment；m1For protectiveness landscape mosaic sum；i1, j1It is the sequence of protectiveness landscape mosaic
Number；Si1It is i-th1The area of a protectiveness landscape mosaic, Sj1For jth1The area of a protectiveness landscape mosaic；TS is nature
The gross area of protecting field；FI is landscape cell fragmentation index, and between 0-1, the value of FI is bigger, protectiveness landscape mosaic
Fragmentation degree is bigger, otherwise smaller；
Utilize formulaTurnover ecosystem Effective exponent is calculated, wherein x2
For turnover ecosystem Effective exponent；Assessment initial stage i-th2Class turnover ecosystem is assessment latter stage jth2Class ecosystem
The area of system；m2、n2The quantity at respectively the first moment and the second moment ecosystem-type；Q is Different ecosystems type
Importance assignment；
When the ecosystem is non-wetlands ecosystems, formula is utilizedEcosystem quality validity is calculated to refer to
It counts, wherein x31Ecosystem quality Effective exponent when for the ecosystem being non-wetlands ecosystems, n3For pixel quantity,
For pixel i3With remote sensing obtain normalized differential vegetation index change rate,
NDVIstIt is the first moment pixel value by the normalized differential vegetation index value of Growing season, NDVIs(t+1)It is raw for the second moment pixel vegetation
The normalized differential vegetation index value in long season；When the ecosystem is wetlands ecosystems, formula is utilizedEcosystem quality Effective exponent is calculated, wherein x32For the ecosystem
Ecosystem quality Effective exponent when for wetlands ecosystems,For the area or water of the first wetland wet season at moment
The i-th of matter index3' a variable,For the second wetland wet season at moment area or water quality indicator i-th3' a variable；
i3' value be 1,2 ..., m3', m3' for the wetland water surface, water quality variable number；
Utilize formulaForest Growth Effective exponent is calculated, wherein x4For Forest Growth Effective exponent；
For Forest Growth Effective exponent i-th4The secondary index change rate of item, secondary index includes canopy density EyWith accumulation Ex, wherein storing
Accumulated amount is high according to the diameter of a cross-section of a tree trunk 1.3 meters above the ground and tree, consults standard volume table and determines；For secondary index weights, w1Take 0.4, w2Take 0.6；
Utilize formulaIt calculates vegetation and stablizes Effective exponent, wherein x5Stablize validity for forest cover
Index, H are the measurement index of species diversity, For species i5Number of individuals Zhan it is always individual
Several ratios, i5=1,2 ..., n5；n5For total Number of Species；
Utilize formulaProtected object species quantity Effective exponent is calculated, wherein x6
For protected object species quantity Effective exponent；For protected natural area i-th6The weight of kind protection species, according to protection species
Relative importance degree determine；n6The species number of species is protected for protected natural area；For the first moment i-th6Kind protection
The population quantity or distribution area of kind；For the second moment i-th6The population quantity or distribution area of kind protection species；
Utilize formulaCalculate interference variations Effective exponent, wherein x7For interference variations Effective exponent；
I is mankind's moving obstacle index, For functional areas coefficient, i.e. Human impact soil
Type is in protected natural area jth7The interference assignment of functional areas；It is i-th7The interference weight of kind Human impact land type,
Protected natural area Human impact land type includes mining lease, industrial land, communications and transportation land used, Tourism land, house use
Ground and Farm Land, weighted value take 0.42,0.22,0.16,0.12,0.05 and 0.03 respectively；It is i-th7The kind mankind are dry
Land type is disturbed in functional areas j7In area.
4. a kind of effectiveness appraisal procedure of ecological class protected natural area according to claim 3, which is characterized in that institute
It states and substitutes into the change rate of every kind of parameter in the corresponding exponential model constructed in advance, obtain the effectiveness of various parameters index
Index specifically includes:
Landscape pattern's Effective exponent is substituted into landscape pattern's Effective exponent model F (x1(the 1+EXP (- 1.635- of)=15/
0.757*x1)) in, obtain landscape pattern's effectiveness index；
The turnover ecosystem Effective exponent is substituted into turnover ecosystem Effective exponent model F (x2(the 1+EXP of)=15/
(-1.635-1.557*x2)) in, obtain turnover ecosystem effectiveness index；
When the ecosystem is non-wetlands ecosystems, the ecosystem quality Effective exponent is substituted into non-wetland ecosystem
Unite quality validation exponential model F (x31(the 1+EXP (- 1.21-0.69*x of)=18/31)) in, obtain non-wetlands ecosystems matter
Measure effectiveness index；When the ecosystem is wetlands ecosystems, the ecosystem quality Effective exponent is substituted into wet
Ground area variability index model F (x32(the 1+EXP (- 1.35-0.65*x of)=10/32)) in, it obtains Wetland Area effectiveness and refers to
Number；
The Forest Growth Effective exponent is substituted into Forest Growth Effective exponent model F (x4(the 1+EXP (- 0.73- of)=10/
0.54*x4)) in, obtain Forest Growth effectiveness index；
The vegetation is stablized into Effective exponent and substitutes into the stable Effective exponent model F (x of vegetation5(the 1+EXP (- 1.83- of)=10/
0.81*x5)) in, it obtains vegetation and stablizes effectiveness index；
The protected object species quantity Effective exponent is substituted into protected object species quantity Effective exponent model F (x6)=
18/(1+EXP(-1.23-0.69*x6)) in, obtain protected object species quantity effectiveness index；
The interference variations Effective exponent is substituted into interference variations Effective exponent model F (x7(the 1+EXP (- 1.135+ of)=14/
0.71*x7)) in, obtain mankind's activity Annoyance Index.
5. a kind of effectiveness assessment system of ecology class protected natural area characterized by comprising
Whole story state acquisition module, for obtaining the state and the state at the second moment at the first moment of the protected natural area；
First moment is the initial stage in monitoring and evaluation period, and second moment is the latter stage in monitoring and evaluation period；The period is
The period for the state that protected natural area is monitored in effectiveness determination process is carried out to protected natural area；
Index selection module, for obtaining the index of the state and the state at second moment that characterize first moment；
Parameters variation computing module calculates the protected natural area effectiveness assessment for the variable according to the index
The change rate of many kinds of parameters；
Effectiveness index computing module, for the change rate of every kind of parameter to be substituted into the corresponding exponential model constructed in advance
In, obtain the effectiveness index of various parameters index；
Summation module, for summing to obtain integrated protection into valid value for the effectiveness index of various parameters.
6. a kind of effectiveness assessment system of ecological class protected natural area according to claim 5, which is characterized in that also
Include:
Semantic contrast module is protected for comparing the integrated protection at valid value and effectiveness standard semantic
The language description of effect.
7. a kind of effectiveness assessment system of ecological class protected natural area according to claim 5, which is characterized in that institute
Stating Parameters variation computing module includes:
Landscape pattern's validity computing unit, for utilizing formulaIt is effective to calculate landscape pattern
Sex index, wherein x1For landscape pattern's Effective exponent；EP is landscape indices,T was the first moment；T+1 was the second moment；m1
For protectiveness landscape mosaic sum；i1, j1It is the serial number of protectiveness landscape mosaic；Si1It is i-th1A protectiveness landscape
The area of mosaics, TS are the gross area of conservation of nature；FI is landscape cell fragmentation index, and between 0-1, the value of FI is got over
Greatly, protectiveness landscape mosaic fragmentation degree is bigger, otherwise smaller；
Turnover ecosystem validity computing unit, for utilizing formulaMeter
Turnover ecosystem Effective exponent is calculated, wherein x2For turnover ecosystem Effective exponent；Assessment initial stage i-th2Class is raw
State is system converting for assessment latter stage jth2The area of the class ecosystem；m2、n2Respectively the first moment and the second moment ecosystem
The quantity of type；Q is the importance assignment of Different ecosystems type；
Ecosystem quality validity computing unit, for utilizing formula when the ecosystem is non-wetlands ecosystemsEcosystem quality Effective exponent is calculated, wherein x31Life when for the ecosystem being non-wetlands ecosystems
State mass of system Effective exponent, n3For pixel quantity,For pixel i3With remote sensing obtain normalized differential vegetation index
Change rate,NDVIstIt is planted for the first moment pixel value by the normalization of Growing season
By index value, NDVIs(t+1)For the normalized differential vegetation index value of the second moment pixel vegetation growing season；When the ecosystem is wetland
When the ecosystem, formula is utilizedEcosystem quality validity is calculated to refer to
It counts, wherein x32Ecosystem quality Effective exponent when for the ecosystem being wetlands ecosystems,For the first moment
The i-th of the area of wetland wet season or water quality indicator3' a variable,For the area or water of the second wetland wet season at moment
The i-th of matter index3' a variable；i3' value be 1,2 ..., m3', m3' for the wetland water surface, water quality variable number；
Forest Growth validity computing unit, for utilizing formulaForest Growth Effective exponent is calculated,
Middle x4For Forest Growth Effective exponent；For Forest Growth Effective exponent i-th4The secondary index change rate of item, secondary index
Including canopy density EyWith accumulation Ex, wherein accumulation is high according to the diameter of a cross-section of a tree trunk 1.3 meters above the ground and tree, consults standard volume table and determines；Refer to for secondary
Mark weight, w1Take 0.4, w2Take 0.6；
Vegetation stablizes validity computing unit, for utilizing formulaThe stable validity of vegetation is calculated to refer to
It counts, wherein x5Stablizing Effective exponent for forest cover, H is the measurement index of species diversity, For species i5Number of individuals account for the ratio of total individual number, i5=1,2 ..., n5；n5For total Number of Species；
Protected object species quantity validity computing unit, for utilizing formulaMeter
Protected object species quantity Effective exponent is calculated, wherein x6For protected object species quantity Effective exponent；For conservation of nature
Ground i-th6The weight of kind protection species, determines according to the relative importance degree of protection species；n6Species are protected for protected natural area
Species number；For the first moment i-th6The population quantity or distribution area of kind protection species；For the second moment i-th6
The population quantity or distribution area of kind protection species；
Interference variations validity computing unit, for utilizing formulaInterference variations Effective exponent is calculated,
Wherein, x7For interference variations Effective exponent；I is mankind's moving obstacle index,
For functional areas coefficient, i.e. Human impact land type is in protected natural area jth7The interference assignment of functional areas；It is i-th7Kind
The interference weight of Human impact land type, protected natural area Human impact land type include mining lease, industrial land, friendship
Logical transport land used, Tourism land, residential land and Farm Land, weighted value takes 0.42 respectively, 0.22,0.16,0.12,
0.05 and 0.03；It is i-th7Kind Human impact land type is in functional areas j7In area.
8. a kind of effectiveness assessment system of ecological class protected natural area according to claim 7, which is characterized in that institute
Stating effectiveness index computing module includes:
Landscape pattern's efficacy data substitutes into unit, for landscape pattern's Effective exponent to be substituted into landscape pattern's validity
Exponential model F (x1(the 1+EXP (- 1.635-0.757*x of)=15/1)) in, obtain landscape pattern's effectiveness index；
Turnover ecosystem efficacy data substitutes into unit, for the turnover ecosystem Effective exponent to be substituted into ecosystem
System conversion Effective exponent model F (x2(the 1+EXP (- 1.635-1.557*x of)=15/2)) in, obtain turnover ecosystem protection
Effect index；
Ecosystem quality efficacy data substitutes into unit, is used for when the ecosystem is non-wetlands ecosystems, by the life
State mass of system Effective exponent substitutes into non-wetland ecosystem quality Effective exponent model F (x31)=18/ (1+EXP (-
1.21-0.69*x31)) in, obtain non-wetland ecosystem quality effectiveness index；When the ecosystem is wetlands ecosystems
When, the ecosystem quality Effective exponent is substituted into Wetland Area variability index model F (x32)=10/ (1+EXP (-
1.35-0.65*x32)) in, obtain Wetland Area effectiveness index；
Forest Growth efficacy data substitutes into unit, for the Forest Growth Effective exponent to be substituted into Forest Growth validity
Exponential model F (x4(the 1+EXP (- 0.73-0.54*x of)=10/4)) in, obtain Forest Growth effectiveness index；
Vegetation stablizes efficacy data and substitutes into unit, stablizes validity for the vegetation to be stablized Effective exponent substitution vegetation
Exponential model F (x5(the 1+EXP (- 1.83-0.81*x of)=10/5)) in, it obtains vegetation and stablizes effectiveness index；
Protected object species quantity efficacy data substitutes into unit, is used for the protected object species quantity Effective exponent generation
Enter protected object species quantity Effective exponent model F (x6(the 1+EXP (- 1.23-0.69*x of)=18/6)) in, obtain protection pair
As species quantity effectiveness index；
Interference variations efficacy data substitutes into unit, for the interference variations Effective exponent to be substituted into interference variations validity
Exponential model F (x7(the 1+EXP (- 1.135+0.71*x of)=14/7)) in, obtain mankind's activity Annoyance Index.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811581119.7A CN109657994B (en) | 2018-12-24 | 2018-12-24 | Method and system for evaluating protection effect of ecological natural protection land |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811581119.7A CN109657994B (en) | 2018-12-24 | 2018-12-24 | Method and system for evaluating protection effect of ecological natural protection land |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109657994A true CN109657994A (en) | 2019-04-19 |
| CN109657994B CN109657994B (en) | 2020-11-10 |
Family
ID=66116427
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811581119.7A Active CN109657994B (en) | 2018-12-24 | 2018-12-24 | Method and system for evaluating protection effect of ecological natural protection land |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109657994B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111260262A (en) * | 2020-03-17 | 2020-06-09 | 国家林业和草原局华东调查规划设计院 | Natural conservation ground protection effect evaluation method and system and application thereof |
| CN111311114A (en) * | 2020-03-11 | 2020-06-19 | 汇明科技(江苏)有限公司 | Method for evaluating land resource restoration effect |
| CN112836969A (en) * | 2021-01-12 | 2021-05-25 | 浙江易智信息技术有限公司 | Remote sensing technology-based river ecological quality evaluation method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100070320A1 (en) * | 2008-09-15 | 2010-03-18 | Parametrix, Inc. | Accounting tool for measuring ecosystem service functional performance at a particular site |
| CN103425892A (en) * | 2013-08-27 | 2013-12-04 | 北京林业大学 | Protective effect computing method for landscape in nature protection area |
| CN103778338A (en) * | 2014-01-26 | 2014-05-07 | 北京林业大学 | Method for calculating landscape protection value of land nature reserves |
| CN108305204A (en) * | 2018-02-12 | 2018-07-20 | 中国科学院东北地理与农业生态研究所 | Wetland Biodiversity effectiveness regional correlation appraisal procedure based on basis of reference |
-
2018
- 2018-12-24 CN CN201811581119.7A patent/CN109657994B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100070320A1 (en) * | 2008-09-15 | 2010-03-18 | Parametrix, Inc. | Accounting tool for measuring ecosystem service functional performance at a particular site |
| CN103425892A (en) * | 2013-08-27 | 2013-12-04 | 北京林业大学 | Protective effect computing method for landscape in nature protection area |
| CN103778338A (en) * | 2014-01-26 | 2014-05-07 | 北京林业大学 | Method for calculating landscape protection value of land nature reserves |
| CN108305204A (en) * | 2018-02-12 | 2018-07-20 | 中国科学院东北地理与农业生态研究所 | Wetland Biodiversity effectiveness regional correlation appraisal procedure based on basis of reference |
Non-Patent Citations (2)
| Title |
|---|
| 郑骁喆 等: "拉市海高原湿地省级自然保护区保护成效评估研究", 《林业资源管理》 * |
| 韦惠兰 等: "秦岭自然保护区保护成效评估", 《生态经济》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111311114A (en) * | 2020-03-11 | 2020-06-19 | 汇明科技(江苏)有限公司 | Method for evaluating land resource restoration effect |
| CN111260262A (en) * | 2020-03-17 | 2020-06-09 | 国家林业和草原局华东调查规划设计院 | Natural conservation ground protection effect evaluation method and system and application thereof |
| CN112836969A (en) * | 2021-01-12 | 2021-05-25 | 浙江易智信息技术有限公司 | Remote sensing technology-based river ecological quality evaluation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109657994B (en) | 2020-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Identification of ecosystem service bundles and driving factors in Beijing and its surrounding areas | |
| Hepinstall-Cymerman et al. | Urban growth patterns and growth management boundaries in the Central Puget Sound, Washington, 1986–2007 | |
| Theobald | Landscape patterns of exurban growth in the USA from 1980 to 2020 | |
| Turner et al. | Wet/dry mapping: using citizen scientists to monitor the extent of perennial surface flow in dryland regions | |
| Geneletti | Ecological evaluation for environmental impact assessment | |
| Duffin et al. | Relative pollen productivity and fall speed estimates for southern African savanna taxa | |
| Joosten et al. | MoorFutures®: integration of additional ecosystem services (including biodiversity) into carbon credits-standard, methodology and transferability to other regions | |
| CN109657994A (en) | A kind of the effectiveness appraisal procedure and system of ecology class protected natural area | |
| Coops et al. | Exploring the relative importance of satellite-derived descriptors of production, topography and land cover for predicting breeding bird species richness over Ontario, Canada | |
| CN111260262A (en) | Natural conservation ground protection effect evaluation method and system and application thereof | |
| Yao et al. | Evaluation of arable land suitability based on the suitability function-A case study of the Qinghai-Tibet Plateau | |
| CN112651548A (en) | Evaluation and identification method for plateau lakeside ecological landscape restoration planning | |
| CN114663782A (en) | Forestry carbon source remittance method for forestry environment data | |
| Yoshikai et al. | Predicting mangrove forest dynamics across a soil salinity gradient using an individual-based vegetation model linked with plant hydraulics | |
| Snelder et al. | Land-use suitability is not an intrinsic property of a land parcel | |
| Isola et al. | Green Infrastructure and Regional Planning: An Operational Framework: An Operational Framework | |
| Weber et al. | Use of LiDAR and supplemental data to estimate forest maturity in Charles County, MD, USA | |
| Hychka et al. | Enhancing a landscape assessment with intensive data: a case study in the Upper Juniata watershed | |
| Ji | Wetland and water resource modeling and assessment: a watershed perspective | |
| Carver et al. | Wildness Study in the Cairngorms National Park | |
| Osipov et al. | Modelling of the suitability of lands to the agrarian use and their resistance to negative processes | |
| Bjorkman et al. | Biomass, Carbon Sequestration, and Avoided Emissions: Assessing the Role of Urban Trees in California | |
| Seid et al. | Spatiotemporal dynamics of land use/land cover changes and its drivers in Bilate watershed, central rift valley, Ethiopia | |
| Xu et al. | Land-Use Change and Driving Force Analysis of Wetland in Poyang Lake Based on Remote Sensing | |
| Niamir | Grazing intensity and ecological change in eastern Senegal: Implications for the monitoring of Sahelian rangelands |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221222 Address after: 518000 211C, Floor 2, West Block, Xincheng Building, No. 1025, Shennan Middle Road, Badeng Community, Nanyuan Street, Futian District, Shenzhen, Guangdong Patentee after: Shenzhen Ring Thumb Electronic Technology Co.,Ltd. Address before: No. 300, Panlong temple, Kunming, Yunnan Province Patentee before: SOUTHWEST FORESTRY University |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230504 Address after: Zone 13-1-101, Zone A, Longxi No.1, Panlong District, Kunming City, Yunnan Province, 650051 Patentee after: Yunnan outdoor Map Technology Co.,Ltd. Address before: 518000 211C, Floor 2, West Block, Xincheng Building, No. 1025, Shennan Middle Road, Badeng Community, Nanyuan Street, Futian District, Shenzhen, Guangdong Patentee before: Shenzhen Ring Thumb Electronic Technology Co.,Ltd. |